流固耦合,换了质量更好的网格流体突然就不收敛了?
-
各位大佬好,我是个openFoam的新手,最近在用preCICE、openFoam、calculix做三维柔性板的流固耦合案例,之前用workbench随便画了个四面体网格(结构用的abaqus画的六面体网格),虽然在耦合界面上没有收敛,但是单相的流体或结构体都收敛了。
我就寻思是网格的问题(流固交界面插值的问题),就重新用ICEM画了个流场的六面体网格,但是这下好了,在相同配置下连带着流体也不收敛
最后的附件是我的计算配置文件和log文件,图片是我的网格文件,我想问下是重新画的网格不太好吗?后续fvSchemes和fvSolution应该按照什么思路修改?
谢谢大佬们!!计算设置文件
controlDict fvSchemes fvSolutionlog文件
log -
@李东岳 在 流固耦合,换了质量更好的网格流体突然就不收敛了? 中说:
所以99%是设置的问题
好的,感谢李老师,我再尝试着调一调设置文件吧
感觉还是fvSchemes和fvSolution的问题,和新的网格不太匹配 -
-
很好奇,现有的OpenFOAM+Calculix耦合框架能用于计算风机类的旋转情况吗?
-
@李东岳 在 流固耦合,换了质量更好的网格流体突然就不收敛了? 中说:
单独openfoam就能计算风机
李老师,我是想算风机的流固耦合。最近利用STAR-CCM+加Abaqus算风机流固耦合,想着能不能用OpenFOAM+Calculix平替了
-
@李东岳 在 流固耦合,换了质量更好的网格流体突然就不收敛了? 中说:
openfoam可以算风机,流固耦合也可以算。所以单独openfoam就可以用。如果你是工程用,那可以用现成的,可能calculix会更适合上手
初学阶段,商软教程多一点,打算先用商软走一遍,学有余力再进OpenFOAM,感觉是个大工程。谢谢李老师回复!
-
你好
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 楼主最开始的是四面体网格吗?
-
@ShuanmingZhao 是的,我基本所有网格都用过了,发现只有纯六面体网格才行,可能就是开源软件的问题
-
lz您好,小弟最近也在尝试类似的case。想请问您这个三维柔性板的计算是基于动网格还是浸没边界法呢?
6/20