Analysis and structure-preserving discretization of the Heat-GLM system
Abstract
In this paper, we study a prototype hyperbolic system arising from the coupling of Generalized-Lagrangian-Multiplier (GLM) curl-cleaning with linear acoustics. We first show that, in the absence of dissipation, the equations can be rigorously derived from an underlying variational principle and that the resulting system is symmetric-hyperbolic. Further analysis shows that it not only conserves the total energy but also admits a set of invariants consisting of quadratic combinations of differential operators applied to the state variables. The addition of a stiff relaxation source term allows the system to be extended to dissipative dynamics, describing, for example, Cattaneo-type heat transfer processes that are asymptotically compatible with the Fourier law in the stiff relaxation limit; in this case, the total energy and invariants are dissipated accordingly. A new semi-implicit compatible numerical scheme on staggered grids is developed to solve this system while exactly preserving its key properties at the discrete level. In particular, we prove that the scheme is asymptotic-preserving in the Fourier limit, with a convergence rate depending on the initial data. We also show that the scheme conserves exactly the total energy as well as all the invariants in the absence of relaxation, and dissipates them consistently in its presence. These findings are demonstrated on a set of representative test cases.
Related
- Surface Stokes-Cahn-Hilliard System: Analysis and Structure-Preserving Discretization
- A dual-field structure-preserving mixed finite element discretization for incompressible Hall MHD equations
- A hyperbolic nonlinear Schrödinger system in multiple spatial dimensions: structure, convergence, and numerical discretization
- High-Order Structure-Preserving SBP Finite Difference Methods for the Vlasov-Maxwell System on Matrix-Free GPUs
- A Structure-Preserving Semi-Lagrangian based bi-fidelity method for the Vlasov--Poisson system with uncertainties
BibTeX
@misc{dhaouadi2026analysis,
title = {{Analysis and structure-preserving discretization of the Heat-GLM system}},
author = {Firas Dhaouadi and Laura Río-Martín and Michael Dumbser},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2607.22151},
eprint = {2607.22151},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.22151},
}
RIS
TY - UNPB
TI - Analysis and structure-preserving discretization of the Heat-GLM system
AU - Firas Dhaouadi
AU - Laura Río-Martín
AU - Michael Dumbser
PY - 2026
DA - 2026/10/27
PB - arXiv
DO - 10.48550/arxiv.2607.22151
UR - https://arxiv.org/abs/2607.22151
AB - In this paper, we study a prototype hyperbolic system arising from the coupling of Generalized-Lagrangian-Multiplier (GLM) curl-cleaning with linear acoustics. We first show that, in the absence of dissipation, the equations can be rigorously derived from an underlying variational principle and that the resulting system is symmetric-hyperbolic. Further analysis shows that it not only conserves the total energy but also admits a set of invariants consisting of quadratic combinations of differential operators applied to the state variables. The addition of a stiff relaxation source term allows the system to be extended to dissipative dynamics, describing, for example, Cattaneo-type heat transfer processes that are asymptotically compatible with the Fourier law in the stiff relaxation limit; in this case, the total energy and invariants are dissipated accordingly. A new semi-implicit compatible numerical scheme on staggered grids is developed to solve this system while exactly preserving its key properties at the discrete level. In particular, we prove that the scheme is asymptotic-preserving in the Fourier limit, with a convergence rate depending on the initial data. We also show that the scheme conserves exactly the total energy as well as all the invariants in the absence of relaxation, and dissipates them consistently in its presence. These findings are demonstrated on a set of representative test cases.
ER -