3、流体力学基础：层流与湍流、连续性方程、纳维-斯托克斯方程简介
各位工程师朋友，大家好。今天我们来聊聊CVD热场仿真里最绕不开的一个话题——流体力学基础。
说实话，我刚开始做CVD仿真那会儿，总觉得流体力学是搞CFD的人的事。我们做热场的，盯着温度场不就行了？结果第一次做垂直流反应器仿真，出来的沉积均匀性一塌糊涂。后来才发现，问题出在气流模式上——我压根没搞懂气体在腔体里是怎么流动的。
所以这一节，咱们把流体力学里最核心的三个概念捋一遍。不搞复杂的推导，只讲你仿真时真正用得上的东西。
3.1 层流与湍流：两种截然不同的流动状态
先问一个问题：CVD腔体里的气体，到底是安静地流过晶圆表面，还是像暴风雨一样乱窜？
答案是：绝大多数CVD工艺都工作在层流区。为什么？因为层流状态下，气体分子以平行层的方式运动，传质均匀，沉积膜厚才能一致。
我见过一个案例，某团队做MOCVD仿真，为了加快计算速度，把进气速度调高了3倍。结果仿真出来的流场完全乱了——从层流变成了湍流，晶圆边缘的沉积速率比中心高了40%。这就是典型的「为了省时间反而浪费时间」。
判断层流还是湍流，看一个无量纲数就够了：雷诺数（Re）。
雷诺数公式：
Re = ρ * v * L / μ其中：ρ = 气体密度 (kg/m³)，v = 流速 (m/s)，L = 特征长度 (m)，μ = 动力粘度 (Pa·s)
判据：Re < 2300 → 层流；2300 < Re < 4000 → 过渡区；Re > 4000 → 湍流
在CVD腔体里，特征长度L通常取反应腔的直径或晶圆直径。我习惯先估算一下Re数——如果Re < 100，那基本可以放心用层流模型。如果Re接近2000，就得小心了，腔体结构稍微有点变化就可能触发转捩。
个人经验：做CVD热场仿真时，我建议始终先按层流跑一遍。如果收敛困难或者结果明显不合理，再检查Re数是否进入了过渡区。不要一上来就开湍流模型——计算量翻倍不说，很多湍流模型对CVD这种低Re数流动其实并不适用。
3.2 连续性方程：质量守恒的数学表达
说白了，连续性方程就是一句话：流进去多少质量，就得流出来多少质量。不能凭空产生，也不能凭空消失。
在CVD腔体里，气体从进气口进入，从排气口排出。如果进气流量是10 slm（标准升每分钟），排气流量也必须是10 slm——除非有化学反应消耗了气体。但即使有反应，总质量还是守恒的，只是组分变了。
连续性方程的微分形式长这样：
∂ρ/∂t + ∇·(ρu) = 0其中：
- ∂ρ/∂t：密度随时间的变化率（稳态仿真时这一项为0）
- ∇·(ρu)：质量通量的散度，说白了就是「净流出量」
稳态情况下，方程简化为：
∇·(ρu) = 0嗯，这里要注意：很多初学者在做CVD仿真时，会忽略密度变化。但CVD腔体里温度变化很大——从室温到1000°C以上，气体密度能差好几倍。如果你还用不可压缩假设，算出来的速度场会严重失真。
我记得有一次帮客户排查问题，他们做的SiC CVD仿真，晶圆表面温度1200°C，进气口温度只有25°C。他们用了不可压缩模型，结果算出来的气流速度在高温区偏高了3倍多。后来改成可压缩模型，结果才跟实验数据对得上。
3.3 纳维-斯托克斯方程：流体力学的核心
纳维-斯托克斯方程（简称N-S方程），是描述流体运动的基本方程。它本质上就是牛顿第二定律在流体上的应用——力 = 质量 × 加速度。
对于CVD热场仿真，我们关心的是N-S方程在直角坐标系下的形式：
ρ(∂u/∂t + u·∇u) = -∇p + μ∇²u + ρg + F每一项的含义：
- ρ(∂u/∂t + u·∇u)：惯性项，描述流体的加速度
- -∇p：压力梯度力，气体从高压区流向低压区
- μ∇²u：粘性力，让流体「粘」在壁面上
- ρg：重力，在水平流CVD中影响不大，但在垂直流中不可忽略
- F：体积力，比如热浮力（自然对流）
避坑指南：我曾经在仿真一个水平管式CVD时，忽略了重力项。结果算出来的流场完全对称，但实际工艺中晶圆上方的气体因为热浮力会向上飘，导致沉积不均匀。后来加上重力项，结果才跟实验吻合。所以——只要腔体里有温度梯度，就别忽略重力。
3.4 知识体系总览
下面这张图，是我自己总结的流体力学基础在CVD热场仿真中的应用框架。你可以把它当作一个「思维导图」来用：
3.5 仿真中的实际应用
好了，理论讲完了。咱们说说在仿真软件里怎么用这些东西。
我个人习惯的流程是这样的：
- 先算Re数——根据工艺条件（压力、温度、流量、腔体尺寸），估算流动状态。如果Re < 500，放心用层流模型。
- 设置边界条件——进气口给质量流量或体积流量，排气口给压力出口。注意：如果用了体积流量，一定要换算到实际工况下的密度。
- 开启能量方程——CVD必然涉及温度场，所以能量方程必须打开。N-S方程和能量方程是耦合的：温度影响密度，密度影响速度，速度又影响对流换热。
- 检查收敛性——连续性方程的残差至少要降到1e-4以下。如果残差震荡不收敛，先检查网格质量，再检查边界条件是否合理。
一个小技巧：做CVD热场仿真时，我习惯先跑一个「等温流场」——假设温度均匀，只算速度场和压力场。等这个收敛了，再打开能量方程，把温度场耦合进来。这样分步走，收敛性会好很多，也容易定位问题。
最后说一句：流体力学这东西，光看书是学不会的。你想想看，Re数、连续性方程、N-S方程，这些概念单独看都不难。但真正放到CVD腔体里，各种物理场耦合在一起，问题就复杂了。我的建议是：先搞懂物理意义，再上手仿真。别一上来就点「Solve」——先问问自己：我这个流动是层流还是湍流？质量守恒了吗？重力考虑了吗？
嗯，这一节就到这里。下一节我们聊聊传热学基础——热传导、对流和辐射，在CVD腔体里是怎么打架的。
公众号:蓝海资料掘金营，微信deep3321