洁净室气流分析
洁净室气流的理论基础
概述
老师！洁净室的气流分析是在半导体工厂等地使用的，对吧？涉及哪些物理学呢？
洁净室气流分析是一种利用CFD来预测为了维持室内洁净度而采用的一方向流（单向流）或乱流置换方式等气流模式的技术。为了达到ISO 14644-1定义的洁净度等级（Class 1～Class 9），需要对FFU（风机滤网组件）吹出的气流如何搬运和排出微粒进行仿真。
也就是说，可以用数值方法验证洁净度等级，这样就清楚了。
支配方程
气流分析用的是Navier-Stokes方程吧？微粒追踪怎么做呢？
首先连续气相用基于RANS的Navier-Stokes方程求解。通常假设不可压缩流动。
连续方程和Navier-Stokes方程如下。
FFU的滤网部分用多孔介质模型表示。包含达西-福尔希海默阻力。
$\alpha$ 是透过率，$C_2$ 是惯性阻力系数吧。能从滤网目录值反算吗？
没错。HEPA滤网的典型值是面风速0.45 m/s时压力损失约250 Pa。从这个和滤网厚度可以推导出 $\alpha$ 和 $C_2$。微粒追踪用DPM（离散相模型）求解微粒运动方程。
还要加入布朗力？亚微米粒子的布朗运动确实很重要。
是的，0.1 μm以下的粒子布朗扩散占主导。还需要Cunningham修正系数 $C_c$。
乱流模型选择
洁净室用什么乱流模型比较合适？
一方向流洁净室存在层流和乱流并存的区域，推荐使用SST $k$-$\omega$ 模型。低雷诺数区域的壁面处理能自然进行。
| 乱流模型 | 推荐度 | 特征 |
|---|---|---|
| SST k-omega | 高 | 低雷诺数壁面处理、分离预测强 |
| Realizable k-epsilon | 中 | 通用性好但需要壁面函数 |
| RNG k-epsilon | 中 | 旋转流略优 |
| LES (Smagorinsky) | 非常高（计算成本大） | 直接解析涡结构 |
在半导体工厂实际案例中SST k-omega用得比较多是吧。LES是研究用途吗？
是的。但最近，人体动作引起的乱流非定常分析中LES也开始被使用。
实务注意事项
请教我现场要注意的要点。
人体发塵模型也要放入CFD，这是洁净室特有的话题，很有意思。
HEPA滤网"99.97%"数字的本质
学习洁净室气流理论时绕不过的是HEPA滤网的捕集效率。"0.3μm微粒捕集率99.97%以上"指的是最坏情况，也就是说这个尺寸最容易穿透。为什么是0.3μm？比这更小的超微粒子受到布朗运动影响，更容易被纤维捕获；比这更大的粒子由于惯性力容易被捕获。正好在0.3μm附近是"惯性小又不被布朗运动支配"的过渡领域，成为最难捕获的尺寸。进行洁净室CFD微粒追踪时，理解这个分布特性来设置粒径范围很重要。
洁净室气流的数值计算方法
数值方法详解
实际求解洁净室CFD时用有限体积法吧？具体用什么离散化格式？
洁净室气流是低马赫数不可压缩流，用基于压力的求解器。用SIMPLE系列算法（SIMPLE、SIMPLEC、PISO）耦合求解压力和速度。
压力速度耦合
SIMPLE和SIMPLEC怎么区分选择？
定常分析用SIMPLEC（不需要压力修正松弛因子，收敛快），非定常分析用PISO（每时间步迭代少）。Coupled Solver也是选项，但内存消耗大。
| 算法 | 定常/非定常 | 特征 |
|---|---|---|
| SIMPLE | 定常 | 基础方法，需要调整松弛系数 |
| SIMPLEC | 定常 | 收敛快，洁净室推荐 |
| PISO | 非定常 | 人体动作非定常分析用 |
| Coupled | 两者 | 健壮但内存2～3倍 |
空间离散化
对流项用什么格式比较好？
洁净室是低速流（0.3～0.5 m/s左右），数值扩散容易成为问题。推荐Second Order Upwind及以上。
QUICK格式不能用于四面体网格吧。
对。QUICK需要结构化网格或六面体网格。多面体网格用Second Order Upwind比较稳妥。
DPM实现细节
请教微粒追踪的具体设置。
DPM用时间积分追踪微粒轨迹。洁净室分析的典型设置如下。
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 粒径分布 | Rosin-Rammler (0.1～10 um) | ISO 14644-1目标粒径 |
| 粒子数 | 10,000以上/注射面 | 统计可信度 |
| 积分方法 | Trapezoidal | 精度和速度平衡 |
| 布朗力 | ON (dp < 1 um) | 亚微米粒子必须 |
| Saffman升力 | ON | 改善壁面附近行为 |
| 壁面条件 | Trap/Reflect | 沉积 vs. 反弹 |
粒子数1万以上，感觉挺多的。对计算时间影响大吗？
DPM是单向耦合的话，只需在气相计算后作为后处理追踪微粒，额外成本全体的10～20%左右。洁净室粒子浓度低，单向耦合就足够。
网格策略
洁净室是大空间，网格数目标多少？
典型半导体工厂一个工作间（10m × 20m × 3m）需要500万～2000万单元。FFU吹出面和晶圆周边需要局部细化，最小单元尺寸约5～10 mm。
地下回风室网格也要弄得相当细，会影响压力损失吗？
地下回风室的开口率约25%的栅板地产生很大压损。有时用多孔质跳跃条件简化，但如果局部偏流成问题需要完整建模。
洁净室CFD的关键是"保持层流"的方法选择
洁净室气流分析的数值方法中，低雷诺数乱流模型（Low-Re k-ε）和层流模型的区分使用是关键。FFU吹出直下雷诺数较低，接近层流，如果用乱流模型求解会产生过大扩散，导致微粒行为非现实。实务中的有效做法是：吹出口附近按层流对待，流向作业区时切换到乱流模型——这种分区切换法。结合拉格朗日微粒追踪（DPM），能定量评估各作业位置是否容易产生污染。
洁净室气流的实务应用
实践指南
老师，请教洁净室CFD的实务分析流程。
整体流程是(1)收集布局信息、(2)构建CAD模型、(3)生成网格、(4)设置边界条件、(5)求解、(6)评估洁净度。
分析流程
CAD模型要详细到什么程度？
1. 形状简化要点
- FFU作为多孔面，给定吹出风速（内部结构不需要）
- 制造装置简化为块状（排气口要精确放置）
- 人体用简易圆柱模型（直径0.4m、高度1.7m）＋发热量75 W
- 管道、电缆架可省略（闭塞度<5%情况）
2. 边界条件设置
- FFU吹出面：Velocity Inlet（0.45 m/s、均匀）＋温度（22℃ 典型）
- 回风口/栅板地：Pressure Outlet（表压0 Pa）
- 装置排气口：Velocity Inlet（负值）或Outflow
- 壁面：绝热No-Slip（天井、墙）、温度指定（发热装置表面）
FFU吹出温度直接用空调设计值就行吗？
基本上是这样，但CFD的目的之一是验证装置发热导致的局部温度上升，所以最好设置空调机出口温度（通常18～20℃）作为吹出温度，把室内温度分布作为计算结果得出。
洁净度评估方法
从CFD结果怎么判定ISO洁净度等级？
从DPM微粒追踪结果推导评估点的粒子浓度。具体步骤如下。
1. 从发塵源（人体、装置）释放微粒
2. 计数经过评估点（晶圆上方300mm等）的微粒数
3. 从通过粒子数和发塵率推导个数浓度 [个/m³]
4. 与ISO 14644-1表对照判定等级
| ISO等级 | 0.1 um [个/m³] | 0.5 um [个/m³] | 应用案例 |
|---|---|---|---|
| Class 1 | 10 | - | 最新半导体光刻 |
| Class 3 | 1,000 | 35 | 半导体前工序 |
| Class 5 | 100,000 | 3,520 | 半导体后工序 |
| Class 7 | - | 352,000 | 药品制造 |
有定量评估能力真棒。不过微粒数的统计误差让人担心。
对。必须确保充足的DPM微粒追踪数。通过将微粒数加倍来验证结果不变，这个"微粒数独立性验证"是必须的。
常见失败与对策
初学者容易犯什么错？
| 失败案例 | 原因 | 对策 |
|---|---|---|
| 面风速分布太不均 | FFU下方网格太粗 | FFU面至少保留10×10网格 |
| 温度在室内均匀 | 数值扩散过大 | 采用Second Order及以上 |
| 微粒评估不稳定 | DPM微粒数不足 | 每个源至少1万微粒 |
| 地板回风偏流 | 栅板开口率设置错 | 确认多孔质跳跃K值 |
| 非定常计算发散 | 时间步长过大 | CFL数 < 1为目标 |
数值扩散导致费力得到的温度分布变得光滑，真的很遗憾。离散化格式的选择很关键吧。
对。洁净室是低速流，一阶风上会由数值扩散导致温度分布几乎均匀。一定要用Second Order Upwind及以上。
"人走一步引起0.3 m/s乱流"——洁净室设计最大的敌人
洁净室气流分析实务中最容易被忽视的是"人体发热和动作产生的乱流"。研究显示，人走动时体周产生0.2～0.5 m/s的乱气流，会扰乱FFU吹出的0.4～0.5 m/s清净向下流。实务中的最好做法是在CFD模型中加入人体形状（热源33～35 W/人），按工作动线模拟污染扩散的"动态洁净室分析"已成为普遍做法。这对入退室的动线设计和防护服义务化的判断都有帮助。
洁净室气流的软件比较
商用工具比较
洁净室CFD用什么工具比较好？听说有专用软件。
除了通用CFD求解器，还有洁净室专用简化工具。各自特点整理如下。
工具列表
| 工具名 | 开发商 | 特征 | 洁净室适用性 |
|---|---|---|---|
| Ansys Fluent | ANSYS Inc. | DPM完整、多孔质完备 | 非常高 |
| Simcenter STAR-CCM+ | Siemens | 多面体网格、自动化 | 非常高 |
| Ansys CFX | ANSYS Inc. | 耦合求解器 | 高 |
| OpenFOAM | 开源 | icoFoam/simpleFoam + DPM | 高（需知识） |
| COMSOL Multiphysics | COMSOL AB | Particle Tracing模块 | 中（小规模） |
| FlowDesigner | Advanced Knowledge Institute | 建筑空调专化，日本产 | 高（空调设计向） |
| Stream | Cradle（MSC） | 建筑、空调CFD | 高（国内实绩丰富） |
FlowDesigner和Stream在日本洁净室设计中实绩多吗？
是的。FlowDesigner的UI针对空调设计人员，FFU和空调机设置很直观。Stream在国内建设公司广泛使用。但如果需要高度的DPM分析或多相流，Fluent和STAR-CCM+更灵活。
工具选择判断标准
选择标准是什么？
按用途选择很重要。
| 用途 | 推荐工具 | 理由 |
|---|---|---|
| 半导体工厂整体气流验证 | Fluent / STAR-CCM+ | 大规模DPM、HPC对应 |
| 空调设备设计检讨 | FlowDesigner / Stream | 设备设计人员向UI |
| 迷你环保间（局部清净空间）详细评估 | Fluent / CFX | 精密DPM ＋ Brownian |
| 研究、参数化研究 | OpenFOAM | 免费许可、自动化 |
| 简易布局检讨 | COMSOL | 多物理场耦合 |
Ansys Fluent设置示例
用Fluent的话具体怎么设置？
1. Viscous Model：SST k-omega、Low-Re Corrections ON
2. Energy Equation：ON（评估温度分布时）
3. Porous Media：FFU面设置Face Porous Jump
4. DPM：Injection → Surface injection（从发塵面）
5. Species Transport：必要时进行示踪气体扩散分析
FFU的Porous Jump设置值目安如下。
| 参数 | HEPA (H13) | ULPA (U15) |
|---|---|---|
| Face Permeability [m²] | 5.0e-10 | 2.0e-10 |
| Pressure Jump Coefficient [1/m] | 200 | 500 |
| Medium Thickness [m] | 0.065 | 0.065 |
滤网等级不同透过率差别这么大。从目录压损曲线反算吗？
对。从厂家目录的"面风速 vs. 压力损失"数据用最小二乘法分离出达西项（粘性阻力）和福尔希海默项（惯性阻力）。
专用工具还是通用CFD——工具选择的分水岭
洁净室气流分析工具选择时经常听到"专用软件还是通用CFD哪个好？"Flomerics（现Mentor Simcenter）的FloVENT和IDEAS/AirPack等专用工具有丰富的洁净室、数据中心特化设置，非专业人士也能快速完成分析。另一方面ANSYS Fluent和OpenFOAM等通用CFD虽自由度高，但边界条件的构建更耗时。实务中常见的做法是"初期设计检讨→专用工具"，"问题调查、详细最优化→通用CFD"这样的分工，单一工具完全解决通常反而低效。
洁净室气流的先端研究
先端课题与研究动向
洁净室CFD最新的趋势是什么？
随着半导体微细化要求的洁净度升高，传统定常RANS分析无法捕捉的现象成为问题。
LES非定常气流评估
听说在洁净室应用LES（大涡模拟）的研究增加了。
对。作业者动作时的发塵本质上是非定常现象，RANS的时间平均无法捕捉。LES用Smagorinsky或WALE（墙面自适应局部涡粘性）模型直接解析涡结构。
Smagorinsky常数 $C_s$ 多少合适？
室内气流一般 $C_s = 0.1$～$0.12$。Dynamic Smagorinsky模型的话自动调整，不用操心常数。
AMC（气载分子污染）分析
不仅粒子，还要评估分子级污染？
最尖端的EUV光刻中，有机分子（放气）沉积在掩膜或镜头上降低良率。这叫AMC分析。用物种输运方程求解有机气体的扩散对流。
$Y_i$ 是化学种质量分率，$D_i$ 是扩散系数吧。源项 $S_i$ 是放气发生率？
对。装置或密封材料的放气发生率 [ng/(cm² hr)] 作为源项。活性炭滤网的去除效率也要建模。
数字孪生与实时监控
洁净室的数字孪生已经实用化了？
部分半导体工厂已经运用粒子计数器和温湿度传感器的实时数据与CFD模型联动的数字孪生。
ROM怎么做？
用POD（正交分解）提取支配性模态，用少数的模态系数来近似气流场。从约100个CFD结果构建的ROM能在毫秒级预测FFU风量或装置发热变化时的气流。
洁净室气流的故障排除
故障排除
洁净室CFD计算不顺利时怎么处理？
按模式整理常见问题。
1. 不收敛、残差振荡
特别是定常计算时残差降不下来。
可能原因与对策：
- Boussinesq浮力与压力插值不匹配：有浮力时用PRESTO!或Body Force Weighted压力插值。Standard插值会残差振荡
- 多孔质跳跃设置值过极端：FFU阻力太大导致压力场不稳定。从小值开始逐步增大
- 低弛系数调整：压力改0.2、动量改0.5试试
压力插值用PRESTO!才行是Boussinesq近似特有的问题吧。
2. DPM微粒全部被壁面捕获
现象：微粒释放直后几乎全部被壁面捕获，评估点无法到达。
对策：
- 确认壁面DPM边界条件不是Trap（改为Reflect或Escape）
- 确认有没有打开布朗力（微小粒子需要ON）
- 确认微粒初速和气流速度一致（为0的话重力落下）
3. 温度分布在室内几乎均匀
装置发热已加上但温度差不出来的情况。
对策：
- 确认离散化格式（一阶风上导致数值扩散温度变平）
- 确认网格不太粗（装置附近最低20 mm以下单元）
- 确认装置发热量单位（W vs. W/m² 搞混）
4. 与实测偏差大
实测风速或微粒计数与CFD不符时怎么办？
确认项：
| 检查项 | 常见问题 | 处理 |
|---|---|---|
| FFU实际风量 | 目录值与实际差异 | 更新为实测值 |
| 泄漏 | 管道接口漏风 | 加入漏风量 |
| 装置排气量 | 实际运转状态不同 | 确认实测值 |
| 地板下回风室 | 模型简化影响 | 提高开口率、管道障碍物精度 |
| 人体模型 | 忽视动态发塵 | 非定常+DPM重新计算 |
FFU实际风量经常与目录值不同吗？
滤网目堵、管道压损、FFU风机老化等原因，实际面风速比新品时低10～20%。尽可能用实测值。
Fluent特有错误
Fluent特有的要点有吗？
Incomplete是微粒在室内一直漂浮的状态吧。加大最大步数或设置超时时间比较好。
对。设置Physical Time Limit后一定时间强制停止，这样计算时间更好预测。从洁净室换气次数推算滞留时间，设为那个3倍左右为目安。
"FFU压力不匀"引发的神秘不均流
洁净室CFD故障排除中经常遇到的是多个FFU间的流量偏差。各FFU独立有电动风机，相邻单元间的静压失衡时有些FFU会接近"逆流模式"。全FFU条件按均匀流量设置但实测不符时，要怀疑这种偏差。把每个FFU的静压—流量特性（P-Q曲线）加入模型做耦合分析，通常就能重现现场的不均流。
相关课题
有益
细些
错误