Radiation Exchange Between Gray Diffuse Surfaces

Understanding radiation exchange is crucial for designing everything from spacecraft thermal control systems to industrial furnaces. While the Stefan-Boltzmann law tells us how much energy a perfect blackbody emits, real engineering surfaces are not perfect, and they exist within enclosures where they constantly exchange energy with other surfaces. The radiosity method, using an electrical network analogy, provides a powerful, systematic way to solve these complex, multi-surface heat transfer problems, transforming a web of radiative interactions into a solvable circuit.

Foundations: Gray, Diffuse Surfaces and Radiosity

To apply the network method, we must first define the behavior of the surfaces involved. A gray surface is one whose emissivity $ϵ$ is independent of wavelength. This is a practical assumption for many engineering materials over a limited temperature range. A diffuse surface emits and reflects radiation uniformly in all directions; its radiative properties are independent of the direction of incoming or outgoing radiation. Most rough engineering surfaces approximate diffuse behavior reasonably well.

The central quantity in this analysis is radiosity, denoted by $J$. Radiosity is defined as the total rate of radiative energy leaving a surface per unit area. It comprises two components: the energy the surface emits directly, plus the portion of incoming radiation that it reflects. For a gray, diffuse surface of area $A$ at temperature $T$, with emissivity $ϵ$ and reflectivity $\rho$, the radiosity is: $$J=ϵ\sigma T^4+\rho G$$ where $G$ is the irradiation (incoming radiative flux) and $\sigma$ is the Stefan-Boltzmann constant. For an opaque surface ($\tau =0$), the reflectivity is $\rho =1-\alpha$. For a gray surface, Kirchhoff's law tells us that absorptivity equals emissivity ($\alpha =ϵ$), so $\rho =1-ϵ$. This simplifies the radiosity equation to a function of emissivity and temperature.

The Electrical Network Analogy: Surface and Space Resistances

The brilliance of the radiosity method is its translation of radiative exchange into an equivalent electrical circuit, where heat flow ($Q$, in Watts) is analogous to current, and the driving potential is a difference in radiosity or emissive power. This circuit is built from two fundamental types of resistors.

The surface resistance arises because a real surface (with $ϵ<1$) does not emit as much energy as a perfect blackbody at the same temperature. The net rate at which radiation leaves a surface due to emission and reflection is $Q=A(J-G)$. Using the relationship between $J$, $G$, $ϵ$, and the blackbody emissive power $E_b=\sigma T^4$, we can derive the surface heat transfer as: $$Q=\frac{E_b-J}{(1-ϵ)/(ϵA)}$$ This equation takes the form of Ohm's Law ($\text{current}=\frac{\text{potential difference}}{\text{resistance}}$). Here, the driving potential is $(E_b-J)$, and the surface resistance is: $$R_{surface}=\frac{1-ϵ}{ϵA}$$ This resistance is a property of the surface itself; a higher emissivity ($ϵ\to 1$) leads to a smaller surface resistance, meaning the surface more readily releases its emissive power.

The space resistance (or geometric resistance) governs the exchange of radiation between surfaces. It depends on how well one surface "sees" another, quantified by the view factor $F_{ij}$ (the fraction of radiation leaving surface $i$ that strikes surface $j$). The radiative exchange between two black surfaces ($ϵ=1$) would be $Q_{i\to j}=A_i{F}_{ij}\sigma (T_i^4-T_j^4)$. In the network analogy for diffuse surfaces, the driving potential is the difference in their radiosities ($J_i-J_j$). The resistance to this exchange is the space resistance: $$R_{space}=\frac{1}{A_i{F}_{ij}}$$ This resistance depends solely on geometry. A larger view factor or a larger area results in a smaller space resistance, allowing for greater radiative coupling.

Constructing and Solving the Radiation Network

For an enclosure consisting of $N$ gray, diffuse surfaces, we construct a complete network. Each surface $i$ is represented by a node for its blackbody emissive power ($E_{b,i}$) and a node for its radiosity ($J_i$). These two nodes are connected by the surface resistance $R_{s,i}=(1-{ϵ}_i)/({ϵ}_i{A}_i)$. Every radiosity node $J_i$ is then connected to every other radiosity node $J_j$ via a space resistance $R_{ij}=1/(A_i{F}_{ij})$.

To solve the network, we apply the foundational principle of circuit analysis: the net current flow into any node must be zero (conservation of energy). For a radiosity node $J_i$, the sum of all heat flows from connected nodes equals zero: $$\frac{E_{b,i}-J_i}{R_{s,i}}+\sum _{j=1}^N\frac{J_j-J_i}{R_{ij}}=0$$ For each node, we have one such equation. The known quantities are typically surface temperatures (giving $E_b$), surface properties ($ϵ$, $A$), and the view factors (which must be determined independently and obey reciprocity and summation rules). The unknowns are the radiosities ($J_i$). Solving this system of $N$ linear equations yields the values for all $J_i$.

