Me305 Chapter 3
Me305 Chapter 3
Me305 Chapter 3
Chapter 3
Introduction to Fluid Dynamics
1
Chapter 3 Introduction to Fluid Dynamics
(a) Material (Lagrangian) Description: identified fluid particles are followed in the course of time
and the variation of their properties are determined.
Path of the fluid particle • As the time passes it moves in the flow field and
rP(0) is located at point B at time=t.
B(xB, yB, zB) • The path followed by this fluid particle from
rP(t)
point A to B can be given by specifying the
position vectors rP as a function of time.
y
x
r r r r
Position vector of P: rP (t) = x P (t) i + y P (t) j + z P (t) k
r
r d rP (t) dx P (t) r dy P (t) r dzP (t) r r r r
Velocity of P: VP (t) = = i + j + k = uP (t ) i + v P (t) j + w P (t)k
dt dt dt dt
2
Material (Lagrangian) Description (cont’d)
• In the Lagrangian description, coordinates of the fluid particle (xP, yP, zP) are dependent variables.
• They depend on time. They are used to locate the particle when the time is known.
• Fluid properties can be expressed as functions of xP, yP, zP and t, such as
r r
VP = VP [ x P ( t ), y P ( t ), z P ( t ), t ]
ρP = ρP [ x P ( t ), y P ( t ), zP ( t ), t ]
etc.
t = t0 t = t1 t = t2
3
Methods for Describing the Fluid Motion (cont’d)
(b) Spatial (Eulerian) Description: attention is focused at fixed points of the flow field and
properties of the fluid particles passing through these points are observed.
• In general, the properties in a flow field can change from point to point and from time to time.
r r
V = V ( x , y , z, t ) ρ = ρ ( x , y , z, t ) etc.
Control volume (open system): is a fixed space with respect to a given coordinate system.
• Fluid particles can enter and leave the control volume through its boundaries, called control surface.
t = t0 t = t1 t = t2
• A flow filed at a certain time can be described by both Lagrangian and Eulerian descriptions,
therefore there must be a relationship between them.
• If at time t = 0, the fluid particle P is located at point A,
r r r
VP = VP ( x P ( t = 0), y P ( t = 0), z P ( t = 0), t = 0 ) = V ( x A , y A , z A , t = 0)
Differential Integral
5
Relations Between the Material and Spatial Descriptions in the Differential Formulation
• One can relate the rate of change of any property of a fluid particle in the material description, to
the rate of chage of the same property in the spatial description at the point of the flow field where
the particle is located.
∂N ∂N ∂N ∂N
• Total differential of this property is DN = dx + dy + dz + dt
∂x ∂y ∂z ∂t
DN ∂ N dx ∂ N dy ∂ N dz ∂N
• Divide by dt = + + +
dt ∂ x dt ∂ y dt ∂ z dt ∂t
DN ∂N ∂N ∂N ∂N
• Use velocity components u, v and w = + u + v + w
dt ∂t ∂x ∂y ∂z
DN ∂N r r
• Put into vectorial notation = + (V ⋅ ∇ ) N
dt ∂t
• D/dt represents the rate of change of a property in the material description. It is called the total
derivative, substantial derivative, material derivative or the derivative following the fluid.
• ∂ /∂t represents the rate of change of a property in the spatial description. It is called the partial
derivative, local derivative or the spatial derivative.
• (V.∇) N represents the change of the property N from one point of the flow field to the other at
a certain time. It is called the convective derivative (read details from the book).
6
Relations Between the Material and Spatial Descriptions in the Integral Formulation
(Reynold’s Transport Theorem)
• To relate the change of any extensive property N of a closed system to the variations of the same
property associated with the control volume, consider the case shown below.
• Initially (at time t0), the system boundary coincides with the boundary of the control volume (CV).
• After a time interval of ∆t (at time t0+∆t), the control volume remains in its original position,
however the system moves to a new position.
• The mass in region A enters into the CV, while the mass in region C leaves the CV.
• Details of the derivation of Reynold’s Transport Theorem will be given during the class.
7
3.5 Classification of Fluid Flow
Steady and Unsteady Flows:
• The flow is referred as steady when the properties of the flow field are independent of time.
• In a steady flow, properties of the flow field are functions of the space coordinates.
r r
V = V ( x , y , z) , p = p ( x , y , z) , ρ = ρ( x , y , z ) , etc.
r
∂V ∂p ∂ρ
or =0 , =0 , =0 , etc.
∂t ∂t ∂t
• Whether the flow is steady or not depends on the choice of the coordinate system with which
the motion is observed.
• An unsteady flow relative to a stationary reference frame can be simplified to a steady flow if
analyzed from a moving reference frame.
Movie Time:
8
Steady and Unsteady Flows (cont’d):
Important: Although the partial derivatives (∂/∂ t) are zero in an steady flow, the total derivatives
(D/dt), which also include the convective derivatives, may be nonzero.
DN ∂N r r
= + ( V ⋅ ∇ )N
dt ∂t
9
Three- , Two- and One-Dimensional Flows:
Three-dimensional flow:
• Generally a fluid flow is a rather complex three-dimensional phenemenon.
• Fluid properties vary in all three mutually perpendicular directions.
• For example in the Cartesian coordinate system
u = u( x , y , z, t ) , v = v ( x , y , z, t ) , w = w ( x , y , z, t ) , p = p(x, y, z, t) , ρ = ρ(x, y, z, t)
• Even in the most simple flows, velocity field generally has all three components nonzero.
• But still, dependence of the fluid properties in one or two directions might be small compared to
the other direction(s).
• In these cases the flow can be simplified as being two- or even one-dimensional.
10
Two-dimensional flow: Properties are functions of time and two space coordinates (say x and y).
u = u( x , y , t ) , v = v ( x , y , t ) , w = w( x , y , t) , p = p(x, y, t) , ρ = ρ(x, y, t)
x u = u( x , y ) , v = 0 , w = 0
p = p(x) , ρ = constant
11
One-dimensional flow: Fluid properties are functions of time and one space coordinate (say x).
u = u( x , t ) , v = v ( x , t ) , w = w( x , t) , p = p(x, t) , ρ = ρ(x, t)
u = u( y ) , v = 0 , w = 0
p = constant , ρ = constant
Uniform flow: For the sake of simplicity, the properties at a cross-section can be approximated by
its average values.
• Then the properties are said to be uniform over that cross-section.
• Streakline: Focus your attention to a fixed point of the flow field and identify all the particles
passing through that point. If you join all these particles with a line you will get a streakline.
• Streaklines are related to the Eulerian description of the flow field.
Movie Time:
• Streamlines: are imaginary lines which are always tangent to the flow direction at a given
instant of time.
• No fluid can pass through a streamline. Therefore solid walls are natural streamlines.
• At the intersection of two streamlines the fluid velocity should be zero.
• For an unsteady flow, streamlines change from time to time, therefore one can talk about
instantaneous streamlines.
Movie Time:
streamlines
V Important:
For a steady flow pathlines, streaklines
and streamlines are the same.
13
Derivation of a streamline:
V = 2x i – 2y j (m/s).
15
3.7 Acceleration of a Fluid Particle
r r
r DV ∂V r r r
a = = + ( V ⋅ ∇) V
dt ∂t
Du ∂u r r ∂u ∂u ∂u ∂u
ax = = + (V ⋅ ∇ )u = + u + v + w
dt ∂t ∂t ∂x ∂y ∂z
Dv ∂v r r ∂v ∂v ∂v ∂v
ay = = + (V ⋅ ∇ )v = + u + v + w
dt ∂t ∂t ∂x ∂y ∂z
Dw ∂w r r ∂w ∂w ∂w ∂w
az = = + (V ⋅ ∇)w = + u + v + w
dt ∂t ∂t ∂x ∂y ∂z
16
Acceleration vector in the cylindrical coordinate system:
v r r r
a = ar ir + a θ i θ + a z i z
• or in more detail
∂Vθ ∂V V ∂Vθ ∂V VV
aθ = + Vr θ + θ + Vz θ + r θ
∂t ∂r r ∂θ ∂z r
∂Vz ∂V V ∂Vz ∂V
az = + Vr z + θ + Vz r
∂t ∂r r ∂θ ∂z
17
Example:
V = x2 i – 2xyt j + 3z k (m/s).
18
Example:
(b) Using the Lagrangian method, determine the equations for the position and acceleration of a
fluid particle as a function of time, which is located at x = 0 at time t = 0.
(c) Show that both expressions for the acceleration give identical results, as this fluid particle
passes through the point at x = L.
∆L
control C
volume
V
B β dAC
A at time t0+∆t
n
system dAC
at time t0 d∀C
20
Mass Flow Rate (cont’d)
mc ∆L r r
m
& = lim
∆t → 0 ∆t
= lim ∫ ρ
∆t → 0 ∆t
cos β dA = ∫ ρ V cos β dA = ∫ ρ V ⋅ n dA
A A A
• Note that the mass in region C flows out through the control surface. Therefore the angle β is
always less than 90o and cosβ is always positive.
• Sign convention for mass flow rate:
• m-dot is positive if the mass goes out of the control volume.
• m-dot is negative if the mass goes into the control volume.
• Note that m-dot can be zero if either V or cos β is zero.
• For incompressible and uniform flow going out of a control volume m-dot = ρ A V
n
V
n n
V V V
n
21
Volumetric Flow Rate (Q) [m3/s]
• For an incompressible fluid the density is constant. Therefore one can define
m
& r r
Q= = ∫ V ⋅ n dA
ρ A
• If the velocity is not uniform over a cross section, then an average velocity, which produces the
same mass flow rate as the original velocity profile, can be defined as
r r u(y) V
m
&
∫ ⋅ n dA
ρ V
A
V= =
∫ ρ dA ∫ ρ dA
A A
r r
Q
∫ V ⋅ n dA
A
• If the fluid is incompressible this simplifies to V= =
A A
22
Example 3.4
Water with a density of 1000 kg/m3 is flowing through a pipe. The velocity of water, which is assumed
to be uniform at the exit of the pipe, is 3 m/s. If the area at the exit is 0.001 m2, determine the mass
flow rate and the volumetric flow rate.
V
n
β=60o
Flow x
Problem 3.13
A viscous and incompressible fluid with a density of ρ is flowing through a wide flat channel with a
heignt of 2h. At a certain section of the channel, the velocity distribution is given by u =UC [1-(y/h)2],
where UC is the centerline velocity. Determine the mass flow rate per unit width, volumetric flow rate
per unit width and the average velocity.
y
u(y) h
x
h
23
Problem 3.14
A viscous and incompressible fluid with a density of ρ is flowing through a pipe with a radius of R. The
velocity distribution in the pipe is given by Vx =UC [1-(r/R)], where UC is the centerline velocity.
Determine the mass flow rate, volumetric flow rate and the average velocity.
r
Vx(r) R
x
24