Introduction To CFD

Introduction of Computational Fluid Dynamics
Wangda Zuo
FAU Erlangen-Nrnberg
JASS 05, St. Petersburg
1. Basics of Computational Fluid Dynamics

1.1.
Concept of Computational Fluid Dynamics
Computational Fluid Dynamics (CFD) is the simulation of fluids engineering systems using
modeling (mathematical physical problem formulation) and numerical methods (discretization
methods, solvers, numerical parameters, and grid generations, etc.). The process is as figure 1.
Figure 1
Process of Computational Fluid Dynamics
Firstly, we have a fluid problem. To solve this problem, we should know the physical properties of
fluid by using Fluid Mechanics. Then we can use mathematical equations to describe these physical
properties. This is Navier-Stokes Equation and it is the governing equation of CFD. As the
Navier-Stokes Equation is analytical, human can understand it and solve them on a piece of paper.
But if we want to solve this equation by computer, we have to translate it to the discretized form.
The translators are numerical discretization methods, such as Finite Difference, Finite Element,
Finite Volume methods. Consequently, we also need to divide our whole problem domain into
many small parts because our discretization is based on them. Then, we can write programs to solve
them. The typical languages are Fortran and C. Normally the programs are run on workstations or
supercomputers. At the end, we can get our simulation results. We can compare and analyze the
simulation results with experiments and the real problem. If the results are not sufficient to solve
the problem, we have to repeat the process until find satisfied solution. This is the process of CFD.
1.2. Importance of Computational Fluid Dynamics

There are three methods in study of Fluid: theory analysis, experiment and simulation (CFD). As a
new method, CFD has many advantages compared to experiments. Please refer table 1.
Simulation (CFD)
Experiment
Cost
Cheap
Expensive
Time
Short
Long
Scale
Any
Small/Middle
Information
All
Measured Point
Repeatable
Yes
Some
Safety
Yes
Some Dangerous
Table 1
Comparison of Simulation and Experiment
1.3. Application of Computational Fluid Dynamics

As CFD has so many advantages, it is already generally used in industry such as aerospace,
automotive, biomedicine, chemical processing, heat ventilation air condition, hydraulics, power
generation, sports and marine etc.
1.4. Physics of Fluid

Fluid is liquid and gas. For example, water and air. Fluid has many important properties, such as
velocity, pressure, temperature, density and viscosity.
The density (1) of a fluid is its mass per unit volume. If the density of fluid is constant (or the
change is very small), we call the fluid is incompressible fluid. If the density of fluid is not constant,
we call the fluid is compressible fluid. Normally, we can treat water and air as incompressible fluid.
If the fluid is incompressible, we can simplify the equations for this type of fluid.
M
V
kg
m 3
(1)
The viscosity (2) is an internal property of a fluid that offers resistance to flow. For example, to stir
water is much easier than to stir honey because the viscosity of water is much smaller than honey.
Ns
= 3 = [Posie]
m
(2)
Table 2 shows the densities and viscosities of air, water and honey.
Substance
Air (18C)
Water (20C)
Honey (20C)
Density (kg/m3)
1.275
1000
1446
Viscosity (P)
1.82e-4
1.002e-2
190
2. 2. Navier-Stokes Equations
2.1. Conservation Law

Navier-Stokes equations are the governing equations of Computational Fluid Dynamics. It is based
on the conservation law of physical properties of fluid. The principle of conservational law is the
change of properties, for example mass, energy, and momentum, in an object is decided by the
input and output.
For example, the change of mass in the object is as follows
dM
= m& in m& out
dt
(3)
dM
=0
dt
(4)
M = const
(5)
& in m& out = 0 , we have

If m
Which means
2.2. Navier-Stokes Equation

