📚 课程目录第2章热力学基础回顾

2. 热力学基础回顾：热传导、热对流、热辐射的基本原理与数学描述

做光机设计的人，往往对力学和光学比较熟，但一碰到热就头疼。我当年刚入行时也是这样，觉得热分析不就是算个温度嘛，后来被现实狠狠教育了一顿——温度分布算不准，镜面变形就全错，光学性能直接崩掉。

所以这一章，咱们把热力学的三大基本传热方式捋一遍。别嫌基础，这些是热光学耦合仿真的根基。

2.1 热传导：固体里的热量搬运工

热传导，说白了就是热量在物体内部从高温区往低温区跑。固体里主要靠晶格振动和自由电子来传热。你想想看，一个镜筒一端被激光加热，另一端摸上去还是凉的，中间那段就是热传导在干活。

2.1.1 傅里叶定律

描述热传导的核心公式就是傅里叶定律：

q = -k · ∇T

其中：

q — 热流密度 (W/m²)，单位面积上流过的热量
k — 导热系数 (W/(m·K))，材料传热能力的度量
∇T — 温度梯度 (K/m)，温度变化的陡峭程度

负号表示热量从高温流向低温，这个方向别搞反了。

我的经验： 做光机仿真时，导热系数k值一定要查准。我曾经用了一个供应商给的典型值，结果仿真和实测差了8°C。后来发现实际材料的k值因为杂质含量不同，比标称值低了15%。从那以后，我习惯在项目初期就要求材料供应商提供实测报告，或者自己拿小块样品去测一下。

2.1.2 热传导方程

实际工程中，我们更关心温度随时间和空间的变化。这就用到热传导方程：

ρ · c_p · ∂T/∂t = ∇ · (k · ∇T) + Q

各项含义：

符号	含义	单位	典型值（铝合金）
ρ	密度	kg/m³	2700
c_p	比热容	J/(kg·K)	900
k	导热系数	W/(m·K)	167
Q	内热源	W/m³	—

左边是温度变化率，右边是热扩散加上内部产热。稳态时左边为0，方程就简化成拉普拉斯方程。

注意： 在光机系统中，内热源Q往往来自激光吸收、电机发热或电子元件。我见过有人把激光加热当成边界条件直接加在表面上，结果忽略了镜体内部的吸收，算出来的温度分布完全不对。正确的做法是把吸收的激光功率按体积分布折算成Q。

2.2 热对流：流体带走的热量

热对流发生在固体表面和流动的流体之间。光机系统里最常见的就是空气自然对流和强迫风冷。你想想看，一个光学平台上的镜架，周围空气在流动，表面热量就被带走了。

2.2.1 牛顿冷却定律

描述对流传热的核心公式：

q = h · (T_s - T_∞)

其中：

h — 对流传热系数 (W/(m²·K))
T_s — 固体表面温度 (K)
T_∞ — 流体主体温度 (K)

h值的大小差别很大，我整理了一个常用范围：

对流类型	h 范围 (W/(m²·K))	典型应用场景
自然对流（空气）	2 ~ 25	静止环境中的镜筒、支架
强迫对流（空气）	25 ~ 250	风扇冷却的电子设备
强迫对流（水）	100 ~ 15000	水冷镜座、激光晶体冷却
沸腾传热	2500 ~ 100000	高热流密度散热（极少用于光机）

  关键点： 自然对流的h值很难精确计算，因为它依赖于表面温度本身（温度影响空气密度，密度影响流动）。我一般先用经验公式估算，再通过实验标定。在热光学耦合仿真中，如果对流边界条件不准，温度场就会偏，镜面变形也跟着偏，最后光学性能分析就白做了。

2.2.2 对流换热系数的工程估算

实际工程中，我们常用无量纲关联式来算h。对于垂直平板的自然对流：

Nu = 0.59 · (Gr · Pr)^1/4   (层流)
Nu = 0.10 · (Gr · Pr)^1/3   (湍流)

其中：
Nu = h · L / k_f   (努塞尔数)
Gr = g · β · ΔT · L³ / ν²   (格拉晓夫数)
Pr = ν / α   (普朗特数)

