Dec 12 – 13, 2016
Institut Henri Poincaré
Europe/Paris timezone

Locally conservative approximation of (conservative) systems written in non conservation form: application to Lagrangian hydrodynamics and multifluid problems

Dec 12, 2016, 5:30 PM
Amphi Hermite (Institut Henri Poincaré)

Amphi Hermite

Institut Henri Poincaré

11 Rue Pierre et Marie Curie, 75005 Paris, France


Prof. Rémi Abgrall (Université de Zürich)


Since the celebrated Lax Wendroff theorem, it is known that the right way of discretising systems of hyperbolic equations written in conservation form is to use a flux formulation. However, in many occasions, the relevant formu- lation, from an engineering point of view, is not to consider this conservative formulation but one non conservative form. For example, with standard notations, a one fluid model writes [équation : voir résumé PDF] (1) but the interesting quantities are the mass, velocity and pressure, which evolution is described by: [équation : voir résumé PDF] (2) Unfortunately this form is not suitable to approximation. In the case of a multi-fluid system, the same problem occurs. In this talk, we will describe a method to overcome this issue. It does not use any flux formulation per se, but can be shown to provide the right solutions. In order to illustrate the method, we will consider several examples in Eulerian and Lagrangian hydrodynamics We will first start from the Residual Distribution (RD) (re-)interpretation of the Dobrev et al. scheme [1] for the numerical solution of the Euler equations in Lagrangian form. The first ingredient of the original scheme is the staggered grid formulation which uses continuous node-based finite element approximations for the kinematic variables and cell-centered discontinuous finite elements for the thermodynamic parameters. The second ingredient of the Dobrev et al. scheme is an artificial viscosity technique applied in order to make possible the computation of strong discontinuities. Using a reformulation in term of RD scheme, we can show that the scheme is indeed locally conservative while the formulation is stricto sensu non conservative. Using this, we can generalise the construction and develop locally conservative artificial viscosity free schemes. To demonstrate the robustness of the proposed RD scheme, we solve several one-dimensional shock tube problems from rather mild to very strong ones: we go from the classical Sod problem, to TNT explosions (with JWL EOS) via the Collela-Woodward blast wave problem. In a second part, we show how to extend this method to the Eulerian framework and give applications on single fluid and multiphase problems via the five equation model. References [1] V. Dobrev, T. Kolev, and R. Rieben. High order curvilinear finite element methods for Lagrangian hydrodynamics. SIAM J. Sci. Comput, 34:B606–B641, 2012.

Presentation materials