Download Governing Equations - Aerodynamics - Lecture Notes and more Study notes Engineering Dynamics in PDF only on Docsity!
II. Governing Equations
In aerodynamics, or fluid mechanics, there are six properties of the flow an
engineer is usually interested in – pressure p, density ρ, the three velocity components
(u,v,w), and the temperature T. For gases and mixtures of gases (e.g. air) the equation of
state links p, ρ and T by:
Equation of State:
p = ρRT
(2.1)
We need to come up five additional equations linking the 6 properties. These five
equations are PDEs and turn out to be:
a) Conservation of Mass or Continuity
b) Conservation of u- momentum
c) Conservation of v-momentum
d) Conservation of w-momentum
e) Conservation of energy
These equations may be derived using a Lagrangean approach, or an Eulerian
approach.
In the Lagrangean approach, we follow a fixed set of fluid particles (e.g. a cloud,
a tornado, tip vortices from an aircraft) and write down equations governing their motion.
This is somewhat like tracking satellites and missiles in space, using equations to
describe their position in space and the forces acting on them.
In the Eulerian approach, we look at a (usually) fixed or (sometimes) moving
volume in space surrounded by permeable boundaries. We develop equations describing
what happens to the fluid inside the control volume as new fluid enters and old fluid
particles leave. Eulerian approach is the preferred approach in most fluid dynamics
applications. This is what we will follow in our derivations.
Conservation of Mass (Continuity):
The conservation of mass stems from the principle that mass can not be created or
destroyed inside the control volume. Obviously, we are situations (e.g. nuclear reactions)
involving the conversion of mass into energy.
Let V be a control volume, a balloon like shape in space. We will assume that it, and its
surface S remain fixed in space. The surface is permeable so that fluid can freely enter in
and leave. The continuity equation says
The time rate of change of Mass within the control volume V =
Rate at which mass enters V through the boundary S
We can assume that the control volume V is made of several infinitely small
(infinitesimal) volume elements dV. The mass of the fluid inside each of these elements
is ρdV, where ρ is the fluid density. The density is free to change from point to point,
from one sub-element dV to another within V. Thus,
Total mass within V = ρdV
V
∫∫∫
In the above equation, the three integral signs simply indicate that we are doing a three-
dimensional integration, or volume integration. The subscript V says that this integration
takes place inside V. This subscript may seem redundant at this time, but may be used to
distinguish three or more control volumes V1, V2, V3 etc. from one another in some
problems.
Then,
Time rate of change of mass within the control volume V =
d
dt
ρ dV V
∫∫∫
In calculus, the order of integration or differentiation may be interchanged so that
these operations do not interfere with each other. Since our control volume is fixed in
space, the limits of the volume integral are not functions of time, and there is no
interaction between the two operations. Thus,
Time rate of change of mass within the control volume V =
∂ρ
∂ t
dV V
∫∫∫
Notice that we are using the partial derivative inside the integral since ρ is a
function of (x,y,z,t) and we are only interested in its variation inside each sub-element dV
with respect to time, while the spatial location (x,y,z) of the sub-element dV remains
fixed.
We next turn our attention to the amount of fluid that enters V through the surface
S. For this purpose, we assume that the surface S is made of many quilt-like infinitesimal
patches dS. At the center of each element is a unit normal vector (i.e. a vector of length
unity, normal to the surface)
n. By common convention, this normal is pointing away
from the surface dS. The normal component of fluid velocity pointing towards the control
volume (entering the control volume) is given by − •
V n. Notice the negative sign. We
are interested in the component of velocity pointing towards the control volume, not
away from it.
We can compute the rate at which mass enters the control volume through dS as
the product of density times normal velocity times area. Thus,
Rate at which mass enters the control voume through dS = - ρV ndS
Then,
When this operator operates on a function g(x,y,z,t) the result is called the gradient of g,
or simply “del g”. The gradient of a scalar in a Cartesian coordinate system is:
(^) ∇ g i + +
g j
g k
g
x y z
Notice that the operation produces a vector.
Since the “del”operator is a vector, it can act on other vectors giving a dot product or
cross product. The dot product between
∇ and
F is called a divergence of
F , and the
cross product between
∇ and F is called a curl of
F. These operations are performed as
follows:
Let
Then,
Divergence of F = F =
F
x
F
y
F
z
Curl of F = F =
i
1 2 3
^
F F i F j F k
j k
x y z F F F
∇ ×
1 2 3
1 2 3
Assignment: We have shown these operations in a Cartesian coordinate system. Read the
relevant sections of Chapter II to find out how the divergence and curl operations are
carried out in the cylindrical coordinate system.
We return back to continuity equation (2.4) and the divergence theorem (2.5) from our
diversion into the land of vector calculus, a mystery for some and a nightmare for many.
Applying equation (2.5) to the surface integral in the continuity equation (2.4) we get:
( )
t
dV V dV V V
∫∫∫ +^ ∫∫∫∇ •^ =
or,
( )
t
V dV V
∫∫∫ (^) =
Consider the above volume integral. It must hold for any arbitrarily shaped
control volume V, at any instance in time for all flows. The only way this can be true is if
the integrand is zero. Therefore,
( )
t
+ ∇ • V =
Equation (2.6) is called the differential form (PDE form) of continuity equation.
For steady flows, the time derivative vanishes. For incompressible flows, r is a constant.
Thus, continuity equation for incompressible flows becomes:
V
or
u x
v y
w z
Body Forces: These are forces acting on every fluid particle within V. Examples of body
forces include gravity, electrostatic forces, and magnetic forces. Let us assume the
symbol a (^) X represents the x- component of all these effects acting on the fluid particles
within V. The quantity a (^) X may vary with x,y,z or t. Then,
Body force along the x - direction acting on all paticles wthin dV = dV
Body force along the x - direction acting on all paticles wthin V = V
a
a dV
x
∫∫∫ x
Surface Pressure Forces: The control volume V is surrounded by the surface S. Pressure
from surrounding fluid (or solid) acts on this surface S. Pressure forces are always
directed towards the fluid within, and is always normal to the surface. The pressure
forces acting on S may therefore be written as − (^) ∫∫ pndS
S
. Notice the negative sign. It is
there because the normal vector
n is pointing outwards, whereas the pressure forces are
acting inwards. The x- component of these pressure forces is found by performing a dot
product of this expression with
i (^).
Thus,
X − component of Surface Forces acting on V= − (^) ∫∫ pn • i dS
Surface Viscous forces: The surrounding fluid can exert an additional type of force on the
surface S, called a “viscous” force. This force may have a component normal to the
surface S, and components tangential to the surface. We will study the effects of viscosity
in more detail later. For now, we will call this contribution Fx-Viscous.
We can finally turn our attention to the rate at which the u- momentum is brought
into V though the surface S. Let dS be an infinitesimal element on S. Then rate at which
u-momentum enters the control volume through dS is − ρ uV • ndS
. Summing up
contributions over the entire surface S, we get:
Rate at which u - momentum enters V through S = - S
ρ uV ndS
∫∫^ •
We have independently found all the terms on the left and right hand side of the
u-momentum equation. Bringing over all the terms, except the viscous forces and the
body-forces, to the left-hand side, we get:
( ) ( )
∂ ρ
∂
ρ υ
u
t
dV pi uV ndS a (^) x Body dV Fx Viscous V S V
∫∫∫ +^ ∫∫ +^ •^ =^ ∫∫∫ − + −
The above equation is the conservation of u- momentum equation in integral
form.
We can similarly derive the v- and w- momentum equations. Or simply replace all
places where u appears with v,
i with j , etc. Then we get
Conservation of v- momentum Equation in Integral Form:
( ) ( )
v
t
dV pj vV ndS a (^) y Body dV Fy Viscous V S V
∫∫∫ +^ ∫∫ +^ •^ =^ ∫∫∫ − + −
Conservation of w- momentum Equation in Integral Form:
( ) ( )
∂ ρ
∂
ρ υ
w
t
dV pk wV ndS a (^) z Body dV Fz Viscous V S V
∫∫∫ +^ ∫∫ +^ •^ =^ ∫∫∫ − + −
In this course, we will henceforth neglect the viscous forces acting on the fluid.
Needless to say, neglecting this effect is an approximation, and all the resulting solutions
only approximate real flows. Viscous forces are small except in regions of high velocity
gradient. In many flows, velocity varies slowly. Neglecting viscous effects turns out to be
an acceptable compromise.
We will also neglect body forces such as gravity, electrical and magnetic effects.
These effects are small in many applications compared to the rest of the terms in the
momentum equations. There are applications such as formation of weather systems where
gravity does play a part. In ionized flows, the ions are attracted and deflected by electrical
and magnetic forces. This is what causes the Aurora Borealis phenomenon! We neglect
these forces because they are small compared to other forces acting on the control
volume.
The effect of these approximations is that the right hand sides of the integral
forms of u- v- and w- momentum equations (2.8, 2.9, 2.10) disappear, leaving us only
with the left-hand side terms.
( ) ( )
( ) ( ) ( ) ( ) ( )
( ) ( ) x
p pi z
k y
j x
pi i
uw z
uv y
u x
u i uvj uwk z
k y
j x
uV i
uV uui vj wk ui uvj uwk
2 2
2
With these operations, equation (2.15) becomes:
( ) ( ) ( ) ( ) 0
2
∂
x
p uw z
uv y
u x
u t
Equation (2.17) is the PDE form of the u- momentum equation, while equation (2.11) in
the integral form.
We can similarly convert the v- and w- momentum equations into PDE forms. Sparing
you the details, here are the final forms:
v- Momentum equation in PDE form:
( ) ( ) ( ) ( ) 0
2
∂
y
p vw z
v y
uv x
v t
w- Momentum Equation in PDE form:
( ) ( ) ( ) ( ) 0
2
∂
z
p w z
vw y
uw x
w t
ρ ρ ρ ρ
You need not memorize these equations. You should, however, be able to recognize them
and identify them if you are given the expressions.
Equations (2.17) - (2.19) may be simplified further. Consider the various terms in
equation (2.17). The derivatives are products of several terms. We can use the product
rule. That is,
( )
( ) ( ) ( ) ( )
( ) ( ) ( )
( ) ( ) ( )
z
w u z
u w z
uw
y
v u y
u v y
uv
x
u u x
u u x
uu
x
u
t
u t
u u t
2
Using these expansions (2.20) in equation (2.17), we get:
( ) ( ) ( ) =^0
z
w
y
v
x
u
t
u x
p
z
u w y
u v x
u u t
u ρ ρ ρ ρ
The term inside the square bracket is the continuity equation and vanishes. The u-
momentum equation in PDE form thus becomes:
u-Momentum Equation: =^0 ∂
x
p
z
u w y
u v x
u u t
u
By a similar process, the v- and w- momentum equations become:
v-Momentum Equation: = 0 ∂
y
p
z
v w y
v v x
v u t
v
w-Momentum equation: =^0 ∂
z
p
z
w w y
w v x
w u t
w
Note that equations (2.22)-(2.24) may be written as:
z
w y
v x
u Dt t
D
where
z
p
Dt
Dw
y
p
Dt
Dv
x
p
Dt
Du
