Sensitivity analysis and numerical diffusion effects for hyperbolic PDE systems with discontinuous solutions. The case of barotropic Euler equations in Lagrangian coordinates - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue SIAM Journal on Scientific Computing Année : 2018

Sensitivity analysis and numerical diffusion effects for hyperbolic PDE systems with discontinuous solutions. The case of barotropic Euler equations in Lagrangian coordinates

Résumé

Sensitivity analysis (SA) is the study of how the output of a mathematical model is affected by changes in the inputs. SA is widely studied, due to its many applications: uncertainty quantifi-cation, quick evaluation of close solutions, and optimization, to name but a few. In this work we show that the classic SA techniques, in particular the continuous sensitivity equation (CSE) method, cannot be used if the mathematical model is a system of hyperbolic partial differential equations (PDEs) with discontinuous solutions. The problem arises from the fact that the CSE method requires the differentiation of the state variable: if the latter is discontinuous, this in turns generates Dirac delta functions in the sensitivity. The focus of the first part of this work is to define a system of sensitivity equations valid also in case of discontinuous state: in order to do that, we add a correction term based on the Rankine-Hugoniot conditions. In the second part of this work we illustrate with some numerical tests how some classical finite volume schemes (an exact Godunov method and a Roe-type method) do not converge to the analytical solution for the sensitivity, due to the numerical diffusion: for this reason, we present an anti-diffusive numerical scheme, which provides the correct results for the sensitivity. In this work we carry out the computation in detail for the barotropic Euler equations in Lagrangian coordinates, but the approach is general and can be applied to any hyperbolic system with discontinuous solutions.
Fichier principal
Vignette du fichier
SIAM_SC_18.pdf (763.52 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-01589337 , version 1 (18-09-2017)
hal-01589337 , version 2 (14-12-2018)

Identifiants

Citer

Christophe Chalons, Régis Duvigneau, Camilla Fiorini. Sensitivity analysis and numerical diffusion effects for hyperbolic PDE systems with discontinuous solutions. The case of barotropic Euler equations in Lagrangian coordinates. SIAM Journal on Scientific Computing, 2018, 40 (6), pp.A3955-A3981. ⟨10.1137/17M1140807⟩. ⟨hal-01589337v2⟩
542 Consultations
520 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More