Fully Discrete Multi-Entropy Stability of High-Order Schemes for Compressible MHD: A Weak-to-Strong Framework
Abstract
We establish a fully discrete weak-to-strong (W2S) multi-entropy stability theory for arbitrarily high-order finite-volume and discontinuous Galerkin (DG) approximations of the ideal compressible magnetohydrodynamics (MHD) equations on general polytopal meshes, whereby a single numerical update simultaneously satisfies discrete entropy inequalities for any prescribed finite family of convex Harten entropy pairs. Because physical entropies are defined only for positive density and pressure, the central analytical difficulty lies in reconciling discrete entropy stability with positivity preservation, the magnetic divergence constraint, and the nonconservative Godunov--Powell coupling. Building upon the provably positivity-preserving MHD framework of Wu and Shu, we develop a weak multi-entropy analysis of the common cell-average evolution of finite-volume and DG methods. Through convex decomposition and relative entropy, we establish two sufficient criteria for this weak stability general polytopal meshes. Compatible magnetic corrections and locally divergence-free approximations cooperatively offset the nonconservative Godunov--Powell coupling and provide the discrete cancellations essential to both positivity-preservation and entropy stability. A W2S lifting then unifies positivity and entropy limiting into a single cellwise scaling limiter, achieving strong multi-entropy stability while strictly preserving cell averages and the locally divergence-free magnetic field. High-order temporal accuracy follows from variable-step strong-stability-preserving multistep methods. Numerical experiments confirm the theoretical findings and demonstrate computational robustness. This framework provides a rigorous, fully discrete foundation for unifying physical admissibility, magnetic divergence control, and multi-entropy stability in high-order MHD approximations.
Related
- Limiter-based fully-discrete entropy stable explicit DG schemes for ideal MHD equations
- Isothermal walls in the compressible Navier-Stokes equations: affine entropy, ballistic energy, and entropy-stable discrete formulations
- Simultaneous energy and entropy stability of volume-penalisation immersed-boundary methods for the compressible Navier--Stokes equations
- A Novel Intuitionistic Fuzzy Entropy with Comprehensive Ranking-Stability Analysis: An Entropy-Weighted IF-TOPSIS Framework
- Skew-Symmetric Formulations and Energy-Stable Schemes for Compressible Reactive Flow Equations
BibTeX
@misc{cao2026fully,
title = {{Fully Discrete Multi-Entropy Stability of High-Order Schemes for Compressible MHD: A Weak-to-Strong Framework}},
author = {Huihui Cao and Kailiang Wu and Caiyou Yuan},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.26791},
eprint = {2609.26791},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.26791},
}
RIS
TY - UNPB
TI - Fully Discrete Multi-Entropy Stability of High-Order Schemes for Compressible MHD: A Weak-to-Strong Framework
AU - Huihui Cao
AU - Kailiang Wu
AU - Caiyou Yuan
PY - 2026
DA - 2026/09/22
PB - arXiv
DO - 10.48550/arxiv.2609.26791
UR - https://arxiv.org/abs/2609.26791
AB - We establish a fully discrete weak-to-strong (W2S) multi-entropy stability theory for arbitrarily high-order finite-volume and discontinuous Galerkin (DG) approximations of the ideal compressible magnetohydrodynamics (MHD) equations on general polytopal meshes, whereby a single numerical update simultaneously satisfies discrete entropy inequalities for any prescribed finite family of convex Harten entropy pairs. Because physical entropies are defined only for positive density and pressure, the central analytical difficulty lies in reconciling discrete entropy stability with positivity preservation, the magnetic divergence constraint, and the nonconservative Godunov--Powell coupling. Building upon the provably positivity-preserving MHD framework of Wu and Shu, we develop a weak multi-entropy analysis of the common cell-average evolution of finite-volume and DG methods. Through convex decomposition and relative entropy, we establish two sufficient criteria for this weak stability general polytopal meshes. Compatible magnetic corrections and locally divergence-free approximations cooperatively offset the nonconservative Godunov--Powell coupling and provide the discrete cancellations essential to both positivity-preservation and entropy stability. A W2S lifting then unifies positivity and entropy limiting into a single cellwise scaling limiter, achieving strong multi-entropy stability while strictly preserving cell averages and the locally divergence-free magnetic field. High-order temporal accuracy follows from variable-step strong-stability-preserving multistep methods. Numerical experiments confirm the theoretical findings and demonstrate computational robustness. This framework provides a rigorous, fully discrete foundation for unifying physical admissibility, magnetic divergence control, and multi-entropy stability in high-order MHD approximations.
ER -