Communications on Applied Mathematics and Computation ›› 2024, Vol. 6 ›› Issue (3): 1832-1859.doi: 10.1007/s42967-023-00334-1

• ORIGINAL PAPERS • Previous Articles     Next Articles

Machine Learning Approaches for the Solution of the Riemann Problem in Fluid Dynamics: a Case Study

Vitaly Gyrya1, Mikhail Shashkov2, Alexei Skurikhin3, Svetlana Tokareva1   

  1. 1 Los Alamos National Laboratory, Applied Mathematics and Plasma Physics (T-5), Los Alamos, USA;
    2 Los Alamos National Laboratory, Continuum Models and Numerical Methods (XCP-4), Los Alamos, USA;
    3 Los Alamos National Laboratory, Space Data Science and Systems (ISR-3), Los Alamos, USA
  • Received:2022-12-26 Revised:2022-12-26 Accepted:2023-09-11 Published:2024-12-20
  • Contact: Svetlana Tokareva,tokareva@lanl.gov;Vitaly Gyrya,vitaliy_gyrya@lanl.gov;Mikhail Shashkov,shashkov@lanl.gov;Alexei Skurikhin,alexei@lanl.gov E-mail:tokareva@lanl.gov;vitaliy_gyrya@lanl.gov;shashkov@lanl.gov;alexei@lanl.gov
  • Supported by:
    This work was performed under the auspices of the National Nuclear Security Administration of the US Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396. The authors gratefully acknowledge the support of the US Department of Energy National Nuclear Security Administration Advanced Simulation and Computing Program. The Los Alamos unlimited release number is LA-UR-19-32257. The authors would like to thank J. Albright for numerous discussions and his valuable comments.

Abstract: We present our results by using a machine learning (ML) approach for the solution of the Riemann problem for the Euler equations of fluid dynamics. The Riemann problem is an initial-value problem with piecewise-constant initial data and it represents a mathematical model of the shock tube. The solution of the Riemann problem is the building block for many numerical algorithms in computational fluid dynamics, such as finite-volume or discontinuous Galerkin methods. Therefore, a fast and accurate approximation of the solution of the Riemann problem and construction of the associated numerical fluxes is of crucial importance. The exact solution of the shock tube problem is fully described by the intermediate pressure and mathematically reduces to finding a solution of a nonlinear equation. Prior to delving into the complexities of ML for the Riemann problem, we consider a much simpler formulation, yet very informative, problem of learning roots of quadratic equations based on their coefficients. We compare two approaches: (i) Gaussian process (GP) regressions, and (ii) neural network (NN) approximations. Among these approaches, NNs prove to be more robust and efficient, although GP can be appreciably more accurate (about 30%). We then use our experience with the quadratic equation to apply the GP and NN approaches to learn the exact solution of the Riemann problem from the initial data or coefficients of the gas equation of state (EOS). We compare GP and NN approximations in both regression and classification analysis and discuss the potential benefits and drawbacks of the ML approach.

Key words: Machine learning (ML), Neural network (NN), Gaussian process (GP), Riemann problem, Numerical fluxes, Finite-volume method

CLC Number: