Abstract:In this paper, we are interested in constructing a scheme solving compressible Navier--Stokes equations, with desired properties including high order spatial accuracy, conservation, and positivity-preserving of density and internal energy under a standard hyperbolic type CFL constraint on the time step size, e.g., $\Delta t=\mathcal O(\Delta x)$. Strang splitting is used to approximate convection and diffusion operators separately. For the convection part, i.e., the compressible Euler equation, the high order accurate postivity-preserving Runge--Kutta discontinuous Galerkin method can be used. For the diffusion part, the equation of internal energy instead of the total energy is considered, and a first order semi-implicit time discretization is used for the ease of achieving positivity. A suitable interior penalty discontinuous Galerkin method for the stress tensor can ensure the conservation of momentum and total energy for any high order polynomial basis. In particular, positivity can be proven with $\Delta t=\mathcal{O}(\Delta x)$ if the Laplacian operator of internal energy is approximated by the $\mathbb{Q}^k$ spectral element method with $k=1,2,3$. So the full scheme with $\mathbb{Q}^k$ ($k=1,2,3$) basis is conservative and positivity-preserving with $\Delta t=\mathcal{O}(\Delta x)$, which is robust for demanding problems such as solutions with low density and low pressure induced by high-speed shock diffraction. Even though the full scheme is only first order accurate in time, numerical tests indicate that higher order polynomial basis produces much better numerical solutions, e.g., better resolution for capturing the roll-ups during shock reflection.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to construct a numerical scheme for solving the compressible Navier - Stokes equations. This scheme has high - order spatial accuracy, conserves quantities, and maintains the positivity of density and internal energy under the standard hyperbolic CFL condition (such as \(\Delta t = O(\Delta x)\)). Specifically, the author focuses on how to design a numerical method that can both maintain the positivity of physical quantities (i.e., density and internal energy are always positive) and conserve mass, momentum, and total energy, while also achieving high spatial accuracy. This is especially important for dealing with complex fluid - dynamic problems involving low - density and low - pressure regions, because these problems are the main sources of instability in high - order methods.
To achieve this goal, the author adopts the Strang splitting method to approximate the convection and diffusion operators separately, and uses the high - order accurate positivity - preserving Runge - Kutta discontinuous Galerkin method to handle the convection part, i.e., the compressible Euler equations. For the diffusion part, the internal energy equation rather than the total energy equation is considered, and a first - order semi - implicit time discretization is adopted to simplify the implementation of positivity. In addition, by choosing an appropriate interior - penalty discontinuous Galerkin method to handle the stress tensor, the conservation of momentum and total energy under any high - order polynomial basis can be ensured. In particular, if the \(Q_k\) spectral - element method (\(k = 1,2,3\)) is used to approximate the Laplacian operator of internal energy, then positivity can be proven under the condition \(\Delta t=O(\Delta x)\).
In summary, this research aims to propose a numerical scheme that is both conservative and positivity - preserving, and can effectively handle complex fluid - dynamic problems, especially those involving low - density and low - pressure regions.