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

Bernhard Müller, Conrad Chan

Research output: Contribution to journalArticleResearchpeer-review

2 Citations (Scopus)

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 3D, the operation count of the new method is competitive with a naive implementation of the truncated spherical harmonics expansion with N ≈ 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.

Original languageEnglish
Article number43
Pages (from-to)1-9
Number of pages9
JournalAstrophysical Journal
Volume870
Issue number1
DOIs
Publication statusPublished - 1 Jan 2019

Keywords

  • gravitation
  • methods: numerical

Cite this

@article{7de6dec9e6d842c49162b9f9022a3610,
title = "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 3D, the operation count of the new method is competitive with a naive implementation of the truncated spherical harmonics expansion with N ℓ ≈ 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.",
keywords = "gravitation, methods: numerical",
author = "Bernhard M{\"u}ller and Conrad Chan",
year = "2019",
month = "1",
day = "1",
doi = "10.3847/1538-4357/aaf100",
language = "English",
volume = "870",
pages = "1--9",
journal = "The Astrophysical Journal",
issn = "0004-637X",
publisher = "American Astronomical Society",
number = "1",

}

An FFT-based Solution Method for the Poisson Equation on 3D Spherical Polar Grids. / Müller, Bernhard; Chan, Conrad.

In: Astrophysical Journal, Vol. 870, No. 1, 43, 01.01.2019, p. 1-9.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

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

AU - Müller, Bernhard

AU - Chan, Conrad

PY - 2019/1/1

Y1 - 2019/1/1

N2 - 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 3D, the operation count of the new method is competitive with a naive implementation of the truncated spherical harmonics expansion with N ℓ ≈ 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.

AB - 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 3D, the operation count of the new method is competitive with a naive implementation of the truncated spherical harmonics expansion with N ℓ ≈ 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.

KW - gravitation

KW - methods: numerical

UR - http://www.scopus.com/inward/record.url?scp=85059838062&partnerID=8YFLogxK

U2 - 10.3847/1538-4357/aaf100

DO - 10.3847/1538-4357/aaf100

M3 - Article

VL - 870

SP - 1

EP - 9

JO - The Astrophysical Journal

JF - The Astrophysical Journal

SN - 0004-637X

IS - 1

M1 - 43

ER -