Parallel implementation of an adaptive and parameter-free N-body integrator

DOI10.1016/J.CPC.2011.01.014MaRDI QIDQ538625zbMATH OpenOpenAlexFDO

Authors C. David Pruett, William H. Ingham, Ralph D. Herman

Publication date 25 May 2011

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1016/j.cpc.2011.01.014

zbMATH Keywords

long time evolutionsolar systemTaylor series method

Mathematics Subject Classification ID

Parallel numerical computation (65Y05) Computational methods for problems pertaining to mechanics of particles and systems (70-08) (n)-body problems (70F10) Celestial mechanics (70F15)

Review text

This paper is concerned with numerical integration of \(N\)-body problems in a parallel implementation, where \(N\) is assumed to be large. The authors consider only \(N\)-body systems under gravitational interactions, and the original \(6N\)-dimensional system in positions \({\mathbf{x}}_i\), \(i=1, \dots ,N\), and velocities \({\mathbf{v}}_i\), \(i=1, \dots ,N\), is extended to an autonomous polynomial system of dimension \( 6N +N(N-1)/2 \) that includes the inverse distances \( S_{ij}= | {\mathbf{x}}_i -{\mathbf{x}}_j |^{-1}\), \( 1 \leq i < j \leq N\). The time evolution of such a system is carried out by a Taylor series expansion of all variables up to order \(m \leq 28\), and variable step sizes are selected dynamically according to a local error estimate algorithm. Parallelization is achieved by a suitable distribution of the bodies in the processors. In particular, a detailed study of the parallelization with 4 processors is given estimating the operation count per processor. In the second part, some numerical experiments are presented to show the advantages and shortcomings of the proposed method. After some introductory numerical experiments with a two-body problem, the authors examine as a test problem long time evolution of the solar system under gravitational forces and some ``swarms of particles with \(N\) between 18 and 480. The results are compared with other standard Runge-Kutta methods.

Recommendations

• The development of a parallel \(N\)-body code for the Edinburgh Concurrent Supercomputer
• N-body simulations
• A Time-Reversible, Regularized, Switching Integrator for the N-Body Problem
• scientific article; zbMATH DE number 663846
• Parallel/vector integration methods for dynamical astronomy

Cites work

• An adaptive N-body algorithm of optimal order.
• FOURTH ORDER SYMPLECTIC INTEGRATION WITH REDUCED PHASE ERROR
• scientific article; zbMATH DE number 1745051 (Why is no real title available?)
• scientific article; zbMATH DE number 2169895 (Why is no real title available?)
• Solution of theN-body problem with recurrent power series

Cited in

(7)

• Parallel \(N\)-body simulation based on the PM and P3M methods using multigrid schemes in conjunction with generic approximate sparse inverses
• Requirements of a package for N-body simulations of the solar system
• An adaptive N-body algorithm of optimal order.
• Parallel/vector integration methods for dynamical astronomy
• Analysis and parallel implementation of a forced N-body problem
• Provably Good Partitioning and Load Balancing Algorithms for Parallel Adaptive N-Body Simulation
• scientific article; zbMATH DE number 970129 (Why is no real title available?)

This page was built for publication: Parallel implementation of an adaptive and parameter-free \(N\)-body integrator

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q538625)