4. 流体力学基础：连续介质假设、纳维-斯托克斯方程、层流与湍流、雷诺数计算
在CVD反应腔室仿真中，流体力学是描述气体输运与化学反应耦合过程的底层理论。本章聚焦于四个核心概念：连续介质假设（建立宏观控制方程的前提）、纳维-斯托克斯方程（动量守恒的数学表达）、层流与湍流（两种基本流态）以及雷诺数（流态判据）。掌握这些基础，是后续构建腔室CFD模型、预测薄膜均匀性的关键。
核心观点： 绝大多数CVD反应腔室（如LPCVD、PECVD）在工艺条件下处于低雷诺数层流状态，但局部几何突变（如喷淋头孔口、晶圆边缘）可能诱发转捩。准确判断流态是选择湍流模型或层流模型的前提。
4.1 连续介质假设
流体力学将流体视为由连续分布的质点组成的介质，忽略分子间隙和个体分子运动，从而用连续函数（如密度、速度、压力）描述流场。该假设成立的条件是：特征长度远大于分子平均自由程。
| 参数 | 定义 | CVD典型值 |
|---|---|---|
| 分子平均自由程 λ | 分子两次碰撞间平均移动距离 | 0.1 ~ 10 μm（低压下） |
| 特征长度 L | 腔室最小几何尺寸（如喷淋头间隙） | 5 ~ 50 mm |
| 克努森数 Kn = λ / L | 判断连续介质是否适用 | Kn < 0.01（连续流） |
仿真提示： 当 Kn > 0.1 时（如超高真空CVD），连续介质假设失效，需改用DSMC（直接模拟蒙特卡洛）方法。常规CVD工艺（压力 > 100 Pa）通常满足连续流条件。
4.2 纳维-斯托克斯方程
纳维-斯托克斯（N-S）方程是描述粘性流体动量守恒的偏微分方程组。在CVD仿真中，通常求解不可压缩或低马赫数可压缩形式的N-S方程，并与能量方程、组分输运方程耦合。
直角坐标系下不可压缩N-S方程（矢量形式）：
ρ (∂u/∂t + u·∇u) = -∇p + μ ∇²u + ρg其中：
- ρ — 流体密度（kg/m³）
- u — 速度矢量（m/s）
- p — 压力（Pa）
- μ — 动力粘度（Pa·s）
- g — 重力加速度（m/s²）
常见误区： 初学者常忽略N-S方程中的对流项 u·∇u。在CVD腔室中，即使入口速度很低（如0.1 m/s），对流项在喷淋头孔口附近仍可能主导动量输运，忽略它将导致流场失真。
4.3 层流与湍流
根据流动状态，CVD腔室内的气流可分为：
| 流态 | 特征 | CVD典型场景 | 仿真模型 |
|---|---|---|---|
| 层流 | 流线平行、无横向掺混、速度分布呈抛物线 | LPCVD炉管、低压平板反应器 | 直接求解N-S方程（层流模型） |
| 湍流 | 随机涡旋、强扩散、速度脉动 | 高流速喷淋头、大尺寸腔室 | RANS（k-ε、k-ω SST）或LES |
| 转捩 | 层流到湍流的过渡区域 | 晶圆边缘、台阶处 | 转捩模型（γ-Reθ） |
工程经验： 对于CVD腔室，当雷诺数 Re < 2000 时，可安全使用层流模型；当 Re > 4000 时，必须使用湍流模型。2000 < Re < 4000 为转捩区，建议使用转捩模型或进行网格敏感性验证。
4.4 雷诺数计算
雷诺数（Reynolds Number, Re）是判断流态的无量纲数，定义为惯性力与粘性力之比：
Re = ρ · U · L / μ其中：
- U — 特征速度（m/s），通常取入口平均速度
- L — 特征长度（m），对于圆管取直径，对于平板取入口到晶圆距离
- ρ — 流体密度（kg/m³）
- μ — 动力粘度（Pa·s）
CVD腔室雷诺数计算示例：
| 参数 | 符号 | 数值 | 单位 |
|---|---|---|---|
| 气体密度（N₂，400 K） | ρ | 0.85 | kg/m³ |
| 入口速度 | U | 0.5 | m/s |
| 喷淋头间隙（特征长度） | L | 0.02 | m |
| 动力粘度（N₂，400 K） | μ | 2.2 × 10⁻⁵ | Pa·s |
| 雷诺数 | Re | 386 | — |
结论： 该算例 Re = 386 < 2000，属于层流。在CVD仿真中，可直接使用层流模型，无需引入湍流模型，从而节省计算资源并避免湍流模型引入的额外误差。
4.5 本章小结
- 连续介质假设是CFD仿真的前提，通过Kn数验证其适用性。
- N-S方程是动量守恒的核心，需与连续性、能量、组分方程耦合求解。
- 层流与湍流的区分直接影响仿真模型选择，CVD腔室多为层流。
- 雷诺数是流态判据，计算公式为 Re = ρUL/μ，需根据实际工艺参数计算。