微流体力学 — 故障排除指南
更详尽的内容请访问 microfluidics.html。
微流体计算的典型故障
请教我微流体CFD中需要特别注意的故障。
与宏观尺度的CFD相比，存在许多不同的故障。让我们按模式进行分析。
1. 寄生电流（Spurious Current）导致液滴崩溃
症状：静止的液滴发生变形和分裂。界面附近产生非物理涡。
原因：表面张力的CSF模型中曲率计算的离散化误差。在Ca << 1的情况下最明显。
定量指标：当寄生电流速度超过物理流速的10%时，结果不可信。
对策：
- 提高界面附近的网格分辨率（每液滴直径最少10个单元）
- 切换到Height Function方法进行曲率计算
- 使用Phase-Field模型（Cahn-Hilliard方程）
- 如果使用OpenFOAM的
interFoam，调整interfaceCompression设置
应该在单相流测试中确认寄生电流的大小，然后再进行主要计算。
2. 接触线被固定，液滴不移动
症状：壁面上的液滴有驱动力但不移动
原因：no-slip条件与界面追踪的矛盾。接触线处速度为零，但界面移动的矛盾（Huh-Scriven悖论）产生数值问题。
对策：
- Navier-slip条件：在壁面引入微小滑移长度 $\beta$。$u_{\text{slip}} = \beta \partial u/\partial n$
- 滑移长度 $\beta = 1\text{--}10\,\text{nm}$ 是物理上正确的，但如果网格无法解析，使用数值较大的值
- Phase-Field模型能自然回避接触线特异性
3. 时间步长变得极小
症状：Courant数限制导致 $\Delta t \sim 10^{-10}\,\text{s}$ 量级，计算进展缓慢
原因：表面张力引起的毛细波传播速度限制。Brackbill限制：
$\Delta x = 1\,\mu\text{m}$，水的物性时，$\Delta t < 1.2 \times 10^{-8}\,\text{s}$。
对策：
- 不要过度细化网格
- 表面张力的隐式处理（某些求解器支持）
- 使用AMR在界面附近细化，远处保持粗
表面张力成为时间步长约束是微流体特有的问题。
4. 两相流质量守恒性差
症状：液滴体积在计算过程中减少（或增加）
原因：Level-Set方法将界面追踪为 $\phi = 0$ 等值面，但再初始化步骤中损失质量。
对策：
- 使用VOF方法（质量守恒性好）
- CLSVOF（Level-Set和VOF的混合）
- 添加质量修正步骤
- Phase-Field模型（Cahn-Hilliard：质量守恒型）
调试步骤
微流体计算的调试应按以下顺序进行。
1. 单相流测试：首先用单相Poiseuille流验证压力损失是否等于 $\Delta p = 12\mu Q L/(w h^3)$（矩形截面）
2. 静止液滴测试：放置球形液滴，验证Laplace压力 $\Delta p = 2\sigma/R$ 是否正确再现。确认寄生电流是否足够小
3. 接触角测试：验证壁面上液滴的平衡接触角是否与设置值一致
4. 动态测试：进行T形汇流处液滴生成等与实验的定量比较
逐步进行分层测试是微流体分析的基本要素。
特别是Laplace压力测试是必须进行的最少实验。如果这不匹配，两相流结果完全不可信。
微流体分析的"数值发散"——网格尺寸是90%的罪魁祸首
微流体分析中解不收敛时，首先怀疑网格的纵横比。通道截面为10μm×100μm时，使用纵横比1:10的单元会使压力梯度数值误差容易膨胀。理想情况是各方向至少5个单元，纵横比3以下。另一个故障源是入口边界条件的"急剧速度变化"，如果不假设充分发展流而设置均匀流速（塞流），入口附近会产生振荡并难以收敛。仅在入口前设置发展长度 L≈0.06·Re·D（微流体中约数十μm）就能改善。
相关主题
详解
错误