灰色体间的辐射交换

灰体辐射换热的理论基础

作为前提的三个假定

面间辐射换热计算的标准模型"灰体漫射面的封闭腔理论"建立在三个假定之上：①灰体（辐射特性与波长无关，适用界限参见基尔霍夫定律）、②漫射（发射与反射与方向无关）、③各面等温且辐射均匀。只要接受这三条，任意面数的辐射交换都能化为线性代数——这是一套性价比极高的工程理论。反过来说，当结果与实测对不上时，这三条也正是应当逐一怀疑的假定清单。

角系数——描述几何的语言

由面 \( i \) 射出的辐射中直接到达面 \( j \) 的份额就是角系数 \( F_{ij} \)。它是纯几何量，并且恒满足两条恒等式。

$$ A_i F_{ij} = A_j F_{ji} \;\;(\text{相对性}), \qquad \sum_j F_{ij} = 1 \;\;(\text{完整性、封闭腔}) $$

这两式可以直接当作角系数计算结果的校核式使用。不满足完整性关系的角系数集合会导致计算中能量泄漏（或凭空产生），因此辐射分析的验证总是从这里开始。

辐射热阻网络——两面系的完整解

灰体面之间的换热用电路类比看得最清楚。给各面配上表面热阻 \( (1-\varepsilon)/(\varepsilon A) \)、面间配上空间热阻 \( 1/(A_1 F_{12}) \)，则两面封闭腔的净换热量为

$$ Q_{12} = \frac{\sigma (T_1^4 - T_2^4)}{\dfrac{1-\varepsilon_1}{\varepsilon_1 A_1} + \dfrac{1}{A_1 F_{12}} + \dfrac{1-\varepsilon_2}{\varepsilon_2 A_2}} $$

大平行平板（\( F_{12}=1, A_1=A_2 \)）时化为 \( q = \sigma(T_1^4 - T_2^4)/(1/\varepsilon_1 + 1/\varepsilon_2 - 1) \)，大空间中的小物体则包络一侧的热阻消失，得 \( Q = \varepsilon_1 A_1 \sigma (T_1^4 - T_2^4) \)——实务中常用的这一系列公式，全都是同一个网络的特例。"ε越小、表面热阻越大"这一电路直觉，原封不动地解释了低发射率面（抛光金属、隔热屏）为何有效。

多面系的计算方法

有效辐射法（辐射度法）——N面系的标准解法

面数一多，就要对各面的有效辐射 \( J_i \)（自身发射＋反射的总射出）求解联立方程组。

$$ J_i = \varepsilon_i \sigma T_i^4 + (1-\varepsilon_i) \sum_j F_{ij} J_j $$

温度已知面与热流已知面（绝热面、重辐射面是 \( q=0 \) 的特例）混杂在一起也照样能解，解完即可得到各面的净热流密度。CFD／FEM工具中的"S2S（surface-to-surface）"模型就是它的实现。数值上的主角不是方程组而是角系数矩阵的计算：面数为N时需要 \( N^2 \) 个几何积分（含遮挡判定），计算成本与精度都由此决定。

角系数的计算方法与验证

无论用哪种方法，都要对算出的矩阵做完整性关系与相对性关系的满足度检查。工具通常提供按完整性关系归一化的选项，但归一化之前误差就很大（例如行和不足0.95）就是光线数或分辨率不够的信号，归一化只是把误差藏起来而已。

面划分的设计——守住"等温面"假定

理论中的假定③（各面等温）翻译成划分准则就是：不要把温度梯度大的面当作一片来处理。把长长的散热肋片当成一个面，根部与顶端的温差就被抹平，辐射量随之算错。反过来，划分过细会让角系数计算爆炸，因此实务上的经验尺度是"按预期温差在数十K以内的单位来划分"。S2S中的聚类设置（把若干网格面并成一个辐射面）正是这一取舍的调节阀。

实务应用指南

手算实例——建立数字的量级感

辐射屏的设计

在两个面之间插入一片低ε的辐射屏，辐射热阻串联增加，热流大致减半（ε相同时N片给出 \( 1/(N+1) \) 倍）。多层隔热材料（MLI）就是这一原理的极致，航天器用几十层实现了0.01量级的有效发射率。设计上有两点要注意：①ε越小屏蔽效果越显著（表面热阻支配）；②实物中间隔件与缝合处的接触导热以及边缘泄漏会大幅侵蚀理论值——MLI的实际性能只有理论值的几分之一是常态，实务上要用实测系数（有效发射率）来设计。

辐射的线性化——并入热网络与瞬态分析

温差不大时，辐射可以线性化，按与对流相同的形式处理：\( h_r = 4\varepsilon\sigma T_m^3 \)（\( T_m \) 为平均温度）。室温附近立刻得到 \( h_r \approx 5\varepsilon \) W/m²K——与自然对流同量级的量级感就此成形。这是把辐射廉价地并入热阻网络、一维分析的常规手法，但温差较大（比值超过1.5倍）时线性化误差不可忽略，此时应保留四次方求解或采用迭代更新。

工具中的处理

各工具的设置要点

开放系统的建模——别忘了"虚拟面"

封闭腔理论以闭空间为前提。对于敞开的机箱、室外设备，惯常做法是在开口处张一个环境温度的虚拟黑体面把空间封闭起来，忘了这一步就会破坏完整性关系、让热量凭空消失。有些工具的"environment temperature"设定会自动承担这一角色，但那个环境温度代表什么（天空、地面、室内壁面）是分析者必须作为物理来判断的事项。在室外，天空有效温度的取法（晴朗夜间相当于比外气低10～20 K）会左右结果。

前沿动向

向镜反射与光谱的扩展

灰体漫射之外，还有部分镜反射（拆分为漫射＋镜反射，进而用BRDF一般化）、波段划分（非灰）乃至含偏振的精细模型。随着光线追踪的GPU化，过去只能"将就用灰体漫射"的大规模模型也开始能做蒙特卡洛分光计算，在前照灯、光学系统、高温炉的设计中正逐步被采用。实务上的常规策略是分层推进：先用灰体漫射把握全局，只对灵敏度高的面做精细化。

与参与性介质（气体、颗粒）的耦合

当燃烧气体（CO₂、H₂O）、烟炱、喷雾介入时，面间交换之上还要叠加介质的吸收、发射与散射，主题就变成辐射输运解法（DO/DOM、P1、蒙特卡洛）与气体非灰模型（WSGG、SLW、k分布）的组合选型。面的S2S理论把"介质透明"这一近似摆到了明处——认识到这一近似会在炉内、发动机缸内失效，正是模型选择的出发点。

近场辐射——理论之外

当面间距离接近热辐射的波长（～10 μm）时，倏逝波的隧穿会带来超出黑体极限若干数量级的热输运，这就是近场辐射。在MEMS、热光伏、纳米间隙器件中实测不断推进，适用于传统角系数理论已经失效的尺度的设计理论正在形成。对宏观设计者而言，目前需要掌握的是这条界线的知识：在微米级间隙里，主宰的是另一套物理。

故障排查

按症状的原因与对策

辐射分析的验证套路

相关文章：基尔霍夫定律、辐射传热文章列表、NAFEMS T4（边界条件验证的套路）。