Modal Discontinuous Galerkin Simulations for Grad's 13 Moment Equations: Application to Riemann Problem in Continuum-Rarefied Flow Regime

Satyvir Singh,Hang Song,Manuel Torrilhon

DOI: https://doi.org/10.1080/23324309.2024.2342947

IF: 0.348

2024-04-24

Journal of Computational and Theoretical Transport

Abstract:For real-world engineering applications, there is a great deal of interest in developing effective models and simulations of continuum-rarefied gas flows. In this study, the numerical simulations of Grad's 13 (G13) moment equations with application to the Riemann problem in a wide range of continuum-rarefied flow regimes are presented. This work emphasizes numerical robustness and wave phenomena in the G13 system to provide a building block for regularized 13 moment systems, and high-order Grad's models. For this purpose, a high-order modal discontinuous Galerkin solver is developed for solving one-dimensional G13 moment equations. For spatial discretization, hierarchical modal basis functions premised on orthogonal-scaled Legendre polynomials are used. The proposed approach reduces the G13 systems into a framework of ordinary differential equations in time, which are addressed by an explicit third-order SSP Runge-Kutta algorithm. Three Riemann test cases, including the shock tube, two shock waves and two rarefaction waves, are examined in continuum-rarefied flow regimes. In the G13 system, the arising characteristic waves and dissipation phenomena are investigated in depth. The numerical results demonstrate that every Riemann problem does not have a solution for the G13 system because of loss of hyperbolicity of the system.

mathematics, applied,physics, mathematical

What problem does this paper attempt to address?

The paper primarily aims to address the problem of effectively modeling and simulating gas flow in the continuum-rarefied flow regime. Specifically, the research focuses on the application of Grad's 13-moment equations (G13) in handling the Riemann problem, particularly concerning flow phenomena in the one-dimensional continuum-rarefied flow regime.

To achieve this goal, the researchers developed a high-order modal discontinuous Galerkin solver to solve the one-dimensional G13 moment equations. By using hierarchical modal basis functions based on orthogonally scaled Legendre polynomials for spatial discretization and employing an explicit third-order strong stability-preserving (SSP) Runge-Kutta scheme for time integration, the G13 equations are transformed into a time framework of ordinary differential equations.

The effectiveness and robustness of this method were validated through the study of three Riemann test cases, including a shock tube, two shock waves, and two rarefaction waves. These tests revealed the characteristic waves and dissipation phenomena present in the G13 system. It is particularly noteworthy that not all Riemann problems have solutions in the G13 system, due to the loss of hyperbolicity in the system.

In summary, this paper aims to investigate the behavior of the G13 moment equations in handling the Riemann problem through numerical simulations, especially in the continuum-rarefied flow regime, providing a theoretical foundation for understanding wave phenomena in such flows.