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

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.

[1]  P. Degond,et al.  All speed scheme for the low Mach number limit of the isentropic Euler equations , 2009, 0908.1929.

[2]  Stéphane Dellacherie,et al.  Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number , 2010, J. Comput. Phys..

[3]  Raphaël Loubère,et al.  A second order all Mach number IMEX finite volume solver for the three dimensional Euler equations , 2020, J. Comput. Phys..

[4]  Ilya Peshkov,et al.  On a pure hyperbolic alternative to the Navier-Stokes equations , 2014 .

[5]  Vincenzo Casulli,et al.  Semi-implicit finite difference methods for the two-dimensional shallow water equation , 1990 .

[6]  Lorenzo Pareschi,et al.  A Unified IMEX Runge-Kutta Approach for Hyperbolic Systems with Multiscale Relaxation , 2017, SIAM J. Numer. Anal..

[7]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[8]  G. Russo,et al.  High Order Asymptotically Strong-Stability-Preserving Methods for Hyperbolic Systems with Stiff Relaxation , 2003 .

[9]  Shi Jin,et al.  Efficient Asymptotic-Preserving (AP) Schemes For Some Multiscale Kinetic Equations , 1999, SIAM J. Sci. Comput..

[10]  Michael Dumbser,et al.  ADER-WENO finite volume schemes with space-time adaptive mesh refinement , 2012, J. Comput. Phys..

[11]  Michael Dumbser,et al.  Semi-implicit discontinuous Galerkin methods for the incompressible Navier–Stokes equations on adaptive staggered Cartesian grids , 2016, 1612.09558.

[12]  Claus-Dieter Munz,et al.  On Godunov-type schemes for Lagrangian gas dynamics , 1994 .

[13]  Vincenzo Casulli,et al.  A SEMI-IMPLICIT FINITE DIFFERENCE METHOD FOR NON-HYDROSTATIC, FREE-SURFACE FLOWS , 1999 .

[14]  Dinshaw S. Balsara,et al.  High order direct Arbitrary-Lagrangian-Eulerian (ALE) PNPM schemes with WENO Adaptive-Order reconstruction on unstructured meshes , 2019, J. Comput. Phys..

[15]  Gabriella Puppo,et al.  CWENO: Uniformly accurate reconstructions for balance laws , 2016, Math. Comput..

[16]  Claus-Dieter Munz,et al.  Coupling of compressible and incompressible flow regions using the multiple pressure variables approach , 2015 .

[17]  Fabrice Deluzet,et al.  An Asymptotic Preserving scheme for the Euler equations in a strong magnetic field , 2008, J. Comput. Phys..

[18]  A. Stroud.  Approximate calculation of multiple integrals , 1973 .

[19]  Michael Dumbser,et al.  Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , 2008, J. Comput. Phys..

[20]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[21]  Lorenzo Pareschi,et al.  Implicit-explicit runge-kutta schemes and applications to hyperbolic systems with relaxation , 2010, 1009.2757.

[22]  Tao Xiong,et al.  A high order semi-implicit IMEX WENO scheme for the all-Mach isentropic Euler system , 2019, J. Comput. Phys..

[23]  Raphaël Loubère,et al.  Second-order implicit-explicit total variation diminishing schemes for the Euler system in the low Mach regime , 2017, J. Comput. Phys..

[24]  E. Hofer,et al.  A Partially Implicit Method for Large Stiff Systems of ODEs with Only Few Equations Introducing Small Time-Constants , 1976 .

[25]  Michael Dumbser,et al.  Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier–Stokes equations , 2010 .

[26]  Michael Dumbser,et al.  Space-time adaptive ADER-DG schemes for dissipative flows: Compressible Navier-Stokes and resistive MHD equations , 2016, Comput. Phys. Commun..

[27]  Eleuterio F. Toro,et al.  Flux splitting schemes for the Euler equations , 2012 .

[28]  Dinshaw S. Balsara,et al.  Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics , 2012, J. Comput. Phys..

