A self-gravity module for the PLUTO code
arXiv:2306.05332 · doi:10.3847/1538-4365/aced0a
Abstract
We present a novel implementation of an iterative solver for the solution of the Poisson equation in the PLUTO code for astrophysical fluid dynamics. Our solver relies on a relaxation method in which convergence is sought as the steady-state solution of a parabolic equation, whose time-discretization is governed by the \textit{Runge-Kutta-Legendre} (RKL) method. Our findings indicate that the RKL-based Poisson solver, which is both fully parallel and rapidly convergent, has the potential to serve as a practical alternative to conventional iterative solvers such as the \textit{Gauss-Seidel} (GS) and \textit{successive over-relaxation} (SOR) methods. Additionally, it can mitigate some of the drawbacks of these traditional techniques. We incorporate our algorithm into a multigrid solver to provide a simple and efficient gravity solver that can be used to obtain the gravitational potentials in self-gravitational hydrodynamics. We test our implementation against a broad range of standard self-gravitating astrophysical problems designed to examine different aspects of the code. We demonstrate that the results match excellently with the analytical predictions (when available), and the findings of similar previous studies.
arXiv commentPublished in ApJS
References in corpus (27)
- Array Programming with NumPy
- The cosmological simulation code GADGET-2
- Theory of Star Formation
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- GADGET: A code for collisionless and gasdynamical cosmological simulations
- PLUTO: a Numerical Code for Computational Astrophysics
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- Density, Velocity, and Magnetic Field Structure in Turbulent Molecular Cloud Models
- Gasoline: An adaptable implementation of TreeSPH
- The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics
- Investigating fragmentation conditions in self-gravitating accretion discs
- CMasher: Scientific colormaps for making accessible, informative and 'cmashing' plots
- CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity
- Dimension as a Key to the Neutrino Mechanism of Core-Collapse Supernova Explosions
- A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
- A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
- An Improved Multipole Approximation for Self-Gravity and Its Importance for Core-Collapse Supernova Simulations
- Resolution requirements for simulating gravitational fragmentation using SPH
- GANDALF - Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids
- Impact of relativistic jets on the star formation rate: a turbulence-regulated framework
- Why does the Jeans Swindle work?
- White Dwarf Mergers on Adaptive Meshes I. Methodology and Code Verification
- A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity
- Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping
- An Extension of the Athena++ Framework for Fully Conservative Self-Gravitating Hydrodynamics
- Accelerating self-gravitating hydrodynamics simulations with adaptive force updates
- Self-gravitational Force Calculation of Second Order Accuracy Using Multigrid Method on Nested Grids