Fluent模拟计算不收敛“Floating Point Exception”常见原因及其解决办法

   在使用Fluent进行计算时，“Floating Point Exception”是一种常见的错误，它意味着浮点溢出，即计算不收敛，通常与数值计算的稳定性或设置相关。这种错误可能会导致计算中断或无法收敛。出现该错误的常见原因及其解决方法如下：

1、网格质量问题

优化网格形状与类型，在边界层或关键区域，可以使用六面体网格来提高精度，或者在流动剧烈变化的区域采用局部加密网格；控制网格单元扭曲度（skewness）在0.8以下；在网格生成后，进行网格光顺可以有效提高网格的质量，光顺操作可以优化网格单元的形状，减少扭曲度并提高正交性。

2、时间步长过大（瞬态）

首先我们要清楚几个参数的含义及关系，比如：时间步长×时间步数=总模拟时间，如果你选择了时间步长为0.5s，时间步数为200步，那么总模拟时间将是100s，这意味着你的模拟将在物理时间上计算100秒的仿真过程。小的时间步长能够提高模拟的精度，但会增加计算量。反之，较大的时间步长可以加快计算速度，但可能导致不收敛或误差增大。

而最大迭代数是指，在每个时间步内，求解器允许的最大迭代次数，该参数用于控制每个时间步的计算量。最大迭代数的设置应确保每个时间步内的解能够充分收敛，避免过度迭代浪费计算资源。

最后我们再回到如何通过调整时间步长来解决浮点溢出的情况，时间步长∆t 我们大致可以通过流速v与最小网格尺寸∆x 估算出来：

时间步长比 ∆t→∆x/ v 所得的time小一到两个数量级或更小。

3、初始条件设置不合理

如果初始条件设置不当，例如过于偏离实际情况，求解器可能无法在早期迭代中找到合理的解，导致数值不稳定。例如，先进行稳态模拟，再使用结果作为瞬态模拟的初始条件。

   

4、边界条件设置错误

对于入口和出口的设置，应确保速度、压力或温度的数值合理，流动湍能、温度、可压缩流体压力保证不为零；或者检查单位是否合理（欧拉多相流模型数值设置稳定性较差，应格外注意）；此外，速度设置为0可能是导致计算不收敛的关键。

5、求解器选择不合理

不同的求解器用于解决不同类型的流动问题和方程，选择合适的求解器可以显著影响模拟的精度、收敛速度和稳定性。下面介绍几种常用的求解器类型，并提供选择指南。（如若不收敛，建议使用SIMPLE算法）

  

(1)SIMPLE

SIMPLE算法是一种半隐式压力-速度耦合算法，用于解决不可压缩流体流动的压力和速度耦合问题，主要用于不可压缩流动和低速流动问题，适合工程应用中的大多数稳态和非稳态流动问题，如管道流动、风洞实验模拟等。

(2)SIMPLEC

SIMPLEC算法是对SIMPLE算法的改进，目的是加快收敛速度。它与SIMPLE的基本求解过程相似，但在压力修正和速度修正的步骤中进行了简化，减少了不必要的迭代次数，特别是在高雷诺数流动中，SIMPLEC能显著提高收敛效率。

SIMPLEC主要用于稳态流动问题，在非定常问题上相比SIMPLE或PISO，优势不明显，特别适用于收敛速度较慢的高雷诺数或湍流流动场，能够显著提升计算效率。

(3)PISO

PISO算法是一种压力-速度耦合方法，最早用于求解不可压缩流动问题。与SIMPLE相比，PISO特别适合用于非稳态流动的求解。PISO算法通过多次修正压力场来避免每个时间步长的多次迭代，从而在非定常流动计算中效率更高。

对于需要高精度压力-速度耦合的瞬态流动问题，PISO是最佳选择。

(4)Coupled

Coupled 求解器是一种全隐式的压力-速度耦合算法，不同于SIMPLE和SIMPLEC等半隐式方法，它在每个时间步长内同时求解动量方程和连续性方程。这种全耦合的方式可以显著提高求解速度，尤其在高马赫数、高雷诺数的复杂流动场中表现优越。

当模拟的流动具有高马赫数（如超音速流动）或涉及显著的压缩效应时，Coupled求解器是最佳选择。

6、处理器选择不合理

合理选择并统一使用单个求解器，并进行适当的迭代和数值精度控制，可以有效降低浮点溢出问题。

在并联计算中，优化处理器之间的通信设置、合理的网格划分、以及对收敛标准的严格控制同样能够有效避免数值发散和浮点溢出。

7、亚松弛因子设置不当

在Fluent中，对于大多数求解，亚松驰因子中的所以变量会默认设置为最佳参数，但对于一些非线性的问题，就会容易导致不收敛。所以我们可以通过减小亚松驰因子来解决浮点溢出的问题。

将压力、密度、动量、湍流动能、湍流耗散率的欠松弛因子改为0.3、0.5、0.5、0.5、0.5，对于SIMPLE求解方法可以不必改变压力参数；如若仍不收敛便将其减少到0.3、0.3、0.3、0.3、0.3。

   

8、物理模型选择问题

模型不适合当前模拟工况（如湍流模型或相互作用模型的选择不当）；物理模型之间的耦合问题（如多相流、反应流模型的交互导致的不收敛）。

在模拟湍流时，尝试从标准“k-epsilon”模型转换为“SST”模型，可以提高收敛性，尤其是在处理边界层时。