Reliability of LES in complex applications

Transcription

1 Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27

2

3 Sample of complex flow simulation DG-FEM treatment of delta-wing: Smagorinsky

4 Connection to Wake-Vortex Hazards Airport throughput limitations: separation up to km

5 LES - Tendency toward complex applications Reliability - Error-bounds? Turbulence model/numerics?? LES error-interactions? Optimal balance computational cost - accuracy?

6 LES - Tendency toward complex applications Reliability - Error-bounds? Turbulence model/numerics?? LES error-interactions? Optimal balance computational cost - accuracy?

7 LES - Tendency toward complex applications Reliability - Error-bounds? Turbulence model/numerics?? LES error-interactions? Optimal balance computational cost - accuracy?

8 Outline Regularization modeling 2 Role of numerics 3 Pragmatic LES 4 Concluding remarks

9 Outline Regularization modeling 2 Role of numerics 3 Pragmatic LES 4 Concluding remarks

10 Energy cascading process in 3D I II III ln(e) k 5/3 k i k d ln(k) I: large-scales stirring at integral length-scale l i /k i II: inviscid nonlinear transfer inertial range E k 5/3 III: viscous dissipation dominant l d /k d

11 DNS and LES in a picture Classical problem: wide dynamic range capture both large and small scales: resolution problem N Re 9/4 ; W Re 3 If Re Re then N 75N and W W Coarsening/mathematical modeling instead: LES

12 Filtering Navier-Stokes equations Convolution-Filtering: filter-kernel G, filter-width u i = L(u i ) = G(x ξ)u(ξ) dξ Application of filter: j u j = t u i + j (u i u j ) + i p Re jju i = j (u i u j u i u j ) Turbulent stress tensor: closure problem τ ij = u i u j u i u j = L(u i u j ) L(u i )L(u j ) = [L,Π ij ](u)

13 Filtering Navier-Stokes equations Convolution-Filtering: filter-kernel G, filter-width u i = L(u i ) = G(x ξ)u(ξ) dξ Application of filter: j u j = t u i + j (u i u j ) + i p Re jju i = j (u i u j u i u j ) Turbulent stress tensor: closure problem τ ij = u i u j u i u j = L(u i u j ) L(u i )L(u j ) = [L,Π ij ](u)

14 Spatial filtering, closure problem Shorthand notation: NS(u) = NS(u) = τ(u, u) M(u) Completes basic LES formulation Find v : NS(v) = M(v) After closure system of PDE s results: dynamic range restricted to scales > Numerics does solution v resemble u? Modeling

15 Spatial filtering, closure problem Shorthand notation: NS(u) = NS(u) = τ(u, u) M(u) Completes basic LES formulation Find v : NS(v) = M(v) After closure system of PDE s results: dynamic range restricted to scales > Numerics does solution v resemble u? Modeling

16 Spatial filtering, closure problem Shorthand notation: NS(u) = NS(u) = τ(u, u) M(u) Completes basic LES formulation Find v : NS(v) = M(v) After closure system of PDE s results: dynamic range restricted to scales > Numerics does solution v resemble u? Modeling

17 Some explicit subgrid models Popular physics-based models: Dissipation: Eddy-viscosity models, e.g., Smagorinsky τ ij ν t S ij = (C S ) 2 S S ij ; effect Similarity: Inertial range, e.g., Bardina ( ) Re Re + ν t τ ij [L,Π ij ](u) = u i u j u i u j Mixed models? m ij = Bardina + C d Smagorinsky C d via dynamic Germano-Lilly procedure Alternative: mathematical, first principles modeling

18 Some explicit subgrid models Popular physics-based models: Dissipation: Eddy-viscosity models, e.g., Smagorinsky τ ij ν t S ij = (C S ) 2 S S ij ; effect Similarity: Inertial range, e.g., Bardina ( ) Re Re + ν t τ ij [L,Π ij ](u) = u i u j u i u j Mixed models? m ij = Bardina + C d Smagorinsky C d via dynamic Germano-Lilly procedure Alternative: mathematical, first principles modeling

19 Example: Leray regularization Alter convective fluxes: t u i + u j j u i + i p Re u i = LES template: t u i + j (u j u i ) + i p ( ) Re u i = j mij L Implied Leray model ( ) mij L = L u j L (u i ) u j u i = u j u i u j u i Uniquely coupled to filter L and its inverse L Rigorous analysis available ( 93s) Provides accurate LES model ( 3)

