Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Майнор / Книги / Chapter-8---Thermal-radiation_2021_Heat-Transfer-Engineering

.pdf
Скачиваний:
0
Добавлен:
08.05.2026
Размер:
4.3 Mб
Скачать

8.9 The radiosity-irradiation method 273

τi = 0,1 ≤ i ≤ N

(8.114)

Let all the surfaces be gray and diffuse, which is a reasonable engineering approximation. In the view of this for all i,

εi = αi

(8.115)

We know that

 

 

αi + ρi = 1

 

 

(8.116)

 

 

ρi = 1 αi = 1 εi

 

 

(8.117)

Substituting for ρi in Eq. (8.118), we have

 

 

 

 

A J

i

= A ε

σ T 4

+ A (1 ε

)G

(8.118)

i

i i

i

i

i

i

It is also evident that if we consider the radiation falling on the ith surface, given by Ai Gi, we get

 

N

N

 

 

AiGi = Aj Fji J j = Ai Fij J j

(8.119)

 

j=1

j=1

 

The irradiation G

on the ith surface is given by

 

i

 

 

 

 

 

N

 

 

Gi

= Fij J j

(8.120)

 

 

j=1

Substituting for Gi in Eq. (8.118), we have

 

 

 

N

 

 

Ji = εi σ Ti4 + (1 εi ) Fij J j

(8.121)

 

 

j =1

Eq. (8.121) is a powerful one and gives us a system of N equations in N unknowns. By solving the system of equations, we can obtain J1 , J2 , J3 ,., JN .

Once we know Ji, we can calculate Gi. From Gi and Ji, we can calculate net radiative heat transfer from any surface by Eq. (8.122)

qi = Ji Gi

(8.122)

or, from Eq. (8.122), we can write qi in terms of Gi alone

qi = εi σ Ti

4

+ (1 εi )Gi Gi

(8.123)

qi = εi σ Ti

4

+ Gi

εiGi

Gi

(8.124)

qi = εi

σ Ti4

εiGi

 

(8.125)

274 CHAPTER 8  Thermal radiation

qi = εi (σTi4 Gi )

or, from Eq. (8.122), we can write qi in terms of Ji alone

Ji = εi σ Ti4 + (1 εi )Gj

Gi =

Ji εiσTi4

 

 

(1 εi )

 

qi = Ji

Ji εiσTi 4

 

(1 εi )

 

 

qi =

Ji εi Ji Ji + εiσTi4

 

(1 εi )

 

 

qi =

εi [σTi 4 Ji ]

 

(1 εi )

(8.126)

(8.127)

(8.128)

(8.129)

(8.130)

(8.131)

Eq. (8.131) requires some caution in its use, as it is not applicable for a black surface.

So, if we know Gi or Ji alone, we can calculate the net radiative heat transfer from any surface.

For a black body, εi = 1, and so Ji is given by

Ji = εiσTi4 + (1 εi )Gi

(8.132)

Ji = σTi 4

(8.133)

Another surface of interest in the radiation problem is a reradiating surface, which is the radiation equivalent of an insulated surface.

For a reradiating surface, qi = 0.

 

qi = Ji Gi = 0

(8.134)

 

 

Ji = Gi

(8.135)

Additionally,

 

 

 

 

 

εi

[σTi 4 Ji ] = 0

 

 

 

 

1 εi

 

 

(8.136)

 

 

J

i

= σT 4

(8.137)

 

 

 

i

In the view of the above the key result is,

Ji ≠ f (εi )

(8.138)

8.9 The radiosity-irradiation method 275

The key question to be answered now is: How do we engineer a reradiating surface? A very thick wall/surface made of a low thermal conductivity medium and with a vacuum or near vacuum will come close to a re-radiating condition. Else, conjugate

heat transfer has to be considered in a wall, that is, qconduction = qradiation and we need to write out the expressions for qconduction and qradiation and we will obtain an equation in which unknown will be the surface temperature.

Example 8.12: Consider two infinitely long parallel plates that are separated by a distance L, as shown in Fig. 8.26. The plates are very deep in the direction perpendicular to the plane of the paper. The space between the two plates is evacuated. The two plates are maintained at temperatures of 810 K and 490 K, respectively. The two plates have emissivities of 0.95 and 0.45 respectively. Determine the net radiation heat transfer between the two plates.

