A Software Framework for the Portable Parallelization of Particle-Mesh Simulations

We present a software framework for the transparent and portable parallelization of simulations using particle-mesh methods. Particles are used to transport physical properties and a mesh is required in order to reinitialize the distorted particle locations, ensuring the convergence of the method. Field quantities are computed on the particles using fast multipole methods or by discretizing and solving the governing equations on the mesh. This combination of meshes and particles presents a challenging set of parallelization issues. The present library addresses these issues for a wide range of applications, and it enables orders of magnitude increase in the number of computational elements employed in particle methods. We demonstrate the performance and scalability of the library on several problems, including the first-ever billion particle simulation of diffusion in real biological cell geometries.

[1]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[2]  P. Degond,et al.  The weighted particle method for convection-diffusion equations , 1989 .

[3]  Petros Koumoutsakos,et al.  Remeshed smoothed particle hydrodynamics for the simulation of viscous and heat conducting flows , 2002 .

[4]  Jens Honore Walther,et al.  An influence matrix particle-particle particle-mesh algorithm with exact particle-particle correction , 2003 .

[5]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[6]  Petros Koumoutsakos,et al.  Effects of organelle shape on fluorescence recovery after photobleaching. , 2005, Biophysical journal.

[7]  Ivo F. Sbalzarini,et al.  PPM - A highly efficient parallel particle-mesh library for the simulation of continuum systems , 2006, J. Comput. Phys..

[8]  J. Monaghan,et al.  Extrapolating B splines for interpolation , 1985 .

[9]  V. Rokhlin,et al.  Rapid Evaluation of Potential Fields in Three Dimensions , 1988 .

[10]  Dongming Wang,et al.  Automated Deduction in Geometry , 1996, Lecture Notes in Computer Science.

[11]  L. Verlet Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[12]  Petros Koumoutsakos,et al.  Three-dimensional vortex methods for particle-laden flows with two-way coupling , 2001 .

[13]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[14]  Tim Colonius,et al.  A general deterministic treatment of derivatives in particle methods , 2002 .

[15]  V. Springel,et al.  GADGET: a code for collisionless and gasdynamical cosmological simulations , 2000, astro-ph/0003162.

[16]  Joel H. Saltz,et al.  Adaptive runtime support for direct simulation Monte Carlo methods on distributed memory architectures , 1994, Proceedings of IEEE Scalable High Performance Computing Conference.

[17]  Michael Bergdorf,et al.  Multilevel Adaptive Particle Methods for Convection-Diffusion Equations , 2005, Multiscale Model. Simul..

[18]  J. Williamson Low-storage Runge-Kutta schemes , 1980 .

[19]  P. Koumoutsakos MULTISCALE FLOW SIMULATIONS USING PARTICLES , 2005 .

[20]  P. Gremaud,et al.  Super-time-stepping acceleration of explicit schemes for parabolic problems , 1996 .

[21]  Petros Koumoutsakos,et al.  Simulations of (an)isotropic diffusion on curved biological surfaces. , 2006, Biophysical journal.

[22]  Petros Koumoutsakos,et al.  On the Water−Carbon Interaction for Use in Molecular Dynamics Simulations of Graphite and Carbon Nanotubes , 2003 .

[23]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .