文章来自——【小白的CFD之旅】26 何为收敛 - 流沙[胡坤] - 博客园
在CFD计算中,每一个网格上都会存储众多物理量,因此每一个网格上的任一个物理量在计算迭代过程中都会存在一个残差,这意味着在一次迭代过程中,同一物理量在不同的计算网格上有不同的计算残差,而实际上我们在进行CFD计算时,每一个迭代步只对应着一个残差值。
CFD中残差分为几种:
最大残差:在一次迭代中,取所有网格中的残差值的最大值作为本次计算的残差。
平均残差:在一次迭代中,计算所有网格中的算术平均值作为本次迭代的残差
均方根残差:在一次计算,计算所有网格中残差值的均方根作为本次迭代计算的残差
在CFD计算中,常常采用均方根残差(RMS)作为残差值。
所谓迭代收敛,简单来说,就是在迭代计算过程中,物理量趋于某一值的情况。CFD计算中判断收敛通常有三种方法。
收敛判断规则之一:残差达到某一设定标准时可以认为迭代计算达到收敛。
这条规则最简单,在实际工程应用中也最常用。通常设定某一标准,当迭代计算过程中残差值低于此标准时则认为计算收敛。这也是几乎所有CFD软件用于判断收敛的基本方法。
然而此规则在实际工程应用中常常无法使用,有些复杂的问题,无论你怎么计算,其残差也不会下降,甚至有时候残差会出现周期性震荡。
残差稳定在某一位置不下降的原因有很多,常见的原因包括:
- 计算区域中存在低质量的网格。低质量的网格会造成计算残差增大及残差震荡
- 边界条件设置有误。错误的边界条件或边界类型搭配都会导致计算残差震荡。
- 利用稳态求解器计算瞬态问题也会造成残差的震荡。
在残差无法达到设定标准的时候该如何判断收敛呢?需要注意的是,此时的收敛并非数学意义上的收敛了,而只是意味着我们可以停止计算。
收敛判断规则之二:进出口物理量通量达到平衡
最常用的是判断进出口质量是否相等,这实际上是判断连续性条件是否达到满足的。实际上还有很多,如计算域中包含化学反应时,判断进出口组分是否守恒;如计算域中包含多相流时,判断进出口各相质量是否守恒等。此规则是一种非常弱的规则,实际上只是收敛的一个必要条件而已。但是在第一条判断规则无法达到时,也常常采用此规则来判断。
收敛判断规则三:计算域中的物理量随迭代进行不再发生变化
这条规则在实际应用中也很常用,甚至比第二条规则更常用。在实际工程中,经常监测某些敏感位置的物理量,当随着迭代进行,监测的量不再发生变化时,基本可以认为计算达到收敛。
注意后两种方法只是在残差无法达到标准时才采用的判断方法,并不意味着计算就收敛了。