辐射度法热辐射传热分析

Stefan-Boltzmann定律 $q = \varepsilon \sigma T^4$ 只能算单个表面向外辐射的热流，处理不了"多个表面互相照射、反射"这样的封闭腔辐射问题。辐射度法（Radiosity Method）在漫射灰体假设下，把每个表面的辐射度 $J$（离开表面的总辐射通量，包括自身热辐射和反射）作为未知量，建立联立方程组求解封闭腔内所有表面的辐射热交换。工业炉膛设计、电子散热、航天器热控都依赖这套方法。

辐射度 $J_i$ 是离开表面 $i$ 的所有辐射，包含两部分：自身热辐射（$\varepsilon_i \sigma T_i^4$）和对入射辐射的反射（$(1-\varepsilon_i)G_i$）：

$\varepsilon_i$ 是表面发射率（漫射灰体），$G_i$ 是到达表面 $i$ 的入射辐射（Irradiation）。$G_i$ 又来自其他表面的辐射度，通过角系数（View Factor）联系：$G_i = \sum_j F_{ij} J_j$。代入后形成 $N$ 个关于 $J_i$ 的联立方程组——这就是辐射度法的核心。

角系数 $F_{ij}$ 表示"从表面 $i$ 出发的总辐射能中，有多少比例直接到达表面 $j$"。数学定义（漫射均匀辐射）：

$\theta_i, \theta_j$ 是各表面法线与两面元连线的夹角，$r$ 是两面元距离。

角系数的两个重要性质（必须验证）：

• 完整性关系（Summation Rule）： $\sum_{j=1}^N F_{ij} = 1$（表面 $i$ 发出的辐射全部落在封闭系统某处）
• 互易关系（Reciprocity）： $A_i F_{ij} = A_j F_{ji}$（利用这个性质，$N$ 个表面只需计算 $N(N-1)/2$ 个独立角系数）

有，部分标准几何有封闭形式。几个常用例子：

• 两个平行的无限大平板（相距 $h$，宽 $L$）： $F_{12} = 1$（无限大），$F_{11} = 0$（平面）
• 两个同心球（半径 $r_1 < r_2$）： $F_{12} = 1$，$F_{21} = (r_1/r_2)^2$
• 有限矩形平板之间： 交叉绳方法（Crossed-String Method，2D）或Hottel代数法（3D），形式复杂
• 凸体对自身： $F_{11} = 0$（凸体不能"看到"自身）

大量角系数解析公式汇编在Howell's "Catalog of View Factors"（网络免费获取）。对于复杂三维几何则需要数值方法（半立方体法、射线追踪）。

将 $G_i = \sum_j F_{ij} J_j$ 代入辐射度定义 $J_i = \varepsilon_i \sigma T_i^4 + (1-\varepsilon_i)G_i$，得到 $N$ 个联立方程：

矩阵形式 $\mathbf{A}\mathbf{J} = \mathbf{b}$，其中：

$\mathbf{A}$ 是 $N \times N$ 满矩阵（原则上每面和其他所有面都有角系数关联）。求解后，表面 $i$ 的净辐射热流（正为离开表面）：

两种边界条件：

① 规定温度边界（最常见）： 给定 $T_i$，方程右侧 $b_i = \varepsilon_i \sigma T_i^4$ 已知，求解方程组得到 $J_i$ 和 $q_i$。黑体（$\varepsilon_i = 1$）时 $J_i = \sigma T_i^4$（辐射度等于黑体辐射力），自动满足，无需列入方程。

② 规定热流边界（绝热/镜面/指定热流）： 给定 $q_i$，用 $\sigma T_i^4 = J_i + q_i(1-\varepsilon_i)/\varepsilon_i$ 替换，将 $T_i$ 作为未知量求解。绝热表面（$q_i = 0$）→ "重辐射（Reradiating）表面"，$J_i = G_i$，在炉膛中的砖墙可视为重辐射表面。

确实是个瓶颈。$N$ 个表面的计算量分析：

实际应对策略：

• 遮挡稀疏化： 互相遮挡（看不见对方）的表面对，角系数精确为零，无需计算存储，利用稀疏性大幅减少内存。光线追踪或半立方体法可自动识别遮挡。
• 分组（Clustering）： 将距离远的小表面聚合为等效大面，减少有效 $N$。
• 蒙特卡洛替代： $N > 10000$ 时，Monte Carlo法内存占用 $O(N)$，比辐射度法更实用。

主要有三种数值方法：

