Sensitivity equation method for Euler equations in presence of shocks applied to uncertainty quantification - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2018

Sensitivity equation method for Euler equations in presence of shocks applied to uncertainty quantification

Résumé

Sensitivity analysis (SA) is the study of how changes in the input of a model affect the outputs. Standard sensitivity analysis techniques, such as the continuous sensitivity equation (CSE) 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 theoretical 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: an HLL-type scheme and an HLLC-type scheme. The effects of the numerical diffusion on the convergence of the schemes with respect to the grid are discussed. Finally, an application of SA to uncertainty quantification is investigated and the different numerical schemes are compared.
Fichier principal
Vignette du fichier
SA_Euler_compl.pdf (1.19 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 1

Citer

Camilla Fiorini, Christophe Chalons, Régis Duvigneau. Sensitivity equation method for Euler equations in presence of shocks applied to uncertainty quantification. 2018. ⟨hal-01817815v1⟩
433 Consultations
633 Téléchargements

Partager

Gmail Facebook X LinkedIn More