①、一般首先是改变初值,尝试不同的初始化,事实上好像初始化很关键,对于收敛。 ②、FLUENT的收敛最基础的是网格的质量,计算的时候看怎样选择CFL数,这个靠经验 ③、首先查找网格问题,如果问题复杂比如多相流问题,与模型、边界、初始条件都有关系。 ④、有时初始条件和边界条件(如很多复杂室内流场pressure outlet 就比outflow边界条件的收敛性要好,因为很多时候室内回风口要符合outflow边界条件的fully developed turbulence 还是很难的)严重影响收敛性,曾经作过一个计算反反复复,通过修改网格,重新定义初始条件,包括具体的选择的模型, 还有老师经常用的方法就是看看哪个因素不收敛,然后寻找和它有关的条件,改变相应参数。就收敛了 ⑤、A.检查是否哪里设定有误:比方用mm的unit建构的mesh,忘了scale;比方给定的边界条件不合理。B从算至发散前几步,看presure分布,看不出来的话,再算几步, 看看问题大概出在那个区域。 C网格,配合第二点作修正,就重建个更漂亮的,或是更粗略的来处理。D再找不出来的话,换个solver。 ⑥、解决的办法是设几个监测点,比如出流或参数变化较大的地方,若这些地方的参数变化很小,就可以认为是收敛了,尽管此时残值曲线还没有降下来。这时再检查下流量是否守恒,若这时流量也守恒,则确实是收敛了------对很多复杂的室内流场计算问题经常出现这个情况 ⑦、调节松弛因子也能影响收敛,不过代价是收敛速度。
看了流量是否平衡吗?在report->flux里面操作,mass flow rate,把所有进出口都选上,compute一下,看看nut flux是什么水平,如果它的值小于总进口流量的1%,并且其他检测量在继续迭代之后不会发生波动,也可以认为你的解是收敛的。 造成连续方程高残差不收敛的原因主要有以下几点: 1.网格质量,主要可能是相邻单元的尺寸大小相差较大,它们的尺寸之比最好控制在1.2以内,不能超过1.4. 2.离散格式及压力速度耦合方法,如果是结构网格,建议使用高阶格式,如2阶迎风格式等,如果是非结构网格,除pressure保持standard格式不变外,其他格式改用高阶格式;压力速度耦合关系,如果使用SIMPLE,SIMPLEC,PISO等segerated solver对联系方程收敛没有提高的话,可以尝试使用coupled solver。另外,对于梯度的计算,不论使用结构或非结构网格,都可以改用node-based来提高计算精度。 |
|