4. 角系数数值计算方法与验证

角系数是热辐射交换模型的核心几何量，它只取决于两个表面的大小、形状、相对位置与朝向，与温度、发射率无关。对于平板、无限长通道、同心球等规则构型，教材给出了解析解；但工程中大量存在的遮挡、有限尺寸、非共面倾斜表面，往往没有闭式解。此时必须借助数值方法。本节围绕三种互补的数值路线展开：直接积分法、蒙特卡洛法与半解析法，并说明如何用解析解或标准数据完成验证。

一、三种方法的数学本质与适用边界

角系数的定义是双重面积分：

F_ij = (1/A_i) ∫_Ai ∫_Aj  cosθ_i · cosθ_j / (π · r²)  dA_j dA_i

其中 r 为两微元连线长度，θ_i、θ_j 分别为连线与各自表面法线的夹角，只有两微元互相可见时被积函数才非零。三种数值方法本质上是对这个积分采取不同的离散策略。

直接积分法把两个表面都剖分为网格，对每对微元计算核函数并做数值求积（如高斯-勒让德求积或中点法则），再对可见性做遮挡判断。它的优点是精度可控、误差随网格加密单调收敛；缺点是计算量为 O(N_i·N_j)，遮挡判断还会进一步增加开销，适合中小规模、几何较规则的构型。

蒙特卡洛法从表面 i 上按面积均匀采样点，在半球方向按余弦加权分布发射射线，追踪射线与场景的交点，统计落在表面 j 上的射线比例即为 F_ij。它天然处理遮挡与复杂几何，实现简单，但收敛速度只有 O(1/√N)，结果带有统计噪声，适合复杂三维场景或作为其他方法的交叉校验。

半解析法利用角系数的可加性与互换性，把复杂构型分解为若干具有解析解的基元（如微元对矩形面、微元对圆盘、平行矩形对），只对剩余的一维或二维做数值积分。典型做法是固定一个表面上的积分点，对另一个表面使用解析公式，从而把二重积分降为一重积分，兼顾精度与效率。

二、直接积分法的实现要点

以两个有限大平行矩形（边长 a×b，间距 c）为例，这是有解析解的经典构型，适合用来标定数值实现。实现时先建立坐标系：表面 i 位于 z=0 平面，法线 +z；表面 j 位于 z=c 平面，法线 −z，两表面中心对齐。对表面 i 用 N×N 高斯点、表面 j 用 M×M 高斯点，则

for p in gauss_points_i:          # 位置 (x_i, y_i, 0)，权重 w_i
    for q in gauss_points_j:      # 位置 (x_j, y_j, c)，权重 w_j
        r_vec = q - p
        r = |r_vec|
        cos_i = (r_vec · n_i) / r
        cos_j = (r_vec · n_j) / r
        if cos_i > 0 and cos_j > 0:      # 互相可见
            K += w_i * w_j * cos_i * cos_j / (π * r^2)
F_ij = K / A_i

这里的关键约定是：法线方向必须由表面指向对方所在半空间，cos 为负说明该微元对不可见，应剔除。若存在第三方遮挡面，还需在 p、q 之间做线段与遮挡面的求交测试。实际调试中，建议先用 8×8 高斯点跑一遍，再加密到 16×16、32×32，观察 F_ij 的变化量是否按预期收敛；若加密后结果反而跳动，多半是遮挡判断的边界容差设置不当。

三、蒙特卡洛法的采样与统计

蒙特卡洛法的核心是“按面积均匀选点、按余弦分布选方向”。在表面 i 上采样时，对矩形可用两个均匀随机数映射到 [0,a]×[0,b]；方向采样则用

u1, u2 ~ U(0,1)
r = sqrt(u1)
phi = 2π·u2
d_local = (r·cosφ, r·sinφ, sqrt(1-u1))   # 局部余弦加权方向
d_world = R · d_local                     # R 为局部到全局的旋转矩阵

随后沿 d_world 发射射线，求与场景所有面的最近交点，判断交点是否落在表面 j 上。重复 N 次，得到命中数 N_hit，则 F_ij ≈ N_hit/N。由于是伯努利试验，标准误差约为 √(F(1−F)/N)。例如 F≈0.3、N=10⁶ 时，标准误差约 4.6×10⁻⁴，通常已满足工程校验需求。实践中应固定随机种子以便复现，并报告置信区间而非单一数值。

四、半解析法的降维思路

半解析法的价值在于“能用解析的地方绝不用数值”。以微元 dA_i 对有限矩形 j 的角系数为例，存在闭式表达式；于是对表面 i 做二维数值积分、对表面 j 用解析式，就把四重积分降为二重积分。更常见的做法是利用可加性：把复杂表面拆成若干矩形块，分别用平行矩形解析公式计算后再叠加。需要注意，可加性只对“接收面”成立，对发射面拆分时必须按面积加权，且拆分后各块之间不能互相遮挡。

五、验证：三类判据缺一不可

数值结果必须经过验证才能进入热网络模型。建议至少执行以下三类判据：

• 与解析解对比：对平行矩形、同心球、无限长同心圆柱等有闭式解的构型，要求相对误差小于 1%（直接积分/半解析）或落在蒙特卡洛 3σ 置信区间内。
• 互换性检验：A_i·F_ij 应等于 A_j·F_ji，相对偏差通常应小于 10⁻³。这是检验核函数与面积归一化是否正确的最灵敏指标。
• 封闭性检验：对封闭腔，Σ_j F_ij = 1。若明显偏离 1，说明存在漏掉的可见面或遮挡判断错误。

三种方法并非互斥：工程实践中常用半解析法快速给出基准值，用直接积分法在规则区域做高精度校核，再用蒙特卡洛法处理含复杂遮挡的三维场景并交叉验证。把每次计算的构型、网格/样本参数、收敛曲线与三类判据结果一并记录，才能让角系数数据在后续热网络求解中真正可信。