Back to Search Start Over

A positivity-preserving implicit-explicit scheme with high order polynomial basis for compressible Navier–Stokes equations.

Authors :
Liu, Chen
Zhang, Xiangxiong
Source :
Journal of Computational Physics. Nov2023, Vol. 493, pN.PAG-N.PAG. 1p.
Publication Year :
2023

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., Δ t = O (Δ 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 Δ t = O (Δ x) if the Laplacian operator of internal energy is approximated by the Q k spectral element method with k = 1 , 2 , 3. So the full scheme with Q k (k = 1 , 2 , 3) basis is conservative and positivity-preserving with Δ t = O (Δ 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. • The hyperbolic subproblem is solved by a positivity-preserving RKDG method. • The parabolic subproblem is solved by classical schemes to achieve positivity. • A proper choice of the IPDG method achieves total conservation. • Numerically a spatially higher order scheme produces better solutions. [ABSTRACT FROM AUTHOR]

Details

Language :
English
ISSN :
00219991
Volume :
493
Database :
Academic Search Index
Journal :
Journal of Computational Physics
Publication Type :
Academic Journal
Accession number :
172809911
Full Text :
https://doi.org/10.1016/j.jcp.2023.112496