Sensitivity Calculus and its Numerical Implementation for Multi-D Hyperbolic Balance Laws
Abstract
We investigate the sensitivity of solutions to multi-dimensional scalar balance laws with respect to perturbations of the initial data analytically and numerically. Reliable first-order sensitivity information is essential in gradient-based optimization and inverse problems constrained by hyperbolic balance laws. In hyperbolic problems, such perturbations affect both the smooth components of the solution and the locations of shocks. Consequently, classical difference quotients of the solution operator generally fail to converge in L1, even in one space dimension. Although generalized tangent-vector techniques have been developed for one-dimensional problems, extending these ideas to multiple space dimensions is considerably more challenging because it requires a geometric description of hypersurfaces of discontinuity. We represent perturbations of hypersurfaces of discontinuity by normal displacements and account for the induced variation of the normal direction. This representation yields evolution equations for the components of a generalized tangent vector. Based on this calculus, we develop a numerical method for computing first-order variations with respect to the initial data. Numerical experiments in two space dimensions confirm the expected first-order accuracy of the resulting approximation.
Related
- A Study of King-Type Operators Connecting Wavelet and Quantum Calculus: Numerical and Graphical Analysis
- Minimal null control time of some 1D linear hyperbolic balance laws with constant coefficients and properties of related kernel equations
- Newton-type inequalities for the class of modified hyperbolic p-convex functions via Riemann-Liouville fractional integrals with numerical justifications
- 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{dreßen2026sensitivity,
title = {{Sensitivity Calculus and its Numerical Implementation for Multi-D Hyperbolic Balance Laws}},
author = {Olivia Dreßen and Michael Herty and Adrian Kolb and Siegfried Müller},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.20770},
eprint = {2609.20770},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.20770},
}
RIS
TY - UNPB
TI - Sensitivity Calculus and its Numerical Implementation for Multi-D Hyperbolic Balance Laws
AU - Olivia Dreßen
AU - Michael Herty
AU - Adrian Kolb
AU - Siegfried Müller
PY - 2026
DA - 2026/09/17
PB - arXiv
DO - 10.48550/arxiv.2609.20770
UR - https://arxiv.org/abs/2609.20770
AB - We investigate the sensitivity of solutions to multi-dimensional scalar balance laws with respect to perturbations of the initial data analytically and numerically. Reliable first-order sensitivity information is essential in gradient-based optimization and inverse problems constrained by hyperbolic balance laws. In hyperbolic problems, such perturbations affect both the smooth components of the solution and the locations of shocks. Consequently, classical difference quotients of the solution operator generally fail to converge in L^1, even in one space dimension. Although generalized tangent-vector techniques have been developed for one-dimensional problems, extending these ideas to multiple space dimensions is considerably more challenging because it requires a geometric description of hypersurfaces of discontinuity. We represent perturbations of hypersurfaces of discontinuity by normal displacements and account for the induced variation of the normal direction. This representation yields evolution equations for the components of a generalized tangent vector. Based on this calculus, we develop a numerical method for computing first-order variations with respect to the initial data. Numerical experiments in two space dimensions confirm the expected first-order accuracy of the resulting approximation.
ER -