Solution:

F11 = 0 and F12 = 1

F21 = 1 and F22 = 0

View factor matrix =

 

0

1

 

 

1

0

 

 

 

 

 

We can solve this problem using the radiosity-irradiation method.

Ji = εiσTi4 + (1 εi )Gi

N

Gi = Fji J j j=1

(8.139)

(8.140)

(8.141)

FIGURE 8.26

Two infinitely long parallel plates undergoing surface radiation exchange and under consideration in example 8.12.

276 CHAPTER 8  Thermal radiation

qi = Ji Gi

(8.142)

Now, we write the radiosity and irradiation expressions for the two surfaces 1 and 2.

(8.143)

(8.144)

J1

= ε1σT14

+ (1 ε1 )J2

(8.145)

J2

= ε2σT24

+ (1 ε2 )J1

(8.146)

Substituting J2 from Eq. (8.146) in Eq. (8.145), we have

J1

= ε1σ T14 + (1 ε1 ) ε2σ T24

+ (1 ε2 )J1

 

 

 

 

 

J1 [1 (1 ε1 )(1 ε2 )] = ε1σ T14 + (1 ε1 )ε2σ T24

J1 =

 

ε1σT14 + ε2σT24 ε1ε2σT24

 

 

 

 

 

1 (1 ε1 )(1 ε2 )

 

q1 = J1 G1 q1 = J1 J2

 

 

ε1ε2

σ T14 σ T24

 

q1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 (1 ε1 )(1 ε2 )

 

 

σ T14

T24

 

 

 

 

q1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

+

 

1

 

1

 

 

 

 

 

 

 

ε1

 

ε2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.67 × 108

8104

490

4

q1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

+

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.95

0.45

 

 

 

 

 

 

 

 

 

 

q1 = 9292.4 W/m2

The expression we derived above is a powerful engineering formula known as the parallel plate formula. This can be used to very swiftly determine the net radiative heat transfer between two infinitely long parallel plates. Double pane windows, radiation shields and several engineering applications find this formula handy.

Example 8.13: Revisit example 8.12. Now a radiation shield with an emissivity of 0.05 on both sides is placed between plates 1 and 2. The details are provided in Fig. 8.27. Determine

1.The net radiation heat transfer in the presence of the shield.

2.The equilibrium temperature of the shield.

8.9 The radiosity-irradiation method 277

FIGURE 8.27

Two infinitely long parallel plates with shield in between under consideration in example 8.13.

Solution:

From the above problem, we have

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

T14

 

T34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q13

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

+

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε1

 

 

 

ε3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

T34

 

T24

 

 

 

 

 

 

 

 

 

 

 

 

q32

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

+

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε3

 

 

 

 

ε2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

At a steady state, q13 = q32 = q12,shield

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

T14

T34

 

 

 

 

 

σ

 

 

T34

T24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

+

 

1

1

 

 

 

 

 

 

 

 

 

+

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε1

 

 

 

ε3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε3

 

ε2

 

 

 

 

 

 

 

8104

 

T34

 

 

 

 

=

 

 

 

 

 

T34 4904

 

 

 

 

 

 

 

1

 

 

+

 

 

1

 

 

 

1

 

 

 

1

 

 

 

 

+

 

1

 

 

1

0.95

 

 

0.05

 

 

0.05

 

 

0.45

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8104

 

T34

 

=

 

 

T34

4904

 

 

 

 

 

 

 

 

 

 

 

 

20.05

 

 

 

 

 

 

 

21.22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21.22 8104

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

+ 20.05

 

490

4

 

T3

=

 

4

 

 

 

 

 

 

 

 

 

21.22

+ 20.05

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T3 706.64K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

T14

T34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

= q13 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

+

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε1

 

ε3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

