Provably Positivity-Preserving Constrained Transport (PPCT) Second-Order Scheme for Ideal Magnetohydrodynamics
This paper proposes and analyzes a robust and efficient second-order positivity-preserving constrained transport (PPCT) scheme for ideal magnetohydrodynamics (MHD) on non-staggered Cartesian meshes. The PPCT scheme ensures two critical physical constraints: a globally discrete divergence-free (DDF)...
Gespeichert in:
Hauptverfasser: | , |
---|---|
Format: | Artikel |
Sprache: | eng |
Schlagworte: | |
Online-Zugang: | Volltext bestellen |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Zusammenfassung: | This paper proposes and analyzes a robust and efficient second-order
positivity-preserving constrained transport (PPCT) scheme for ideal
magnetohydrodynamics (MHD) on non-staggered Cartesian meshes. The PPCT scheme
ensures two critical physical constraints: a globally discrete divergence-free
(DDF) condition on the magnetic field and the positivity of density and
pressure. The method is inspired by a novel splitting technique from [T.A. Dao,
M. Nazarov and I. Tomas, J. Comput. Phys., 508:113009, 2024], which divides the
MHD system into an Euler subsystem with steady magnetic fields and a magnetic
subsystem with steady density and internal energy. To achieve these
structure-preserving properties, the PPCT scheme combines a
positivity-preserving (PP) finite volume method for the Euler subsystem with a
finite difference constrained transport (CT) method for the magnetic subsystem
via Strang splitting. The finite volume method employs a new PP limiter that
retains second-order accuracy and enforces the positivity of density and
pressure, with rigorous proof provided using the geometric quasilinearization
(GQL) approach [K. Wu and C.-W. Shu, SIAM Review, 65:1031-1073, 2023]. For the
magnetic subsystem, we develop an implicit finite difference CT method that
conserves energy and maintains a globally DDF constraint. This nonlinear system
is efficiently solved to machine precision using an iterative algorithm. Since
the CT method is unconditionally energy-stable and conserves steady density and
internal energy, the PPCT scheme requires only a mild CFL condition for the
finite volume method to ensure stability and the PP property. While the focus
is on 2D cases for clarity, the extension to 3D is discussed. Several
challenging numerical experiments, including highly magnetized MHD jets with
high Mach numbers, validate the PPCT scheme's accuracy, robustness, and high
resolution. |
---|---|
DOI: | 10.48550/arxiv.2410.05173 |