AN ADAPTIVE PARTICLE-MESH GRAVITY SOLVER FOR ENZO
- Argelander Institute für Astronomie, Bonn Universität, D-53121 Bonn (Germany)
- Department of Astronomy, Columbia University, New York, NY 10027 (United States)
We describe and implement an adaptive particle-mesh algorithm to solve the Poisson equation for grid-based hydrodynamics codes with nested grids. The algorithm is implemented and extensively tested within the astrophysical code Enzo against the multigrid solver available by default. We find that while both algorithms show similar accuracy for smooth mass distributions, the adaptive particle-mesh algorithm is more accurate for the case of point masses, and is generally less noisy. We also demonstrate that the two-body problem can be solved accurately in a configuration with nested grids. In addition, we discuss the effect of subcycling, and demonstrate that evolving all the levels with the same timestep yields even greater precision.
- OSTI ID:
- 22340105
- Journal Information:
- Astrophysical Journal, Supplement Series, Vol. 215, Issue 1; Other Information: Country of input: International Atomic Energy Agency (IAEA); ISSN 0067-0049
- Country of Publication:
- United States
- Language:
- English
Similar Records
Elliptic Solvers for Adaptive Mesh Refinement Grids
Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway