On some implicit and semi-implicit staggered schemes for the shallow water and Euler equations

¹ Aix-Marseille Université, CNRS, Centrale Marseille, I2M, UMR 7373, 13453 Marseille, France.
raphaele.herbin@univ-amu.fr
² Institut de Radioprotection et de Sûreté Nucléaire (IRSN), BP 3, 13115 Saint-Paul-lez-Durance cedex, France.
kheriji.walid@gmail.com; jean-claude.latche@irsn.fr

Received: 1 January 2013
Revised: 16 January 2014

Abstract

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.