Uncertainty quantifiction in numerical aerodynamics

18
Uncertainty Quantification in numerical Aerodynamics Center for Uncertainty Quantification Dr. A. Litvinenko http://sri-uq.kaust.edu.sa/ September 29, 2013

Transcript of Uncertainty quantifiction in numerical aerodynamics

Page 1: Uncertainty quantifiction in numerical aerodynamics

Uncertainty Quantification in numerical Aerodynamics

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

Dr. A. Litvinenkohttp://sri-uq.kaust.edu.sa/

September 29, 2013

Page 2: Uncertainty quantifiction in numerical aerodynamics

4*

Introduction to uncertainties

Uncertain Input:Variables (α, Ma), geometry of airfoil

Uncertain solution:

1. statistical moments of (v , p, ρ), exceedance probabilitiesP(ρ > ρ∗)

2. probability density functions of CL, CD , position of shock.

Our aims:

1. Sparse representation of the input data (random fields)

2. Computing process in a reasonable time

3. Use the deterministic solver as a black box

4. A sparse format for the solution

5. Postprocessing in the sparse formatCenter for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

2 / 21

Page 3: Uncertainty quantifiction in numerical aerodynamics

4*

Stochastical Methods Overview

1. Monte Carlo Simulations (easy to implement, parallelisable,expensive, dim. indepen.).

2. Stoch. collocation methods with global or local polynomials(easy to implement, parallelisable, cheaper than MC, dim.depen.).

3. Stochastic Galerkin (difficult to implement, non-trivialparallelisation, the cheapest from all, dim. depen.)

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

3 / 21

Page 4: Uncertainty quantifiction in numerical aerodynamics

4*

Example 1:

500 MC realisations of cp in dependence on αi and Mai

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

4 / 21

Page 5: Uncertainty quantifiction in numerical aerodynamics

4*

Example 2:

500 MC realisations of cf in dependence on αi and Mai

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

5 / 21

Page 6: Uncertainty quantifiction in numerical aerodynamics

4*

Example 3:

5% and 95% quantiles for cp from 500 MC realisations.

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

6 / 21

Page 7: Uncertainty quantifiction in numerical aerodynamics

4*

Example 4:

5% and 95% quantiles for cf from 500 MC realisations.

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

7 / 21

Page 8: Uncertainty quantifiction in numerical aerodynamics

4*

Quantification of uncertainties of the solution

Compute difference with the deterministic solution (α = 2.74,Ma = 0.73)

(a) det.turb. viscos. ∈ (2e − 5, 3.6e −4)

(b) det. density ∈ (0.45, 1.25)

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

8 / 21

Page 9: Uncertainty quantifiction in numerical aerodynamics

4*

Comparison of determ. and uncert. solutions

(c) det. pressure ∈ (0.5, 1.4) (d) det.kin.energy ∈ (2e− 3, 3.8e− 2)

(e) det. velocity. in x ∈ (0, 1.4) (f) det. velocity in z ∈ (−0.2, 0.7)

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

9 / 21

Page 10: Uncertainty quantifiction in numerical aerodynamics

4*

Conclusion 1:

I It is important to model and to quantify uncertainties.

I One can see significant influence of uncertainties on thesolution.

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

10 / 21

Page 11: Uncertainty quantifiction in numerical aerodynamics

4*

Part II. Compressing of data

I Compressing of N snapshots

[q(x ,θ1), ..., q(x ,θN)] ≈ ABT .

influence on the mean value and variance of q(x ,θ) ?

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

11 / 21

Page 12: Uncertainty quantifiction in numerical aerodynamics

4*

Approximate mean of the pressure

(a) (k=5) (b) (k=10)

(c) (k=30) (d) (k=50)

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

12 / 21

Page 13: Uncertainty quantifiction in numerical aerodynamics

4*

Approximate variance of the pressure

(a) (k=5) (b) (k=10)

(c) (k=30) (d) (k=50)

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

13 / 21

Page 14: Uncertainty quantifiction in numerical aerodynamics

4*

Case 1, density

(a) |ρ̄PCE137 − ρ̄MC |F = 0.5534relative error: 0.00222

(b) |ρ̄PCE281 − ρ̄MC |F = 0.5670relative error: 0.00228

(c) |var(ρPCE137) − var(ρMC )|F =0.2503

(d) |var(ρPCE281) − var(ρMC )|F =0.2502

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

14 / 21

Page 15: Uncertainty quantifiction in numerical aerodynamics

4*

Case 1, Pressure

(a) |p̄PCE137 − p̄MC |F = 0.7114relative error: 0.00282

(b) |p̄PCE281 − p̄MC |F = 0.7240relative error: 0.00287

(c) |var(pPCE137) − var(pMC )|F =0.4714

(d) |var(pPCE281) − var(pMC )|F =0.4713

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

15 / 21

Page 16: Uncertainty quantifiction in numerical aerodynamics

4*

Case 9, density

(a) |ρ̄PCE29 − ρ̄MC |F = 7.0045relative error: 0.01431

(b) |ρ̄PCE201 − ρ̄MC |F = 7.2286relative error: 0.01477

(c) |var(ρPCE29) − var(ρMC )|F =0.7885

(d) |var(ρPCE201) − var(ρMC )|F =0.6205

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

16 / 21

Page 17: Uncertainty quantifiction in numerical aerodynamics

4*

Case 9, Pressure

(a) |p̄PCE29 − p̄MC |F = 9.0311relative error: 0.01820

(b) |p̄PCE201 − p̄MC |F = 9.3225relative error: 0.01878

(c) |var(pPCE29) − var(pMC )|F =1.3030

(d) |var(pPCE201) − var(pMC )|F =1.0294

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

17 / 21

Page 18: Uncertainty quantifiction in numerical aerodynamics

4*

Results are published

total 390 pages, our results on pp. 265-281

Center for UncertaintyQuantification

Center for UncertaintyQuantification

Center for Uncertainty Quantification Logo Lock-up

18 / 21