[29]  Rémi Abgrall,et al.  An adaptive shock-capturing algorithm for solving unsteady reactive flows , 2003 .

[30]  W. Boscheri.  A space-time semi-Lagrangian advection scheme on staggered Voronoi meshes applied to free surface flows , 2020 .

[31]  Michael Dumbser,et al.  A conservative, weakly nonlinear semi-implicit finite volume scheme for the compressible Navier-Stokes equations with general equation of state , 2016, Appl. Math. Comput..

[32]  Philip M. Gresho,et al.  On the theory of semi‐implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. Part 1: Theory , 1990 .

[33]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[34]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[35]  Claus-Dieter Munz,et al.  A contribution to the construction of diffusion fluxes for finite volume and discontinuous Galerkin schemes , 2007, J. Comput. Phys..

[36]  R. Klein.  Semi-implicit extension of a Godunov-type scheme based on low Mach number asymptotics , 1995 .

[37]  E. Turkel,et al.  Preconditioned methods for solving the incompressible low speed compressible equations , 1987 .

[38]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[39]  G. Russo,et al.  Central WENO schemes for hyperbolic systems of conservation laws , 1999 .

[40]  P. Roe,et al.  On Godunov-type methods near low densities , 1991 .

[41]  Samuel Kokh,et al.  Large Time Step and Asymptotic Preserving Numerical Schemes for the Gas Dynamics Equations with Source Terms , 2013, SIAM J. Sci. Comput..

[42]  S. Orszag,et al.  Small-scale structure of the Taylor–Green vortex , 1983, Journal of Fluid Mechanics.

[43]  Francis Filbet,et al.  High Order Semi-implicit Schemes for Time Dependent Partial Differential Equations , 2016, Journal of Scientific Computing.

[44]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[45]  Pierre Degond,et al.  An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations , 2011, J. Comput. Phys..

[46]  Florian Bernard,et al.  Linearly implicit all Mach number shock capturing schemes for the Euler equations , 2019, J. Comput. Phys..

[47]  H. Guillard,et al.  On the behavior of upwind schemes in the low Mach number limit: II. Godunov type schemes , 2004 .

[48]  M. Dumbser,et al.  An efficient semi-implicit finite volume method for axially symmetric compressible flows in compliant tubes , 2015 .

[49]  Santiago Badia,et al.  Assessment of variational multiscale models for the large eddy simulation of turbulent incompressible flows , 2015 .

[50]  Michael Dumbser,et al.  A staggered space-time discontinuous Galerkin method for the three-dimensional incompressible Navier-Stokes equations on unstructured tetrahedral meshes , 2016, J. Comput. Phys..

[51]  M. Dumbser,et al.  Semi-implicit staggered discontinuous Galerkin schemes for axially symmetric viscous compressible flows in elastic tubes , 2018 .

[52]  Luigi Brugnano,et al.  Iterative Solution of Piecewise Linear Systems , 2007, SIAM J. Sci. Comput..

[53]  Michael Dumbser,et al.  Central Weighted ENO Schemes for Hyperbolic Conservation Laws on Fixed and Moving Unstructured Meshes , 2017, SIAM J. Sci. Comput..

[54]  Axel Klar.  An Asymptotic Preserving Numerical Scheme for Kinetic Equations in the Low Mach Number Limit , 1999 .

[55]  T. Alazard,et al.  Low Mach Number Limit of the Full Navier-Stokes Equations , 2005, math/0501386.

[56]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[57]  Michael Dumbser,et al.  Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems , 2007, J. Comput. Phys..

[58]  With Invariant Submanifolds,et al.  Systems of Conservation Laws , 2009 .

[59]  Song Jiang,et al.  Low Mach number limit of full Navier–Stokes equations in a 3D bounded domain , 2015 .

[60]  Giovanni Russo,et al.  On a Class of Uniformly Accurate IMEX Runge--Kutta Schemes and Applications to Hyperbolic Systems with Relaxation , 2009, SIAM J. Sci. Comput..

[61]  Katja Bachmeier,et al.  Numerical Heat Transfer And Fluid Flow , 2016 .