A Provable Positivity-Preserving Local Discontinuous Galerkin Method for the Viscous and Resistive MHD Equations

Abstract

In this paper, we construct a high-order discontinuous Galerkin (DG) method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics (VRMHD). To control the divergence error in the magnetic field, both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD. Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes, respectively, showing that the scheme can maintain the positivity-preserving (PP) property under some CFL conditions when combined with the strong-stability-preserving time discretization. Then, general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes. Numerical tests demonstrate the effectiveness of the proposed schemes.

Keywords

Viscous and resistive MHD equations / Positivity-preserving / Discontinuous Galerkin (DG) method / High order accuracy

References

[1.] Balsara DS. Second-order-accurate schemes for magnetohydrodynamics with divergence-free reconstruction. Astrophys. J. Suppl. Ser., 2004, 151(1): 149

[2.] Balsara DS. Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics. J. Comput. Phys., 2012, 231(22): 7504-7517

[3.] Balsara DS, Spicer DS. A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. J. Comput. Phys., 1999, 149(2): 270-292

[4.] Barth T. Arnold DN, Bochev PB, Lehoucq RB, Nicolaides RA, Shashkov M. On the role of involutions in the discontinuous Galerkin discretization of Maxwell and magnetohydrodynamic systems. Compatible Spatial Discretizations, 2006 New York, NY Springer 69-88

[5.] Barth T. Kröner D, Ohlberger M, Rohde C. Numerical methods for gasdynamic systems on unstructured meshes. An Introduction to Recent Developments in Theory and Numerics for Conservation Laws, 1999 Berlin, Heidelberg Springer 195-285

[6.] Brackbill JU, Barnes DC. The effect of nonzero ∇ · B \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nabla \cdot {B}$$\end{document} on the numerical solution of the magnetohydrodynamic equations. J. Comput. Phys., 1980, 35(3): 426-430

[7.] Chandrashekar P, Gallego-Valencia JP, Klingenberg C. Klingenberg C, Westdickenberg M. A Runge-Kutta discontinuous Galerkin scheme for the ideal magnetohydrodynamical model. Theory, Numerics and Applications of Hyperbolic Problems I. HYP 2016, 2016 Cham Springer 335-344

[8.] Chandrashekar P, Klingenberg C. Entropy stable finite volume scheme for ideal compressible MHD on 2-D Cartesian meshes. SIAM J. Numer. Anal., 2016, 54(2): 1313-1340

[9.] Cheng Y, Li F, Qiu J, Xu LW. Positivity-preserving DG and central DG methods for ideal MHD equations. J. Comput. Phys., 2013, 238: 255-280

[10.] Christlieb AJ, Feng X, Jiang Y, Tang Q. A high-order finite difference WENO scheme for ideal magnetohydrodynamics on curvilinear meshes. SIAM J. Sci. Comput., 2018, 40(4): A2631-A2666

[11.] Christlieb AJ, Liu Y, Tang Q, Xu ZF. Positivity-preserving finite difference weighted ENO schemes with constrained transport for ideal magnetohydrodynamic equations. SIAM J. Sci. Comput., 2015, 37(4): A1825-A1845

[12.] Christlieb AJ, Rossmanith JA, Tang Q. Finite difference weighted essentially non-oscillatory schemes with constrained transport for ideal magnetohydrodynamics. J. Comput. Phys., 2014, 268: 302-325

[13.] Dai W, Woodward PR. A simple finite difference scheme for multidimensional magnetohydrodynamical equations. J. Comput. Phys., 1998, 142(2): 331-369

[14.] Dedner A, Kemm F, Króner D, Munz CD, Schnitzer T, Wesenberg M. Hyperbolic divergence cleaning for the MHD equations. J. Comput. Phys., 2002, 175(2): 645-673

[15.] Evans CR, Hawley JF. Simulation of magnetohydrodynamic flows: a constrained transport method. Astrophys. J., 1988, 332: 659-677

[16.] Godunov SK. Symmetric form of the equations of magnetohydrodynamics. Numer. Methods Mech. Cont. Medium., 1972, 1: 26-34

[17.] Hu XY, Adams NA, Shu C-W. Positivity-preserving method for high-order conservative schemes solving compressible Euler equations. J. Comput. Phys., 2013, 242: 169-180

[18.] Li F, Shu C-W. Locally divergence-free discontinuous Galerkin methods for MHD equations. J. Sci. Comput., 2005, 22(1): 413-442

[19.] Li F, Xu LW. Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations. J. Comput. Phys., 2012, 231(6): 2655-2675

[20.] Li F, Xu LW, Yakovlev S. Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field. J. Comput. Phys., 2011, 230(12): 4828-4847

[21.] Liang C, Xu Z. Parametrized maximum principle preserving flux limiters for high order schemes solving multi-dimensional scalar hyperbolic conservation laws. J. Sci. Comput., 2014, 58(1): 41-60

[22.] Liu Y, Shu C-W, Zhang M. Entropy stable high order discontinuous Galerkin methods for ideal compressible MHD on structured meshes. J. Comput. Phys., 2018, 354: 163-178

