• Fluent计算出现“Floating point exception”时的解决办法


            “Floating point exception”代表的含义是浮点溢出,一般在运行某一算法的时候无法继续计算时出现。通常来说,浮点溢出就代表着计算不收敛。下面提供几种解决这个问题的方法:

    1.检查边界条件的设置:

            检查以下各项,保证其数值不为0:流动湍能(turbulent kinetic energy)、温度、可压缩流体仿真中的压力。

    2.调整时间步长:

            当时间步长过大时也会造成计算结果不收敛,尤其是进行第一次迭代时,流体速度通常会比较小,CFL准则会导致时间步长设置过大。

    3.改变求解器:尽可能使用耦合显式求解器,也可以选择GMAG或者平滑求解器。平滑求解器在矩阵条件不好的时候有效。

    4.处理器:

            串联或并联计算时,使用单个求解器可能能够解决浮点溢出的问题。

    5.提高网格质量:

            调整网格参数,重新划分网格,保证网格歪斜度(skewness)在0.8以下,或尝试使用conformal网格。

    6.使用亚松弛(under-relaxation)方法:

            Fluent默认的更新方法是利用新计算的数值代替旧值。在进行仿真计算时,一般会根据初始条件估算出下一步的结果,但大多数情况下估算结果会与真实结果差距较大,因此通过迭代的方法进行计算,收敛得到真实结果。

            亚松弛方法的原理是当计算出新指后,不以新值作为下一步迭代的初始值,而是按权重组合旧值和新值作为下一步的初始值。

            如果亚松弛因子为1,则代表新值的权重为1,旧值为0,即完全使用新值作为下一步的初始值;反之,若松弛因子为0,则代表新值的权重为0,旧值为1,即以上一步的旧值作为下一步的初始值。

            一般来说,对于速度的亚松弛因子取0.7,对于压力及其他标量为0.3.如果计算结果不收敛,可以尝试使用更低的亚松弛因子,例如速度亚松弛因子取0.3,压力及其他标量取0.1。

    7.使用upwind convection schemes:

            upwind convection schemes引入了数植扩散,可以确保在特定的情况下计算收敛。对于稳态仿真,使用bounded Gauss upwind schemes更加合适,因为这个能够确保计算满足边界条件。

            另外,还要避免使用高阶项,因为高阶项更容易导致计算结果不收敛。

  • 相关阅读:
    相同切入点的抽取
    【BOOST C++ 5 】通信(02 可扩展性和多线程 )
    ElementUI中使用ECharts
    C++--第一个代码hello world
    【C++初阶7-string】真方便,真舒服
    Design Compiler工具学习笔记(6)
    MHA实验
    Linux从入门到入土③(GCC、静态与动态链接库、Git、Make、Makefile)
    Shiro授权&&Shiro的注解式开发
    [SWPUCTF 2022 新生赛]善哉善哉(隐写,新佛曰,MD5)
  • 原文地址:https://blog.csdn.net/qq_50981222/article/details/126231585