An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations
arXiv:1007.2606 · doi:10.1016/j.jcp.2011.02.009
Abstract
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [SIAM J. Sci. Comp. 28, 1766 (2006)], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J. Comp. Phys. 165, 126 (2000)]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The resulting scheme is applied to several numerical test cases.
46 pages, 12 figures
References in corpus (1)
Cited by in corpus (18)
- A Solution Accurate, Efficient and Stable Unsplit Staggered Mesh Scheme for Three Dimensional Magnetohydrodynamics
- Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics
- Improved Accuracy of High-Order WENO Finite Volume Methods on Cartesian Grids
- DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
- A high-order positivity-preserving single-stage single-step method for the ideal magnetohydrodynamic equations
- High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
- Systematic construction of upwind constrained transport schemes for MHD
- A -order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions
- A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations
- High order asymptotic preserving finite difference WENO schemes with constrained transport for MHD equations in all sonic Mach numbers
- WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
- Meshless Methods for Magnetohydrodynamics with Vector Potential
- High-Order Discontinuous Galerkin Finite Element Methods with Globally Divergence-Free Constrained Transport for Ideal MHD
- A high-order unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations based on the method of lines
- A high-order finite difference WENO scheme for ideal magnetohydrodynamics on curvilinear meshes
- A Kernel Based High Order "Explicit" Unconditionally Stable Constrained Transport Method for Ideal Magnetohydrodynamics
- Stable Finite Element Methods Preserving Exactly for MHD Models
- Free-stream preserving finite difference schemes for ideal magnetohydrodynamics on curvilinear meshes