📚 课程目录第3章角系数定义、性质与计算几何基础

3. 角系数定义、性质与计算几何基础

在热辐射交换模型中，角系数（View Factor，又称形状因子、交换因子）是连接几何与辐射能量传递的核心桥梁。它回答一个非常具体的问题：从一个表面发出的漫辐射能量中，有多大比例直接到达另一个表面？本节从积分定义出发，梳理互换性、可加性与封闭性三条基本性质，并说明计算几何在角系数求解中的支撑作用。

一、角系数的积分定义

考虑两个有限面积表面 A₁ 与 A₂，均为漫射、均匀辐射的灰体表面。角系数 F₁₂ 定义为：从 A₁ 发出的辐射能中，直接投射到 A₂ 上的份额。其严格积分表达式为：

F₁₂ = (1/A₁) ∬_A₁ ∬_A₂ [cosθ₁ · cosθ₂ / (π · r²)] dA₂ dA₁

式中：θ₁、θ₂ 分别为连线 r 与两表面法线的夹角；r 为两微元面之间的距离；π 来自朗伯漫射假设下的辐射强度与辐射力关系。这个四重积分（两个面积分）揭示了角系数的本质：它完全由几何形状、相对位置和朝向决定，与温度、发射率、波长等物理量无关。这是角系数最重要的前提——它只描述“几何可见性”，不描述“能量大小”。

理解定义时需注意三个约定：其一，表面必须是漫射的，否则 cosθ 的权重关系不成立；其二，积分中隐含遮挡判断，若连线被其他表面阻断，则该微元对不贡献；其三，F₁₂ 与 A₁ 的大小有关，但 F₁₂·A₁ 是对称量，这直接引出互换性。

二、三条基本性质

互换性（Reciprocity）：A₁F₁₂ = A₂F₂₁。这条性质来自积分表达式中被积函数对两表面完全对称。它的工程价值在于：当两个表面面积悬殊时，可先求小表面对大表面的角系数，再用面积比换算，避免直接计算大表面的积分。例如一块 0.5 m² 的传感器表面对 20 m² 的炉壁，若已知 F_壁→传感器 = 0.02，则 F_传感器→壁 = 20/0.5 × 0.02 = 0.8。

可加性（Additivity）：若 A₂ 可分解为互不重叠的 A₂ₐ 与 A₂ᵦ，则 F₁₂ = F₁₂ₐ + F₁₂ᵦ。可加性要求分解后的子表面之间无遮挡、无重叠，且各自仍满足漫射假设。它使得复杂几何可以拆成若干简单几何的组合，是工程计算中最常用的化简手段。

封闭性（Enclosure / Summation）：对于构成封闭腔的 N 个表面，任一表面 i 满足 Σⱼ Fᵢⱼ = 1（j = 1…N）。封闭性本质是能量守恒：从表面 i 发出的全部辐射必然落在腔内某个表面上。对于非封闭的开放场景，需引入“虚拟表面”补全封闭腔，否则求和小于 1，剩余部分代表逃逸到环境中的能量。

三条性质常联合使用：用可加性拆分复杂面，用互换性转换计算方向，用封闭性校验结果是否自洽。若某表面角系数之和明显偏离 1（封闭腔情形），通常意味着几何分解有重叠、遗漏或遮挡判断错误。

三、计算几何的支撑作用

解析积分仅对少数规则几何（平行矩形、垂直矩形、同轴圆盘、无限长平行平板等）有闭式解。实际工程中的曲面、遮挡、非规则边界，几乎无法手工积分。计算几何在此承担三类支撑：

可见性判定：通过射线与三角形求交，判断两微元面之间是否被遮挡，这是积分中“直接到达”的几何实现。
表面离散与数值积分：将表面剖分为网格微元，用高斯积分或蒙特卡洛采样近似四重积分，网格越细、采样越多，精度越高但代价越大。
几何预处理：法线方向统一、面积计算、坐标变换、对称性识别，为角系数矩阵的构建提供一致的数据基础。

下面这张图概括了从几何输入到角系数矩阵、再到性质校验的完整支撑链路，并标出了三条性质在链路中的作用位置。

四、一个可复现的计算示例

设一个封闭矩形腔由 6 个面组成，其中底面 A₁ = 2 m × 3 m，顶面 A₂ 与底面平行且正对，尺寸相同，间距 1.5 m。其余四面为侧壁。要求 F₁₂。

平行矩形角系数的解析式（Hottel 交叉线法）为：

F₁₂ = (2/(πXY)) · [ ln√((1+X²)(1+Y²)/(1+X²+Y²))
       + X√(1+Y²)·arctan(X/√(1+Y²))
       + Y√(1+X²)·arctan(Y/√(1+X²))
       - X·arctan X - Y·arctan Y ]

其中 X = a/c，Y = b/c，a、b 为矩形边长，c 为间距。代入 a=2，b=3，c=1.5，得 X=1.333，Y=2.0。按上式逐项计算（示例数值，供教学演示）：ln 项约 0.427，X 项约 1.333×2.236×0.538≈1.604，Y 项约 2.0×1.667×0.927≈3.090，减去 X·arctanX≈1.333×0.927≈1.236，减去 Y·arctanY≈2.0×1.107≈2.214，合计约 1.671，乘以 2/(π×1.333×2.0)≈0.2387，得 F₁₂ ≈ 0.399。

校验：由对称性 F₂₁ = F₁₂ ≈ 0.399。底面到四个侧壁的角系数之和应为 1 − 0.399 = 0.601，若四个侧壁几何对称，则每个约 0.150。用封闭性 ΣF₁ⱼ = 0.399 + 4×0.150 ≈ 0.999，与 1 的偏差来自四舍五入，说明结果自洽。这一“解析计算 + 封闭性校验”的组合，正是工程中验证数值角系数程序的标准做法。

当表面间存在遮挡或表面为曲面时，解析式不再适用，必须依赖计算几何的射线求交与数值积分。此时封闭性校验尤为关键：若 ΣFᵢⱼ 与 1 偏差超过 1%～2%，应优先排查网格法线朝向、遮挡判定和虚拟表面补全是否正确，而不是盲目加密网格。

综上，角系数的积分定义给出几何本质，互换性、可加性、封闭性提供化简与校验工具，计算几何则把不可积的复杂场景转化为可执行的数值流程。三者结合，才能在实际热辐射交换模型中稳定地构建角系数矩阵。