20 Example: Leray regularization Alter convective fluxes: t u i + u j j u i + i p Re u i = LES template: t u i + j (u j u i ) + i p ( ) Re u i = j mij L Implied Leray model ( ) mij L = L u j L (u i ) u j u i = u j u i u j u i Uniquely coupled to filter L and its inverse L Rigorous analysis available ( 93s) Provides accurate LES model ( 3)

21 Example: Leray regularization Alter convective fluxes: t u i + u j j u i + i p Re u i = LES template: t u i + j (u j u i ) + i p ( ) Re u i = j mij L Implied Leray model ( ) mij L = L u j L (u i ) u j u i = u j u i u j u i Uniquely coupled to filter L and its inverse L Rigorous analysis available ( 93s) Provides accurate LES model ( 3)

22 Kelvin s circulation theorem d ( ) u j dx j dt Γ(u) Re NS-α regularization Γ(u) kk u j dx j = NS eqs

23 Kelvin s circulation theorem d ( ) u j dx j dt Γ(u) Re NS-α regularization Γ(u) kk u j dx j = NS eqs Filtered Kelvin theorem: yields extended Leray model u u t t 2

24 Cascade-dynamics computability E M k 5/3 / k 3/3 k L k 3 k k NSa k DNS NS-α,Leray are dispersive Regularization alters spectrum controllable cross-over as k / : steeper than 5/3

25 Mixing layer: testing ground for models t = 2 t = 4 t = 8 Temporal at different t Spatial at different x

26 Some basic flow properties 7 6 momentum thickness A(k) time 9 Momentum thickness and energy spectrum Smagorinsky too dissipative, Bardina not enough dynamic models quite accurate problems, e.g., intermediate and smallest scales k

27 Instantaneous snapshots of spanwise vorticity x2 x x (a) x (b) x2 x x (c) x a: DNS, b: Bardina, c: Smagorinsky, d: dynamic Accuracy limited: regularization models better? (d)

28 Leray and NS-α predictions: Re = 5, = l/6 (DNS) (DNS) (Leray) Snapshot u 2 : red (blue) corresponds to up/down (NS-α)

29 Momentum thickness θ as = l/6 7 6 θ t NS-α (solid), Leray (dash), dynamic (dash-dotted), DNS approximately grid-independent at 96 3 NS-α accurate but converges slowest

30 Application to complex physics Intuitive subgrid-scale modeling: collection of assumptions/parameters increases consistency becomes problem predictive power? Can first principles approach be extended? Wall-bounded flow, Combustion, Particle-laden flow,... But how much focus should be put on SGS modeling?

31 Application to complex physics Intuitive subgrid-scale modeling: collection of assumptions/parameters increases consistency becomes problem predictive power? Can first principles approach be extended? Wall-bounded flow, Combustion, Particle-laden flow,... But how much focus should be put on SGS modeling?

32 Application to complex physics Intuitive subgrid-scale modeling: collection of assumptions/parameters increases consistency becomes problem predictive power? Can first principles approach be extended? Wall-bounded flow, Combustion, Particle-laden flow,... But how much focus should be put on SGS modeling?

33 Application to complex physics Intuitive subgrid-scale modeling: collection of assumptions/parameters increases consistency becomes problem predictive power? Can first principles approach be extended? Wall-bounded flow, Combustion, Particle-laden flow,... But how much focus should be put on SGS modeling?

34 Outline Regularization modeling 2 Role of numerics 3 Pragmatic LES 4 Concluding remarks

35 Numerics in academic LES setting Goal: approximate the unique solution to system of PDE s resulting after adopting explicit closure model General requirements: Filter separates scales > from scales < Computational grid provides additional length-scale h Require /h to be sufficiently large ( /h ) Good numerics: v(x, t :, h) v(x, t :, ) rapidly However: computational costs N 4 : implies modest /h potentially large role of numerical method in computational dynamics because of marginal resolution

36 Numerics in academic LES setting Goal: approximate the unique solution to system of PDE s resulting after adopting explicit closure model General requirements: Filter separates scales > from scales < Computational grid provides additional length-scale h Require /h to be sufficiently large ( /h ) Good numerics: v(x, t :, h) v(x, t :, ) rapidly However: computational costs N 4 : implies modest /h potentially large role of numerical method in computational dynamics because of marginal resolution

