Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics

Abstract

In this work we develop a class of high-order finite difference weighted essentially non-oscillatory (FD-WENO) schemes for solving the ideal magnetohydrodynamic (MHD) equations in 2D and 3D. The philosophy of this work is to use efficient high-order WENO spatial discretizations with high-order strong stability-preserving Runge-Kutta (SSP-RK) time-stepping schemes. Numerical results have shown that with such methods we are able to resolve solution structures that are only visible at much higher grid resolutions with lower-order schemes. The key challenge in applying such methods to ideal MHD is to control divergence errors in the magnetic field. We achieve this by augmenting the base scheme with a novel high-order constrained transport approach that updates the magnetic vector potential. The predicted magnetic field from the base scheme is replaced by a divergence-free magnetic field that is obtained from the curl of this magnetic potential. The non-conservative weakly hyperbolic system that the magnetic vector potential satisfies is solved using a version of FD-WENO developed for Hamilton-Jacobi equations. The resulting numerical method is endowed with several important properties: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as point values on the same mesh (i.e., there is no mesh staggering); (2) both the spatial and temporal orders of accuracy are fourth-order; (3) no spatial integration or multidimensional reconstructions are needed in any step; and (4) special limiters in the magnetic vector potential update are used to control unphysical oscillations in the magnetic field. Several 2D and 3D numerical examples are presented to verify the order of accuracy on smooth test problems and to show high-resolution on test problems that involve shocks.

39 pages, 9 figures, 4 tables

References in corpus (3)

• Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
• An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations
• High-order multiderivative time integrators for hyperbolic conservation laws

Cited by in corpus (15)

• Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes
• Provably Positive Discontinuous Galerkin Methods for Multidimensional Ideal Magnetohydrodynamics
• Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics
• A high-order positivity-preserving single-stage single-step method for the ideal magnetohydrodynamic equations
• A $4^{\mathrm{th}}$-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions
• A high-order weighted finite difference scheme with a multi-state approximate Riemann solver for divergence-free magnetohydrodynamic simulations
• High order asymptotic preserving finite difference WENO schemes with constrained transport for MHD equations in all sonic Mach numbers
• An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD
• HOW-MHD: A High-Order WENO-Based Magnetohydrodynamic Code with a High-Order Constrained Transport Algorithm for Astrophysical Applications
• High-Order Discontinuous Galerkin Finite Element Methods with Globally Divergence-Free Constrained Transport for Ideal MHD
• Mimetic Explicit Time Discretiztions
• A Kernel Based High Order "Explicit" Unconditionally Stable Constrained Transport Method for Ideal Magnetohydrodynamics
• A high-order finite difference WENO scheme for ideal magnetohydrodynamics on curvilinear meshes
• Explicit Time Mimetic Discretiztions Of Wave Equations
• Free-stream preserving finite difference schemes for ideal magnetohydrodynamics on curvilinear meshes