[23.] Londrillo P, Del Zanna L. High-order upwind schemes for multidimensional magnetohydrodynamics. Astrophys. J., 2000, 530(1): 508

[24.] Powell KG. Hussaini MY, van Leer B, Van Rosendale J. An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension). Upwind and High-Resolution Schemes, 1997 Berlin, Heidelberg Springer

[25.] Powell, K.G., Roe, P., Myong, R., Gombosi, T.: An upwind scheme for magnetohydrodynamics. In: AIAA 12thComputational Fluid Dynamics Conference, San Diego, CA, June 19–22 1995, AIAA-95-1704-CP, AIAA (1995)

[26.] Qiu J, Shu C-W. Runge-Kutta discontinuous Galerkin method using WENO limiters. SIAM J. Sci. Comput., 2005, 26(3): 907-929

[27.] Rossmanith JA. An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows. SIAM J. Sci. Comput., 2006, 28(5): 1766-1797

[28.] Ryu D, Miniati F, Jones TW, Frank A. A divergence-free upwind code for multidimensional magnetohydrodynamic flows. Astrophys. J., 1998, 509(1): 244-255

[29.] Torrilhon M. Locally divergence-preserving upwind finite volume schemes for magnetohydrodynamic equations. SIAM J. Sci. Comput., 2005, 26(4): 1166-1191

[30.] Tóth G. The ∇ · B \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nabla \cdot \textbf{B}$$\end{document}= 0 constraint in shock-capturing magnetohydrodynamics codes. J. Comput. Phys., 2000, 161(2): 605-652

[31.] Warburton TC, Karniadakis GE. A discontinuous Galerkin method for the viscous MHD equations. J. Comput. Phys., 1999, 152(2): 608-641

[32.] Wu KL. Design of provably physical-constraint-preserving methods for general relativistic hydrodynamics. Phys. Rev. D, 2017, 95(10)

[33.] Wu KL. Positivity-preserving analysis of numerical schemes for ideal magnetohydrodynamics. SIAM J. Numer. Anal., 2018, 56(4): 2124-2147

[34.] Wu, K.L., Jiang, H., Shu, C.-W.: Provably positive central DG schemes via geometric quasilinearization for ideal MHD equations. arXiv:2203.14853 (2022)

[35.] Wu KL, Shu C-W. A provably positive discontinuous Galerkin method for multidimensional ideal magnetohydrodynamics. SIAM J. Sci. Comput., 2018, 40(5): B1302-B1329

[36.] Wu KL, Shu C-W. Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes. Numer. Math., 2019, 142(4): 995-1047

[37.] Wu, K.L., Shu, C.-W.: Geometric quasilinearization framework for analysis and design of bound-preserving schemes. arXiv:2111.04722 (2021)

[38.] Wu KL, Shu C-W. Provably physical-constraint-preserving discontinuous Galerkin methods for multidimensional relativistic MHD equations. Numer. Math., 2021, 148(3): 699-741

[39.] Wu KL, Tang H. Admissible states and physical-constraints-preserving schemes for relativistic magnetohydrodynamic equations. Math. Models Methods Appl. Sci., 2017, 27(10): 1871-1928

[40.] Xing Y, Zhang X, Shu C-W. Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations. Adv. Water Resour., 2010, 33(12): 1476-1493

[41.] Xu Z. Parametrized maximum principle preserving flux limiters for high order schemes solving hyperbolic conservation laws: one-dimensional scalar problem. Math. Comput., 2014, 83(289): 2213-2238

[42.] Xu Z, Zhang X. Abgrall R, Shu C-W. Bound-preserving high-order schemes. Handbook of Numerical Analysis, 2017 North-Holland, Amsterdam Elsevier 81-102

[43.] Yakovlev S, Xu LW, Li F. Locally divergence-free central discontinuous Galerkin methods for ideal MHD equations. J. Comput. Sci., 2013, 4(1/2): 80-91

[44.] Yu YZ, Jiang Y, Zhang M. Free-stream preserving finite difference schemes for ideal magnetohydrodynamics on curvilinear meshes. J. Sci. Comput., 2020, 82(1): 1-26

[45.] Zhang X. On positivity-preserving high order discontinuous Galerkin schemes for compressible Navier-Stokes equations. J. Comput. Phys., 2017, 328: 301-343

[46.] Zhang X, Shu C-W. On maximum-principle-satisfying high order schemes for scalar conservation laws. J. Comput. Phys., 2010, 229(9): 3091-3120

[47.] Zhang X, Shu C-W. On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes. J. Comput. Phys., 2010, 229(23): 8918-8934

[48.] Zhang X, Shu C-W. Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: survey and new developments. Proc. R. Soc. A Math. Phys. Eng. Sci., 2011, 467(2134): 2752-2776

[49.] Zhao J, Tang H. Runge-Kutta discontinuous Galerkin methods for the special relativistic magnetohydrodynamics. J. Comput. Phys., 2017, 343: 33-72

Funding

National Natural Science Foundation of China(12126604); Cyrus Tang Foundation