37 Numerics in academic LES setting Goal: approximate the unique solution to system of PDE s resulting after adopting explicit closure model General requirements: Filter separates scales > from scales < Computational grid provides additional length-scale h Require /h to be sufficiently large ( /h ) Good numerics: v(x, t :, h) v(x, t :, ) rapidly However: computational costs N 4 : implies modest /h potentially large role of numerical method in computational dynamics because of marginal resolution

38 Numerics in academic LES setting Goal: approximate the unique solution to system of PDE s resulting after adopting explicit closure model General requirements: Filter separates scales > from scales < Computational grid provides additional length-scale h Require /h to be sufficiently large ( /h ) Good numerics: v(x, t :, h) v(x, t :, ) rapidly However: computational costs N 4 : implies modest /h potentially large role of numerical method in computational dynamics because of marginal resolution

39 Numerics in academic LES setting Goal: approximate the unique solution to system of PDE s resulting after adopting explicit closure model General requirements: Filter separates scales > from scales < Computational grid provides additional length-scale h Require /h to be sufficiently large ( /h ) Good numerics: v(x, t :, h) v(x, t :, ) rapidly However: computational costs N 4 : implies modest /h potentially large role of numerical method in computational dynamics because of marginal resolution

40 Numerics in academic LES setting Goal: approximate the unique solution to system of PDE s resulting after adopting explicit closure model General requirements: Filter separates scales > from scales < Computational grid provides additional length-scale h Require /h to be sufficiently large ( /h ) Good numerics: v(x, t :, h) v(x, t :, ) rapidly However: computational costs N 4 : implies modest /h potentially large role of numerical method in computational dynamics because of marginal resolution

41 Modified closure problem high pass filter Discretization induces spatial filter: δ x f(x) = x ( f) Convective contribution: [ ] x (u 2 ) = δ x (u 2 ) + x (u 2 ) δ x (u 2 ) = δ x (u 2 ) + x (u 2 û2 ) = x (û2 ) + x (ξ) Modified mean flux Computational Turbulent Stress Tensor ξ = u 2 û2 = (u 2 u 2 ) + (u 2 û2 ) = τ + H(u 2 ) Numerically induced high-pass filter: H(f) = f f as r = /h

42 Modified closure problem high pass filter Discretization induces spatial filter: δ x f(x) = x ( f) Convective contribution: [ ] x (u 2 ) = δ x (u 2 ) + x (u 2 ) δ x (u 2 ) = δ x (u 2 ) + x (u 2 û2 ) = x (û2 ) + x (ξ) Modified mean flux Computational Turbulent Stress Tensor ξ = u 2 û2 = (u 2 u 2 ) + (u 2 û2 ) = τ + H(u 2 ) Numerically induced high-pass filter: H(f) = f f as r = /h

43 Modified closure problem high pass filter Discretization induces spatial filter: δ x f(x) = x ( f) Convective contribution: [ ] x (u 2 ) = δ x (u 2 ) + x (u 2 ) δ x (u 2 ) = δ x (u 2 ) + x (u 2 û2 ) = x (û2 ) + x (ξ) Modified mean flux Computational Turbulent Stress Tensor ξ = u 2 û2 = (u 2 u 2 ) + (u 2 û2 ) = τ + H(u 2 ) Numerically induced high-pass filter: H(f) = f f as r = /h

44 Modified closure problem high pass filter Discretization induces spatial filter: δ x f(x) = x ( f) Convective contribution: [ ] x (u 2 ) = δ x (u 2 ) + x (u 2 ) δ x (u 2 ) = δ x (u 2 ) + x (u 2 û2 ) = x (û2 ) + x (ξ) Modified mean flux Computational Turbulent Stress Tensor ξ = u 2 û2 = (u 2 u 2 ) + (u 2 û2 ) = τ + H(u 2 ) Numerically induced high-pass filter: H(f) = f f as r = /h

45 Dynamic importance - subgrid resolution Contributions associated with u = e ıkx : τ = A τ (k )e 2ıkx ; H(u 2 ) = A H (k, r)e 2ıkx

46 Dynamic importance - subgrid resolution Contributions associated with u = e ıkx : τ = A τ (k )e 2ıkx ; H(u 2 ) = A H (k, r)e 2ıkx A τ A H A τ : solid k

47 Dynamic importance - subgrid resolution Contributions associated with u = e ıkx : τ = A τ (k )e 2ıkx ; H(u 2 ) = A H (k, r)e 2ıkx A τ A H A τ : solid A H at r = (- -) r = 2 (-.-) and r = 4 ( ) k

