Module 3 - FLUID FLOW BASIC EQUATIONS-printable
Module 3 - FLUID FLOW BASIC EQUATIONS-printable
This lecture will introduce the student to some of the important aspects of fluid flow (momentum transfer). Later,
other subjects will be introduced to give the fledgling engineer the competence required to meaningfully deal with
this overall area of momentum transport and fluid flow.
The subjects in this lecture will include:
Fluid statics
Fluid flow phenomena,
Categories of fluid flow behavior,
Equations of change relating the momentum transport,
Macroscopic approach to fluid flow.
FLUID STATICS:
The behavior of fluids at rest
This can be involved in process operations in a number of ways. One such case is in the measurement of
pressure differentials in a system.
¦
Fluid Static 1
Variation of pressure with depth
The pressure in a fluid at rest does not change in the horizontal direction. This can be shown easily
by considering a thin horizontal layer of fluid and doing a force balance in any horizontal direction.
However, in the vertical direction in a gravity field, pressure in a fluid increases with depth
because more fluid rests on deeper layers, and the effect of this “extra weight” on a deeper layer
is balanced by an increase in pressure as shown in the figure below:
o The pressure of a fluid at rest increases with depth (as a result of added weight)
o Free-body diagram of a rectangular fluid element in equilibrium.
Let:
P1 = Patm , be the pressure at Pt. 1 (free surface of a liquid open to atmosphere)
Therefore:
𝐠
𝐏𝟐 = 𝐏𝐚𝐭𝐦 + 𝛒𝐡
𝐠𝐜
𝐠
𝐨𝐫 𝐏𝐠𝐚𝐠𝐞 = 𝐏𝟐 − 𝐏𝐚𝐭𝐦 = 𝛒𝐡
𝐠𝐜
Pressure in a liquid at rest increases linearly with distance from the free surface.
For gases at rest and for small to moderate distances, the variation of pressure with height is negligible because
of their low density.
𝐠
𝐏𝐛𝐨𝐭𝐭𝐨𝐦 = 𝐏𝐭𝐨𝐩 + 𝛒𝐡 𝐠 -
𝐜
𝐠
𝐨𝐫 𝐏𝐠𝐚𝐠𝐞 = 𝐏𝐛𝐨𝐭𝐭𝐨𝐦 − 𝐏𝐭𝐨𝐩 = 𝛒𝐡 ≅𝟎
𝐠𝐜
The pressure in tank or in a room containing a gas can be considered to be uniform since the weight of the
gas is too small to make significant difference
𝐠
𝐏𝐀 = 𝐏𝐁 = 𝐏𝐂 = 𝐏𝐃 = 𝐏𝐄 = 𝐏𝐅 = 𝐏𝐆 = 𝐏𝐚𝐭𝐦 + 𝛒𝐡
𝐠𝐜
𝐏𝐇 ≠ 𝐏𝐈
Pascal Law:
A consequence of the pressure in a fluid remaining constant in the horizontal direction is that that
the pressure applied to a confined fluid increases the pressure throughout by the same amount.
o Lifting of a large weight by a small force by the application of Pascal’s law (Pascal’s
Machine).
Since ∶ P1 = P2
F1 F2 F2 A 2
Also, = → =
A1 A2 F1 A1
A
∗ 2 is known as "ideal mechanical advantage of the hydraulic lift
A1
A2
∗ Using a hydraulic car jack with a piston area ratio of = 10
A1
∗ A person can lift a 1000 − kg car by applying a force of just 100 kgf (908N)
Barometic Equations:
A. Directly Vertical
An important equation relating to fluid statics is the barometric equation:
𝐝𝐏 𝐠
=𝛒 = (𝐃𝐢𝐫𝐞𝐜𝐭𝐥𝐲 𝐕𝐞𝐫𝐭𝐢𝐜𝐚𝐥)
𝐝𝐙 𝐠𝐜
𝐠
or : ∆𝐏 = 𝐏𝟐 − 𝐏𝟏 = 𝛒 ∆𝐳
𝐠𝐜
Where:
P = pressure; Z = vertical distance; 𝛒 = fluid density;
𝐟𝐭 𝐦
g = acceleration of gravity = 𝟑𝟐. 𝟐 𝟐 𝐨𝐫 𝟗. 𝟖𝟏 𝟐
𝐬𝐞𝐜 𝐬𝐞𝐜
𝐥𝐛 𝐦𝐚𝐬𝐬 − 𝐟𝐞𝐞𝐭 𝐤𝐠 𝐦𝐚𝐬𝐬 − 𝐦
𝐠 𝐜 = 𝐠𝐫𝐚𝐯𝐢𝐭𝐚𝐭𝐢𝐨𝐧𝐚𝐥 𝐜𝐨𝐧𝐬𝐭𝐚𝐧𝐭 = 𝟑𝟐. 𝟐 𝟐
𝐨𝐫 𝟏
𝐬𝐞𝐜 − 𝐥𝐛 𝐟𝐨𝐫𝐜𝐞 𝐬𝐞𝐜 𝟐 − 𝐧𝐞𝐰𝐭𝐨𝐧
Hydrostatic Equilibrium
For engineering calculations and for incompressible fluid with density constant and between the two
definite heights 𝑍𝑎 𝑎𝑛𝑑 𝑍𝑏 ,
𝐏𝐛 𝐏𝐚 𝐠
− = (𝐙𝐚 − 𝐙𝐛 )
𝛒 𝛒 𝐠𝐜
𝐏𝐛 − 𝐏𝐚 𝐠
= (𝐙𝐚 − 𝐙𝐛 )
𝛒 𝐠𝐜
∆𝐏 𝐠
= ∆𝐙
𝛒 𝐠𝐜
o “Consider the barometric equation with respect to the world itself. If, for example, we find
ourselves on the surface of the Pacific Ocean, then we know that the pressure is atmospheric. On
the other hand, if we dive to the bottom of the ocean until you reach the ocean floor, we would find
that the pressure is many times atmospheric with the difference due to the effect predicted by the
barometric equation for ocean water.
o Likewise, if we would go to California, we would find that the pressure in Death Valley (282 feet
below sea level) is higher than that at the top of Mount Apo (elevation 14,494 feet). Once again the
difference in pressure would be governed by the barometric equation.”
B. Directly Not Vertical direction
𝐝𝐏 𝐠
=𝛒 𝐂𝐨𝐬 𝛉 (𝐃𝐢𝐫𝐞𝐜𝐭𝐥𝐲 𝐧𝐨𝐭 𝐕𝐞𝐫𝐭𝐢𝐜𝐚𝐥)
𝐝𝐥 𝐠𝐜
Where:
𝜃 = angle between the vertical Z and the dimension 𝒍
*Note:
In using the barometric equation, we also must consider the usage of the terms absolute, gauge, and
atmospheric pressure.
𝐏𝐀𝐛𝐬𝐨𝐥𝐮𝐭𝐞 = 𝐏𝐠 + 𝐏𝐀𝐭𝐦
Where:
𝑃𝐴𝑡𝑚 = ambient pressure at that point where we make the reading.
𝑃𝑔 = gauge pressure read by some measuring device for a vessel or a container.
C. Barometric equation for ideal gas:
𝐝𝐏 𝐠
=𝛒 (𝐃𝐢𝐫𝐞𝐜𝐭𝐥𝐲 𝐕𝐞𝐫𝐭𝐢𝐜𝐚𝐥)
𝐝𝐙 𝐠𝐜
𝐏𝐌
𝐁𝐮𝐭 ∶ 𝛒 =
𝐑𝐓
𝐠 𝐏𝐌 𝐠
𝐝𝐏 = 𝛒 𝐝𝐙 = ( ) 𝐝𝐙
𝐠𝐜 𝐑𝐓 𝐠 𝐜
𝐝𝐏 𝐠𝐌 𝐝𝐏 𝐠𝐌
= 𝐝𝐙 ; 𝐀𝐥𝐬𝐨, − 𝐝𝐙 = 𝟎
𝐏 𝐠 𝐜 𝐑𝐓 𝐏 𝐠 𝐜 𝐑𝐓
𝐏𝐛 𝐠𝐌
𝐥𝐧 + (𝐙 − 𝐙𝐚 ) = 𝟎 (𝐖𝐡𝐞𝐫𝐞 𝐙𝐚 > 𝐙𝐛 )
𝐏𝐚 𝐠 𝐜 𝐑𝐓 𝐛
𝐏𝐛 𝐠𝐌 𝐏𝐛 𝐠𝐌(𝐙𝐛−𝐙𝐚 )
[− ]
𝐥𝐧 =− (𝐙 − 𝐙𝐚 ) 𝐨𝐫 =𝐞 𝐠𝐜 𝐑𝐓
𝐏𝐚 𝐠 𝐜 𝐑𝐓 𝐛 𝐏𝐚
Fluid Static 2
Another aspect of fluid statics that will be utilized later in this text is the Principle of Archimedes relating
to buoyant force.
Buoyant force for an object floating in a liquid is:
𝐠 𝐠
𝐅𝐁 = [𝛒𝐋𝐢𝐪𝐮𝐢𝐝 𝐕𝐋𝐢𝐪𝐮𝐢𝐝 + 𝛒𝐀𝐢𝐫 𝐕𝐀𝐢𝐫 ]
𝐠𝐜 𝐠𝐜
Where:
𝑉𝐿𝑖𝑞𝑢𝑖𝑑 = 𝑡ℎ𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 𝑜𝑓 𝑙𝑖𝑞𝑢𝑖𝑑 𝑑𝑖𝑠𝑝𝑙𝑎𝑐𝑒𝑑 𝑏𝑦 𝑡ℎ𝑒 𝑜𝑏𝑗𝑒𝑐𝑡
𝑉𝐴𝑖𝑟 = 𝑡ℎ𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 𝑜𝑓 𝐴𝑖𝑟 𝑑𝑖𝑠𝑝𝑙𝑎𝑐𝑒𝑑
Hydrostatic Equilibrium in a centrifugal Field:
F = ma = mass x acceleration
𝐹𝑜𝑟 𝑟𝑜𝑡𝑎𝑡𝑖𝑛𝑔 𝑙𝑖𝑞𝑢𝑖𝑑 𝑜𝑟 𝑖𝑛 𝑐𝑖𝑟𝑐𝑢𝑙𝑎𝑟 𝑚𝑜𝑡𝑖𝑜𝑛, 𝑎 = 𝑟𝜔2
𝐝𝐅 = 𝐫𝛚𝟐 𝐝𝐦
Where: dF = centrifugal force
dm = mass of liquid
ω = angular velocity, rad/s
r = radial distance
𝐝𝐅 = 𝐫𝛚𝟐 𝐝𝐦 = 𝐫𝛚𝟐 [𝟐𝛑𝛒𝐫𝐛𝐝𝐫]
dF
= rω2 [ρdr] = dP
2πrb
P2 r2
∫ dP = ω2 ρ ∫ rdr
P1 r1
𝛚𝟐 𝛒(𝐫𝟐 𝟐 − 𝐫𝟏 𝟐 )
𝐏𝟐 − 𝐏𝟏 = (𝐰𝐡𝐞𝐫𝐞 𝐫𝟐 − 𝐫𝟏 ≈ 𝐬𝐦𝐚𝐥𝐥 𝐝𝐢𝐟𝐟𝐞𝐫𝐞𝐧𝐜𝐞)
𝟐
Manometry:
Suppose a manometer is used to measure a pressure differential in a pipe with a flowing fluid ‘‘x” at room
temperature as shown in figure below The manometer reads a differential height of 0.3323 m. The liquid in the pipe
has a density of 1260 kg mass/ m3 . Mercury (density of 13,600 kg/ m3) is the manometer fluid. What is the pressure
measured?
𝐦 𝑘𝑔 𝑚𝑎𝑠𝑠−𝑚
For metric, 𝑔 = 𝟗. 𝟖𝟏 𝒂𝒏𝒅 𝒈𝒄 = 1
𝐬𝐞𝐜 𝟐 𝑠𝑒𝑐 2 −𝑛𝑒𝑤𝑡𝑜𝑛
𝐦
𝟗. 𝟖𝟏 𝟐 kgm
Therefore, P5 − P1 = 𝐬𝐞𝐜 (0.3323 m)(1260 − 13,600) 3
𝑘𝑔 𝑚𝑎𝑠𝑠 − 𝑚 m
1
𝑠𝑒𝑐 2 − 𝑛𝑒𝑤𝑡𝑜𝑛
𝑁 𝑁
P5 − P1 = −40,226.709 or P1 − P5 = 40,226.709 ≈ 40,266.709 Pa
m2 m2
𝑃1 𝑃2
𝐠
∆𝐏 = 𝐏𝟐 − 𝐏𝟏 = 𝐑 𝐦 (𝛒𝐦 − 𝛒𝐟 )
𝐠𝐜
Where:
P2 = Pressure on the depressed side
P1 = Pressure on the Elevated side
∆h = R m = depressed distance
ρm = higher density or manometer fluid
ρf = lighter density fluid
g
∆P = R m ρm
gc
For Well-type manometer,
𝑨𝒄 𝒈
∆𝑷 = 𝑷𝟐 − 𝑷𝟏 = 𝝆𝒇 𝑹𝒎 (𝟏 + )
𝑨𝒘 𝒈𝒄
Where:
𝑃2 = 𝑃𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑜𝑛 𝑡ℎ𝑒 𝑑𝑒𝑝𝑟𝑒𝑠𝑠𝑒𝑑 𝑠𝑖𝑑𝑒
𝑃1 = 𝑃𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑜𝑛 𝑡ℎ𝑒 𝐸𝑙𝑒𝑣𝑎𝑡𝑒𝑑 𝑠𝑖𝑑𝑒
𝑦 = 𝑅𝑚 = 𝑑𝑒𝑝𝑟𝑒𝑠𝑠𝑒𝑑 𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑒
𝜌𝑓 = 𝑙𝑖𝑔ℎ𝑡𝑒𝑟 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑓𝑙𝑢𝑖𝑑
𝐴𝑐 = 𝑐𝑟𝑜𝑠𝑠 𝑠𝑒𝑐𝑡𝑖𝑜𝑛 𝑎𝑟𝑒𝑎 𝑜𝑓 𝑐𝑜𝑙𝑢𝑚𝑛
𝐴𝑤 = 𝑐𝑟𝑜𝑠𝑠 𝑠𝑒𝑐𝑡𝑖𝑜𝑛 𝑎𝑟𝑒𝑎 𝑜𝑓 𝑤𝑒𝑙𝑙
If the well large in comparison to the column,
𝒈
∆𝐏 = 𝑷𝟐 − 𝑷𝟏 = 𝑹𝒎 𝝆𝒇
𝒈𝒄
𝒈
𝑷𝒂 − 𝑷𝒃 = 𝑹𝟏 (𝝆𝑨 − 𝝆𝑩 )𝒔𝒊𝒏 𝜶
𝒈𝒄
𝑔
∆𝑷 = (𝑷𝟐 − 𝑷𝟏 ) = 𝒉𝒎 (𝝆𝒎𝟏 − 𝝆𝒎𝟐 )
𝑔𝑐
The water in a tank is pressured by air, and the pressure is measured by a multi-fluid manometer as shown
above. The tank is located on a mountain at an altitude of 1400 meters where the atmospheric pressure
is 85.6 kPa. Determine the air pressure in the tank if ℎ1 = 0.10 𝑚 , ℎ2 = 0.20 𝑚 𝑎𝑛𝑑 ℎ3 = 0.35 𝑚. Take
𝑘𝑔 𝑘𝑔 𝑘𝑔
the densities of water, oil and mercury to be 1000 3 , 850 3 , 𝑎𝑛𝑑 13,600 3 𝑟𝑒𝑠𝑝𝑒𝑐𝑡𝑖𝑣𝑒𝑙𝑦.
𝑚 𝑚 𝑚
Solution:
The air pressure in the pressurized water tank is measured by a multi-fluid manometer
g g g
P1 + ρwater h1 + ρoil h2 − ρmercury h3 = Patm
gc gc gc
g
P1 = Patm − (ρwater h1 + ρoil h2 − ρmercury h3 )
gc
m
9.81 2 𝑘𝑔 𝑘𝑔 𝑘𝑔
𝑃1 = 85.6 𝑘𝑃𝑎 − ( sec ) [(1000 3 ) (0.1𝑚) + (850 3 ) (0.2𝑚) − (13,600 3 ) (0.35𝑚)]
𝑘𝑔 𝑚𝑎𝑠𝑠 − 𝑚 𝑚 𝑚 𝑚
1
𝑠𝑒𝑐 2 − 𝑛𝑒𝑤𝑡𝑜𝑛
𝑵 𝟏 𝒌𝑷𝒂
𝑷𝟏 = 𝟖𝟓. 𝟔 𝒌𝑷𝒂 − [𝟗𝟖𝟏 + 𝟏𝟔𝟔𝟕. 𝟕 − 𝟒𝟔, 𝟔𝟗𝟓. 𝟔] ( ) = 𝟏𝟐𝟗. 𝟔𝟒𝟕 𝒌𝑷𝒂
𝒎 𝟏𝟎𝟎𝟎 𝑵
𝟐
𝒎𝟐
FLUID DYNAMICS -A PHENOMONOLOGICAL APPROACH
The fluid is put into motion,
Fluid behavior is determined by its physical nature, the flow geometry, and its velocity.
Osborne Reynolds’ Experiment on Newtonian fluids (those that obey Newton’s Law of Viscosity).
Conclusion:
These behavior patterns led to conclusion that such flows were streamline (i.e., the dye stream showed a
straight line behavior with a given velocity at a circumference). Furthermore, since the velocity moved from
a maximum at the tube center to a minimum at the wall, the fluid itself moved in shells or lamina. Because
of these patterns of behavior, it was termed such flows as streamline and laminar
𝐈𝐧𝐞𝐫𝐭𝐢𝐚 𝐟𝐨𝐫𝐜𝐞𝐬 ̅𝟐
𝝆𝝊 𝑫𝝊̅ 𝝆
= = = 𝒅𝒊𝒎𝒆𝒏𝒔𝒊𝒐𝒏𝒍𝒆𝒔𝒔 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕 = 𝑹𝒆𝒚𝒏𝒐𝒍𝒅𝒔 𝒏𝒖𝒎𝒃𝒆𝒓 = 𝑹𝒆
𝐕𝐢𝐬𝐜𝐨𝐮𝐬 𝐅𝐨𝐫𝐜𝐞𝐬 𝝁𝝊̅ 𝝁
𝑫
Flow characteristics as a function of the shape of relative velocity profiles/distribution in a tube or pipe:
*Note: The Ultimate beyond Laminar and Turbulent is “Ideal” or Plug flow where all velocities
across the tube cross section are the same just like in the flow reactors or flow-through packed bed.
CLASSIFICATION OF FLUID BEHAVIOR
Many fluids obey Newton’s Law of Viscosity. In these fluids, which are called Newtonian, the viscosity is
a property of the system.
As such, it depends on the substance or substances in the system, temperature, and pressure but not on the
velocity gradient, which is the rate of shear or on the time parameter.
In nature there are many other fluid systems other than the Newtonian fluids. Therefore, the general
classification of fluid system according to its behavior are:
1. Newtonian
Newtonian fluids are named after Sir Issac Newton (1642 - 1726) who described the flow
behavior of fluids with a simple linear relation between shear stress [mPa] and shear rate [1/s].
Some examples of Newtonian fluids include water, organic solvents, and honey. For those
fluids viscosity is only dependent on temperature.
This means that the viscosity of Newtonian fluids will remain a constant no matter how fast
they are forced to flow through a pipe or channel (i.e. viscosity is independent of the rate of
shear).
An exception to the rule is Bingham plastics, which are fluids that require a minimum stress to
be applied before they flow. These are strictly non-Newtonian, but once the flow starts they
behave essentially as Newtonian fluids (i.e. shear stress is linear with shear rate). A great
example of this kind of behavior is mayonnaise.
2. “Simple” non-Newtonian (the viscosity is a function of shear rate )
In reality most fluids are non-Newtonian, which means that their viscosity is dependent on
shear rate (Shear Thinning or Thickening) or the deformation history (Thixotropic fluids). In
contrast to Newtonian fluids, non-Newtonian fluids display either a non-linear relation between
shear stress and shear rate, have a yield stress, or viscosity that is dependent on time or
deformation history (or a combination of all the above).
3. “Complex” non-Newtonian (the behavior is a function of both rate of shear and the time
parameter)
4. Fluids influenced by external force fields
5. Fluids that are non-continuous
6. Relativistic fluids
“Simple” non-Newtonian (the viscosity is a function of shear rate)
Those fluids in which the rate of shear influences the flow behavior
𝒅𝑽𝑿
Comparison shear stress (𝝉) vs. shear rate ( )
𝒅𝒚
*Newtonian (A), shear-thinning (B), and shear-thickening (C) fluids.
Using Newton’s Law of Viscosity so that we have an apparent viscosity (𝝁𝒂𝒑𝒑 ) :
𝒅𝑽𝑿
𝝉𝒚𝒙 = −𝝁𝒂𝒑𝒑 (𝑵𝒆𝒘𝒕𝒐𝒏’𝒔 𝑳𝒂𝒘 𝒐𝒇 𝑽𝒊𝒔𝒄𝒐𝒔𝒊𝒕𝒚)
𝒅𝒚
𝝉𝒚𝒙
𝝁𝒂𝒑𝒑 =
𝒅𝑽
− 𝑿
𝒅𝒚
If we plot the behavior of apparent viscosity, 𝝁𝒂𝒑𝒑 on a log-log plot, we obtain Figure below:
𝑑𝑉𝑋
Log of apparent viscosity (𝜇𝑎𝑝𝑝 ) vs. log of shear rate ( )
𝑑𝑦
*Newtonian (A), shear-thinning (B), and shear-thickening (C) fluids.
Figure Interpretation:
The apparent viscosity, 𝝁𝒂𝒑𝒑 for a Newtonian is a constant.
The other fluids show a decreasing apparent viscosity with increasing shear rate (shear-thinning fluid) or an
increasing apparent viscosity with an increasing shear rate (shear-thickening fluid).
Shear-thinning fluids are also called pseudoplastic, whereas shear-thickening fluids are termed dilatant.
The overall science that considers flow and deformation of fluids (as well as solids) is termed “rheology”.
The velocity profile of a pseudoplastic fluid flowing in a circular tube is not a parabola (not even in
laminar flow) but rather a blunted profile (same as that of a Newtonian turbulent flow).
Thixotropy is a time-dependent shear thinning property. Certain gels or fluids that are thick (viscous) under
static conditions will flow (become thin, less viscous) over time when shaken, agitated, sheared or
otherwise stressed (time dependent viscosity). They then take a fixed time to return to a more viscous state.
In more technical language: some non-Newtonian pseudoplastic fluids show a time-dependent change
in viscosity; the longer the fluid undergoes shear stress, the lower its viscosity. A thixotropic fluid is a fluid
which takes a finite time to attain equilibrium viscosity when introduced to a step change in shear rate.
Some thixotropic fluids return to a gel state almost instantly, such as ketchup, and are
called pseudoplastic fluids. Others such as yogurt take much longer and can become nearly solid.
Many gels and colloids are thixotropic materials, exhibiting a stable form at rest but becoming fluid when
agitated.
Some fluids are anti-thixotropic: constant shear stress for a time causes an increase in viscosity or even
solidification. Constant shear stress can be applied by shaking or mixing. Fluids which exhibit this property
are usually called rheopectic. They are much less common.
“Complex” non-Newtonian:
Fluids in which the time parameter becomes a factor.
𝒅𝜸 𝒅𝜸
𝝉 = 𝝓 ( , 𝒕) 𝑰𝒇 𝒇𝒐𝒓 𝑵𝒆𝒘𝒕𝒐𝒏𝒊𝒂𝒏 𝒇𝒍𝒖𝒊𝒅 ∶ 𝝉 = 𝝁
𝒅𝒕 𝒅𝒕
MULTIDIMENSIONAL FLUID FLOW - THE EQUATIONS OF CHANGE
Navier-Stoke Equations:
In physics, the Navier–Stokes equations describe the motion of viscous fluid substances.
These balance equations arise from applying Newton's second law to fluid motion, together with the
assumption that the stress in the fluid is the sum of a diffusing viscous term (proportional to the gradient of
velocity) and a pressure term—hence describing viscous flow.
Navier–Stokes equations are useful because they describe the physics of many phenomena
of scientific and engineering interest. They may be used to model the weather, ocean currents, water flow
in a pipe and air flow around a wing.
The Navier–Stokes equations in their full and simplified forms help with the design of aircraft and cars, the
study of blood flow, the design of power stations, the analysis of pollution, and many other things.
The solution of the Navier–Stokes equations is a flow velocity. It is a field, since it is defined at every point
in a region of space and an interval of time. Once the velocity field is calculated, other quantities of
interest, such as pressure or temperature, may be found.
The Navier-Stokes equations, developed by Claude-Louis Navier and George Gabriel Stokes in 1822, are
equations which can be used to determine the velocity vector field that applies to a fluid, given some
initial conditions. They arise from the application of Newton’s second law in combination with a fluid
stress (due to viscosity) and a pressure term.
General Equation: ∑𝑭 = 𝒎𝒂
For fluid infinitesimally small
Consider the conservation of mass, using “stationary Volume element for mass balance” or on a per unit
volume analysis,
Mass rate in - Mass rate out = Mass accumulation rate
Based on the above figure , by letting ∆𝑥∆𝑦∆𝑧 approach zero, the differential equation will be:
𝝏𝝆 𝝏 𝝏 𝝏
= − ( (𝝆𝝊𝒙 ) + (𝝆𝝊𝒚 ) + (𝝆𝝊𝒛 )) (𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑓𝑜𝑟 𝑟𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑐𝑜𝑜𝑟𝑑𝑖𝑛𝑎𝑡𝑒𝑠)
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏𝝆 𝝏 𝝏 𝝏
𝒐𝒓 + (𝝆𝝊𝒙 ) + (𝝆𝝊𝒚 ) + (𝝆𝝊𝒛 ) = 𝟎 (𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑓𝑜𝑟 𝑟𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑐𝑜𝑜𝑟𝑑𝑖𝑛𝑎𝑡𝑒𝑠)
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛
Using differential and separated into density and velocity derivatives:
𝜕𝜌 𝜕𝜐𝑥 𝜕𝜐𝑦 𝜕𝜐𝑧 𝜕𝜌 𝜕𝜌 𝜕𝜌
= −𝜌 ( + + ) = − [𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 ]
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝝏𝝆 𝝏𝝆 𝝏𝝆 𝝏𝝆
= −𝝆(𝛁 ∙ 𝝊) = − [𝝊𝒙 + 𝝊𝒚 + 𝝊𝒛 ]
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛
For a fluid of constant density (i.e., incompressible, usually being assumed in the fluid transport):
𝜕𝜌 𝜕𝜌 𝜕𝜌 𝜕𝜌
𝑆𝑖𝑛𝑐𝑒, = −𝜌(∇ ∙ 𝜐) = − [𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 ] = 0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
(Momentum rate in) - (Momentum rate out) + (Sum of all forces acting on the system)
= (Momentum accumulation rate)
Note:
𝜕𝜐𝑥
The momentum rate expressions would include both a convective (for example, 𝜌𝜐𝑥 ) and
𝜕𝑥
a molecular transport term (involving shear stress).
The force term would include both pressure and gravity forces.
Using the approach for the equation of continuity and letting ∆𝑥∆𝑦∆𝑧 approach zero,
𝜕𝜌 𝜕 𝜕 𝜕
= − ( (𝜌𝜐𝑥 ) + (𝜌𝜐𝑦 ) + (𝜌𝜐𝑧 )) (𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑓𝑜𝑟 𝑟𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑐𝑜𝑜𝑟𝑑𝑖𝑛𝑎𝑡𝑒𝑠)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝜌 𝜕 𝜕 𝜕
𝑜𝑟 + (𝜌𝜐𝑥 ) + (𝜌𝜐𝑦 ) + (𝜌𝜐𝑧 ) = 0 (𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑓𝑜𝑟 𝑟𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑐𝑜𝑜𝑟𝑑𝑖𝑛𝑎𝑡𝑒𝑠)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
a. For Simple non-Newtonian Fluids (∑ 𝑭 = 𝒎 ∙ 𝒂)
For x-component:
For y-component:
𝜕𝜐𝑦 𝜕𝜐𝑦 𝜕𝜐𝑦 𝜕𝜐𝑦 𝜕𝑃 𝜕𝜏𝑥𝑦 𝜕𝜏𝑦𝑦 𝜕𝜏𝑧𝑦
𝜌( + 𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 ) = −𝜌𝑔𝑦 − −( + + )
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧
For z-component:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 𝜕𝜏𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜏𝑧𝑧
𝜌( + 𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 ) = −𝜌𝑔𝑧 − −( + + )
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
For x-component: ∑ 𝑭𝒙 = 𝒎 ∙ 𝒂𝒙
𝜕𝑃 𝜕 2 𝜐𝑥 𝜕 2 𝜐𝑥 𝜕 2 𝜐𝑥 𝜕𝜐𝑥 𝜕𝜐𝑥 𝜕𝜐𝑥 𝜕𝜐𝑥
+𝜌𝑔𝑥 − +𝜇( 2 + 2
+ 2
) = 𝜌( + 𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 )
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝜐𝑥 𝜕𝜐𝑥 𝜕𝜐𝑥 𝜕𝜐𝑥
∗∗ 𝑤ℎ𝑒𝑟𝑒 ∶ = 𝑙𝑜𝑐𝑎𝑙 𝑎𝑐𝑐𝑒𝑙𝑒𝑟𝑎𝑡𝑖𝑜𝑛 ; 𝜐𝑥 , 𝜐𝑦 𝑎𝑛𝑑 𝜐𝑧 = 𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑣𝑒 𝑎𝑐𝑐𝑒𝑙𝑒𝑟𝑎𝑡𝑖𝑜𝑛
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
For y-component: ∑ 𝑭𝒚 = 𝒎 ∙ 𝒂𝒚
𝜕𝑃 𝜕 2 𝜐𝑦 𝜕 2 𝜐𝑦 𝜕 2 𝜐𝑦 𝜕𝜐𝑦 𝜕𝜐𝑦 𝜕𝜐𝑦 𝜕𝜐𝑦
+𝜌𝑔𝑦 − +𝜇( 2 + 2
+ 2
) = 𝜌( + 𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 )
𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
For z-component: ∑ 𝑭𝒛 = 𝒎 ∙ 𝒂𝒛
𝜕𝑃 𝜕 2 𝜐𝑧 𝜕 2 𝜐𝑧 𝜕 2 𝜐𝑧 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝜐𝑧
+𝜌𝑔𝑧 − +𝜇( 2 + 2
+ 2
) = 𝜌( + 𝜐𝑥 + 𝜐𝑦 + 𝜐𝑧 )
𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
Mass Balance:
𝜕𝜌 1 𝜕 1 𝜕 𝜕
+ (𝜌𝑟𝜐𝑟 ) + (𝜌𝜐𝜃 ) + (𝜌𝜐𝑧 ) = 0 (𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑓𝑜𝑟 𝑐𝑦𝑙𝑖𝑛𝑑𝑟𝑖𝑐𝑎𝑙)
𝜕𝑡 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
For r-component:
𝜕𝜐𝑟 𝜕𝜐𝑟 𝜐𝜃 𝜕𝜐𝑟 𝜐𝜃 2 𝜕𝜐𝑟 𝜕𝑃 𝜕 1 𝜕(𝑟𝜐𝑟 ) 1 𝜕 2 𝜐𝑟 2 𝜕𝜐𝜃 𝜕 2 𝜐𝑟
𝜌( + 𝜐𝑟 + − + 𝜐𝑧 )= − +𝜇[ ( )+ 2 − + ] + 𝜌𝑔𝑟
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
For 𝜃-component:
𝜕𝜐𝜃 𝜕𝜐𝜃 𝜐𝜃 𝜕𝜐𝜃 𝜐𝑟 𝜐𝜃 𝜕𝜐𝜃 1 𝜕𝑃 𝜕 1 𝜕(𝑟𝜐𝜃 ) 1 𝜕 2 𝜐𝜃 2 𝜕𝜐𝑟 𝜕 2 𝜐𝜃
𝜌( + 𝜐𝑟 + + + 𝜐𝑧 )= − +𝜇[ ( )+ 2 + + ] + 𝜌𝑔𝜃
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝑟 𝜕𝜃 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
For z-component:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 1𝜕 𝜕𝜐𝑧 1 𝜕 2 𝜐𝑧 𝜕 2 𝜐𝑧
𝜌( + 𝜐𝑟 + + 𝜐𝑧 )= − +𝜇[ (𝑟 )+ 2 + ] + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2
Mass Balance:
For 𝜽-component:
𝜕𝜐𝜃 𝜕𝜐𝜃 𝜐𝜃 𝜕𝜐𝜃 𝜐𝜙 𝜕𝜐𝜃 𝜐𝑟 𝜐𝜃 𝜐𝜙 2 cot 𝜃
𝜌( + 𝜐𝑟 + − + − )
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟𝑠𝑖𝑛 𝜃 𝜕𝜙 𝑟 𝑟
2
1 𝜕𝑃 1 𝜕(𝑟 𝜏𝑟𝜃 ) 1 𝜕(𝜏𝜃𝜃 𝑆𝑖𝑛 𝜃) 1 𝜕𝜏𝑟𝜙 𝜏𝑟𝜃 𝐶𝑜𝑡 𝜃
= − −( 2 + + + − 𝜏𝜙𝜙 ) + 𝜌𝑔𝜃
𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟𝑆𝑖𝑛𝜃 𝜕𝜃 𝑟 𝑆𝑖𝑛𝜃 𝜕𝜙 𝑟 𝑟
For 𝝓-component:
𝜕𝜐𝜙 𝜕𝜐𝜙 𝜐𝜃 𝜕𝜐𝜙 𝜐𝜙 𝜕𝜐𝜙 𝜐𝜙 𝜐𝑟 𝜐𝜃 𝜐𝜙
𝜌( + 𝜐𝑟 + + + + cot 𝜃)
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟𝑠𝑖𝑛 𝜃 𝜕𝜙 𝑟 𝑟
1 𝜕𝑃 1 𝜕(𝑟 2 𝜏𝑟𝜙 ) 1 𝜕(𝜏𝜃𝜙 ) 1 𝜕𝜏𝜙𝜙 𝜏𝑟𝜙 2 cot 𝜃
= − −( 2 + + + + 𝜏𝜃𝜙 ) + 𝜌𝑔𝜙
𝑟𝑆𝑖𝑛𝜃 𝜕𝜙 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝑆𝑖𝑛𝜃 𝜕𝜙 𝑟 𝑟
For r-component:
𝜕𝜐𝑟 𝜕𝜐𝑟 𝜐𝜃 𝜕𝜐𝑟 𝜐𝜙 𝜕𝜐𝑟 𝜐𝜃 2 + 𝜐𝜙 2
𝜌( + 𝜐𝑟 + + − )
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟𝑠𝑖𝑛 𝜃 𝜕𝜙 𝑟
𝜕𝑃 1 𝜕2 1 𝜕 𝜕𝜐𝑟 1 𝜕 2 𝜐𝑟
= − + 𝜇 [ 2 2 (𝑟 2 𝜐𝑟 ) + 2 (𝑠𝑖𝑛 𝜃 )+ 2 2
( 2 )] + 𝜌𝑔𝑟
𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝑆𝑖𝑛 𝜃 𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑖𝑛 𝜃 𝜕𝜙
For 𝜽-component:
𝜕𝜐𝜃 𝜕𝜐𝜃 𝜐𝜃 𝜕𝜐𝜃 𝜐𝜙 𝜕𝜐𝜃 𝜐𝑟 𝜐𝜃 𝜐𝜙 2
𝜌( + 𝜐𝑟 + + + − cot 𝜃)
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟𝑠𝑖𝑛 𝜃 𝜕𝜙 𝑟 𝑟
1 𝜕𝑃 1 𝜕 𝜕𝜐𝜃 1 𝜕 1 𝜕 𝜐𝜃 sin 𝜃 1 𝜕 2 𝜐𝜃 2 𝜕𝜐𝑟 2 cos 𝜃 𝜕𝜐𝜙
= − + 𝜇 [ 2 (𝑟 2 )+ 2 ( )+ 2 ( )+ 2 − ] + 𝜌𝑔𝜃
𝑟 𝜕𝜃 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝑠𝑖𝑛 𝜃 𝜕𝜃 𝑟 𝑠𝑖𝑛2 𝜃 𝜕𝜙 2 𝑟 𝜕𝜃 𝑟 2 𝑠𝑖𝑛2 𝜃 𝜕𝜙
For 𝝓-component:
𝜕𝜐𝜙 𝜕𝜐𝜙 𝜐𝜃 𝜕𝜐𝜙 𝜐𝜙 𝜕𝜐𝜙 𝜐𝜙 𝜐𝑟 𝜐𝜃 𝜐𝜙
𝜌( + 𝜐𝑟 + + + − cot 𝜃)
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟𝑠𝑖𝑛 𝜃 𝜕𝜙 𝑟 𝑟
1 𝜕𝑃 1 𝜕 𝜕𝜐𝜙 1 𝜕 1 𝜕 𝜐𝜙 sin 𝜃 1 𝜕 2 𝜐𝜙 2 𝜕𝜐𝑟 2 cos 𝜃 𝜕𝜐𝜃
= − + 𝜇 [ 2 (𝑟 2 )+ 2 ( )+ 2 ( )+ 2 + ]
𝑟 sin 𝜃 𝜕𝜙 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝑠𝑖𝑛 𝜃 𝜕𝜃 𝑟 𝑠𝑖𝑛2 𝜃 𝜕𝜙 2 𝑟 sin 𝜃 𝜕𝜙 𝑟 2 𝑠𝑖𝑛2 𝜃 𝜕𝜙
+ 𝜌𝑔𝜙
Sample Problem:
Consider a Newtonian fluid flowing in a circular tube at constant temperature as shown on figure below. The fluid
(in laminar flow) is in steady-state flow and has a fully developed velocity profile.
a. What is the velocity profile across the tube?
b. Also, find the shear stress profile
*Since there is only a 𝜐𝑍 (velocity in the axial direction) and no pressure gradients for both components r and 𝜃,
component r and 𝜃 = 0
For r-component:
𝜕𝜐𝑟 𝜕𝜐𝑟 𝜐𝜃 𝜕𝜐𝑟 𝜐𝜃 2 𝜕𝜐𝑟 𝜕𝑃 𝜕 1 𝜕(𝑟𝜐𝑟 ) 1 𝜕 2 𝜐𝑟 2 𝜕𝜐𝜃 𝜕 2 𝜐𝑟
𝜌( + 𝜐𝑟 + − + 𝜐𝑧 )= − +𝜇[ ( )+ 2 − + ] + 𝜌𝑔𝑟 = 0
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
For 𝜃-component:
𝜕𝜐𝜃 𝜕𝜐𝜃 𝜐𝜃 𝜕𝜐𝜃 𝜐𝑟 𝜐𝜃 𝜕𝜐𝜃 1 𝜕𝑃 𝜕 1 𝜕(𝑟𝜐𝜃 ) 1 𝜕 2 𝜐𝜃 2 𝜕𝜐𝑟 𝜕 2 𝜐𝜃
𝜌( + 𝜐𝑟 + + + 𝜐𝑧 )=− +𝜇[ ( )+ 2 + + ] + 𝜌𝑔𝜃 = 0
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝑟 𝜕𝜃 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
For z-component:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 1 𝜕 𝜕𝜐𝑧 1 𝜕 2 𝜐𝑧 𝜕 2 𝜐𝑧
𝜌( + 𝜐𝑟 + + 𝜐𝑧 )= − +𝜇[ (𝑟 )+ 2 + ] + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2
𝜕𝜌 1 𝜕 1 𝜕 𝜕
+ (𝜌𝑟𝜐𝑟 ) + (𝜌𝜐𝜃 ) + (𝜌𝜐𝑧 ) = 0
𝜕𝑡 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
𝜕𝜌
= 0 (𝑆𝑡𝑒𝑎𝑑𝑦 𝑠𝑡𝑎𝑡𝑒 𝑎𝑛𝑑 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡)
𝜕𝑡
1 𝜕
(𝜌𝑟𝜐𝑟 ) = 0 (𝑁𝑜 𝜐𝑟 )
𝑟 𝜕𝑟
1 𝜕
(𝜌𝜐𝜃 ) = 0 (𝑁𝑜 𝜐𝜃 )
𝑟 𝜕𝜃
𝜕
Therefore, (𝜌𝜐𝑧 ) = 0
𝜕𝑧
𝜕𝜐 𝜕2 𝜐
Since 𝜌 is constant, 𝜕𝑧𝑧 = 0 𝑎𝑙𝑠𝑜, 𝜕𝑧 𝑧 = 0
Using the z-component of Equation of Motion in Cylindrical Coordinates in terms of velocity gradient,
𝜕𝜐𝑧
= 0 (𝑆𝑡𝑒𝑎𝑑𝑦 𝑠𝑡𝑎𝑡𝑒)
𝜕𝑡
𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧
𝜐𝑟 𝑎𝑛𝑑 = 0 ( 𝑁𝑜 𝜐𝑟 𝑎𝑛𝑑 𝜐𝜃 )
𝜕𝑟 𝑟 𝜕𝜃
𝜕𝜐𝑧
𝜐𝑧 =0 (𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦)
𝜕𝑧
1 𝜕 2 𝜐𝑧
= 0 (𝑛𝑜 𝜃 𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑒)
𝑟 2 𝜕𝜃 2
𝜕 2 𝜐𝑧
= 0 (𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦)
𝜕𝑧 2
Therefore,
𝜕𝑃 𝜇 𝜕 𝜕𝜐𝑧
0=− + (𝑟 ) + 𝜌𝑔𝑧
𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟
𝜕𝑃 𝜇 𝜕 𝜕𝜐𝑧
− 𝜌𝑔𝑧 = (𝑟 )
𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟
𝜕𝑃 ∆𝑃
𝐿𝑒𝑡 𝑧 = 𝐿 , 𝑡ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒 =
𝜕𝑧 𝐿
∆𝑃 𝜇 𝜕 𝜕𝜐𝑧
− 𝜌𝑔𝑧 = (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
∆𝑃 − 𝜌𝑔𝑧 𝜇 𝜕 𝜕𝜐𝑧
= (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
𝐿𝑒𝑡: ∆𝑃 = 𝑃 − 𝑃̅𝑜
̅ 𝑜 ) − 𝜌𝑔𝑧
(𝑃 − 𝑃 𝜇 𝜕 𝜕𝜐𝑧
= (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
𝐿𝑒𝑡: 𝑃̅𝐿 = 𝑃 − 𝜌𝑔𝑧
̅𝐿 − 𝑃̅𝑜 ) 𝜇 𝜕 𝜕𝜐𝑧
(𝑃
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, = (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
̅0 − 𝑃̅𝐿 )
−(𝑃 𝜇 𝜕 𝜕𝜐𝑧
𝑜𝑟 = (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
̅0 − 𝑃̅𝐿 )𝑟
−(𝑃 𝜕𝜐𝑧
𝜕𝑟 = 𝜕 (𝑟 )
𝜇𝐿 𝜕𝑟
*Since only r is involved, we can write the above equation as ordinary differentials
̅0 − 𝑃̅𝐿 )𝑟
−(𝑃 𝜕𝜐𝑧
𝑑𝑟 = 𝑑 (𝑟 )
𝜇𝐿 𝑑𝑟
𝜕𝜐𝑧 𝜕𝜐𝑧
𝑟=𝑟; 𝑟 = 𝑟
𝑑𝑟 𝑑𝑟
𝜕𝜐𝑧 𝜕𝜐
𝑟=𝑟 ̅0 − 𝑃̅𝐿 )𝑟 𝑟 =𝑟 𝑧
−(𝑃 𝑑𝑟 𝑑𝑟 𝜕𝜐𝑧
∫ [ 𝑑𝑟] = ∫ [𝑑 (𝑟 )]
𝑟=0 𝜇𝐿 𝜕𝜐
𝑟 𝑧 =0 𝑑𝑟
𝑑𝑟
̅0 − 𝑃̅𝐿 )𝑟2
−(𝑃 𝜕𝜐𝑧
=𝑟
2𝜇𝐿 𝑑𝑟
̅0 − 𝑃̅𝐿 )𝑟
−(𝑃 𝜕𝜐𝑧
=
2𝜇𝐿 𝑑𝑟
Again integrating and using the boundary conditions,
̅0 − 𝑃̅𝐿 )𝑟
−(𝑃
𝑑𝑟 = 𝜕𝜐𝑧
𝜇𝐿
𝑟=𝑟
−(𝑃̅0 − 𝑃̅𝐿 )𝑟 𝜐𝑧 =𝜐𝑧
∫ [ 𝑑𝑟] = ∫ 𝜕𝜐𝑧
𝑟=𝑅 2𝜇𝐿 𝜐𝑧 =0
̅0 − 𝑃̅𝐿 )
−(𝑃
[𝑟2 − 𝑅2 ] = [𝜐𝑧 − 0]
4𝜇𝐿
̅0 − 𝑃̅𝐿 )
(𝑃
[−𝑟2 + 𝑅2 ] = 𝜐𝑧
4𝜇𝐿
̅𝐋 )𝑹𝟐
̅𝟎 − 𝐏
(𝐏 𝐫 𝟐
Or: 𝛖𝐳 = [𝟏 − (𝐑) ] (𝑽𝒆𝒍𝒐𝒄𝒊𝒕𝒚 𝑫𝒊𝒔𝒕𝒓𝒊𝒃𝒖𝒕𝒊𝒐𝒏)
𝟒𝛍𝐋
̅𝟎 − 𝐏
(𝐏 ̅𝐋 )𝑹𝟐
𝐓𝐡𝐞𝐫𝐞𝐟𝐨𝐫𝐞, 𝐌𝐚𝐱𝐢𝐦𝐮𝐦 𝐕𝐞𝐥𝐨𝐜𝐢𝐭𝐲 (𝛖𝐳 )𝐦𝐚𝐱 =
𝟒𝛍𝐋
𝐫 𝟐
𝐇𝐞𝐧𝐜𝐞, 𝛖𝐳 = (𝛖𝐳 )𝐦𝐚𝐱 [𝟏 − ( ) ] (𝐏𝐚𝐫𝐚𝐛𝐨𝐥𝐢𝐜 𝐞𝐪𝐮𝐚𝐭𝐢𝐨𝐧 𝐟𝐫𝐨𝐦 𝐩𝐫𝐞𝐯𝐢𝐨𝐮𝐬 𝐟𝐢𝐠𝐮𝐫𝐞)
𝐑
(𝐏𝟎 − 𝐏𝐋 )𝑹𝟐 𝟏
𝑨𝒍𝒔𝒐, 𝑨𝒗𝒆𝒓𝒂𝒈𝒆 𝒗𝒆𝒍𝒐𝒄𝒊𝒕𝒚, 𝝊𝒁 𝒂𝒗𝒆 = = (𝝊𝒁 𝒎𝒂𝒙 )
𝟖𝛍𝐋 𝟐
𝐫 𝟐
∗ 𝛖𝐳 = 𝟐 (𝝊𝒁 𝒂𝒗𝒆 ) [𝟏 − ( ) ]
𝐑
𝝅(𝐏𝟎 − 𝐏𝐋 )𝑹𝟒 𝝆
𝑴𝒂𝒔𝒔 𝒇𝒍𝒐𝒘 𝒓𝒂𝒕𝒆 (𝒎̇) =
𝟖𝛍𝐋
b. Shear stress profile2
As in the previous computation , only the Equation of Continuity and the z component of the Equation of Motion
apply.
From the Equation of Continuity,
𝝏𝝊𝒛
=𝟎 (𝑬𝒒𝒖𝒂𝒕𝒊𝒐𝒏 𝒐𝒇 𝒄𝒐𝒏𝒕𝒊𝒏𝒖𝒊𝒕𝒚)
𝝏𝒛
The z-component of Equation of Motion in Cylindrical Coordinates in terms of velocity gradient:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 1 𝜕 𝜕𝜐𝑧 1 𝜕 2 𝜐𝑧 𝜕 2 𝜐𝑧
𝜌( + 𝜐𝑟 + + 𝜐𝑧 )= − +𝜇[ (𝑟 )+ 2 + ] + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2
Using the z-component of Equation of Motion in Cylindrical Coordinates in terms of shear stress, 𝜏:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 ) 1 𝜕𝜏𝜃𝑧 𝜕𝜏𝑧𝑧
𝜌( + 𝜐𝑟 + + 𝜐𝑧 )= − −( + + ) + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
𝜕𝜐𝑧
= 0 (𝑆𝑡𝑒𝑎𝑑𝑦 𝑠𝑡𝑎𝑡𝑒)
𝜕𝑡
𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧
𝜐𝑟 𝑎𝑛𝑑 = 0 ( 𝑁𝑜 𝜐𝑟 𝑎𝑛𝑑 𝜐𝜃 )
𝜕𝑟 𝑟 𝜕𝜃
𝜕𝜐𝑧
𝜐𝑧 =0 (𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦)
𝜕𝑧
𝜕𝜏𝜃𝑧
= 0 (𝑛𝑜 𝜃 𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑒)
𝜕𝜃
𝜕𝜏𝑧𝑧
= 0 (𝑛𝑜 𝑧 𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑒)
𝜕𝑧
Therefore,
𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
0= − − + 𝜌𝑔𝑧
𝜕𝑧 𝑟 𝜕𝑟
𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
𝑜𝑟 0 = + − 𝜌𝑔𝑧
𝜕𝑧 𝑟 𝜕𝑟
Therefore,
𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
− 𝜌𝑔𝑧 = −
𝜕𝑧 𝑟 𝜕𝑟
𝜕𝑃 ∆𝑃
𝐿𝑒𝑡 𝑧 = 𝐿 , 𝑡ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒 =
𝜕𝑧 𝐿
∆𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
− 𝜌𝑔𝑧 = −
𝐿 𝑟 𝜕𝑟
∆𝑃 − 𝜌𝑔𝑧 1 𝜕(𝑟𝜏𝑟𝑧 )
=−
𝐿 𝑟 𝜕𝑟
𝐿𝑒𝑡: ∆𝑃 = 𝑃 − 𝑃̅𝑜
̅ 𝑜 ) − 𝜌𝑔𝑧
(𝑃 − 𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
=−
𝐿 𝑟 𝜕𝑟
𝐿𝑒𝑡: 𝑃̅𝐿 = 𝑃 − 𝜌𝑔𝑧
̅𝐿 − 𝑃̅𝑜 )
(𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, =−
𝐿 𝑟 𝜕𝑟
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜)
𝑟𝑑𝑟 = −𝑑(𝑟𝜏𝑟𝑧 )
𝐿
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜 )𝑟 2
= − (𝑟𝜏𝑟𝑧 + 𝐶1 )
2𝐿
(𝑃̅𝐿 − 𝑃̅𝑜 )𝑟 𝐶1
÷ 𝑟, = − (𝜏𝑟𝑧 + )
2𝐿 𝑟
̅𝑳 − 𝑷
(𝑷 ̅ 𝒐 )𝒓 ̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 )𝒓
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, 𝝉𝒓𝒛 = − =
𝟐𝑳 𝟐𝑳
Exercise:
An incompressible fluid is flowing at steady state in the annular region (i.e., torus or ring between two concentric
cylinders). The coaxial cylinders have an outside radius of 𝑟2 and inner radius of 𝑟1 Find:
a) Shear stress profile
b) Velocity profile
c) Maximum and average velocities
Using the z-component of Equation of Motion in Cylindrical Coordinates in terms of shear stress, 𝜏:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 ) 1 𝜕𝜏𝜃𝑧 𝜕𝜏𝑧𝑧
𝜌( + 𝜐𝑟 + + 𝜐𝑧 )= − −( + + ) + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
𝜕𝜐𝑧
= 0 (𝑆𝑡𝑒𝑎𝑑𝑦 𝑠𝑡𝑎𝑡𝑒)
𝜕𝑡
𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧
𝜐𝑟 𝑎𝑛𝑑 = 0 ( 𝑁𝑜 𝜐𝑟 𝑎𝑛𝑑 𝜐𝜃 )
𝜕𝑟 𝑟 𝜕𝜃
𝜕𝜐𝑧
𝜐𝑧 =0 (𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦)
𝜕𝑧
𝜕𝜏𝜃𝑧
= 0 (𝑛𝑜 𝜃 𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑒)
𝜕𝜃
𝜕𝜏𝑧𝑧
= 0 (𝑛𝑜 𝑧 𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑒)
𝜕𝑧
Therefore,
𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
0= − − + 𝜌𝑔𝑧
𝜕𝑧 𝑟 𝜕𝑟
𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
𝑜𝑟 0 = + − 𝜌𝑔𝑧
𝜕𝑧 𝑟 𝜕𝑟
Therefore,
𝜕𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
− 𝜌𝑔𝑧 = −
𝜕𝑧 𝑟 𝜕𝑟
𝜕𝑃 ∆𝑃
𝐿𝑒𝑡 𝑧 = 𝐿 , 𝑡ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒 =
𝜕𝑧 𝐿
∆𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
− 𝜌𝑔𝑧 = −
𝐿 𝑟 𝜕𝑟
∆𝑃 − 𝜌𝑔𝑧 1 𝜕(𝑟𝜏𝑟𝑧 )
=−
𝐿 𝑟 𝜕𝑟
𝐿𝑒𝑡: ∆𝑃 = 𝑃 − 𝑃̅𝑜
̅ 𝑜 ) − 𝜌𝑔𝑧
(𝑃 − 𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
=−
𝐿 𝑟 𝜕𝑟
𝐿𝑒𝑡: 𝑃̅𝐿 = 𝑃 − 𝜌𝑔𝑧
̅𝐿 − 𝑃̅𝑜 )
(𝑃 1 𝜕(𝑟𝜏𝑟𝑧 )
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, =−
𝐿 𝑟 𝜕𝑟
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜)
𝑟𝑑𝑟 = −𝑑(𝑟𝜏𝑟𝑧 )
𝐿
Integrating and using the boundary conditions,
𝑟 = 𝑟1 ; 𝑟𝜏𝑟𝑧 = 𝑟1 𝜏𝑟𝑧
𝑟 = 𝑟2 ; 𝑟𝜏𝑟𝑧 = 𝑟2 𝜏𝑟𝑧
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜 ) 𝑟=𝑟2 𝑟𝜏𝑟𝑧 =𝑟2 𝜏𝑟𝑧
∫ 𝑟𝑑𝑟 = − ∫ 𝑑(𝑟𝜏𝑟𝑧 )
𝐿 𝑟=𝑟1 𝑟𝜏𝑟𝑧 =𝑟1 𝜏𝑟𝑧
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜 )𝑟2 2 (𝑃 ̅𝐿 − 𝑃 ̅ 𝑜 )𝑟1 2 𝐶1
− = − (𝑟2 𝜏𝑟𝑧 − 𝑟1 𝜏𝑟𝑧 + )
2𝐿 2𝐿 𝑟
𝐵𝑢𝑡 𝐶1 = 0 𝑏𝑒𝑐𝑎𝑢𝑠𝑒 𝝉𝒓𝒛 would be infinite at 𝑟 = 𝑟1 𝑜𝑟 𝑟2 .
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜 )(𝑟2 2 − 𝑟1 2 )
= −(𝑟2 − 𝑟1 )(𝜏𝑟𝑧 )
2𝐿
̅𝐿 − 𝑃
(𝑃 ̅ 𝑜 )(𝑟2 − 𝑟1 )(𝑟2 + 𝑟1 )
= −(𝑟2 − 𝑟1 )(𝜏𝑟𝑧 )
2𝐿
̅ ̅
(𝑃𝐿 − 𝑃𝑜 )(𝑟2 + 𝑟1 )
= −(𝜏𝑟𝑧 )
2𝐿
̅𝑳 − 𝑷
(𝑷 ̅ 𝒐 )(𝒓𝟐 + 𝒓𝟏 ) ̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 )(𝒓𝟐 + 𝒓𝟏 )
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, 𝝉𝒓𝒛 = − =
𝟐𝑳 𝟐𝑳
b. Velocity profile across the tube
For cylindrical coordinates,
Let: z = axial dimension
r = radial dimension
Assume constant 𝜌 𝑎𝑛𝑑 𝜇
*Since there is only a 𝜐𝑍 (velocity in the axial direction) and no pressure gradients for both components r and 𝜃,
component r and 𝜃 = 0
For r-component:
𝜕𝜐𝑟 𝜕𝜐𝑟 𝜐𝜃 𝜕𝜐𝑟 𝜐𝜃 2 𝜕𝜐𝑟 𝜕𝑃 𝜕 1 𝜕(𝑟𝜐𝑟 ) 1 𝜕 2 𝜐𝑟 2 𝜕𝜐𝜃 𝜕 2 𝜐𝑟
𝜌( + 𝜐𝑟 + − + 𝜐𝑧 )= − +𝜇[ ( )+ 2 − + ] + 𝜌𝑔𝑟 = 0
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
For 𝜃-component:
𝜕𝜐𝜃 𝜕𝜐𝜃 𝜐𝜃 𝜕𝜐𝜃 𝜐𝑟 𝜐𝜃 𝜕𝜐𝜃 1 𝜕𝑃 𝜕 1 𝜕(𝑟𝜐𝜃 ) 1 𝜕 2 𝜐𝜃 2 𝜕𝜐𝑟 𝜕 2 𝜐𝜃
𝜌( + 𝜐𝑟 + + + 𝜐𝑧 )=− +𝜇[ ( )+ 2 + + ] + 𝜌𝑔𝜃 = 0
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝑟 𝜕𝜃 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
For z-component:
𝜕𝜐𝑧 𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧 𝜕𝜐𝑧 𝜕𝑃 1 𝜕 𝜕𝜐𝑧 1 𝜕 2 𝜐𝑧 𝜕 2 𝜐𝑧
𝜌( + 𝜐𝑟 + + 𝜐𝑧 )= − +𝜇[ (𝑟 )+ 2 + ] + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2
𝜕𝜌 1 𝜕 1 𝜕 𝜕
+ (𝜌𝑟𝜐𝑟 ) + (𝜌𝜐𝜃 ) + (𝜌𝜐𝑧 ) = 0
𝜕𝑡 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
𝜕𝜌
= 0 (𝑆𝑡𝑒𝑎𝑑𝑦 𝑠𝑡𝑎𝑡𝑒 𝑎𝑛𝑑 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡)
𝜕𝑡
1 𝜕
(𝜌𝑟𝜐𝑟 ) = 0 (𝑁𝑜 𝜐𝑟 )
𝑟 𝜕𝑟
1 𝜕
(𝜌𝜐𝜃 ) = 0 (𝑁𝑜 𝜐𝜃 )
𝑟 𝜕𝜃
𝜕
Therefore, (𝜌𝜐𝑧 ) = 0
𝜕𝑧
𝜕𝜐 𝜕2 𝜐
Since 𝜌 is constant, 𝜕𝑧𝑧 = 0 𝑎𝑙𝑠𝑜, 𝜕𝑧 𝑧 = 0
Using the z-component of Equation of Motion in Cylindrical Coordinates in terms of velocity gradient,
𝜕𝜐𝑧
= 0 (𝑆𝑡𝑒𝑎𝑑𝑦 𝑠𝑡𝑎𝑡𝑒)
𝜕𝑡
𝜕𝜐𝑧 𝜐𝜃 𝜕𝜐𝑧
𝜐𝑟 𝑎𝑛𝑑 = 0 ( 𝑁𝑜 𝜐𝑟 𝑎𝑛𝑑 𝜐𝜃 )
𝜕𝑟 𝑟 𝜕𝜃
𝜕𝜐𝑧
𝜐𝑧 =0 (𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦)
𝜕𝑧
2
1 𝜕 𝜐𝑧
= 0 (𝑛𝑜 𝜃 𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑒)
𝑟 2 𝜕𝜃 2
𝜕 2 𝜐𝑧
= 0 (𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦)
𝜕𝑧 2
Therefore,
𝜕𝑃 𝜇 𝜕 𝜕𝜐𝑧
0=− + (𝑟 ) + 𝜌𝑔𝑧
𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟
𝜕𝑃 𝜇 𝜕 𝜕𝜐𝑧
− 𝜌𝑔𝑧 = (𝑟 )
𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟
𝜕𝑃 ∆𝑃
𝐿𝑒𝑡 𝑧 = 𝐿 , 𝑡ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒 =
𝜕𝑧 𝐿
∆𝑃 𝜇 𝜕 𝜕𝜐𝑧
− 𝜌𝑔𝑧 = (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
∆𝑃 − 𝜌𝑔𝑧 𝜇 𝜕 𝜕𝜐𝑧
= (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
𝐿𝑒𝑡: ∆𝑃 = 𝑃 − 𝑃̅𝑜
̅ 𝑜 ) − 𝜌𝑔𝑧
(𝑃 − 𝑃 𝜇 𝜕 𝜕𝜐𝑧
= (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
𝐿𝑒𝑡: 𝑃̅𝐿 = 𝑃 − 𝜌𝑔𝑧
̅𝐿 − 𝑃̅𝑜 ) 𝜇 𝜕 𝜕𝜐𝑧
(𝑃
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, = (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
̅0 − 𝑃̅𝐿 )
−(𝑃 𝜇 𝜕 𝜕𝜐𝑧
𝑜𝑟 = (𝑟 )
𝐿 𝑟 𝜕𝑟 𝜕𝑟
̅0 − 𝑃̅𝐿 )𝑟
−(𝑃 𝜕𝜐𝑧
𝜕𝑟 = 𝜕 (𝑟 )
𝜇𝐿 𝜕𝑟
*Since only r is involved, we can write the above equation as ordinary differentials
̅𝟎 − 𝑷
−(𝑷 ̅ 𝑳 )𝒓 𝝏𝝊𝒛
𝒅𝒓 = 𝒅 (𝒓 )
𝝁𝑳 𝒅𝒓
**The above equation still holds, but the velocity in the annulus will reach a maximum at some radius r = 𝒓𝒎𝒂𝒙 which is
between 𝒓𝟏 𝒂𝒏𝒅 𝒓𝟐 as shown in the above figure.
̅𝟎 − 𝑷
−(𝑷 ̅ 𝑳 )(𝒓𝒎𝒂𝒙 𝟐 − 𝒓𝟐 ) 𝝏𝝊𝒛
=𝒓
𝟐𝝁𝑳 𝒅𝒓
̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 )(𝒓𝟐 − 𝒓𝒎𝒂𝒙 𝟐 ) 𝝏𝝊𝒛
𝒐𝒓 =𝒓
𝟐𝝁𝑳 𝒅𝒓
̅0 − 𝑃
(𝑃 ̅ 𝐿 )(𝑟 2 − 𝑟𝑚𝑎𝑥 2 ) 𝑑𝑟
𝜕𝜐𝑧 = ( )
2𝜇𝐿 𝑟
𝜐𝑍 =𝜐𝑍 ̅0 − 𝑃
(𝑃 ̅ 𝐿 ) 𝑟=𝑟1 𝑑𝑟
∫ 𝜕𝜐𝑧 = ∫ (𝑟 2 − 𝑟𝑚𝑎𝑥 2 ) ( )
𝜐𝑍 =0 2𝜇𝐿 𝑟=𝑟 𝑟
𝜐𝑍 =𝜐𝑍 ̅0 − 𝑃
(𝑃 ̅ 𝐿 ) 𝑟=𝑟1 𝑑𝑟 𝑟=𝑟1
𝑑𝑟
∫ 𝜕𝜐𝑧 = [∫ (𝑟 2 ) − ∫ (𝑟𝑚𝑎𝑥 2 )]
𝜐𝑍 =0 2𝜇𝐿 𝑟=𝑟 𝑟 𝑟=𝑟 𝑟
𝜐𝑍 =𝜐𝑍 ̅0 − 𝑃
(𝑃 ̅ 𝐿 ) 𝑟=𝑟1 𝑟=𝑟1
𝑑𝑟
∫ 𝜕𝜐𝑧 = [∫ (𝑟𝑑𝑟) − ∫ (𝑟𝑚𝑎𝑥 2 )]
𝜐𝑍 =0 2𝜇𝐿 𝑟=𝑟 𝑟=𝑟 𝑟
̅0 − 𝑃
(𝑃 ̅ 𝐿) 𝑟2 𝑟1
𝜐𝑍 = [ − 𝑟𝑚𝑎𝑥 2 ln 𝑟 ]
2𝜇𝐿 2 𝑟
̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 ) (𝒓𝟏 𝟐 − 𝒓𝟐 ) 𝒓𝟏
𝝊𝒁 = [ − 𝒓𝒎𝒂𝒙 𝟐 𝐥𝐧 ]
𝟐𝝁𝑳 𝟐 𝒓
Re- Integrating and using the boundary conditions, (for outer radius, 𝒓𝟐 ),
𝑟=𝑟; 𝜐𝑍 = 0
𝑟 = 𝑟2 ; 𝜐𝑍 = 𝜐𝑍
̅0 − 𝑃
(𝑃 ̅ 𝐿 )(𝑟 2 − 𝑟𝑚𝑎𝑥 2 ) 𝑑𝑟
𝜕𝜐𝑧 = ( )
2𝜇𝐿 𝑟
𝜐𝑍 =𝜐𝑍 ̅0 − 𝑃
(𝑃 ̅ 𝐿 ) 𝑟=𝑟2 𝑑𝑟
∫ 𝜕𝜐𝑧 = ∫ (𝑟 2 − 𝑟𝑚𝑎𝑥 2 ) ( )
𝜐𝑍 =0 2𝜇𝐿 𝑟=𝑟 𝑟
𝜐𝑍 =𝜐𝑍 ̅0 − 𝑃
(𝑃 ̅ 𝐿 ) 𝑟=𝑟2 𝑑𝑟 𝑟=𝑟2
𝑑𝑟
∫ 𝜕𝜐𝑧 = [∫ (𝑟 2 ) − ∫ (𝑟𝑚𝑎𝑥 2 )]
𝜐𝑍 =0 2𝜇𝐿 𝑟=𝑟 𝑟 𝑟=𝑟 𝑟
𝜐𝑍 =𝜐𝑍 ̅0 − 𝑃
(𝑃 ̅ 𝐿 ) 𝑟=𝑟2 𝑟=𝑟2
𝑑𝑟
∫ 𝜕𝜐𝑧 = [∫ (𝑟𝑑𝑟) − ∫ (𝑟𝑚𝑎𝑥 2 )]
𝜐𝑍 =0 2𝜇𝐿 𝑟=𝑟 𝑟=𝑟 𝑟
̅0 − 𝑃
(𝑃 ̅ 𝐿) 𝑟2 𝑟2
𝜐𝑍 = [ − 𝑟𝑚𝑎𝑥 2 ln 𝑟 ]
2𝜇𝐿 2 𝑟
̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 ) (𝒓𝟐 𝟐 − 𝒓𝟐 ) 𝒓𝟐
𝝊𝒁 = [ − 𝒓𝒎𝒂𝒙 𝟐 𝐥𝐧 ]
𝟐𝝁𝑳 𝟐 𝒓
Combining the Re-integrated equation for inner radius, 𝑟1 𝑎𝑛𝑑 𝑜𝑢𝑡𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠, 𝑟2
̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 ) (𝒓𝟏 𝟐 − 𝒓𝟐 ) 𝒓𝟏
𝝊𝒁 = [ − 𝒓𝒎𝒂𝒙 𝟐 𝐥𝐧 ]
𝟐𝝁𝑳 𝟐 𝒓
̅𝟎 − 𝑷
(𝑷 ̅ 𝑳 ) (𝒓𝟐 𝟐 − 𝒓𝟐 ) 𝒓𝟐
𝝊𝒁 = [ − 𝒓𝒎𝒂𝒙 𝟐 𝐥𝐧 ]
𝟐𝝁𝑳 𝟐 𝒓
𝟏 𝒓𝟐 𝟐 − 𝒓𝟏 𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, 𝒓𝒎𝒂𝒙 = √ 𝒓𝟐 ( )
𝐥𝐧 ( ) 𝟐
𝒓𝟏
̅0 − 𝑃̅𝐿 )
(𝑃 𝑟 2
𝐹𝑜𝑟 𝑐𝑦𝑙𝑖𝑛𝑑𝑒𝑟, 𝜐𝑧 = [1 − ( ) ]
4𝜇𝐿 𝑅
̅𝟎 − 𝑷
(𝑷 ̅ 𝑳) 𝒓𝒎𝒂𝒙 𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, (𝝊𝒛 )𝒎𝒂𝒙 = [𝟏 − ( ) ]
𝟒𝝁𝑳 𝒓𝟐 − 𝒓𝟏
1 𝑟2 2 − 𝑟1 2
𝑤ℎ𝑒𝑟𝑒: 𝑟𝑚𝑎𝑥 = √ 𝑟 ( )
ln ( 2 ) 2
𝑟1
𝟐
1 𝑟2 2 − 𝑟1 2
√ 𝑟2 ( 2
)
̅𝟎 − 𝑷
̅ 𝑳) ln ( )
(𝑷 𝑟1
(𝝊𝒛 )𝒎𝒂𝒙 = 𝟏−
𝟒𝝁𝑳 (𝒓𝟐 − 𝒓𝟏 )
{ [ ] }
Mass Balance:
In a steady flow/ streamline flow, the rate of mass entering the flow system equals that leaving, as
mass can neither accumulated nor depleted within a flow system under steady conditions
𝒎̇ = 𝝆𝝊𝑺 = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕
Therefore, the mass of fluid entering and leaving the tube per unit time ,
𝒎̇ = 𝝆𝒂 𝝊𝒂 𝑺𝒂 = 𝝆𝒃 𝝊𝒃 𝑺𝒃
Average Velocity:
If the flow through the stream tube lies wholly or in part in which shear stresses exist, the
velocity (𝝊), will vary from point to point across the area 𝑺𝒂 and also 𝝊𝒃 will vary from
point to point across area 𝑺𝒃 . At this instance, local velocity differs with the average
velocity.
𝒔𝒃
𝒎̇ = 𝝆 ∫ 𝝊𝒅𝑺
𝒔𝒂
̅ of the entire stream flowing through cross sectional area S will be,
The average velocity 𝝊
𝒎̇ 𝟏 𝒔𝒃
̅=
𝝊 = ∫ 𝝊𝒅𝑺
𝝆𝑺 𝑺 𝒔𝒂
Therefore,
𝒒
̅=
𝒗 𝒘𝒉𝒆𝒓𝒆: 𝒒 = 𝒗𝒐𝒍𝒖𝒎𝒆𝒕𝒓𝒊𝒄 𝒇𝒍𝒐𝒘 𝒓𝒂𝒕𝒆
𝑺
The First Law of Thermodynamics for a closed system: Constant Mass
∆𝑼 = 𝑸 − 𝑾
In general, ∆𝑬 = 𝑸 − 𝑾
Where: E = all the energy types involved = (internal, kinetic, potential, electric, etc.)
The First Law of Thermodynamics for an open system : Constant Volume and steady state situation
The sum of the kinetic, potential, and flow energies of a fluid particle is constant along a
streamline during steady flow when the compressibility and frictional effects are negligible
𝑷 𝒈 𝝊𝟐 −𝒅𝑾𝒔
∆( + 𝒛 + ) = − 𝑭𝒉 (𝑩𝒆𝒓𝒏𝒐𝒖𝒍𝒍𝒊 𝑬𝒒𝒖𝒂𝒕𝒊𝒐𝒏 𝒘𝒊𝒕𝒉 𝒇𝒓𝒊𝒄𝒕𝒊𝒐𝒏𝒂𝒍 𝒉𝒆𝒂𝒕𝒊𝒏𝒈)
𝝆 𝒈𝒄 𝟐 𝒅𝒎
or
𝒈 𝑷𝟏 𝝊𝟏 𝟐 𝒈 𝑷𝟐 𝝊𝟐 𝟐
𝒁𝟏 ( ) + + + 𝑾𝒔𝒉𝒂𝒇𝒕 = 𝒁𝟐 ( ) + + − 𝑭𝒉
𝒈𝒄 𝝆 𝟐 𝒈𝒄 𝝆 𝟐
𝑷𝟏 𝒈 𝝊𝟏 𝟐 𝑷𝟐 𝒈 𝝊𝟐 𝟐
+ 𝒛 + = + 𝒛 + = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕 (𝒂𝒍𝒐𝒏𝒈 𝒂 𝒔𝒕𝒓𝒆𝒂𝒎𝒍𝒊𝒏𝒆)
𝝆 𝒈𝒄 𝟏 𝟐 𝝆 𝒈𝒄 𝟐 𝟐
𝒐𝒓
𝝊𝟏 𝟐 𝒈 𝝊𝟐 𝟐 𝒈
𝑷𝟏 + 𝝆 + 𝝆 𝒛𝟏 = 𝑷𝟐 + 𝝆 + 𝝆 𝒛𝟐 = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕 (𝒂𝒍𝒐𝒏𝒈 𝒂 𝒔𝒕𝒓𝒆𝒂𝒎𝒍𝒊𝒏𝒆)
𝟐 𝒈𝒄 𝟐 𝒈𝒄
Where:
𝑷 = 𝒔𝒕𝒂𝒕𝒊𝒄 𝒑𝒓𝒆𝒔𝒔𝒖𝒓𝒆 = represents the actual thermodynamic pressure of the fluid. This is the
same as the pressure used in thermodynamics and property tables.
𝝊𝟐
𝝆 𝟐 = 𝒅𝒚𝒏𝒂𝒎𝒊𝒄 𝒑𝒓𝒆𝒔𝒔𝒖𝒓𝒆 =represents the pressure rise when the fluid in motion is brought
to a stop isentropically.
𝒈
𝝆 𝒈 𝒛 = 𝒉𝒚𝒅𝒓𝒐𝒔𝒕𝒂𝒕𝒊𝒄 𝒑𝒓𝒆𝒔𝒔𝒖𝒓𝒆 =which is not pressure in a real sense since its value depends
𝒄
on the reference level selected; it accounts for the elevation effects, i.e., of fluid weight on
pressure.
The sum of the static, dynamic, and hydrostatic pressures is called the total pressure
Injection work + potential energy+kinetic energy = - Shaft work - internal energy + heat.
𝑷 𝑔 𝝊𝟐 −𝒅𝑾𝒔 𝒅𝑸
∆( + 𝒛+ ) = − (∆𝒖 − )
𝝆 𝑔𝑐 𝟐 𝒅𝒎 𝒅𝒎
Where: P= pressure
𝜌 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦
𝑧 = ℎ𝑒𝑖𝑔ℎ𝑡 𝑎𝑏𝑜𝑣𝑒 𝑎 𝑑𝑎𝑡𝑢𝑚 𝑝𝑙𝑎𝑛𝑒
𝜐 = 𝑣𝑒𝑙𝑜𝑐𝑖𝑡𝑦 𝑜𝑓 𝑓𝑙𝑜𝑤
𝑊𝑠 = 𝑠ℎ𝑎𝑓𝑡 𝑤𝑜𝑟𝑘 or pump work
* 𝑢 = 𝑖𝑛𝑡𝑒𝑟𝑛𝑎𝑙 𝑒𝑛𝑒𝑟𝑔𝑦 𝑐ℎ𝑎𝑟𝑔𝑒; 𝑄 = ℎ𝑒𝑎𝑡; 𝑚 = 𝑚𝑎𝑠𝑠
Δ𝑃
= 𝑖𝑛𝑗𝑒𝑐𝑡𝑖𝑜𝑛 𝑤𝑜𝑟𝑘
𝜌
𝑔
Δ𝑧 = 𝑃𝑜𝑡𝑒𝑛𝑡𝑖𝑎𝑙 𝐸𝑛𝑒𝑟𝑔𝑦
𝑔𝑐
Δ𝜐 2
= 𝑘𝑖𝑛𝑒𝑡𝑖𝑐 𝑒𝑛𝑒𝑟𝑔𝑦
2
*All of the quantities in equation (pressure, density, etc.) except internal energy (∆𝒖)
𝒅𝑸
and heat ( ) can be directly measured.
𝒅𝒎
𝒅𝑸
Evaluation of (∆𝒖 − ):
𝒅𝒎
For incompressible fluid (good approximation for most liquids and also for gases under certain conditions):
𝒅𝑸
(∆𝒖 − ) = 𝑭𝒉 ~ 𝐹𝑟𝑖𝑐𝑡𝑖𝑜𝑛 ℎ𝑒𝑎𝑡𝑖𝑛𝑔 (𝒏𝒐𝒏 − 𝒖𝒔𝒆𝒇𝒖𝒍 𝒇𝒐𝒓𝒎)
𝒅𝒎
𝑷 𝒈 𝝊𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, ∆( + 𝒛+ ) = 𝟎 (𝑾𝒊𝒕𝒉𝒐𝒖𝒕 𝒇𝒓𝒊𝒄𝒕𝒊𝒐𝒏𝒂𝒍 𝒉𝒆𝒂𝒕𝒊𝒏𝒈 & 𝒔𝒉𝒂𝒇𝒕 𝒘𝒐𝒓𝒌)
𝝆 𝒈𝒄 𝟐
𝑷 𝒈 𝝊𝟐 −𝒅𝑾𝒔
∆( + 𝒛+ ) = − 𝑭𝒉
𝝆 𝒈𝒄 𝟐 𝒅𝒎
Boundary layer
Bernoulli equation is valid for fluid flow where viscous effects are negligibly small
compared to inertial, gravitational, and pressure effects.
Bernoulli’s Equation cannot be valid for an entire flow field of practical interest. In other
words, we cannot apply the Bernoulli equation everywhere in a flow, no matter how small
the fluid’s viscosity.
It turns out that the equation is reasonable in certain regions of many practical flows. We
refer to such regions as “inviscid regions of flow”. They are regions where net viscous or
frictional forces are negligibly small compared to other forces acting on fluid particles.
It is typically useful in flow regions outside of boundary layers (i.e. very close to solid
walls), where the fluid motion is governed by the combined effects of pressure and gravity
forces.
Limitations of Bernoulli’s Equation:
Bernoulli equation is applicable to steady flow. It should not be used during the transient start-
up and shut-down periods, or during periods of change in the flow conditions.
Frictional/Frictionless flow .Every flow involves some friction, no matter how small, and frictional
effects may or may not be negligible. The situation is complicated even more by the amount of
error that can be tolerated. In general, frictional effects are negligible for short flow sections with
large cross sections, especially at low flow velocities. Frictional effects are usually significant in
long and narrow flow passages, in the wake region downstream of an object, and in diverging flow
sections such as diffusers because of the increased possibility of the fluid separating from the
walls in such geometries.
The Bernoulli equation was derived from a force balance on a particle moving along a
“streamline”. Therefore, the Bernoulli equation is not applicable in a flow section that involves a
pump, turbine, fan, or any other machine or impeller since such devices destroy the streamlines
and carry out energy interactions with the fluid particles. When the flow section considered
involves any of these devices, the energy equation should be used instead to account for the shaft
work input or output.
Incompressible fluid. One of the assumptions used in the derivation of the Bernoulli equation is
that density is constant and thus the flow is incompressible.
No heat transfer. The density of a gas is inversely proportional to temperature, and thus the
Bernoulli equation should not be used for flow sections that involve significant temperature
change such as heating or cooling sections.
Application of Bernoulli’s Equation (Without Frictional heating and shaft work) :
A. Pitot Tube: A tube with a pressure tap at the stagnation point that measures stagnation pressure.
𝑷𝟐 − 𝑷𝟏 𝝊𝟐 𝟐 − 𝝊𝟏 𝟐
+ =𝟎
𝝆 𝟐
𝒈
𝜐1 = √2 ℎ1
𝒈𝒄
B. Pitot-Static Tube : A tube that has both a stagnation pressure tap and several circumferential static
pressure taps and it measures both stagnation and static pressures
𝑷𝟐 − 𝑷𝟏 𝝊𝟐 𝟐 − 𝝊𝟏 𝟐
+ =𝟎
𝝆 𝟐
√𝟐∆𝑷
𝝊𝟏 =
𝝆
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝝊=√
𝝆
𝟏/𝟐 𝟏/𝟐
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝝆 𝟐(𝑷𝟏 − 𝑷𝟐 )
𝝊𝟐 = 𝒐𝒓 𝝊𝟐 = 𝑪𝑽𝒗
𝑺 𝟐 𝑺𝟐 𝟐
𝟏 − 𝟐𝟐 𝝆 (𝟏 − )
[ 𝑺𝟏 ] [ 𝑺𝟏 𝟐 ]
Where: 𝐶𝑉𝑣 = coefficient of discharge for venture meter to compensate for friction heating & non-uniform
flow
= 𝑎 𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛 𝑜𝑓 (𝑅𝑒)
𝑷𝟐 − 𝑷𝟏 𝝊𝟐 𝟐 − 𝝊𝟏 𝟐
+ =𝟎
𝝆 𝟐
𝟏/𝟐
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝝊𝟐 = 𝑪𝑽𝒐
𝑺 𝟐
𝝆 (𝟏 − 𝟐 𝟐 )
[ 𝑺𝟏 ]
Where:
𝐶𝑉𝑜 = coefficient of discharge for orifice meter to compensate for friction heating & non-uniform flow
𝐷2
= 𝑎 𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛 𝑜𝑓 (𝑅𝑒) 𝑎𝑛𝑑
𝐷1
𝑆2 = 𝑜𝑟𝑖𝑓𝑖𝑐𝑒 𝑐𝑟𝑜𝑠𝑠 𝑠𝑒𝑐𝑡𝑖𝑜𝑛𝑎𝑙 𝑎𝑟𝑒𝑎 , 𝑚2
𝑆1 = 𝑝𝑖𝑝𝑒 𝑐𝑟𝑜𝑠𝑠 𝑠𝑒𝑐𝑡𝑖𝑜𝑛𝑎𝑙 𝑎𝑟𝑒𝑎, 𝑚2
𝑁
(𝑃1 − 𝑃2 ) = 𝐷𝑖𝑓𝑓𝑒𝑟𝑒𝑛𝑐𝑒 𝑖𝑛 𝑚𝑎𝑛𝑜𝑚𝑒𝑡𝑒𝑟 𝑓𝑙𝑢𝑖𝑑 𝑙𝑒𝑣𝑒𝑙 ,
𝑚2
𝟐(𝑷𝟏 − 𝑷𝟐 ) 𝒅𝒐
𝝊𝟐 = 𝑪𝑽𝒐 √ 𝒘𝒉𝒆𝒓𝒆: 𝜷 = (𝒅𝒊𝒂𝒎𝒆𝒕𝒆𝒓 𝒓𝒂𝒕𝒊𝒐)
𝝆(𝟏 − 𝜷𝟒 ) 𝑫𝒕
𝐹𝑙𝑜𝑤 𝑟𝑎𝑡𝑒, 𝑞 = 𝜐1 𝑆1 = 𝜐2 𝑆2
The velocity is obtained by assuming no loss, and thus it is the maximum velocity that can occur
at the constriction site. In reality, some pressure losses due to frictional effects are inevitable,
and thus the velocity will be less.
Also, the fluid stream will continue to contract past the obstruction, and the vena contracta area
is less than the flow area of the obstruction. Both losses can be accounted for by incorporating a
correction factor called the discharge coefficient, (𝑪𝑽 )
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝑭𝒍𝒐𝒘𝒓𝒂𝒕𝒆 𝑬𝒒𝒖𝒂𝒕𝒊𝒐𝒏: 𝒒 = 𝑺𝒐 𝑪𝑽𝒐 √
𝝆(𝟏 − 𝜷𝟒 )
𝒅𝒐 𝝅𝒅𝟎 𝟐
𝑾𝒉𝒆𝒓𝒆: 𝜷 = ; 𝑪𝑽 = 𝑫𝒊𝒔𝒄𝒉𝒂𝒓𝒈𝒆 𝑪𝒐𝒆𝒇𝒇𝒊𝒄𝒊𝒆𝒏𝒕 ; 𝑺𝟎 = 𝑺 𝟐 =
𝑫𝒕 𝟒
Venturi Meter
𝑊ℎ𝑒𝑟𝑒: 𝐷2 = 𝑑𝑜 𝑎𝑛𝑑 𝐷1 = 𝐷𝑡
Sample Problem:
The pitot-static tube uses mercury as a manometer fluid. When it is placed in a water flow, the
manometer height reading is h = 21.30 cm. Neglecting yaw and other errors, what is the flow
velocity in m/sec?
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝝊=√
𝝆
𝐠
∆𝐏 = 𝐏𝟐 − 𝐏𝟏 = 𝐑 𝐦 (𝛒𝐦 − 𝛒𝐟 )
𝐠𝐜
m
g 9.81 kg 𝑛𝑒𝑤𝑡𝑜𝑛
∆𝑃 = (𝑃1 − 𝑃2 ) = Rm (ρm − ρf ) = ( sec 2 ) (0.213 𝑚)(13600 − 1000) 3 = 26,328.10
gc 1
𝑘𝑔 𝑚𝑎𝑠𝑠 − 𝑚 m 𝑚2
𝑠𝑒𝑐 2 − 𝑛𝑒𝑤𝑡𝑜𝑛
𝑛𝑒𝑤𝑡𝑜𝑛
𝟐(𝑷𝟏 − 𝑷𝟐 ) 𝟐 (26,328.10 𝑚2 ) 𝒎
𝝊=√ = √ kg
= 𝟕. 𝟐𝟓𝟔
𝝆 1000 3 𝒔𝒆𝒄
m
Sample Problem:
What is the water velocity as measured by an orifice plate (0.06-m diameter) in a 0.305-m-
𝑛
diameter pipe with the measured pressure drop being 75,150 𝑚2 ?
For Orifice meter:
𝟏/𝟐
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝝊𝟐 = 𝑪𝑽𝒐
𝑺 𝟐
𝝆 (𝟏 − 𝟐 𝟐 )
[ 𝑺𝟏 ]
𝑛
∆𝑃 = 𝑃1 − 𝑃2 = 75,150 2
𝑚
𝜋𝐷𝑝𝑖𝑝𝑒 2 𝜋(0.305𝑚)2
𝑆1 = = = 0.0731𝑚2
4 4
𝜋𝐷𝑜𝑟𝑖𝑓𝑖𝑐𝑒 2 𝜋(0.06𝑚)2
𝑆2 = = = 2.827 𝑥 10−3 𝑚2
4 4
𝑘𝑔
𝜌𝑤𝑎𝑡𝑒𝑟 = 1000 3
𝑚
𝐷
𝐸𝑣𝑎𝑙𝑢𝑎𝑡𝑖𝑜𝑛 𝑜𝑓 𝐶𝑉𝑜 : as a 𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛 𝑜𝑓 (𝑅𝑒) 𝑎𝑛𝑑 𝐷2
1
𝐷1 𝜐1 𝜌
𝑅𝑒 =
𝜇
*𝑉1 cannot be determined unless 𝑉2 is known (by using the volumetric flow rate, (𝜐1 𝑆1 = 𝜐2 𝑆2 ). However,
neither 𝜐1 𝑛𝑜𝑟 𝜐2 is known, which gives the initial aspect of an unsolvable problem.
𝐷 𝐷 0.06 𝑚
* However, if we consider the behavior of the curve, where 𝐷2 = 𝐷𝑜𝑟𝑖𝑓𝑖𝑐𝑒 = 0.305 𝑚 = 0.20
1 𝑝𝑖𝑝𝑒
From the graph, Reynolds number (Re) is about 40 to 40,000 with 𝐶𝑉𝑜 𝑟𝑎𝑛𝑔𝑖𝑛𝑔 𝑓𝑟𝑜𝑚 0.60 𝑡𝑜 0.68,
By trial and error / Iteration,
Assume 𝑪𝑽𝒐 = 𝟎. 𝟔𝟓
1 1
2
𝑁 2
2(𝑃1 − 𝑃2 ) 2 (75,150 2 )
𝜐2 = 𝐶𝑉𝑜 = 0.65 [ 𝑚 ]
𝑆 2 𝑘𝑔 (2.827 𝑥 10−3 𝑚2 )2
𝜌 (1 − 2 2 ) 1000 3 [1 − ]
[ 𝑆1 ] 𝑚 (0.0731𝑚2 )2
𝑁 1/2
150,300
𝑉𝜐2 = 0.65 [ 𝑚2 ] = 7.97 𝑚
998.504 𝑠𝑒𝑐
𝐵𝑢𝑡 𝜐1 𝑆1 = 𝜐2 𝑆2
𝑆2 𝑚 2.827 𝑥 10−3 𝑚2 𝑚
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, 𝜐̅1 = 𝜐̅2 ( ) = 7.97 ( 2 ) = 0.308
𝑆1 𝑠𝑒𝑐 0.0731𝑚 𝑠𝑒𝑐
𝑚 𝑘𝑔
𝐷1 𝜐1 𝜌 0.305 𝑚(0.308 𝑠𝑒𝑐 )(1000 3 )
𝑅𝑒 = = 𝑚 = 93,940
𝜇 𝑘𝑔
0.001 𝑚 − 𝑠𝑒𝑐
𝑭𝒓𝒐𝒎 𝒈𝒓𝒂𝒑𝒉, 𝑪𝒗𝒐 ≈ 𝟎. 𝟔𝟐 ≠ 𝟎. 𝟔𝟓
Assume 𝑪𝑽𝒐 = 𝟎. 𝟔𝟐
1 1
2
𝑁 2
2(𝑃1 − 𝑃2 ) 2 (75,150 2 )
𝜐2 = 𝐶𝑉𝑜 = 0.62 [ 𝑚 ]
𝑆 2 𝑘𝑔 (2.827 𝑥 10−3 𝑚2 )2
𝜌 (1 − 2 2 ) 1000 3 [1 − ]
[ 𝑆1 ] 𝑚 (0.0731𝑚2 )2
𝑁 1/2
150,300
𝜐2 = 0.62 [ 𝑚2 ] = 7.61 𝑚
998.504 𝑠𝑒𝑐
𝐵𝑢𝑡 𝜐1 𝑆1 = 𝜐2 𝑆2
𝑆2 𝑚 2.827 𝑥 10−3 𝑚2 𝑚
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, 𝜐̅1 = 𝜐̅2 ( ) = 7.61 ( 2 ) = 0.294
𝑆1 𝑠𝑒𝑐 0.0731𝑚 𝑠𝑒𝑐
𝑚 𝑘𝑔
𝐷1 𝜐1 𝜌 0.305 𝑚(0.294 𝑠𝑒𝑐 )(1000 3 )
𝑅𝑒 = = 𝑚 = 89670
𝜇 𝑘𝑔
0.001 𝑚 − 𝑠𝑒𝑐
𝑭𝒓𝒐𝒎 𝒈𝒓𝒂𝒑𝒉, 𝑪𝒗𝒐 ≈ 𝟎. 𝟔𝟐
𝒎
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, 𝝊𝟏 = 𝟎. 𝟐𝟗𝟒
𝒔𝒆𝒄
Alternative Solution:
∗ 𝑰𝒏 𝒕𝒉𝒆 𝒂𝒃𝒔𝒆𝒏𝒄𝒆 𝒐𝒇 𝒅𝒂𝒕𝒂, 𝒇𝒐𝒓 𝒉𝒊𝒈𝒉 𝑹𝒆, 𝑪𝑽 ≈ 𝟎. 𝟔𝟏
𝑨𝒔𝒔𝒖𝒎𝒆: 𝑪𝑽 ≈ 𝟎. 𝟔𝟏
𝑁
2(𝑃1 − 𝑃2 ) 2 (75,150 ) 𝑚3
𝑞 = 𝑆𝑜 𝐶𝑉 √
−3
= (2.827 𝑥 10 𝑚2 ) (0.61)√
𝑚2 = 2.116𝑥10−2
𝜌(1 − 𝛽 4 ) kg 𝑠𝑒𝑐
1000 (1 − 0.204 )
m3
𝐹𝑙𝑜𝑤 𝑟𝑎𝑡𝑒, 𝑞 = 𝜐1 𝑆1 = 𝜐2 𝑆2
𝑚3
𝑞 2.116𝑥10−2 𝑚
𝑠𝑒𝑐
𝑞 = 𝜐1 𝑆1 ; 𝜐1 = = 2
= 0.29
𝑆1 𝜋(0.305𝑚) 𝑠𝑒𝑐
4
Validate: 𝐶𝑉 = 0.61
𝑚 kg
𝐷𝜐1 𝜌 (0.305𝑚) (0.29 𝑠𝑒𝑐 ) (1000 m3 )
𝑅𝑒 = = = 8.845𝑥104
𝜇 𝑘𝑔
0.001 𝑚 − 𝑠𝑒𝑐
𝟗𝟏. 𝟕𝟏𝜷𝟐.𝟓
𝑪𝑽 = 𝟎. 𝟓𝟗𝟓𝟗 + 𝟎𝟑𝟏𝟐𝜷𝟐.𝟏 − 𝟎. 𝟏𝟖𝟒𝜷𝟖 +
𝑹𝒆𝟎.𝟕𝟓
(𝑽𝒂𝒍𝒊𝒅 𝒇𝒐𝒓 𝟎. 𝟐𝟓 < 𝜷 < 𝟎. 𝟕𝟓 𝒂𝒏𝒅 𝟏𝟎𝟒 < 𝑹𝒆 < 𝟏𝟎𝟕 )
Sample Problem:
kg kg
The flow rate of methanol at 20℃ (ρ = 788.4 m3 and μ = 5.857 x10−4 m∙s) through a 4 −
cm diameter pipe is to be measured with a 3 − cm diameter orifice meter equipped with a
mercury manometer across the orifice plate. If the differential height of the manometer is read to
be 11 cm, determine the flow rate of methanol through the pipe and the average flow velocity
Given: 𝐷𝑡 = 4 𝑐𝑚 (0.04 𝑚)
Methanol (20℃)
𝒅𝟎 = 𝟑 𝒄𝒎 (𝟎. 𝟎𝟑𝒎)
Required:
a. Flowrate of methanol, q
b. Average velocity, 𝜐
Solution:
𝟐(𝑷𝟏 − 𝑷𝟐 )
𝑭𝒍𝒐𝒘𝒓𝒂𝒕𝒆 𝑬𝒒𝒖𝒂𝒕𝒊𝒐𝒏: 𝒒 = 𝑺𝒐 𝑪𝑽 √
𝝆(𝟏 − 𝜷𝟒 )
𝑑𝑜 0.03
𝑊ℎ𝑒𝑟𝑒: 𝛽 = = = 0.75
𝐷𝑡 0.04
m
𝑔 9.81 kg 𝑛𝑒𝑤𝑡𝑜𝑛
∆𝑃 = (𝑃1 − 𝑃2 ) = ℎ𝑚 (𝜌𝑚2 − 𝜌𝑚1 ) = ( sec 2 ) (0.11𝑚)(13600 − 788.4) 3 = 13,825
𝑔𝑐 𝑘𝑔 𝑚𝑎𝑠𝑠 − 𝑚 m 𝑚2
1
𝑠𝑒𝑐 2 − 𝑛𝑒𝑤𝑡𝑜𝑛
𝜋𝑑0 2 𝜋(0.03𝑚)2
𝑆0 = 𝑆2 = = = 7.07𝑥10−4 𝑚2
4 4
𝑪𝑽 = 𝑫𝒊𝒔𝒄𝒉𝒂𝒓𝒈𝒆 𝑪𝒐𝒆𝒇𝒇𝒊𝒄𝒊𝒆𝒏𝒕
∗ 𝑰𝒏 𝒕𝒉𝒆 𝒂𝒃𝒔𝒆𝒏𝒄𝒆 𝒐𝒇 𝒅𝒂𝒕𝒂, 𝒇𝒐𝒓 𝒉𝒊𝒈𝒉 𝑹𝒆, 𝑪𝑽 ≈ 𝟎. 𝟔𝟏
𝑨𝒔𝒔𝒖𝒎𝒆: 𝑪𝑽 ≈ 𝟎. 𝟔𝟏
𝑛𝑒𝑤𝑡𝑜𝑛
2(𝑃1 − 𝑃2 ) 2 (13,825 ) 𝑚3
𝑞 = 𝑆𝑜 𝐶𝑉 √ = (7.07𝑥10 −4
𝑚 2 )(0.61) √ 𝑚2 = 3.09𝑥10 −3
𝜌(1 − 𝛽 4 ) kg 𝑠𝑒𝑐
788.4 3 (1 − 0.754 )
m
𝐹𝑙𝑜𝑤 𝑟𝑎𝑡𝑒, 𝑞 = 𝜐1 𝑆1 = 𝜐2 𝑆2
𝑚3
𝑞 3.09𝑥10−3 𝑚
𝑠𝑒𝑐
𝑞 = 𝜐1 𝑆1 ; 𝜐1 = = 2
= 2.46
𝑆1 𝜋(0.04𝑚) 𝑠𝑒𝑐
4
Validate: 𝐶𝑉 = 0.61
𝑚 kg
𝐷𝜐1 𝜌 (0.04𝑚) (2.46 𝑠𝑒𝑐 ) (788.4 m3 )
𝑅𝑒 = = = 1.3245𝑥105
𝜇 kg
5.857 x10−4 m ∙ s
𝟗𝟏. 𝟕𝟏𝜷𝟐.𝟓
𝑪𝑽 = 𝟎. 𝟓𝟗𝟓𝟗 + 𝟎𝟑𝟏𝟐𝜷𝟐.𝟏 − 𝟎. 𝟏𝟖𝟒𝜷𝟖 +
𝑹𝒆𝟎.𝟕𝟓
(𝑽𝒂𝒍𝒊𝒅 𝒇𝒐𝒓 𝟎. 𝟐𝟓 < 𝜷 < 𝟎. 𝟕𝟓 𝒂𝒏𝒅 𝟏𝟎𝟒 < 𝑹𝒆 < 𝟏𝟎𝟕 )
Sample Problems:
What is the volumetric flow rate for the tank (h = 3.048 m; exit cross-sectional area of 0.279 m2) shown
in figure below if 𝐹ℎ is taken to be zero.
Given:
ℎ = 3.048; 𝑆2 = 0.279 𝑚2 𝐹ℎ = 0
Required: Volumetric flow rate, q
Solution:
Using Bernoulli Equation,
(Using orderly and systematic approach),
Inlet 1: Liquid surface in the tank
Outlet 2: Exit
𝑷 𝒈 𝝊𝟐 −𝒅𝑾𝒔
∆( + 𝒛 + ) = − 𝑭𝒉 (𝐵𝑒𝑟𝑛𝑜𝑢𝑙𝑙𝑖 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑤𝑖𝑡ℎ 𝑓𝑟𝑖𝑐𝑡𝑖𝑜𝑛𝑎𝑙 ℎ𝑒𝑎𝑡𝑖𝑛𝑔)
𝝆 𝒈𝒄 𝟐 𝒅𝒎
Δ𝑃
= 0 ( 𝐼𝑛𝑙𝑒𝑡 𝑎𝑛𝑑 𝑜𝑢𝑡𝑙𝑒𝑡 𝑎𝑡 𝑎𝑡𝑚𝑜𝑠𝑝ℎ𝑒𝑟𝑖𝑐 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒)
𝜌
−𝑑𝑊𝑠
= 0 (𝑁𝑜 𝑠ℎ𝑎𝑓𝑡 𝑤𝑜𝑟𝑘 𝑝𝑒𝑟𝑓𝑜𝑟𝑚𝑒𝑑)
𝑑𝑚
𝐹ℎ = 0 (𝐹𝑟𝑖𝑐𝑡𝑖𝑜𝑛𝑙𝑒𝑠𝑠)
𝒈 𝝊𝟐
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, ∆ ( 𝒛+ )=𝟎
𝒈𝒄 𝟐
𝒈 𝝊𝟐 𝟐 − 𝝊𝟏 𝟐
(𝒛𝟐 − 𝒛𝟏 ) + =𝟎
𝒈𝒄 𝟐
*Neglect υ1 , since the average velocity from the exit υ2 will be much greater than the velocity υ1 ≈0
𝒈 ̅𝟐 𝟐
𝝊
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, (𝟎 − 𝒉) + =𝟎
𝒈𝒄 𝟐
𝒈 𝒎 𝒎
̅𝟐 = √𝟐 𝒉 = √𝟐(𝟗. 𝟖
𝝊 𝟐
)(𝟑. 𝟎𝟒𝟖 𝒎) = 𝟕. 𝟕𝟕
𝒈𝒄 𝒔𝒆𝒄 𝒔𝒆𝒄
Volumetric Flowrate, q:
𝒒 = 𝝊𝟏 𝑺𝟏 = 𝝊𝟐 𝑺𝟐
𝐦 𝐦𝟑
𝒒 = ( 𝟕. 𝟕𝟕 𝐬𝐞𝐜) (𝟎. 𝟐𝟕𝟗 𝐦𝟐 ) = 𝟐. 𝟏𝟕 𝐬𝐞𝐜
Sample Problem:
Water issues from the nozzle of a horizontal hose-pipe. The hose has an internal diameter of 60 mm and
the nozzle tapers to an exit diameter of 20 mm. If the gauge pressure at the connection between the nozzle
𝑘𝑔
and the pipe is 200 kPa, what is the flowrate? The density of water is 1000 𝑚3 .
Given:
60-mm 20-mm
𝑷 𝝊𝟐
∆( + ) = 𝟎 𝒘𝒉𝒆𝒓𝒆: 𝚫𝝊𝟐 = 𝝊𝑩 𝟐 − 𝝊𝑨 𝟐
𝝆 𝟐
𝑁
where: ∆P = 200,000
𝑚2
𝑁
Δ𝑃 200,000 2 𝑚2
2 𝟐 𝟐
Δ𝜐 = 𝝊𝑩 − 𝝊𝑨 = −2 ( ) = 2 ( 𝑚 ) = 400 ;
𝜌 𝑘𝑔 𝑠𝑒𝑐 2
1000 3
𝑚
𝑈𝑠𝑖𝑛𝑔 𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡𝑦 𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛,
𝑚̇ = 𝜌𝑎 𝜐𝑎 𝑆𝑎 = 𝜌𝑏 𝜐𝑏 𝑆𝑏
𝜋(0.06𝑚)2 𝜋(0.02𝑚)2
𝑆𝐴 = = 2.827𝑥10−3 𝑚2 𝑎𝑛𝑑 𝑆𝐵 = = 3.142𝑥10−4 𝑚2
4 4
𝜐𝐴 𝑆𝐴 𝜐𝐴 (2.827𝑥10−3 𝑚2 )
𝜐𝐵 = = = 8.997 𝜐𝐴
𝑆𝐵 3.142𝑥10−4 𝑚2
𝑚2
𝑆𝑖𝑛𝑐𝑒, 𝝊𝑩 𝟐 − 𝝊𝑨 𝟐 = 400
𝑠𝑒𝑐 2
𝑚2
(8.997 𝜐𝐴 )2 − 𝝊𝑨 𝟐 = 400
𝑠𝑒𝑐 2
𝑚2
400
𝜐𝐴 =
√ 𝑠𝑒𝑐 2 = 2.237 𝑚 ; 𝜐𝐵 = 8.997 (2.237
𝑚
) = 20.126
𝑚
79.946 𝑠𝑒𝑐 𝑠𝑒𝑐 𝑠𝑒𝑐
𝑚 𝑚3
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒, 𝑞 = 𝜐𝑎 𝑆𝑎 = 2.237 (2.827𝑥10−3 𝑚2 ) = 6.324 𝑥10−3
𝑠𝑒𝑐 𝑠𝑒𝑐
Sample Problem:
𝑘𝑔
Water with a density of 998 𝑚3 enters a 50 mm pipe fitting horizontally as shown below at a steady velocity
𝑚 𝑘𝑁
of 1 and a gauge pressure of 100 2. It leaves the fitting horizontally, at the same elevation, at an angle
𝑠 𝑚
of 45° with the entrance direction. The diameter at the outlet is 20 mm. Assuming the fluid density is
constant, the kinetic energy and momentum correction factors at both entrance and exit are unity, and the
friction loss in the fitting is negligible, calculate the gauge pressure at the exit of the fitting
Given:
𝒎
̅𝒂 = 𝟏
𝝊𝟏 = 𝝊
𝒔𝒆𝒄
𝑫𝟏 = 𝟓𝟎 𝒎𝒎 ; 𝑫𝟐 = 𝟐𝟎 𝒎𝒎
𝒌𝒈
𝝆 = 𝟗𝟗𝟖 = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕 𝒅𝒆𝒏𝒔𝒊𝒕𝒚
𝒎𝟑
𝒌𝑵
𝑷𝟏 = 𝟏𝟎𝟎 𝟐
𝒎
Required: 𝑃2𝑔𝑎𝑢𝑔𝑒
Solution:
𝑷 𝑔 𝝊𝟐 −𝒅𝑾𝒔
∆( + 𝒛+ ) = − 𝑭𝒉 (𝐵𝑒𝑟𝑛𝑜𝑢𝑙𝑙𝑖 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑤𝑖𝑡ℎ 𝑓𝑟𝑖𝑐𝑡𝑖𝑜𝑛𝑎𝑙 ℎ𝑒𝑎𝑡𝑖𝑛𝑔)
𝝆 𝑔𝑐 𝟐 𝒅𝒎
−𝑑𝑊𝑠
= 0 (𝑁𝑜 𝑠ℎ𝑎𝑓𝑡 𝑤𝑜𝑟𝑘 𝑝𝑒𝑟𝑓𝑜𝑟𝑚𝑒𝑑)
𝑑𝑚
𝐹ℎ = 0 (𝐹𝑟𝑖𝑐𝑡𝑖𝑜𝑛𝑙𝑒𝑠𝑠)
𝑔 𝑔
∆ ( 𝒛) = 𝟎 or 𝑔𝑐
(𝒛𝟐 − 𝒛𝟏 ) = 𝟎
𝑔𝑐
(Kinetic energy and momentum correction factors at both entrance and exit are unity and fluid
enters and leaves at the same elevation)
𝑷 𝝊𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, ∆ ( + ) = 𝟎
𝝆 𝟐
𝑷𝟐 − 𝑷𝟏 𝝊𝟐 𝟐 − 𝝊𝟏 𝟐
𝐆𝐞𝐧𝐞𝐫𝐚𝐥 𝐄𝐪𝐮𝐚𝐭𝐢𝐨𝐧: + =𝟎
𝝆 𝟐
a. Using the volumetric flow rate equation:
𝒒 = 𝝊𝟏 𝑺𝟏 = 𝝊𝟐 𝑺𝟐
𝜋𝐷1 2 𝜋(0.05𝑚)2
𝑆1 == = 1.963 𝑥 10−3 𝑚2
4 4
𝜋𝐷2 2 𝜋(0.02𝑚)2
𝑆2 = = = 3.1416 𝑥 10−4 𝑚2
4 4
𝒎
𝝊𝟏 𝑺𝟏 (𝟏 𝒔𝒆𝒄)(1.963 𝑥 10−3 𝑚2 ) 𝒎
̅𝒃 =
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆, 𝝊𝟐 = 𝝊 = = 𝟔. 𝟐𝟒𝟖
𝑺𝟐 3.1416 𝑥 10−4 𝑚2 𝒔𝒆𝒄
𝑷𝟐 − 𝑷𝟏 𝝊𝟐 𝟐 − 𝝊𝟏 𝟐
+ =𝟎
𝝆 𝟐
𝑘𝑔 𝑚 2 𝑚 2
2
𝜌(𝜐2 − 𝜐1 2)
𝑘𝑁 998 3 [(6.248 𝑠𝑒𝑐 ) − (1 𝑠𝑒𝑐 ) ] 1𝑘𝑁
𝑚
𝑃2 = 𝑃𝑏 = 𝑃1 − = 100 2 − { 𝑥( )}
2 𝑚 2 1000𝑁
𝑘𝑁 𝑘𝑁 𝑘𝑁
𝑃2 = 𝑃𝑏 = 100 2
− 18.981 2 = 81.019 2
𝑚 𝑚 𝑚
2. A large tank open to the atmosphere is filled with water to a height of 5 m from the outlet tap (See figure below).
A tap near the bottom of the tank is now opened, and water flows out from the smooth and rounded outlet. Determine
the water velocity at the outlet.
Given:
Large tank open to atmosphere
Water height (Z1 ) = 5m
∆𝐏
P2 = Patm ; Therefore, ΔP = (P2 − Patm ) = 0 → =𝟎
𝛒
−𝐝𝐖𝐬
Also, − 𝐅𝐡 = 𝟎 (smooth and rounded outlet and no pump work)
𝐝𝐦
𝒈 ∆𝝊𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆: ∆𝒁 + =𝟎
𝒈𝒄 𝟐
Using Point 2 as the reference point
z2 = 0 z1 = 5 → Δz = (𝑧2 − 𝑧1 ) = (0 − 5)m = −5 m
𝒈 ∆𝝊𝟐
( ∆𝒛 + )=𝟎
𝒈𝒄 𝟐
𝑔 𝑣2 2 − 𝜐1 2
𝑇ℎ𝑒𝑟𝑒𝑓𝑜𝑟𝑒: (𝑧2 − 𝑧1 ) + =0 𝑧2 = 0 𝑎𝑛𝑑 𝜐1 = 0
𝑔𝑐 2
𝐠 𝛖𝟐 𝟐 − 𝟎
(𝟎 − 𝟓)𝐦 + =𝟎
𝐠𝐜 𝟐
𝐠 𝛖𝟐 𝟐
(−𝟓𝐦) +
𝐠𝐜 𝟐
𝐦
𝐠 𝐠 𝟗. 𝟖𝟏 𝐦
𝛖𝟐 𝟐 = 𝟐 (𝟓 𝐦) → 𝛖𝟐 = √𝟐 (𝟓𝐦) = 𝟐 (𝟓𝐦) 𝐬𝐞𝐜𝟐 = 𝟗. 𝟗
𝐠𝐜 𝐠𝐜 𝐤𝐠 𝐦𝐚𝐬𝐬 − 𝐦 𝐬𝐞𝐜
𝟏 𝐤𝐠 𝐦𝐚𝐬𝐬 − 𝐦
𝟐
√ ( 𝐬𝐞𝐜 − 𝐬𝐞𝐜𝟐 )
Sample Problem:
During a trip to the beach (𝑃𝑎𝑡𝑚 = 1 𝑎𝑡𝑚) car runs out of gasoline, and it becomes necessary to siphon
gas out of the car of a Good Samaritan (As shown in figure below). The siphon is a small-diameter hose,
and to start the siphon it is necessary to insert one siphon end in the full gas tank, fill the hose with
gasoline via suction, and then place the other end in a gas can container below the level of the gas
tank. The difference in pressure between point 1 (at the free surface of the gasoline in the tank) and
point 2 (at the outlet of the tube) causes the liquid to flow from the higher to the lower elevation. Point
2 is located 0.75 m below point 1 in this case, and point 3 is located 2 m above point 1. The siphon
diameter is 5 mm, and frictional heating and losses in the siphon is negligible. If the density of gasoline
𝑘𝑔
is 750 , determine:
𝑚3
a. Minimum time to withdraw 4 liters of gasoline from the tank to the gas can container
b. The pressure at point 3
Given:
∆𝑷 𝒈 ∆𝝊𝟐 −𝒅𝑾𝒔
( + ∆𝒛 + )= − 𝑭𝒉
𝝆 𝒈𝒄 𝟐 𝒅𝒎
−𝑑𝑊𝑠
𝐹𝑜𝑟 ( − 𝐹ℎ ) = 0 (𝑛𝑜 𝑝𝑢𝑚𝑝 𝑤𝑜𝑟𝑘 𝑎𝑛𝑑 𝑓𝑟𝑖𝑐𝑡𝑖𝑜𝑛𝑎𝑙 ℎ𝑒𝑎𝑡𝑖𝑛𝑔)
𝑑𝑚
∆𝑃
𝐹𝑜𝑟 ( ) = 0
𝜌
𝑃1 = 𝑃𝑎𝑡𝑚 (𝑂𝑝𝑒𝑛 𝑡𝑜 𝑎𝑡𝑚𝑜𝑠𝑝ℎ𝑒𝑟𝑒) 𝑎𝑛𝑑 𝑃2 = 𝑃𝑎𝑡𝑚 (𝐷𝑖𝑠𝑐ℎ𝑎𝑟𝑔𝑒 𝑡𝑜 𝑡ℎ𝑒 𝑜𝑝𝑒𝑛 𝑔𝑎𝑠 𝑐𝑎𝑛 𝑐𝑜𝑛𝑡𝑎𝑖𝑛𝑒𝑟 )
𝒈 ∆𝝊𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆: ∆𝒛 + =𝟎
𝒈𝒄 𝟐
𝒈 𝒗𝟐 𝟐 − 𝝊𝟏 𝟐
𝑻𝒉𝒆𝒓𝒆𝒇𝒐𝒓𝒆: (𝒛𝟐 − 𝒛𝟏 ) + =𝟎 𝒛𝟐 = 𝟎 𝒂𝒏𝒅 𝝊𝟏 = 𝟎
𝒈𝒄 𝟐
𝒈 𝒗𝟐 𝟐 − 𝟎
(𝟎 − 𝒛𝟏 ) + =𝟎
𝒈𝒄 𝟐
𝒗𝟐 𝟐 𝒈
= (𝒛𝟏 )
𝟐 𝒈𝒄
m
𝑔 9.81 𝑚
𝜐2 = √2 (𝑧1 ) = 2 sec2 (0.75𝑚) = 3.84
𝑔𝑐 𝑘𝑔 𝑚𝑎𝑠𝑠 − 𝑚 𝑠𝑒𝑐
1 𝑘𝑔 𝑚𝑎𝑠𝑠 − 𝑚
2
√ ( 𝑠𝑒𝑐 − 𝑠𝑒𝑐2 )
The cross-sectional area of the tube and the flow rate of gasoline are:
𝜋𝐷𝑡𝑢𝑏𝑒 2 𝜋(0.005𝑚)2
𝑆= = = 1.96 𝑥10−5 𝑚2
4 4
𝑚 𝑚3 𝑙𝑖
𝑞 = 𝜐2 𝑆 = 3.84 (1.96 𝑥10−5 𝑚2 ) = 7.526 𝑥 10−5 (0.07526 )
𝑠𝑒𝑐 𝑠𝑒𝑐 𝑠𝑒𝑐
END OF LECTURE