English

Solving Random Hyperbolic Conservation Laws Using Linear Programming

Numerical Analysis 2025-10-29 v3 Numerical Analysis

Abstract

A novel structure-preserving numerical method to solve random hyperbolic systems of conservation laws is presented. The method uses a concept of generalized, measure-valued solutions to random conservation laws. This yields a linear partial differential equation with respect to the Young measure and allows to compute the approximation based on linear programming problems. We analyze structure-preserving properties of the derived numerical method and discuss its advantages and disadvantages. We numerically demonstrate the approach on the one-dimensional Burgers and isentropic Euler equations and compare with stochastic collocation. In addition, we introduce a discontinuous-flux test in which different entropies used in the linear-program objective select different weak entropy solutions, and we report the corresponding changes in the moments and supports of the Young measure.

Keywords

Cite

@article{arxiv.2501.10104,
  title  = {Solving Random Hyperbolic Conservation Laws Using Linear Programming},
  author = {Shaoshuai Chu and Michael Herty and Maria Lukacova-Medvidova and Yizhou Zhou},
  journal= {arXiv preprint arXiv:2501.10104},
  year   = {2025}
}