48 Dynamic importance - subgrid resolution Contributions associated with u = e ıkx : τ = A τ (k )e 2ıkx ; H(u 2 ) = A H (k, r)e 2ıkx A τ A H A τ : solid A H at r = (- -) r = 2 (-.-) and r = 4 ( ) 2nd order (thin) 4th order (thick) k

49 Dynamic importance - subgrid resolution Contributions associated with u = e ıkx : τ = A τ (k )e 2ıkx ; H(u 2 ) = A H (k, r)e 2ıkx A τ A H k A τ : solid A H at r = (- -) r = 2 (-.-) and r = 4 ( ) 2nd order (thin) 4th order (thick) Strong effect r = 2; reduction as r 4

50 Observe: Numerics or modeling or both? At marginal resolution the numerics strongly modifies the equations that should be solved (alteration of mathematical nature) Likewise, the introduction of a subgrid model modifies these equations (alteration of flow-physics nature) Dilemma: unclear which is to be preferred Pragmatic guideline: minimal total error at given computational costs Governed by error-interactions: nonlinear error-accumulation... not just best model/best numerics

51 Observe: Numerics or modeling or both? At marginal resolution the numerics strongly modifies the equations that should be solved (alteration of mathematical nature) Likewise, the introduction of a subgrid model modifies these equations (alteration of flow-physics nature) Dilemma: unclear which is to be preferred Pragmatic guideline: minimal total error at given computational costs Governed by error-interactions: nonlinear error-accumulation... not just best model/best numerics

52 Observe: Numerics or modeling or both? At marginal resolution the numerics strongly modifies the equations that should be solved (alteration of mathematical nature) Likewise, the introduction of a subgrid model modifies these equations (alteration of flow-physics nature) Dilemma: unclear which is to be preferred Pragmatic guideline: minimal total error at given computational costs Governed by error-interactions: nonlinear error-accumulation... not just best model/best numerics

53 Observe: Numerics or modeling or both? At marginal resolution the numerics strongly modifies the equations that should be solved (alteration of mathematical nature) Likewise, the introduction of a subgrid model modifies these equations (alteration of flow-physics nature) Dilemma: unclear which is to be preferred Pragmatic guideline: minimal total error at given computational costs Governed by error-interactions: nonlinear error-accumulation... not just best model/best numerics

54 Observe: Numerics or modeling or both? At marginal resolution the numerics strongly modifies the equations that should be solved (alteration of mathematical nature) Likewise, the introduction of a subgrid model modifies these equations (alteration of flow-physics nature) Dilemma: unclear which is to be preferred Pragmatic guideline: minimal total error at given computational costs Governed by error-interactions: nonlinear error-accumulation... not just best model/best numerics

55 Counter-acting errors: LES-paradoxes Decaying turbulence: discretization, modeling and total-error Bernard J. Geurts: Reliability of LES in complex applications

56 Counter-acting errors: LES-paradoxes time error in kinetic energy Decaying turbulence: discretization, modeling and total-error Bernard J. Geurts: Reliability of LES in complex applications

57 Counter-acting errors: LES-paradoxes time error in kinetic energy Decaying turbulence: discretization, modeling and total-error better model may result in worse predictions higher order discretization may result in worse predictions optimal setting and refinement? Bernard J. Geurts: Reliability of LES in complex applications

58 Counter-acting errors: LES-paradoxes time error in kinetic energy Decaying turbulence: discretization, modeling and total-error better model may result in worse predictions higher order discretization may result in worse predictions optimal setting and refinement? Bernard J. Geurts: Reliability of LES in complex applications

59 Counter-acting errors: LES-paradoxes time error in kinetic energy Decaying turbulence: discretization, modeling and total-error better model may result in worse predictions higher order discretization may result in worse predictions optimal setting and refinement? Bernard J. Geurts: Reliability of LES in complex applications

60 Counter-acting errors: LES-paradoxes time error in kinetic energy Decaying turbulence: discretization, modeling and total-error better model may result in worse predictions higher order discretization may result in worse predictions optimal setting and refinement? Bernard J. Geurts: Reliability of LES in complex applications

61 Outline Regularization modeling 2 Role of numerics 3 Pragmatic LES 4 Concluding remarks