Once radiosities are known, the net radiation heat transfer from each surface $i$ is easily found by calculating the current through its surface resistor: $$Q_i=\frac{E_{b,i}-J_i}{R_{s,i}}=\frac{E_{b,i}-J_i}{(1-{ϵ}_i)/({ϵ}_i{A}_i)}$$ This $Q_i$ is the final answer for the surface's heat loss or gain. The network solution yields heat transfer between all surface pairs indirectly; the exchange between any two surfaces $i$ and $j$ can be found as $Q_{ij}=(J_i-J_j)/R_{ij}$.

Worked Example: A Three-Surface Enclosure

Consider a long triangular duct, a common simplification. The three surfaces are: Surface 1 (hot floor, $T_1$, ${ϵ}_1$), Surface 2 (side wall, $T_2$, ${ϵ}_2$), and Surface 3 (side wall/insulated, $T_3$ unknown, ${ϵ}_3$). Surface 3 is adiabatic, meaning $Q_3=0$.

Step 1: Determine View Factors. For an enclosure, the summation rule holds: $F_{11}+F_{12}+F_{13}=1$. For long ducts, $F_{11}=0$. From geometry, we can often find $F_{12}$ and $F_{13}$ (e.g., using crossed-string method or tables). Reciprocity ($A_1{F}_{12}=A_2{F}_{21}$) gives the others.

Step 2: Draw the Network. Create nodes $E_{b1}$, $J_1$, $E_{b2}$, $J_2$, $E_{b3}$, $J_3$. Connect each $E_b$ to its $J$ with $R_s$. Connect $J_1$ to $J_2$ with $R_{12}=1/(A_1{F}_{12})$, $J_1$ to $J_3$ with $R_{13}$, and $J_2$ to $J_3$ with $R_{23}$.

Step 3: Apply Nodal Equations. Write energy balance at $J_1$, $J_2$, and $J_3$. For the adiabatic surface 3, note that $Q_3=0$ implies $E_{b3}=J_3$, so its surface resistance has no current flow. This simplifies the network.

Step 4: Solve. Insert known values ($T_1$, $T_2$, $ϵ$, $A$, $F_{ij}$) into the three nodal equations. Solve the linear system for $J_1$, $J_2$, and $J_3$ (and thus $E_{b3}=J_3$). Finally, calculate $Q_1=-Q_2=(E_{b1}-J_1)/R_{s1}$.

Common Pitfalls

1. Ignoring the Enclosure Assumption: The radiosity network method strictly requires that all surfaces form a complete enclosure. A "surface" must account for all radiation leaving another surface. A common mistake is analyzing a single object in open space without defining an artificial "surroundings" surface to close the system.
2. Misapplying Surface and Space Resistances: Confusing when to use each resistance is frequent. Remember: Surface resistance ($(1-ϵ)/(ϵA)$) connects a surface's temperature potential ($E_b$) to its radiosity potential ($J$). Space resistance ($1/(A_i{F}_{ij})$) connects the radiosity potentials of two different surfaces. Using the wrong resistance leads to physically meaningless results.
3. Forgetting that Radiosity is a Flux Potential: Treating radiosity $J$ as energy can cause dimensional errors. $J$ is a flux (W/m²), while $Q$ is power (W). The resistances in the network have units of m⁻² or W⁻¹, making $Q=\Delta J/R$ dimensionally consistent in watts.
4. Neglecting the Adiabatic Surface Condition: An insulated (reradiating) surface has a net heat transfer of zero, but it is not inactive. It absorbs and re-emits radiation, influencing the heat exchange between other surfaces. The correct treatment is to set $Q=0$, which forces $E_b=J$ for that surface, effectively shorting its surface resistance in the network.

Summary

• The radiosity method models complex radiative exchange in an enclosure of gray, diffuse surfaces by converting it into an equivalent electrical circuit problem, where heat flow is current and radiosities or emissive powers are potentials.
• The two key resistances are the surface resistance $R_{surface}=(1-ϵ)/(ϵA)$, which depends on a surface's emissivity, and the space resistance $R_{space}=1/(A{F}_{ij})$, which depends solely on the geometry and view factor between surfaces.
• Constructing the network involves creating a node for the blackbody emissive power $E_b$ and radiosity $J$ for each surface, connecting them with their surface resistances, and linking radiosity nodes with space resistances.
• Solving the system of linear nodal equations (energy balances) at each radiosity node yields all unknown radiosities, from which the net heat transfer for every surface is calculated.
• This method provides a clear, systematic, and computationally manageable framework for solving real-world radiation heat transfer problems in engineering design, from electronic cooling to thermal processing.