Introductory Notes and Fluid Properties
Introductory Notes and Fluid Properties
1. History of hydraulics, by Hunter Rouse and Simon Ince, Dover publications, Inc., New York, 1963.
INTRODUCTORY NOTES AND FLUID PROPERTIES 3
I.1.2. 18th and 19th Century Developments. Eighteenth and nineteenth centuries saw a
tremendous advancement in the field of hydraulics and hydrodynamics. It was during this period
that most of the experimental hydraulics was developed at the hands of the French, the Italian and
the German engineers and hydraulicians. Noteworthy amongst them were Poleni, Pitot, Chezy,
Borda, Venturi, Weber, Fourneyron, Belanger, Russel, Reech, Hagen, Poiseuille, Weisbach, Darcy,
Bazin, Kutter, Manning and Froude. It was in the early 18th century when Daniel Bernoulli along
with Leonhard Euler and d Alembert founded the mathematical science of hydrodynamics. Others
who significantly contributed to the development of hydrodynamics are Lagrange, Laplace, Navier,
Stokes, Helmholtz, Kirchhoff, Boussinesq, Reynolds, Thomson (Lord Kelvin), Strutt (Lord Rayleigh)
and Joukowsky. The important contributions of some of these distinguished investigators have
been described in the following paragraphs.
Daniel Bernoulli (1700 1782) was a Swiss mathematician who published a treatise named
Hydrodynamica dealing with various aspects of hydrostatics and hydraulics. He introduced for
the first time the word Hydrodynamics to encompass various topics of fluid statics and dynamics.
He was first to use piezometer openings in the walls of conduits for pressure indication. His energy
principle utilized only two terms, namely, the pressure and the velocity. Bernoulli s principle used
for evaluating pressure seemed to indicate the constancy of pressure and velocity heads. He showed
both experimentally and analytically that the pressure would become negative if the velocity
increased sufficiently.
Leonhard Euler (1707 1783) was also a Swiss mathematician who worked mainly on
hydrodynamics and hydraulic machinery. He was the first to explain the role of fluid pressure in
fluid flow. Euler investigated the motion of a fluid under the action of an external force, rightly
regarding the isotropic pressure as a function of space only ; formulated basic equations of motion
(now known as Euler s equations of acceleration or motion). Assuming the fluid to be incompressible,
flow to be steady and irrotational and utilizing the concept of a force potential he combined his
three equations of acceleration to yield a single relationship involving pressure, velocity and elevation
heads. This equation of Euler in its present familiar form is attributed to Bernoulli. Euler also
contributed significantly to the hydrodynamics theory of centrifugal machinery. Analysing the
performance of reaction turbines, he expressed for the first time the basic relationship by equating
the torque to the change in the moment of momentum as the fluid is passed through the rotating
part.
Jean le Rond d Alembert (1717 1783), was a French mathematician first to introduce
concepts of components of fluid velocity and acceleration and also the differential expression of
continuity. Assuming similar condition in the rear of a body as in the front, the summation of the
elementary pressures exerted on each part of the body surface led to the paradoxical result of a
zero longitudinal force on the body. d Alembert, however, left it on of future investigators to explain
this anomaly between the theory and the reality. This is the paradox of zero resistance to steady
non-uniform motion, known as d Alembert s paradox.
Joseph Louis Lagrange (1736 1813) was a self-trained French mathematician who
succeeded Euler as the world s leading mathematician. In the course of analysis of fluid motion, he
introduced concepts of velocity potential and stream function which became of fundamental
importance in describing the flow pattern. Lagrange was first to derive an equation for velocity of
propagation of a wave of infinitesimal height in channel of finite depth (V = gy ).
Merquis Giovanni Poleni (1683 1761) is the only Italian who deserves a special mention
for his contribution to the experimental hydraulics of the early 18th century. His major contributions
were three-fold : (1) Based on measurements he obtained the coefficient of contraction for a sharp-
edged orifice to be 0.62, a significant improvement over the value of 1/ 2 proposed by Newton. (2)
He also experimented with short tubes (mouth-pieces) attached to the orifice and found that a
maximum rate of discharge was obtained with an intermediate length of tube. (3) His last
contribution was his treatment of the discharge through rectangular sharp- crested weirs. He
considered the discharge as occurring in a series of horizontal elements, the velocity of each being
4 FLUID MECHANICS
assumed proportional to the square root of its depth below the liquid surface. The same approach
was later used for deriving the discharge relationship for sharp-crested weirs, and hence the basic
weir equation is often named after Poleni.
Henri de Pitot (1695 1771) was a French engineer whose principal contribution is the
invention of a device (known after him as the Pitot-tube) for measuring fluid velocity. The original
Pitot tube consisted of two parallel tubes mounted on a slender frame containing a scale and four
petcocks ; one of the tubes being straight and the other bent through 90° at its lower end.
Antoine Chezy (1718 1798) was also a French engineer who contributed significantly to
the understanding of the resistance in uniform open-channel flow. Credit goes to Chezy for not only
presenting the first but also the most lasting resistance formula later to be known by his name.
Jean Charles Borda (1733 1799) was a French military engineer who devoted himself to
experiments in hydraulics and hydraulic machinery. His resistance studies verified the prevalent
theory that the drag of an immersed body varied with the square of the relative velocity and showed
that it would depend upon a still higher power if surface wave were produced. Borda was first to
introduce the concept of elementary stream-tubes. He showed that not only the contraction of jet
but also loss of energy must be taken into account in obtaining an expression for discharge. Making
use of the momentum principle, he found that for the particular case of a re-entrant tube (i.e. the
Borda mouth-piece) the coefficient of contraction has a value of 0.5.
Venturi s published work in Paris in 1797 reported his findings on various forms of
mouthpieces fitted to the orifice. He demonstrated the effect of eddies formed at abrupt changes in
section and, incidentally the change in discharge which would result from their elimination. Venturi
observed that the replacement of the cylindrical tube with two conical sections essentially eliminated
the eddies and increased the rate of flow, but it still produced the local reduction in pressure. This
form of boundary is now used to measure flowrate in pipes and is rightly known as the venturimeter.
Giorgio Bidone (1781 1839) was an Italian hydraulician who is credited with having
discovered the hydraulic jump. He was the first to study it systematically and to attempt its analysis.
Giuseppe Venturoli (1768 1816) also an Italian hydraulician, derived the elementary
back water equation for rectangular channels. Through graphical integration, Venturoli succeeded
in plotting various branches of the surface profile.
Claude Burdin (1790 1893), a French engineer, coined the word turbine and developed
one with free efflux of water.
Benoit Fourneyron (1802 1867) improved upon Burdin s original device and developed a
successful hydraulic turbine. More than hundred similar turbines were built by him for various
parts of the world.
John Scott Russel (1808 1882), a Scottish engineer, was the first to study the problems
of unsteady, non-uniform open-channel flow without discontinuity. He studied the effects of waves
on the resistance of ships and proposed a reverse-curve form of bow which he believed would reduce
the wave effect to a minimum.
A French contemporary of Russel, Ferdinand Reech (1805 1880) advocated the
practicability of model tests and developed similitude principles based on Newton s laws of motion.
Reech was first to express what is now known as the Froude criterion of similitude.
Gotthilf Heinrich Ludwig Hagen (1707 1884), a German hydraulic engineer made
original contributions to the resistance of pipe flow. Based on extensive and accurate experiments
on flow through small diameter tubes, he reasoned that the flow took place in series of cylindrical
layers, the velocity of which varied linearly (an assumption which later proved wrong) from zero at
the wall to a maximum at the centre. Hagen proposed an expression for the resistance to flow in
small diameter tubes (i.e. laminar flow) based on the above assumption. He also carried out
experiments on resistance of pipes in turbulent flow and correlated his measurements by means of
a resistance equation.
INTRODUCTORY NOTES AND FLUID PROPERTIES 5
Jean Louis Poiseuille (1799 1869), a French physician and not an engineer, was interested
in experimental physiology. He conducted research on pumping power of the heart, the movement
of blood in the veins and capillary vessels and the resistance to flow through tubes. He carried out
accurate experiments on very small diameter tubes and presented an empirical relationship for
the discharge in terms of the head loss, the tube diameter and the tube length. The resistance law
for laminar flow was later named after Poiseuille rather than Hagen, and it still continues to be
known by his name.
Julius Weisbach (1806 1871), a German hydraulician, wrote a treatise on hydraulics for
engineering application which even now can be considered a textbook on hydraulics. Weisbach not
only incorporated the best available experimental information, but in many cases supplemented it
with the results of his own experiments. He advocated the use of the non-dimensional coefficient,
and was first to express the resistance equation for pipes in the form h = fLV2/2gD. Weisbach also
modified the weir equation to include the velocity of approach and to eliminate the successive
approximations involved in the determination of discharge he proposed instead and empirical
equation of the type later adopted by Bazin.
Antoine Charles Bresse (1822 1883), a French engineer, accomplished integration of
the equation of gradually varied open-channel flow and prepared tables of function, now known as
Bresse s back water function, for general use. He also presented a correct formulation of the
momentum characteristics of the hydraulic jump.
Henry Philibert Gaspard Darcy (1803 1858), also a French engineer, conducted studies
on the flow of water in both pipes and permeable soils. His experiments included pipes of different
sizes, materials and in various states of deterioration. His greatest contribution was his conclusive
demonstration of the fact that the resistance depended upon the type and the condition of the
boundary material. As a result of his studies on pipe flow, Darcy s name is commonly associated
with that of Weisbach in designating the present-day resistance equation first formulated by
Weisbach. On the basis of his filteration studies, Darcy concluded that the loss of head through a
filter bed was proportional to the rate of flow rather than to its square root as was then generally
believed.
Henri Emile Bazin (1829 1917), a French engineer and an associate of Darcy, conducted
extensive experiments on open-channel resistance, propagation of waves and flow over weirs. Based
upon his experiments in canals of various materials and shapes, Bazin proposed a formula of
resistance in open-channel flow. He also carried out experiments for measurements of velocity
distribution at various cross-sections with different linings. From these tests, he noted that the
depth of point of maximum velocity varied with the relative width of cross-section, approaching
zero as the width-depth ratio exceeded 5. Bazin s subsequent studies on discharge over vertical
and inclined weirs included the precise determination of nappe profile and the distribution of velocity
and pressure through the nappe. As a result of these studies he introduced a new dimensional term
in the weir-discharge equation proposed by Weisbach.
Two Swiss engineers Emile Oscar Ganguillet (1818 1894) and Wilhelm Rudolf Kutter
(1818 1888) are known for their contribution to the open-channel resistance. On the basis of
several hundred experiments, they expressed the Chezy s coefficient C as a function of a roughness
factor n, hydraulic radius R and the channel slope S.
Robert Manning (1816 1897), an Irish engineer, proposed in 1889 a relationship for open-
channel flow of the form, V = K R2/3 S1/2, which was in better agreement with the available data
than the earlier ones. The present-day formula named after Manning was neither recommended
nor ever devised in full by Manning himself. He also did not suggest use of Kutter s n a coefficient
now associated with the Manning s formula.
William Froude (1810 1879), an English engineer, developed and perfected the towing-
tank techniques for testing of model ships. Froude believed that Experiments duly conducted on
small scale model will give results truly indicating of the performance of the full size ships . He
considered the total resistance to be made up of the skin friction resistance and the resistance due
6 FLUID MECHANICS
to other factors such as waves. It is interesting to note that Froude s name has been inseparably
associated with a law of similarity and a non dimensional number, which were neither originated
nor even used by him.
Lousi Marie Henry Navier (1785 1836), a French engineer, derived through a purely
mathematical analysis equations of motion for a viscous flow. His name is, therefore, most frequently
associated with the present day equations for viscous flow. He however, did not identify the fluid
viscosity as a variable affecting the flow but instead considered molecular spacing.
George Gabriel Stokes (1819 1903), a British mathematician, was the one whose name
becomes finally associated with that of Navier in designating the equations of motion for a viscous
fluid. His paper On the Theories of Internal Friction of Fluids in motion published in 1845 contained
the derivations of what are now known as the Navier-Stokes equations. The general coefficient ε
appearing in Navier s equations was replaced by the dynamic viscosity µ. Stokes also derived an
expression for the terminal velocity of fall of spheres, which is now known as the Stokes law.
Osborne Reynolds (1842 1912), a British engineer, was : (1) the first to demonstrate the
phenomenon of cavitation and attribute the accompanying noise to the collapse of vapour bubbles,
just as in a kettle beginning to boil ; (2) the first to correlate the length and time scales in the study
of distorted models ; and (3) the first to introduce the viscosity into a parameter now bearing his
name, demarcating the limit between laminar and turbulent flows. Reynolds demonstrated
experimentally in 1888 that the velocity at which eddy motion (i.e. turbulent flow) began did indeed
vary with the tube diameter and fluid characteristics in such a manner as to yield a fairly definite
value of this parameter (now known as the Reynolds number). The most important contribution of
Reynolds was his application and extension of Navier-Stokes equations to turbulent flow. His lasting
contribution, however, was the derivation of the equations for motion now known as the Reynolds
equations for turbulent flow.
I.1.3. 20th Century Developments Advent of Fluid Mechanics. Until the early 20th
century the two distinctly divergent approaches namely, the experimental hydraulics and the
theoretical hydrodynamics had developed to such an extent that the apparent gap between theory
and the practical reality was bothering the genius amongst the hydraulicians and the
hydrodynamists. The emergence of aeronautics at a rapid pace also hastened the activities of
engineering talent towards bridging the gap between theory and the fact. This was brought about
by a new concept of boundary layer originated by Prandtl of Germany.
Ludwig Prandtl (1875 1953) is regarded as the founder of the present-day fluid mechanics.
He realized the need for better correlation between theory and experiment in problems of fluid
flow. Prandtl conducted his first experiments on the flow of air while he was employed as a engineer
in a large machinery firm. Later he joined as a faculty member of Polytechnic Institute, Hannover
and continued his research on flow of air. Within three years he presented a paper containing his
findings before the Third International Congress of Mathematicians in 1904. In the eight-page
paper he introduced the concept of boundary layer according to which the motion of fluids of low
internal resistance (i.e. low viscosity) can be divided into two interdependent zones : (i) very close
to the fixed boundaries, there is a transition layer in which the fluid velocity changes from zero at
the boundary to practically the same value as in the free-stream at the edge of this layer. It is in
this layer that effects of viscosity are predominant, and (ii) away from the boundary lies the zone of
free-stream across which there is hardly any change in velocity and within which viscous effects
are negligible. The flow in this part can be analysed on the basis of potential flow theory. The layer
in the vicinity of the fixed boundary and to which the viscous effects are confined was given the
name of boundary layer .
The emergence of the theory of boundary layer has come to play a vital role not only in
aeronautics but also in hydraulics and other related fields. This has enabled a much rapid
development of the science of motion of fluids known as the fluid mechanics.
As a result of his pioneering research on flow of low-viscosity fluids, Prandlt was invited as
Professor and Director of Research Institute at the University of Göttingen. At Göttingen, he and
INTRODUCTORY NOTES AND FLUID PROPERTIES 7
his students further worked on the boundary layer theory, analyzed the phenomena of turbulence
and drag and evolved dynamic principles of aerofoil behaviour. The first of his students to achieve
recognition was Paul Richard Heinrich Blasius (1833 ) of Berlin. He published in 1908 an
analytical solution for velocity distribution and resistance of laminar boundary layer wherein
Pandtl s qualitative theory was quantitatively verified by laboratory experiments.
Theoder von Karman (1881 1963) was perhaps the most illustrious alumnus of the
Göttingen Institute where he worked under Prandtl. Karman was gifted with a combination of
rare physical insight and mathematical ability. His primary contributions to fluid mechanics have
been in the fields of form drag turbulence and surface resistance and the analogy between the
sound and gravity waves. He investigated the problem of eddy formation behind circular cylinders
and provided an analytical solution of what has since been known as the Karman vortex trail.
Karman and Prandtl contributed to the analysis of velocity distribution and resistance to turbulent
flow in pipes as well as long flat surfaces. The resulting logarithmic equations for resistance and
velocity distributions are now known by their joint names.
Other notable investigators associated with the Göttingen Institute were : Walter Ludwig
Christian Schiller (1882 ) whose primary interest was in the field of pipe resistance. Walter
Gustav Johannes Tollmein (1900 ) distinguished himself by his analysis of flow stability and
turbulent diffusion ; Hermann Schlichting (1907 ) contributed greatly to the analysis of stability
and boundary layer development. Carl Wieselberger (1887 1941) contributed significantly on
the phenomena of drag ; Otty Flachspart (1898 1957) contributed on drag and particularly on
wind pressures on buildings and Johann Nikuradse (1894 ) contributed in the field of pipe
resistance and is known for his famous experiments on artificially roughened pipes.
Credit goes to Mortiz Weber (1871 1951) of Berlin to put the general principles of similitude
in their present form. It was he who specifically named the Frode and Reynolds numbers associated
Cauchy s name with elastic similarity and introduced a capillarity parameter in his paper presented
in 1919. The capillarity parameter W is subsequently named as the Weber number.
Geoffrey Ingram Taylor (1886 ), a British physicist while employed as a meteorologist
at the University of Cambridge, studied, eddy motion in the atmosphere. He published a series of
papers dealing with the fundamental analysis of fluid turbulence by methods of statistics. Taylor
presented his theory of diffusion by continuous movements and related the diffusive and dissipative
characteristics of turbulent motion.
component of force in the direction of the acceleration is proportional to the product of the mass
and the acceleration. Thus
F ∝ m.a
1
= m.a ...(1)
Cn
where F = force in the direction of acceleration,
m = mass
a = acceleration, and
Cn = constant of proportionality known as Newton s constant.
If we select the units such that one unit of force acting on one unit of mass produces one unit
of acceleration, the proportionality constant Cn will be unity. For this case we can write Newton s
second law in the conventional form, F = m.a, expressing that the force equals mass multiplied by
acceleration. The proportionality constant Cn will be unity only when the absolute system of units
is used, that is when the system of measuring units is coherent.
A coherent system of units may be defined as the one in which the product or quotient of any
two unit quantities involved in the phenomenon is the unit of the resultant quantity. As examples
of this statement, in any coherent system the quotient of unit length and unit time gives unit
velocity and similarly the product of unit mass and unit acceleration gives unit force. Any absolute
system of units is essentially a coherent system.
I.2.2. SI Units2. In the International System (SI) of units the base quantities are mass,
length, time and the thermodynamic temperature and the corresponding base units are kilogram
(kg), metre (m), second (s) and Kelvin (K) respectively. From these base units the derived unit of
force is newton (N) which is the force required to accelerate 1 kg mass at the rate of 1 m/s2, and is
obtained from Eq. (1).
1
1N= × 1 kg (mass) × 1 m/s2
Cn
from which the Newton s constant is
1 kg (mass) m
Cn = ...(2)
N.s2
I.2.3. MKS Gravitational Units. This system has the base quantities as force, length,
time and thermodynamic temperature and the corresponding base units are kilogram (kg), metre
(m), second (s) and Celsius (C). In the gravitational system, unit force is proportional to the product
of unit mass and acceleration due to gravity. Unit force of one kilogram is defined as the force
required to accelerate one kg mass at the rate of 9.81 m/s2. Thus from Eq. (1), we obtain
1
1 kg (force) = × 1 kg (mass) × 9.81 m/s2.
Cn
kg (mass) . m
Cn = 9.81 ...(3)
kg (force) . s2
Comparing Eqs. (2) and (3), 1 kg (force) = 9.81 newton (N) ...(4)
The derived unit of mass is called metric slug in analogy to slug the British gravitational
unit of mass. It is defined as the mass which will be accelerated at the rate of 1 m/s2 when acted
upon by a force of 1 kg. Equation (1) may be used to determine the Newton s constant.
2. The International System of Units Editors Chester H. Page and Paul Vigoureux, London, Her
Majesty s Stationery Office 1973.
INTRODUCTORY NOTES AND FLUID PROPERTIES 9
1
1 kg (force) = × 1 metric slug × 1 m/s2.
Cn
metric slug m
Cn = ...(5)
kg (force) s2
Comparing Eqs. (3) and (5),
1 metric slug = 9.81 kg (mass) ...(6)
Table 1 shows SI and MKS (gravitational) units of various quantities used in fluid mechanics.
System of units
Quantity SI MKS (gravitational)
Unit Symbol Unit Symbol
The conversion factors from MKS to SI, MKS to FPS and MKS to CGS units appear in
Appendix E.
I.3. System
In the most general terms, a system may be defined as that region of space occupied by the
quantity of fluid under consideration. The fluid contained within the system is separated from the
surroundings by a boundary. The system may contain either a constant or a variable mass. Its
boundaries may be fixed or deformable. It may be in motion or at rest with respect to a chosen
co-ordinate system. The region outside the system s boundary is known as the surroundings. A
10 FLUID MECHANICS
system diagram shows the system s boundary and significant interactions between the system and
its surroundings. If the surroundings of importance are forces then the system diagram reduces to
the familiar free-body diagram .
I.3.1. Closed and Open Systems. A system is said to be a closed system if the same body
of fluid remains within the system during a process. The process is then called a non-flow process.
In such a process, work and heat may be transferred across the boundary but no fluid crosses the
boundary. In Fig. 1 (a) dotted lines indicate the systems boundary and only the mass transfer has
been considered.
(1)
m1
Condenser
Reaction
turbine
(1)
Expansion
valve
Compressor System’s Draft
boundary tube
Evaporator
(2) (2)
m2
Refrigeration unit Water turbine
(a) Closed system (a) Open system
Tank wall
D C
FH
2
W E
H
D A
FH
1
A¢ F¢V
F¢H
G
Resultant
reaction F¢ B¢ B
PROPERTIES OF FLUID
1.0. Introduction
Matter can be distinguished by the physical form of its existence. These forms known as
phases, are solid, liquid and gas.
Solid, Liquid and Gas
The liquid and gaseous phases are usually combined and given a common name of fluid,
because of the common characteristics exhibited by liquids and gases. Solids differ from liquids
and liquids from gases on account of their molecular structure (spacing of molecules and the ease
with which they can move). The spacing of molecules is large in a gas, smaller in a liquid and
extremely small in a solid. Very strong intermolecular attractive forces exist in solids which give
them the property of rigidity. These forces are weaker in liquids and extremely small in gases.
Definition of a Fluid
The word fluid means a substance having particles which readily change their relative
positions. A fluid may be defined as substance which deforms continuously under the action of
shear stress, regardless of its magnitude.
Distinction between a Liquid and a Gas
A fluid may be either a liquid or a gas. The molecules of a liquid are very closely spaced as
compared to those of a gas. While a liquid has a free surface and occupies a certain volume in a
container, a gas does not possess a free surface and fills the entire space of the container regardless
of its size. For all practical purposes, a liquid is incompressible while a gas is compressible and
expands unless contained or enclosed in a container. In the words of Sir Oliver Lodge A solid has
volume and shape, a liquid has volume but no shape, a gas has neither. A vapour is a gas whose
temperature and pressure are such that it is very near to the liquid phase, and hence the steam is
12 FLUID MECHANICS
considered as a vapour. A gas may be defined as a superheated vapour. Air is regarded as a gas on
account of its state being normally very far from that of liquid air.
1.2. Density
The density of a substance is defined as the mass per unit volume and is denoted by the
symbol ρ (Greek letter rho). It has the dimension [ML 3]. The expansion or contraction of the
substance results in a change in density. The density of liquids may be considered as constant
while that of gases will be subjected to changes depending on the pressure and temperature. The
fluid density at a point is defined by
∆m
ρ= Lim
∆V → 0 ∆V
where ∆m = mass contained in a small volume ∆V.
Specific Gravity. The specific gravity represents a numerical ratio of two densities, and
water is commonly taken as a reference substance. Thus
Density of the substance
Specific gravity of a substance = .
Density of water
It is also called the relative density.
Mass 174.63/9.81
(ii) Density of mixture = = = 1128 kg/m3.
Volume 15.76 × 10 3
(iii) Volume of mixture = 15.76 litre.
(iv) Weight of mixture = 174.63 N.
3
Example 1.3. If 5 m of a certain oil weighs 40 kN, calculate the specific weight, mass density,
specific volume and relative density of the oil.
Weight 40 (kN)
Solution. Specific weight γ = = = 8 kN/m3
Volume 5 (m 3 )
1 1 3
Specific volume, v = = = 1.23 × 10 m3/kg
ρ 815.49
1.4. Viscosity
Among all the fluid properties, viscosity is the most important and is recognised as the only
single property which influences the fluid motion to a great extent. The viscosity is the property by
virtue of which a fluid offers resistance to deformation under the influence of a shear force. For a
given fluid, the rate of deformation is dependent upon the magnitude of shear force. The molecular
friction or shear resistance within the fluid opposes such continuous deformation.
Let us consider a fluid contained in between two parallel plates as shown in Fig. 1.2, the
bottom one being kept stationary while the top one moves at a constant speed U under influence of
the applied shearing force F.
14 FLUID MECHANICS
Moving plate
1 2 3 4 U
F
Velocity
distribution
b (Linear)
y dq
Stationary
plate
0
Fig. 1.2. Fluid motion between parallel plates, bottom plate stationary, top one moving.
Initially when the top plate is about to start moving all the fluid particles are at rest. The
position of fluid particles lying along the vertical line 0 1 changes with time once the motion
starts. The fluid particles sticking to the moving plate move with the same velocity U while those
adhering to the bottom stationary plate are at rest. The velocity of the intermediate particles vary
from 0 to U. If the gap separating the two plates is small, the velocity distribution will be linear
(straight line) as shown, and the fluid particles originally lying on line 0 1 after a certain time,
say t = ∆t, will occupy the positions indicated by the line 0 2, and at times t = 2∆t and 3∆t they will
lie along the lines 0 3 and 0 4 respectively. The maximum deformation of fluid takes place at
y = b, the magnitude of which in unit time is U, and the zero deformation occurs at y = 0, the bottom
plate being stationary. The time rate of deformation is, therefore, equal to U/b. The rate of angular
deformation dθ/dt is given by
dθ dθ
b =U or y =u
dt dt
dθ U u
∴ = =
dt b y
where u is the velocity at a distance y from the stationary plate.
If A is the area of the moving surface and F is the force required to move the surface at a
constant velocity, it has been established that the shear stress F/A is directly proportional to the
time rate of deformation, thus,
F dθ
∝
A dt
F dθ U u
or =µ =µ =µ
A dt b y
where µ is the proportionality constant and is called coefficient of viscosity. It is also known as the
dynamic viscosity or simply viscosity. The shear stress which is usually denoted by the symbol τ
(Greek letter tau) may be expressed as
U dθ u
τ=µ =µ =µ ...(1.2)
b dt y
Equation (1.2) states that the rate of angular deformation is proportional to the shear stress.
Equation (1.2) in a differential form is expressed as
du
τ=µ ...(1.3)
dy
For a linear velocity distribution, du/dy, is a constant. But if the gap separating the parallel
surfaces is large, the velocity distribution can no longer be assumed linear.
INTRODUCTORY NOTES AND FLUID PROPERTIES 15
Newtonian Non-Newtonian
Time-independent Time-dependent
(i) Dilatant (i) Thixotropic
(ii) Bingham plastic (ii) Rheopectic
(iii) Pseudoplastic
A general relationship between shear stress and velocity gradient (rate of shear strain) for
non-Newtonian fluid may be written as :
τ=A
FG du IJ n
+B ...(1.4)
H dy K
where A and B are constants which depend upon the type of fluid and conditions imposed on the
flow (shear stress).
The fluids which obey Eq. (1.4) are called power-law fluids. The additive constant B is zero
for the fluids except Bingham plastic. Based on the value of power index n in Eq. (1.4), the non-
Newtonian fluids are classified as :
16 FLUID MECHANICS
nt
ing t
ta
B
ila
id
D
flu
Dilatant
an
ty
N
B
Ideal fluid
Pseudoplastic
du/dy
O Velocity gradient log t
du
Fig. 1.4. Plot of τ versus . Fig. 1.5. Plot of log µ versus log τ.
dy
Figs. 1.4 and 1.5 illustrate the shear or viscous characteristics of different fluids. Fig. 1.4
illustrates the shear stress-velocity gradient relationship for various type of fluids. The Newtonian
fluids are characterised by linear relationship between the shear stress and the velocity gradient.
They are represented by a straight line like OA passing through the origin and inclined at an angle
α with the horizontal such that µ = tan α.
The Bingham plastic fluids require a certain minimum shear stress τy known as the yield
stress before they start flowing and exhibit a linear relationship between the shear stress and the
velocity gradient as shown by the straight line PQ. The dilatant and pseudoplastic fluids are shown
by curves marked OC and OB respectively.
Fig. 1.5 illustrates the viscosity µ and shear stress τ relationship for different fluids on log-
log scale. For Newtonian fluids this relationship is a straight line parallel to τ-axis indicating that
viscosity is independent of stress. For non-Newtonian fluids, the fact that the viscosity is a function
of the shear stress can be noticed.
Thixotropic fluids are those which show an increase in apparent viscosity with time. Lipstic
and certain paints and enamels exhibit thixotropic behaviour. The apparent viscosity may be defined
τ
as µapp = . Those fluids which show a decrease in the apparent viscosity with time are called
du/dy
rheopectic. Rheopectic fluids are much less common than thixotropic fluids. Gypsum suspensions
in water and bentonite solutions are examples of rheopectic fluids. Thixotropy is an important
property of paints and enamels. When subjected to high shear by the brush during application of
paint, the apparent viscosity is reduced so that the paint covers the surface smoothly, and brush
marks disappear subsequently.
INTRODUCTORY NOTES AND FLUID PROPERTIES 17
The unit of measurement of viscosity µ, may be obtained using Eq. (1.4 b).
18 FLUID MECHANICS
m2/s. These units can be converted from one system to the other using the same procedure as
above :
2
1 ft2/sec = (30.48) cm2/s = (30.48)2 stokes = 930 stokes.
1
1 stoke = 1
cm2
= FG IJ m /s = 10
1 2 4 m2/s.
sec H K
100
ν= 0.0179
1 + 0.0337 T + 0.000221 T 2
in which ν is the kinematic viscosity in cm2/s and T is the temperature in degree centigrade.
In case of gases, the viscosity is mainly due to transfer of molecular momentum in the
transverse direction brought about by the molecular agitation. The contribution to the viscosity by
the intermolecular cohesive force being negligible due to large spacing of molecules. As the molecular
agitation increases with the rise of temperature the viscosity of gases also increases with temperature
rise.
The viscosity of a fluid may thus be considered to be composed of two parts :
1. that due to intermolecular cohesion, and
2. that due to transfer of molecular momentum, and thus
INTRODUCTORY NOTES AND FLUID PROPERTIES 19
The viscosity of gases, like that of liquids, changes with temperature but is practically
unaffected by pressure. The kinematic viscosity, depending as it does on density, varies with both
temperature and pressure. The following equation, given by Holman, may be used for determining
the value of µ at different temperatures :
µ = 1.7150 × 10 4 (1 + 0.00275 T 3.4 × 10 7 T2) (poise)
in which µ is in poise and T in degrees centigrade. The viscosity in kg s/m2 and Ns/m2 are given by :
µ = 1.7500 × 10 6 (1 + 0.00265 T 3.4 × 10 7 T2) (kg-s/m2
5
µ = 1.715 × 10 (1 + 0.00275 T 3.4 × 10 T ) 7 2 (N-s/m2).
Example 1.4. Calculate the velocity gradient at distances of 0, 100, 150 mm from the boundary
if the velocity profile is a parabola with the vertex 150 mm from the boundary, where the velocity is
1 m/s. Also calculate the shear stresses at these points if the fluid has a viscosity of 0.804 N-s/m2.
Solution. Let the equation of the parabolic velocity profile be
u = Ay2 + By + C ...(1) y
1 m/s
where A, B and C are constants to be determined from the following
boundary conditions :
(i) u = 0 at y = 0
(ii) u = 1 m/s, at y = 0.15 m 150
mm
(iii) du/dy = 0 at the vertex, i.e., y = 0.15 m
Boundary condition (i) gives C = 0 and from (ii), we obtain
1 = A(0.15)2 + B(0.15) ...(2)
u
du O
and from (iii), = 2Ay + B,
dy
or 0 = 2A(0.15) + B ...(3)
Solving Eqs. (2) and (3), A = 44.4 and B = 13.33
Eq. (1) for the velocity profile now becomes
u = 44.4y2 + 13.33y.
20 FLUID MECHANICS
The velocity gradients and the shear stresses the desired points may be obtained as below :
du 1
(a) At y = 0 mm, = 2 × 44.4 × 0 + 13.33 = 13.33 sec
dy
du
Shear stress τ=µ = 0.804 (13.33) = 10.8 N/m2
dy
du 1
(b) At y = 100 mm, = 2 × (44.4) × (0.1) + 13.33 = 4.45 sec
dy
τ = 0.804 (4.45)
= 3.575 N/m2
du
(c) At y = 150 mm, = 2 × (44.4)(0.15) × 15 + 13.33 = 0
dy
du
τ=µ = 0.
dy
Example 1.5. Two horizontal plates are placed 12.5 mm apart, the space between them
being filled with oil of viscosity 14 poise. Calculate the shear stress in the oil if the upper plate moves
with a velocity of 2.5 m/s.
1
Solution. 1 poise = 1 dyne-s/cm2 = kg-s/m2 = 0.1 N-s/m2
98.1
du
Shear stress τ=µ
dy
Relative velocity between the plates, du = 2.5 m/s.
Distance between the plates, dy = 1.25 cm = 0.0125 m
Viscosity of the oil, µ = 14 poise
14
= = 0.143 kg-s/m2 = 1.4 N-s/m2
98.1
(since 1 Poise = 0.1 N-s/m2)
Substituting in the formula,
2.5
τ = 0.143 × = 28.55 kg/m2 = 280.0 N/m2.
0.0125
Example 1.6. A rectangular plate 1.2 m × 0.4 m, weighing 970 N slides down a 45° inclined
surface at a uniform velocity of 2.25 m/s. If the 2 mm gap between the plate and the inclined surface
is filled with oil, determine its viscosity.
Solution. The sliding plate will attain the uniform velocity when the fore causing the motion
(i.e. the component of the plates weight along the inclined surface) balances the fluid resistance
offered by the oil filled in the gap.
du
W sin 45° = µ .A
dy
V Plate 2 mm gap
=µ × 1.2 × 0.4 filled with oil
2 / (10 × 100)
2.25 × 103
or 970 sin 45° = µ × 1.2 × 0.4
/s
W = 970 N
m
2
25
2.
=
970 × 1 / 2 45°
V
Example 1.7. A piston 796 mm diameter and 200 mm long works in a cylinder of 800 mm
diameter. If the annular space is filled with a lubricating oil of viscosity 5 cP, calculate the speed of
descent of piston in vertical position. The weight of piston and the axial load are 9.8 N.
Solution. Viscosity of oil
= 5 cP = 5 × 10 2 Poise
= 0.5 × 10 2 N-s/m2 = 5 × 10 3 N-s/m2.
From Newton s law of viscosity, the shear stress
τ = µ du ...(i) 9.8 N
dy
796
Shear stress developed by the vertical descent of the piston mm
800
= Shear force (i. e. weight of piston etc. ) mm
Piston area in contact with oil
9.8
= = 19.61 N/m2.
π × 0.796 × 0.20
Velocity gradient when the piston attains a constant velocity of V (m/s) in the annular gap of
2 mm
du = V
.
dy 2 × 10 3
Substituting in the Newton s law, Eq. (i)
3 V
19.61 = 5.1 × 10 × 3
2 × 10
from which, V = 7.841 m/s.
Example 1.8. A cylinder of 150 mm radius rotates concentrically inside a fixed cylinder of
155 mm radius. Both cylinders are 300 mm long. Determine the viscosity of the liquid which fills the
space between the cylinders if a torque of 0.98 N-m is required to maintain an angular velocity of
60 r.p.m.
Solution. The torque is transmitted through the fluid layers to the
outer cylinder.
Tangential velocity of the inner cylinder = rω w
2πn
=r
60
2π × 60
= 0.15 × = 0.943 m/s.
60
For the small space between the cylinders, the velocity profile may be assumed to be a
straight line, then
du 0.943
= = 188.6 per sec.
dy (15.5 15.0)
100
Torque applied = Torque resisted
0.98 = τ × Area × Lever arm = τ × (2π × 0.15 × 0.30) × 0.15
0.98
τ= = 23.15 N/m2
2π × 0.045 × 0.15
Example 1.9. A circular disc of a diameter d is slowly rotated in a liquid of large viscosity
µ at a small distance h from a fixed surface. Derive an expression for torque T necessary to maintain
an angular velocity ω .
Solution. Consider an element of disc at a radius r and having a width dr.
Linear velocity at this radius = rω
du
Shear stress τ=µ w
dy
Stationary
Torque = Shear stress × Area × r surface r
= τ × 2πr dr × r dr
= µ du . 2πr2 dr. d
dy
Assuming the gap h to be small so that the velocity distribution may be assumed linear.
du = rω
dy h
∴ Torque, dT, on the element
rω 2πµω 3
dT = µ . 2πr2 dr = r dr
h h
z
d/2
d/2 2πµω 3 2πµω r 4 µπd 4ω
Total torque, T= r dr = = .
0 h h 4 32 h
9
Example 1.10. A space 25 mm wide between two large plane surfaces is filled with glycerine.
What force is required to drag a very thin plate 0.75 sq metre in area between the surfaces at a speed
of 0.5 m/s (i) if this plate remains equidistant from the two surfaces, (ii) if it is at a distance of 10
mm from one of the surfaces ? Take µ = 0.785 N-s/m2.
Solution. Total force required to drag the plate Velocity distribution
=µ
FG du IJ × 0.75 = 0.785 ×
FG 0.5 IJ × 0.75 = 23.5 N
H dy K 1
H 2.5 / 2 × 100 K
The force on the bottom side of the plate
F2 = µ
FG du IJ × Area of plate = 0.785
FG 0.5 IJ × 0.75 = 23.5 N.
H dy K 2
H 2.5/2 × 100 K
∴ The total resistance force (which is equal to the force required to drive it) experienced by
the plate = F1 + F2 = 47.0 N.
Case II. When the plate is located at a distance of 10 mm from one of the surfaces.
Force on the upper side of the plate
FG du IJ × 0.75 = 0.785 FG 0.5 IJ × 5 = 19.60 N.
F1 = µ
H dy K H 1.5 / 100 K
INTRODUCTORY NOTES AND FLUID PROPERTIES 23
τ=µ
FG du IJ = 0.245 × 1.13 × 10 4
H dy K
Clearance
0.0125 cm
= 2.77 × 103 N/m2 = 2.77 kN/m2
Shear force on the shaft
15.2 25
F = τ . 2πrl = 2.77 × 103 × 2 × π × × = 327.5 N
2 × 100 100
Torque to be overcome by the shaft = Shear force × Radius of shaft
T=F.r
15
= 327.5 × = 24.6 N-m.
2 × 100
Power corresponding to this torque at a speed of n = 180 r.p.m.
2πn 2π × 180
P= ×T= × 24.6
60 60
= 463.0 N-m/s = 463.0 W = 0.463 kW.
Example 1.12. A cylinder 0.25 m in radius and 2 m length rotates coaxially inside a fixed
cylinder of the same length and 0.30 m radius. Olive oil of viscosity 4.9 × 10 2 Ns/m2 fills the space
between the cylinders. A torque 4.9 N-m is applied to the inner cylinder. After constant velocity is
attained, calculate the velocity gradient at
T(Resisting)
the cylinder walls, the resulting r.p.m., and Fixed
the power dissipated by fluid resistance cylinder
ignoring end effects. Velocity rad 0.3 m
distribution V
Solution. The surface area of the V
T (Driving)
outer cylinder is larger than that of the inner
one, since the former has a larger radius. X X X X 0.25
m
Accordingly the shear force and the velocity h
gradient at the outer cylinder will be less
than the respective quantities on the inner Olive oil
one. The velocity profile through the fluid
24 FLUID MECHANICS
will be non-linear as indicated in the figure, since the gap between the inner and outer cylinders is
comparatively larger.
The torque of 4.9 Nm is transmitted from inner cylinder to the outer one through fluid
friction (viscous effect). Let r be the radial distance of any fluid layer.
Then 4.9 = τ × (2πrl) × r = τ.2πr × 2 × r
= τ.4πr2 = µ du . 4πr2
dy
2 du
= 4.9 × 10 × 4πr2
dy
du 100
∴ = = 7.95 .
dy 4 πr 2 r2
The velocity gradients at the inner at outer cylinders are :
FG du IJ = 7.95 = 127.2 sec 1
H dy K i (0.25)2
and
FG du IJ 7.95
= 88.3 sec 1.
=
H dy K
0 (0.30)2
Substituting ( dr) for dy in the equation for du/dy since velocity decreases as r increases.
Integrating,
∴
z 0
V
du =
V = 7.95
LM 1OP 0.25
= 5.30 m/s
NrQ 0.30
V 5.3
speed of inner cylinder, ω= = = 21.2 rad/sec.
r 0.25
60ω 60 × 21.2
n= = = 202.4 r.p.m.
2π 2π
Assuming the velocity profile to be linear for an approximate calculation
V = 127.1 × 0.05 = 6.35 m/s
and n = 242.5 r.p.m.
Since this result differs from the former by nearly 20%, the approximation is not satisfactory
in this case.
The power dissipated in fluid friction
2πnT 2π × 202.4 × 4.9
= = = 104.0 Nm/s = 104.0 W.
60 60
Example 1.13. The lower end of a vertical shaft of diameter 10 cm rests in a foot step bearing
(length 100 mm). The clearance between the lower end of the shaft and the bearing surface is 0.5
mm. If the shaft has to run at 750 rpm, find the torque required to keep the shaft in motion. Find
also the power required. Take dynamic viscosity as 1.5 poise. (RGPV, 2013 June)
Solution. Torque required to keep the shaft in motion
= shear force on shaft × its radius
dµ
Shear stress, τ = µ
dy
INTRODUCTORY NOTES AND FLUID PROPERTIES 25
T= zz
2π R
0 0
dT = 2µω × 103 zz
2π R
0 0
r 3 dr dθ
26 FLUID MECHANICS
= 2 × 103 µω
z
0
2π r4
4
R
0
dθ = 2 × 103 µω z0
2π R4
4
dθ
2π R4
= 2 × 103 µω 0 0 4
2µωR4
= . 2π = πR4µω
4
Power reqd. by the shaft bottom = T.ω
= πR4µω2
FG 2π × 750 IJ 2
= π × (5 × 10 2)4 × 0.15 ×
H 60 K Nm/s
= 1.728 N.m
Power dissipated in overcoming viscous friction
2π NT 2π × 1200 × 1.728
P= =
60 60
= 217.06 N.m/s = 217.06 W.
INTRODUCTORY NOTES AND FLUID PROPERTIES 27
Example 1.15. A thin plate of large area is placed midway in a gap of height h filled with oil
of viscosity µ0 and the plate is pulled at a constant velocity V. If a lighter oil of viscosity µ1 is then
substituted in the gap, it is found that for the same velocity V, the drag force will be the same as
before if the plate is located unsymmetrically in the gap but parallel to the walls. Find µ1 in terms of
µ0 and the distance from the nearer wall to the plate.
Solution. Case (i) When the liquid of viscosity µ0
fills the gap. Since the plate is placed midway in the gap, 1
the velocity profile on both the sides of the plate will be 2
h V
symmetrical, thus h V
du V 2V 1 m0
= = . 2
h
dy h/2 h
The shear force on the upper and the bottom side
of the plate will be same, and hence the drag force on the
plate
LM F du I FG du IJ OP A = Lµ 2V 2V OP
4Vµ 0 A
F1 = µ 0
MN GH dy JK H dy K PQ MN A=
Q
+ µ0 0 + µ0
u b
h h h
where A is the area of plate.
Case (ii) When the liquid of viscosity µ1 fills the gap. Let the plate be placed at a distance of
y from the bottom wall as shown.
The velocity gradients for the upper and bottom
sides of the plate are :
FG du IJ =
V
(h–y)
V
m1
H dy K u
(h y) h
V
and
FG du IJ =
V
.
y
H dy K b
y
Now, the drag force on the plate
LM F du I FG du IJ OP A = Lµ V V OP A = µ VhA .
MN GH dy JK PQ MN
1
F2 = µ 1
u
+ µ1
H dy K b
1
(h y)
+ µ1
y Q y(h y)
But since the drag forces F1 and F2 are equal, we have
4Vµ 0 A
= 1
µ VhA
or µ1 = 4µ0
y
1
y LM OP
h y(h y) h h N Q
Example 1.16. A vertical gap 23.5 mm wide of infinite extent contains oil 0.1 m/s
of specific gravity 0.95 and viscosity 2.45 N-s/m2. A metal plate 1.5 m × 1.5 m × 1.5
mm weighing 49 N is to be lifted through the gap at a constant speed of 0.1 m/s.
Estimate the force required.
Solution. Let the plate be placed midway in the gap. The velocity gradient
du 0.1
= = 9.09 sec 1.
dy 1.10 × 100
Viscous resistance to be overcome by the plate 1.175
cm
= Shear force on the plate
= Sum of the shear force acting on each face
2.35 cm
28 FLUID MECHANICS
du
=2.µ .A
dy
= 2 × 2.45 × 9.09 × (1.5 × 1.5)
= 100.1 N
Total force required
= Immersed weight of the plate + Viscous resistance
= (49.0 0.95 × 9810 × 1.5 × 1.5 × 0.0015) + 100.1
= 118.0 N.
Example 1.17. Through a very narrow gap of height h, a thin plate of a very large extent is
being pulled at constant velocity V. On one side of the plate is oil of viscosity µ and on the other side
oil of viscosity Kµ. Calculate the position of the plate so that drag force on it will be a minimum.
Solution. Let the thin plate be placed at a distance y from one of the surfaces as shown. The
drag force per unit area of the plate
F = Sum of shear forces per unit area on both the faces of plate
=µ
FG du IJ + Kµ
FG du IJ
H dy K u
H dy K b
m
where the subscripts u and b refer to the upper and bottom V
sides of the plate, thus h
V
V V y km
F=µ + Kµ
(h y) y
For the drag force to be minimum,
dF
=0
dy
dF 1 V
= µV Kµ =0
dy (h y)2 y2
µ h2 + y2 2 hy h2 h
or = 2
= 2 +1 2
Kµ y y y
or
h2
2
h
+ 1
FG µ IJ = 0.
y2 y H Kµ K
Solving the quadratic equation for h/y
2± 4 4
FG µ IJ
h
=
µH Kµ K =1±
y 2 Kµ
Since h/y cannot be less than unity, using the plus sign
h µ h
=1+ or y=
y Kµ µ
1+
Kµ
Example 1.18. Calculate the approximate viscosity of the oil for the following case :
Solution. When the constant velocity of 0.5 m/s is attained, the viscous resistance to the
motion is equal to the component of the weight of plate along the slope.
INTRODUCTORY NOTES AND FLUID PROPERTIES 29
1 m × 1 m Square plate
ght 13
Wei 50 N 5
/s 1
5m W= q
V = 0. 12
gravity and the surface tension may be taken as 13.55 and 0.53 kg/m (5.2 N/m) respectively. With
h and d both expressed in mm, Eq. (1.6) for glass-air contact becomes
9.6
h= mm ...(1.8)
d
Glass tubes are commonly used for measuring pressure of flow, and in order that they give
correct pressure observations, it is necessary that the rise of water or any other liquid should not
be influenced by the capillary action. To ensure this, the diameter of tube should be large enough
so that the capillary rise is negligible.
Table 1.2 gives the values of h for glass tubes of different diameter for water and mercury as
calculated from the above equation.
Table 1.2
From this table it is obvious that if a smaller tube is used for pressure measurement, the
height of liquid in the tube which indicates pressure of flow will not represent the correct pressure,
as it includes the capillary rise. For this reason, the diameter of tube should never be less than 1
cm.
Liquid drop, jet and soap bubble :
In case of a liquid drop or inside a jet, the action of the surface tension is to increase the
pressure inside, in relation to the outside pressure. In a liquid drop of diameter d, if ∆p is the
difference of pressure between the inside and outside of the drop, then using Fig. 1.8 (a).
∆p . πd2/4 = σπd
4σ
∆p = ...(1.9)
d
and in case of a liquid jet of diameter d and of unit length, we have from Fig. 1.8 (b).
2σ = ∆p . d
2σ
∆p = ...(1.10)
d
A soap bubble in air has two surfaces in contact with air, one inside and the other outside,
Fig. 1.8 (c). The forces that act on the hemispherical section are same as those for the drop, but the
surface tension force is twice as great. The pressure difference is given by
s
Outer
surface
p
Dp
Inner
surface
spd s
(a) Forces on hemispherical (b) Process on half-cylindrical (c) Two surfaces of a
section of liquid drop section of liquid jet soap bubble
Fig. 1.8
32 FLUID MECHANICS
πd 2
∆p . = 2πdσ
4
8σ
or ∆p = ...(1.11)
d
A soap solution has a high value of σ, which causes a soap bubble to be larger in diameter for
small pressure of blowing.
Example 1.19. A soap bubble 25 mm in diameter has inside pressure of 20.0 N/m2 above
atmosphere. Calculate the tension in the soap film.
8σ
Solution. Using Eq. (1.11), ∆p =
d
∆p . d
∴ Surface tension in soap film σ =
8
20
= (25/1000) = 0.0625 N/m.
8
Example 1.20. If the surface tension of water in contact with air is 0.075 N/m, what correction
need be applied toward capillary rise in the manometric reading in tube of 3 mm diameter.
Solution. Assuming the manometer tube made of glass, for water-air-glass contact gibson
has determined the angle of contact θ as 25° 32′, and surface tension σ = 0.0075 kg/m. The capillary
rise is given by Eq. (1.6) as
(4 cos θ) σ
h=
γd
4 cos 25.533° × 0.075
= = 0.009198 m = 9.198 mm
9810 × (3/1000)
However for pure water and clean glass, θ = 0°, and the capillary rise,
4σ
h=
γd
4 × 0.075
= = 0.01019 m = 10.19 mm
9810 × (3/1000)
Note. From the above two values of h, it can be seen that the effect of impure water and
unclean glass is reflected in decrease in the value of capillary rise.
Example 1.21. Calculate the capillary effect in millimetres in a glass tube of 4 mm diameter,
when immersed in (i) water and (ii) mercury. The temperature of the liquid is 20°C and the values of
surface tension of water and mercury at 20°C in contact with air are 0.0736 N/m and 0.51 N/m
respectively. The angle of contact for water is zero and that for mercury 130°.
Solution. The capillary effect is given by Eq. (1.6),
(4 cos θ) σ
h=
γd
For water : σ = 0.0736 N/m γ at 20°C = 9790 N/m3
θ = 0°, d = 4 × 10 3 m
4 cos 0° × 0.0736
h= = 7.51 × 10 3 m
9790 × 4 × 10 3
= 7.51 mm (rise of water)
For mercury : σ = 0.51 N/m γ = 13.6 × 9790 = 133 kN/m3
θ = 130°, d = 4 × 10 3 m
INTRODUCTORY NOTES AND FLUID PROPERTIES 33
= 2.46 mm (depression).
Example 1.22. A U-tube is made of two capillaries of diameter 1.0 mm 1.5 mm respectively.
The tube is kept vertically and partially filled with water of surface tension 0.0736 N/m and zero
contact angle. Calculate the difference in the levels of miniscii caused by the capillarity.
Solution. Capillary rise in a circular tube is given by Eq. (1.6),
4σ cos θ
h= 1 mm 1.5 mm
γd
According to the data given y
θ = 0°, σ = 0.0736 N/m
(i) Capillary rise in 1.0 mm tube
4 × 0.0736 × 1
h1 = = 0.030 m
9810 × 1.0 × 10 2
= 30 mm
(ii) Capillary rise in 1.5 mm tube
4 × 0.0736 × 1
h2 = = 0.020 m
9810 × 1.5 × 10 2
= 20 mm.
∴ The difference in levels of water in the two limbs caused by the surface tension effect
y = h1 h2 = 30 20 = 10 mm.
Example 1.23. The diameters of the two glass limbs of a differential U-tube manometer were
found to be 5 mm and 6 mm respectively. In an experiment the differential pressure readings of 50,
100, 250, 400 and 500 mm were indicated by the manometer. Determine the percentage error caused
by the capillary effect. Take surface tension of water as 0.0736 N/m and the angle of contact as zero.
Solution. Capillary rise in
(i) 5 mm tube (ii) 6 mm tube
4σ cos θ 4σ cos θ 4 × 0.0736 × 1
h1 = h2 = =
γd1 γd2 9810 × 6 × 10 3
4 × 0.0736 × 1
= = 0.006 = 0.005 = 5 mm
9810 × 5 × 10 3
= 6 mm
∴ Capillary effect = h1 h 2
= 6 5 = 1 mm.
Percentage Errors :
1
(i) When pressure difference is 50 mm, % error =± × 100 = ± 2%
50
1
(ii) When pressure difference is 100 mm, % error = ± × 100 = ± 1%
100
1
(iii) When pressure difference is 250 mm, % error = ± × 100 = ± 0.4%
250
1
(iv) When pressure difference is 400 mm, % error = ± × 100 = ± 0.25%
400
1
(v) When pressure difference is 500 mm, % error = ± × 100 = ± 0.20%.
500
34 FLUID MECHANICS
Example 1.24. Determine the absolute pressure and the gauge pressure that would exist
within :
(i) a spherical droplet of water 5 mm in diameter
(ii) a jet of water 5 mm in diameter.
Surface tension of water at the prevalent temperature is 0.0736 N/m and the barometer read-
ing stands at 750 mm of mercury. Take specific gravity of mercury as 13.55 and specific weight of
water as 9810 N/m3.
Solution. Case (i) For a spherical droplet of water, Eq. (1.9) gives differential pressure
within the drop as compared to the outside atmosphere as
∆p = 4σ/d = 4 × 0.0736/5 × 10 3 = 58.86 N/m2
This being the pressure measured above the local atmospheric one, it represents the gauge
pressure.
Local atmospheric pressure = 13.55 × 750 mm water = 9810 × 13.55 × 750 × 10 3 N/m2
= 99.67 kN/m2 = 99.67 kPa.
Absolute pressure inside the droplet
= Gauge pressure + Local atmospheric pressure
= 58.86 + 99.67 × 103 = 99.73 kPa
Case (ii) For the liquid jet, the differential pressure is
2σ 2 × 0.0736
∆p = = = 29.43 N/m2 = 29.43 Pa
d 5 × 10 3
Absolute pressure inside the liquid jet
= 29.43 + 99.67 × 103 N/m2 = 99.7 kN/m2.
Example 1.25. In measuring the unit surface energy of a mineral oil (sp. gr. = 0.85) by the
bubble method, a tube having an internal diameter of 1.5 mm is immersed to a depth of 12.5 mm in
the oil. Air is forced through the tube forming a bubble at the lower end. What magnitude of unit
surface energy will be indicated by a maximum bubble pressure intensity of 147.15 Pa ?
Solution. Specific weight of the mineral oil
γ = 0.85 × 9810 = 8338.5 N/m3
12.5
Pressure at a depth of 12.5 mm = γh = 8338.5 × = 104.2 N/m2 = 104.2 Pa
1000
Pressure difference between inside and outside of the bubble
∆p = 147.15 104.2 = 42.95 Pa
Using Eq. (1.9), the unit surface energy
∆p . d 1.5
σ= =4× = 16.1 N/m.
4 103
Example 1.26. Air is introduced through a nozzle into a tank of water to form a stream of
bubbles. If the bubbles are intended to have a diameter of 2 mm, calculate by how much the pressure
of the air at the nozzle must exceed that of the surrounding water. Assume that surface tension is
71.6 mN/m.
Solution. The pressure at the tip of the nozzle must be
the same as the excess pressure inside the air bubble. The
difference in pressure between inside and outside of the bubble Water
is given by Eq. (1.9),
Bubbles
3
4σ 4 × 71.6 × 10 of 2 mm dia.
∆p = =
d 2 × 10 3 Air
The low vapour pressure of mercury, as evidenced from Table 1.3 (along with its high den-
sity) makes it very suitable for use in barometers and other pressure measuring devices. The vapour
pressure of mercury is so low that there is an almost perfect vacuum above the mercury column in
a barometer.
The problem of cavitation, encountered in hydraulic structures like spillways and sluice
gates, and hydromachinery such as turbines and pumps, is a direct result of local pressures being
equivalent to or less than the vapour pressure of the liquid. When such a situation develops, vapour
bubbles or cavities are formed in the flow. The unsteady nature of these bubbles and their ultimate
collapse is responsible for the high pressure which often leads to vibrations, noise, pitting and
erosion of metal parts of machines and concrete surfaces of hydraulic structures.
Tables 1.5 and 1.6 show the properties of some common fluids and that of water respec-
tively. Table 1.7 deals with the properties of standard atmosphere.
Table 1.5. Properties of some common fluids at 20°C and atmospheric pressure
p = ρRT ...(1.13)
which may be used to compute the density if the gas constant R is known. The absolute tempera-
ture (T) is measured in °K (Kelvin), and is related to temperature in Celcius by the equation
T = 273 + Tc, where Tc = Temperature in °C.
Another fundamental equation for a perfect gas is
pvn = p1v1n = p2v2n = constant ...(1.14)
in which n may have any value from zero to infinity depending upon the process to which the gas is
subjected. By combining Eqs. (1.12), (1.13) and (1.14) the following useful relationship can be ob-
tained :
T2 F I
p
(n 1)/n
Fv I n 1
Fρ I n 1
T1
GH JK
= 2
p1
= GH v JK
1
2
= GH ρ JK
1
2
...(1.15)
Example 1.27. One cubic metre of air at 40°C and pressure 0.105 × 106 N/m2 is compressed
adiabatically to 0.5 m3. What are the temperature and pressure of the gas ? If the process had been
isothermal, what would be the temperature and pressure ?
Solution. The gas constant for air, R = 287 N-m/kg °K
From the equation of state, pv = RT
Substituting the given data,
0.105 × 106 v1 = 287(273 + 40), from which v1 = 1.17 m3/kg
For adiabatic process, pvk = constant
p2(0.5)1.4 = 0.105 × 106(1.17)1.4
p2 = 0.277 × 106 N/m2.
Also, 0.277 × 106(1.17 × 0.5) = 287(273 + T2)
T2 = 292°K.
If the process is isothermal, T2 = 40°C.
pv = constant
0.105 × 106 × 1
and p= = 52.5 × 103 N/m2.
0.5
1.8.1. Compressibility
The compressibility is the measure of change in volume (or density) when a substance is
subjected to pressure. The reciprocal of coefficient of compressibility β is known as the bulk modu-
lus of elasticity. Thus
Coefficient of compressibility = Percentage change in volume for a given change in pressure.
or β = 1/Ev
A fluid may be compressed by the application of pressure, thereby reducing its volume and
giving rise to a volumetric strain. Such a compressed fluid will expand to its original volume when
the applied pressure is withdrawn. This property of compressibility of a fluid is expressed by the
bulk modulus of elasticity. If by applying a pressure dp, the decrease in the fluid volume is dV,
then the bulk modulus of elasticity is defined as
dp
Ev = ...(1.16)
dV /V
where V = original fluid volume. The negative sign indicates a decrease in volume with the increase
in pressure. Since most of the liquids have a comparatively high value of bulk modulus of elasticity,
the compressibility is very close to zero and hence the liquids are considered practically
incompressible under ordinary conditions. The bulk modulus of elasticity of fluid is not a constant
but increases with increasing pressure. Tables 1.5 and 1.6 indicate Ev-values for some common
fluids and water respectively.
As the density is equal to the mass divided by volume, we have
m
ρ=
V
Since the mass m of a certain volume V is constant, differentiating ρ,
dρ = d
FG m IJ = mdV
=
dV
HVK V2
ρ
V
dV dρ
= ...(1.17)
V ρ
From Eqs. (1.16) and (1.17),
dp
Ev = ρ ...(1.18)
dρ
40 FLUID MECHANICS
Example 1.29. A liquid with a volume of 0.2 m3 at 300 kPa is subjected to a pressure of 3000
kPa and its volume is found to decrease by 0.2%. Calculate the bulk modulus of eleasticity of the
liquid.
Solution. Bulk modulus of elasticity of a fluid is defined by Eq. (1.16),
dp
Ev =
dV /V
Pressure increase, dp = 3000 300 = 2700 kPa
0.2 ( 2)
Resulting dcrease in volume, dV/V = 0.2% = =
100 1000
∴ Bulk modulus of the liquid
2700
EV = = 13.5 × 105 kPa.
( 2/1000)
Example 1.30. A cylinder contains 0.75 m3 of gas at 20°C and 2.5 bar pressure. After
compression, the volume gets reduced to 0.15 m3. Determine final pressure and bulk modulus of
compressed gas if compression takes place under :
(i) Isothermal conditions
(ii) Adiabatic conditions (n = 1.4) (RGPV, June 2013)
Solution. (i) Isothermal conditions :
Equation of state, for a perfect gas Eq. (1.12) gives
pv = RT, which for isothermal conditions, T is constunt, hence = constant
or p1v1 = p2v2, where v1 and v2 are volumes per unit weight pressure after compression, i.e.,
v 0.75
p2 = p1 1 = 2.5 × = 12.5 bar
v2 0.15
∴ Final pressure is 12.5 bar.
Bulk modulus of compressed gas, Ev is given by Eq. (1.16) :
dp p2 − p1 12.5 − 0.75
Ev = = = = 12.5 bar
dv/v (v2 − v1)/v1 (0.15 − 0.75)/0.75
1 bar = 105 N/m2 = 105 Pascal
Ev = 12.5 × 105 N/m2
(ii) Adiabatic compression :
For a perfect gas, Eq. (1.14) gives another relationship
pvn = constant
Fv I n
or p 1v 1n = p 2v 2n, ∴ p2 = p1 GH v JK
1
2
0.75 FG IJ 1.4
= 2.5
0.15 H K
= 23.79 bar
12. Match the kinematic viscosity of each liquid choosing the correct value from column B
Col. (A) Col. (B)
Name of liquid Kinematic viscosity in m2/s
(i) Water (a) 150 × 10 7
(ii) Glycerine (b) 10 × 10 7
(iii) Air (c) 6.3 × 10 4.
13. The following table shows the dynamic viscosity of air in N s/m2 at atmospheric pressure at
different temperatures
(i) 2 × 10 5 (ii) 1.90 × 10 5
(iii) 1.81 × 10 5 (iv) 1.715 × 10 5
If the temperatures involved are 0, 20, 40 and 60 degree C, identify each viscosity value
with its temperature.
14. MKS unit of 1 kg force is equal to how many newtons
(a) 0.981 (b) 98.1 (c) 9.81 (d) 981.
2
15. To convert the MKS unit of dynamic viscosity (kg s/m ) into poise, the multiplying factor
is
(i) 89.1 (ii) 981
(iii) 98.1 (iv) 9.81.
16. 1 poise is equal to
(a) 1 dyne s/cm2 (b) 98.1 dyne s/cm2
(c) 1 dyne s/m 2 (d) 1 kg s/m2
(e) 1 kg s/cm2.
17. To convert the MKS unit of viscosity (kg s/m2) into its SI equivalent (N s/m2) multiply by
(i) 98.1 (ii) 981
(iii) 9.81 (iv) 0.981.
18. One stoke is equal to
(a) 1 cm2/s (b) 1 m2/s
(c) 1 ft2/s (d) 1 mm2/s.
2
19. The FPS unit of kinematic viscosity, ft /s, is equal to how many stokes
(i) 93 (ii) 930
(iii) 9.30 (iv) 9300.
20. Surface tension is a phenomenon due to
(a) cohesion only (b) viscous force
(c) adhesion between liquid and solid molecules
(d) difference in magnitude between the forces due to adhesion and cohesion.
21. Weight of liquid that rises in a capillary tube is supported by
(i) the friction between the tube wall and the liquid
(ii) the atmospheric pressure
(iii) the vertical component of force due to surface tension
(iv) the curvature of the miniscus.
22. The capillary depression in mercury is on account of
(a) adhesion being greater than cohesion
(b) surface tension being larger than the viscosity
(c) cohesion being greater than the adhesion
(d) vapour pressure being small.
46 FLUID MECHANICS
REVIEW QUESTIONS
13. Cite examples where surface tension effects play a prominent role.
14. Answer briefly the following questions : (i) define a fluid from mechanics point of view, (ii)
can a fluid sustain (a) tension (b) volumetric compression (c) shear ? In each case state the
relationship between the applied stress and the corresponding shear developed.
15. To what type of flow is the Newton s law of viscosity applicable ? Can it be used to determine
shear stress in turbulent flow ?
16. Classify the fluids on the basis of existence of interface (Free-surface).
17. What are the characteristic fluid properties to which the following phenomena are
attributable, viz :
(i) rise of sap in a trees, (ii) spherical shape of a drop of liquid,
(iii) cavitation, (iv) flow of a jet of oil in an unbroken stream,
(v) water hammer ?
Express the quantities involved in the metric system and also in terms of fundamental
units.
18. Explain the property of fluids on the basis of molecular motion.
19. Define Newtonians fluids.
PROBLEMS
1.1. If the equation of velocity profile is u = 3y2/3 (u in cm/s, y in cm), what is the velocity gradient
at the boundary and at 10 cm from it ?
1.2. A plate 2.5 × 10 5 m distant from a fixed plate moves at 0.60 m/s and requires a force of 1.96
N/m2 to maintain this speed. Determine the fluid viscosity of the substance that fills the
space between the plates.
1.3. A piston 12 cm dia and 15 cm long moves down in a 12.04 cm dia cylinder. The oil filling the
annular space has a viscosity of 8.0 × 10 2 Ns/m2 and the weight of the piston is 9.81 N. Find
the speed with which the piston slides down.
1.4. If two coaxial cylinders 10 cm and 9.75 cm in dia and 25 cm high have a certain liquid filled
in between, find the viscosity of the liquid which produces a torque of 0.98 Nm upon the
inner cylinder when the outer one rotates at the rate of 90 rpm.
2
1.5. If the velocity distribution over a plate is given by u = y y2
3
in which u is the velocity in m/s at a distance y metres above the plate, determine the shear
stress at y = 0 and y = 0.15 m. Take µ = 0.863 Ns/m2.
1.6. Derive the dimensions of :
(i) dynamic viscosity, and (ii) kinematic viscosity
and hence obtain their units in (i) MKS, (iii) FPS and (iii) SI systems.
1.7. A piece of pipe 0.5 m long weighing 9.81 N and having internal diameter of 5.25 cm is
slipped over a vertical shaft 5.0 cm in diameter and allowed to fall. Calculate the approximate
velocity attained by the pipe if a film of oil of viscosity 0.196 Ns/m 2 is maintained between
pipe and shaft.
1.8. A piece of pipe of 5.25 cm internal diameter and 15 cm long slides down a vertical shaft of
5.0 cm diameter at a constant speed of 0.1 m/s. A vertical force 14.7 N is required to pull the
pipe back up the shaft at the same constant speed. Calculate the approximate viscosity of oil
which fills the small gap between the pipe and shaft.
1.9. A Newtonian fluid is filled in the clearance between a shaft and concentric sleeve. When a
force of 490 N is applied to the sleeve parallel to the shaft the sleeve attains a speed of 70
50 FLUID MECHANICS
cm/sec. If 2450 N force is applied, what speed will the sleeve attain ? The temperature of the
sleeve remains constant.
1.10. A very large thin plate is centred in a gap of width 6 cm with different oils of unknown
viscosities above and below, the viscosity of one being twice that of the other. When the
plate is pulled at a velocity of 30 cm/sec, the resulting force on one square metre of plate due
to viscous shear on both sides is 29.4 N. Assuming viscous flow and neglecting all end effects,
calculate viscosities of the oils.
1.11. Two coaxial cylinders with the gap in between completely filled with a viscous fluid, one
cylinder rotating while the other one remains stationary, show that the shear stress on the
inner cylinder is always greater than that on the outer one. Indicate the velocity profiles for
(a) inner cylinder rotating while outer one is stationary, and (b) outer cylinder rotating
while the inner one is kept stationary.
1.12. A rotating viscometer consists of a disc pivoted above a stationary boundary, the fluid to be
tested filling the very small space between the parallel surfaces. Means are available to
measure the driving torque and the rotational speed of the disc. Determine by integration of
the expression for shear stress over the lower surface of the disc., the torque-speed ratio
which would be obtained for a liquid having a viscosity of 0.15 Ns/m2 with a disc having 20
cm diameter and a boundary spacing of 0.19 cm.
1.13. (a) On a plot of shear stress versus velocity gradient represent a Newtonian and non-
Newtonian fluid.
(b) Two large parallel flat plates are placed 1.25 cm apart. A 0.25 cm thick plate of 0.2 m2
area is being towed in glycerine filled between the above plates with a constant force of 9.81
N. Calculate the towing speed of the plate when it is held equidistant from the two parallel
plates. Take µ = 0.01 cm (mass) per cm sec.
1.14. A shaft of diameter 74.9 mm rotates in a bearing of diameter 75 mm and of length 75 mm.
The annular space between the shaft and the bearing is filled with oil having a coefficient of
viscosity of 0.2 stokes and the specific gravity 0.94. Determine the power in overcoming
viscous resistance in this bearing at 1400 rpm.
1.15. Calculate the maximum capillary rise of water (20°C) to be expected in a vertical glass tube
1 mm in diameter. The surface tension at 20°C is 0.0718 N/m.
1.16. Derive an equation for theoretical capillary rise between vertical parallel plates.
1.17. Calculate the maximum rise of water (20°C) to be expected between two vertical, clean glass
plates spaced 1 mm apart.
1.18. What force is necessary to lift a thin wire ring 2.5 cm in diameter from a water surface at
20°C ? Neglect the weight of ring.
1.19. A soap bubble 5 cm in diameter contains a pressure (in excess of atmosphere) of 20.07 N/m2.
Calculate the tension in the soap film.
1.20. Determine the velocity of sound in air at (i) 20°C and 101.2 kN/m2, and (ii) 267°C and 706.0
kN/m2 and find out the ratio of the two velocities.
1.21. Calculate the dynamic viscosity of standard air using table 1.7. Comment upon the effect of
pressure on the dynamic viscosity.
1.22. Compute kinematic viscosity of air at atmospheric pressure using table 1.1. What conclusions
can be drawn regarding the effect of temperature on the viscosity (both dynamic and
kinematic) ?
1.23. A cubical block weighing 196.2 N and having a 20 cm edge is allowed to slide down on an
inclined plane surface making an angle of 20° with the horizontal on which there is a thin
film of oil having a viscosity of 2.16 × 10 3 Ns/m2. What terminal velocity will be attained by
the block, if the film thickness is estimated to be 0.025 mm ?
INTRODUCTORY NOTES AND FLUID PROPERTIES 51
1.24. In a stream of glycerine in motion at a certain point the velocity gradient is 0.25 per second.
The mass density of fluid is 1268 kg per cubic metre and the kinematic viscosity is 6.30 ×
10 4 m2/s. Calculate the shear stress at the point.
1.25. Calculate the capillary effect in millimeter in a glass tube of 4 mm diameter when immersed
in (i) water and (ii) in mercury. The temperature of the liquids is 20°C, and the value of
surface tension of water and mercury at 20°C in contact with air respectively 0.07357 N/m
and 0.490 N/m. The contact angle for θ = 0°, and for mercury θ = 130° 24′.
ANSWERS
1.1. ∞, 0.9275 s 1. 1.2. 81.8 × 10 6 Ns/m2.
1.3. 0.434 m/s. 1.4. 0.7. Ns/m2.
1.5. 0.574 N/m2, 0.316 N/m2.
1.6. (i) FTL 2, kgs/m2, lbs/ft2, Ns/m2. (ii) L2T 1,
m2/s, ft2/s, m2/s,
1.7. 1.515 m/s. 1.8. 3.72 Ns/m2
1.9. 3.5 m/s. 1.10. 0.98 and 1.96 Ns/m2
1.12. 2/3978. 1.13. 1.275 m/s.
1.14. 199.75 kW 1.15. 2.932 cm.
2σ cos θ
1.16. h= 1.17. 1.466 cm.
rwd
1.18. 11.28 × 10 2 N. 1.19. 0.1293 N/m.
1.20. 343 m/s, 465.5 m/s, 1.354.
1.21. 1.18 × 10 6 msl/m-s (1.775 × 10 5 kg/n-s), 1.50 × 10 6 msl/ms (1.50 × 10 6 kg/m-s).
1.22. 13.3 × 10 6 m2/s, 7.85 × 10 6 m2/s. 1.23. 19.4 m/s.
1.24. 19.94 × 10 2 N/m2. 1.25. 7.5 mm in water, and 2.47 mm in mercury.