基尔霍夫定律
基尔霍夫定律的理论基础
定律的准确表述——它在"光谱、方向"层面成立
热辐射的基尔霍夫定律指出,处于热平衡的表面发射率与吸收率相等。但其严格表述针对的是波长 \( \lambda \)、方向 \( \theta \)、温度 \( T \) 都对齐的光谱定向量。
$$ \varepsilon(\lambda, \theta, T) = \alpha(\lambda, \theta, T) $$
其来源是细致平衡:在平衡的黑体辐射场中,若存在某个表面在某一波长、某一方向上吸收多于发射,就能不做功而制造温差,违反热力学第二定律。这正是"善于吸收的表面也善于发射"这一直觉的依据,也是黑体(\( \alpha = 1 \))成为最强发射体(\( \varepsilon = 1 \))的原因。
最重要的陷阱——全量上一般 \( \varepsilon \ne \alpha \)
等一下。既然有基尔霍夫定律,材料表里出现"发射率0.9、吸收率0.3"这样的数字不是很奇怪吗?
这正是这条定律最要命的考点。等号成立的前提是同一波长之间。实务中用的"全发射率""全吸收率"都是光谱量的加权平均,而且权重不同:发射率按表面自身温度的普朗克分布加权,吸收率则按入射辐射的光谱加权。想想刷成白色的屋顶——太阳光(5800 K的光谱、以可见光为中心)被很好地反射掉,所以 \( \alpha_{solar} \approx 0.3 \);但屋顶自身是300 K,在红外波段发射,所以 \( \varepsilon_{IR} \approx 0.9 \)。任何波长上 \( \varepsilon(\lambda) = \alpha(\lambda) \) 都依然成立,而全量上0.3与0.9却可以并存。这不是矛盾,而是平均权重的差别——想通这一点,基尔霍夫这一关就算过了。
灰体近似的意义与适用界限
假定辐射特性与波长无关(\( \varepsilon(\lambda) = \) 常数)就是灰体近似;只有在这一前提下,全量上 \( \varepsilon = \alpha \) 才成立,计算也随之大幅简化。判断能否适用的准则是"在所涉辐射源光谱相互重叠的波段内,实际的光谱特性是否大致平坦"。温度相近的表面之间的辐射交换(炉内、电子设备内)用灰体近似基本妥当。而太阳光对常温物体、高温火焰对低温壁面这类光谱相距很远的组合,灰体近似就会成为系统误差的主因。
计算中的处理
带平均的做法
这里把由光谱数据生成实用值的公式明确写出来。温度为 \( T_s \) 的表面的全发射率,以及对光谱为 \( G(\lambda) \) 的入射的全吸收率为
$$ \varepsilon = \frac{\int \varepsilon(\lambda)\, E_{b\lambda}(T_s)\, d\lambda}{\int E_{b\lambda}(T_s)\, d\lambda}, \qquad \alpha = \frac{\int \alpha(\lambda)\, G(\lambda)\, d\lambda}{\int G(\lambda)\, d\lambda} $$
其中 \( E_{b\lambda} \) 是普朗克分布。实务上往往不做连续积分,而采用双波段(半灰)近似——太阳波段(~2.5 μm)与红外波段(2.5 μm~)分别取不同常数——它在精度与工作量之间取得了很好的平衡,已成为室外设备与航天器热设计的标准做法。航天热设计中"\( \alpha_s \)(太阳吸收率)与 \( \varepsilon \)(红外发射率)"这一经典搭配,正是这种双波段近似的产物。
并入辐射换热计算
面间辐射的辐射度法(S2S、基于角系数)中,各面的 \( \varepsilon \) 是在灰体、漫射的前提下使用的。与基尔霍夫定律有关、需要确认的有三点:①各面的 \( \varepsilon \) 是否是在该面的温度区间下评价的,②存在太阳等外部光谱源时,吸收一侧是否用 \( \alpha_s \) 单独处理(工具的"solar load"功能内置了这种分离),③是否把以漫射为前提的S2S用在了镜反射很强的表面(抛光金属)上。在涉及气体辐射(CO₂、H₂O的带状吸收)的燃烧系统中,不仅表面、连介质也具有很强的非灰性,因此带模型(WSGG、SLW等)的选择就成了支配因素。
发射率数据的现实——表面状态决定一切
手册里的 \( \varepsilon \) 是"在该表面状态下的值"。铝在抛光面是0.05、阳极氧化后是0.8——同一材料相差一个数量级以上,这正是辐射领域的常态:氧化、粗糙度、污垢、涂装的影响都比母材本身更大。分析输入的 \( \varepsilon \) 应当附上:①标明了表面处理的出处,②对老化(氧化推进)的设想,③条件允许时的实测(发射率仪、FTIR)。辐射计算的不确定度收支,通常由ε的输入不确定度而非模型所支配。
实务用法与误用防范
典型误用模式
| 误用 | 会发生什么 | 正确做法 |
|---|---|---|
| 对阳光下的物体套用 \( \alpha = \varepsilon_{IR} \) | 白色面、选择性表面上的吸热量被高估/低估数倍 | 双波段近似(\( \alpha_s \) 与 \( \varepsilon_{IR} \) 分别输入) |
| 抛光金属的 \( \varepsilon \) 用了氧化面的值 | 辐射损失被高估一个数量级 | 改用明确了表面处理的值,老化用敏感性工况处理 |
| 辐射温度计设定的ε与对象不一致 | 低ε面(金属光泽)上出现数十°C量级的误读,并叠加环境映射 | 按对象设定ε,并用黑体胶带辅助校准 |
| 高温炉壁对低温被加热物用单一灰体ε计算 | 光谱不匹配所致的系统误差 | 分波段的ε,或采用非灰模型 |
把"ε≠α"变成武器——选择性表面
全量上能做到 \( \varepsilon \ne \alpha \),这不是约束,而是设计自由度。太阳能集热器的选择性吸收膜(\( \alpha_s \approx 0.95, \varepsilon_{IR} \approx 0.1 \):吸进来的热不再靠再辐射跑掉)、航天器的白色涂层与OSR(\( \alpha_s \) 小、\( \varepsilon_{IR} \) 大:挡住日照、把排热做到最大)、近年的日间辐射冷却材料(只在大气窗口8~13 μm强烈发射)——它们都是在恪守光谱基尔霍夫定律的同时去设计带平均的技术。就分析一侧而言,一旦要处理这类表面,放弃灰体近似的决定也就自动作出了。
辐射分析评审的检查清单
- 是否确认了外部光谱源(太阳、高温火焰、灯具)的有无,若有是否做了波段分离
- 各面ε的出处是否明确标注了表面处理与温度
- 是否考察了ε的不确定度(±0.05~0.1很常见)会让结果变动多少
- 是否检查了把S2S(以漫射为前提)用于强镜反射面的情况
- 是否确认了辐射温度计所得验证数据本身的ε设定与环境映射
工具的支持情况
各工具的辐射模型与非灰支持
| 工具 | 辐射模型 | 非灰(分波段)支持 |
|---|---|---|
| Ansys Fluent | S2S、DO、Monte Carlo等+Solar Load模型 | DO支持非灰波段划分。太阳波段由Solar Load分离 |
| STAR-CCM+ | S2S、DOM+太阳辐射 | 多波段(按波段设定表面特性) |
| Ansys Mechanical/通用FEM | 面间辐射(灰体漫射) | 波段分离能力有限。实务上把太阳输入手工分离为热流密度 |
| TAITherm、ESATAN-TMS等(航天、车辆热) | 辐射专业软件。α_s/ε_IR双波段是标准设计变量 | 轨道日照、反照率、行星红外一并涵盖 |
| COMSOL | 面对面辐射+半透明介质 | 覆盖多个波段的多光谱分波段 |
设置时的确认点
无论用哪个工具,都要确认三点:①"emissivity"一栏是只用于红外侧(自身发射),还是同时用于吸收,②太阳吸收率能否单独输入(不能输入的工具,就用 \( \alpha_s \times \) 日照手算出太阳热输入,再换算为热流密度边界),③半透明材料(玻璃、树脂)透射波段的处理方式。把手册的辐射章节精读一遍,弄清其中"gray""spectral""band"各词是如何定义的,就能避免大部分设置事故。
前沿动向
日间辐射冷却与超表面的光谱设计
只在大气透射窗口(8~13 μm)提高发射率、同时抑制太阳波段吸收的表面,即使在白天也能实现比环境空气低数°C的无源冷却。用光子结构与超表面自由设计辐射光谱是材料研究的最前沿,在建筑外围护、车辆、纺织品上的应用正迅速走向实用。热分析一侧,能处理这类强非灰表面的带模型建设,则是需要跟上的课题。
非互易材料与"基尔霍夫定律的破缺"
在磁光材料、时间调制结构等打破互易性的体系中,理论与实验都表明:即便同一波长、同一方向,也能做到 \( \varepsilon(\lambda,\theta) \ne \alpha(\lambda,\theta) \)(因为细致平衡所依赖的互易性不再成立,所以与热力学并不矛盾)。若能独立控制发射与吸收,就有望改写热光伏发电的效率极限,因此这在基础与应用两方面都很活跃。而普通工业材料满足互易性,实务中基尔霍夫定律的用法并不会因此改变——这条界线也应当把握住。
光谱测量与原位辨识的进展
基于FTIR的光谱发射率测量、结合热像仪对ε与温度的同时反演、用机器学习从光谱推断表面状态等,正在推动从"假定ε"转向"测量并辨识ε"。在辐射占主导的高温工艺(锻造、烧结、半导体RTP)中,含ε辨识的温度控制直接关系到产品质量,因此产业需求正牵引着这方面的研究。
故障排查
按症状的原因与对策
| 症状 | 可能原因 | 对策 |
|---|---|---|
| 室外设备的温度预测与实测大幅偏离 | α_s与ε_IR混淆(用单一ε计算太阳吸收) | 改为双波段输入。使用工具的Solar Load功能 |
| 辐射温度计与热电偶对不上 | ε设定有误、低ε面的环境映射 | 按对象设定ε,在贴黑体胶带处校准,斜视以避开映射 |
| 金属零件的辐射损失大于实测 | 把氧化面的ε值用在了光亮面上 | 改用明确了表面状态的值。用敏感性分析给出取值范围 |
| 炉内分析中被加热物的升温对不上 | 高温源对低温面的非灰性、忽略了气体辐射的波段特性 | 改用分波段ε与WSGG等非灰气体模型 |
| 低发射率设定的面在工具中呈现异常高温/低温 | ε→0时辐射平衡在数值上变得敏感、多次反射计算不足 | 把ε的下限取为现实值(≥0.02~0.05),加强反射计算的收敛设置 |
| 设计评审上因"不是应该ε=α吗"而争执不下 | 把光谱定律与全量平均混为一谈 | 用本文白漆屋顶的例子说明:光谱上恒成立,全量则取决于权重 |
记忆方法的终极版
说到底,基尔霍夫定律该怎么记,才能在实务中不出错?
"同一波长下永远ε=α;全量的等号,要先看清'对方的光谱'"——只要这两句就够了。落到操作上就是:开始做辐射计算之前,先数一数"这个问题里出现的辐射光谱有几种"。若只有自身温度的红外,就用灰体近似取单一ε;一旦太阳或高温源进来,立刻分波段。养成数光谱种类的习惯——这就是在实务中正确使用基尔霍夫定律的全部。
帮助
更多
错误