In this paper, we propose implicit and semi-implicit in time finite volume schemes for the barotropic Euler equations (hence, as a particular case, for the shallow water equations) and for the full Euler equations, based on staggered discretizations. For structured meshes, we use the MAC finite volume scheme, and, for general mixed quadrangular/hexahedral and simplicial meshes, we use the discrete unknowns of the Rannacher-Turek or Crouzeix-Raviart finite elements. We first show that a solution to each of these schemes satisfies a discrete kinetic energy equation. In the barotropic case, a solution also satisfies a discrete elastic potential balance; integrating these equations over the domain readily yields discrete counterparts of the stability estimates which are known for the continuous problem. In the case of the full Euler equations, the scheme relies on the discretization of the internal energy balance equation, which offers two main advantages: first, we avoid the space discretization of the total energy, which involves cell-centered and face-centered variables; second, we obtain an algorithm which boils down to a usual pressure correction scheme in the incompressible limit. Consistency (in a weak sense) with the original total energy conservative equation is obtained thanks to corrective terms in the internal energy balance, designed to compensate numerical dissipation terms appearing in the discrete kinetic energy inequality. It is then shown in the 1D case, that, supposing the convergence of a sequence of solutions, the limit is an entropy weak solution of the continuous problem in the barotropic case, and a weak solution in the full Euler case. Finally, we present numerical results which confirm this theory.
Mots clés : finite volumes, finite elements, staggered, pressure correction, Euler equations, shallow-water equations, compressible flows, analysis
@article{M2AN_2014__48_6_1807_0, author = {Herbin, R. and Kheriji, W. and Latch\'e, J.-C.}, title = {On some implicit and semi-implicit staggered schemes for the shallow water and {Euler} equations}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {1807--1857}, publisher = {EDP-Sciences}, volume = {48}, number = {6}, year = {2014}, doi = {10.1051/m2an/2014021}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2014021/} }
TY - JOUR AU - Herbin, R. AU - Kheriji, W. AU - Latché, J.-C. TI - On some implicit and semi-implicit staggered schemes for the shallow water and Euler equations JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2014 SP - 1807 EP - 1857 VL - 48 IS - 6 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2014021/ DO - 10.1051/m2an/2014021 LA - en ID - M2AN_2014__48_6_1807_0 ER -
%0 Journal Article %A Herbin, R. %A Kheriji, W. %A Latché, J.-C. %T On some implicit and semi-implicit staggered schemes for the shallow water and Euler equations %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2014 %P 1807-1857 %V 48 %N 6 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2014021/ %R 10.1051/m2an/2014021 %G en %F M2AN_2014__48_6_1807_0
Herbin, R.; Kheriji, W.; Latché, J.-C. On some implicit and semi-implicit staggered schemes for the shallow water and Euler equations. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 48 (2014) no. 6, pp. 1807-1857. doi : 10.1051/m2an/2014021. http://www.numdam.org/articles/10.1051/m2an/2014021/
[1] An L2-stable approximation of the Navier-Stokes convection operator for low-order non-conforming finite elements. Int. J. Numer. Methods Fluids 66 (2011) 555-580. | MR
, , and ,[2] Comparative study of pressure-correction and Godunov-type schemes on unsteady compressible cases. Comput. Fluids 38 (2009) 1495-1509. | MR | Zbl
, and ,[3] Notes on PCICE method: simplification, generalization and compressibility properties. J. Comput. Phys. 215 (2006) 6-11. | MR | Zbl
,[4] A unified method for computing incompressible and compressible flows in boundary-fitted coordinates. J. Comput. Phys. 141 (1998) 153-173. | MR | Zbl
and ,[5] CALIF3S. A software components library for the computation of reactive turbulent flows. Available on https://gforge.irsn.fr/gf/project/isis.
[6] Pressure method for the numerical solution of transient, compressible fluid flows. Int. J. Numer. Methods Fluids 4 (1984) 1001-1012. | Zbl
and ,[7] Numerical solution of the Navier-Stokes equations. Math. Comput. 22 (1968) 745-762. | MR | Zbl
,[8] Basic error estimates for elliptic problems, in vol. II of Handb. Numer. Anal. Edited by P. Ciarlet and J. Lions. North Holland (1991) 17-351. | MR | Zbl
,[9] Conforming and nonconforming finite element methods for solving the stationary Stokes equations. RAIRO Série Rouge 7 (1973) 33-75. | Numdam | MR | Zbl
and ,[10] I. Demirdžić, v. Lilek and M. Perić, A collocated finite volume method for predicting flows at all speeds. Int. J. Numer. Methods Fluids 16 (1993) 1029-1050. | Zbl
[11] Finite volume methods, in vol. VII of Handb. Numer. Anal. Edited by P. Ciarlet and J. Lions. North Holland (2000) 713-1020. | Zbl
, and ,[12] Convergence of the MAC scheme for the compressible Stokes equations. SIAM J. Numer. Anal. 48 (2010) 2218-2246. | MR
, , and ,[13] Dynamics of Viscous Compressible Flows. In vol. 26 of Oxford Lect. Ser. Math. Appl. Oxford University Press (2004). | MR | Zbl
,[14] An unconditionally stable pressure correction scheme for compressible barotropic Navier-Stokes equations. Math. Model. Numer. Anal. 42 (2008) 303-331. | Numdam | MR | Zbl
, , and ,[15] Staggered discretizations, pressure correction schemes and all speed barotropic flows, in Finite Volumes for Complex Applications VI − Problems and Perspectives Vol. 2, − Prague, Czech Republic (2011) 39-56. | MR | Zbl
, , , and ,[16] A discretization of phase mass balance in fractional step algorithms for the drift-flux model. IMA J. Numer. Anal. 3 (2011) 116-146. | MR
, and ,[17] Explicit high order staggered schemes for the Euler equations (2014).
, , and ,[18] An unconditionally stable pressure correction scheme for the compressible Navier-Stokes equations. Submitted (2014). | Numdam | Zbl
, , and ,[19] An overview of projection methods for incompressible flows. Comput. Methods Appl. Mech. Engrg. 195 (2006) 6011-6045. | MR | Zbl
, and ,[20] Entropy-based nonlinear viscosity for Fourier approximations of conservation laws. C.R. Acad. Sci. Paris - Série I - Analyse Numérique 346 (2008) 801-806. | MR | Zbl
and ,[21] Entropy viscosity method for nonlinear conservation laws. J. Comput. Phys. 230 (2011) 4248-4267. | MR | Zbl
, and ,[22] A projection FEM for variable density incompressible flows. J. Comput. Phys. 165 (2000) 167-188. | MR | Zbl
and ,[23] Numerical calculation of almost incompressible flow. J. Comput. Phys. 3 (1968) 80-93. | Zbl
and ,[24] A numerical fluid dynamics calculation method for all flow speeds. J. Comput. Phys. 8 (1971) 197-213. | Zbl
and ,[25] Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface. Phys. Fluids 8 (1965) 2182-2189. | Zbl
and ,[26] Staggered schemes for all speed flows. ESAIM Proc. 35 (2012) 22-150. | MR
, and ,[27] Pressure correction staggered schemes for barotropic monophasic and two-phase flows. Comput. Fluids 88 (2013) 524-542. | MR
, and ,[28] Kinetic energy control in the MAC discretization of the compressible Navier-Stokes equations. Int. J. Finites Volumes 7 (2010). | MR
and ,[29] Convergence of the MAC scheme for the steady-state incompressible Navier-Stokes equations on non-uniform grids. Proc. of Finite Volumes for Complex Applications VII − Problems and Perspectives, Berlin, Germany (2014). | MR
, and ,[30] An explicit staggered scheme for the shallow water and Euler equations. Submitted (2013).
, and ,[31] Explicit staggered schemes for the compressible euler equations. ESAIM Proc. 40 (2013) 83-102. | MR
, and ,[32] Computer simulation of reactive gas dynamics. Vol. 5 of Modeling, Identification and Control (1985) 211-236.
,[33] A robust, colocated, implicit algorithm for direct numerical simulation of compressible, turbulent flows. J. Comput. Phys. 205 (2005) 205-221. | Zbl
and ,[34] Solution of the implicitly discretised fluid flow equations by operator splitting. J. Comput. Phys. 62 (1985) 40-65. | MR | Zbl
,[35] The computation of compressible and incompressible recirculating flows by a non-iterative implicit scheme. J. Comput. Phys. 62 (1986) 66-82. | MR | Zbl
, and ,[36] Pressure-based compressible calculation method utilizing total variation diminishing schemes. AIAA J. 36 (1998) 1652-1657.
and ,[37] A second order primitive preconditioner for solving all speed multi-phase flows. J. Comput. Phys. 209 (2005) 477-503. | MR | Zbl
, , , and ,[38] Pressure based calculation procedure for viscous flows at all speeds in arbitrary configurations. AIAA J. 27 (1989) 1167-1174.
and ,[39] Characteristic-based pressure correction at all speeds. AIAA J. 34 (1996) 272-280. | Zbl
and .[40] New adaptative artificial viscosity method for hyperbolic systems of conservation laws. J. Comput. Phys. 231 (2012) 8114-8132. | MR | Zbl
and ,[41] A method for avoiding the acoustic time step restriction in compressible flow. J. Comput. Phys. 228 (2009) 4146-4161. | MR | Zbl
, , and ,[42] A convergent staggered scheme for variable density incompressible Navier-Stokes equations. Submitted (2014).
and ,[43] A pressure-based unstructured grid method for all-speed flows. Int. J. Numer. Methods Fluids 33 (2000) 355-374. | Zbl
,[44] Mathematical Topics in Fluid Mechanics - Volume 2 - Compressible Models. Vol. 10 of Oxford Lect. Ser. Math. Appl. Oxford University Press (1998). | MR | Zbl
,[45] The derivation and numerical solution of the equations for zero Mach number solution. Combust. Sci. Techn. 42 (1985) 185-205.
and .[46] The pressure-corrected ICE finite element method for compressible flows on unstructured meshes. J. Comput. Phys. 198 (2004) 659-685. | Zbl
and ,[47] Shock capturing using a pressure-correction method. AIAA J. 28 (1990) 1751-1757.
and ,[48] A high-resolution pressure-based algorithm for fluid flow at all speeds. J. Comput. Phys. 168 (2001) 101-133. | MR | Zbl
and ,[49] An efficient semi-implicit compressible solver for large-eddy simulations. J. Comput. Phys. 226 (2007) 1256-1270. | MR | Zbl
, and ,[50] Mach-uniformity through the coupled pressure and temperature correction algorithm. J. Comput. Phys. 206 (2005) 597-623. | Zbl
, and ,[51] A Mach-uniform algorithm: coupled versus segregated approach. J. Comput. Phys. 224 (2007) 314-331. | MR | Zbl
, and .[52] The Characteristic-Based Split (CBS) scheme - A unified approach to fluid dynamics. Int. J. Numer. Methods Engrg. 66 (2006) 1514-1546. | MR | Zbl
, and ,[53] Introduction to the Mathematical Theory of Compressible Flow. Vol. 27 of Oxford Lect. Ser. Math. Appl. Oxford University Press (2004). | Zbl
and ,[54] A barely implicit correction for flux-corrected transport. J. Comput. Phys. 71 (1987) 1-20. | Zbl
, , and ,[55] PELICANS, Collaborative development environment. Available on https://gforge.irsn.fr/gf/project/pelicans.
[56] A formally second order cell centered scheme for convection-diffusion equations on unstructured nonconforming grids. Int. J. Numer. Methods Fluids 71 (2013) 873-890. | MR
, , and ,[57] A pressure-based algorithm for high-speed turbomachinery flows. Int. J. Numer. Methods Fluids 25 (1997) 63-80. | Zbl
and ,[58] Simple nonconforming quadrilateral Stokes element. Numer. Methods Partial Differ. Equ. 8 (1992) 97-111. | MR | Zbl
and .[59] A time-accurate variable property algorithm for calculating flows with large temperature variations. Comput. Fluids 37 (2008) 51-63. | Zbl
and ,[60] Sur l'approximation de la solution des équations de Navier-Stokes par la méthode des pas fractionnaires II. Arch. Rat. Mech. Anal. 33 (1969) 377-385. | MR | Zbl
,[61] A multiblock operator-splitting algorithm for unsteady flows at all speeds in complex geometries. Int. J. Numer. Methods Fluids 46 (2004) 383-413. | MR | Zbl
and ,[62] Comparison of consistent explicit schemes on staggered and colocated meshes (2014).
and ,[63] Riemann solvers and numerical methods for fluid dynamics - A practical introduction, 3rd edition. Springer (2009). | MR | Zbl
,[64] Stability analysis of segregated solution methods for compressible flow. Appl. Numer. Math. 38 (2001) 257-274. | MR | Zbl
, and ,[65] A conservative pressure-correction method for flow at all speeds. Comput. Fluids 32 (2003) 1113-1132. | MR | Zbl
, and .[66] The segregated approach to predicting viscous compressible fluid flows. Trans. ASME 109 (1987) 268-277.
, and ,[67] A superlinearly convergent Mach-uniform finite volume method for the Euler equations on staggered unstructured grids. J. Comput. Phys. 217 (2006) 277-294. | MR | Zbl
, and ,[68] A semi-implicit method for resolution of acoustic waves in low Mach number flows. J. Comput. Phys. 181 (2002) 545-563. | MR | Zbl
, and ,[69] A Mach-uniform unstructured staggered grid method. Int. J. Numer. Methods Fluids 40 (2002) 1209-1235. | MR | Zbl
, and ,[70] A hybrid pressure-density-based algorithm for the Euler equations at all Mach number regimes. Int. J. Numer. Methods Fluids, online (2011).
, , and ,[71] The unified simulation for incompressible and compressible flow by the predictor-corrector scheme based on the CIP method. Comput. Phys. Commun. 119 (1999) 149-158. | Zbl
and ,[72] A general algorithm for compressible and incompressible flow - Part I. The split characteristic-based scheme. Int. J. Numer. Methods Fluids 20 (1995) 869-885. | MR | Zbl
and ,Cité par Sources :