看着公式多，其实现在仿真软件里都内置了这些关联式。但我建议你至少理解背后的物理意义——Gr数代表浮升力与粘性力的比值，Pr数代表动量扩散与热扩散的比值。理解了这些，你才能判断软件算出来的结果靠不靠谱。

避坑指南： 我曾经在一个项目中，用软件默认的自然对流关联式算一个水平放置的镜筒，结果温度分布明显不对称。后来发现，软件默认用的是垂直平板的关联式，而水平圆柱的关联式完全不同。改过来之后，结果就合理了。所以，一定要确认你用的关联式适用于你的几何和流动状态。

2.3 热辐射：不需要介质的传热

热辐射是唯一不需要介质的传热方式。在真空中工作的光学系统，比如太空望远镜，辐射是主要的散热途径。即使在地面系统中，高温部件（比如激光晶体）的辐射散热也不能忽略。

2.3.1 斯特藩-玻尔兹曼定律

一个黑体表面辐射的总功率：

E = σ · T⁴

其中σ = 5.67 × 10⁻⁸ W/(m²·K⁴)，是斯特藩-玻尔兹曼常数。注意温度是四次方，所以温度一高，辐射功率增长极快。

实际物体不是黑体，需要乘以发射率ε：

E = ε · σ · T⁴

2.3.2 辐射换热的基本公式

两个表面之间的净辐射换热量：

Q₁₂ = ε_eff · σ · A₁ · F₁₂ · (T₁⁴ - T₂⁴)

其中：

ε_eff — 有效发射率，取决于两个表面的发射率和几何关系
F₁₂ — 角系数，表示表面1发出的辐射到达表面2的比例
A₁ — 表面1的面积

注意： 角系数的计算是辐射换热中最容易出错的地方。在光机系统中，镜面、镜筒、遮光罩之间往往存在复杂的遮挡关系。我见过有人用软件默认的"对所有表面辐射"设置，结果算出来的温度比实测低了20%，因为忽略了镜筒内壁之间的多次反射。正确的做法是仔细定义辐射视角，必要时做光线追踪。

2.3.3 材料表面特性对辐射的影响

不同材料表面的发射率差别很大：

材料/表面处理	发射率 ε (常温)	备注
抛光铝	0.04 ~ 0.06	低发射率，反射率高
阳极氧化铝（黑）	0.80 ~ 0.90	高发射率，常用作散热面
光学玻璃	0.85 ~ 0.95	对红外波段不透明
镀金表面	0.02 ~ 0.03	极低发射率，用于热控
黑漆	0.85 ~ 0.95	标准高发射率涂层

我的习惯： 在做热光学耦合仿真时，我会把辐射边界条件单独列一个检查清单：①所有表面的发射率是否赋值？②角系数计算是否考虑了遮挡？③是否开启了辐射与对流的耦合？④对于透明介质（如光学窗口），是否考虑了辐射穿透？这个清单帮我避免过很多低级错误。

2.4 三种传热方式的对比与耦合

在实际光机系统中，三种传热方式往往是同时存在的。比如一个被激光加热的镜片：

镜片内部是热传导
镜片表面与周围空气是热对流
镜片表面与周围结构是热辐射

做仿真时，不能只考虑其中一种。我见过有人只算传导，忽略了对流和辐射，结果温度比实测高了30%。反过来，也有人过度强调辐射，忽略了自然对流在低温差下的主导作用。

  判断准则： 在常温附近（300K左右），自然对流的换热能力通常比辐射强。但当温度超过400K时，辐射的T⁴效应开始显现，辐射换热量会迅速超过自然对流。在真空环境中，辐射是唯一的散热途径，必须精确建模。

下面这张图总结了三种传热方式的核心要素和数学描述：

做热光学耦合仿真时，我的做法是：先判断系统的主要传热路径，然后对次要路径做简化处理。比如一个在地面实验室测试的光学系统，如果温差不大（<20°C），辐射可以先用等效对流系数近似，等初步结果出来后再决定要不要精确建模。这样能节省大量计算时间。