流固耦合，换了质量更好的网格流体突然就不收敛了？

[Jimmy Yao]

@李东岳 这个是照着油管上的一个算例做的设置，那个老师说多迭代几步能提高耦合界面的收敛性

[李东岳]

那个老师说多迭代几步能提高耦合界面的收敛性

crazy...

最近在用preCICE、openFoam、calculix做三维柔性板的流固耦合案例

目前我这面只搞过openfoam，没有更建设性的意见。不过这么简单的工况，并且这几个软件都比较成熟，所以99%是设置的问题

[Jimmy Yao]

@李东岳 在 流固耦合，换了质量更好的网格流体突然就不收敛了？ 中说：

所以99%是设置的问题

好的，感谢李老师，我再尝试着调一调设置文件吧
感觉还是fvSchemes和fvSolution的问题，和新的网格不太匹配

[coolhhh]

猜测跟网格问题比较大，ICEM这种o-block网格质量其实很差，之前画过一个类似的网格，用fluent的PISO算法算就是发散，这网格连fluent都发散。后面改用fluentmeshing或者snappyHexMesh画的网格就质量很好，并且能够方便设置不同加密区的网格也是更加适用研究问题

[Jimmy Yao]

@coolhhh 在 流固耦合，换了质量更好的网格流体突然就不收敛了？ 中说：

改

是的，大佬
昨天重新画了下网格，结果收敛了
猜测是O切的问题，最后直接用线切分块，流场也顺利收敛了

[HITSC30]

很好奇，现有的OpenFOAM+Calculix耦合框架能用于计算风机类的旋转情况吗？

[李东岳]

单独openfoam就能计算风机

[HITSC30]

@李东岳 在 流固耦合，换了质量更好的网格流体突然就不收敛了？ 中说：

单独openfoam就能计算风机

李老师，我是想算风机的流固耦合。最近利用STAR-CCM+加Abaqus算风机流固耦合，想着能不能用OpenFOAM+Calculix平替了

[李东岳]

openfoam可以算风机，流固耦合也可以算。所以单独openfoam就可以用。如果你是工程用，那可以用现成的，可能calculix会更适合上手

[HITSC30]

@李东岳 在 流固耦合，换了质量更好的网格流体突然就不收敛了？ 中说：

openfoam可以算风机，流固耦合也可以算。所以单独openfoam就可以用。如果你是工程用，那可以用现成的，可能calculix会更适合上手

初学阶段，商软教程多一点，打算先用商软走一遍，学有余力再进OpenFOAM，感觉是个大工程。谢谢李老师回复！

[Mina Lee]

你好

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次就可以收敛了。流固耦合交界面前几个时间步不收敛很正常，但是你这个是流体不收敛，这就不对了，建议你先单独跑一下流体看看，你这个流体设置不太好。
祝好~