圆柱周围的流动 — 故障排除指南
本文已迁移至统合版
更充实的内容请见 cylinder-flow.html。
更充实的内容请见 cylinder-flow.html。
圆柱周围的流动 — 故障排除指南
常见问题和对策
请教圆柱周围流动计算中经常遇到的要点。
让我们整理一下实务中频繁出现的问题。
1. 没有出现涡（收敛到定常解）
原因和对策:
- 使用一阶精度迎风格式 → 数值扩散消除涡。至少改为二阶精度
- 计算域过窄 → 堵塞效应改变实际Re。侧面保持 $10D$ 以上距离
- 初始条件对称 → 对称性无法破裂。在初始条件中加入微小的非对称扰动
- 使用定常求解器 → 切换到非定常求解器（PISO/PIMPLE）
- 时间步过大 → 确认 CFL < 1 。每个涡放出周期至少200个时间步
在初始条件中加入扰动，这在物理上是正确的吗？
实际流动中总是存在微小扰动。数值计算中没有这些扰动，所以即使不稳定也可能不会发生分岔。在$y$方向速度上加入 $O(10^{-3})$ 量级的随机噪声就足够了。
2. 斯特劳哈尔数不匹配
确认要点:
- 时间平均统计是否充分 → 放弃初始暂态（至少20个周期）后再统计
- 网格分辨率 → 尾流域网格过粗会导致涡快速耗散，周频率改变
- 堵塞效应 → 计算域过窄时St会升高。存在Behr等人的修正公式
- FFT窗函数 → 使用足够长的时间数据进行FFT。零填充以保证频率分辨率
3. $C_D$ 与实验值不匹配
$C_D$ 算出来是 $1.5$ 这样的值。
Re=100时如果 $C_D > 1.4$，请确认以下几点。
- 是否用2D计算与3D实验值比较: 2D计算的 $C_D$ 比3D实验值高 $5\text{--}10\%$。这在物理上是正确的（3D效应降低阻力）
- 壁面网格分辨率: 确认壁面压力分布是否得到正确分辨
- 出口边界影响: 出口过近会改变背压，从而影响 $C_D$
4. 计算发散
按症状分别对策:
| 症状 | 原因 | 对策 |
|---|---|---|
| 库朗数爆炸 | 时间步过大 | 将 $\Delta t$ 减半。使用adjustableRunTime |
| 压力振荡 | 棋盘现象 | 确认Rhie-Chow插值。非结构网格则确认单元心配置 |
| 速度在壁面附近发散 | $y^+$ 超出模型适用范围 | 壁函数情况：$30 < y^+ < 300$；壁面分辨情况：$y^+ < 1$ |
| LES 中出现非物理涡 | SGS模型不当 | Smagorinsky ($C_s = 0.1$)在壁面耗散过大。改用WALE |
5. 3D计算中跨向周期性不出现
3D计算中跨向只出现均匀的涡。
可能是跨向计算长度不足，或跨向网格过粗。要捕捉Mode A需要 $L_z > 4D$，捕捉Mode B需要 $\Delta z < 0.5D$。同时，跨向端面必须设置周期边界条件。若设为壁面条件，端部效应会占主导地位。
原来如此。先用2D在低Re验证，再进行3D，是安全的方法呢。
正确。首先用2D·低Re验证 $C_D$, St → 用3D低Re确认Mode A/B → 逐步升高到目标Re，这是最稳妥的方法。
Coffee Break 闲谈小角
"涡呈现左右对称"——暂态解析的初始化问题
圆柱周围非定常计算中"卡门涡无法发展，涡呈现左右对称地流出"是初学者容易踩到的坑。原因是"初始场完全对称，没有数值的对称性破裂扰动"。Navier-Stokes方程在左右对称的边界条件和初始条件下会继续输出对称解。要引起涡放出需要微小的非对称扰动。解决方法是"在圆柱略微非对称的位置放置网格节点"，或"仅在开始几个时间步加入微小的横向速度扰动"。"物理上应该发生的现象在CFD中没有发生"通常源于初始条件的对称性，这是很好的例子。
相关主题
本文评价
谢谢您的回答！
有参考
价值
价值
更详细
的内容
的内容
报告
错误
错误