Turbulence Models: Dr.N.V.Mahalakshmi Asst Professor, Department of Mechanical Engineering
Turbulence Models: Dr.N.V.Mahalakshmi Asst Professor, Department of Mechanical Engineering
Turbulence Models: Dr.N.V.Mahalakshmi Asst Professor, Department of Mechanical Engineering
Dr.N.V.Mahalakshmi
Asst Professor,
Department of Mechanical Engineering
1
Turbulence models
• A turbulence model is a computational procedure to close the
system of mean flow equations.
• For most engineering applications it is unnecessary to resolve the
details of the turbulent fluctuations.
• Turbulence models allow the calculation of the mean flow without
first calculating the full time-dependent flow field.
• We only need to know how turbulence affected the mean flow.
• In particular we need expressions for the Reynolds stresses.
• For a turbulence model to be useful it:
– must have wide applicability,
– be accurate,
– simple,
– and economical to run.
2
Common turbulence models
• Classical models. Based on Reynolds Averaged Navier-Stokes
(RANS) equations (time averaged):
– 1. Zero equation model: mixing length model.
– 2. One equation model:
– 3. Two equation models: k-ε style models (standard, RNG,
realizable), k-ω model, and ASM.
– 4. Seven equation model: Reynolds stress model.
• The number of equations denotes the number of additional PDEs
that are being solved.
• Large eddy simulation. Based on space-filtered equations. Time
dependent calculations are performed. Large eddies are explicitly
calculated. For small eddies, their effect on the flow pattern is
taken into account with a “subgrid model” of which many styles are
available.
3
Boussinesq hypothesis
• Many turbulence models are based upon the Boussinesq hypothesis.
– It was experimentally observed that turbulence decays unless there is
shear in isothermal incompressible flows.
– Turbulence was found to increase as the mean rate of deformation
increases.
– Boussinesq proposed in 1877 that the Reynolds stresses could be linked
to the mean rate of deformation.
• Using the suffix notation where i, j, and k denote the x-, y-, and z-
directions respectively, viscous stresses are given by:
∂ui ∂u j
• Similarly, link Reynolds eij = µ to the
τ ij = µstresses rate of deformation:
+ mean
∂x j ∂xi
∂U i ∂U j
τ ij = − ρ ui ' u j ' = µ t +
∂x j ∂xi
4
Turbulent viscosity
• A new quantity appears: the turbulent viscosity µ t.
∂U i ∂U j
τ ij = − ρ ui ' u j ' = µ t +
∂x j ∂xi
5
Turbulent Schmidt number
• The turbulent viscosity is used to close the momentum equations.
• We can use a similar assumption for the turbulent fluctuation terms that appear
in the scalar transport equations.
• For a scalar property φ (t) = Φ + ϕ ’(t):
∂Φ
− ρui 'ϕ ' = Γt
• Here Γ t is the turbulent diffusivity.
∂xi
• The turbulent diffusivity is calculated from the turbulent viscosity, using a
model constant called the turbulent Schmidt number (AKA Prandtl number) σ t:
µ
• Experiments have shown that theσturbulent
= t
Schmidt number is nearly constant
with typical values between 0.7 and 1. Γt
t
6
Predicting the turbulent viscosity
• The following models can be used to predict the turbulent
viscosity:
– Mixing length model.
– Standard k-ε model.
– Realizable k-ε model.
– k-ω model.
• We will discuss these one by one.
7
Mixing length model
• On dimensional grounds one can express the kinematic turbulent
viscosity as the product of a velocity scale and a length scale:
ν t (m 2 / s ) ∝ ϑ (m / s ) (m)
9
One-equation model
• In the one-equation model or Norris and Reynolds the eddy
viscosity concept is used as well, however with a difference
function fµ. This model calculates only the turbulent kinetic energy
K from a transport equation, whereas the dissipation rate is now
obtained from
– Є = k3/2 / L (1 + CЄ / (k1/2 L/v) )
– Where L is prescribed algebraically in the form of a ramp distribution
like in the mixing – length approach
– L = CD min(ky,λ )
• is the boundary layer thickness.
10
The k-ε model
• The k-ε model focuses on the mechanisms that affect the turbulent kinetic energy (per unit mass) k.
• The instantaneous kinetic energy k(t) of a turbulent flow is the sum of mean kinetic energy K and
turbulent kinetic energy k:
K = 12 U 2 + V 2 + W 2
1
• ε is the dissipation rate of k. k = 2 u ' + v' +
2 2
w'2
• If k and ε are known, we can model the turbulent viscosity as:
k (t ) = K + k
k 3/ 2 k2
νt ∝ ϑ ∝ k 1/ 2
=
ε ε
11
Mean flow kinetic energy K
• The equation for the mean kinetic energy is as follows:
∂ ( ρK )
+ div( ρKU ) = div( − PU + 2 µUEij − ρUui ' u j ') − 2 µEij .Eij − ( − ρ ui ' u j '.Eij )
∂t
(I ) ( II ) ( III ) ( IV ) (V ) (VI ) (VII )
• Here Eij is the mean rate of deformation tensor.
• This equation can be read as:
– (I) the rate of change of K, plus
– (II) transport of K by convection, equals
– (III) transport of K by pressure, plus
– (IV) transport of K by viscous stresses, plus
– (V) transport of K by Reynolds stresses, minus
– (VI) rate of dissipation of K, minus
– (VII) turbulence production.
12
Turbulent kinetic energy k
• The equation for the turbulent kinetic energy k is as follows:
∂ ( ρk )
+ div ( ρkU) = div( − p ' u' + 2 µ u' eij ' − ρ 12 ui '.ui ' u j ') − 2 µ eij '.eij ' + (− ρ ui ' u j '.Eij )
∂t
(I ) ( II ) ( III ) ( IV ) (V ) (VI ) (VII )
• Here eij ’ is fluctuating component of rate of deformation tensor.
• This equation can be read as:
– (I) the rate of change of k, plus
– (II) transport of k by convection, equals
– (III) transport of k by pressure, plus
– (IV) transport of k by viscous stresses, plus
– (V) transport of k by Reynolds stresses, minus
– (VI) rate of dissipation of k, plus
– (VII) turbulence production.
13
Model equation for k
• The equation for k contains additional turbulent fluctuation terms,
that are unknown. Again using the Boussinesq assumption, these
fluctuation terms can be linked to the mean flow.
• The following (simplified) model equation for k is commonly used.
∂ ( ρk ) µ
+ div( ρkU) = div t grad k + 2 µ t Eij .Eij − ρ ε
∂t σ k
Rate of Rate of
Convective Rate of
increase Diffusive destruction
transport production
transport
14
Turbulent dissipation
• The equations look quite similar.
• However, the k equation mainly contains primed quantities,
indicating that changes in k are mainly governed by turbulent
interactions.
• Furthermore, term (VII) is equal in both equations. But it is
actually negative in the K equation (destruction) and positive in
the k equation: energy transfers from the mean flow to the
turbulence.
• The viscous dissipation term (VI) in the k equation − 2 µ eij '.eij '
describes the dissipation of k because of the work done by the
smallest eddies against the viscous stresses.
• We can now define the rate of dissipation per unit mass ε as:
ε = 2ν eij '.eij '
15
Dissipation rate - analytical equation
• The analytical equation for ε is shown below. Because of the
many unknown higher order terms, this equation can not be
solved, and simplified model equations need to be derived.
16
Model equation for ε
• A model equation for ε is derived by multiplying the k equation by (ε/k) and
introducing model constants.
• The following (simplified) model equation for ε is commonly used.
∂ ( ρ ε) µt ε ε2
+ div( ρ U
ε ) = div grad ε + C1ε 2 µ t Eij .Eij − C2ε ρ
∂t σ ε k k
17
Calculating the Reynolds stresses from k and ε
• The turbulent viscosity is calculated from:
k2
µt = Cµ C µ = 0.09
ε
• The Reynolds stresses are then calculated as follows:
∂U i ∂U j 2 2
− ρ ui ' u j ' = µ t + − ρkδ ij = 2 µ t Eij − ρkδ ij
∂x j ∂xi 3 3
δ ij = 1 if i = j and δ ij = 0 if i ≠ j
18
k-ε model discussion
• Advantages:
– Relatively simple to implement.
– Leads to stable calculations that converge relatively easily.
– Reasonable predictions for many flows.
• Disadvantages:
– Poor predictions for:
• swirling and rotating flows,
• flows with strong separation,
• axisymmetric jets,
• certain unconfined flows, and
• fully developed flows in non-circular ducts.
– Valid only for fully turbulent flows.
– Simplistic ε equation.
19
k-ω model
• This is another two equation model. In this model ω is an inverse
time scale that is associated with the turbulence.
• This model solves two additional PDEs:
– A modified version of the k equation used in the k-ε model.
– A transport equation for ω.
• The turbulent viscosity is then calculated as follows:
k
µt = ρ
ω
• Its numerical behavior is similar to that of the k-ε models.
• It suffers from some of the same drawbacks, such as the assumption
that μt is isotropic.
20
Algebraic stress model
• The same k and ε equations are solved as with the standard k-ε model.
• However, the Boussinesq assumption is not used.
• The full Reynolds stress equations are first derived, and then some
simplifying assumptions are made that allow the derivation of algebraic
equations for the Reynolds stresses.
• Thus fewer PDEs have to be solved than with the full RSM and it is much
easier to implement.
• The algebraic equations themselves are not very stable, however, and
computer time is significantly more than with the standard k-ε model.
• This model was used in the 1980s and early 1990s. Research continues
but this model is rarely used in industry anymore now that most
commercial CFD codes have full RSM implementations available.
21
Reynolds stress model
• RSM closes the Reynolds-Averaged Navier-Stokes equations by solving
additional transport equations for the six independent Reynolds stresses.
– Transport equations derived by Reynolds averaging the product of the
momentum equations with a fluctuating property.
– Closure also requires one equation for turbulent dissipation.
– Isotropic eddy viscosity assumption is avoided.
• Resulting equations contain terms that need to be modeled.
• RSM is good for accurately predicting complex flows.
– Accounts for streamline curvature, swirl, rotation and high strain rates.
• Cyclone flows, swirling combustor flows.
• Rotating flow passages, secondary flows.
• Flows involving separation.
22
Reynolds stress transport equation
• The exact equation for the transport of the Reynolds stress Rij :
DRij
= Pij + Dij − ε ij + Π ij + Ω ij
Dt
• This equation can be read as:
– rate of change of R = u ' uplus
ij i j'
– transport of Rij by convection, equals
– rate of production Pij , plus
– transport by diffusion Dij , minus
– rate of dissipation εij , plus
– transport due to turbulent pressure-strain interactions πij , plus
– transport due to rotation Ωij .
• This equation describes six partial differential equations, one for the
transport of each of the six independent Reynolds stresses.
23
Reynolds stress transport equation
• The various terms are modeled as follows:
– Production Pij is retained in its exact form.
– Diffusive transport Dij is modeled using a gradient diffusion assumption.
– The dissipation εij , is related to ε as calculated from the standard ε
equation, although more advanced ε models are available also.
– Pressure strain interactions πij , are very important. These include
pressure fluctuations due to eddies interacting with each other, and due
to interactions between eddies and regions of the flow with a different
mean velocity. The overall effect is to make the normal stresses more
isotropic and to decrease shear stresses. It does not change the total
turbulent kinetic energy. This is a difficult to model term, and various
models are available. Common is the Launder model. Improved, non-
equilibrium models are available also.
– Transport due to rotation Ωij is retained in its exact form.
24