A Positivity-Preserving Semi-Implicit Discontinuous Galerkin Scheme for Solving Extended Magnetohydrodynamics Equations
Document Type
Article
Publication Date
12-2014
Department
Department of Mathematical Sciences
Abstract
A positivity-preserving discontinuous Galerkin (DG) scheme [42] is used to solve the Extended Magnetohydrodynamics (XMHD) model, which is a two-fluid model expressed with a center-of-mass formulation. We prove that DG scheme with a positivity-preserving limiter is stable for the system governed by the XMHD model or the resistive MHD model. We use the relaxation system formulation [28] for describing the XMHD model, and solve the equations using a split level implicit-explicit time advance scheme, stepping over the time step constraint imposed by the stiff source terms. The magnetic field is represented in an exact locally divergence-free form of DG [23], which greatly improves the accuracy and stability of MHD simulations. As presently constructed, the method is able to handle a wide range of density variations, solve XMHD model on MHD time scales, and provide greatly improved accuracy over a Finite Volume implementation of the same model.
Publication Title
Journal of Computational Physics
Recommended Citation
Zhao, X.,
Yang, Y.,
&
Seyler, C.
(2014).
A Positivity-Preserving Semi-Implicit Discontinuous Galerkin Scheme for Solving Extended Magnetohydrodynamics Equations.
Journal of Computational Physics,
278, 400-415.
http://doi.org/10.1016/j.jcp.2014.08.044
Retrieved from: https://digitalcommons.mtu.edu/michigantech-p/6664
Publisher's Statement
© 2014 The Authors.