流固耦合,换了质量更好的网格流体突然就不收敛了?
-
各位大佬好,我是个openFoam的新手,最近在用preCICE、openFoam、calculix做三维柔性板的流固耦合案例,之前用workbench随便画了个四面体网格(结构用的abaqus画的六面体网格),虽然在耦合界面上没有收敛,但是单相的流体或结构体都收敛了。
我就寻思是网格的问题(流固交界面插值的问题),就重新用ICEM画了个流场的六面体网格,但是这下好了,在相同配置下连带着流体也不收敛
最后的附件是我的计算配置文件和log文件,图片是我的网格文件,我想问下是重新画的网格不太好吗?后续fvSchemes和fvSolution应该按照什么思路修改?
谢谢大佬们!!这是计算模型
这是之前随便画的网格
这是后来用ICEM画的
计算设置文件
controlDict fvSchemes fvSolutionlog文件
log -
你好
PIMPLE: iteration 40 **time step continuity errors : sum local = 20238.39271, global = -1326.825454, cumulative = -1286.292262** **time step continuity errors : sum local = 1965.022984, global = -493.7266286, cumulative = -1780.018891** PIMPLE: iteration 22 GAMG: Solving for p, Initial residual = 0.1839153372, Final residual = 1.265686901e-06, No Iterations 20 GAMG: Solving for p, Initial residual = 0.1398796637, Final residual = 1.397575026e-06, No Iterations 17 GAMG: Solving for p, Initial residual = 0.1463017811, Final residual = 1.437462198e-06, No Iterations 16 GAMG: Solving for p, Initial residual = 0.121320625, Final residual = 6.97291136e-08, No Iterations 20 PIMPLE: iteration 23 .GAMG: Solving for p, Initial residual = 0.1600126604, Final residual = 1.256444652e-06, No Iterations 23 GAMG: Solving for p, Initial residual = 0.09022545056, Final residual = 7.317342264e-07, No Iterations 20 GAMG: Solving for p, Initial residual = 0.08084307981, Final residual = 7.297894683e-07, No Iterations 18 GAMG: Solving for p, Initial residual = 0.06296048891, Final residual = 7.500379861e-08, No Iterations 22
你的问题在于:
1、瞬态项上,你的瞬态项很难收敛,残差爆掉了,这个似乎与你的瞬态项离散方式有关,可以试一下二阶向后差分。ddtSchemes { default CrankNicolson 0.9; }
2、你的泊松方程收敛性很差,外循环很难收敛,残差很高,你定义的外循环收敛条件根本达不到,泊松方程在非正交修正后,迭代次数反而变多,证明你的网格质量可能较差。
建议:1.更换瞬态项离散方式。2.增加非正交修正次数,提高你的网格质量。
注意:外循环定义50次合理,一般难收敛问题,最多10次就可以收敛了。流固耦合交界面前几个时间步不收敛很正常,但是你这个是流体不收敛,这就不对了,建议你先单独跑一下流体看看,你这个流体设置不太好。
祝好~ -
@Jimmy-Yao 楼主最开始的是四面体网格吗?