Gravitational Self-force Errors of Poisson Solvers on Adaptively Refined Meshes
arXiv:2104.10170 · doi:10.3847/1538-4365/abe86f
Abstract
An error in the gravitational force that the source of gravity induces on itself (a self-force error) violates both the conservation of linear momentum and the conservation of energy. If such errors are present in a self-gravitating system and are not sufficiently random to average out, the obtained numerical solution will become progressively more unphysical with time: the system will acquire or lose momentum and energy due to numerical effects. In this paper, we demonstrate how self-force errors can arise in the case where self-gravity is solved on an adaptively refined mesh when the refinement is nonuniform. We provide the analytical expression for the self-force error and numerical examples that demonstrate such self-force errors in idealized settings. We also show how these errors can be corrected to an arbitrary order by straightforward addition of correction terms at the refinement boundaries.
arXiv comment8 pages, 5 figures, ApJS accepted
References in corpus (10)
- The cosmological simulation code GADGET-2
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
- Enzo: An Adaptive Mesh Refinement Code for Astrophysics
- The Arepo public code release
- Effects of Baryons and Dissipation on the Matter Power Spectrum
- Constrained Simulations of the Real Universe: II. Observational Signatures of Intergalactic Gas in the Local Supercluster Region
- A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
- A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
- An adaptive particle-mesh gravity solver for ENZO