A data-driven physics-informed finite-volume scheme for nonclassical undercompressive shocks

Abstract

We propose a data-driven physics-informed finite-volume scheme for the approximation of small-scale dependent shocks. Nonlinear hyperbolic conservation laws with non-convex fluxes allow nonclassical shock wave solutions. In this work, we consider the cubic scalar conservation law as representative of such systems. As standard numerical schemes fail to approximate nonclassical shocks, schemes with controlled dissipation and schemes with well-controlled dissipation have been introduced by LeFloch and Mohammadian and by Ernest and coworkers, respectively. Emphasis has been placed on matching the truncation error of the numerical scheme with physically relevant small-scale mechanisms. However, aforementioned schemes can introduce oscillations as well as excessive dissipation around shocks. In our approach, a convolutional neural network is used for an adaptive nonlinear flux reconstruction. Based on the local flow field, the network combines local interpolation polynomials with a regularization term to form the numerical flux. This allows to modify the discretization error by nonlinear terms. In a supervised learning task, the model is trained to predict the time evolution of exact solutions to Riemann problems. The model is physics-informed as it respects the underlying conservation law. Numerical experiments for the cubic scalar conservation law show that the resulting method is able to approximate nonclassical shocks very well. The adaptive reconstruction suppresses oscillations and enables sharp shock capturing. Generalization to unseen shock configurations, smooth initial value problems, and shock interactions is robust and shows very good results.