FR38 · 热输运

薄膜里的弹道与扩散导热

从一侧加热一层薄膜，热是像墨水一样慢慢扩散过去，还是直接飞到对面的壁？ 厚膜里，声子在穿过的路上被散射很多次，热慢慢扩散过去，傅里叶定律成立。膜厚接近一个平均自由程时，许多声子一口气直接飞到对面，冷侧很早就开始升温，但稳定后的热流反而比傅里叶定律预测的小。膜再薄，热流也超不过由声子速度和热容决定的弹道极限。

从一侧加热类硅薄膜，看声子穿过它。薄膜导的热比傅里叶定律预测的少，厚膜又回到傅里叶定律。

克努森数：平均自由程 / 厚度 $0.943$

此刻冷壁流出的热量 / 弹道极限 $0.276$

稳态薄膜热导率 / 体材料值 $0.429$

热壁开启后的时间 $20.0\,\mathrm{ps}$

设置薄膜

薄膜厚度（nm）

40

L

声子平均自由程（nm）

37.7

\Lambda

拖动蓝色冷壁可改变厚度。滑块按对数刻度变化。

求解器分辨率、声子包与种子

求解器每个半球的方向数

12

稳态薄膜中的蒙特卡罗声子包数

6000

蒙特卡罗的随机种子

3807

模型、单位与来源

$\frac{\partial e}{\partial t}+\mu\frac{\partial e}{\partial x}=\frac{\theta-e}{\mathrm{Kn}},\qquad \frac{q}{q_B}=2\int_{-1}^{1}\mu\,e\,d\mu,\qquad q_B=\frac{Cv\,\Delta T}{4}$

长度以膜厚为单位，时间以穿越膜厚的飞行时间为单位。求解器用离散方向、带限制器的迎风格式和精确弛豫；蒙特卡罗用同一模型，声子包在指数分布的自由时间后散射。绿色稳态值来自另一个积分方程，不使用任何一个模拟。

Ioffe NSM, Si thermal properties · Ioffe NSM, Si mechanical properties · Maznev, Johnson and Nelson (2011)

教程章节：薄膜里的弹道与扩散导热

模型假设

灰体声子玻尔兹曼方程，弛豫时间近似：所有声子用同一个速度和同一个平均自由程，散射各向同性，两侧为黑壁，吸收所有到达的声子并按自身温度发射。约 300 K 下硅的数据取自 Ioffe NSM 数据表；单一速度（[100] 方向声速的德拜平均）和由动理学公式反推的单一平均自由程是模型选择。真实硅的平均自由程分布很宽，所以灰体模型会低估薄膜要多厚才回到傅里叶定律。温度是相对冷壁的小偏差，方程保持线性。每个声子包的竖直位置只为让舞台易读而画出。

此刻膜内的温度分布

黄色：求解器。粉色：蒙特卡罗分层统计（误差线为一倍标准误差）。蓝色虚线：傅里叶。紫色点线：自由飞行。

冷壁流出的热量

黄色：求解器。粉色：蒙特卡罗。蓝色虚线：傅里叶。紫色点线：自由飞行。绿色虚线：精确稳态值。

稳态热导与克努森数

绿色：灰体精确解（圆点为当前薄膜）。黄色圆圈：求解器当前值。蓝色虚线：傅里叶。紫色点线：弹道极限。

更多数值与检验方法

薄膜厚度: $0$
平均自由程: $0$
穿越薄膜的飞行时间: $0$
稳态热流 / 弹道极限（积分方程）: $0$
傅里叶定律给出的同一比值: $0$
此刻热壁流入的热量 / 弹道极限: $0$
储存能量减去净流入能量（求解器检验）: $0$
膜内蒙特卡罗声子包数: $0$
求解器的网格数: $0$
求解器方向数: $0$

求解器和蒙特卡罗共用一个模型，但不共用代码：两者在误差线内一致，说明没有只出现在其中一方的计算错误。稳态值来自单独求解的积分方程，求解器在一次运行结束时逼近它。

物理教程

飞行与扩散的赛跑

背景在绝缘晶体里，热由声子携带。声子是晶格振动的波包，以接近声速飞行，直到被散射。两次散射之间平均走一个平均自由程；按单一速度估算，室温硅的平均自由程是几十纳米。

为什么值得理解傅里叶定律默认声子在穿过材料的路上要被散射很多次。如今芯片里的硅层只有几个平均自由程厚，这个前提不再成立，薄膜导的热比傅里叶定律给的少。

先抓重点

核心问题: 从一侧加热一层薄膜，热是像墨水一样慢慢扩散过去，还是直接飞到对面的壁？
一句话直觉: 决定结果的是一个数：克努森数 $\mathrm{Kn}=\Lambda/L$ 。它小时，热慢慢到达，傅里叶定律成立；它大时，声子一次飞过，对面的壁很早就开始升温，稳态热流被弹道极限 $q_B=Cv\,\Delta T/4$ 封顶，薄膜再薄也不会超过它。

核心数学模型

灰体玻尔兹曼方程

\frac{\partial e}{\partial t}+v\mu\frac{\partial e}{\partial x}=\frac{\theta-e}{\tau},\qquad \theta=\frac12\int_{-1}^{1}e\,d\mu

这里 $e$ 是沿方向余弦 $\mu$ 携带的能量，散射在平均时间 $\tau=\Lambda/v$ 后把它拉回局部温度 $\theta$ 。

傅里叶与弹道的对比

\frac{q_F}{q_B}=\frac{\kappa\,\Delta T/L}{Cv\,\Delta T/4}=\frac{4}{3}\,\mathrm{Kn}

代入 $\kappa=Cv\Lambda/3$ ，当 $\mathrm{Kn}>3/4$ 时傅里叶定律给出的热流会超过弹道极限，而灰体模型的精确解永远不会。

常见难点

不是越薄导热越好

容易误解膜越薄路程越短，所以导热应该越来越快。

正确理解只有声子还在膜内散射时，厚度减半才会让傅里叶热流翻倍。一旦声子一口气飞过，热流就不再增长，而是饱和在弹道极限；薄膜就像两面互相辐射的壁之间的空隙。

动手实验

01
看这场赛跑
从“厚度接近平均自由程”开始，看冷壁热流图中的黄色曲线。
你会看到： 热流比傅里叶虚线更早升起，最后稳定在明显低于傅里叶值的位置；粉色蒙特卡罗点沿着同一条曲线。
02
先变薄，再变厚
选择“很薄：声子直接飞过”，再选“厚膜：回到傅里叶定律”。
你会看到： 薄膜沿自由飞行曲线走，稳态接近弹道极限；厚膜沿傅里叶曲线走，稳态点落在傅里叶直线上。
03
换材料，不换尺寸
选择“同一薄膜，自由程很短”。
你会看到： 同样 40 nm 的薄膜现在几乎是扩散的：起作用的是平均自由程与厚度之比，而不是厚度本身。

薄膜里的弹道与扩散导热

飞行与扩散的赛跑

核心数学模型

灰体玻尔兹曼方程

傅里叶与弹道的对比

常见难点

不是越薄导热越好

动手实验

看这场赛跑

先变薄，再变厚

换材料，不换尺寸

相关实验与文章

从声子飞行到第二声、界面和热电