涡度方程式 — 故障排除指南
更丰富的内容请查看 vorticity-equation.html。
涡度分析常见问题
涡度分析结果不正确时，应检查哪些内容？
整理涡度分析特有的故障及对策。
1. 涡在下游消失（数值扩散）
我解了圆柱的卡门涡，但涡在下游很快就消失了。
这是最常见的问题。原因是数值扩散（numerical diffusion），需要检查以下内容。
| 检查项目 | 常见原因 | 对策 |
|---|---|---|
| 离散化方案 | 使用一阶迎风差分 | 改为二阶或更高精度（QUICK, MUSCL, Central等） |
| 网格方向 | 流向与网格线斜向 | 将网格与涡的移流方向对齐 |
| 网格密度 | 涡核相对于单元较粗 | 涡核直径至少10个单元 |
| 网格类型 | 仅使用四面体网格 | 将涡通过区域改为六面体网格 |
对于Fluent，应将Spatial Discretization > Momentum设为Second Order Upwind或更高。理想情况下，使用Bounded Central Differencing（用于LES）可将数值扩散最小化。
2. 涡脱离未开始
在圆柱周围的非稳态分析中，涡始终没有出现。
可能的原因和对策如下。
- 完全对称网格：上下对称的网格不会产生数值扰动。应在网格中引入微小的非对称性，或在初始条件中添加微小扰动
- 使用稳态分析：应切换到非迭代时间推进，或明确选择瞬态分析
- 时间步长过大：应减小 $\Delta t$ 使 CFL < 1
- 雷诺数过低：Re < 47 时圆柱后流是稳定的，不会产生涡脱离
3. 壁面涡度振荡
壁面上的涡度分布呈锯齿状，这是正常的吗？
壁面涡度为 $\omega_{wall} = -\frac{\partial u_{tangential}}{\partial n}$，所以当壁面附近的速度梯度不能精确解析时就会出现振荡。对策为
- 确认壁面第一单元的 $y^+$：使用壁面函数时 $30 < y^+ < 300$，壁面解析时 $y^+ < 1$
- 网格膨胀率：壁面向外扩展方向的膨胀率应控制在1.2以下
- Fluent特定注意：使用基于节点的梯度（Node-Based Gradient）可改进壁面附近的梯度评估
4. 拟涡度不守恒
在非粘性流的涡度分析中，拟涡度在增长。
在非粘性（$\nu=0$）的二维流中，拟涡度应该守恒。如果增长，则说明数值方案的耗散和色散特性有问题。
- 守恒型方案：检查是否使用守恒型（divergence form）对流项离散化
- 能量守恒方案：使用严格保守动能的离散化（Morinishi et al., 1998）
- 伪谱法：对于周期边界条件，基于FFT的谱法对拟涡度守恒最有优势
OpenFOAM中的涡度调试
OpenFOAM中出现涡度相关错误时的处理方法请教一下。
总结代表性的故障及对策。
| 错误/症状 | 原因 | 对策 |
|---|---|---|
Maximum Courant Number: xxx 过大 | 违反CFL条件 | 设置 maxCo 0.5，启用 adjustTimeStep yes |
| 涡度后处理失败 | functionObjects配置错误 | 在 system/controlDict 中添加 functions { vorticity { type vorticity; libs (fieldFunctionObjects); } } |
| 涡度场中的尖刺 | 网格质量不良 | 用 checkMesh -allGeometry -allTopology 检查质量，修复 skewness > 4 的单元 |
Fluent中涡度计算结果有时会出现NaN。
NaN出现的主要原因是发散。以下是涡度分析中容易发散的情况及对策。
- 高Re数初期过渡：从稳态解开始切换到瞬态。或者只在前几十步使用一阶精度，稳定后再切换到二阶精度
- 低松弛因子：将压力设为0.2~0.3，动量设为0.5~0.7
- Fluent的TUI命令：
/solve/set/equations/flow no可冻结流场，先收敛湍流变量再求解流场。这也是一个有效技巧
涡度分析是速度场的导数，所以速度场本身的精度很关键呢。
完全正确。涡度是速度的一阶导数，变形速率张量也是一阶导数。因此速度场的二阶精度只能得到涡度的一阶精度。如果想提高涡度精度，需要用三阶或更高精度求解速度场，或者确保充分的网格分辨率。
涡度的"发散"是真实故障还是数值误差
在直接求解涡度方程的（ψ-ω法）CFD代码中，经常出现涡度急剧增长的问题。在慌张之前，应该确认这是真实流动（例如：涡旋增强的分离）还是数值误差积累。特别是壁面涡度边界条件的设置错误比较明显，用ω=−∂²ψ/∂n²近似时精度低会导致涡度异常蓄积在壁面附近，整体看起来像发散。正确的处理方法是细致地可视化壁面附近的涡度分布。
相关主题
更详细
错误