High order pressure-based semi-implicit IMEX schemes for the 3D Navier-Stokes equations at all Mach numbers

• Published in: Journal of computational physics Vol. 434; p. 110206
• Main Authors: Boscheri, Walter, Pareschi, Lorenzo
• Format: Journal Article
• Language: English
• Published: Cambridge Elsevier Inc 01.06.2021; Elsevier Science Ltd
• Subjects: 3D compressible Euler and Navier-Stokes equations; All Mach number flow solver; Approximation; Asymptotic preserving methods; Compressibility; Computational fluid dynamics; Computational physics; Discretization; Energy dissipation; Enthalpy; Equations of state; Fluid flow; Fluxes; General equation of state (EOS); High Mach number; Ideal gas; Iterative methods; Kinetic energy; Mach number; Mathematical analysis; Navier-Stokes equations; Newton methods; Nonlinear systems; Quadrature-free WENO; Quadratures; Robustness (mathematics); Semi-implicit IMEX schemes; Solvers; Viscous flow; Quadrature-free WENO; Semi-implicit IMEX schemes; Asymptotic preserving methods; General equation of state (EOS); All Mach number flow solver; 3D compressible Euler and Navier-Stokes equations
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2021.110206

This article aims at developing a high order pressure-based solver for the solution of the 3D compressible Navier-Stokes system at all Mach numbers. We propose a cell-centered discretization of the governing equations that splits the fluxes into a fast and a slow scale part, that are treated implicitly and explicitly, respectively. A novel semi-implicit discretization is proposed for the kinetic energy as well as the enthalpy fluxes in the energy equation, hence avoiding any need of iterative solvers. The implicit discretization yields an elliptic equation on the pressure that can be solved for both ideal gas and general equation of state (EOS). A nested Newton method is used to solve the mildly nonlinear system for the pressure in case of nonlinear EOS. High order in time is granted by implicit-explicit (IMEX) time stepping, whereas a novel CWENO technique efficiently implemented in a dimension-by-dimension manner is developed for achieving high order in space for the discretization of explicit convective and viscous fluxes. A quadrature-free finite volume solver is then derived for the high order approximation of numerical fluxes. Central schemes with no dissipation of suitable order of accuracy are finally employed for the numerical approximation of the implicit terms. Consequently, the CFL-type stability condition on the maximum admissible time step is based only on the fluid velocity and not on the sound speed, so that the novel schemes work uniformly for all Mach numbers. Convergence and robustness of the proposed method are assessed through a wide set of benchmark problems involving low and high Mach number regimes, as well as inviscid and viscous flows. •All Mach solver for the Navier-Stokes equations on collocated 3D Cartesian grid.•Asymptotic Preserving discretization.•High order space-time method.•General EOS (equation of states).•Quadrature-free CWENO finite volume scheme for the explicit sub-systems.