Integral Invariants of Vasudeva Murthy's Relaxation Systems: Analysis and Numerical Validation
Abstract
Vasudeva Murthy's relaxation approach [A.S. Vasudeva Murthy, J. Comput. Appl. Math., 203(2), pp. 437-443, 2007], originally proposed for the Jin-Xin relaxation model, provides an alternative formulation with invariant properties that is consistent and retains the semilinear structure incomparison to the standard one. In this work, Vasudeva Murthy's relaxation approach for various relaxation systems are proposed such as the shallow water equations, the Broadwell model, the Euler equations with heat transfer and two-dimensional Jin-Xin model. For proposed relaxation models, the associated integral invariants are rigorously established at the theoretical level. The main advantage of the integral invariant is that it provides a conserved quantity for the relaxation system by incorporating the coupled contributions of the solution variables in vector form. To validate the analytical results, numerical simulations are carried out for each model using three second-order numerical schemes: CS-EBT2, a semi-implicit second-order central finite-volume scheme for hyperbolic systems with relaxation source terms [S. Sahu, E. Macca, and R. Samala, J. Comput. Phys., 563, 115100, 2026]; UCS2, a finite-volume central relaxation-type scheme [S. F. Liotta, V. Romano, and G. Russo, SIAM J. Numer. Anal., 38(4), 1337-1356, 2000]; and IMEX-RK2, a second-order Implicit-Explicit Runge-Kutta scheme [Pareschi and Russo, J. Sci. Comput., 25, 129-155, 2005]. Numerical results, compared to exact or finely resolved reference solutions, confirm that the models preserve integral invariants, remain stable under CFL restrictions, and exhibit robust and accurate behavior across all benchmark systems tested.
Related
- Identification of generic polygonal domains by integral-geometric invariants
- Grothendieck ring of varieties and at-the-characteristic numerical invariants
- Numerical invariants for three-generated ideals
- A Unified System of Tensor Invariants and Their Integral Reduction Theorems
- Numerical Invariants and Parameter Recovery for Quasi-Homogeneous Submodules [zk−wℓ]
BibTeX
@misc{sahu2026integral,
title = {{Integral Invariants of Vasudeva Murthy's Relaxation Systems: Analysis and Numerical Validation}},
author = {Sudipta Sahu and Rathan Samala},
howpublished = {arXiv},
year = {2026},
month = aug,
doi = {10.48550/arxiv.2608.30770},
eprint = {2608.30770},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.30770},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Integral Invariants of Vasudeva Murthy's Relaxation Systems: Analysis and Numerical Validation
AU - Sudipta Sahu
AU - Rathan Samala
PY - 2026
DA - 2026/08/31
PB - arXiv
DO - 10.48550/arxiv.2608.30770
UR - https://arxiv.org/abs/2608.30770
AB - Vasudeva Murthy's relaxation approach [A.S. Vasudeva Murthy, J. Comput. Appl. Math., 203(2), pp. 437-443, 2007], originally proposed for the Jin-Xin relaxation model, provides an alternative formulation with invariant properties that is consistent and retains the semilinear structure incomparison to the standard one. In this work, Vasudeva Murthy's relaxation approach for various relaxation systems are proposed such as the shallow water equations, the Broadwell model, the Euler equations with heat transfer and two-dimensional Jin-Xin model. For proposed relaxation models, the associated integral invariants are rigorously established at the theoretical level. The main advantage of the integral invariant is that it provides a conserved quantity for the relaxation system by incorporating the coupled contributions of the solution variables in vector form. To validate the analytical results, numerical simulations are carried out for each model using three second-order numerical schemes: CS-EBT2, a semi-implicit second-order central finite-volume scheme for hyperbolic systems with relaxation source terms [S. Sahu, E. Macca, and R. Samala, J. Comput. Phys., 563, 115100, 2026]; UCS2, a finite-volume central relaxation-type scheme [S. F. Liotta, V. Romano, and G. Russo, SIAM J. Numer. Anal., 38(4), 1337-1356, 2000]; and IMEX-RK2, a second-order Implicit-Explicit Runge-Kutta scheme [Pareschi and Russo, J. Sci. Comput., 25, 129-155, 2005]. Numerical results, compared to exact or finely resolved reference solutions, confirm that the models preserve integral invariants, remain stable under CFL restrictions, and exhibit robust and accurate behavior across all benchmark systems tested.
ER -