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腔体里,这三种传热方式不是孤立的。举个例子:加热器通过热辐射把热量传给晶圆,晶圆内部通过热传导把热量均匀分布,同时反应气体流过晶圆表面,通过对流带走一部分热量。这三者相互影响,共同决定了晶圆表面的温度分布。
为了让你更直观地理解,我画了一张图:
从这张图里,你可以看到:加热器通过热辐射加热晶圆,晶圆内部通过热传导把热量传开,同时反应气体流过晶圆表面,通过对流带走热量。腔体壁也在和晶圆、加热器进行辐射换热。这就是一个典型的耦合传热问题。
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腔体建模中,看看实际怎么操作。