Lecture Note 1 - Introduction and Fluid Properties
Lecture Note 1 - Introduction and Fluid Properties
Lecture Note 1 - Introduction and Fluid Properties
1
The Course: ME 231
2
University of Tabuk, Mechanical Engineering Department
Assessment Methods: Detailed
3
University of Tabuk, Mechanical Engineering Department
Course References
• Lecture Notes
Other References:
• Engineering Fluid Mechanics by D.F. Elgar, B.C. Williams, C.T. Crowe, J.A.
Roberson, 11th edition.
4
University of Tabuk, Mechanical Engineering Department
ME 231 Course Contents
• Introduction to fluid mechanics and units, properties of fluids, Fluid statics;
Pressure measurements; Buoyancy and stability
• Boundary layer Concepts. Pipe flows, Pipe network and water hammer analysis
Note: These are covered in chapters 1-9 of the reference book
5
University of Tabuk, Mechanical Engineering Department
ME 231 Course Learning Outcomes (CLOs)
• Demonstrate the knowledge of fluid mechanics and units, properties of fluids, Fluid
statics; Pressure measurements and boundary layer Concepts, pipe flows
• Solve static pressure , buoyancy and stability and Fluid flow kinematics problems
6
University of Tabuk, Mechanical Engineering Department
Course Motivation- Application of Fluid Mechanics
Fluid Mechanics have an
unlimited number of
applications the flow in human
body; water flow in households;
household appliances like ACs,
vacuum cleaners, refrigerators;
airplanes; ships; cars; power
plants; are all impossible to
design and operate without the
knowledge of fluid mechanics
7
University of Tabuk, Mechanical Engineering Department
Course Motivation- Application of Fluid Mechanics
The heart is constantly pumping blood through the arteries and veins to all parts of the body (5.6 L total
blood circulating through the body 3 times in one minute; blood in our body travels 19,000 km a day!!),
Lungs are the sites of airflow in alternating directions. All artificial hearts, breathing machines, and
dialysis systems are designed using fluid dynamics
8
University of Tabuk, Mechanical Engineering Department
Fluid Mechanics History
One of the first engineering problems humankind faced as cities were developed was the supply of
water for domestic use and irrigation of crops. Hence fluid mechanics is as old as human civilization.
9
University of Tabuk, Mechanical Engineering Department
What is a Fluid?
Distinction between a solid and a fluid is made on the basis of the
substance’s ability to resist an applied shear (or tangential) stress that tends
to change its shape.
If this experiment were repeated with a fluid (with two large parallel plates
placed in a large body of water, for example), the fluid layer in contact with
the upper plate would move with the plate continuously at the velocity of
the plate no matter how small the force F. The fluid velocity would decrease
with depth because of friction between fluid layers, reaching zero at the
lower plate (no slip Boundary Condition).
Gas molecules are widely spaces with very small cohesive force; hence, gas expands until it encounters
the walls of the container and fills the entire available space.
Any practical fluid system consists of a large number of molecules, and the properties of the system naturally depend
on the behavior of these molecules. Macroscopic approach does not require a knowledge of the behavior of individual
molecules and provides a direct and easy way to analyze engineering problems. Microscopic or statistical approach, is
based on the average behavior of large groups of individual molecules.
11
University of Tabuk, Mechanical Engineering Department
What is Mechanics?
Mechanics is the oldest physical science that deals with both stationary and moving
bodies under the influence of forces. The branch of mechanics that deals with bodies
at rest is called statics, while the branch that deals with bodies in motion is called
dynamics.
Fluid mechanics is defined as the science that deals with the behavior of fluids at rest
(fluid statics) or in motion (fluid dynamics), and the interaction of fluids with solids or
other fluids at the boundaries.
12
University of Tabuk, Mechanical Engineering Department
Branches of Fluid Mechanics
The study of the motion of fluids that can be approximated as incompressible (such as liquids, especially
water, and gases at low speeds) is usually referred to as hydrodynamics.
A subcategory of hydrodynamics is hydraulics, which deals with liquid flows in pipes and open channels.
Gas dynamics deals with the flow of fluids that undergo significant density changes, such as the flow of
gases through nozzles at high speeds.
The category aerodynamics deals with the flow of gases (especially air) over bodies such as aircraft,
rockets, and automobiles at high or low speeds.
Some other specialized categories such as meteorology, oceanography, and hydrology deal with naturally
occurring flows
Note: In all of the above categorizations, we study fluid mechanics by studying either the forces on the
fluid at rest (Hydrostatics), the forces that cause the fluid motion (Dynamics) , or the geometry of the
fluid motion without regard to the forces acting or causing it (Kinematics).
13
University of Tabuk, Mechanical Engineering Department
Classifications of Fluid Flow
There is a wide variety of fluid flow problems encountered in practice, and it is usually convenient
to classify them on the basis of some common characteristics to make it feasible to study them in
groups, some general categories are as follows with bolded ones being simpler to analyze:
• Viscous and Inviscid flow; The internal resistance to flow is quantified by the fluid property
viscosity, inviscid flow have negligible viscosity relative to the flow inertia. Viscosity is caused
by cohesive forces between the molecules in liquids and by molecular collisions in gases.
• Laminar and Turbulent Flow; Some flows are smooth and orderly while others are chaotic. The
highly ordered fluid motion characterized by smooth layers of fluid is called laminar while high
velocity flows characterized by fluctuations are called turbulent flow. Reynold’s number is used
to characterize this flow categories.
• Steady and Unsteady flow; The term steady implies no change of properties with time. The
opposite of steady is unsteady. Uniform flow implies no change of properties with location.
14
University of Tabuk, Mechanical Engineering Department
Classifications of Fluid Flow
• Compressible and incompressible flow; flow is said to be incompressible if the density
remains nearly constant throughout. Subsonic flows at low Mach numbers (V/c <0.3), air can
be treated as incompressible (velocity less than 100 m/s compared to sound speed of 340 m/s)
• Internal and External flow; The flow of an unbounded fluid over a surface such as a plate, a
wire, or a pipe is external flow. The flow in a pipe or duct is internal flow if the fluid is
completely bounded by solid surfaces. Open-channels flow is bounded but with a free surface
• Natural and forced flow; In forced flow, a fluid is forced to flow over a surface or in a pipe by
external means such as a pump or a fan. In natural flows, fluid motion is due to natural means
such as the buoyancy effect (density variation).
• One-, two-, and three-dimensional flow; A flow field is best characterized by its velocity
distribution, and thus a flow is said to be one-, two-, or three-dimensional if the flow velocity
varies in one, two, or three primary dimensions, respectively
15
University of Tabuk, Mechanical Engineering Department
Can you identify these classes of Flow?
16
University of Tabuk, Mechanical Engineering Department
Primary Dimensions and Units
17
University of Tabuk, Mechanical Engineering Department
Derived Dimensions and Units
𝐹 = 𝑚 𝑎
𝐿
𝐹=𝑀 2
𝑇
18
University of Tabuk, Mechanical Engineering Department
Modeling in Engineering
19
University of Tabuk, Mechanical Engineering Department
Problem solving Techniques
• Problem Statement
21
Properties
For the purpose of this course, we can also categorize properties as fluid properties
(the properties of the fluid like density and viscosity) or flow properties (the
characteristics of the flow like velocity and flow rate).
22
University of Tabuk, Mechanical Engineering Department
The Continuum Assumption
The continuum assumption treats a fluid as a continuous, homogeneous matter with no holes
by disregarding the spaces between the individual atoms or molecules that makes up the fluid
(even for gases where the molecules may be widely spaced apart).
The continuum idealization allows us to treat properties as point functions and to assume that
the properties vary continually in space with no jump discontinuities
The continuum idealization is valid as long as the characteristic length of the system (e.g.
diameter) we deal with is large relative to the space between the molecules (mean free path,
MFP). This holds for almost all of the physical flow systems we are interested in.
As an example, the diameter of an oxygen molecule is about 3 x 10-10 m and its mass is 5.3 x
10-26 kg. The MFP of oxygen at 1 atm pressure and 20°C is 6.3 x 10-8 m. Hence, an oxygen
molecule travels, on average, a distance of 6.3 x 10-8 m (about 200 times its diameter) before
it collides with another molecule.
The MFP is so small compared to the size of any real, practical flow!!!!!!!!
23
University of Tabuk, Mechanical Engineering Department
Fluid Properties
Density is defined as mass per unit volume. That is,
𝒎
Density: 𝝆 = Ѵ (kg/m3)
The reciprocal of density is the specific volume v, which is defined as volume per unit mass. That is,
Ѵ 1
ѵ= =
𝑚 𝜌
Sometimes the density of a substance is given relative to the density of a well-known substance. Then it is called
specific gravity, or relative density, and is defined as the ratio of the density of a substance to the density of
some standard substance at a specified temperature (for liquid; usually water at 4°C, for which 𝜌𝐻2𝑂 =
1000 𝑘𝑔/𝑚3 and for gases we use air density at atmospheric pressure and 20oC 𝜌𝑎𝑖𝑟 = 1.205 𝑘𝑔/𝑚3 ) That is;
𝜌
𝑆𝐺 =
𝜌𝐻2𝑂
The density of a substance, in general, depends on temperature and pressure. The density of most gases is
proportional to pressure and inversely proportional to temperature. Liquids and solids densities, on the other
hand does not change much with pressure but changes with temperature.
The weight per unit volume of a substance is called specific weight or weight density and is expressed as; 𝛶𝑠 =
𝜌𝑔 where g is the gravitational acceleration.
24
University of Tabuk, Mechanical Engineering Department
Density and Specific Volume
25
University of Tabuk, Mechanical Engineering Department
Density and Specific Gravity
26
University of Tabuk, Mechanical Engineering Department
Fluid Properties
The vapor pressure Pv of a pure substance is defined as the pressure exerted by its vapor in phase
equilibrium with its liquid at a given temperature. Pv is a property of the pure substance and turns out to
be identical to the saturation pressure Psat of the liquid (Pv = Psat).
The rate of evaporation from open water bodies such as lakes is controlled by the difference between the
vapor pressure and the partial pressure of H2O in air.
27
University of Tabuk, Mechanical Engineering Department
Fluid Properties
For example, the vapor pressure of water vapor in air at 20°C is 2.34 kPa (N2 = 78
kpa, O2 = 21 kpa). Therefore, a bucket of water at 20°C left in a room with dry air (1
atm:101.325 kpa) will continue evaporating until one of two things happens:
1. The water evaporates away (there is not enough water to establish phase
equilibrium in the room), or
2. The evaporation stops when the partial pressure of the water vapor in the room
rises to 2.34 kPa at which point phase equilibrium is established
This fact is used for vacuum cooling!!!
28
University of Tabuk, Mechanical Engineering Department
Vapor pressure and Cavitation
• At lower pressures, liquids vaporizes at lower temperatures. Cavitation occurs as a result of
liquid vaporization in lower pressure (below Pv or Psat) regions of the flow (e.g. accelerated
flow due to sudden pipe contraction or other flow obstructions that reduces pressure), so
that vapor bubbles are generated withing the liquid flow.
• The collapse of those bubbles at regions of higher pressure generates highly destructive,
extremely high-pressure waves called cavitation; which can cause drop in performance
(causing vibrations) and even the erosion of pump impeller blades, pipes or other flow
passages. Cavitation is a limiting factor for maximum possible lift
29
University of Tabuk, Mechanical Engineering Department
Vapor pressure and Cavitation
Cavitation can be predicted using cavitation index (CI, σ, or KT), which is the ratio of potential to resist
vapor formation to the potential for causing the formation of the vapors.
For an upstream pressure, P1 and downstream pressure, P2,
the Pressure reduction due to flow acceleration;
𝟏
𝑷𝑺 = 𝑷𝟐 − 𝑷𝟏 = 𝝆(𝒗𝟐𝟐 − 𝒗𝟐𝟏 )
𝟐
Cavitation index, KT
𝑷𝟏 − 𝑷𝒗
𝑲𝑻 =
𝑷𝑺
The Table gives example for the range of KT values for valves
To avoid cavitation, gradual decrease in pressure is key (not sudden large drops).
30
University of Tabuk, Mechanical Engineering Department
Coefficient of Compressibility
We know from experience that the volume (or density) of a fluid changes with a change in
its temperature or pressure. Fluids usually expand as they are heated or depressurized and
contract as they are cooled or pressurized.
Two properties that relate volume changes to the changes in pressure and temperature are;
The coefficient of compressibility (also called the bulk modulus of compressibility or bulk
modulus of elasticity) for fluids, 𝜅, and the coefficient of volume expansion β.
𝜕𝑃 𝜕𝑃 ∆𝑃 ∆𝑃
𝜅 = −Ѵ =𝜌 ≅− ≅ (𝑝𝑎)
𝜕Ѵ 𝑇 𝜕𝜌 𝑇 ∆ѴΤѴ ∆𝜌Τ𝜌
A large value of 𝜅 indicates that a large change in pressure is needed to cause a small
fractional change in volume, and thus a fluid with a large 𝜅 is essentially incompressible.
For example, 𝜿 for water at normal atmospheric conditions is 21, 000 atm, which means
that the pressure of water must be raised to 210 atm to compress it 1 percent.
The acoustic waves (amplitude, ∆𝑃 = ρ𝑐𝑣 where 𝑐 = 𝜅Τ𝜌 is the pressure wave speed and 𝑣 is flow
velocity) that are produced strike the pipe surfaces, bends, and valves as they propagate and reflect along
the pipe, causing the pipe to vibrate and produce the familiar sound like the pipe is being hammered.
A second important parameter in the analysis of compressible fluid flow is the Mach
number Ma, named after the Austrian physicist Ernst Mach (1838–1916). It is the
ratio of the actual speed of the fluid (or an object in still fluid) to the speed of sound
in the same fluid at the same state:
𝑀𝑎 = 𝑉 Τ𝑐
Fluid flow regimes are often described in terms of the flow Mach number. The flow
is called sonic when Ma = 1, subsonic when Ma < 1, supersonic when Ma > 1,
hypersonic when Ma >> 1
33
University of Tabuk, Mechanical Engineering Department
Coefficient of Volume Expansion
The density of a fluid, in general, depends more strongly on temperature than it does on
pressure, and the variation of density with temperature is responsible for numerous
natural phenomena such as winds, currents in oceans, rise of plumes in chimneys, the
operation of hot-air balloons, and heat transfer by natural convection.
1 𝜕Ѵ 1 𝜕𝜌 ∆ѴΤѴ ∆𝜌Τ𝜌
β= =− ≅ ≅− (1Τ𝐾)
Ѵ 𝜕𝑇 𝑃
𝜌 𝜕𝑇 𝑃
∆𝑇 ∆𝑇
A large value expansion β means a large change in density with temperature, and the
product β∆T represents the fraction of volume change of a fluid that corresponds to a
temperature change of ∆T at constant pressure.
For ideal gases, β𝒊𝒅𝒆𝒂𝒍 𝒈𝒂𝒔 = 𝟏Τ𝑻 where T is the absolute temperature
34
University of Tabuk, Mechanical Engineering Department
PRESSURE
Absolute pressure is measured relative to absolute vacuum. Most pressure-measuring devices,
however, are calibrated to read zero in the atmosphere, and so they indicate the difference
between the absolute pressure and the local atmospheric pressure. This difference is called
the gage pressure. Pressures below atmospheric pressure are called vacuum pressures.
In thermodynamic relations and tables, absolute pressure is almost always used. In these
notes, the pressure P will denote absolute pressure unless specified otherwise
35
University of Tabuk, Mechanical Engineering Department
PRESSURE
What is the atmospheric pressure on a planet if the absolute pressure is 100 kPa and the gage
pressure is 10 kPa?
36
University of Tabuk, Mechanical Engineering Department
Viscosity
VISCOSITY refers to the internal resistance to flow. High internal resistance to flow means a higher
viscosity such as honey.
The force a flowing fluid exerts on a body in the flow direction is called the drag force,
and the magnitude of this force depends, in part, on viscosity.
37
University of Tabuk, Mechanical Engineering Department
Viscosity
Most common fluids like water, oxygen, nitrogen, and glycerin are
Newtonian
Viscosity of fluids can be determined with a sliding plate viscometer test. Consider two plates
of area A separated by a fluid with thickness, δ. The bottom plate is fixed, and the top plate is
kept in motion at a constant velocity, V, by a force, F.
40
University of Tabuk, Mechanical Engineering Department
Determination of Viscosity
A sliding-plate viscometer is used to measure the viscosity of a Newtonian fluid. A force of 25
N is required to keep the top plate moving at a constant velocity of 5 m/s. What is the
viscosity of the fluid?
41
University of Tabuk, Mechanical Engineering Department
Determination of Viscosity
The shear force acting on a Newtonian fluid layer in-between two parallel plates (or, by
Newton’s third law, the force acting on the plate) is:
𝑑𝑢
𝐹 = 𝜏𝐴 = 𝜇𝐴
𝑑𝑦
where again A is the contact area between the plate and the fluid. Then the force F required to
move the upper at a constant speed of V while the lower plate remains stationary is:
𝑉
𝐹 = 𝜏𝐴 = 𝜇𝐴
𝑙
This relation can alternately be used to calculate μ when the force F is measured. Therefore,
fluid confined in-between two concentric cylinders with the outer cylinder stationary while the
inner one is rotating can be used to measure the viscosity of fluid (a viscometer) in-between
the cylinders. For cylinders of length L and fluid layer with thickness ℓ,
Where torque; 𝑻 = 𝑭𝑹, Tangential velocity 𝑽 = 𝝎𝑹, 𝑨 = 𝟐𝝅𝑹𝑳 and 𝝎 = 𝟐𝝅𝑵Τ𝟔𝟎 N is the
rpm in revolution per minutes and 𝝎 is angular velocity in radians per second
42
University of Tabuk, Mechanical Engineering Department
Viscosity
43
University of Tabuk, Mechanical Engineering Department
Viscosity
44
University of Tabuk, Mechanical Engineering Department
Viscosity - Example
A 50-cm x 30-cm x 20-cm block weighing 150 N is to be moved at a constant velocity of 0.80
m/s on an inclined surface with a friction coefficient of 0.27. (a) Determine the force F that
needs to be applied in the horizontal direction. (b) If a 0.40-mm-thick oil film with a dynamic
viscosity of 0.012 Pa.s is applied between the block and inclined surface, determine the
percent reduction in the required force.
45
University of Tabuk, Mechanical Engineering Department
Viscosity - Example
46
University of Tabuk, Mechanical Engineering Department
Viscosity - Examples
47
University of Tabuk, Mechanical Engineering Department
Viscosity
48
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
It is often observed that a drop of blood forms a hump on a horizontal glass; a drop of
mercury forms a near-perfect sphere and can be rolled just like a steel ball over a
smooth surface; water droplets from rain or dew hang from branches or leaves of
trees.
In these and other observances, liquid droplets behave like small balloons filled with
the liquid, and the surface of the liquid acts like a stretched elastic membrane under
tension.
The pulling force that causes this tension acts parallel to the surface and is due to the
attractive forces between the molecules of the liquid. The magnitude of this force per
unit length is called surface tension or coefficient of surface tension σs and is usually
expressed in the unit N/m (or lbf/ft in English units).
49
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
𝑭
Surface tension 𝝈𝒔 =
𝒃
Water droplets, soap bubbles, needle floating on water, and filled cup are
examples of the effect of surface tension 50
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
Typical range of values for σs : 0.001 – 0.1 N/m. For water @ 20 °C, σs = 0.073 N/m.
Surface tension:
• It decreases with temperature.
• Strongly affected by impurities. e.g. Adding soap to water
51
University of Tabuk, Mechanical Engineering Department
Surface Tension Measurement
Wire frame Apparatus: Surface tension is determined by measuring the force
necessary to keep the sliding side stationary against the surface tension pull of the
film. Since the film has two surfaces (i.e two surface tensions), the surface tension is:
𝑭
𝝈=
𝟐𝑳
52
University of Tabuk, Mechanical Engineering Department
Surface Tension
The film width in a surface tension experiment is 10 cm. If mercury is the fluid (surface
tension = 0.52 N/m), what is the maximum force that can be applied without breaking the
membrane? Neglect gravitational force.
53
University of Tabuk, Mechanical Engineering Department
Surface Tension Measurement
Du Nouy Ring Apparatus: Alternatively, surface tension can also be measured by
measuring the force required to a Du Nuoy wire ring out of the liquid as shown in the
figure. Because the ring's inner and outer sides are in contact with the liquid, the
wetted perimeter is twice the circumference, the surface tension is:
𝑭
𝝈=
𝟒𝝅𝒓
54
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
Another interesting consequence of surface tension is the capillary effect, which is the
rise or fall of a liquid in a small-diameter tube inserted into the liquid. Such narrow
tubes or confined flow channels are called capillaries. The curved free surface
of a liquid in a capillary tube is called the meniscus.
Liquid is said to wet the surface when ϕ < 90° and not to wet the surface when ϕ > 90°
55
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
56
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
57
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
To derive an equation for capillary rise,
define a system comprised of the water
inside the capillary tube. Then, draw a free
body diagram (FBD). As shown, the pull of
surface tension lifts the column of water.
58
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
59
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
60
University of Tabuk, Mechanical Engineering Department
Surface Tension and Capillary Action
Notes:
• The capillary rise is inversely proportional to the radius of the tube. Therefore, the thinner the tube is,
the greater the rise (or fall) of the liquid in the tube.
• In practice, the capillary effect for water is negligible in tubes whose diameter is greater than 1cm.
• When pressure measurements are made using manometers and barometers, it is important to use
sufficiently large tubes to minimize the capillary effect.
• The capillary rise is also inversely proportional to the density of the liquid, as expected. Therefore, in
general, lighter liquids experience greater capillary rises.
• In atmospheric air, the contact angle of water (and most other organic liquids) with glass is nearly zero.
The contact angle is 130° for mercury–glass and 26° for kerosene–glass in air. The contact angle, in
general, is different in different environments.
• Keep in mind that the equation for column rise (h) is derived for constant-diameter tubes and should
not be used for tubes of variable cross section
61
University of Tabuk, Mechanical Engineering Department
Summary – Fluid Properties
Fluid Properties
Density or
Specific Vapor Modulus of Volume Surface
Specific Viscosity
Weight Pressure elasticity Expansivity Tension
Gravity
62
University of Tabuk, Mechanical Engineering Department
Summary - Forces
63
University of Tabuk, Mechanical Engineering Department
Exercises
2-81- A thin 30-cm x 30-cm flat plate is pulled at 3 m/s horizontally through a 3.6-mm-thick oil
layer sandwiched between two plates, one stationary and the other moving at a constant velocity
of 0.3 m/s, as shown in Fig. P2–81. The dynamic viscosity of the oil is 0.027 Pa.s. Assuming the
velocity in each oil layer to vary linearly, (a) plot the velocity profile and find the location where
the oil velocity is zero and (b) determine the force that needs to be applied on the plate to
maintain this motion.
64
University of Tabuk, Mechanical Engineering Department
Exercises
65
University of Tabuk, Mechanical Engineering Department
Exercises
2-87- For flow over a plate, the variation of velocity with vertical distance y from the plate is given
as u(y) = ay - by2 where a and b are constants. Obtain a relation for the wall shear stress in terms of
a, b, and µ.
66
University of Tabuk, Mechanical Engineering Department
Exercises
2-94-A thin plate moves between two parallel, horizontal, stationary flat surfaces at a constant
velocity of 5 m/s. The two stationary surfaces are spaced 4 cm apart, and the medium between
them is filled with oil whose viscosity is 0.9 N.s/m2. The part of the plate immersed in oil at any
given time is 2-m long and 0.5-m wide. If the plate moves through the mid-plane between the
surfaces, determine the force required to maintain this motion. What would your response be if the
plate was 1 cm from the bottom surface (h2) and 3 cm from the top surface (h1)?
67
University of Tabuk, Mechanical Engineering Department
Exercises
68
University of Tabuk, Mechanical Engineering Department
Exercises
2–103- A 1.2-mm-diameter tube is inserted into an unknown liquid whose density is 960 kg/m3,
and it is observed that the liquid rises 5 mm in the tube, making a contact angle of 15°.
Determine the surface tension of the liquid.
69
University of Tabuk, Mechanical Engineering Department
Exercises
70
University of Tabuk, Mechanical Engineering Department