SIMPLE法 — 故障排除

SIMPLE法 — 故障排除

常见问题及解决方案

1. 残差振荡而无法收敛

原因和对策:

• 松弛系数过大 → 将压力降至0.2，速度降至0.5
• 网格非正交性高 → 将 nNonOrthogonalCorrectors 增加到2～3（OpenFOAM）
• 离散化格式过高 → 暂时降低为一阶迎风，使其稳定

2. 残差停滞而不下降

原因和对策:

• 网格品质不良 → 检查单元歪度、长宽比
• 边界条件不匹配 → 确认入口/出口质量流量平衡
• 湍流模型初值不当 → 根据流速重新估算 $k$ 和 $\epsilon$ 的初值

3. 残差发散

原因和对策:

• CFL数过大（使用伪时间推进法时）→ 减小伪时间步长
• 边界条件出现反流 → 设置出口为压力出口，确认反流条件
• 网格存在负体积单元 → 用网格生成器重新生成

4. 连续方程残差无法下降

• 增加AMG前处理循环数（V循环 → W循环）
• 增加压力线性求解器迭代次数（50→200）
• 检查网格粗化比（coarsening ratio）

求解器特定的典型错误

• Divergence detected in AMG solver: pressure correction → 降低压力松弛，确认网格品质
• Reversed flow in N faces → 将出口边界远离流动区域

OpenFOAM:

• GAMG: Solving for p, Initial residual = 1, Final residual = 1 → 压力方程无法求解。确认边界条件类型（fixedValue vs zeroGradient）
• Floating point exception → 负体积单元、零除法。运行 checkMesh

STAR-CCM+:

• Diverging solution - Loss of fluid cells → 检查单元品质，降低时间步长因子

调试步骤

1. 检查网格品质: checkMesh（OpenFOAM）、网格品质报告（Fluent）

2. 简化边界条件: 首先用全壁面进行测试

3. 一阶精度稳定计算: 收敛后切换到二阶精度

4. 保守设置松弛系数: 从压力0.2、速度0.5开始

5. 观察残差变量推移: 确定哪个变量存在问题

6. 可视化物理量场: 确认没有非物理区域