Computing the gravitational potential on nested meshes using the convolution method

Abstract

Aims. Our aim is to derive a fast and accurate method for computing the gravitational potential of astrophysical objects with high contrasts in density, for which nested or adaptive meshes are required. Methods. We present an extension of the convolution method for computing the gravitational potential to the nested Cartesian grids. The method makes use of the convolution theorem to compute the gravitational potential using its integral form. Results. A comparison of our method with the iterative outside-in conjugate gradient and generalized minimal residual methods for solving the Poisson equation using nonspherically symmetric density configurations has shown a comparable performance in terms of the errors relative to the analytic solutions. However, the convolution method is characterized by several advantages and outperforms the considered iterative methods by factors 10--200 in terms of the runtime, especially when graphics processor units are utilized. The convolution method also shows an overall second-order convergence, except for the errors at the grid interfaces where the convergence is linear. Conclusions. High computational speed and ease in implementation can make the convolution method a preferred choice when using a large number of nested grids. The convolution method, however, becomes more computationally costly if the dipole moments of tightly spaced gravitating objects are to be considered at coarser grids.

arXiv commentAccepted for publication in Astronomy & Astrophysics

References in corpus (14)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• The Physics of Star Formation
• Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Non-Ideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks
• A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
• Treating gravity in thin disk simulations
• Collapse and Fragmentation of Rotating Magnetized Clouds. I. Magnetic Flux - Spin Relation
• Self-gravitational Magnetohydrodynamics with Adaptive Mesh Refinement for Protostellar Collapse
• Using the PPML approach for constructing a low-dissipation, operator-splitting scheme for numerical simulations of hydrodynamic flows
• Protostellar collapse: the conditions to form dust rich protoplanetary disks
• A Fast Algorithm for Solving the Poisson Equation on a Nested Grid
• The FRIGG project: From intermediate galactic scales to self-gravitating cores
• Accretion bursts in magnetized gas-dust protoplanetary disks
• Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
• Self-gravitational Force Calculation of Second Order Accuracy Using Multigrid Method on Nested Grids

Cited by in corpus (3)

• Accelerating Fortran Codes: A Method for Integrating Coarray Fortran with CUDA Fortran and OpenMP
• Primordial dust rings, hidden dust mass, and the first generation of planetesimals in gravitationally unstable protoplanetary disks
• Dust enrichment and growth in the earliest stages of protoplanetary disk formation