5.67 × 108 (8104 706.644 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

+

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.95

 

 

 

0.05

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q = 512.15W/m2

278 CHAPTER 8  Thermal radiation

We see that q has drastically reduced (with respect to the solution in example 8.12) by placing just one low emissivity shield. If more are placed, this could help us obtain a “super insulation’. If all the surfaces and “N” shields placed in between all have the same emissivity, it can be shown that, qNshields = qnoshield/(N + 1).

8.10  Introduction to gas radiation

The study of heat transfer through media, which can absorb, emit, and scatter radiation, has been receiving increasing attention in the last few decades. Such media are called radiatively participating media or simply participating media. Applications include emission from rocket nozzles, combustion chambers, radiation to and from the Earth and other planets, glass making, IC engines so on.

Consider solar radiation incident on the Earth’s surface, as shown in Fig. 8.28. The curve with the dotted line corresponds to the intensity of radiation emitted by a black body at 5800 K. This corresponds to the temperature of the photosphere or the outer layer of the Sun. The actual solar irradiation at the top of the atmosphere is shown as a solid line. The inner curve shows the solar irradiation after passing through the Earth’s atmosphere. Here, we can see regions having lower irradiation

FIGURE 8.28

Approximation of variation of incident solar radiation and attenuation by the atmosphere (Adapted from Balaji C., 2014).

8.11 Equation of transfer or radiative transfer equation (RTE) 279

due to high absorption by CO2 and H2O. This figure clearly shows that gas absorption is spectrally dependent and that the atmosphere is truly a participating medium. There are also transparent windows where there is little absorption.

8.11  Equation of transfer or radiative transfer equation (RTE)

Consider a gas volume with a cross sectional area dA and thickness dx as shown in Fig. 8.29. The incoming spectral radiation in the direction x is Iλ,x while Iλ, x+dx is the spectral radiation exiting the gas volume. The area dA is normal to the direction x such that the radiation is traveling in a direction normal to the cross-sectional area. When we neglect scattering, the only two other phenomena involved in this problem are emission and absorption.

The energy absorbed by the gas in the interval dλ is given by

κ λ Iλ dAdx

(8.147)

where kλ is the monochromatic or spectral absorption coefficient (m1). The energy emitted by the gas volume is given by

ελ Ib,λ (Tg )dAdx

(8.148)

where ελ is the monochromatic or spectral emission coefficient (m1). By energy balance, we have

dIλ = ελ Ibλ (Tg ) κ λ Iλ dx

On rearranging, we get the following equation

dIλ +κ λ Iλ = ελ Ib,λ (Tg ) dx

(8.149)

(8.150)

Eq. (8.150) is known as the RTE, or radiative transfer equation. It is deceptively innocuous as ελ and αλ invariably have strong spectral dependence. Additionally, if “reflection ” or scattering from gas molecules is concerned, the problem gets even more formidable, with scattering being a function of wavelength and angle in many cases.

FIGURE 8.29

Gas volume used in the derivation of the RTE.

280 CHAPTER 8  Thermal radiation

Consider the gas under question in an isothermal cavity at Tg. Let the temperature of the gas be Tg. For this situation we have isotropic radiation; that is, radiation intensity is the same everywhere and so dIλ /dx is 0 everywhere within the isothermal cavity. This can be written as

 

 

Iλ = Ib,λ (Tg )

(8.151)

Substituting in the above Eq. (8.150)

 

0 +κ λ Ib,λ (Tg ) = ελ Ib,λ (Tg )

(8.152)

 

 

κ λ = ελ

(8.153)

The above is actually Kirchhoff’s law.

 

Substituting for ελ in Eq. (8.150)

 

 

dIλ

+κ λ Iλ = κ λ Ib,λ (Tg )

(8.154)

 

 

 

dx

 

Let us consider a situation where the absorption is much larger compared to emission. Consider a plane wall at temperature Tw (corresponding to x = 0) surrounded by the gas that is emitting and absorbing. In view of the foregoing assumption, the RTE becomes

dIλ +κ λ Iλ = 0 dx

or

dIλ = −κ λ dx Iλ

If kλ f (λ) , that is, the gas is gray,

dII = −κ dx

If we integrate this, we get I = Cekx . At x = 0, I = Iw, i.e. the wall intensity.

C = Iw, Now the solution becomes I = Iw ekx.

(8.155)

(8.156)

(8.157)

The above solution clearly shows that radiation decreases exponentially with x. This is known as Beer-Lambert’s law. A typical example is solar radiation decreasing exponentially from the surface as we go deeper into the ocean waters.

8.11.1  Determination of heat fluxes

The heat flux q is of fundamental importance to engineers. Hence, the intensity I or Iλ, to be more general, must be converted into a flux form.

Now let us consider the full spectrum of electromagnetic radiation, that is, λ = 0 to λ = . The RTE becomes

8.11 Equation of transfer or radiative transfer equation (RTE) 281

 

dI

+ kI

 

= kI

 

(T )

 

 

λ

b

 

dx

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

Now Ib

=

 

EbT

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

=

σ Tg4

 

 

b

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

The solution to Eq. (8.158) is given by

 

+

 

 

 

4

 

kx

 

σ T 4

 

 

kx

 

 

 

 

 

 

 

 

kx

 

σ Tw

 

cos θ

 

 

 

I

 

 

 

=

 

e

 

+

 

1

− ecos θ

 

 

 

 

 

cosθ

 

π

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.158)

