CFD Basics
CFD Basics
Ladzinski Horvat
Collaboration – Innovation – Productivity - Quality
THE INTERNATIONAL ASSOCIATION
FOR THE ENGINEERING ANALYSIS
COMMUNITY
Matthew Ladzinski
NAFEMS North America
Silver Sponsors
For more information, please visit: www.nafems.org Collaboration – Innovation – Productivity - Quality
Fire Modelling
in Computational Fluid
Dynamics (CFD)
13
Introduction IFS
Fire modelling is an area of computational modelling which aims to
predict fire behaviour in different environmental conditions.
Therefore, these computational models need to take into account
fluid dynamics, combustion and heat transfer processes.
The complexity of the fire modelling arises from significantly
different time scales of the modelled processes. Also, not
completely understood physics and chemistry of fire adds the
uncertainty to the modelling process.
14
IFS
15
Transport equations IFS
Eulerian and Lagrangian description
17
Transport equations IFS
18
Transport equations IFS
The following physical laws and terms also need to be
included
}
- Newton's viscosity law diffusive terms -
- Fourier's law of heat conduction flux is a linear
function of a gradient
- Fick's law of mass transfer
- Sources and sinks due to thermal radiation, chemical
reactions etc.
19
Transport equations IFS
Transport of mass and composition
∂ t ρ + ∂ i (vi ρ) = M ∂ t (ρξ j ) + ∂ i (vi ρξ j ) = ∂ i (ρD∂ i ξ j ) + M j
Transport of momentum
Transport of energy
volumetric
term
change diffusion
flux difference
in a control vol.
(convection) 20
Transport equations IFS
Lagrangian formulation is simpler
r r
- particle location equation dx dt = u
23
Averaging and simplification of
transport equations
IFS
The presented set of transport equations is analytically
unsolvable for majority of cases
ρ, vi , p, h
w χ,τ
25
Averaging and simplification of
transport equations
IFS
Kolmogorov scale is (for most cases) too small to be
captured with a numerical grid
26
Averaging and simplification of
transport equations
IFS
Transport equation variables can be decomposed onto
a filtered (averaged) part and a residual (fluctuation)
ρ = ρ + ρ' p = p + p' ~ ~
vi = ρvi ρ + v*i = v~i + v*i ξi = ξi + ξ* hi = h + h*
29
Averaging and simplification of
transport equations
IFS
a) b)
CFD Modelling
32
Turbulence models IFS
laminar flow
transitional flow
turbulent flow
34
Turbulence models IFS
In general, 2 kinds of averaging (filtering) exist, which
leads to 2 families of turbulence models:
- filtering over a spatial interval → Large Eddy Simulation
(LES) models
- filtering over a time interval → Reynolds Averaged
Navier-Stokes (RANS) models: k-epsilon model, SST
model, Reynolds Stress models etc
For RANS models, size of the averaging time interval is
not known or given (statistical average of experimental
data)
For LES models, size of the filter or the spatial averaging
interval is a basic input parameter (in most cases, it is
equal to grid spacing)
35
Turbulence models IFS
Reynolds Averaged Navier-Stokes (RANS) models
For two-equation models (e.g. k-epsilon, k-omega or
SST), 2 additional transport equations need to be solved:
or
- for frequency of turbulent fluctuations ω~ ε k
37
Turbulence models IFS
Large Eddy Simulation (LES) models
- Large Eddy Simulation (LES) models are based on spatial
filtering (averaging)
- many different forms of the filter exist, but the most
common is "top hat" filter (simple geometrical averaging)
- size of the filter is based on a grid node spacing
Basic assumption of LES methodology:
Size of the used filter is so small that the averaged flow
structures do no influence large structures, which do
contain most of the energy.
These small structures are being deformed, disintegrated
onto even smaller structures until they do not dissipate
due to viscosity (kinetic energy → thermal energy). 38
Turbulence models IFS
Large Eddy Simulation (LES) models
- eddy (turbulent) viscosity is defined as
µ t ~ ρ l 4 / 3ε1 / 3 where l ~ Cs ∆
grid spacing
- using the definition of turbulence (subgrid) stresses
2 2
ρ Π ji − ρ k δ ji = −2µ t S ij + µ t (∂ l ~
vl ) δ ji
3 3
µt ~
ρΩ j = − ∂ jh
and turbulence fluxes Prt
40
Turbulence models IFS
Comparison of turbulence models
k-e (Gr = 10 10 ) k-e (Gr = 10 10 )
k-e &&B (Gr = 10 1 0 ) k-e &&B (Gr = 10 10 )
8 R&G k-e (Gr = 1 0 10 ) R&G k-e (Gr = 1 0 10 )
7 S S T (Gr = 10 1 0 ) S S T (Gr = 1 0 10 )
10 10
6 S S G (Gr = 10 ) S S G (Gr = 10 )
10 10
5 LES (Gr = 10 ) LES (Gr = 10 )
Rous e e t al. (195 2) Rous e e t al. (195 2)
4 S habbir and Ge o rge (1 994) S habbir and Ge org e (1994)
-1
10
3
b c (R 5 /F20 )1/3
w c (R/F 0 )1/3
1 10 -2
25 50 75 100 25 50 75 100
z/R z/R
a) b)
Buoyant flow over a heat source: a) velocity, b) temperature*
41
Combustion models IFS
Grinstein,
Chen et al., 1988 Kailasanath, 1992
42
Combustion models IFS
Combustion can be modelled with heat sources
- information on chemical composition is lost
- thermal loading is usually under-estimated
43
Combustion models IFS
Modelling of composition requires solving n-1 transport
equations for mixture components – mass or molar
(volume) fractions
µ µ ~
( ) (
~
∂t ρ ξ j + ∂i ρ~
~
)
vi ξ j = ∂ i + t ∂ i ξ j + M j
Sc Sct
or
, ,
∑ν' I ↔ ∑ν" I
I
I = A,B ,C.... I = A,B ,C....
I
44
Combustion models IFS
Reaction source term is defined as
M j = W j (ν" j −ν' j )R
or for multiple M j = W j ∑ (ν" k , j −ν ' k , j )Rk
reactions k
45
Combustion models IFS
Constant burning velocity sL
ρc
sF = sL speed of flame front propagation is
ρh
larger due to expansion
outflow of comb.
products
fireball
49
Combustion models IFS
Flamelet model
- describes interaction of reaction kinetics with turbulent
structures for a fast reaction (high Damköhler number)
- basic assumption is that combustion is taking place in
thin sheets - flamelets
- turbulent flame is an ensemble of laminar flamelets
- the model gives a detailed picture of the chemical
composition - resolution of small length and time scales
of the flow is not needed
- the model is also known as "Mixed-is-burnt" - large
difference between various implementations of the
model
50
Combustion models IFS
Flamelet model
- it is based on definition of a mixture fraction
Z kg/s fuel
A 1 kg/s mixture
mixing
1-Z kg/s oxidiser M
proces
B
βM − βB
Z β A + (1 − Z )β B = β M or Z= where β = ξ f − ξo i
β A − βB
dI ν (Ω ) K
= −(K aν + K sν )I ν (Ω ) + K aν I eν + sν ∫ I (Ω' )P (Ω' → Ω)dΩ'
sν ν
ds 4π 4π
2πν 2 n 2 hν
Eν (T ) = πI ν (T ) = 2 [Wm-2 Hz -1 ]
c exp(hν k BT ) − 1
ν - frequency
c - speed of light
n - refraction index
h - Planck's constant
kB - Boltzmann's constant
∞
, j = ∫ I ν (r , s ) cos ϕ j cos θ dΩ
qνrad
4π
57
Thermal radiation IFS
Flashover simulation
wooden targets
propane burner 58
IFS
Conclusions
59
Conclusions IFS
The webinar gave a short (but demanding) overview of
fluid mechanics and heat transfer theory that is relevant
for fire simulations
All current commercial CFD software packages (ANSYS-
CFX, ANSYS-Fluent, Star-CD, Flow3D, CFDRC, AVL Fire)
contain most of the shown models and methods:
- they are based on the finite volume or the finite element
method and they use transport equations in their
conservative form
- numerical grid is unstructured for greater geometrical
flexibility
- open-source computational packages exist and are freely
accessible (FDS, OpenFoam, SmartFire, Sophie)
60
Contact information IFS
Andrej Horvat
Intelligent Fluid Solutions Ltd.
99 Milton Park, Abingdon, Oxfordshire, OX14 4RY, UK
Tel.: +44 (0)1235 841 505
Fax: +44 (0)1235 854 001
Mobile: +44 (0)78 33 55 63 73
Skype: a.horvat
E-mail: [email protected]
Web: www.intelligentfluidsolutions.co.uk
61
Questions
Website: www.nafems.org
Collaboration – Innovation – Productivity - Quality
THE INTERNATIONAL ASSOCIATION
FOR THE ENGINEERING ANALYSIS COMMUNITY
Thank you!