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.
35 pages, 32 figures
References in corpus (8)
- PLUTO: a Numerical Code for Computational Astrophysics
- Athena: A New Code for Astrophysical MHD
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
- Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
- Cholla : A New Massively-Parallel Hydrodynamics Code For Astrophysical Simulation
- New numerical solver for flows at various Mach numbers
- WOMBAT: A Scalable and High Performance Astrophysical MHD Code
- APSARA: A multi-dimensional unsplit fourth-order explicit Eulerian hydrodynamics code for arbitrary curvilinear grids
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