We develop a new second-order unstaggered semidiscrete path-conservative central- upwind (PCCU) scheme for ideal and shallow water magnetohydrodynamics (MHD) equations. The new scheme possesses several important properties: it locally preserves the divergence-free constraint, it does not rely on any (approximate) Riemann problem solver, and it robustly produces high- resolution and nonoscillatory results. The derivation of the scheme is based on the Godunov-Powell nonconservative modifications of the studied MHD systems. The local divergence-free property is enforced by augmenting the modified systems with the evolution equations for the corresponding derivatives of the magnetic field components. These derivatives are then used to design a special piecewise linear reconstruction of the magnetic field, which guarantees a nonoscillatory nature of the resulting scheme. In addition, the proposed PCCU discretization accounts for the jump of the nonconservative product terms across cell interfaces, thereby ensuring stability. We test the proposed PCCU scheme on several benchmarks for both ideal and shallow water MHD systems. The obtained numerical results illustrate the performance of the new scheme, its robustness, and its ability not only to achieve high resolution, but also to preserve the positivity of computed quantities such as density, pressure, and water depth.
more »
« less
A divergence-free semi-implicit finite volume scheme for ideal, viscous, and resistive magnetohydrodynamics: A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics
- Award ID(s):
- 1713765
- PAR ID:
- 10108025
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Fluids
- Volume:
- 89
- Issue:
- 1-2
- ISSN:
- 0271-2091
- Page Range / eLocation ID:
- 16 to 42
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We construct a structure-preserving finite element method and time-stepping scheme for compressible barotropic magnetohydrodynamics both in the ideal and resistive cases, and in the presence of viscosity. The method is deduced from the geometric variational formulation of the equations. It preserves the balance laws governing the evolution of total energy and magnetic helicity, and preserves mass and the constraint $$\text {div}B = 0$$ to machine precision, both at the spatially and temporally discrete levels. In particular, conservation of energy and magnetic helicity hold at the discrete levels in the ideal case. It is observed that cross-helicity is well conserved in our simulation in the ideal case.more » « less
An official website of the United States government

