Module 3 Lesson 1-3
Module 3 Lesson 1-3
Module III
84
MODULE III
HEAT TRANSFER BY CONVECTION
INTRODUCTION
OBJECTIVES
There are three lessons in the module. Read each lesson carefully then answer
the exercises/activities to find out how much you have benefited from it. Work on
these exercises carefully and submit your output to your instructor.
In case you encounter difficulty, discuss this with your instructor during the
face-to-face meeting.
Lesson 1
Introduction to Convection
Introduction to Convection
It is well known that a hot plate of metal will cool faster when placed in front of a
fan than when exposed to still air.
For example, we know that the velocity at which the air blows over the hot plate
obviously influences the heat transfer rate. But does it influence the cooling in a
linear way? i.e. if the velocity is doubled, will the heat transfer rate have doubled?
Thus, we might compute the heat transfer using Fourier’s equation of conduction
i.e. 𝑞 = −k𝐴 d𝑡/d𝑥 with the thermal conductivity of fluid and the fluid temperature
gradient at wall.
Why then, if the heat flows by conduction in this layer, do we speak of “Convection”
heat transfer and need to consider the velocity of the fluid?
The answer is that the temperature gradient is dependent on the rate at which the
fluid carries the heat away; a high velocity produces a large temperature gradient,
and so on.
It must be remembered that the physical mechanism of heat transfer at the wall is
a conduction process.
q = hA ( ts - t ) → Equation 3.1
Where:
𝑄 = Convective heat flow rate
𝐴 = Surface area exposed to heat transfer
𝑡𝑠 = Surface temperature of solid and
𝑡∞ = Temperature of the fluid (Stagnant or Undisturbed)
ℎ = The convective heat transfer coefficient or the film
coefficient or the surface conductance
W W Cal Btu
; ; ;
m −K
2
m −C
2
m − hr − C hr − ft 2 − o F
2 o
The viscosity influences the velocity profile and correspondingly the energy
transfer rate in the region near the wall.
With respect to the cause of fluid flow, two types of convection are
distinguished:
1. Free Convection or Natural Convection and
2. Forced Convection.
A density difference is created between this adjacent layer and the still fluid
surrounding it.
The density difference introduces a buoyant force causing flow of fluid near the
surface.
Thus, “Free or Natural convection is the process of heat transfer which occurs
due to movement of the fluid particles by density changes associated with
temperature differential in a fluid.”
This mode of heat transfer occurs very commonly, some of the examples are:
1. House heating system
2. The cooling of transmission lines, electric transformers and rectifiers.
Forced Convection
Examples of forced convection are: cooling of I.C. Engines, Air conditioner, Heat
exchangers, etc.
Solved Problems:
Problem 16:
A motor cycle cylinder consists of ten fins, each 150 mm outside diameter and
75 mm inside diameter. The average fin temperature is 500oC and the
surrounding air is at 20oC temperature. Make calculations for the rate of heat
dissipation from the cylinder fins by convection when (i) motor cycle is
stationery and convection coefficient h:6 W/m2K (ii) motor cycle is moving at
60 km/h? and h = 75 W/m2K.
Solution:
Since both sides of each fin are exposed to the surroundings air, forced airflows over
a convective heat exchanger in a room heater, resulting in a convective heat
A = 10 2 ( 0.152 − 0.0752 ) = 0.265m2
4
The heat flow from the surface of the fluid is given by Newton’s Law of cooling:
Q = hA ( ts − t )
Q = hA ( ts − t )
(
Q = 6W
m2 − K ) ( 0.265m ) (500 − 20) C
2 o
Q = 763.2W
Q = hA ( ts − t )
(
Q = 75W
m2 − K ) (0.265m ) (500 − 20) C
2 o
Q = 9, 6540W
Problem 17:
Forced air flows over a convective heat exchanger in a room heater, resulting
in a convective heat transfer coefficient 1.136 kW/m2-K. The surface
temperature of heat exchanger may be considered constant at 65oC, and the air
is at 20oC. Determine the heat exchanger surface area required for 8.8 kW of
heating.
Solution:
The convective heat flow from a solid surface to the surrounding fluid is given by
Newton’s Law of Cooling
Q = hA ( ts − t )
kW
8.8kW = 1.136 2 ( A )( 65 − 20 ) C
o
m −K
A = 0.172m 2
Lesson 2
Natural Convection
Natural Convection
Consider a hot object exposed to cold air. The temperature of the outside of
the object will drop (as a result of heat transfer with cold air), and the temperature
of adjacent air to the object will rise. Consequently, the object is surrounded with a
thin layer of warmer air and heat will be transferred from this layer to the outer
layers of air.
The temperature of the air adjacent to the hot object is higher, thus its density is
lower. As a result, the heated air rises. This movement is called the natural
convection current. Note that in the absence of this movement, heat transfer would
be by conduction only and its rate would be much lower.
In a gravitational field, there is a net force that pushes a light fluid placed in a heavier
fluid upwards. This force is called the buoyancy force.
The magnitude of the buoyancy force is the weight of the fluid displaced by the body.
where Vbody is the volume of the portion of the body immersed in the fluid.
Note that the net force is proportional to the difference in the densities of the fluid
and the body. This is known as Archimedes’ principle.
We all encounter the feeling of “weight loss” in water which is caused by the
buoyancy force. Other examples are hot balloon rising, and the chimney effect.
Note that the buoyancy force needs the gravity field, thus in space (where no gravity
exists) the buoyancy effects does not exist.
where T is the absolute temperature. Note that the parameter βΔT represents the
fraction of volume change of a fluid that corresponds to a temperature change ΔT at
constant pressure.
Since the buoyancy force is proportional to the density difference, the larger the
temperature difference between the fluid and the body, the larger the buoyancy
force will be.
Whenever two bodies in contact move relative to each other, a friction force develops
at the contact surface in the direction opposite to that of the motion. Under steady
conditions, the air flow rate driven by buoyancy is established by balancing the
buoyancy force with the frictional force.
Grashof Number
It is also expressed as
where
g = gravitational acceleration, m/s2
β = coefficient of volume expansion, 1/K
δ = characteristic length of the geometry, m
ν = kinematics viscosity of the fluid, m2/s
It indicates the relative importance of the inertial and viscous effects in a fluid
motion.
At low Reynolds number, the viscous effect dominates and the fluid motion is
laminar.
For highly viscous oils, Pr is quite large (100 to 10000) and that indicates rapid
diffusion of momentum by viscous action compared to the diffusion of energy.
For gases, Pr is about 1, which indicates that both momentum and heat dissipate
through the field at about the same rate.
The liquid metals (liquid sodium or liquid potassium) have Pr = 0.003 to 0.01 and
that indicates more rapid diffusion of energy compared to the momentum diffusion
rate.
The Prandtl number is connecting link between the velocity field and the
temperature field, and its value strongly influences relative growth of velocity and
thermal boundary layers.
The role played by Reynolds number in forced convection is played by the Grashoff
number in natural convection.
The critical Grashof number is observed to be about 109 for vertical plates. Thus,
the flow regime on a vertical plate becomes turbulent at Grashof number greater
than 109. The heat transfer rate in natural convection is expressed by Newton’s law
of cooling as:
Q = h A (Ts ‐ T∞)
Figure 3.4 Velocity and temperature profile for natural convection flow over a
hot vertical plate. Grcritical = 109
Natural Convection over Surfaces
The velocity and temperature distribution for natural convection over a hot vertical
plate are shown in Figure 3.4.
Note that the velocity at the edge of the boundary layer becomes zero. It is
expected since the fluid beyond the boundary layer is stationary.
The shape of the velocity and temperature profiles, in the cold plate case, remains
the same but their direction is reversed.
The complexities of the fluid flow make it very difficult to obtain simple analytical
relations for natural convection. Thus, most of the relationships in natural convection
are based on experimental correlations.
The Rayleigh number is defined as the product of the Grashof and Prandtl numbers:
where the constants C and n depend on the geometry of the surface and the flow.
Table 14‐1 in Cengel book lists these constants for a variety of geometries.
These relationships are for isothermal surfaces, but could be used approximately for
the case of non‐isothermal surfaces by assuming surface temperature to be constant
at some average value.
The characteristics length is A/p where the surface area is A, and perimeter is p.
Solved Problems:
Solution:
Finned surfaces of various shapes (heat sinks) are used in microelectronics cooling.
One of most crucial parameters in designing heat sinks is the fin spacing, S. Closely
packed fins will have greater surface area for heat transfer, but a smaller heat
transfer coefficient (due to extra resistance of additional fins). A heat sink with
widely spaced fins will have a higher heat transfer coefficient but smaller surface
area. Thus, an optimum spacing exists that maximizes the natural convection from
the heat sink.
Consider a heat sink with base dimension W (width) and L (length) in which the fins
are assumed to be isothermal and the fin thickness t is small relative to fin spacing
S. The optimum fin spacing for a vertical heat sink is given by Rohsenow and Bar‐
Cohen as
where L is the characteristic length in Ra number. All the fluid properties are
determined at the film temperature. The heat transfer coefficient for the optimum
spacing can be found from
Assumptions:
The fin thickness t is much smaller than the fin spacing S.
Solution:
The properties of air are evaluated at the film temperature:
The number of fins and the heat transfer coefficient for the optimum fin
spacing case are
Problem 20.
A hot plate of 400mm x 400mm at 100°C is exposed to air at 20°C. Calculate heat
loss from both the surfaces of the plate if (a) the plate is kept vertical (b) plate
is kept horizontal. Air properties at mean temperature are ρ = 1.06 kg/m3, k =
0.028 W/m-k, Cp= 1.008 KJ/kg-k, and ν = 18.97 x 10-6 m2/s.
Use following correlations:
𝑁𝑢 = 0.125(𝐺𝑟. 𝑃𝑟)0.33 𝑓𝑜𝑟 𝑣𝑒𝑟𝑡𝑖𝑐𝑎𝑙 𝑝𝑙𝑎𝑡𝑒
𝑁𝑢 = 0.72(𝐺𝑟. 𝑃𝑟)0.25 𝑓𝑜𝑟 𝑢𝑝𝑝𝑒𝑟 𝑠𝑢𝑟𝑓𝑎𝑐𝑒
𝑁𝑢 = 0.35(𝐺𝑟. 𝑃𝑟)0.25 𝑓𝑜𝑟 𝑙𝑜𝑤𝑒𝑟 𝑠𝑢𝑟𝑓𝑎𝑐𝑒
Lesson 3
Forced Convection
Convection heat transfer is complicated since it involves fluid motion as well as heat
conduction. The fluid motion enhances heat transfer (the higher the velocity the
higher the heat transfer rate).
q = hA (Ts − T ) → W
m2
Q = h (Ts − T ) → W
The convective heat transfer coefficient h strongly depends on the fluid properties
and roughness of the solid surface, and the type of the fluid flow (laminar or
turbulent).
It is assumed that the velocity of the fluid is zero at the wall, this assumption is called
noslip condition. As a result, the heat transfer from the solid surface to the fluid
layer adjacent to the surface is by pure conduction, since the fluid is motionless.
Thus,
The convection heat transfer coefficient, in general, varies along the flow direction.
The mean or average convection heat transfer coefficient for a surface is determined
by (properly) averaging the local heat transfer coefficient over the entire surface.
Consider the flow of a fluid over a flat plate, the velocity and the temperature of the
fluid approaching the plate is uniform at U∞ and T∞. The fluid can be considered as
adjacent layers on top of each other.
Assuming no‐slip condition at the wall, the velocity of the fluid layer at the wall is
zero. The motionless layer slows down the particles of the neighboring fluid layers as
a result of friction between the two adjacent layers. The presence of the plate is felt
up to some distance from the plate beyond which the fluid velocity U∞ remains
unchanged. This region is called velocity boundary layer.
Boundary layer region is the region where the viscous effects and the velocity changes
are significant and the inviscid region is the region in which the frictional effects are
negligible and the velocity remains essentially constant.
The friction between two adjacent layers between two layers acts similar to a drag
force (friction force). The drag force per unit area is called the shear stress:
➢ The flow in boundary layer starts as smooth and streamlined which is called
laminar flow.
➢ At some distance from the leading edge, the flow turns chaotic, which is called
turbulent and it is characterized by velocity fluctuations and highly
disordered motion.
➢ The transition from laminar to turbulent flow occurs over some region which
is called transition region.
➢ The intense mixing of the fluid in turbulent flow enhances heat and
momentum transfer between fluid particles, which in turn increases the
friction force and the convection heat
transfer coefficient.
Non‐dimensional Groups
where δ is the characteristic length, i.e. D for the tube and L for the flat plate.
Nusselt number represents the enhancement of heat transfer through a fluid as a
result of convection relative to conduction across the same fluid layer.
At large Re numbers, the inertia forces, which are proportional to the density and
the velocity of the fluid, are large relative to the viscous forces; thus, the viscous
forces cannot prevent the random and rapid fluctuations of the fluid (turbulent
regime).
The Reynolds number at which the flow becomes turbulent is called the critical
Reynolds number. For flat plate the critical Re is experimentally determined to be
approximately Re critical = 5 x105.
➢ The relative thickness of the velocity and the thermal boundary layers is
described by the Prandtl number.
➢ For low Prandtl number fluids, i.e. liquid metals, heat diffuses much faster
than momentum flow (remember Pr = ν/α<<1) and the velocity boundary layer
is fully contained within the thermal boundary layer. On the other hand, for
high Prandtl number fluids, i.e. oils, heat diffuses much slower than the
momentum and the thermal boundary layer is contained within the velocity
boundary layer.
The friction and heat transfer coefficient for a flat plate can be determined
by solving the conservation of mass, momentum, and energy equations (either
approximately or numerically). They can also be measured experimentally. It is found
that the Nusselt number can be expressed as:
where C, m, and n are constants and L is the length of the flat plate. The properties
of the fluid are usually evaluated at the film temperature defined as:
Laminar Flow
The local friction coefficient and the Nusselt number at the location x for
laminar flow over a flat plate are
where x is the distant from the leading edge of the plate and Rex = ρV∞x/μ.
The averaged friction coefficient and the Nusselt number over the entire isothermal
plate for laminar regime are:
Taking the critical Reynolds number to be 5 x105, the length of the plate xcr over
which the flow is laminar can be determined from
Turbulent Flow
The local friction coefficient and the Nusselt number at location x for
turbulent flow over a flat isothermal plate are:
The averaged friction coefficient and Nusselt number over the isothermal
plate in turbulent region are:
If the plate is sufficiently long for the flow to become turbulent (and not long enough
to disregard the laminar flow region), we should use the average values for friction
coefficient and the Nusselt number.
where the critical Reynolds number is assumed to be 5x105. After performing the
integrals and simplifications, one obtains:
The above relationships have been obtained for the case of isothermal surfaces, but
could also be used approximately for the case of non‐isothermal surfaces. In such
cases assume the surface temperature be constant at some average value.
For isoflux (uniform heat flux) plates, the local Nusselt number for laminar and
turbulent flow can be found from:
Note the isoflux relationships give values that are 36% higher for laminar and 4% for
turbulent flows relative to isothermal plate case.
Problem 21:
Engine oil at 60°C flows over a 5 m long flat plate whose temperature is 20°C
with a velocity of 2 m/s. Determine the total drag force and the rate of heat
transfer per unit width of the entire plate.
Solution:
We assume the critical Reynolds number is 5x105. The properties of the oil at the
film temperature are:
which is less than the critical Re. Thus we have laminar flow. The friction coefficient
and the drag force can be found from:
The critical Re for the flow across spheres or tubes is 2x105. The approaching
fluid to the cylinder (a sphere) will branch out and encircle the body, forming a
boundary layer.
Figure 3.9 Typical flow patterns over sphere and streamlined body and drag forces.
At low Re (Re < 4) numbers the fluid completely wraps around the body. At higher Re
numbers, the fluid is too fast to remain attached to the surface as it approaches the
top of the cylinder. Thus, the boundary layer detaches from the surface, forming a
wake behind the body. This point is called the separation point.
To reduce the drag coefficient, streamlined bodies are more suitable, e.g. airplanes
are built to resemble birds and submarine to resemble fish, Figure 3.9.
In flow past cylinder or spheres, flow separation occurs around 80° for laminar flow
and 140° for turbulent flow.
The drag force acting on a body is caused by two effects: the friction drag (due to
the shear stress at the surface) and the pressure drag which is due to pressure
differential between the front and rear side of the body.
As a result of transition to turbulent flow, which moves the separation point further
to the rear of the body, a large reduction in the drag coefficient occurs. As a result,
the surface of golf balls is intentionally roughened to induce turbulent at a lower Re
number, see Figure 3.10.
The average heat transfer coefficient for cross‐flow over a cylinder can be found from
the correlation presented by Churchill and Bernstein:
Tf = (Ts + T∞)/2
which is valid for 3.5 < Re < 80,000 and 0.7 < Pr < 380. The fluid properties are
evaluated at the free‐stream temperature T∞, except for μs which is evaluated at
surface temperature.
The average Nusselt number for flow across circular and non‐circular cylinders can
be found from Table 10‐3 Cengel book.
Problem 22:
The decorative plastic film on a copper sphere of 10‐mm diameter is cured in
an oven at 75°C. Upon removal from the oven, the sphere is subjected to an
air stream at 1 atm and23°C having a velocity of 10 m/s, estimate how long it
will take to cool the sphere to 35°C.
Assumptions:
1. Negligible thermal resistance and capacitance for the plastic layer.
2. Spatially isothermal sphere.
3. Negligible Radiation.
The time required to complete the cooling process may be obtained from the
results for a lumped capacitance.
Whitaker relationship can be used to find h for the flow over sphere:
Hence,
MODULE SUMMARY
The rate of heat transfer by convection from the surface of a conducting body
to the surrounding fluid is given by Q = hA ( t s − t ) where ts is the temperature of the
surface, t∞ is the temperature of fluid and A is the surface area of the body. The
parameter h is the convective heat transfer coefficient. The convective coefficient
is a complicated function of (i) geometry of the system, (ii) thermal properties of the
fluid, and (iii) characteristics of fluid flow.
LEARNING TASK 3
Practice Problems
1. The forced convective heat transfer coefficient for a hot fluid flowing over
a cool surface has been estimated to be 816 kJ/m2-oC. The fluid temperature
upstream of the cool surface is 120oC, and the surface is held at 10oC.
Determine the heat transfer, per unit surface area, from the fluid to the
surface.
References:
1. Principles of Heat Transfer by Frank Kreith, Cengage Learning. 2018
2. Heat and Mass Transfer, 3G eLearning. 2014
3. Heat and Mass Transfer by D. S. Kumar, S K Kataria and Sons Publications.2014
4. Foundation of Heat Transfer, Incropera, Frank. John Wiley. 2013
5. Basics of Heat and Mass Transfer. Kumar, D.S. S.K. Kataria & Sons.2011