① 半立方体法（Hemi-Cube Method）： 在每个表面法线方向放置虚拟半立方体，将其他表面投影到半立方体各面，统计投影面积比例。Ansys Fluent S2S（Surface-to-Surface）默认方法。分辨率（HCU_RESOLUTION参数）越高，精度越高但越慢。对于有遮挡的复杂几何，遮挡判断用Z-buffer算法。

② 射线追踪法（Ray Tracing）： 从每个表面随机发射 $M$ 条射线，统计落在各表面的比例。Ansys Mechanical默认方法。精度正比于 $1/\sqrt{M}$，$M = 10^4$ 时误差约1%。

③ 高斯积分法： 对角系数积分式用数值积分（Gauss-Legendre），适合简单几何（矩形、圆形），精度高但只对规则形状有效。

质量检查：完成计算后验证 $\sum_j F_{ij} = 1$ 是否成立（允许误差<1%），不满足说明封闭性假设被违反或计算精度不足。

对于大多数电子散热场景（温度<300℃、空气介质、表面近似漫射灰体），辐射度法完全够用。实践注意点：

1. 发射率设置（最重要！）： 裸铝约0.03～0.1（极低！），氧化铝约0.8，喷漆铝约0.85～0.95。同一面积，发射率从0.05提升到0.9，辐射换热量增加18倍——发射率对结果影响远大于几何精度。
2. 封闭腔构建： 辐射度法要求所有表面构成封闭腔（或加虚拟开口面）。忘记某个开口面会导致完整性关系 $\sum F_{ij} \ne 1$，计算发散。
3. 与对流耦合： 电子散热中辐射和自然对流量级相当，必须同时考虑（联合求解）。单独只算辐射或只算对流均误差很大。
4. 角系数完整性验证： 每行 $\sum_j F_{ij}$ 应为1.00±0.01，超出说明封闭性定义有问题。

辐射度法假设介质透明（非参与性），不适合处理CO₂/H₂O气体辐射。对于含参与性介质的问题（工业炉膛、燃烧室），需要其他方法：

• 离散坐标法（DOM，SN法）： 把辐射传输方程离散为 $N$ 个方向，Ansys Fluent默认的工业燃烧辐射方法，精度和效率均好，支持谱带气体模型（WSGG）
• P-N近似（P1最常用）： 将辐射强度展开为球谐函数，P1方程扩散形式简单、易实现，对光学厚介质精度好
• 蒙特卡洛法（MCM）： 追踪单个光子路径，可处理最复杂的几何和谱带，但计算量大

经验判断：工业空气炉膛（1000℃以上，含CO₂/H₂O）→DOM+WSGG；真空电阻炉（不含气体）→辐射度法；LED灯光学设计（复杂散射）→蒙特卡洛。

主要工具的热辐射功能对比：

几个容易踩坑的地方：

• 封闭性： S2S要求模型是封闭腔。如果有开口（如通风孔），需要创建"开口面"并设为黑体温度边界，否则完整性关系不满足
• Hemi-Cube分辨率： 默认10往往不够，对有小间隙或高纵横比几何推荐设为50～100；分辨率太低会导致角系数误差>5%
• 集群（Clustering）设置： Fluent默认启用表面集群，会将角系数光滑化，对细节精度有影响。精细分析时关闭集群
• 计算顺序： 角系数必须在第一次迭代前计算完毕。更改表面定义后必须重算角系数
• 半透明材料： S2S不支持半透明介质（如玻璃），需用DO模型

辐射度法的局限在于漫射灰体假设，现实中两个条件常常不满足：

• 谱带辐射度法（Band Radiosity）： 将波长划分为多个谱带（例如每1000 nm一段），各谱带独立求解，谱带内取均匀发射率。适合Low-E玻璃（$\varepsilon$ 随波长急剧变化）、太阳能集热器（区分太阳光谱和热红外）。
• 双向反射分布函数（BRDF）集成： 将辐射度法扩展到处理方向性反射（非漫射表面），在太阳能抛物槽收集器、LED反射镜、激光热处理中有应用。
• 机器学习角系数预测： 用图神经网络（GNN）或卷积神经网络（CNN）从几何形状直接预测角系数，避免大规模积分计算，速度提升100倍以上（2022-2024年多篇论文）。
• GPU加速蒙特卡洛辐射： 百万表面级别的辐射计算在GPU上可实现实时仿真（NVIDIA CUDA MCM），用于数字孪生的在线热监控，已在半导体制造炉管控制中实验应用。
• 量子辐射效应： 当两表面间距<微米级（MEMS器件），近场辐射（Nearfield Radiation）超越黑体Planck定律上限，经典辐射度法完全失效，需要基于涨落-耗散定理的近场辐射理论。