62 Experimental error-assessment Pragmatic: minimal total error at given computational costs Discuss: error-landscape/optimal refinement strategy optimality of MILES in DG-FEM LES-specific error-minimization: SIPI

63 Experimental error-assessment Pragmatic: minimal total error at given computational costs Discuss: error-landscape/optimal refinement strategy optimality of MILES in DG-FEM LES-specific error-minimization: SIPI

64 Smagorinsky fluid Homogeneous decaying turbulence at Re λ = 5, Smagorinsky fluid subgrid model: m S ij = 2(C S ) 2 S S ij = 2l 2 S S S ij introduces Smagorinsky-length l S Also dynamic Smagorinsky fluid with length-scale l d

65 Smagorinsky fluid Homogeneous decaying turbulence at Re λ = 5, Smagorinsky fluid subgrid model: m S ij = 2(C S ) 2 S S ij = 2l 2 S S S ij introduces Smagorinsky-length l S Also dynamic Smagorinsky fluid with length-scale l d

66 Accuracy measures Monitor resolved kinetic energy E = Ω 2 u u dx = u u 2 Ω Measure relative error: top-hat filter, grid h = /r E δ E (, r) = LES (, r) E DNS (, r) E DNS (, r) with error integrated over time f 2 = t t t t f 2 (t)dt each simulation represented by single number concise representation facilitates comparison

67 Accuracy measures Monitor resolved kinetic energy E = Ω 2 u u dx = u u 2 Ω Measure relative error: top-hat filter, grid h = /r E δ E (, r) = LES (, r) E DNS (, r) E DNS (, r) with error integrated over time f 2 = t t t t f 2 (t)dt each simulation represented by single number concise representation facilitates comparison

68 Error-landscape: Definition Framework for collecting error information: l S h δ E N Each Smagorinsky LES corresponds to single point: ( N, l ) S ; error : δ h E Contours of δ E fingerprint of LES

69 Error-landscape: Definition Framework for collecting error information: l S h δ E N Each Smagorinsky LES corresponds to single point: ( N, l ) S ; error : δ h E Contours of δ E fingerprint of LES

70 Total error-landscape combination of central discretization and Smagorinsky optimum at C S > : SGS modeling is viable here

71 Total error-landscape combination of central discretization and Smagorinsky optimum at C S > : SGS modeling is viable here

72 Error-landscape: optimal refinement (a) (b) Optimal trajectory for Re λ = 5 (a) and Re λ = (b) Under-resolution leads to strong error-increase Dynamic model over-estimates viscosity reduces error well with increasing resolution

73 Error-landscape: optimal refinement (a) (b) Optimal trajectory for Re λ = 5 (a) and Re λ = (b) Under-resolution leads to strong error-increase Dynamic model over-estimates viscosity reduces error well with increasing resolution

74 Observation: MILES philosophy practical LES implies marginal resolution which implies large role of specific numerical discretization next to dynamics due to subgrid model and leads to strong interactions and complex error-accumulation Proposal: obtain smoothing via appropriate numerical method alone accept that there is no grid-independent solution, other than DNS accept that predictions become discretization dependent Is no-model/just numerics option optimal/viable? Consider example: DG-FEM and homogeneous turbulence

75 Observation: MILES philosophy practical LES implies marginal resolution which implies large role of specific numerical discretization next to dynamics due to subgrid model and leads to strong interactions and complex error-accumulation Proposal: obtain smoothing via appropriate numerical method alone accept that there is no grid-independent solution, other than DNS accept that predictions become discretization dependent Is no-model/just numerics option optimal/viable? Consider example: DG-FEM and homogeneous turbulence

76 Observation: MILES philosophy practical LES implies marginal resolution which implies large role of specific numerical discretization next to dynamics due to subgrid model and leads to strong interactions and complex error-accumulation Proposal: obtain smoothing via appropriate numerical method alone accept that there is no grid-independent solution, other than DNS accept that predictions become discretization dependent Is no-model/just numerics option optimal/viable? Consider example: DG-FEM and homogeneous turbulence

77 Observation: MILES philosophy practical LES implies marginal resolution which implies large role of specific numerical discretization next to dynamics due to subgrid model and leads to strong interactions and complex error-accumulation Proposal: obtain smoothing via appropriate numerical method alone accept that there is no grid-independent solution, other than DNS accept that predictions become discretization dependent Is no-model/just numerics option optimal/viable? Consider example: DG-FEM and homogeneous turbulence