Applying the mass, momentum and energy conservation, we can derive the continuity
equation, momentum equation and energy equation as follows.
Continuity Equation
D
U i
+
=0
Dt
xi
(7)
Momentum Equation
U j
P ij
+ U i
=
+ g j
3
x j xi {
t
xi
12
{
V
1424
3
{
I
U j
II
III
(8)
IV
Where
U j
ij =
xi
I:
II:
III:
IV:
V:
U i 2
U k
+ ij
x j 3
xk
(9)
Local change with time

Momentum convection
Surface force
Molecular-dependent momentum exchange (diffusion)
Mass force
Energy Equation
U j
T
T
U
2T
+ cU i
= P i + 2 ij
t
xi
xi
xi
xi
123
1424
3
123
123
123
I:
II:
III:
IV:
V:
II
III
IV
(10)
Local energy change with time

Convective term
Pressure work
Heat flux (diffusion)
Irreversible transfer of mechanical energy into heat
If the fluid id compressible, we can simplify the continuity equation and momentum equation as
follows.
Continuity Equation
U i
=0
xi
Momentum Equation
(11)
2U j
P
+ g j
=
+ U i
xi2
xi
x j
t
U j
U j
(12)
2.3. General Form of Navier-Stokes Equation

To simplify the Navier-Stokes equations, we can rewrite them as the general form.
( )
+
t
xi

= q
U i
x
i
(13)
When = 1, U j , T , we can respectively get continuity equation, momentum equation and energy
equation.
3. Finite Volume Method

The Navier-Stokes equations are analytical equations. Human can understand and solve them, but if
we want to solve them by computer, we have to transfer them into discretized form. This process is
discretization. The typical discretization methods are finite difference, finite element and finite
volume methods. Here we introduce finite volume method.
3.1. The Approach of Finite Volume Method

Integrate the general form of Navier-Stokes equation over a control volume and apply Gauss
Theory
dV = ni dS
(14)
We can get the integral form of Navier-Stoke equation

( )
ni dS = q dV
dV + U i
xi
t
S
V
(15)
To approximate the the volume integral, we can multiply the volume and the value at the center of
the control volume. For example, we have a 2D domain as fig 2. To approximate the mass and
momentum of control volume P, we have
m = dV pV ,
Vi
mu = i ui dV P u PV
(16)
Vi
To approximate the surface integral, for example pressure force, we have
Si
P dS Pk S k
k = n, s, e, w
(17)
Normally we store our variables at the center of control volume, so we need to interpolate them to
get Pk , which are located at the surface of control volume.
Typically, we have two types of interpolations, one is upwind interpolation, and the other one is
central interpolation.
UE
UP
Ue
Figure 2
2D Structured Grid Domain
Upwind Interpolation
U P
Ue =
U E
r r
if (U n ) e > 0
r r
if (U n ) e < 0
Central Interpolation
U e = U E e + U P (1 e )
e =
xe xP
xE xP
3.2. Conservation of Finite Volume Method

If we use finite difference and finite element approach to discretized Navier-Stokes equation, we
have to manually control the conservation of mass, momentum and energy. But with finite volume
method, we can easily find out that, if the Navier-Stokes equation is satisfied in every control
volume, it will automatically be satisfied for the whole domain. In another words, if the
conservation is satisfied in every control volume, it will be automatically satisfied in whole domain.
That is the reason why finite volume is preferred in computational fluid dynamics.
4. Grids
There are three types of grids: structured grids, unstructured grids and block structured grids.
The simplest one is structured grid (fig 3). This type of grids, all nodes have the same number of
elements around it. We can describe and store them easily. But this type of grid is only for the
simple domain.
If we have a complex domain, we can use unstructured grid. For example, fig 4 is an airfoil. The
structure of airfoil is very complex. The flow near the object is very important and complex, we
need very fine grid at this region. Far away from the airfoil, the flow is comparably simple, so we
can use coarse grid. Generally, unstructured grid is suitable for all geometries. It is very popular in
CFD. The disadvantage is that because the data structure is irregular, it is more difficult to describe
and store them.
Block structure grid is a compromising of structured and unstructured grid. The idea is, firstly,
divide the domain into several blocks, then use different structured grids in different blocks.
Figure 3
Structured Grids
Figure 4
Unstructured Grids
5. Boundary Conditions
To solve the equation system, we also need boundary conditions. The typical boundary conditions
in CFD are No-slip boundary condition, Axisymmetric boundary condition, Inlet, outlet boundary
condition and Periodic boundary condition.
For example, fig 5 is a pipe, the fluid flows from left to right. We can use inlet at left side, which
means we can set the velocity manually. At the right side, we use outlet boundary condition to keep
all the properties constant, which means all the gradients are zero.
At the wall of pipe, we can set the velocity to zero. This is no-slip boundary condition.
At the center of pipe, we can use axisymmetric boundary condition.
No-Slip B.C
Outlet B.C.
Inlet B.C.
Axisymmetric B.C.
Figure 5
Boundary Conditions of Pipe Flow

Introduction To CFD

Uploaded by

Copyright:

Available Formats

Introduction To CFD

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Introduction To CFD

Uploaded by

Copyright:

Available Formats

Introduction of Computational Fluid Dynamics

1. Basics of Computational Fluid Dynamics

Concept of Computational Fluid Dynamics

Process of Computational Fluid Dynamics

1.2. Importance of Computational Fluid Dynamics

Comparison of Simulation and Experiment

1.3. Application of Computational Fluid Dynamics

1.4. Physics of Fluid

2.1. Conservation Law

& in m& out = 0 , we have

2.2. Navier-Stokes Equation

Local change with time

Local energy change with time

2.3. General Form of Navier-Stokes Equation

3. Finite Volume Method

3.1. The Approach of Finite Volume Method

We can get the integral form of Navier-Stoke equation

To approximate the surface integral, for example pressure force, we have

2D Structured Grid Domain

3.2. Conservation of Finite Volume Method

Boundary Conditions of Pipe Flow

You might also like