×

Well-balanced central scheme for the system of MHD equations with gravitational source term. (English) Zbl 1496.65138

Summary: A well-balanced second order finite volume central scheme for the magnetohydrodynamic (MHD) equations with gravitational source term is developed in this paper. The scheme is an unstaggered central scheme that evolves the numerical solution on a single grid and avoids solving Riemann problems at the cell interfaces using ghost staggered cells. A subtraction technique is used on the conservative variables with the support of a known steady state in order to manifest the well-balanced property of the scheme. The divergence-free constraint of the magnetic field is satisfied after applying the constrained transport method (CTM) for unstaggered central schemes at the end of each time-step by correcting the components of the magnetic field. The robustness of the proposed scheme is verified on a list of numerical test cases from the literature.

MSC:

65M08 Finite volume methods for initial value and initial-boundary value problems involving PDEs
76M12 Finite volume methods applied to problems in fluid mechanics
65M22 Numerical solution of discretized equations for initial value and initial-boundary value problems involving PDEs

References:

[1] P. Arminjon and R. Touma, Central finite volume methods with constrained transport diver-gence treatment for ideal MHD, Journal of Computational Physics, 204:737-759, 2005. · Zbl 1329.76196
[2] J. Berberich, P. Chandrashekar and C. Klingenberg, High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws, Computers & Fluids, 209, 2021 (https://doi.org/10.1016/j.compfluid.2021.104858). · Zbl 1521.76393 · doi:10.1016/j.compfluid.2021.104858
[3] T.J. Bogdan et al., Waves in the magnetized solar atmosphere II: Waves from localized sources in magnetic flux concentrations, Astrophys.J., 599:626-660, 2003.
[4] J.U. Brackbill and D.C. Barnes, The effect of nonzero ∇.B on the numerical solution of the magnetohydrodynamic equations, J. Comp. Phys., 201:261-285, 2004.
[5] J.U. Brackbill and D.C. Barnes, The effect of nonzero divb on the numerical solution of the magnetohydrodynamic equations, J. Comput. Phys., 35:426-430, 1980. · Zbl 0429.76079
[6] C.R. Evans and J.F. Hawley, Simulation of magnetohydrodynamic flows: A constrained transport method, Astrophys. J. Lett., 332:659, 1988.
[7] F.G. Fuchs, A.D. McMurry, S. Mishra, N.H. Risbro and K. Waagan, High order well-balanced finite volume schemes for simulating wave propagation in stratified magnetic atmospheres, Journal of Computational Physics, 229:4033-4058, 2010. · Zbl 1190.76153
[8] S.K. Godunov, The symmetric form of magnetohydrodynamics equation, Numer. Methods Mech. Contin. Media 1, 26, 1972.
[9] G.S. Jiang, D. Levy, C.T. Lin, S. Osher and E. Tadmor, High-resolution nonoscillatory cen-tral schemes with nonstaggered grids for hyperbolic conservation laws, SIAM Journal on Numerical Analysis, 35(6):2147-2168, 1998. · Zbl 0920.65053
[10] F. Kanbar, R. Touma and C. Klingenberg, Well-balanced central schemes for the one and two-dimensional Euler systems with gravity, Applied Numerical Mathematics, 156:608-626, 2020. · Zbl 1442.65199
[11] H. Nessyahu and E. Tadmor, Non-oscillatory central differencing for hyperbolic conserva-tion laws, Journal of Computational Physics, 87(2):408-463, 1990. · Zbl 0697.65068
[12] K.G. Powell, An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension), ICASE-Report 94-24 (NASA CR-194902) (NASA Langley Re-search Center, Hampton, VA 23681-0001, 8. April 1994).
[13] K.G. Powell, An approximate riemann solver for magnetohydrodynamics (that works in more than one space dimension), Technical report, 94-24, ICASE, Langley, VA, 1994.
[14] K.G. Powell, P.L. Roe, R.S. Myong, T. Gombosi and D. De Zeeuw, An upwind scheme for the magnetohydrodynamics, AIAA Paper 95-1704-CP, 1995.
[15] K.G. Powell, P.L. Roe, T.J. Linde, T.I. Gombosi and D.L. De Zeeuw, A solution adaptive upwind scheme for ideal MHD, J. Comput. Phys., 154(2):284-309, 1999. · Zbl 0952.76045
[16] C.S. Rosenthal et al., Waves in the magnetized solar atmosphere I: Basic processes and inter-network oscillations, Astrophys. J., 564:508-524, 2002.
[17] G. Toth, The divb=0 constraint in shock capturing magnetohydrodynamics codes, J. Com-put. Phys., 161:605-652, 2000. · Zbl 0980.76051
[18] R. Touma, Central unstaggered finite volume schemes for hyperbolic systems: Applications to unsteady shallow water equations, Applied Mathematics and Computation, 213(1):47-59, 2009. · Zbl 1165.76034
[19] R. Touma, Unstaggered central schemes with constrained transport treatment for ideal and shallow water magnetohydrodynamics, Applied Numerical Mathematics, 60:752-766, 2010. · Zbl 1425.76194
[20] R. Touma and S. Khankan, Well-balanced unstaggered central schemes for one and two-dimensional shallow water equation systems, Applied Mathematics and Computation, 218(10):5948-5960, 2012. · Zbl 1426.76433
[21] R. Touma and C. Klingenberg, Well-balanced central finite volume methods for the Ripa system, Applied Numerical Mathematics, 97:42-68, 2015. · Zbl 1329.76217
[22] R. Touma, U. Koley, and C. Klingenberg, Well-balanced unstaggered central schemes for the Euler equations with gravitation, SIAM J. Sci. Comput., 38(5):773-807, 2016. · Zbl 1346.76100
[23] Kailiang Wu and Chi-Wang Shu, Provably positive discontinuous Galerkin methods for multidimensional magnetohydrodynamics, SIAM J. Sci. Comput., 40:B1302-1329, 2018. · Zbl 1404.65184
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.