Efficient, Divergence-Free, High Order MHD on 3D Spherical Meshes with Optimal Geodesic Meshing

Abstract

There is a great need in several areas of astrophysics and space-physics to carry out high order of accuracy, divergence-free MHD simulations on spherical meshes. This requires us to pay careful attention to the interplay between mesh quality and numerical algorithms. Methods have been designed that fundamentally integrate high order isoparametric mappings with the other high accuracy algorithms that are needed for divergence-free MHD simulations on geodesic meshes. The goal of this paper is to document such algorithms that are implemented in the geodesic mesh version of the RIEMANN code. The fluid variables are reconstructed using a special kind of WENO-AO algorithm that integrates the mesh geometry into the reconstruction process from the ground-up. A novel divergence-free reconstruction strategy for the magnetic field that performs efficiently at all orders, even on isoparametrically mapped meshes, is then presented. The MHD equations are evolved in space and time using a novel ADER predictor algorithm that is efficiently adapted to the isoparametrically mapped geometry. The application of one-dimensional and multidimensional Riemann solvers at suitable locations on the mesh then provides the corrector step. The corrector step for the magnetic field uses a Yee-type staggering of magnetic fields. This results in a scheme with divergence-free update for the magnetic field. The use of ADER enables a one-step update which only requires one messaging operation per complete timestep. This is very beneficial for parallel processing. Several accuracy tests are presented as are stringent test problems. PetaScale performance is also demonstrated on the largest available supercomputers.

arXiv commentAccepted for publication in MNRAS

References in corpus (17)

• Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
• Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
• On the Divergence-Free Condition in Godunov-Type Schemes for Ideal Magnetohydrodynamics: the Upwind Constrained Transport Method
• Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
• Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
• A Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows
• Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
• ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement
• A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
• Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics
• On the accuracy of WENO and CWENO reconstructions of third order on nonuniform meshes
• Central WENO schemes for hyperbolic conservation laws on fixed and moving unstructured meshes
• Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
• Three Dimensional MHD Simulations of Accretion to an Inclined Rotator: The "Cubed Sphere" Method
• Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes
• A New Code for Nonlinear Force-Free Field Extrapolation of the Global Corona
• MHD modeling on geodesic grids

Cited by in corpus (4)

• On GLM curl cleaning for a first order reduction of the CCZ4 formulation of the Einstein field equations
• Systematic construction of upwind constrained transport schemes for MHD
• Modeling Magnetically Channeled Winds in 3D: I. Isothermal Simulations of a Magnetic O Supergiant
• Technologies for supporting high-order geodesic mesh frameworks for computational astrophysics and space sciences