John M. Hong

Reyna Quita, Yu-Shuo Chen, Hsin-Yi Lee Alex C. Hu et al.

In this paper, a modified version of conservative Physics-informed Neural Networks (cPINN for short) is provided to construct the weak solutions of Riemann problem for the hyperbolic scalar conservation laws in non-conservative form. To demonstrate the results, we use the model of generalized Buckley-Leverett equation (GBL equation for short) with discontinuous porosity in porous media. By inventing a new unknown, the GBL equation is transformed into a two-by-two resonant hyperbolic conservation laws in conservative form. The modified method of cPINN is invented to overcome the difficulties due to the discontinuity of the porosity and the appearance of the critical states (near vacuum) in the Riemann data. We experiment with our idea by using a deep learning algorithm to solve the GBL equation in both conservative and non-conservative forms, as well as the cases of critical and non-critical states. This method provides a combination of two different neural networks and corresponding loss functions, one is for the two-by-two resonant hyperbolic system, and the other is for the scalar conservation law with a discontinuous perturbation term in the non-convex flux. The technique of re-scaling to the unknowns is adopted to avoid the oscillation of the Riemann solutions in the cases of critical Riemann data. The solutions constructed by the modified cPINN match the exact solutions constructed by the theoretical analysis for hyperbolic conservation laws. In addition, the solutions are identical in both conservative and non-conservative cases. Finally, we compare the performance of the modified cPINN with numerical method called WENO5. Whereas WENO5 struggles with the highly oscillation of approximate solutions for the Riemann problems of GBL equation in non-conservative form, cPINN works admirably.

John M. Hong, Philippe G. LeFloch

We introduce a generalization of Glimm's random choice method, which provides us with an approximation of entropy solutions to quasilinear hyperbolic system of balance laws. The flux-function and the source term of the equations may depend on the unknown as well as on the time and space variables. The method is based on local approximate solutions of the generalized Riemann problem, which form building blocks in our scheme and allow us to take into account naturally the effects of the flux and source terms. To establish the nonlinear stability of these approximations, we investigate nonlinear interactions between generalized wave patterns. This analysis leads us to a global existence result for quasilinear hyperbolic systems with source-term, and applies, for instance, to the compressible Euler equations in general geometries and to hyperbolic systems posed on a Lorentzian manifold.