Shock Structure for Hyperbolic Balance Laws with Evolutionary Constraints
Abstract
We study travelling shock profiles for dissipative quasilinear hyperbolic balance laws endowed with a convex entropy principle and subject to involutive differential constraints. For systems without constraints, the classical Boillat--Ruggeri result excludes continuous shock profiles whose speed exceeds the largest characteristic velocity at the equilibrium state ahead of the shock. We show that this bound remains unchanged for homogeneous constraints and for nonhomogeneous constraints whose source Jacobian has full row rank at equilibrium; a longitudinal two-species Gaussian ten-moment plasma with Landau collisions and Gauss's law provides a physical example of the latter case. For rank-deficient constraints, travelling-wave compatibility alone leads, after restriction to the joint conservation--constraint manifold, to a quadratic correction determined by a finite-dimensional Lyapunov equation. For 2×2 systems with one conservation law, one genuine balance law and one scalar involutive constraint, we prove that a regular noncharacteristic rank-deficient tail cannot sustain a nonzero correction, so that the original Boillat--Ruggeri bound is recovered. The remaining degenerate case is illustrated by an exact propagated two-field model in which the constraint selects the travelling speed and the smooth profile approaches equilibrium algebraically. The same model also admits entropy-admissible Lax composite profiles containing a sub-shock, while a completely smooth connection exists between the same equilibrium states. Numerical entropy evolutions further show that, depending on the smooth compatible initial datum, the first-order model may either remain smooth or dynamically develop a sub-shock.
Related
- Well-posedness for a class of n× n hyperbolic systems and shock stability for the viscous approximation
- Minimal null control time of some 1D linear hyperbolic balance laws with constant coefficients and properties of related kernel equations
- Sensitivity Calculus and its Numerical Implementation for Multi-D Hyperbolic Balance Laws
- Boundary stabilization of a moving interface in a 2 × 2 weakly hyperbolic system of balance laws: application to continuous sedimentation
- Boundary stabilization of a moving interface in a 2 x 2 weakly hyperbolic system of balance laws: application to continuous sedimentation
BibTeX
@misc{manganaro2026shock,
title = {{Shock Structure for Hyperbolic Balance Laws with Evolutionary Constraints}},
author = {Natale Manganaro and Tommaso Ruggeri},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.29533},
eprint = {2609.29533},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.29533},
}
RIS
TY - UNPB
TI - Shock Structure for Hyperbolic Balance Laws with Evolutionary Constraints
AU - Natale Manganaro
AU - Tommaso Ruggeri
PY - 2026
DA - 2026/09/01
PB - arXiv
DO - 10.48550/arxiv.2609.29533
UR - https://arxiv.org/abs/2609.29533
AB - We study travelling shock profiles for dissipative quasilinear hyperbolic balance laws endowed with a convex entropy principle and subject to involutive differential constraints. For systems without constraints, the classical Boillat--Ruggeri result excludes continuous shock profiles whose speed exceeds the largest characteristic velocity at the equilibrium state ahead of the shock. We show that this bound remains unchanged for homogeneous constraints and for nonhomogeneous constraints whose source Jacobian has full row rank at equilibrium; a longitudinal two-species Gaussian ten-moment plasma with Landau collisions and Gauss's law provides a physical example of the latter case. For rank-deficient constraints, travelling-wave compatibility alone leads, after restriction to the joint conservation--constraint manifold, to a quadratic correction determined by a finite-dimensional Lyapunov equation. For 2×2 systems with one conservation law, one genuine balance law and one scalar involutive constraint, we prove that a regular noncharacteristic rank-deficient tail cannot sustain a nonzero correction, so that the original Boillat--Ruggeri bound is recovered. The remaining degenerate case is illustrated by an exact propagated two-field model in which the constraint selects the travelling speed and the smooth profile approaches equilibrium algebraically. The same model also admits entropy-admissible Lax composite profiles containing a sub-shock, while a completely smooth connection exists between the same equilibrium states. Numerical entropy evolutions further show that, depending on the smooth compatible initial datum, the first-order model may either remain smooth or dynamically develop a sub-shock.
ER -