78 DG-FEM of homogeneous turbulence Discretization: Approximate Riemann solver F = F central + γf dissipative ; HLLC flux c ILES plane N ELES plane c s Three-dimensional accuracy charts

79 LES with DG-FEM: dissipative numerics γ c c s N Third order DG-FEM at Re λ = : red symbols - optimal setting

80 Optimal refinement strategies: 2nd order From dissipative to central: c s.. c s..5 c s N N N left-to-right: γ c =, γ c =., γ c =. C S = at γ c = : MILES best option γ c < implies CS : MILES sub-optimal decrease γ c implies increase CS : exchange of dissipation

81 Optimal refinement strategies: 2nd order From dissipative to central: c s.. c s..5 c s N N N left-to-right: γ c =, γ c =., γ c =. C S = at γ c = : MILES best option γ c < implies CS : MILES sub-optimal decrease γ c implies increase CS : exchange of dissipation

82 Optimal refinement strategies: 2nd order From dissipative to central: c s.. c s..5 c s N N N left-to-right: γ c =, γ c =., γ c =. C S = at γ c = : MILES best option γ c < implies CS : MILES sub-optimal decrease γ c implies increase CS : exchange of dissipation

83 Optimal refinement strategies: 2nd order From dissipative to central: c s.. c s..5 c s N N N left-to-right: γ c =, γ c =., γ c =. C S = at γ c = : MILES best option γ c < implies CS : MILES sub-optimal decrease γ c implies increase CS : exchange of dissipation

84 Optimal refinement strategies: 3rd order From dissipative to central: c s c s c s N N N left-to-right: γ c =, γ c =., γ c =. for all γ c [, ] find CS : MILES sub-optimal optimal C S is less sensitive to γ c value than 2nd order

85 Optimal refinement strategies: 3rd order From dissipative to central: c s c s c s N N N left-to-right: γ c =, γ c =., γ c =. for all γ c [, ] find CS : MILES sub-optimal optimal C S is less sensitive to γ c value than 2nd order

86 Optimal refinement strategies: 3rd order From dissipative to central: c s c s c s N N N left-to-right: γ c =, γ c =., γ c =. for all γ c [, ] find CS : MILES sub-optimal optimal C S is less sensitive to γ c value than 2nd order

87 Optimality of MILES (γ c = )? δ L2. δ L N (a) N (b) (a): 2nd order ; (b): 3rd order γ c =. (dot), γ c =. (dash) and γ c =. ( ) 2nd: MILES-error larger than with explicit SGS model 3rd: optimum requires explicit SGS model Option: directly optimize total simulation error - SIPI

88 Optimality of MILES (γ c = )? δ L2. δ L N (a) N (b) (a): 2nd order ; (b): 3rd order γ c =. (dot), γ c =. (dash) and γ c =. ( ) 2nd: MILES-error larger than with explicit SGS model 3rd: optimum requires explicit SGS model Option: directly optimize total simulation error - SIPI

89 Optimality of MILES (γ c = )? δ L2. δ L N (a) N (b) (a): 2nd order ; (b): 3rd order γ c =. (dot), γ c =. (dash) and γ c =. ( ) 2nd: MILES-error larger than with explicit SGS model 3rd: optimum requires explicit SGS model Option: directly optimize total simulation error - SIPI

90 Optimality of MILES (γ c = )? δ L2. δ L N (a) N (b) (a): 2nd order ; (b): 3rd order γ c =. (dot), γ c =. (dash) and γ c =. ( ) 2nd: MILES-error larger than with explicit SGS model 3rd: optimum requires explicit SGS model Option: directly optimize total simulation error - SIPI

91 SIPI - basic algorithm Goal: minimize total error at given N δe C S a b d c Initial triplet: no-model, dynamic and half-way New iterand CS d = b (b a) 2 [δ E (b) δ E (c)] (b c) 2 [δ E (b) δ E (a)] 2 (b a)[δ E (b) δ E (c)] (b c)[δ E (b) δ E (a)]

92 SIPI applied to homogeneous turbulence Each iteration = separate simulation (a) (b) Re λ = 5 (a) and Re λ = (b). Resolutions N = 24 (solid), N = 32 (dashed) and N = 48 (dash-dotted) Iterations: +

