Reliability of LES in complex applications
|
|
- Luke Wilson
- 5 years ago
- Views:
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)
Regularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationDatabase analysis of errors in large-eddy simulation
PHYSICS OF FLUIDS VOLUME 15, NUMBER 9 SEPTEMBER 2003 Johan Meyers a) Department of Mechanical Engineering, Katholieke Universiteit Leuven, Celestijnenlaan 300, B3001 Leuven, Belgium Bernard J. Geurts b)
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More information43rd AIAA Aerospace Sciences Meeting and Exhibit, Jan 2005, Reno, Nevada
43rd AIAA Aerospace Sciences Meeting and Exhibit, 10-13 Jan 2005, Reno, Nevada A Dynamic Procedure for the Lagrangian Averaged Navier-Stokes-α Model of Turbulent Flows Kamran Mohseni and Hongwu Zhao Aerospace
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationMULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
More informationEnsemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3
Center for Turbulence Research Proceedings of the Summer Program 1996 237 Ensemble averaged dynamic modeling By D. Carati 1,A.Wray 2 AND W. Cabot 3 The possibility of using the information from simultaneous
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationA priori analysis of an anisotropic subfilter model for heavy particle dispersion
Workshop Particles in Turbulence, COST Action MP0806 Leiden, May 14-16, 2012 A priori analysis of an anisotropic subfilter model for heavy particle dispersion Maria Knorps & Jacek Pozorski Institute of
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
More informationCOMMUTATION ERRORS IN PITM SIMULATION
COMMUTATION ERRORS IN PITM SIMULATION B. Chaouat ONERA, 93 Châtillon, France Bruno.Chaouat@onera.fr Introduction Large eddy simulation is a promising route. This approach has been largely developed in
More informationL.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133, ORSAY CEDEX, FRANCE fax number:
Large Eddy Simulations of a spatially developing incompressible 3D mixing layer using the v-ω formulation. C. TENAUD, S. PELLERIN, A. DULIEU and L. TA PHUOC L.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133,
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationThree-dimensional wall filtering formulation for large-eddy simulation
Center for Turbulence Research Annual Research Briefs 6 55 Three-dimensional wall filtering formulation for large-eddy simulation By M. Shoeybi AND J. A. Templeton 1. Motivation and objectives Large-eddy
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationOn the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows
PHYSICS OF FLUIDS VOLUME 11, NUMBER 5 MAY 1999 On the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows L. Shao a) Laboratoire de Mécanique des Fluides
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationDYNAMIC SUBGRID-SCALE MODELING FOR LARGE-EDDY SIMULATION OF TURBULENT FLOWS WITH A STABILIZED FINITE ELEMENT METHOD
DYNAMIC SUBGRID-SCALE MODELING FOR LARGE-EDDY SIMULATION OF TURBULENT FLOWS WITH A STABILIZED FINITE ELEMENT METHOD By Andrés E. Tejada-Martínez A Thesis Submitted to the Graduate Faculty of Rensselaer
More informationLES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2
MAGNETOHYDRODYNAMICS Vol. 00 (1964), No. 00, pp. 1 5 LES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2 1 Institute for Electrothermal
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationSubgrid-Scale Models for Compressible Large-Eddy Simulations
Theoret. Comput. Fluid Dynamics (000 13: 361 376 Theoretical and Computational Fluid Dynamics Springer-Verlag 000 Subgrid-Scale Models for Compressible Large-Eddy Simulations M. Pino Martín Department
More informationJ.-L. Guermond 1 FOR TURBULENT FLOWS LARGE EDDY SIMULATION MODEL A HYPERVISCOSITY SPECTRAL
A HYPERVISCOSITY SPECTRAL LARGE EDDY SIMULATION MODEL FOR TURBULENT FLOWS J.-L. Guermond 1 Collaborator: S. Prudhomme 2 Mathematical Aspects of Computational Fluid Dynamics Oberwolfach November 9th to
More informationA scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric boundary layer
J. Fluid Mech. (2000), vol. 415, pp. 261 284. Printed in the United Kingdom c 2000 Cambridge University Press 261 A scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric
More informationMulti-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah
Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationcfl Copyright by Fotini V. Katopodes 2000 All rights reserved.
A theory for the subfilter-scale model in large-eddy simulation Fotini V. Katopodes, Robert L. Street, Joel H. Ferziger March, 2000 Technical Report 2000-K1 Environmental Fluid Mechanics Laboratory Stanford,
More informationDynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh
40th Fluid Dynamics Conference and Exhibit 8 June - July 00, Chicago, Illinois AIAA 00-506 Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows Xiaochuan Chai and Krishnan Mahesh University
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationModeling of turbulence in stirred vessels using large eddy simulation
Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationEstimation of Turbulent Dissipation Rate Using 2D Data in Channel Flows
Proceedings of the 3 rd World Congress on Mechanical, Chemical, and Material Engineering (MCM'17) Rome, Italy June 8 10, 2017 Paper No. HTFF 140 ISSN: 2369-8136 DOI: 10.11159/htff17.140 Estimation of Turbulent
More informationPredicting natural transition using large eddy simulation
Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact
More informationA simple subgrid-scale model for astrophysical turbulence
CfCA User Meeting NAOJ, 29-30 November 2016 A simple subgrid-scale model for astrophysical turbulence Nobumitsu Yokoi Institute of Industrial Science (IIS), Univ. of Tokyo Collaborators Axel Brandenburg
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationProject Topic. Simulation of turbulent flow laden with finite-size particles using LBM. Leila Jahanshaloo
Project Topic Simulation of turbulent flow laden with finite-size particles using LBM Leila Jahanshaloo Project Details Turbulent flow modeling Lattice Boltzmann Method All I know about my project Solid-liquid
More informationThe lattice Boltzmann equation (LBE) has become an alternative method for solving various fluid dynamic
36th AIAA Fluid Dynamics Conference and Exhibit 5-8 June 2006, San Francisco, California AIAA 2006-3904 Direct and Large-Eddy Simulation of Decaying and Forced Isotropic Turbulence Using Lattice Boltzmann
More informationA NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS
June 30 - July 3, 2015 Melbourne, Australia 9 7B-4 A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS C.-Y. Chang, S. Jakirlić, B. Krumbein and C. Tropea Institute of Fluid Mechanics and Aerodynamics /
More informationSUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION
SUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION By Sergei G. Chumakov A dissertation submitted in partial fulfillment of the requirements for the degree
More informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationA NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS
June 3 - July 3, 25 Melbourne, Australia 9 B-3 A NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS Ryan N. King Department of Mechanical Engineering University of Colorado Boulder, CO, 839, USA ryan.n.king@colorado.edu
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal PI: Sarma L. Rani Department of Mechanical and Aerospace Engineering The
More informationOn the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows
Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More informationEulerian models. 2.1 Basic equations
2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing
More informationLecture 4: The Navier-Stokes Equations: Turbulence
Lecture 4: The Navier-Stokes Equations: Turbulence September 23, 2015 1 Goal In this Lecture, we shall present the main ideas behind the simulation of fluid turbulence. We firts discuss the case of the
More informationMassimo GERMANO Politecnico di Torino
Hybrid Massimo GERMANO Politecnico di Torino Martín SÁNCHEZ-ROCHA Dassault Systèmes SIMULIA Corporation Suresh MENON Georgia Institute of Technology 64th Annual APS-DFD Meeting Baltimore, Maryland November
More informationMULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS
MULTISCALE MODEL. SIMUL. Vol. 6, No. 4, pp. 1317 1346 c 2008 Society for Industrial and Applied Mathematics MULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS
More informationA priori tests of one-equation LES modeling of rotating turbulence
Journal of Turbulence Volume 8, No. 37, 2007 A priori tests of one-equation LES modeling of rotating turbulence HAO LU, CHRISTOPHER J. RUTLAND and LESLIE M. SMITH, Engine Research Center, Department of
More informationCharacteristics of Linearly-Forced Scalar Mixing in Homogeneous, Isotropic Turbulence
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-1103 Characteristics of Linearly-Forced Scalar Mixing in Homogeneous, Isotropic Turbulence
More informationVelocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step
Copyright c 2004 Tech Science Press CMC, vol.1, no.3, pp.275-288, 2004 Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step B. Wang 1, H.Q. Zhang 1, C.K. Chan 2 and X.L.
More informationCOMPARISON OF DIFFERENT SUBGRID TURBULENCE MODELS AND BOUNDARY CONDITIONS FOR LARGE-EDDY-SIMULATIONS OF ROOM AIR FLOWS.
7 TH INTRNATINAL CNFRNC N AIR DISTRIBTIN IN RMS, RMVNT 2 pp. 31-36 CMPARISN F DIFFRNT SBGRID TRBLNC MDLS AND BNDARY CNDITINS FR LARG-DDY-SIMLATINS F RM AIR FLWS. D. Müller 1, L. Davidson 2 1 Lehrstuhl
More informationAn improved velocity increment model based on Kolmogorov equation of filtered velocity
An improved velocity increment model based on Kolmogorov equation of filtered velocity Le Fang, Liang Shao, Jean-Pierre Bertoglio, Guixiang X. Cui, Chun-Xiao Xu, Zhaoshun Zhang To cite this version: Le
More informationKey words. turbulence models, sub-grid scale models, large eddy simulations, global attractors, inviscid regularization of Euler equations.
COMM. MATH. SCI. Vol. 4, No. 4, pp. 823 848 c 2006 International Press GLOBAL WELL-POSEDNESS OF THE THREE-DIMENSIONAL VISCOUS AND INVISCID SIMPLIFIED BARDINA TURBULENCE MODELS YANPING CAO, EVELYN M. LUNASIN,
More informationLES ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT
2th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT Wang T.*, Gao S.F., Liu Y.W., Lu Z.H. and Hu H.P. *Author
More informationTurbulence modelling. Sørensen, Niels N. Publication date: Link back to DTU Orbit
Downloaded from orbit.dtu.dk on: Dec 19, 2017 Turbulence modelling Sørensen, Niels N. Publication date: 2010 Link back to DTU Orbit Citation (APA): Sørensen, N. N. (2010). Turbulence modelling. Paper presented
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal and Vijaya Rani PI: Sarma L. Rani Department of Mechanical and Aerospace
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationDynamics of the Coarse-Grained Vorticity
(B) Dynamics of the Coarse-Grained Vorticity See T & L, Section 3.3 Since our interest is mainly in inertial-range turbulence dynamics, we shall consider first the equations of the coarse-grained vorticity.
More informationSimulations for Enhancing Aerodynamic Designs
Simulations for Enhancing Aerodynamic Designs 2. Governing Equations and Turbulence Models by Dr. KANNAN B T, M.E (Aero), M.B.A (Airline & Airport), PhD (Aerospace Engg), Grad.Ae.S.I, M.I.E, M.I.A.Eng,
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationA SEAMLESS HYBRID RANS/LES MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS
A SEAMS HYBRID RANS/ MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS NUMBER FLOWS ON COARSE GRIDS P. Nguyen 1, J. Uribe 2, I. Afgan 1 and D. Laurence 1 1 School of Mechanical, Aerospace
More informationarxiv: v1 [physics.flu-dyn] 16 Nov 2018
Turbulence collapses at a threshold particle loading in a dilute particle-gas suspension. V. Kumaran, 1 P. Muramalla, 2 A. Tyagi, 1 and P. S. Goswami 2 arxiv:1811.06694v1 [physics.flu-dyn] 16 Nov 2018
More information+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2):
International Conference on Methods of Aerophysical Research, ICMAR 008 LARGE EDDY SIMULATION OF TURBULENT ROUND IMPINGING JET B.B. Ilyushin, D.V. Krasinsky Kutateladze Institute of Thermophysics SB RAS
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationOn the Euler rotation angle and axis of a subgrid-scale stress model
Center for Turbulence Research Proceedings of the Summer Program 212 117 On the Euler rotation angle and axis of a subgrid-scale stress model By M. Saeedi, B.-C. Wang AND S. T. Bose In this research, geometrical
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationSubgrid-scale modeling of helicity and energy dissipation in helical turbulence
Subgrid-scale modeling of helicity and energy dissipation in helical turbulence Yi Li, 1 Charles Meneveau, 1 Shiyi Chen, 1,2 and Gregory L. Eyink 3 1 Department of Mechanical Engineering and Center for
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More information