An FFT-based Solution Method for the Poisson Equation on 3D Spherical Polar Grids

Abstract

The solution of the Poisson equation is a ubiquitous problem in computational astrophysics. Most notably, the treatment of self-gravitating flows involves the Poisson equation for the gravitational field. In hydrodynamics codes using spherical polar grids, one often resorts to a truncated spherical harmonics expansion for an approximate solution. Here we present a non-iterative method that is similar in spirit, but uses the full set of eigenfunctions of the discretized Laplacian to obtain an exact solution of the discretized Poisson equation. This allows the solver to handle density distributions for which the truncated multipole expansion fails, such as off-center point masses. In three dimensions, the operation count of the new method is competitive with a naive implementation of the truncated spherical harmonics expansion with $N_{\ell }\approx 15$ multipoles. We also discuss the parallel implementation of the algorithm. The serial code and a template for the parallel solver are made publicly available.

9 pages, 4 figures, accepted for publication in ApJ. Added convergence tests and included some other minor revisions

References in corpus (8)

• Delayed neutrino-driven supernova explosions aided by the standing accretion-shock instability
• The isotropic diffusion source approximation for supernova neutrino transport
• The Dynamics of Neutrino-Driven Supernova Explosions after Shock Revival in 2D and 3D
• Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue
• Crucial Physical Dependencies of the Core-Collapse Supernova Mechanism
• An axis-free overset grid in spherical polar coordinates for simulating 3D self-gravitating flows
• APSARA: A multi-dimensional unsplit fourth-order explicit Eulerian hydrodynamics code for arbitrary curvilinear grids
• Parallelized Solution Method of the Three-dimensional Gravitational Potential on the Yin-Yang Grid

Cited by in corpus (3)

• Hydrodynamics of core-collapse supernovae and their progenitors
• A Fast Poisson Solver of Second-Order Accuracy for Isolated Systems in Three-Dimensional Cartesian and Cylindrical Coordinates
• A Generalized Solution Method for Parallelized Computation of the Three-dimensional Gravitational Potential on a Multi-patch grid in Spherical Geometry