93 Convergence example Re λ = 5 Re λ = n C (n) S (24) C(n) S (48) C(n) S (24) C(n) S (48) Computational overhead SIPI: CPU-time T N 4 implies approximate optimization at N can be (almost) completed within cost of one simulation at 3N/2

94 Outline Regularization modeling 2 Role of numerics 3 Pragmatic LES 4 Concluding remarks

95 Concluding remarks Regularization modeling LES filtering approach allows grid-independent LES closure problem: intuitive or first principles modeling Leray and LANS-α: accuracy and extension to complex physics Leray robust, LANS-α (slightly) more accurate turbulent combustion: Leray not (yet) more accurate but more systematic

96 Concluding remarks Regularization modeling LES filtering approach allows grid-independent LES closure problem: intuitive or first principles modeling Leray and LANS-α: accuracy and extension to complex physics Leray robust, LANS-α (slightly) more accurate turbulent combustion: Leray not (yet) more accurate but more systematic

97 Concluding remarks Regularization modeling LES filtering approach allows grid-independent LES closure problem: intuitive or first principles modeling Leray and LANS-α: accuracy and extension to complex physics Leray robust, LANS-α (slightly) more accurate turbulent combustion: Leray not (yet) more accurate but more systematic

98 Concluding remarks Regularization modeling LES filtering approach allows grid-independent LES closure problem: intuitive or first principles modeling Leray and LANS-α: accuracy and extension to complex physics Leray robust, LANS-α (slightly) more accurate turbulent combustion: Leray not (yet) more accurate but more systematic

99 Concluding remarks Regularization modeling LES filtering approach allows grid-independent LES closure problem: intuitive or first principles modeling Leray and LANS-α: accuracy and extension to complex physics Leray robust, LANS-α (slightly) more accurate turbulent combustion: Leray not (yet) more accurate but more systematic

100 Concluding remarks Grid-independent LES for computational error-assessment error-decomposition: modeling, discretization effects LES-paradoxes and interacting errors: better models/numerics may not lead to better predictions error-landscape optimal refinement strategy dynamic procedure efficient error-reduction Error-interaction and a priori error-bounds hard to include: SIPI to account for modeling and numerics Thanks: Johan Meyers (Leuven), Fedderik van der Bos (Munich), Darryl Holm (London)

101 Concluding remarks Grid-independent LES for computational error-assessment error-decomposition: modeling, discretization effects LES-paradoxes and interacting errors: better models/numerics may not lead to better predictions error-landscape optimal refinement strategy dynamic procedure efficient error-reduction Error-interaction and a priori error-bounds hard to include: SIPI to account for modeling and numerics Thanks: Johan Meyers (Leuven), Fedderik van der Bos (Munich), Darryl Holm (London)

102 Concluding remarks Grid-independent LES for computational error-assessment error-decomposition: modeling, discretization effects LES-paradoxes and interacting errors: better models/numerics may not lead to better predictions error-landscape optimal refinement strategy dynamic procedure efficient error-reduction Error-interaction and a priori error-bounds hard to include: SIPI to account for modeling and numerics Thanks: Johan Meyers (Leuven), Fedderik van der Bos (Munich), Darryl Holm (London)

103 Concluding remarks Grid-independent LES for computational error-assessment error-decomposition: modeling, discretization effects LES-paradoxes and interacting errors: better models/numerics may not lead to better predictions error-landscape optimal refinement strategy dynamic procedure efficient error-reduction Error-interaction and a priori error-bounds hard to include: SIPI to account for modeling and numerics Thanks: Johan Meyers (Leuven), Fedderik van der Bos (Munich), Darryl Holm (London)

104 Concluding remarks Grid-independent LES for computational error-assessment error-decomposition: modeling, discretization effects LES-paradoxes and interacting errors: better models/numerics may not lead to better predictions error-landscape optimal refinement strategy dynamic procedure efficient error-reduction Error-interaction and a priori error-bounds hard to include: SIPI to account for modeling and numerics Thanks: Johan Meyers (Leuven), Fedderik van der Bos (Munich), Darryl Holm (London)

105 Concluding remarks Grid-independent LES for computational error-assessment error-decomposition: modeling, discretization effects LES-paradoxes and interacting errors: better models/numerics may not lead to better predictions error-landscape optimal refinement strategy dynamic procedure efficient error-reduction Error-interaction and a priori error-bounds hard to include: SIPI to account for modeling and numerics Thanks: Johan Meyers (Leuven), Fedderik van der Bos (Munich), Darryl Holm (London)