Commun. Comput. Phys., 10 (2011), pp. 940-978.


Staggered Lagrangian Discretization Based on Cell-Centered Riemann Solver and Associated Hydrodynamics Scheme

Pierre-Henri Maire 1, Raphael Loubere 2*, Pavel Vachal 3

1 CEA-CESTA, BP 2, 33114, Le Barp, France.
2 Institut de Mathematiques de Toulouse (IMT) Universite Paul-Sabatier and CNRS, Toulouse, France.
3 FNSPE, Czech Technical University (CTU) in Prague, Czech Republic.

Received 17 March 2010; Accepted (in revised version) 25 November 2010
Available online 22 June 2011
doi:10.4208/cicp.170310.251110a

Abstract

The aim of the present work is to develop a general formalism to derive staggered discretizations for Lagrangian hydrodynamics on two-dimensional unstructured grids. To this end, we make use of the compatible discretization that has been initially introduced by E. J. Caramana et al., in J. Comput. Phys., 146 (1998). Namely, momentum equation is discretized by means of subcell forces and specific internal energy equation is obtained using total energy conservation. The main contribution of this work lies in the fact that the subcell force is derived invoking Galilean invariance and thermodynamic consistency. That is, we deduce a general form of the sub-cell force so that a cell entropy inequality is satisfied. The subcell force writes as a pressure contribution plus a tensorial viscous contribution which is proportional to the difference between the nodal velocity and the cell-centered velocity. This cell-centered velocity is a supplementary degree of freedom that is solved by means of a cell-centered approximate Riemann solver. To satisfy the second law of thermodynamics, the local subcell tensor involved in the viscous part of the subcell force must be symmetric positive definite. This subcell tensor is the cornerstone of the scheme. One particular expression of this tensor is given. A high-order extension of this discretization is provided. Numerical tests are presented in order to assess the efficiency of this approach. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of this scheme.

AMS subject classifications: 52B10, 65D18, 68U05, 68U07

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Lagrangian scheme, Riemann solver, artificial viscosity, hydrodynamics, compressible flow.

*Corresponding author.
Email: maire@celia.u-bordeaux1.fr (P.-H. Maire), raphael.loubere@math.univ-toulouse.fr (R. Loubere), vachal@galileo.fjfi.cvut.cz (P. Vachal)
 

The Global Science Journal