Module 3 - Thermal Radiation
Module 3 - Thermal Radiation
Thermal Radiation
Thermal radiation is a process by which energy, in the form of electromagnetic radiation,
emitted by a heated surface in all directions.
Solid and liquid surfaces at all temperatures continuously emit energy in the form of radiation.
The emitted radiation is proportional to the temperature of the body. For radiation, no material or
medium is required and can be transported even in the space and vacuum.
Thermal radiation is similar to other electromagnetic waves like x-rays and 𝛾-rays. The energy of
radiation is emitted from the body at all the wave lengths (frequencies), but the heating effect
occurs only over a limited range of wave lengths. The radiation in this wave length range is
known as thermal radiation.
𝛾-rays: Used in medicine for killing Cancer cells
X-rays: Used to view inside of bodies or objects
UV rays used in fluorescent tubes; absorbed by the skin
Visible light: make the things able to be seen
Infrared: transmits the heat from sun, fires, radiators etc
Micro wave: Used in cooking, radar, telephone and other signals
Radio wave: used to broadcast radio and television
The thermal radiation lies in the range of wave lengths from 0.1 μm to 100μm. The light
radiation which is visible lies in the range of 0.4 μm to 0.7μm. It can be seen that the thermal
radiation overlaps and includes the visible light range. The electric filament bulb is a common
example.
Radiation properties
Absobed radiation
Absorptivity (𝛼) =
Incident radiation
Reflected radiation
Reflectivity (𝜌)=
Incident radiation
Tramsmitted radiation
Transmitivity (𝜏) = Incident radiation
QA + QR + QT =Qi
𝑄𝐴 𝑄𝑅 𝑄𝑇
+ + =1
𝑄𝑖 𝑄𝑖 𝑄𝑖
𝛼+𝜌+𝜏=1
Gray body: If a body absorbs a definite percentage of incident radiation falling on it, irrespective of
their wave length is known as Gray body.
White body: If all the incident radiation falling on the body is reflected, then it is known as White
body.
∴𝛼+𝜌=1
Transparent body: It is one, which transmits a part of the radiation falling it.
Black body: A black body is an object which absorbs all the radiation falling on it.
A black body is an idealized physical body that absorbs all incident electromagnetic radiation,
regardless of frequency or angle of incidence. The name "black body" is given because it absorbs
radiation in all frequencies, not because it only absorbs: a black body can emit radiation.
It is necessary to define a theoretical black body to compare the performances of other bodies. Black
body is one which absorbs all the radiation and none is reflected nor transmitted. In other words a
black body is one whose Absorptivity is 100%.
In practice a perfect black body does not exist, but the concept of a black body is very important in
radiation heat transfer. A black body is also a perfect emitter.
Cavity acts as a black body. A cavity with a small hole is shown in figure. An incident ray through
opening is emitted many times on the inside surface before escaping finally from the opening. With
each emission, there is a fraction of energy absorbed. After many such absorption, practically all the
incident radiation is absorbed.
It is the total radiant energy emitted by the surface in all directions over the entire wave length per
unit surface area per unit time (W/m2). The total emissive power depends on the surface temperature.
Monochromatic (Spectral) Emissive Power (Eλb)
The energy emitted by a surface in all directions at a given wave length is called monochromatic
Emissive power.
Emissivity
It is the ratio of emissive power of a surface to emissive power of a black body at the same
temperature. It is denoted by Ԑ.
𝐄
Ԑ=
𝐄𝐛
Emissivity is the ability of a surface to radiate heat. It depends on the nature of the surface and
independent wave length.
1. Stefan-Boltzmann law
2. Plank’s law
3. Wien’s displacement law
4. Kirchhoff’s law
Stefan-Boltzmann law
This law states that the emissive power of a black body is directly proportional to the fourth power of
its absolute temperature.
Eb = 𝝈 T4
= 2898 𝜇𝑚 − 𝑘
Where, λmax is the wave length at which maximum value of monochromatic power occurs.
Plank’s law
This law states that the total monochromatic emissive power (energy) emitted by the black body is in
all directions per unit area per unit wavelength interval around λ.
C1 λ−5
Ebλ = 𝐶
𝐸𝑥𝑝( 2 )−1
λT
𝜆= Wavelength (𝜇𝑚)
T= Absolute temperature in K
Kirchhoff’s law
It states that at any temperature, the ratio of total emissive power (E) to the total Absorptivity ( α ) is a
constant for all the substances which are in thermal equilibrium with each other.
E
= Constant
α
Consider a large radiating body which encloses a small body (1) of surface area A1. Let Eb be the
energy falling on the surface of body (1). Of this energy, a fraction 𝛂 will be absorbed by the small
body.
Thus at equilibrium, A1 E1 = 𝛼1 A1 Eb
E1
∴ Eb =
α1
E2
Replacing the small real body (1) with another small body (2), having Absorptivity α2, then Eb =
α2
E1 E2 E3 E
Eb = = = -------= α eq(1)
α1 α2 α3
𝐄 E
Also, by definition of emissivity Ԑ = i.e Eb = eq (2)
𝐄𝐛 Ԑ
E E
=
Ԑ α
∴Ԑ= α
Kirchhoff’s law also states that under thermal equilibrium, emissivity of a body is equal to
Absorptivity.
Where, A=Area of the surface of the sphere enclosed by the conical surface forming the angle
R= Radius of sphere
Intensity of radiation: It is the energy emitted by a surface of unit area at all wave lengths in a
particular direction through a unit solid angle.
𝐸𝑏 𝜎𝑇 4
I= =
𝜋 𝜋
𝜃 = Zenith angle ∅= Azimuth angle In= Normal intensity Iθ= Intensity at an angle θ
This law states that the rate of energy radiated by a diffuse surface in any particular direction is
proportional to the cosine of the angle between directions under considerations and normal to the
surface.
∴Iθ= In Cos θ
The energy radiated out decreases with increase in θ and becomes zero at 90°
Example:
The body having area of 1000 cm2 has effective temperature of 900 K. Find
3 According to Lamberts’ cosine law, I𝜃 = In Cos θ = 11841.4 Cos 60° = 592.07 W/m2- Steradian
Example:
An industrial furnace in the form of a black body emits radiation at 3000 K. calculate the following:
1. Monochromatic emissive power at 1 𝜇𝑚 wave length.
2. Maximum emissive power
3. Total emissive power
4. Compare the total emissive power of the furnace, if it is assumed as a real surface having
emissivity 0.8
5. Wave length at which emission is maximum
Solution:
𝐶1 𝜆− 5
1. From Planks law, monochromatic emissive power (Eλ) b = 2 −1]𝐶
𝐸𝑥𝑝[𝜆𝑇
2 Maximum emissive power, (Eλb) max = 1.287 *10-5* T5 = 3.127 * 1012 W/m2
3 Total Emissive Power = 𝜎 T4 = 5.67*10- 8 * 30004 = 4.59 *106 W/m2
4 Total emissive power of a real surface E= Ԑ 𝜎 T4 = 0.8 * 67*10- 8 * 30004 = 3.67 *106 W/m2
5 According to Wien’s displacement law, 𝜆maxT= 2.898 *10-3 = 𝜆max 3000= 2.898 *10-3
Example:
Determine the radiant energy emission between wave lengths 2 to 5 𝜇𝑚 at 1400 k by a black body of
area 1m2.
λ1T = 2800 𝜇𝑚 − 𝐾
λ2T = 7000 𝜇𝑚 − 𝐾
From DHB P-83, the black body radiation functions at these values are
Eb0−2 Eb0−5
4 = 0.2278 and = 0.8081
𝜎T 𝜎 T4
Eb(5−2)
∴ = 0.8081- 0.2278 = 0.5802
𝜎 T4
Example: The monochromatic emissivity [spectral hemispherical emissivity] of fire brick at 750 K as
a function of wave length is as follows:
Ԑ=0.1 for 0 ≤ 𝜆 ≤ 2 𝜇𝑚
Ԑ=0.6 for 2 ≤ 𝜆 ≤ 14 𝜇𝑚
Ԑ=0.3 for 𝜆 ≥ 14 𝜇𝑚
Calculate the hemispherical emissivity over entire wave lengths and overall radiant energy emission.
λ1T = 1500 𝜇𝑚 − 𝐾
λ2T = 10500 𝜇𝑚 − 𝐾
The black body radiation functions from DHB, at λ1T = 1500, is 0.013754 and at λ2T = 10500 is
0.923710 respectively.
E
= Ԑ = 0.1(0.013754 – 0) + 0.6(0.923710 - 0.013754) + 0.3(1- 0.923710)
𝜎 T4
Ԑ = 0.57
E E
Also, Ԑ = 4 i.e. 0.57 = (5.67 x 10−8 x 7504 )
𝜎T
Overall radiant energy emission, E= 10.223 W/m2
Show that the total emissive power is equal to 𝝅 times the intensity of radiation
𝜃 = Zenith angle ∅= Azimuth angle In= Normal intensity Iθ= Intensity at an angle θ
When the collector is located at an angle θ from the normal and subtending a solid angle dω, the rate
of energy emitted through the surface is given by
(dEb)𝜃 = Iθ dω dA = In Cos θ dω dA
Eb = In dA 𝜋 … (i)
Also the total emissive power of the emitter with area dA and temperature T is given by
Eb = 𝜎 T4 dA …. (ii)
Comparing the equations (i) & (ii)
𝜎 T4 dA = In dA 𝜋
Eb = In 𝝅
Example:
A surface with area 2 cm2 emits radiation as a black body at a temperature 1000 K. Calculate the
π
radiation emitted into a solid angle subtended by 0 ≤ ∅ ≤ 2𝜋 and 0 ≤ θ ≤ .
6
Solution: The radiation energy Eb emitted by an area a through the solid angle dω= Sinθ dθd∅ is given
by Eb = In Sinθ Cos θdθ d∅ A
π
The energy into solid angle subtended by the angles ∅ = 0 to 2𝜋 and θ = 0 to is
6
𝜋
2𝜋
Eb= In A ∫06 Sinθ Cos θdθ ∫0 d∅
𝜋
1 2𝜋
= In A ∫06 2 Sin2θ dθ ∫0 d∅
𝜋
6
1 −Cos 2θ
= 𝝈 T4 A 2 [ 2
]0 ∅]2𝜋
0
1
= (5.67 x 10-8 x10004) (2 x 10-4)(− 4 ){𝐶𝑜𝑠 60 − 𝐶𝑜𝑠0} 2𝜋 Cos 60°= 0.5
1
= (5.67 x 10-8 x10004) (2 x 10-4)(− 4 ){−0.5} 2𝜋
= 8.9 W
1. The views of the surfaces have of each other, i.e how they see each other
2. The emitting and absorbing characteristics
Radiation shape factor is defined as “the fraction of the radiative energy that is diffused from one
surface and strikes directly at the other surface without intermediate reflection”.
1. By Summation rule: If N surface, makeup an enclosure then all the radiation streaming
from one surface has to be received by the remaining(n-1) surfaces and by the surface
itself.
Normally, F1-1 = 0
2. By reciprocity theorem: the shape factors F1-2 and F2-1 between two surfaces 1 and 2 are
related by
F1-2 A1= F2-1 A2
3. By Geometry
The shape factor of convex surface or Flat surface with the other surface enclosing the
first is always unity. This is because all the radiation coming out from the convex surface
is received by the enclosing surface.
F1-2 =1 and F1-1 = 0
The concave surface has a shape factor with itself because the radiant Energy coming out
from one part of the surface is received by the part of the same surface.
F1-1 ≠ 0
Example: A bulb hanging at the centre of the roof of the room is measuring (6x4x3) meters. The bulb
is assumed to be spheres of diameter 50 mm. determine the shape factor of the room with respect to
bulb.
Solution:
Body 2: Room,
Area A2= (6x4x2) + (4x3x2) + (6x3x2) = 108 m2
∴ F21=2.9 x 10-4
Example: Determine the shape factor from the floor of the furnace of 1 m x 2 m x 1.5 m size to the
side surfaces and the roof.
Solution:
Solution:
∴ F1-2+ F1-3 = 1
(i) Shape factor from the base of a cylinder to the curved face (F1-3)
d/L= 1
From DHB (P91), at d/L=1 and referring curve 1, shape factor F1-3 = 0.18
F1-2+ 0.18 = 1
F1-2 = 0.82
Example: Determine the shape factor F1-2 and F2-1 for the following cases
Solution:
Case (1): Surface 1 is the base of hemisphere. All the radiation from the surface1, reaches surface 2
∴ F1-2=1
Using Reciprocal theorem, F1-2 A1 = F2-1 A2
𝑑2 𝑑2
1x𝜋 = F2-1x 2 𝜋 r2 Area of circular plate, A1= 𝜋
4 4
∴ F1-2 =1
Using Reciprocal theorem, F1-2 A1 = F2-1 A2
F2-1= 0.424
Considering the surface 2,
F2-1 + F2-2 = 1
0.424 + F2-2 =1
∴ F2-2 =0.575
Assignment: Considering the following long grooves of width 2R determine the shape factors
from the grooves to the outside and also groove itself.
Radiation heat exchange between parallel plates
Consider two surfaces of area A1and A2 (A1=A2) maintained at temperature T1 and T2 as shown in
fig.
Net Radiation heat exchange between two long parallel plates (Real surfaces)
Consider two long parallel plates of surface area A1 and A2 as shown in figure. The surface
emissivities are Ԑ1 and Ԑ2 respectively. Since the plates are very long and parallel to each other and is
equal area (i.e A1 = A2 =A), the amount of heat leaving one surface is completely reaching the other
surface, the shape factor F12 = 1.
Radiation energy emitted by plate 1= E1= Ԑ1σ T14
Similarly energy absorbed by the plate 2 on third incidence = E1 𝜶𝟐 [1- 𝜶𝟐 ]2 [𝟏 − 𝜶𝟏 ]2 ……. (3)
1
Q1= E1𝛼2 [1−[1−𝛼 ]
1 ][1−𝛼2 ]
𝛼2
∴ Q1= E1 [ ] By Kirchhoff’s law, Absorptivity = emissivity i.e 𝛼 = Ԑ
1−[1−𝛼1 ][1−𝛼2 ]
Ԑ
= E [1−[1−Ԑ 2][1−Ԑ ]]
1
1 2
Ԑ2 Ԑ2
= E1 [ ] = E1 [ ]
1−1+Ԑ1 +Ԑ2 −Ԑ1 Ԑ2 ] Ԑ1 +Ԑ2 −Ԑ1 Ԑ2 ]
Ԑ1
Similarly, Q2 = E2 [ ]
Ԑ1 +Ԑ2 −Ԑ1 Ԑ2 ]
E1 Ԑ2 − E2 Ԑ1
∴ The net heat flow from plate 1 to plate 2 is, Q12 = Q1- Q2 =
Ԑ1 +Ԑ2 −Ԑ1 Ԑ2
1
Q12 = 1 1 σ [ T14 - T24 ]
+ −1
Ԑ1 Ԑ2
1
Q12 = Ԑ̅ σ A (𝑻𝟒𝟏 - 𝑻𝟒𝟐 ) where, Ԑ̅ = equivalent emissivity = 1 1
+ −1
Ԑ1 Ԑ2
Net Radiation heat exchange between two long concentric cylinders (Real
surfaces)
Consider two long concentric cylinders as shown in fig.
Let,
Since the two cylinders are long and concentric, the shape factor F12 is unity and the net heat
exchange between them is given by
1
Q12 = Ԑ̅ σA (𝑻𝟒𝟏 - 𝑻𝟒𝟐 ) where, Ԑ̅ = equivalent emissivity = 1 A1 1
+ ( − 1)
Ԑ1 A Ԑ2 2
Example:
Determine the net radiation heat exchange between two parallel black rectangular plates (60x80) cm,
separated by a distance of 40 cm, if their surface temperatures are 500℃ and 200℃.
Example: Two black discs of diameter 60 cm are arranged directly opposite to each other and
separated by a distance of 120 cm. The temperatures of the discs are 1200 K and 600 k. Estimate the
heat flow by radiation between the discs for the following cases.
1. When no other surfaces are present
2. When the discs are connected by a black right cylindrical non flux surface. (Refractory).
𝜋 𝑑2
Area A= = 0.2827 m2
4
Ratio d/L= 0.6/1.2= 0.5
From DHB (P 91) referring curve 1 for discs not connected by other walls, F12= 0.07
Ԑ1 = Ԑ2 =1 (Black body)
Ԑ̅= Ԑ1 Ԑ2=1
Heat exchange between two bodies is given by, Q12 = Ԑ̅ σ F12 A [T14 - T24 ]
Q12 = 1 x (5.67 x 10-8) x 0.06 x 0.2827 x [12004 – 6004]
Q12 =1879 W
2 When the discs are connected by a black right cylindrical non flux surface. (Refractory)
From DHB (P 91) referring curve 1 for discs not connected by other walls, F12=0.34
Heat exchange between two bodies is given by, Q12 = Ԑ̅ σ F12 A [T14 - T24 ]
Q12 = 1 x (5.67 x 10-8) x 0.34 x 0.2827 x [12004 – 6004]
Q12 =10650 W
Example: Two closely spaced parallel plates radiating only to each other are maintained at a
temperature of 900 K and 450 K respectively. Each plate has an area of 5 m 2. Compare the net heat
exchange between the plates for the following cases:
Ԑ1 = Ԑ2 =1 (Black body)
Ԑ̅= 1
Heat exchange between two bodies is given by, Q12 = Ԑ̅ σ F12 A [T14 - T24 ]
Q12 = 1 x (5.67 x 10-8) x 1 x 5 x [9004 – 4504]
Q12 =174379 W
Ԑ1 = Ԑ2 =0.4
1 1
Equivalent emissivity, Ԑ̅ = 1 1 = 1 1 = 0.25
+ −1
Ԑ1 Ԑ2
+ −1
0.4 0.4
Heat exchange between two bodies is given by, Q12 = Ԑ̅ σ F12 A [T14 - T24 ]
Q12 = 1 x (5.67 x 10-8) x 0.25 x 5 x [9004 – 4504]
Q12 =43594.75 W
Example: Two concentric spheres 25 cm and 35 cm in diameter with the space between them
evacuated are used to store liquid air at -150℃ in a room at 20℃. The surfaces of the sphere are
flushed with aluminum (Ԑ=0.04). Calculate the rate of evaporation of liquid air if the latent heat of
vaporization of liquid air is 218 KJ/kg.
Solution: d1= 25x 10-2 m, d2= 35 x 10-2 m, T1= -150+273= 123 K, T2= 20+273= 293 K
1 1
Equivalent emissivity Ԑ̅ = 1 A1 1 = 1 2
4 𝜋 𝑟1 1
= 0.0268
+ ( − 1)
Ԑ1 A2 Ԑ2 + 2 ( − 1)
0.04 4 𝜋 𝑟2 0.04
Also, Q= mf x hfg
2.13 = mf x 218
Rate of evaporation, mf =9.7 x10-6 kg/sec
Radiation shields
In many engineering problems it is necessary to reduce the radiant heat transfer between surfaces.
One method is use materials which are highly reflective. Another method is to use radiation shields,
between the heat exchange surfaces. These shields do not deliver or remove any heat from the overall
system; but they only place another resistance in the heat flow path, resulting in decreasing in overall
heat transfer.
Example: Two black discs of diameter 60 cm are arranged directly opposite to each other and
separated by a distance of 120 cm. The temperatures of the discs are 1200 K and 600 k. Estimate the
heat flow by radiation between the discs for the following cases.
1. When no other surfaces are present
2. When the discs are connected by a black right cylindrical non flux surface. (Refractory).
Heat transfer from plate 1 to plate 2, without radiation shield (shown in fig.1.)
𝟏
Q12 = 𝟏 𝟏 σ A (𝑻𝟒𝟏 -𝑻𝟒𝟐 ) …….(1)
+
Ԑ𝟏 Ԑ𝟐
− 𝟏
When the radiation shield is placed (refer fig.2), the transfer
𝟏
Q13 = 𝟏 𝟏 σ A (𝑻𝟒𝟏 - 𝑻𝟒𝟑 ) ……..(2)
+
Ԑ𝟏 Ԑ𝟑
− 𝟏
𝟏
Heat transfer from plate 3 to plate 2, Q32 = 𝟏 𝟏 σ A (𝑻𝟒𝟑 - 𝑻𝟒𝟐 ) ……..(3)
+
Ԑ𝟑 Ԑ𝟐
− 𝟏
Since, the sheet does not deliver or remove the heat from the system, Q13 = Q32
𝟏 𝟏
𝟏 𝟏 σ A (𝑻𝟒𝟏 - 𝑻𝟒𝟑 ) = 𝟏 𝟏 σ A (𝑻𝟒𝟑 - 𝑻𝟒𝟐 )
+ −𝟏
Ԑ𝟏 Ԑ𝟑
+
Ԑ𝟑 Ԑ𝟐
− 𝟏
1 T41 + T42
Q13 = 1 1
+ −1
[
σ A 𝑇14 - (
2
) ]
Ԑ1 Ԑ3
1 T41 − T42
Q13 = 1 1
+ −1
σA[
2
] but, Ԑ1 = Ԑ3 = Ԑ
Ԑ1 Ԑ3
4 4
σ A [T1 − T2 ]
Q13 = 2 .......(4) and
2 ( Ԑ −1 )
4 4
σ A [T1 − T2 ]
Equation (1) becomes, Q12 = 2 .......(4)
( Ԑ −1 )
Dividing Q13 by Q12
𝑄1−3 1
=
𝑄1−2 2
𝟏
Heat transfer with shield =
𝟐 heat transfer without shield
1
If there are N shield, QN = Q0
N+1
𝟏
Q with N shield = Q without shield
𝐍+𝟏
Example: Two parallel plates at T1=900 K and T2= 500 K, have emissivities Ԑ1=0.6 and Ԑ2=0.9
respectively. A radiation shield having an emissivity Ԑ31=0.15 on one side and emissivity Ԑ32=0.06 on
the other side is placed between the plates. Calculate the heat transfer by radiation per sq meter with
and without radiation shield.
Solution: T1=900 K, T2= 500 K, Ԑ1=0.6, Ԑ2=0.9, Ԑ13=0.15, Ԑ32=0.06, Q12 w/o shield =?
Q12 with shield =?
𝟏
Q12 w/o shield= 𝟏 𝟏 σ A (𝑻𝟒𝟏 -𝑻𝟒𝟐 )
+ −𝟏
Ԑ𝟏 Ԑ𝟐
1
= 1 1 (5.67 x 10-8) 1(9004- 5004)
+ −1
0.6 0.9
σ A (𝑇41 −𝑇42 )
Heat transfer with shield, Q12 with shield= 1 1 1 1
[ + − 1]+ [Ԑ32 + Ԑ2 − 1]
Ԑ1 Ԑ13
(5.67 x 10−8 ) 1[ 9004 − 5004 ]
= 1 1 1 1
[ + − 1]+ [0.06 + 0.9 − 1]
0.6 0.15
Example: Two large brass plates are maintained parallel to each other and one at 4000C and other at
400C. Find the % reduction in the heat transfer when polished aluminum sheet is placed between
them. Also find the temperature of aluminum sheet.
For brass, Ԑ=0.1 For aluminum, Ԑ=0.06 (facing higher temperature)
Ԑ= 0.03 (facing lower temperature)
1
Q12 w/o shield= 1 1 σ A (T14 -T24 )
+ −1
Ԑ1 Ԑ2
1
= 1 1 (5.67 x 10-8) 1(6734- 3134)
+ −1
0.1 0.1
σ A (𝑇41 −𝑇42 )
Heat transfer with shield, Q12 with shield= 1 1 1 1
[ + − 1]+ [Ԑ32 + Ԑ2 − 1]
Ԑ1 Ԑ13
(5.67 x 10−8 ) 1[ 6734 − 3134 ]
= 1 1 1 1
[ + − 1]+ [0.03 + 0.1 − 1]
0.1 0.06
T3=601.9 K
Example: Two large parallel plates with Ԑ=0.5 each, maintained at two different temperatures are
exchanging heat only by radiation. Two equally large radiation shields with surface emissivity Ԑ=0.05
are introduced in parallel to the plates. Find the percentage heat reduction in net radiation heat
exchange.
Solution:
Heat transfer without shield
σ A (𝑇41 −𝑇42 )
Q12 w/o shield= 1 1 let, C= σ A (𝑇14 − 𝑇24 )
[ + − 1]
Ԑ1 Ԑ2
C
= 1 1
[
0.5
+ 0.5− 1]
Q12 w/o shield = 0.33 C
σ A (𝑇41 −𝑇42 ) C
Q12 with shield= 1 1 1 1 1 1 = 1 1 1 1 1 1
[ + − 1]+ [Ԑ3 + Ԑ4 − 1]+[Ԑ4 + Ԑ2 − 1] [ + − 1]+ [0.05 + 0.05− 1]+[0.05 + 0.5 − 1]
Ԑ1 Ԑ3 0.5 0.05
C
= [ ]+ [ ]+[ ]
21 39 21