Lecture 02 Conduction
Lecture 02 Conduction
Conduction
by
Dr Raj Kumar Arya
[PhD(IITB), M.Tech.(IITD), B.Tech.(HBTIK)]
Associate Professor
Department of Chemical Engineering
Dr. B. R. Ambedkar National Institute of Technology, Jalandhar
Webpage: https://www.nitj.ac.in/index.php/nitj_cinfo/Faculty/180
https://www.youtube.com/channel/UCbbCr7YE-vDAtpn3sivdsoQ/
www.drrajkumararya.com
aryark@nitj.ac.in
Example
A 10-cm diameter copper ball is to be heated from 100°C to an
average temperature of 150°C in 30 minutes. Taking the average
density and specific heat of copper in this temperature range to be
8950 kg/m3 and Cp 0.395 kJ/kg · °C, respectively, determine (a)
the total amount of heat transfer to the copper ball, (b) the average
rate of heat transfer to the ball, and (c) the average heat flux.
3:29 AM 2
SOLUTION The copper ball is to be heated from 100°C to 150°C. The total heat
transfer, the average rate of heat transfer, and the average heat flux are to
be determined.
Analysis (a) The amount of heat transferred to the copper ball is simply the
change in its internal energy, and is determined from
3:29 AM 3
where
Substituting,
(b) The rate of heat transfer normally changes during a process with time.
However, we can determine the average rate of heat transfer by dividing the total
amount of heat transfer by the time interval. Therefore,
3:29 AM 4
(c) Heat flux is defined as the heat transfer per unit time per unit area, or the
rate of heat transfer per unit area. Therefore, the average heat flux in this case
is
Discussion Note that heat flux may vary with location on a surface. The value
calculated above is the average heat flux over the entire surface of the ball.
3:29 AM 5
Example
1.2 kg of liquid water initially at 15°C is to be heated to 95°C in a teapot equipped
with a 1200-W electric heating element inside. The teapot is 0.5 kg and has an
average specific heat of 0.7 kJ/kg · °C. Taking the specific heat of water to be 4.18
kJ/kg · °C and disregarding any heat loss from the teapot, determine how long it
will take for the water to be heated.
3:29 AM 6
SOLUTION
Liquid water is to be heated in an electric teapot. The heating time
is to be determined.
Assumptions
1. Heat loss from the teapot is negligible.
2. Constant properties can be used for both the teapot and the water.
Properties The average specific heats are given to be 0.7 kJ/kg · °C for the
teapot and 4.18 kJ/kg · °C for water.
Analysis We take the teapot and the water in it as the system, which is a closed
system (fixed mass). The energy balance in this case can be expressed as
3:29 AM 7
Then the amount of energy needed to raise the temperature of water and the
teapot from 15°C to 95°C is
The 1200-W electric heating unit will supply energy at a rate of 1.2 kW or 1.2
kJ per second. Therefore, the time needed for this heater to supply 429.3 kJ
of heat is determined from
3:29 AM 8
Discussion
In reality, it will take more than 6 minutes to accomplish this heating process
since some heat loss is inevitable during heating.
3:29 AM 9
HEAT TRANSFER MODES
Conduction
- needs matter
- molecular phenomenon (diffusion process)
- without bulk motion of matter
Convection
- heat carried away by bulk motion of fluid
- needs fluid matter
Radiation
- does not needs matter
- transmission of energy by electromagnetic waves
3:29 AM 10
CONDUCTION
Conduction: The transfer of energy from the more energetic
particles of a substance to the adjacent less energetic ones as a
result of interactions between the particles.
In gases and liquids, conduction is due to the collisions and
diffusion of the molecules during their random motion.
In solids, it is due to the combination of vibrations of the molecules
in a lattice and the energy transport by free electrons.
3:29 AM 11
The rate of heat conduction through a plane layer is proportional to the
temperature difference across the layer and the heat transfer area, but is
inversely proportional to the thickness of the layer.
Heat conduction through a large plane wall of
thickness x and area A.
Hence,
3:29 AM 12
When x→ 0
3:29 AM 13
The rate of heat conduction through a
solid is directly proportional to its
thermal conductivity.
3:29 AM 16
Statement of Fourier’s Law
The (mod of) heat flux, q (the flow of heat per unit area and per unit time), at a
point in a medium is directly proportional to the temperature gradient at the
point.
dT
Qx = − kA
&
dx
3:29 AM 18
A rate equation that allows determination of the conduction
heat flux from knowledge of the temperature distribution in a
medium
→ →
q = −k T
3:29 AM 19
Fourier’s Law serves to define the thermal conductivity of the
medium
→ →
k − q / T
3:29 AM 20
Mathematical Description
• Temperature is a scalar quantity.
• Gradient of temperature (▼T) will be a vector.
• Heat flux is defined with direction and Magnitude : A
Vector.
• Mathematically it is possible to represent heat flux as:
r
q = qxiˆ + q y ˆj + qz kˆ
Using the principles of vector calculus:
r
q = − k T
r T ˆ T ˆ T ˆ
q = −k i + j+ k
x y z
3:29 AM 21
Thermal Conductivity
Thermal conductivity: The rate of heat transfer through a unit
thickness of the material per unit area per unit temperature
difference.
3:29 AM 22
A simple experimental setup to determine the thermal conductivity of a
material.
23
The range of thermal conductivity of various
materials at room temperature.
24
The thermal conductivities of gases such as air
vary by a factor of 104 from those of pure
metals such as copper.
Pure crystals and metals have the highest
thermal conductivities, and gases and insulating
materials the lowest.
3:29 AM 26
Mathematical Fantasy to Physical Existence
• Will k be same in all directions?
• Why k cannot be different each direction?
• Why k cannot be a vector?
k = k xiˆ + k y ˆj + k z kˆ
Will mathematics approve this ?
3:29 AM 27
Most General form of Fourier Law of Conduction
dT
qx = − k
dx
r
q = − k T
k = k xiˆ + k y ˆj + k z kˆ
3:29 AM 28
Physically -mathematically Feasible Model
• Taking both physics and mathematics into consideration,
the most feasible model for Fourier’s Law of conduction is:
r %
q = − k%
.T
Thermal conductivity of a general material is a tensor.
T
x
qx k xx k xy k xz
T
q y = − k yx k yy k yz
y
qz k zx k zy k zz T
3:29 AM
z 29
Surprising Inventions !!!
T T T
qx = − k xx + k xy + k xz
x y z
T T T
q y = − k yx + k yy + k yz
x y z
T T T
qz = − k zx + k zy + k zz
x y z
3:29 AM 30
Thermal Diffusivity
3:29 AM 31
A material that has a high thermal
conductivity or a low heat capacity
will obviously have a large thermal
diffusivity.
3:29 AM 33
Solution
Assumptions
1. Steady-state one-dimensional conduction
2. Thermal conductivity is constant for the temperature range of interest
3. The heat loss through the edge side surface is insignificant
4. The layers are in perfect thermal contact
Given,
3:29 AM 34
Example
The roof of an electrically heated home is 6 m long, 8 m wide, and 0.25 m thick,
and is made of a flat layer of concrete whose thermal conductivity is k 0.8 W/m ·
°C . The temperatures of the inner and the outer surfaces of the roof one night
are measured to be 15°C and 4°C, respectively, for a period of 10 hours.
Determine
(a) the rate of heat loss through the roof that
night and
(b) the cost of that heat loss to the home owner if the cost of electricity
is $0.08/kWh.
3:29 AM 35
SOLUTION
The inner and outer surfaces of the flat concrete roof of an electrically heated
home are maintained at specified temperatures during a night. The heat loss
through the roof and its cost that night are to be determined.
Assumptions
1. Steady operating conditions exist during the entire night since the surface
temperatures of the roof remain constant at the specified values.
2. Constant properties can be used for the roof.
Properties
The thermal conductivity of the roof is given to be k = 0.8 W/m · °C.
3:29 AM 36
Analysis
(a) Noting that heat transfer through the roof is by conduction and the area of
the roof is A = 6 m × 8 m = 48 m2, the steady rate of heat transfer through the
roof is determined to be
(b) The amount of heat lost through the roof during a 10-hour period and its
cost are determined from
Discussion The cost to the home owner of the heat loss through the roof that night was
$1.35. The total heating bill of the house will be much larger since the heat losses through the
walls are not considered in these calculations.
3:29 AM 37
Coordinate systems
3:29 AM 38
Steady vs Transient Heat Transfer
3:29 AM 39
Lumped systems
Most heat transfer problems encountered in practice are
transient in nature, but they are usually analyzed under some
presumed steady conditions since steady processes are easier
to analyze, and they provide the answers to our questions.
3:29 AM 40
Multidimensional Heat Transfer
Heat transfer problems are also classified as being
One-dimensional, 𝑻 = 𝒇(𝒙)
Two-dimensional, 𝑻 = 𝒇(𝒙, 𝒚)
or
Three-dimensional, 𝑻 = 𝒇(𝒙, 𝒚, 𝒛)
3:29 AM 41
The temperature gradient dT/dx is simply the slope of the
temperature curve on a T-x diagram.
3:29 AM 42
Heat Generation
A medium through which heat is conducted may involve the
conversion of electrical, nuclear, or chemical energy into heat (or
thermal) energy. In heat conduction analysis, such conversion
processes are characterized as heat generation.
Oven
3:29 AM 43
➢Heat generation is a volumetric phenomenon.
➢The rate of heat generation in a medium is usually specified
per unit volume and is denoted by 𝒈,ሶ whose unit is 𝐖/𝐦𝟑 or
Btu/h.ft3.
➢The rate of heat generation in a medium may vary with time
as well as position within the medium.
➢The total rate of heat generation, 𝑮ሶ in a medium of volume V
can be determined from
3:29 AM 44
EXAMPLE
The resistance wire of a 1200-W hair dryer is 80 cm long and has a
diameter of D = 0.3 cm. Determine the rate of heat generation in
the wire per unit volume, in W/cm3, and the heat flux on the outer
surface of the wire as a result of this heat generation.
3:29 AM 45
SOLUTION
The power consumed by the resistance wire of a hair dryer is given.
The heat generation and the heat flux are to be determined.
Assumptions
Heat is generated uniformly in the resistance wire.
Analysis
A 1200-W hair dryer will convert electrical energy into heat in the wire at a
rate of 1200 W. Therefore, the rate of heat generation in a resistance wire
is equal to the power consumption of a resistance heater. Then the rate of
heat generation in the wire per unit volume is determined by dividing the
total rate of heat generation by the volume of the wire,
3:29 AM 46
Similarly, heat flux on the outer surface of the wire as a
result of this heat generation is determined by dividing the
total rate of heat generation by the surface area of the wire,
3:29 AM 47
ONE-DIMENSIONAL HEAT
CONDUCTION
Heat Conduction Equation in a Large Plane Wall
3:29 AM 49
An energy balance on this thin element during a small time
interval Δt can be expressed as
3:29 AM 50
But the change in the energy content of the element and the
rate of heat generation within the element can be expressed as
Substituting
3:29 AM 51
since, from the definition of the derivative and Fourier’s law
of heat conduction and taking the limit as Δx → 0
Hence,
3:29 AM
It is for variable thermal conductivity 52
for constant thermal conductivity
𝒌
We know that thermal diffusivity, 𝜶 =
𝝆𝑪
3:29 AM 53
For constant k
Steady-state:
𝝏𝑻
=𝟎
𝝏𝒕
For constant k
Transient and no
heat generation:
𝒈ሶ = 𝟎
For constant k
Steady-state and no
heat generation:
𝝏𝑻
𝒈ሶ = 𝟎; =𝟎
𝝏𝒕
3:29 AM 54
Heat Conduction Equation in a Long Cylinder
3:29 AM 55
consider a thin cylindrical shell element of thickness Δr in a long
this case, and thus it varies with location. An energy balance on this
expressed as
3:29 AM 56
……………..(1)
The change in the energy content of the element and the rate
of heat generation within the element can be expressed as
……………..(2)
3:29 AM 57
Substituting values in eq. (1) from eq. (2)
……………..(3)
……………..(4)
3:29 AM 58
Taking the limit as Δr → 0 and Δ t → 0 and
applying Fourier's law yields
……………..(5)
Hence,
……………..(6)
Noting that the heat transfer area in this case is A =2πrL, the
one-dimensional transient heat conduction equation in a
cylinder becomes
3:29 AM 59
……………..(7)
……………..(8)
3:29 AM 60
For constant k
Steady-state:
𝝏𝑻
=𝟎
𝝏𝒕
For constant k
Transient and no
heat generation:
𝒈ሶ = 𝟎
For constant k
Steady-state and no
heat generation:
𝝏𝑻
𝒈ሶ = 𝟎; =𝟎
𝝏𝒕
3:29 AM 61
Heat Conduction Equation in a Sphere
3:29 AM 62
Considering a sphere with constant density ρ, specific heat C,
and outer radius R. The area of the sphere normal to the
direction of heat transfer at any location is A= 4πr2, where r is
the value of the radius at that location. Note that the heat
transfer area A depends on r in this case also, and thus it
varies with location.
By considering a thin spherical shell element of thickness Δr
The one-dimensional transient heat conduction equation for
a sphere is determined to be
3:29 AM 63
For variable thermal conductivity, k
3:29 AM 64
𝝏𝑻
For constant k ; Steady-state: =𝟎
𝝏𝒕
or
3:29 AM 65
Combined 1-D Heat Conduction Equation
An examination of the one-dimensional transient heat
conduction equations for the plane wall, cylinder, and sphere
reveals that all three equations can be expressed in a
compact form as
3:29 AM 66
GENERAL HEAT CONDUCTION EQUATION
Rectangular Coordinates
3:29 AM 67
or
………..(1)
3:29 AM 68
………..(2)
Substituting into Eq. (1), we get
………..(3)
………..(4)
3:29 AM 69
Noting that the heat transfer areas of the element for heat
conduction in the x, y, and z directions are A 𝒙= ∆y∆z, A 𝒚 = ∆x∆z,
and A 𝒛 = ∆y∆x, respectively, and taking the limit as Δx, Δ y, Δ z
and Δ t → 0 yields
3:29 AM 70
Hence, the general heat conduction equation in rectangular
coordinates is
For constant k,
3:29 AM 71
Steady-state (called the Poisson equation)
3:29 AM 72
The general heat conduction equation in rectangular
coordinates is
For constant k,
3:29 AM 73
Steady-state (called the Poisson equation)
3:29 AM 74
General heat conduction equation in Cylindrical Coordinates
𝒙 = 𝒓 𝐜𝐨𝐬 ∅
𝒚 = 𝒓 𝐬𝐢𝐧 ∅
𝒛=𝒛
3:29 AM 75
General heat conduction equation in Spherical Coordinates
𝒙 = 𝒓 𝐜𝐨𝐬 ∅ sin 𝜽
𝒚 = 𝒓 𝐬𝐢𝐧 ∅ sin 𝜽
𝒛 = 𝒛 cos 𝜽
3:29 AM 76
BOUNDARY AND INITIAL CONDITIONS
……….
3:29 AM 77
3:29 AM 78
Reconsidering the brick wall. The temperature at any point
on the wall at a specified time also depends on the condition
of the wall at the beginning of the heat conduction process.
Such a condition, which is usually specified at time t = 0, is
called the initial condition
3:29 AM 79
3:29 AM 80
Specified Heat Flux Boundary Condition
3:29 AM 81
For a plate of thickness L subjected to heat flux of 50 W/m2
into the medium from both sides, the specified heat flux
boundary conditions can be expressed as
3:29 AM 82
Insulated Boundary
3:29 AM 83
Thermal Symmetry
3:29 AM 84
Convection Boundary Condition
3:29 AM 85
where h1 and h2 are the convection heat transfer coefficients
and T∞1 and T∞2 are the temperatures of the surrounding
mediums on the two sides of the plate
3:29 AM 86
Radiation Boundary Condition
3:29 AM 87
where ε1 and ε2 are the emissivities of the boundary surfaces,
σ=5.67×108 W/m2 · K4 is the Stefan–Boltzmann constant, and
Tsurr, 1 and Tsurr, 2 are the average temperatures of the surfaces
surrounding the two sides of the plate, respectively.
3:29 AM 88
Interface Boundary Conditions
3:29 AM 90
Basic steps involved in the solution of heat transfer
problems
3:29 AM 91
Example
Consider a large plane wall of
thickness L = 0.2 m, thermal
conductivity k = 1.2 W/m°C, and
surface area A= 15 m2. The two
sides of the wall are maintained at
constant temperatures of T1 =120°C
and T2 =50°C, respectively.
Determine:
(a) the variation of temperature
within the wall and the value of
temperature at x=0.1 m
(b) the rate of heat conduction
through the wall under steady
conditions.
3:29 AM 92
Solution
Assumptions
1. Heat conduction is steady.
2. Heat conduction is one dimensional since the wall is large
relative to its thickness and the thermal conditions on both
sides are uniform.
3. Thermal conductivity is constant.
4. There is no heat generation.
3:29 AM 93
Properties
The thermal conductivity is given to be k =1.2 W/m · °C.
Analysis
(a) Taking the direction normal to the surface of the wall
to be the x-direction, the differential equation for this
problem can be expressed as
0 0
Boundary Conditions
Hence,
3:29 AM 94
3:29 AM 95
3:29 AM 96
The first boundary condition can be interpreted as in the
general solution, replace all the x’s by zero and T (x ) by T1
𝑻 𝟎 = 𝑪𝟏 × 𝟎 + 𝑪𝟐 → 𝑪𝟐 = 𝑻𝟏
3:29 AM 97
Substituting the C1 and C2 expressions into the general
solution, we obtain
3:29 AM 98
(b) The rate of heat conduction anywhere in the wall is
determined from Fourier’s law to be
3:29 AM 99
HEAT GENERATION IN A SOLID
Heat generation is usually expressed per unit volume of the
medium, and is denoted by 𝒈ሶ , whose unit is W/m3.
3:29 AM 100
Consider a solid medium of surface area As, volume V, and
constant thermal conductivity k, where heat is generated at a
constant rate of 𝒈ሶ per unit volume. Heat is transferred from
the solid to the surrounding medium at T, with a constant
heat transfer coefficient of h. All the surfaces of the solid are
maintained at a common temperature Ts. Under steady
conditions, the energy balance for this solid can be expressed
3:29 AM 101
Disregarding radiation (or incorporating it in the heat transfer
coefficient h), the heat transfer rate can also be expressed from
Newton’s law of cooling as
3:29 AM 102
a long solid cylinder of radius ro
3:29 AM 103
For a long solid cylinder of radius ro
3:29 AM 104
where To is the centerline
temperature of the cylinder, which
is the maximum temperature, and
ΔTmax is the difference between the
centerline and the surface
temperatures of the cylinder,
which is the maximum temperature
rise in the cylinder above the
surface temperature.
3:29 AM 105
3:29 AM 106
EXAMPLE
Consider a long resistance wire of radius r1 = 0.2 cm and thermal
conductivity kwire=15W/m·°C in which heat is generated uniformly as a
result of resistance heating at a constant rate of 50 W/cm3. The wire is
embedded in a 0.5-cm-thick layer of ceramic whose thermal conductivity is
kceramic=1.2 W/m·°C. If the outer surface temperature of the ceramic layer is
measured to be Ts=45°C, determine the temperatures at the center of the
resistance wire and the interface of the wire and the ceramic layer under
steady conditions.
3:29 AM 107
SOLUTION
The surface and interface temperatures of a resistance wire
covered with a ceramic layer are to be determined.
Assumptions
1. Heat transfer is steady since there is no change with time.
2. Heat transfer is one-dimensional since this two-layer heat
transfer problem possesses symmetry about the centerline
and involves no change in the axial direction, and thus T =T
(r ).
3. Thermal conductivities are constant.
4. Heat generation in the wire is uniform.
3:29 AM 108
3:29 AM 109
3:29 AM 110
3:29 AM 111
T T T T T T T
k + k + k + e&gen + v = cˆ p +u +v +w
x x y y z z t x y z
1 T 1 T T T T v T T
k r + 2 k + k +
gen
&
e + = ˆ
c p + v + + v
r r r r z z t r r z
v r z
T T v T v T 1 2 T 1 T 1 2T
cp
ˆ + vr + + = k r + k sin +
t r r r sin r 2 r r r 2 sin r 2 sin 2
+ e&gen + v
3:29 AM 112
3:29 AM 113
3:29 AM 114
3:29 AM 115
3:29 AM 116
VARIABLE THERMAL CONDUCTIVITY, k (T )
3:29 AM 117
STEADY HEAT CONDUCTION
we have dT/dx = constant, which means that the temperature through the wall varies
linearly with x. That is, the temperature distribution in the wall under steady conditions is
a straight line
Separating the variables in the above equation and integrating from x =0, T(0)= T1, to x = L,
where T(L) = T2, we get
3:29 AM 118
Analogy between thermal and electrical resistance concepts
3:29 AM 119
In case of convection heat transfer
3:29 AM 120
COMBINED HEAT TRANSFER COEFFICIENT
3:29 AM 121
3:29 AM 122
3:29 AM 123
where U is the overall heat transfer coefficient.
3:29 AM 124
Interfacial temperature
3:29 AM 125
Example
3:29 AM 126
3:29 AM 127
3:29 AM 128
THERMAL CONTACT RESISTANCE
3:29 AM 129
A typical experimental setup for the determination of thermal contact
resistance
3:29 AM 130
Equivalent Thickness for Contact Resistance
3:29 AM 131
SOLUTION
The thickness of the aluminium plate whose thermal resistance is equal to the
thermal contact resistance is to be determined.
Properties
The thermal conductivity of aluminium at room temperature is k = 237 W/m·°C
3:29 AM 132
HEAT CONDUCTION IN CYLINDERS
3:29 AM 133
HEAT CONDUCTION IN SPHERES
3:29 AM 134
CRITICAL RADIUS OF INSULATION
(a) CYLINDER
𝑸ሶ 𝜶 𝑨 ;
𝐀𝜶𝒓
3:29 AM 135
3:29 AM 136
3:29 AM 137
Example
3:29 AM 138
3:29 AM 139
the interface temperature can be determined as
3:29 AM 140
HEAT TRANSFER FROM FINNED SURFACES
3:29 AM 141
In the special case of constant cross section and constant thermal
conductivity, the above differential equation can reduces to
Where,
3:29 AM 142
Therefore, the general solution of the differential equation
3:29 AM 143
Negligible Heat Loss from the Fin Tip (Adiabatic fin tip)
(Insulated fin tip, Q fin tip = 0)
3:29 AM 144
Fin Efficiency
3:29 AM 145
Fin Effectiveness
3:29 AM 146
We can draw several important conclusions from the fin effectiveness
relation above for consideration in the design and selection of the fins:
• The ratio of the perimeter to the cross-sectional area of the fin p/Ac should
be as high as possible.
3:29 AM 147