Turbulence Modeling in Phoenics
Turbulence Modeling in Phoenics
Turbulence Modeling in Phoenics
• The eddies’ smallest scale lk ratio for Re of 400 and 20000 is of 19:1
Prediction Methods
l h = l/ReL3/4
Direct numerical simulation (DNS)
ui u j
ij eij
x j xi
9
Predicting the turbulent viscosity
• The turbulent viscosity can be estimated by algebraic models or by
solving a system of PDE.
• Algebraic models (mixing length and LVEL) employs expressions
similar to: T = C(dU/dy). They are computationally cheap but limited
to very simple flows (pipe and boundary layers)
• The 2nd category usually solves two equations to get T . So they are
called by two equation models: k-, k-, k-l, k- RNG, etc. For k-:
T = Ck2/ ; T derived by other 2 eq. models similarly follows k-.
• There is also the Reynolds stress model which instead of estimate T it
solves a system of six turbulent stress equations in order to get the
Reynolds turbulent stress tensor.
• The two equation models are largely applied to solve industrial
problems. Very often the constants are tuned to specific cases to
increase the accuracy.
10
Comparison of RANS turbulence models
11
Wall bounded flows
• The presence of a wall damps the velocity and the turbulence near the
wall region.
• To capture the wall effect the turbulence models have to be modified
near the wall.
• The embodied wall models on the turbulence models are grid
dependent.
• The success of turbulent flow simulation starts on a careful choice of
the grid spacing near the wall.
• This presentation focus on near wall grid spacing selection to a
successful simulation. To beginners this a potential area to make
mistakes
Wall bounded flows and velocity damping near the wall
y v*
y
u
u *
u u
u
v* where
u* W
is the eddy
velocity scale.
5 40 200
y v*
y
Wall
layers Inner Overlap Log Outter
Overlap: sensitivity to adverse pressure gradient
Phoenics turbulence models
• Algebraic models: constant turbulent viscosity, mixing length, lvel.
Lvel is the most popular, it can be applied starting from the inner layer
or from the log layer. It is an algebraic expression which fits the near
wall flow. These models are computationally cheap and delivers good
results for channel or boundary layer flows.
• Two equation models (standard) : refers to the k-e and variants. The
distance from the wall of the first volume center has to lay within the
log layer, 40 < y+ < 200.
• Two equation models (low Re) : turbulent low Reynolds flows usually
have the first node distance for 40 < y+ < 200 corresponding to a
significant percentage of the domain (> 15%). To avoid numerical
errors induced by the coarse grid it is preferred to starting integrating
the flow starting from the wall. The first node must lay at y+ ~ 1 and at
least 3 more volumes up to y+ < 5. This assures a smooth integration
along the inner, overlap and log layer.
How to get good results from turbulence models
• All models have limitations, one can not ask more than the model
can deliver.
• Algebraic models are the simplest and the most limited. Usually
work well for pipe flows and boundary layers.
• 2 equation models are more accurate than algebraic models but
their constants do not have universally, it is possible for certain types
of flow a better matching by tuning the constants.
• 2 equation models considers flow in equilibrium: all the
production of k is dissipated at the same point. Flow with sudden
changes in direction are prone to be out of equilibrium.
• In general, turbulent models demand a special care with the near
wall mesh. If you get this wrong your simulation will fail.
WKSP #1-Turbulent flow develop. 2D channel
• This WKSH deals with a turbulent flow development along a 2D
channel flow.
• The flow simulations will be done using Parabolic model (visit (1)
and (2) links to further information).
• Parabolic model apply only to one-way flows (no recirculation
zones present). It will used along the WKSH because it is far more
efficient than Elliptic model.
• Along the WKSH will be supplied specific hints to set up a problem
employing Parabolic model.
2D Channel dimensions & properties
y nwall (plate) (b) nooutlet Uniform GRID (c)
Inlet specification NX=1
NY=28
H/2=0.05m NZ=80
Center line 2H U
z (a) Re
L = 10m
S E T T I N G PA R A B O L I C M O D E L :
Model: choose parabolic (confined because is a channel)
Ground: Parabolic model visits only one slab at a time and stores only the last slab.
To provide full field storage go to VR-EDITOR box GRND and set IDSPA=1,
IDSPB=1 and IDSPC=80 (always equal to NZ), IDSPD=1
One may reproduce this figure using autoplot or get similar results
employing ‘Ploting variable’ in VRVIEWER.
2D Channel Results – W center-line velocity
overshoot
experimental
overshoot,
see link
STRS = w/
STRS fully developed = 0.192
STRS Colebrook-White = 0.195 (1.5% off)
ENUT
KE
EP
KE – Low Re
KE – Low Re
KE - Standard
KE - Standard
Workshop#1 – adjusting grids
• Based on the table below select between ‘standard’ or ‘low Re’ KE
models to cases where W = 1 and 50 m/s, also define a new NY. Do
these changes on the previous Q1.
• For reference, the previous case velocity is W = 10m/s or Re 105,
• STRS Cole is the w/ value estimated by Colebrook-White for a fully
developed flow.
• Y+ = 1 thru 100 are the estimated first node distance to the wall.
• Case W = 1m/s – ‘Low Re’ because the 1st volume for y+ = 40 would
be at 7mm from the wall which is nearly 14% of the Y length. One can
use a uniform grid or a more economical 2 region grid, the near wall
grid extends up to y+ = 40 (~7mm)
• Case W = 50m/s – Standard,
Workshop#2 – Flow in backward face step
• A simple backward step flow challenges the 2 equation turbulence
models because it has a shear layer, a recirculation zone, null wall
shear stress point and a redeveloping boundary layer!
x (m) x/H
k-e
Chen-Kim k-e
RNG k-e
k-omega
LVEL
Workshop#2 – Flow in backward face step
values
(a) For parabolic model the main flow direction has to be aligned with the Z axis.
(b) Parabolic flow has w component always along z > 0 direction, if happens to have
w along z<0 it will fail. Due this feature it does not need specification of an outlet.
(c) Parabolic model visits only one slab at a time and stores only the last slab. To
provide full field storage go to VR-EDITOR box GRND and set IDSPA=1, IDSPB=1
and IDSPC=80 (always equal to NZ), IDSPD=1