Thermal RadiationIntroduction and contextThermal radiation is energy transported by electromagnetic waves or photons as a consequence of matter being at a temperature above absolute zero. Unlike conduction and convection, radiation requires neither a material medium nor bulk fluid motion. It is therefore the dominant heat-transfer mode in furnaces, combustion chambers, re-entry vehicles, solar-thermal systems, infrared instruments, and cryogenic environments. Radiation is also central to thermodynamics. A surface emits, absorbs, reflects, and transmits radiation according to its temperature and spectral characteristics. Engineering analysis must consequently combine energy conservation with geometric relationships, wavelength dependence, surface properties, and the laws of blackbody radiation. Learning objectivesAfter studying this chapter, the reader should be able to: - explain the physical origin and distinctive features of thermal radiation;
- apply Planck’s, Wien’s, and Stefan–Boltzmann laws;
- distinguish black, gray, diffuse, specular, opaque, and transparent surfaces;
- calculate radiation exchange between surfaces using view factors;
- analyze radiation shields, enclosures, and participating gases;
- formulate engineering solutions using radiosity and irradiation;
- assess assumptions, limitations, uncertainty, and practical design implications.
1. Physical basis of thermal radiationMatter contains charged particles in continual thermal motion. Accelerated charges emit electromagnetic radiation. Thermal radiation occupies a broad range of wavelengths, although engineering heat transfer is commonly concentrated in the infrared region. The wavelength λ and frequency ν are related by c = λν where c is the speed of light in vacuum, approximately 2.998 × 10⁸ m/s. A photon has energy Eph = hν = hc/λ where h is Planck’s constant, 6.626 × 10⁻³⁴ J·s. Radiation is emitted over a spectrum; it is not generally confined to one wavelength. The net radiative heat transfer between bodies depends strongly on the fourth power of absolute temperature. This makes radiation particularly important at high temperature. In contrast with conduction and convection, radiative transfer can occur through a vacuum, as demonstrated by energy received from the Sun. Radiation incident on a surface is divided into three fractions: α + ρ + τ = 1 Here α is absorptivity, ρ is reflectivity, and τ is transmissivity. For an opaque surface, τ = 0 and therefore α + ρ = 1. These quantities may depend on wavelength and direction. 2. Blackbody and spectral emissionA blackbody is an ideal surface that absorbs all incident radiation at every wavelength and direction. Thus α = 1, ρ = 0, and τ = 0 for an opaque blackbody. It is also the maximum possible emitter at a specified temperature. Planck’s lawThe spectral emissive power of a blackbody is Eλ,b(λ,T) = [C₁/λ⁵] / [exp(C₂/(λT)) − 1] where E_{λ,b} is the emitted power per unit area per unit wavelength, T is absolute temperature in kelvin, λ is wavelength, C₁ = 2πhc², and C₂ = hc/k, with k being Boltzmann’s constant. The wavelength must be expressed consistently with the chosen units of the constants. Planck’s law explains both the magnitude and spectral distribution of blackbody emission. Integration over all wavelengths gives the Stefan–Boltzmann law: Eb = σT⁴ where E_b is total blackbody emissive power in W/m² and σ = 5.670 × 10⁻⁸ W/(m²·K⁴). This law applies to total emission from an ideal blackbody. Wien’s displacement lawThe wavelength at which blackbody emission is maximum is λmax T = b where b = 2.898 × 10⁻³ m·K. Increasing temperature shifts the peak toward shorter wavelengths. For a surface at 1000 K, λ_max ≈ 2.90 μm; for the Sun at approximately 5800 K, the peak is near 0.50 μm. The fraction of blackbody radiation emitted between selected wavelengths is obtained from tabulated or computational blackbody radiation functions. Such fractions are essential when filters, selective coatings, or spectral gases are present. 3. Emissivity and real surfacesThe emissivity of a surface is the ratio of its emission to blackbody emission at the same temperature: ε = E/Eb For total hemispherical emission, 0 ≤ ε ≤ 1. A real surface may have spectral emissivity ε_λ and directional emissivity ε(θ, φ). A polished metal often has low emissivity, whereas an oxidized or non-metallic surface commonly has higher emissivity. A gray surface is an idealization whose emissivity is independent of wavelength over the relevant spectrum. A diffuse surface has direction-independent intensity. Real engineering surfaces are neither perfectly gray nor perfectly diffuse, but these assumptions often provide useful first-order models. Kirchhoff’s lawAt thermal equilibrium, for a specified wavelength and direction, ε_λ(θ,φ) = α_λ(θ,φ) For a gray, opaque surface this is often written ε = α. Kirchhoff’s law does not state that emissivity equals absorptivity under arbitrary non-equilibrium conditions unless the same spectral and directional conditions are being compared. This distinction matters for selective surfaces and laser or narrow-band radiation. For a diffuse-gray opaque surface, the emitted power is E = εσT⁴ The net exchange with a large black surrounding at temperature T_sur is q̇ = εσA(Ts⁴ − Tsur⁴) where q̇ is net heat transfer from the surface, A is area, and T_s and T_sur are absolute temperatures. The surroundings must be sufficiently large that the surface views them completely and their radiation is approximated as blackbody radiation. 4. Radiation intensity and geometryRadiation intensity I is the rate of energy leaving a surface per unit projected area per unit solid angle. For a diffuse emitter, intensity is independent of direction. The relation between emissive power and intensity is E = πI The cosine dependence of projected area is called Lambert’s cosine law. It explains why a diffuse surface appears equally bright from different directions when brightness is measured appropriately, even though its emitted power per unit solid angle varies with viewing angle. View factorsA view factor F_{i→j} is the fraction of radiation leaving surface i that directly reaches surface j. For surfaces A_i and A_j, Fi→j = (radiation from i reaching j)/(total radiation leaving i) For an enclosure, the summation rule is ∑j Fi→j = 1 and the reciprocity relation is AiFi→j = AjFj→i These relations follow from conservation of energy and geometric reciprocity. For a convex surface, F_{i→i} = 0. View factors depend only on geometry, not on temperature or material properties. Standard charts, analytical integrations, and numerical ray-tracing methods are used to determine them. 5. Exchange between black surfacesFor two black surfaces exchanging radiation directly, the net transfer from surface 1 to surface 2 is q̇1→2 = A₁F1→2σ(T₁⁴ − T₂⁴) This equation assumes diffuse emission, black surfaces, and no participating medium or intervening radiation source. It is especially useful for large parallel plates, concentric cylinders, and enclosure components. For two large, parallel black plates, F_{1→2} = 1, giving q̇ = Aσ(T₁⁴ − T₂⁴) The fourth-power dependence demonstrates why replacing temperature in kelvin by degrees Celsius produces incorrect results. 6. Radiosity method for diffuse-gray enclosuresFor a real diffuse-gray opaque surface, the irradiation G is the incident radiation per unit area and the radiosity J is the total radiation leaving the surface. Radiosity consists of emitted and reflected components: Ji = εiσTi⁴ + ρiGi For an opaque gray surface, ρ_i = 1 − ε_i, so Ji = εiσTi⁴ + (1 − εi)Gi The net heat leaving surface i is q̇i = Ai(Ji − Gi) For an enclosure, Gi = ∑j Fi→jJj These equations form a simultaneous linear system in the radiosities. They are applicable to diffuse-gray, opaque surfaces with a nonparticipating medium. Once radiosities are known, surface heat rates follow directly. An equivalent resistance network is often useful. For two diffuse-gray surfaces with direct exchange, q̇ = [σ(T₁⁴ − T₂⁴)] / [ (1−ε₁)/(A₁ε₁) + 1/(A₁F1→2) + (1−ε₂)/(A₂ε₂) ] The denominator contains surface resistances and a space resistance. The expression is valid when the two-surface enclosure representation is appropriate and areas, view factors, and temperatures are consistently defined. Worked example: two large platesTwo large parallel plates have T₁ = 800 K, T₂ = 500 K, ε₁ = 0.80, ε₂ = 0.60, and A = 2 m². For large parallel plates, F_{1→2} = 1. Therefore, q̇ = Aσ(T₁⁴ − T₂⁴) / [1/ε₁ + 1/ε₂ − 1] The temperature term is T₁⁴ − T₂⁴ = 800⁴ − 500⁴ = 3.59 × 10¹¹ K⁴. The denominator is 1/0.80 + 1/0.60 − 1 = 2.9167. Hence, q̇ ≈ (2)(5.670 × 10⁻⁸)(3.59 × 10¹¹)/2.9167 ≈ 1.40 × 10⁴ W. Thus approximately 14.0 kW is transferred from the hotter plate to the colder plate. Convection and conduction through the intervening space have been neglected. 7. Radiation shieldsA radiation shield is a low-emissivity surface inserted between hot and cold surfaces. It reduces exchange by introducing additional surface and space resistances. For two large parallel gray plates with N identical, negligible-thickness shields, all having emissivity ε_s, the total resistance is found by adding the resistances of the surfaces and gaps. For a particularly useful special case in which both original plates and all shields have the same emissivity ε, q̇ = Aσ(T₁⁴ − T₂⁴) / [(N + 1)(2/ε − 1)] This result assumes large parallel surfaces, equal emissivities, diffuse-gray behavior, vacuum or nonparticipating gaps, and negligible conduction through the shields. In practice, shield supports can create parasitic conduction paths, and shield temperature distributions must be determined if shield temperatures are not prescribed. 8. Radiation in gases and participating mediaA nonparticipating medium neither absorbs nor emits appreciable radiation; dry air in many short-path laboratory problems is approximately nonparticipating. Combustion gases containing CO₂ and H₂O can be strongly participating, especially over long optical paths and at high temperatures. For a participating medium, radiative transfer depends on absorption coefficient, emission, scattering, composition, pressure, temperature, and path length. The radiative transfer equation (RTE) describes the change of intensity along a path. In a non-scattering medium, a one-dimensional form is dI_λ/ds = −κ_λ I_λ + κ_λ Iλ,b where I_λ is spectral intensity, s is path length, κ_λ is the spectral absorption coefficient, and I_{λ,b} is blackbody spectral intensity. The first term represents attenuation; the second represents thermal emission. Spectral-line and band models are used because gas properties vary strongly with wavelength. Engineering methods include mean-beam-length correlations, weighted-sum-of-gray-gases models, narrow-band models, and computational solutions of the RTE. Soot in flames can make radiation approximately gray over selected ranges but also introduces strong uncertainties in emissivity and concentration. 9. Radiation combined with convection and conductionReal systems commonly involve all three heat-transfer modes. At a surface, a steady energy balance may be written q̇in − q̇out = 0 with convective and radiative terms combined, for example, q̇ = hA(T_∞ − Ts) + εσA(Tsur⁴ − Ts⁴) Here h is the convection coefficient and T_∞ is the fluid temperature. The equation is nonlinear in T_s because of the fourth-power radiation term. It should be solved iteratively unless a linearized radiation coefficient is acceptable. A useful linearization is hr = εσ(Ts + Tsur)(Ts² + Tsur²) so that q̇rad = hrA(Tsur − Ts) This is an exact algebraic rearrangement for two specified temperatures, not a constant material property. The value of h_r changes with temperature and emissivity. 10. Applied case: thermal design of a vacuum furnaceBackgroundAn illustrative industrial case concerns a vacuum furnace used to heat metal components to 900 K. Because the chamber pressure is low, convection is negligible. The components are surrounded by stainless-steel walls near 450 K. ProblemThe initial design uses a single wall with emissivity 0.70. The components have area 1.5 m² and emissivity 0.80. Excessive radiative loss increases heater power and lengthens the heating cycle. Analysis and approachThe chamber is approximated as two large facing gray surfaces. With component temperature T₁ = 900 K and wall temperature T₂ = 450 K, q̇ = Aσ(T₁⁴ − T₂⁴) / (1/ε₁ + 1/ε₂ − 1) The design team evaluates a polished metallic shield with ε_s = 0.10 placed between the component and wall. The shield is supported by ceramic spacers, and the support conduction is estimated separately. The radiative model predicts a major reduction in exchange because the shield introduces two low-emissivity surface resistances. A thermal balance on the shield is then required; its temperature cannot simply be assumed to equal either adjacent surface. Outcome and engineering judgementThe shield reduces heater demand, but the final design must include conduction through supports, radiation from openings, nonuniform temperatures, and the effect of oxidation. A coating that initially has ε = 0.10 may become substantially more emissive after repeated thermal cycling. Instrumentation should therefore measure heater power, wall temperature, component temperature, and shield temperature during qualification tests. Learning points- Vacuum eliminates convection but not radiation.
- Surface emissivity can be more influential than material bulk conductivity.
- Radiation shields require a complete resistance and support-conduction assessment.
- Long-term surface degradation can invalidate an idealized design.
- Model validation should include energy balances and calibrated temperature measurements.
11. Practical design, limitations, and contemporary methodsImportant sources of error include using Celsius in fourth-power expressions, neglecting view factors, treating a selective surface as gray, ignoring gas radiation, assuming a shield has a prescribed temperature, and confusing emissivity with absorptivity under non-equilibrium spectral conditions. Good practice includes: - use kelvin for all absolute-temperature radiation laws;
- document whether properties are spectral, directional, or total hemispherical;
- calculate view factors independently and check summation and reciprocity;
- perform sensitivity studies for emissivity, temperature, and gas composition;
- include conduction through supports and convection through residual gas where relevant;
- validate models with calorimetry, infrared thermography, or calibrated heat-flux sensors;
- report uncertainty in surface temperature and emissivity, not only numerical precision.
Contemporary thermal systems increasingly use selective solar absorbers, multilayer insulation, metamaterial-inspired coatings, thermal camouflage, radiative cooling surfaces, and coupled CFD–radiation solvers. Monte Carlo ray tracing is effective for complex geometries and participating media, while discrete ordinates, finite-volume, and spherical-harmonics methods solve approximate forms of the RTE. Machine-learning surrogates can accelerate repeated radiation calculations, but they must be trained over physically valid ranges and checked against conservation laws. AI should assist—not replace—verification of energy balance, boundary conditions, and property data. Radiation design also has safety implications. High-temperature furnaces can produce intense infrared exposure; laser-like sources require wavelength-specific controls; and reflective surfaces may redirect hazardous energy. Relevant occupational, pressure-vessel, optical-radiation, and fire-safety regulations must be incorporated into professional design. Chapter SummaryThermal radiation is electromagnetic energy emitted by matter because of its temperature and can transfer heat through a vacuum. Blackbody behavior is described by Planck’s law, while integration gives the Stefan–Boltzmann law and the emission peak is located using Wien’s law. Real surfaces are characterized by absorptivity, reflectivity, transmissivity, and emissivity; Kirchhoff’s law connects absorptivity and emissivity under identical spectral and directional equilibrium conditions. View factors describe geometric coupling between surfaces and obey summation and reciprocity relations. For real diffuse-gray enclosures, the radiosity method combines emission and reflection into solvable energy balances. Radiation shields reduce heat transfer by adding surface and space resistances. Participating gases require spectral radiative-transfer models rather than simple surface-to-surface formulas. Practical design requires coupling radiation with convection and conduction, checking units and assumptions, accounting for degradation and parasitic paths, and validating models experimentally. Key TermsBlackbody: Ideal surface that absorbs all incident radiation and emits the maximum possible radiation at a given temperature. Emissive power: Radiation emitted per unit surface area. Emissivity: Ratio of real-surface emission to blackbody emission at the same temperature. Absorptivity: Fraction of incident radiation absorbed by a surface. Reflectivity: Fraction of incident radiation reflected by a surface. Transmissivity: Fraction of incident radiation transmitted through a surface. Planck’s law: Spectral distribution law for blackbody emission. Stefan–Boltzmann law: Total blackbody emission proportional to absolute temperature to the fourth power. Wien’s displacement law: Relation between emission-peak wavelength and absolute temperature. Gray surface: Ideal surface with wavelength-independent properties over the relevant range. Diffuse surface: Surface whose radiation is direction-independent in intensity. Irradiation: Radiation incident on a surface per unit area. Radiosity: Total radiation leaving a surface per unit area. View factor: Fraction of radiation leaving one surface that directly reaches another. Kirchhoff’s law: Equality of spectral directional emissivity and absorptivity at thermal equilibrium. Radiation shield: Low-emissivity intervening surface that reduces radiative exchange. Participating medium: Medium that absorbs, emits, or scatters radiation. Radiative transfer equation: Differential equation governing intensity through an emitting and absorbing medium. Radiation resistance: Equivalent resistance representing surface or geometric opposition to radiative exchange. Review QuestionsShort Answer Questions- Why can thermal radiation transfer energy through a vacuum?
- Define blackbody emissive power and state the Stefan–Boltzmann law.
- Distinguish irradiation from radiosity.
- State the view-factor reciprocity relation.
- Why is emissivity important in radiation-shield design?
Long Answer Questions- Explain Planck’s law, Wien’s displacement law, and the Stefan–Boltzmann law, including their engineering significance.
- Derive the radiosity relation for an opaque diffuse-gray surface and explain each term.
- Discuss the assumptions and limitations of the two-surface radiation-resistance method.
- Explain how participating gases alter radiation analysis in combustion systems.
- Describe a systematic procedure for designing and validating a vacuum-furnace radiation model.
Analytical / Application Questions- Two large parallel plates at 750 K and 400 K have emissivities 0.75 and 0.50. Derive the net radiative heat-transfer expression and identify every resistance.
- A radiation shield is installed in a cryogenic enclosure. Determine which additional heat-transfer paths must be evaluated before claiming the predicted reduction in heat leak.
- A furnace model predicts a much smaller heat loss than experiment. Develop a diagnostic sequence using view-factor checks, emissivity measurements, gas participation, support conduction, and energy-balance closure.
Multiple Choice Questions with AnswersQuestion 1Which law gives the total emissive power of a blackbody? A. Kirchhoff’s law B. Stefan–Boltzmann law C. Wien’s displacement law D. Fourier’s law Correct Answer: B. Stefan–Boltzmann law Question 2For an opaque surface, which relation is correct? A. α + ρ = 1 B. α + τ = 1 C. ρ + τ = 1 D. α = ρ = τ Correct Answer: A. α + ρ = 1 Question 3The wavelength of maximum blackbody emission varies with temperature according to A. λ_max/T = constant B. λ_maxT = constant C. λ_max + T = constant D. λ_maxT⁴ = constant Correct Answer: B. λ_maxT = constant Question 4For a convex surface in an enclosure, the self-view factor is A. 1 B. 0.5 C. 0 D. Dependent only on emissivity Correct Answer: C. 0 Question 5Radiosity is best defined as A. Radiation incident on a surface B. Total radiation leaving a surface C. Radiation absorbed by a gas D. Radiation transmitted through a vacuum Correct Answer: B. Total radiation leaving a surface Question 6A radiation shield primarily reduces heat transfer by A. Increasing convection B. Adding radiative resistance through additional low-emissivity surfaces C. Increasing gas pressure D. Increasing the fourth power of temperature Correct Answer: B. Adding radiative resistance through additional low-emissivity surfaces Question 7A surface at 1000 K has its blackbody emission peak approximately at A. 0.29 μm B. 2.90 μm C. 29.0 μm D. 290 μm Correct Answer: B. 2.90 μm Question 8For two large parallel gray plates, the view factor from plate 1 to plate 2 is A. 0 B. 0.25 C. 0.5 D. 1 Correct Answer: D. 1 Question 9Which medium requires a radiative-transfer equation rather than only a surface-exchange model? A. A vacuum B. A nonparticipating transparent gas C. A strongly absorbing and emitting combustion gas D. An ideal empty enclosure Correct Answer: C. A strongly absorbing and emitting combustion gas Question 10In the expression h_r = εσ(T_s + T_sur)(T_s² + T_sur²), temperatures must be expressed in A. degrees Celsius B. degrees Fahrenheit C. kelvin D. any temperature scale Correct Answer: C. kelvin
|