2、传热学基础:热传导、热对流、热辐射的基本概念与数学描述

各位好,我是老张。在半导体设备这行摸爬滚打了十几年,今天咱们来聊聊传热学。你可能会问:“搞CVD仿真,为啥非得啃传热学?” 说白了,CVD腔体里那点事儿,本质上就是个热场游戏。晶圆表面的温度均匀性,直接决定了薄膜的厚度和成分。我见过太多因为热场没调好,整批晶圆报废的案例。所以,打好这个基础,比啥都重要。

传热有三种基本方式:热传导、热对流、热辐射。在CVD腔体里,这三种方式往往同时存在,互相纠缠。咱们一个一个拆开看。

2.1 热传导:固体里的热量传递

热传导,就是热量从高温区往低温区跑,靠的是分子、原子之间的碰撞和振动。在固体里,这是最主要的传热方式。比如晶圆本身、加热器、腔体壁,热量都是这么传的。

数学描述:傅里叶定律

这个定律是热传导的基石。公式很简单:

q = -k * (dT/dx)

其中:

  • q:热流密度,单位是 W/m²。说白了就是单位面积上,每秒钟流过多少热量。
  • k:导热系数,单位是 W/(m·K)。这是材料本身的属性。铜的导热系数高,陶瓷的导热系数低。我在选加热器材料时,特别看重这个值。
  • dT/dx:温度梯度。温度变化越剧烈,传热越快。

避坑指南: 我曾经在仿真一个SiC加热器时,忽略了导热系数随温度的变化。结果仿真结果和实测差了20%。记住,k值不是常数,尤其是高温下,变化很大。一定要查材料手册,或者用实验数据拟合。

对于三维问题,傅里叶定律可以写成更一般的形式:

q = -k * ∇T

这里的∇T是温度梯度向量。嗯,看着复杂,但核心思想没变:热量沿着温度下降最快的方向跑。

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

热对流,是流体(气体或液体)流动时带走热量的过程。在CVD腔体里,反应气体流过晶圆表面,一边反应,一边带走热量。这直接影响晶圆表面的温度分布。

数学描述:牛顿冷却定律

这个定律描述的是流体和固体表面之间的换热:

q = h * (T_s - T_f)

其中:

  • q:对流换热热流密度,W/m²。
  • h:对流换热系数,W/(m²·K)。这个系数很关键,它受流体速度、物性、表面形状等多种因素影响。说白了,h值越大,换热越猛。
  • T_s:固体表面温度。
  • T_f:流体主流温度。

个人经验: 我建议你在做CVD仿真时,对流换热系数h不要随便拍脑袋。最好通过CFD(计算流体动力学)仿真或者实验关联式来估算。比如,对于平板层流,h ≈ 0.664 * (k/L) * Re^(1/2) * Pr^(1/3)。Re是雷诺数,Pr是普朗特数。这些参数决定了流体的流动状态和换热能力。

你想想看,如果气体流速太快,晶圆表面热量被迅速带走,温度就上不去。如果流速太慢,反应气体供应不足,沉积速率又不够。所以,找到那个平衡点,就是咱们工程师的活儿。

2.3 热辐射:不需要介质的能量传递

热辐射,是物体通过电磁波向外发射能量的过程。在高温CVD腔体里(比如800°C以上),热辐射是主要的传热方式。它不需要介质,在真空中也能传热。加热器、晶圆、腔体壁之间,都在互相辐射。

数学描述:斯特藩-玻尔兹曼定律

这个定律描述了黑体(理想辐射体)的辐射能力:

E = σ * T^4

其中:

  • E:辐射力,单位是 W/m²。
  • σ:斯特藩-玻尔兹曼常数,5.67×10⁻⁸ W/(m²·K⁴)。
  • T:绝对温度,单位是K。

注意,这里是T的四次方!温度稍微一高,辐射能量就急剧增加。这也是为什么高温下辐射占主导的原因。

实际物体不是黑体,所以需要引入发射率ε:

E = ε * σ * T^4

ε在0到1之间。抛光金属的ε很小(比如0.05),而粗糙表面的ε接近1。我在仿真中,经常需要查材料在不同温度下的发射率。这个参数对结果影响很大。

警告: 热辐射是高度非线性的。因为它是T的四次方,所以温度场计算时,辐射项会让方程变得很难收敛。我曾经在调试一个辐射模型时,迭代了上千步才稳定下来。建议你使用合适的辐射模型(比如离散坐标法或蒙特卡洛法),并注意网格质量。

2.4 三种传热方式的耦合

在CVD腔体里,这三种传热方式不是孤立的。举个例子:加热器通过热辐射把热量传给晶圆,晶圆内部通过热传导把热量均匀分布,同时反应气体流过晶圆表面,通过对流带走一部分热量。这三者相互影响,共同决定了晶圆表面的温度分布。

为了让你更直观地理解,我画了一张图:

CVD腔体传热方式示意图 腔体壁(热传导 + 热辐射) 加热器(热辐射源) 晶圆(热传导) 反应气体(热对流) 热辐射 热辐射 热传导 热对流

从这张图里,你可以看到:加热器通过热辐射加热晶圆,晶圆内部通过热传导把热量传开,同时反应气体流过晶圆表面,通过对流带走热量。腔体壁也在和晶圆、加热器进行辐射换热。这就是一个典型的耦合传热问题。

2.5 传热方程:统一描述

在仿真中,我们通常用能量守恒方程来统一描述这三种传热方式。对于固体区域,方程是:

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

其中:

  • ρ:密度
  • Cp:比热容
  • ∂T/∂t:温度随时间的变化(瞬态项)
  • ∇ · (k ∇T):热传导项
  • Q:内热源(比如加热器的焦耳热)

对于流体区域,还需要加上对流项:

ρ * Cp * (∂T/∂t + u · ∇T) = ∇ · (k ∇T) + Q

这里的u是流体速度矢量。u · ∇T就是对流项,它描述了流体流动对温度场的影响。

小技巧: 在CVD仿真中,我通常先做流场计算(求解Navier-Stokes方程),得到速度场u,然后再代入能量方程求解温度场。这样能避免同时求解带来的收敛困难。当然,如果你用商业软件,这些都已经封装好了,但理解背后的逻辑,能帮你更好地设置边界条件和判断结果是否合理。

好了,传热学的基础就聊到这儿。这三种传热方式,你搞清楚了,CVD热场仿真的大门就算迈进去了。下一节,咱们会把这些知识用到具体的CVD腔体建模中,看看实际怎么操作。


专注资料整理