Study Material - Heat and Mass Transfer MODULE 5-MODULE - 55 PDF
Study Material - Heat and Mass Transfer MODULE 5-MODULE - 55 PDF
Study Material - Heat and Mass Transfer MODULE 5-MODULE - 55 PDF
Module V: Radiation: Nature of thermal radiation-definitions and concepts- monochromatic and total emissive power-
Intensity of radiation- solid angle- absorptivity, reflectivity and transmissivity-Concept of black body- Planck’ law- Kirchoff’s
law- Wein’s displacement law-Stefan Boltzmann’s law- black, gray and real surfaces-Configuration factor (derivation for
simple geometries only)- Electrical analogy- Heat exchange between black/gray surfaces- infinite parallel plates, equal and
parallel opposite plates-perpendicular rectangles having common edge- parallel discs (simple problems using charts and
tables). Radiation shields (no derivation).
RADIATION
1.1 INTRODUCTION
Radiation heat transfer is defined as “the transfer of energy across a system boundary by means of
an electromagnetic mechanism which is caused solely by a temperature difference. Whereas the heat,
transfer by conduction and convection takes place only in the presence of a medium, radiation heat
transfer does not require a medium. Radiation exchange, in fact, occurs most effectively in vacuum.
Further, the rate of heat transfer by conduction and convection varies as the temperature difference to the
first power. Whereas the radiant heat exchange between two bodies depends on the difference between
their temperature to the 'forth power'. Both the amount of radiation and the quality of radiation depend
upon temperature.
Examples: The dissipation from the filament of a vacuum tube or the heat leakage through the
evacuated walls of a thermos flask are some familiar examples of heat transfer by radiation. The
contribution of radiation to heat transfer is very significant at high absolute temperature levels such as
those prevailing in furnaces, combustion chambers, nuclear explosions and in space applications. The
solar energy incident upon the earth is also governed by the laws of radiation.
The energy which a radiating surface releases is not continuous but is in the form of successive
and separate (discrete) packet or quanta of energy called photons. The photons are propagated through
space as rays; the movement of groups of photons is described as electromagnetic waves. All types of
electromagnetic waves are classified in terms of wavelength and are propagated at the speed of light (c)
i.e.3 × 108 m/s. The electromagnetic spectrum in shown in Fig. 1.1 The distinction between one form of
radiation and another lies only in its frequency (f) and wavelength (λ) which are related by,
c = λ× f (1 . 1)
The emission of thermal radiation (range
lies between wavelength of 10 -7 m and 10 -4 m)
depends upon the nature, temperature and state of
the emitting surface; however, with gases the
dependence is also upon the thickness of the
emitting layer and the gas pressure. Thermal
radiations exhibit characteristics similar to those of
visible light, and follow optical laws. These can
reflect, refracted and are subject to scattering and Fig. 1.1. Spectrum of electromagnetic radiation.
absorption when they pass through a media. They get polarized and weakened in strength with Inverse
Square of radial distance from the radiating surface.
1.2 SURFACE EMISSION PROPERTIES
The rate of emission of radiation by a body depends upon (i) The temperature of the surface, (ii)
The nature of the surface and (iii) The wavelength or frequency of radiation.
1.1.1. The parameters which deal with the surface emission properties:
(i) Total emissive power (E): The emissive power is defined as the total amount of radiation emitted
by a body per unit area and time. It is expressed in W/m2. The emissive power of a black body,
according to Stefan-Boltzmann, is proportional to absolute temperature to the fourth power.
Eb = σ T 4 W / m2 or E b = σ AT 4 W (1 . 2)
Where σ = Stefan-Boltzmann constant = 5.67 × 10-8W/m2K4
(iii) Emission from real surface-emissivity: The emissive power from a real surface is given by,
E = ε σ AT 4 W (1 . 4)
Where ε = Emissivity of the material.
Emissivity (ε): It is defined as the ability of the surface of a body to radiate heat. It is also defined
as the ratio of the emissive power of anybody to the emissive power of a black body of equal
temperature (i.e., ε = E/Eb). Its values vary for different substances ranging from 0 to 1. For a
black body ε = 1, for a white body surface ε = 0 and for gray bodies it lies between 0 and 1. It
may vary with temperature or wavelength.
(iv) Intensity of radiation or Irradiation (G): it is defined as the total incident radiation on a surface
from all directions per unit time and per unit area of surface: expressed W/m2.
(v) Radiation density and pressure.
(vi) Radiosity (J).It refers to all of the radiant energy leaving a surface.
(vii) Interrelationship between surface emission and irradiation properties.
1.3ABSORPTIVITY, REFLECTIVITY AND TRANSMISSIVITY
When incident radiation Impinges on a surface, three things happens; a part is reflected back (Gr)
apart a transmitted through (Gt) and the remainder is absorbed (Ga), depending upon the characteristics
of the body, as shown in Fig. 1.2.
fraction of incident radiation transmitted). When the incident radiation is absorbed, it is converted into
internal energy.
1.2.1. Black body
For perfectly absorbing body, α=1, ρ = 0, τ =0. Such a body is called a ‘black body’. Black body
is one which neither reflects nor transmits any part of the incident radiation but absorbs all of it. In
practice, a perfect black body (α = 1) does not exists. However its concept is very important.
1.2.2. Opaque body:
When no incident radiation is transmitted through the body, it is called an opaque body. For the
opaque body τ = 0, and Eqn (1.4) reduces to
α + ρ = 1 (1 . 6)
A1E1=α 1 A1E b (1 . 9)
Now we remove body (1) and replace it by body (2) having absorptivity α2. The radiative energy
impinging on the surface of this body is again Eb. In this case, we may write
A2 E 2=α 2 A2 E b (1 . 1 0)
By considering generality of bodies, we obtain
A1E1 A2 E 2 A E E1 E2 E
Eb = = = or Eb = = = (1 . 1 1)
α 1 A1 α 2 A2 α A α1 α2 α
Also, as per definition of emissivity ε,
E1 E1
ε= r Eb = (1 . 1 2)
Eb ε
By comparing Eqns. (1.11) and (1.12), ε = α (1 . 1 3)
α is always smaller than 1. Therefore, the emissive power E is always smaller than the emissive
power of a black body at equal temperature).
Thus, Kirchhoff’s law also states that the emissivity of a body is equal to its absorptivity when the
body remains in thermal equilibrium with its surroundings
1.7 PLANCK’S LAW
In 1900 Max Planck showed by quantum arguments that the spectral distribution of the radiation
intensity of a black body is given by
2π c 2 hλ − 5
( Eλ )b = (1 . 1 4)
ch
e x p −1
λk T
Where (Eλ)b = monochromatic (single wavelength) emissive power of a black body in
W/m2.µm
c = velocity of light is vacuum, 2.998 × 108≈ 3 × 108m/s
Dept. of Mechanical Engineering S.N.M.I.M.T. Page | 4
h = Planck’s constant = 6.625 × 10-34j.s
λ = wavelength, µm [1 µm (micro meter) = 10-6 m]
k = Boltzmann constant = 1.3805 × 10-23J/K
T = absolute temperature, K
Hence the unit of (Eλ)b is W/m2.µm
Quite often the Planck’s law is written as
C1λ − 5
( Eλ )b = (1 . 1 5)
C
e x p 2 −1
λT
Where, C1 = 2π c 2 h = 3 . 7 4 2 × 1 0 8 W .µ m 4 / m 2
ch
and C2 = = 1 . 4 3 8 8 ×1 0 4 µ m K
k
Equation (1.14) is of great importance as it provides from a black body. The quantity (Eλ)b,
monochromatic emissive power, is defined as the energy emitted by the black surface in all directions at a
given wavelength λ per unit wavelength interval around λ; that is, the rate of energy emission in the
interval dλ is equal to (Eλ)b dλ. The total emissive power and monochromatic emissive power are related
by the equation,
∞
Eb = ∫ ( E λ ) bd λ (1 . 1 6)
0
This integral represents the total emissive power per unit area radiated from a black body.
1.8 WIEN’S DISPLACEMENT LAW
In 1893 Wien’s established a relationship between the temperature of a black body and the
wavelength at which the maximum value of monochromatic emissive power occurs. A peak
monochromatic emissive power occurs at a particular wavelength. Wien’s displacement law states that
the product of λmax and T is constant, i.e.,
λ m a xT = C o n s t a n t (1 . 1 7)
θ =π 2 θ =π 2
Q = 2π I d A1 ∫ s i n θ .c o s θ .dθ = π I d A1 ∫ 2s i n θ .c o s θ .dθ
θ =0 θ =0
θ =π 2
Q = π I d A1 ∫ s i n 2θ .dθ = π I d A1 (1 . 2 2)
θ =0
d A2 .c o s θ 2
d Q1− 2 = I b1.d A1.c o s θ 1.
r2
I c o s θ 1.c o s θ 2 .d A1.d A2 .
d Q1−2 = b1 (2 . 2)
r2
Dept. of Mechanical Engineering S.N.M.I.M.T. Page | 8
This energy is absorbed by the elementary area dA2, since both the surfaces are black. The
quantity of energy radiated by dA2and absorbed by dA1is given by,
d A1.c o s θ 1
d Q 2−1 = I b 2 .d A2 .c o s θ 2 .d ω 2 = I b 2 .d A2 .c o s θ 2 .
r2
I b 2 .c o s θ 1.c o s θ 2 .d A1.d A2
d Q 2−1 = (2 . 3)
r2
Net heat exchange between dA1 and dA2
d Q1 2 = d Q1− 2 − d Q 2−1
I b1.c o s θ 1.c o s θ 2 .d A1.d A2 . I b 2 .c o s θ 1.c o s θ 2 .d A1.d A2 .
d Q1 2 = −
r2 r2
d A .d A2 .c o s θ 1.c o s θ 2
= ( I b1 − I b 2 ) 1
r2
σ .T14
E b1 σ .T 24
E b2
But, I b1 = = and I b2 = = (2 . 4)
π π π π
(
d Q1 2 = σ T1 4 − T 2 4 ) d A .d A .cπor s θ .c o s θ
1 2
2
1 2
(2 . 5)
The net heat transfer is between the Area A1 and A2 is given by,
(
Q1 2 = ∫ d Q1 2 = σ T1 4 − T 2 4 ) ∫ ∫ c o sθ .c oπsrθ
1
2
2 .d A1.d A2
(2 . 6)
A1 A2
The rate of radiant energy emitted by A1 that falls on A2, from Eqn. (12.2), is given by
c o s θ 1.c o s θ 2 .d A1.d A2 c o s θ 1.c o s θ 2 .d A1.d A2
Q1− 2 = ∫ d Q1− 2 = I b1 ∫ ∫ 2
= σ T1 4 ∫ ∫ (2 . 7)
A1 A2 r A 1 A2
πr2
Q1− 2
= F1− 2 (2 . 8 a )
Q1
F1-2 is known as ‘configuration factor or surface factor or ‘view factor between the two radiating
surfaces and is a function of geometry only. Thus, the shape factor may be defined as “The fraction of
radiative energy that is diffused from one surface element and strikes the other surface directly with no
intervening reflections.”
Q1 2 = Q1F1 2 or Q1 2 = F1 2σ . A1.T14 (2 . 9)
Similarly, the rate of radiant energy by A2 that falls on A1from Eqn. (12.3), is given by
c o s θ 1.c o s θ 2 .d A1.d A2 c o s θ 1.c o s θ 2 .d A1.d A2
Q 2−1 = ∫ d Q 2−1 = I b 2 ∫ ∫ 2
= σ T2 4 ∫ ∫
A1 A2 r A 1 A2
πr2
Q 2−1
= F2−1 and Q 2−1 = F2−1.Q 2 = F2−1. A2 .σ .T2 4 (2 . 1 1)
Q2
F2-1is the shape factor of A2with respect to A1
From (2.8) and (2.10),
A1.Q1−2 A .Q A1.F1− 2 .Q1 A2 .F2−1.Q 2
= 2 2−1 or =
Q1 Q2 Q1 Q2
A1.F1− 2 = A2 .F2−1 (2 . 1 2)
The above result is known as reciprocity theorem. It indicates that the net radiant interchange may
be evaluated by computing one way configuration factor from either surface to the other. Thus the net
rate of heat transfer between two surfaces A1 and A2is given by
( ) (
Q1 2 = A1.F1−2 .σ . T1 4 − T2 4 = A2 .F2−1.σ . T1 4 − T2 4 ) (2 . 1 3)
It may be noted that Eqn. (12.13) is applicable to black surfaces only and must not be used for
surfaces having emissivities very different from unity. The evaluation of the integral equation (12.8) for
determining the shape factor is rather complex and cumbersome. Therefore, results have been obtained
and presented in graphical form for the geometries normally encountered in engineering practice.
Geometrical factors for parallel planes (discs and rectangles) directly opposed and those for radiation
between perpendicular rectangles with a common edge are shown in Figs. 2.2.
Fig. 2.2. Shape factor for perpendicular rectangles with a common edge.
2.3 SHAPE FACTOR ALGEBRA AND SALIENT FEATURES OF THE SHAPE FACTOR
In order to compute the shape factor for certain geometric arrangements for which shape factors
or equations are not available, the concept of shape factor as fraction of intercepted energy, and
reciprocity theorem can be used. The shape factors for these geometries can be derived in term of known
shape factors of other geometries. The interrelation between various factors is called shape factor
algebra.
For the calculation of shape factors for specific geometries and for the analysis of radiant heat
exchange between surfaces, the following facts and properties will the useful:
1. The shape factor is purely a function of geometric parameters only.
In general, Ai .Fi − j = A j .F j −i
This reciprocal relation is particularly useful when one of the shape factors is unity.
3. When all the radiation emanating from a convex surface 1 is intercepted by the enclosing surface
2, the shape factor of convex surface with respect to the enclosure F1-2 is unity. Then in
conformity with reciprocity theorem, the shape factor F2-1is merely the ratio of areas.
ie. When surface A1is entirely convex, say a sphere, completely enclosed by A2, then according
to reciprocity relation, we have
A1.F1−2 = A2 .F2−1 and A1 = A2 .F2−1
[∵ F1−2 = 1, a s s u r f a c e1 c o m p l e t e l y s e e s s u r f a c e 2]
A1
Or F2−1 = (i e. , r a t i o o f a r e a s ) , a n d F2−1 + F1− 2 = 1
A2
(
I n t h i s c a s e , t h e b l a c k b o d y r a d i a t i o n e x c h a n g e i s , Q1 2 = A1.σ . T1 4 − T2 4 )
4. A concave surface has a shape factor with itself because the radiant energy coming out from one
part of the surface is intercepted by another part of the same surface. The shape factor of a
surface with respect to itself is F1-1
5. For a flat or convex surface, the shape factor with respect to itself is zero (i.e., F1-1= 0).This is
due to the fact that for any part of flat or convex surface, one cannot see/view any other part of
the same surface.
Q1 = E1 − α 1 (1 − α 2 ) E1[1 + (1 − α 1 ) (1 − α 2 ) + (1 − α 1 ) 2 (1 − α 2 ) 2 + . . . . ]
Q1 = E1 − α 1 (1 − α 2 ) E1[1 + Z + Z 2 + . . . . ]
W h e r e , Z = (1 − α 1 ) (1 − α 2 )
ε1 ε 2 1
W h e r e , f 1−2 = Or f 1− 2 =
ε 1 + ε 2 − ε 1ε 2 (1 ε 1 ) + (1 ε 2 ) − 1
E1ε 2 E 2ε 1. ( A1 A2 )
Q1 2 = A1 − A2
( A1 A2 ) ε 1 + ε 2 − ( A1 A2 ) ε 1ε 2 ( A1 A2 ) ε 1 + ε 2 − ( A1 A2 ) ε 1ε 2
A1E1ε 2 − A1E 2ε 1
Q1 2 =
( A1 A2 ) ε 1 + ε 2 − ( A1 A2 ) ε 1ε 2
Q1 2 =
A1ε 1ε 2σ T1 4 − A1ε 1ε 2σ T 2 4
or Q1 2 =
(
ε 1ε 2 A1σ T1 4 − T 2 4 )
( A1 A2 ) ε 1 + ε 2 − ( A1 A2 ) ε 1ε 2 1
( A1 A2 ) ε 1ε 2 − 1 + ε 2
ε2
Q1 2 =
(
A1σ T1 4 − T2 4 ) (2 . 2 4)
1 A 1
+ 1 − 1
ε 1 A2 ε 2
(
Q1 2 = f 1− 2 A1σ T1 4 − T 2 4 ) (2 . 2 4 a)
1
W h e r e , f 1− 2 =
1 A1 1
+ − 1
ε 1 A2 ε 2
f1-2 is called interchange factor for the radiant heat exchange between infinite long concentric
cylinders. The equation (2.24) is equally applicable to concentric spheres except that for concentric
cylinders of equal length.
2.4.3 Small Gray Bodies
Consider two small gray bodies (the gray bodies are said to be small if their size is very small as
compared to the distance between them) 1 and 2 having emissivities ε1andε2, or absorptivities α1and α2.
The radiant energy emitted by body 1 is absorbed by body 2. The portion of radiation unabsorbed and
thus reflected on the first incidence is considered to be lost in space (due to surface being small) ie.
Nothing is return back to surface 1. Similar in this case with surface 2,
The net heat energy transfer between two bodies, Q1 2 = ε 1ε 2 F1− 2 A1σ T1 4 − ε 1ε 2 F2−1 A2σ T 2 4
But A1F1− 2 = A2 F2−1 ……..By reciprocal theorem
(
Q1 2 = ε 1ε 2 A1F1− 2σ T1 4 − T2 4 ) or (
Q1 2 = f 1− 2 A1F1−2σ T1 4 − T 2 4 ) (2 . 2 5)
Where, f1-2 is called equivalent emissivity or interchange factor for the radiant heat exchange
between two small gray bodies.
The net radiant heat exchange between the small body 1 and the outer large enclosure 2,
Q1− 2 =ε 1 A1σ T1 4 − ε 1ε 2 A2 F2−1σ T 2 4
I f T1 = T 2 a n d Q1 2 = 0, t h e n A1 = A2ε 2 F2−1
(
Q1− 2 =ε 1 A1σ T1 4 − T 2 4 ) or (
Q1− 2 = f 1-2 A1σ T1 4 − T2 4 ) (2 . 2 6)
(
Q n e t = f 1-2 F1− 2 A1σ T1 4 − T2 4 ) (2 . 2 7)
2.5 ELECTRICAL NETWORK ANALOGY FOR THERMAL RADIATION SYSTEMS
An electrical network analogy is an alternative approach for analysing radiation heat exchange
between gray or black surfaces. In this approach the two terms commonly used are irradiation and
radiosity.
2.5.1 Irradiation (G)
It is defined as the total radiation incident upon a surface per unit time per unit area; it is
expressed in W/m2.
2.5.2Radiosity (J)
This term is used to indicate the total radiation leaving a surface per unit time per unit area. It is
also expressed in W/m2.
The radiosity comprises the original emittance φ from the surface plus the reflected portion of
any radiation incident upon it [Refer Fig. 2.7(a)].
i e . J = E + ρG
Or J = ε E b + ρG (2 . 2 8)
WhereEb= emissive power of a perfect blackbody at
the same temperature.
Also α +ρ+τ = 1
Or α +ρ = 1 (τ = 0, t h e s u r f a c eb e i n g o p a q u e)
Or ρ = 1 −α
∴ J = ε E b + (1 − α )G
But α =ε . . . .B y K i r c h h o f f ' s l a w
∴ J = ε E b + (1 − ε )G (2 . 2 9)
Fig. 2.7. Irradiation and radiosity
Or J − ε E b = (1 − ε )G
J − ε Eb
Or G= (2 . 3 0)
(1 − ε )
The net energy leaving a surface is the difference between its radiosity and irradiation. Thus,
Qn e t Qn e t J − ε Eb
= J −G T h e n, =J− (2 . 3 1)
A A (1 − ε )
Q n e t J (1 − ε ) − ( J − ε E b ) Qn e t J − ε J − J + ε Eb ε ( E b − J )
= Or = =
A (1 − ε ) A 1− ε 1− ε
Qn e t J − ε J − J + ε Eb ε ( E b − J )
Or = =
A 1− ε 1− ε
Aε ( E b − J )
∴ Qn e t = (2 . 3 2)
1− ε
∴ Q1 2 =
( J1 − J 2 ) (2 . 3 4)
1 A1F1−2
This eqn. (2.34) can be represented in the form of electric network as shown in Fig. 2.8. The
quantity 1 A1F1−2 is called the space resistance because it is due to the distance and geometry of the
radiating bodies.
If the surface resistances of the two bodies and spaceresistance between them are considered then
the net heat flow can be represented by an electric circuit as shown in Fig. 2.9. The net heat exchange
between the two gray surfaces is given by,
E b1 − E b 2 A1 σ (T1 4 − T 2 4 )
( Q1 2 ) n e t =
1− ε1 1 1− ε 2
or ( Q1 2 ) n e t =
1− ε1 1 1 − ε 2 A1
+ + + + .
A1ε 1 A1F1−2 A2ε 2 ε1 F1−2 ε 2 A2
or ( Q1 2 ) n e t = ( Fg ) 1−2 A1 σ (T1 4 − T2 4 ) (2 . 3 5)
1
Where, ( Fg )1−2 = 1 − ε 1 1 1 − ε 2 A1
+ + .
ε1 F1− 2 ε 2 A2
(Fg)1-2is known as gray body factor. When exchange of heat takes place between two black
surfaces, the surface resistance becomes zero as ε1 = ε2 = 1; (Fg)1-2changes to F1-2(i.e., the configuration
factor) and the eqn.(2.35) reduces to
Q n e t = F1− 2 A1 σ (T1 4 − T 2 4 ) . . . . .F o r b l a c k s u r f a c e s
( Fg )1−2 = 1 − ε11 (2 . 3 7)
+1
ε1
Examples: (i) A pipe carrying steam in a large room, and (ii) A thermocouple bead located inside
a duct to measure temperature of the fluid.
With no radiation shields, the net heat exchange between the parallel planes is given by:
A σ (T1 4 − T2 4 )
( Q1 2 ) n e t =
1 1
(2 . 2 1)
+ −1
ε1 ε2
If the emissivity of the radiation shield is ε3, the equation is to find heat exchange between
surfaces 1, 3 and 3, 2.
A σ (T1 4 − T3 4 )
( Q1 3 ) n e t =
1 1
(2 . 3 8)
+ −1
ε1 ε3
1 1 1 1
T1 4 + − 1 + T 2 4 + − 1
T3 4 = ε3 ε2 ε1 ε 3 (2 . 4 1)
1 1 1 1
+ − 1 + + − 1
ε 3 ε 2 ε1 ε 3
Substituting the value of T3in the left hand side of eqn. (2.40), we get
A σ (T1 4 − T 2 4 )
( Q1 2 ) n e t = (2 . 4 2)
1 1 1 1
+ − 1 + + − 1
ε1 ε 3 ε 3 ε 2
Dividing eqn. 2.42 by eqn2.21, we get
1 1
( Q1 2 ) + −1
net withshield ε1 ε 2
= (2 . 4 3)
( Q1 2 ) 1 1 1 1
net withoutshield
+ − 1 + + − 1
ε1 ε 3 ε 3 ε 2
In the general case where there are n shields, all the surface resistances would be the same, since
the emissivities are equal. There will be two surface resistances for each shield and one for each heat
transfer surface. There will also be (n + 1) ‘space resistances' but the configuration factor is unity for each
infinite parallel plane.
1− ε
T o t a l R e s i s t a n c e ( R ) n - s h i e l d = ( 2 n + 2 ) + ( n + 1)(1) A (2 . 4 6)
ε
2
= ( n + 1) − 1 + ( n + 1)(1) A (2 . 4 7)
ε
The radiant heat transfer rate between two infinitely large parallel plates separated by n-shields is,
therefore,
Then ( Q ) w i t h o u t s h i e l d =
(
Aσ T1 4 − T2 4 ) (2 . 5 0)
( ( 2 ε ) − 1)
From eqns. 2.48 to 2.50,
( Q ) n− s h i e l d ( R) w i t h o u t s h i e l d 1
= = (2 . 5 1)
( R) w i t h o u t s h i e l d ( R ) n− s h i e l d n +1
From this provided that the reduction in heat flow does not affect the wall temperatures.
Arranging the eqn. (2.43) in a slightly different form,
A σ (T1 4 − T2 4 )
( Q1 2 ) n e t = (2 . 4 3)
1 1 1
+ + 2 − 2
ε1 ε 2 ε3
This equation can be generalized for a system of two parallel plates separated by n-shields of
emissivity εs1, εs2….εsn as
Aσ (T1 4 − T2 4 )
( Q1 2 ) n e t = n
(2 . 4 4)
1 1 1
+ + 2∑ − ( n + 1)
ε1 ε2 i =1 ε si