arrow
Return

The Athena plus plus Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-gravity

delete2023-04-27
delete9
delete
OA
AI
K
Kengo Tomida *
J
James M. Stone
DOI:10.3847/1538-4365/acc2c0delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
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 Message Passing Interface and OpenMP is adopted, and we present results of tests demonstrating the accuracy and scaling of the methods. On a uniform grid, we show that multigrid significantly outperforms methods based on fast Fourier transforms, and requires only a small fraction of the computing 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.
Keywords:
POISSON SOLVER
MAGNETOHYDRODYNAMICS
CODE
TREE
SIMULATIONS
ACCRETION
EXTENSION
TRANSPORT
COLLAPSE
TURBULENCE

Journal

Astrophysical Journal Supplement Series cover
Astrophysical Journal Supplement Series
IF:
8.5
Papers:
5.7K
Citations:
4.2W

Organization

T
tohoku university
Scholars:
4.3W
Papers: 3.6W
Citations: 31
I
institute for advanced study - usa
Scholars:
933
Papers: 1.4K
Citations: 15