斯托克斯流（低雷诺数） — 故障排除指南
本文已移至综合版
更丰富的内容请访问 stokes-flow.html。
更丰富的内容请访问 stokes-flow.html。
斯托克斯流（低雷诺数） — 故障排除指南
故障排除
低Re数分析中容易出现什么故障吗？
与高Re数不同，会出现不同类型的问题。主要是精度和物理模型的问题。
常见问题和对策
1. 压力场不稳定（FEM情况）
症状: 压力呈棋盘式振荡，或压力接近零的异常值。
原因: 速度-压力的要素组合不满足inf-sup条件（LBB条件）。
对策:
- 使用Taylor-Hood要素（P2/P1）
- 应用稳定化方法（PSPG, GLS）
- COMSOL中，Creeping Flow物理会自动选择稳定的公式
2. 粒子的阻力不符合斯托克斯定律
用CFD计算球的阻力，结果与 $6\pi\mu RU$ 相差10%以上…
| 原因 | 对策 |
|---|---|
| 计算域太小 | 确保球径50倍以上的外部区域 |
| 网格太粗 | 球表面至少30个单元（周向） |
| 壁面影响 | 无穷远条件（远场边界处均匀流） |
| Re > 0.1 | 考虑Oseen修正: $C_D = (24/Re)(1 + 3Re/16)$ |
斯托克斯解假定无穷域，所以有限的计算域中壁面效应不可避免。将外部边界放置得足够远很重要。
3. 微流道流量不符合理论值
症状: 哈根-泊肃叶的理论流量偏差很大。
确认事项:
- 流道截面不是圆管时，需要根据截面形状应用修正系数
- 入口助流区（$L_e \approx 0.06 Re D$，斯托克斯流中几乎为零）
- 壁面滑移条件: 在微观尺度下，克努森数 $Kn = \lambda/L$ 有限，no-slip可能不成立
4. 非牛顿流体收敛困难
症状: 幂律流体中粘度在局部非常大（或很小），不收敛。
对策:
- 设置粘度的上下限（$\mu_{\min} < \mu_{\text{eff}} < \mu_{\max}$）
- 使用Carreau模型（由零剪切粘度和无穷剪切粘度规定上下界）
- 降低松弛因子
- 使用牛顿流体（平均粘度）的解作为初始值
5. 多粒子问题中相互作用被低估
粒子很多的情况下，逐个计算不行吗？
孤立粒子的假设仅在粒子间距足够大时有效。体积分率超过0.1时，粒子间流体动力学相互作用不能忽视。
对策:
即使是斯托克斯流，实际工作中也有各种问题呢。
即使方程是线性的，分析也不一定简单。需要理解尺度问题、壁面效应、非牛顿性等低Re特有的课题。
Coffee Break 杂谈
斯托克斯流计算"不符合"—— 怀疑Re数
斯托克斯流（Re≪1）分析中最常见的失败是"低速所以可以用斯托克斯近似"的误解。Re数不仅取决于流速，还取决于代表长度。微通道内的1μm粒子Re≈0.001，完全处于斯托克斯区域，但相同流速下1mm气泡的Re≈1接近，惯性项开始生效。当分析与实验值不符时，首先要重新检查问题的代表长度——这是斯托克斯流故障排除的铁则。
相关话题
本文评价
感谢您的回答！
有参考
价值
价值
需要
更详细
更详细
报告
错误
错误