CFD分析的发散

从数值的角度看，残差单调增加，解变成物理上没有意义的极端值（例如，局部压力达到10^30 Pa）的现象。从物理角度看，这意味着设定的条件超出了支配方程（纳维-斯托克斯方程）能够稳定保持解的范围。例如，流速太快导致流动分离，计算网格无法捕捉到这种急剧变化。

是的，有一个重要的概念叫CFL（Courant-Friedrichs-Lewy）条件。这是对显式方法的时间步长

理论上隐式方法被称为无条件稳定，但在实际的非线性问题中并不那么简单。如果时间步长过大，非线性求解器（例如牛顿法）的收敛性会恶化，最终导致发散。Ansys Fluent的默认时间步长设置需要考虑所模拟现象的时间尺度。例如，在求解圆柱体周围卡门涡的情况下，应根据Strouhal数从主流速度和圆柱直径估算特征时间，将初始

在初期计算或流动复杂的区域，应从"一阶迎风差分"开始，而不是高阶精度格式。这是因为数值扩散较大，易于稳定解。在OpenFOAM中，`divSchemes`设置可以选择`Gauss linearUpwind`等。待解稳定后，可切换到二阶精度的`Gauss linear`或`limitedLinear`来提高精度。如果从一开始就使用`Gauss QUICK`这样的高阶格式，发散风险会大大增加。

有很大的区别。SIMPLE方法稳定性高但收敛慢，是定常分析的标准方法。PISO方法用于非定常分析，如果不增加每个时间步的修正次数，容易发散。Ansys Fluent的"耦合"求解器同时求解压力和速度，在复杂流动或浮力流中强大，但内存消耗大，初始条件不好时可能立即发散。实际应用中的一个策略是，先用SIMPLE处理复杂的内部流动得到初始场，再切换到耦合求解器。

在发散趋势或发散直后，应从默认值降低。例如在Fluent中，压力的松弛因子默认为0.3，在问题时可降至0.1或0.05。动量松弛因子从0.7降至0.3。这做的是减少更新量，抑制解的急剧变动。但收敛速度会极大地放慢，所以待稳定后要逐步恢复原值。反之，如果只是收敛太慢而不是发散，可提高松弛因子（但上限为0.9）来加速。

1) 立即停止计算。2) 从最后保存的数据文件（例如Fluent的`.dat`文件）恢复。3) 检查"残差历史"和"区域内最小/最大值（例如最小压力）"。定位到压力为-10^5 Pa等物理上不可能的值的单元。4) 检查该单元周围的网格质量（畸变率、长宽比）。这是最常见的初始操作。

一般最佳实践是：歪斜度（Skewness）小于0.9（理想值为0.5以下），长宽比（Aspect Ratio）100以下（边界层可超过1000但需注意），直交性（Orthogonal Quality）大于0.1（越接近1越好）。Ansys Meshing或Pointwise可以生成这些指标的报告。特别需要严格检查流速或压力梯度大的区域的网格质量。

有几种方法。第一，如前所述降低松弛因子。第二，启用梯度限制（Gradient Limiting）。Fluent中可在"Solution Limits"中设置最小绝对压力（例如-1e5 Pa）或最大温度，使解超出范围时被截断。第三，临时切换到更稳健的离散化格式。这些不是根本解决方案，但有助于继续计算并定位问题区域。

区别很大。Fluent在GUI上将残差显示为红色"Diverged"，易于监控最小/最大值。求解器还会自动发出"Turbulent Viscosity Limited"等警告。相比之下，OpenFOAM在命令行输出中显示`time step continuity errors`急剧增加（例如1e+10）或因`solution singularity`而停止。OpenFOAM有一个`controlDict`中的`maxCo`（最大Courant数）设置，超过时自动削减时间步，对防止发散很有效。

与其说强弱，不如说方法不同。COMSOL默认使用"分离式求解器"，按顺序求解各物理场。稳定性很高，但对有强耦合的问题（如伴有化学反应的高速流动）可能收敛困难。此时需切换到"完全耦合求解器"。COMSOL的优势在于有求解器的"稳定化"设置，通过自动阻尼来防止发散。但网格依赖性与其他软件相同。

有的，在"Implicit Unsteady"求解器的"Robustness Settings"中有特色功能。包括"Pressure Variation Damping"和"Momentum Damping"，当发散趋势出现时求解器自动增加阻尼系数。还有"Accuracy/Stability Optimizer"功能，对网格质量差的区域自动添加数值扩散。早期的基准测试显示，即使刻意给予劣质网格，Star-CCM+的计算比其他软件更不容易失败。

这是非常普遍的问题。原因是流入边界或壁面附近设定的乱流强度I或乱流粘度比μ_t/μ过大，或初始场的k和ε不平衡。具体来说，如果初值不满足

界面急剧变形会产生体积分率的陡梯度，对流项离散化容易破坏。对策是：1) 在界面附近充分细化网格（至少用5-10个单元来解析界面）。2) 界面捕捉格式从高阶但不稳定的`Compressive`系改为稳定性优先的`Geometric Reconstruction`或`High Resolution Interface Capturing`（Fluent）开始。3) 严格控制时间步，使界面的Courant数保持在0.25以下。在OpenFOAM的`interFoam`中，增加`nAlphaCorr`（迭代次数）对`alpha.water`场也有效。

原因有重叠之处，但对应方式不同。停滞是求解器被困在解的局部"凹陷"无法逃脱的状态。对策包括：1) 改变初始猜测值（从另一个稳定案例读取初始场）。2) 更换求解器（从SIMPLE改为耦合型）。3) 临时降低离散化格式精度以"平滑化"。相比之下，发散是解失控暴走的状态，需要更强的稳定化措施，如降低松弛因子、限制解的范围或修改网格。可以这样理解：停滞是"卡住"，发散是"失控"。