2. 热辐射基本定律与交换模型物理基础

在建立任何热辐射交换模型之前，必须先厘清四条基本定律各自回答什么问题，以及发射率、吸收率、反射率在模型中扮演什么角色。很多数值验证失败，并非求解器出错，而是物理量用错了位置：把光谱量当成全量、把方向量当成半球量、把发射率当成吸收率。本节围绕这一主线展开。

四条定律的分工

普朗克定律给出黑体光谱辐射力随波长与温度分布的完整描述，是辐射模型中最底层的“源项”。在波长 λ、温度 T 下，黑体光谱辐射力为

E_bλ(λ,T) = C1 / (λ^5 · (exp(C2/(λT)) − 1))
C1 = 3.7418e-16 W·m²      C2 = 1.4388e-2 m·K

它的工程意义在于：任何真实表面的光谱发射都由 E_bλ 乘以光谱发射率得到，任何气体或颗粒的参与性介质吸收也以它为权重。交换模型中若需要分波段处理（如太阳短波与常温长波分开），普朗克定律就是划分波段并计算波段内黑体份额的依据。

斯特藩-玻尔兹曼定律是普朗克定律在全波段上的积分，给出黑体全辐射力 E_b = σT⁴，σ = 5.670×10⁻⁸ W/(m²·K⁴)。它把复杂的谱分布压缩成一个温度的四次方关系，是表面辐射换热网络（如热阻网络、角系数法）中计算节点辐射出射的核心公式。需要注意，σT⁴ 只对黑体严格成立；对灰体要乘发射率，对选择性表面则必须回到分波段积分。

基尔霍夫定律指出，在热平衡条件下，任何表面对某一波长、某一方向的发射率等于其同波长、同方向的吸收率：ε_λθ = α_λθ。它的价值在于把“发射能力”和“吸收能力”这两个看似独立的物性绑定起来。工程上常用的灰体假设（ε_λ = α_λ = 常数）正是基尔霍夫定律在谱平均意义上的简化。但必须警惕：全发射率与全吸收率一般不相等，因为二者对应的入射谱不同——太阳辐射峰值约 0.5 μm，而 300 K 表面辐射峰值约 10 μm，这正是选择性涂层（低太阳吸收、高红外发射）能够工作的物理根源。

兰贝特余弦定律描述漫射表面的方向特性：定向辐射强度与方向无关，而定向辐射力随天顶角按 cosθ 衰减。它决定了半球空间积分时 cosθ 权重的来源，也是角系数计算中“漫射-灰”假设的基础。若表面不是漫射而是镜反射或强方向性，兰贝特假设失效，交换模型必须改用双向反射分布函数描述。

三类物性在交换模型中的角色辨析

把上述定律落到交换模型，表面辐射行为可归结为一个能量分配问题。设入射辐射为 G，则表面吸收 αG、反射 ρG、透射 τG，且 α + ρ + τ = 1；对不透明表面 τ = 0，于是 α + ρ = 1。表面自身的出射为 εσT⁴（灰体近似下）。净辐射换热量为

q_net = εσT⁴ − αG

这个式子清楚显示了两个物性的不同角色：ε 只出现在出射项，α 只出现在入射项。只有在灰体假设下才能令 α = ε 合并简化；一旦入射谱与表面发射谱显著错位（太阳能集热、航天热控、红外隐身），就必须分别取值。反射率 ρ 则通过“多次反射”进入模型：在两表面构成的封闭腔中，A 面反射的能量会再次打到 B 面，形成无穷级数，这正是辐射网络法中表面热阻与空间热阻并联结构的来源。

一个可手算的示例

设两块平行不透明灰板，面积均为 1 m²、间距远小于板宽（角系数 F₁₂ = F₂₁ ≈ 1），T₁ = 500 K、ε₁ = 0.8，T₂ = 300 K、ε₂ = 0.6，无其他表面参与。按灰体假设 α = ε，用两表面网络法：

空间热阻      R_s = 1/(A1·F12) = 1.0
表面热阻      R_1 = (1−ε1)/(A1·ε1) = 0.25
              R_2 = (1−ε2)/(A2·ε2) = 0.667
总热阻        R = 0.25 + 1.0 + 0.667 = 1.917
黑体辐射力差  Eb1 − Eb2 = σ(500⁴ − 300⁴) = 5.670e-8 ×(6.25e10 − 8.1e9)
              ≈ 3084 W/m²
q = (Eb1 − Eb2)/R ≈ 3084 / 1.917 ≈ 1609 W

若误把 ε₁ 当作 α₁ 用于入射项，或忽略反射的多次耦合而只算一次吸收，结果会偏离 10% 以上。验证时可做两个自检：其一，令 ε₁ = ε₂ = 1，结果应退化为 σ(T₁⁴ − T₂⁴) ≈ 3084 W；其二，令 T₁ = T₂，净热流必须为零，否则说明 α 与 ε 的取值破坏了热平衡一致性。