In:
The Astrophysical Journal Supplement Series, American Astronomical Society, Vol. 252, No. 2 ( 2021-02-01), p. 14-
Abstract:
Self-gravity calculations for 3D are expensive in terms of computational time. Several methods exist for this computation, for example multigrid and spectral methods. Unfortunately, these approaches require the imposition of boundary conditions, which can be either numerically expensive (direct Newtonian sums), artificial (periodicity assumptions), or potentially imprecise (multipolar expansions). In this work we present a novel direct numerical method to calculate the gravitational potential and forces by solving the Poisson equation without the need to prescribe artificial boundary conditions; this method, despite being direct, turns out to be efficient due to the possibility of using a fast Fourier transform for its implementation. For a grid having N zones in each dimension, the computational complexity of the method presented here is , which is comparable with multigrid methods under no consideration of boundary settings. Finally, a numerical study shows this proposed method can achieve second order for calculations of both potential and forces.
Type of Medium:
Online Resource
ISSN:
0067-0049
,
1538-4365
DOI:
10.3847/1538-4365/abca97
Language:
Unknown
Publisher:
American Astronomical Society
Publication Date:
2021
detail.hit.zdb_id:
2961-0
detail.hit.zdb_id:
2006860-8
SSG:
16,12
Permalink