Gray Body Radiation Calculation

Reference ID: MET-972E | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

Gray-body radiation is the dominant mode of heat loss from hot, uninsulated process equipment. Unlike ideal black bodies, whose behavior is described in our black body radiation calculation, real metallic surfaces emit only a fraction of the theoretical maximum; this fraction is the emissivity \( \varepsilon \). Accurate prediction of radiative flux \( E \) is essential for:

Designing fired heaters, reformers and cracking furnaces
Estimating heat losses from reactors, distillation columns and flare stacks
Setting safe touch temperatures on piping and guarding against insulation degradation
Balancing energy budgets in high-temperature drying, calcination and sintering operations

Methodology & Formulas

Convert practical temperature to absolute scale
\[ T[\text{K}] = T[^{\circ}\text{C}] + 273.15 \]
Compute hemispherical emissive power
\[ E = \varepsilon \, \sigma \, T^{4} \] where
\( \sigma = 5.670 \times 10^{-8} \ \text{W m}^{-2}\text{ K}^{-4} \) (Stefan–Boltzmann constant)
\( \varepsilon \) = surface emissivity, dimensionless

Validity Criteria
Parameter	Lower Limit	Upper Limit	Remark
Temperature \( T \)	\( > 0\ \text{K} \)	—	Absolute scale required
Emissivity \( \varepsilon \)	\( 0 \)	\( 1 \)	Gray-body assumption

The resulting flux \( E \) is expressed in W m^-2; divide by 1000 to obtain kW m^-2 for plant-level energy balances.

A black body absorbs all incident radiation and emits the maximum possible thermal radiation at every wavelength.

A gray body absorbs and emits a constant fraction of that ideal amount; this fraction is the emissivity ε (0 < ε < 1).

In practice, most industrial surfaces are treated as gray bodies to simplify heat-transfer models while still giving reliable results.

Use the gray-body radiation equation:

Q = σ · F₁₂ · A₁ · (T₁⁴ – T₂⁴)

where

σ = 5.670 × 10⁻⁸ W m⁻² K⁻⁴ (Stefan–Boltzmann constant)
F₁₂ = view factor times gray-body factor [1/(1/ε₁ + 1/ε₂ – 1)] for two infinite parallel plates
A₁ = area of surface 1 (m²)
T₁, T₂ = absolute temperatures of the two surfaces (K)

Adjust F₁₂ for other geometries using the appropriate view-factor relation.

Standard references: Perry’s Chemical Engineers’ Handbook, ASME Heat Transfer Data Book, and Touloukian & Ho thermophysical property series.
Manufacturer data sheets for refractories, oxidized steels, and specialty alloys.
In-house measurements with IR thermography or emissometers when coatings or surface conditions vary.

Always use values measured at the wavelength and temperature range of your process; emissivity can change with oxidation, roughness, and temperature.

Pure gray-body assumptions are inadequate for CO₂ and H₂O because these gases emit and absorb only in discrete infrared bands.

For furnace design, use the weighted-sum-of-gray-gases model or spectral (band) radiation models to account for the wavelength-dependent absorption.

Process simulators such as Aspen Plus, Fluent, or HSC Chemistry include these models and will give far more accurate heat-flux predictions than a simple ε·σ·T⁴ approach.

Worked Example: Radiation Heat Loss from a Stainless-Steel Reactor Vessel

A process engineer is evaluating the heat loss from a 0.5 m diameter horizontal reactor operating at 350 °C. The vessel’s outer surface is currently polished stainless steel, but after six months of service it will become heavily oxidised. Compare the radiative heat loss per square metre for both surface conditions.

Knowns

Stefan–Boltzmann constant, \(\sigma = 5.670\times10^{-8}\) W m\(^{-2}\) K\(^{-4}\)
Surface temperature, \(T = 350\) °C
Emissivity of polished stainless steel, \(\varepsilon_{\text{polished}} = 0.07\)
Emissivity of oxidised stainless steel, \(\varepsilon_{\text{oxidised}} = 0.79\)
Absolute temperature, \(T = 350 + 273.15 = 623.15\) K

Step-by-Step Calculation

Convert the operating temperature to kelvin: \[ T = 350 + 273.15 = 623.15\ \text{K} \]
Compute the black-body emissive power: \[ E_b = \sigma T^4 = 5.670\times10^{-8} \times (623.15)^4 = 8550.465\ \text{W m}^{-2} \]
Calculate the gray-body emissive power for the polished surface: \[ E_{\text{polished}} = \varepsilon_{\text{polished}} E_b = 0.07 \times 8550.465 = 598.483\ \text{W m}^{-2} \]
Calculate the gray-body emissive power for the oxidised surface: \[ E_{\text{oxidised}} = \varepsilon_{\text{oxidised}} E_b = 0.79 \times 8550.465 = 6754.303\ \text{W m}^{-2} \]

Final Answer

The radiative heat loss per unit area is 598 W m\(^{-2}\) for the polished surface and 6754 W m\(^{-2}\) for the oxidised surface—more than an eleven-fold increase.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle

Gray Body Radiation Calculation

Introduction & Context

🚀 Skip the Manual Math!

Methodology & Formulas

Worked Example: Radiation Heat Loss from a Stainless-Steel Reactor Vessel