A modified sensitivity equation method for the Euler equations in presence of shocks - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue Numerical Methods for Partial Differential Equations Année : 2019

A modified sensitivity equation method for the Euler equations in presence of shocks

Résumé

The Continuous Sensitivity Equation (CSE) method allows to quantify how changes in the input of a Partial Differential Equation (PDE) model affect the outputs, by solving additional PDEs ob- tained by differentiating the model. However, this method cannot be used directly in the framework of hyperbolic PDE systems with discontinuous solution, because it yields Dirac delta functions in the sensitivity solution at the location of state discontinuities. This difficulty is well known from theoret- ical viewpoint, but only a few works can be found in the literature regarding the possible numerical treatment. Therefore, we investigate in this study how classical numerical schemes for compressible Euler equations can be modified to account for shocks when computing the sensitivity solution. In particular, we propose the introduction of a source term, that allows to remove the spikes associated to the Dirac delta functions in the numerical solution. Numerical studies exhibit a strong impact of the numerical diffusion on the accuracy of this strategy. Therefore, we propose an anti-diffusive numerical scheme coupled with the approximate Riemann solver of Roe for the state problem. For the sensitivity problem, two different numerical schemes are implemented and compared: one which takes into account the contact wave and another that neglects it. The effects of the numerical diffu- sion on the convergence of the schemes with respect to the grid are discussed. Finally, an application to uncertainty propagation is investigated and the different numerical schemes are compared.
Fichier principal
Vignette du fichier
NumPDE19.pdf (1.62 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-01817815 , version 1 (18-06-2018)
hal-01817815 , version 2 (09-12-2019)
hal-01817815 , version 3 (11-12-2019)

Identifiants

  • HAL Id : hal-01817815 , version 2

Citer

Camilla Fiorini, Christophe Chalons, Régis Duvigneau. A modified sensitivity equation method for the Euler equations in presence of shocks. Numerical Methods for Partial Differential Equations, In press. ⟨hal-01817815v2⟩
433 Consultations
635 Téléchargements

Partager

Gmail Facebook X LinkedIn More