TU Darmstadt / ULB / TUbiblio

Parallel summation of symmetric inter-particle forces in smoothed particle hydrodynamics

Willkomm, J. and Bücker, H. M.
Griebel, M. and Schweitzer, M. A. (eds.) (2011):
Parallel summation of symmetric inter-particle forces in smoothed particle hydrodynamics.
In: Meshfree Methods for Partial Differential Equations V, Berlin, Springer, pp. 235-248, [Book Section]

Abstract

In the smoothed particle hydrodynamics (SPH) method, the forces between all particles are efficiently summed up in a serial environment as follows. Each pair of particles is considered once. The resulting inter-particle force is computed and then the contributions to both particles are updated, taking into account the symmetry of the problem. This algorithm is difficult to parallelise when concurrently accessing the same memory location in a multi-threaded process. We develop a parallel 1D summation algorithm consisting of two passes on a Cartesian grid of cells in which particles move freely. In a first pass, we consider all cells with an even index n. We compute the inter-particle forces for all pairs of particles located in cell n and those pairs where a particle is located in cell n and another particle is located in cell n + 1. Each cell n is handled by a different thread since no data race can occur. In a second pass, we do the same for cells with odd n. This way, all interparticle forces are computed. We generalise this algorithm to 2D, 3D, and arbitrarily high dimensions and report performance results on three different shared-memory platforms using the OpenMP programming paradigm.

Item Type: Book Section
Erschienen: 2011
Editors: Griebel, M. and Schweitzer, M. A.
Creators: Willkomm, J. and Bücker, H. M.
Title: Parallel summation of symmetric inter-particle forces in smoothed particle hydrodynamics
Language: English
Abstract:

In the smoothed particle hydrodynamics (SPH) method, the forces between all particles are efficiently summed up in a serial environment as follows. Each pair of particles is considered once. The resulting inter-particle force is computed and then the contributions to both particles are updated, taking into account the symmetry of the problem. This algorithm is difficult to parallelise when concurrently accessing the same memory location in a multi-threaded process. We develop a parallel 1D summation algorithm consisting of two passes on a Cartesian grid of cells in which particles move freely. In a first pass, we consider all cells with an even index n. We compute the inter-particle forces for all pairs of particles located in cell n and those pairs where a particle is located in cell n and another particle is located in cell n + 1. Each cell n is handled by a different thread since no data race can occur. In a second pass, we do the same for cells with odd n. This way, all interparticle forces are computed. We generalise this algorithm to 2D, 3D, and arbitrarily high dimensions and report performance results on three different shared-memory platforms using the OpenMP programming paradigm.

Title of Book: Meshfree Methods for Partial Differential Equations V
Series Name: Lecture Notes in Computational Science and Engineering
Volume: 79
Place of Publication: Berlin
Publisher: Springer
Uncontrolled Keywords: Smoothed particle hydrodynamics; SPH; Parallelization; OpenMP
Divisions: 20 Department of Computer Science > Scientific Computing
20 Department of Computer Science
Date Deposited: 16 May 2013 11:28
Export:
Suche nach Titel in: TUfind oder in Google

Optionen (nur für Redakteure)

View Item View Item