(8.159)

(See Balaji, 2014), for a derivation of the above result. In the above equation θ is the angle between a straight path of radiation and a path.) The heat flux at x = 0, going out in the positive direction of x is given by

qL+ = 2π

π

 

2

IL+ (cosθ )cosθsinθ dθ dφ

φ =0

θ =0

(8.160)

The integral appearing on the right hand side of Eq. (8.160) is hard to integrate. We can, of course, easily do it on a computer this day and age. Historically, this was solved by an introducing an exponential integral En(t), as defined below. We can introduce the exponential function En(t) here as

En (t) =

1

 

 

t

 

0

n2 exp

 

d

 

 

 

 

 

 

 

(8.161)

 

 

 

 

 

 

The exponential integral of third order is given by

 

(t) =

1

 

 

t

 

E3

0

1exp

 

d

 

 

 

 

 

 

 

 

(8.162)

 

 

 

 

 

 

 

In Eq. (8.162) µ is basically a dummy variable.

E3(x) satisfies the following asymptotic properties.

lim E3

 

1

 

 

(x) =

 

− x

(8.163)

t0

2

 

E3 (0) = 0.5 . The x in Eq. (8.162) corresponds to k.x. The product of k and x is known as the optical depth and is given by τ or τx.

When the optical depth is very small, E3(x) approaches 0. In gas radiation, this is known as the optically thin limit.

qL+ can now be written as

qL+ = σTg4 + 2σ (TW4 Tg4 )E3 (τ L )

(8.164)

 

 

qL+ = 2E3 (τ L )σTW4 + σTg4 [1 2E3

(τ L )]

(8.165)

 

 

282 CHAPTER 8  Thermal radiation

For an optically thin gas,

1

(8.166)

E3 (kL ) = 2 kL

 

and 1 2E3 (kL ) = τ L

(8.167)

qL+ = 2τ LσTg4 + (1 2τ L )σTw4 (8.168)

Therefore the radiation arriving in a plane gas layer of thickness L with the wall at x = 0 being at a temperature TW at x = L consists of two parts. These are

1.Radiation coming directly from the gas (first term).

2.Radiation coming from the wall and that is attenuated by the gas (second term).

We can now introduce two key quantities, namely (1) gas emissivity εg and (2) gas absorptivity τg as

εg = 2τ L = (2L)(K)

(8.169)

 

τ g = 1 2τ L = (1 εg )

(8.170)

 

qL+ = εgσTg4 + τ gσTw4 (8.171)

Finally we are in a position to define the gas emissivity as 2.L.k, where k is the gas absorptivity and 2L represents a kind of mean length traveled by the radiation. This mean length is frequently referred to as the mean beam length.

Consider a hemispherical gas volume whose radius is Ro (see Fig. 8.30). Consider an elemental area dA at the center of the bottom surface. Let the gas volume emit and absorb with no scattering.

Consider an optically thin gas at a temperature Tg. For this situation,

 

σ TW4

(κ x)

 

σ Tg4

(κ x)

 

(8.172)

I(x) =

π

e

 

+

π

(1e

 

)

 

 

 

FIGURE 8.30

Hemispherical gas volume used for elucidating the concept of mean beam length.