The Athena++ Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-Gravity

Abstract

We describe the implementation of multigrid solvers in the Athena++ adaptive mesh refinement (AMR) framework and their application to the solution of the Poisson equation for self-gravity. The new solvers are built on top of the AMR hierarchy and TaskList framework of Athena++ for efficient parallelization. We adopt a conservative formulation for the Laplacian operator that avoids artificial accelerations at level boundaries. Periodic, fixed, and zero-gradient boundary conditions are implemented, as well as open boundary conditions based on a multipole expansion. Hybrid parallelization using both MPI and OpenMP is adopted, and we present results of tests demonstrating the accuracy and scaling of the methods. On a uniform grid we show multigrid significantly outperforms methods based on FFTs, and requires only a small fraction of the compute time required by the (highly optimized) magnetohydrodynamic solver in Athena++. As a demonstration of the capabilities of the methods, we present the results of a test calculation of magnetized protostellar collapse on an adaptive mesh.

28 pages, 13 figures, submitted to AAS Journals

References in corpus (13)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• Magnetic processes in a collapsing dense core. I Accretion and Ejection
• Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Non-Ideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks
• A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• Grand Design Spiral Arms in A Young Forming Circumstellar Disk
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity
• An Extension of the Athena++ Framework for Fully Conservative Self-Gravitating Hydrodynamics
• Unravelling cosmic velocity flows: a Helmholtz-Hodge decomposition algorithm for cosmological simulations
• Self-gravitational Force Calculation of Second Order Accuracy Using Multigrid Method on Nested Grids