WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
arXiv:1812.04278 · doi:10.3847/1538-4365/ab09fb
Abstract
Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine "Hazel Hen", with a sustained performance of about 5 percent of peak at scale.
arXiv comment35 pages, 32 figures
References in corpus (40)
- 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
- PLUTO: a Numerical Code for Computational Astrophysics
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- Athena: A New Code for Astrophysical MHD
- Enzo: An Adaptive Mesh Refinement Code for Astrophysics
- Adaptive Refinement Tree - a new high-resolution N-body code for cosmological simulations
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport
- A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
- Nyx: A Massively Parallel AMR Code for Computational Cosmology
- 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
- High Order Upwind Schemes for Multidimensional Magnetohydrodynamics
- Accurate, Meshless Methods for Magneto-Hydrodynamics
- A Solution Accurate, Efficient and Stable Unsplit Staggered Mesh Scheme for Three Dimensional Magnetohydrodynamics
- Moving mesh cosmology: the hydrodynamics of galaxy formation
- High-order conservative finite difference GLM-MHD schemes for cell-centered MHD
- A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
- The MHD Kelvin-Helmholtz Instability: A Two-Dimensional Numerical Study
- Computational Eulerian Hydrodynamics and Galilean Invariance
- A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
- Cholla : A New Massively-Parallel Hydrodynamics Code For Astrophysical Simulation
- GAMER-2: a GPU-accelerated adaptive mesh refinement code -- accuracy, performance, and scalability
- New numerical solver for flows at various Mach numbers
- An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations
- A fourth-order accurate finite volume method for ideal MHD via upwind constrained transport
- Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
- A Hybrid Cosmological Hydrodynamic/N-body Code Based on a Weighted Essentially Non-Oscillatory Scheme
- A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
- The Recent Growth History of the Fornax Cluster Derived from Simultaneous Sloshing and Gas Stripping: Simulating the Infall of NGC 1404
- High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
- DISPATCH: A Numerical Simulation Framework for the Exa-scale Era. I. Fundamentals
- Multi-Dimensional, Compressible Viscous Flow on a Moving Voronoi Mesh
- WOMBAT: A Scalable and High Performance Astrophysical MHD Code
- APSARA: A multi-dimensional unsplit fourth-order explicit Eulerian hydrodynamics code for arbitrary curvilinear grids
- Planet-disc interactions with Discontinuous Galerkin Methods using GPUs
- Higher Order Central Schemes for Multi-dimensional Hyperbolic Problems
Cited by in corpus (3)
- HOW-MHD: A High-Order WENO-Based Magnetohydrodynamic Code with a High-Order Constrained Transport Algorithm for Astrophysical Applications
- Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows
- An Alternative Finite Difference WENO-like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems