Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Me305 Chapter 3

Download as pdf or txt
Download as pdf or txt
You are on page 1of 24

ME 305 Fluid Mechanics I

Chapter 3
Introduction to Fluid Dynamics

These presentations are prepared by


Dr. Cüneyt Sert
Mechanical Engineering Department
Middle East Technical University
Ankara, Turkey
cuneyt@me.metu.edu.tr

They can not be used without the permission of the author.

1
Chapter 3 Introduction to Fluid Dynamics

Methods for Describing the Fluid Motion

(a) Material (Lagrangian) Description: identified fluid particles are followed in the course of time
and the variation of their properties are determined.

z A(xA, yA, zA) • Fluid particle P is located at point A at time=0.

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

where xP, yP, zP are the coordinates of the fluid particle.

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.

Material Volume (Closed system): A collection of identified fluid particles.


• Its shape and position may change, but it contains always the same fluid particles.
• It is surrounded by the system boundaries, through which no particle can pass.

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)

ρ P = ρ P ( x P ( t = 0), y P ( t = 0), z P ( t = 0), t = 0 ) = ρ ( x A , y A , z A , t = 0)


4
Methods for the Mathematical Formulation of the Fluid Flow

(a) Differential Formulation:


• When used with the Lagrangian description, basic laws are applied to all fluid particles.
• When used with the Eulerian description, basic laws are applied to an infinitesimal fluid element.
• In either case the final product is a differential equation.
• When this differential equation is solved, detailed information about all fluid particles or about all
points in the flow field is obtained.

(b) Integral Formulation:


• In many engineering problems, one does not require a very detailed information of the flow.
• Then basic laws are applied to a closed system (Lagrangian) or a control volume (Eulerian).
• Overall behavior of the fluid motion is obtained.

Differential Integral

Solve the flow


Flow Thrust
equations inside
in out
the jet engine

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.

• In material description any property N can be expressed as N = N ( x ( t ), y ( t ), z(t) , t )

∂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.

control flow field control


volume volume C
B
A
system
system
time t0 time t0 + ∆t

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

• In an unsteady flow, one or more properties change with time.

• 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

Steady state uniform flow in a converging


diverging nozzle. ∂u/∂t of any fluid particle is
zero, but Du/dt is not.
Steady state operation of a water
heater. ∂T/∂t of any fluid particle
is zero, but DT/dt is not.

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)

Flow visualization of the complex


three-dimensional flow field past
a model airfoil.
Movie Time:

• 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)

• Variation of the properties in the third direction (z) is assumed to be negligible.

Two-dimensional flow between two


y diverging plates.

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.

Uniform one-dimensional flow Uniform one-dimensional flow


in a constant diameter pipe. between two diverging plates.
12
3.6 Visualization of the Fluid Flow
• Pathline: is the path of a fluid particle moving in a flow field.
• To determine a pathline, identify a fluid particle by using dye and take a film of its motion.
• Pathlines are related to the Lagrangian description of the flow field.

• 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:

z • Consider a fluid particle with a position vector r at a certain


dr time.
V instantaneous
• After an infinitesimal time interval, the particle will move to a
streamline
r different location.
r+dr
• During this infinitesimal time interval, the position change of
y
the particle, dr, coincides with the instantaneous streamline.
x

• That is, the velocity vector V=ui+vj+wk is parallel to dr = dx i + dy j + dz k.


dx dy dz
• In other words V x dr = 0. This gives the equation of a streamline as = =
u v w

Streamtube: is formed by a bundle of streamlines passing


through a closed curve in space. streamlines
• For a steady flow, streamtube behaves like a real tube,
because no particles can pass through the boundaries of a
streamtube (they are formed by streamlines).
• For an unsteady flow, streamtube changes its position from
time to time and the fluid particles that form the streamtube
change with time.
14
Example:

• The velocity field for a fluid flow is given as

V = 2x i – 2y j (m/s).

(a) Determine whether the flow field is three-, two- or one-dimensional.


(b) Determine the components of the velocity vector at point P(2,2,0).
(c) Determine the magnitude and direction of the velocity vector at point P.
(d) Determine the equation of the streamline passing through point P.
(e) Sketch the velocity vector and the streamline.

Solution will be given during the class.

15
3.7 Acceleration of a Fluid Particle
r r
r DV ∂V r r r
a = = + ( V ⋅ ∇) V
dt ∂t

• ∂V/∂t is the local acceleration.


• (V.∇)V is the convective acceleration.
• This is a vector equation and can be expressed in the Cartesian coordinate system as follows
r
v r r r DV Du r Dv r Dw r
a = a x i + a y j + a zk = = i + j + k
dt dt dt dt
• or in more detail

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

∂Vr ∂Vr Vθ ∂Vr ∂Vr ∂Vθ2


ar = + Vr + + Vz -
∂t ∂r r ∂θ ∂z r

∂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

• Think about where the additional terms are coming from.


• Similarly acceleration components in the spherical coordinate system () can be written.

17
Example:

• A velocity field is given as

V = x2 i – 2xyt j + 3z k (m/s).

(a) Determine the acceleration field a .


(b) Determine the acceleration at point P(1,2,3) and time t=10 s.

Solution will be given during the class.

18
Example:

• Consider one-dimensional, steady, incompressible flow through a converging channel.

• The velocity field is given by V = Uo (1 + x/L ) i

(a) Determine the acceleration field a by using the Eulerian method.

(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.

Solution will be given during the class.


19
3.8 Flow Rates and Average Velocity
• Flow rate is the quantity of fluid flowing per unit time across any cross section.
• Mass flow rate is commonly used for compressible fluids (mainly gasses).
• Volumetric flow rate is commonly used for incompressible fluids (mainly liquids).

Mass Flow Rate (m-dot) [kg/s]


• Remember the surface integrals in the derivation of the Reynold’s Transport Theorem.

∆L
control C
volume
V
B β dAC
A at time t0+∆t
n
system dAC
at time t0 d∀C

• Mass in the differential volume d∀C is dm = ρ ∆L cosβ dA

• Mass of the fluid in the region C is mC = ∫ ρ ∆L cos β dA


Ac

20
Mass Flow Rate (cont’d)

• The mass flow rate going out of the control volume is

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

• Sign convention for Q is the same as m-dot.


• For uniform flow going out of a control volume Q=AV

Average Velocity (V-bar) [m/s]

• 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

You might also like