计算圆管绕流，的流固耦合fsiFoam的Floating point exception (core dumped)

[炎冰灬灬]

@bestucan 我的耦合求解器设置的是e-6，流体也是e-6，固体是e-9。是不是固体这边设置的残差太小了呢？下面的截图依次是流体、耦合和固体的残差设置。

[bestucan]

@炎冰灬灬 我没做过，所以不清除多大合适。但是你说每次都计算到30次才进入下一个时间步是因为这个设置的。设大了，就不会要计算到30步才停了。或者把最大迭代步设大点。

判断当前的残差还需不需要算，需要看残差图，看看随着迭代的增加残差是否还明显下降，如果残差线都平了，再算意义也不大了。如果还有很大的下降空间，再算算稳一些。

[bestucan]

OpenFOAM生成残差图，好像之前有人在论坛里讨论过

[炎冰灬灬]

@bestucan 我设置的计算时间是6秒，前面都还好，算得很快，每次算到2秒的时候，就开始每个时间步迭代到30次了

[bestucan]

@炎冰灬灬 八成是算到两秒的时候，物理量开始急剧变化了（比如撞倒边界突然减速之类的），所以追踪起来比较费劲。那就让他算到30步。等后面这种变化平稳了，就不会每一步都要算到30步了

[炎冰灬灬]

@bestucan 我算的是柔性翼形。我看自带算例里面，大概在两秒的时候，自带算例里面的固体开始出现变形。所以我在考虑两秒的时候是不是就是固体开始变形的时候呢，然后这个时候有什么东西导致了迭代出现30次

[炎冰灬灬]

@bestucan 自带算例里面也是，从某一时刻开始一直到计算结束，每个时间步迭代为16次

[bestucan]

@炎冰灬灬 反正不崩了，能算下去就是希望₍₍ (ง ˘ω˘ )ว⁾⁾
先过一遍再慢慢改善

[炎冰灬灬]

@bestucan 有没有可能是差分格式或者是求解格式的问题呢

[bestucan]

@炎冰灬灬 判断差分和求解格式需要的理论深度比较深，我说不来。

我目前的认知：如果是差分的问题，可能会造成特殊的结果，不符合实际（比如中心差分的棋盘压力，比如迎风算带回流的）。如果是求解格式，它会影响能不能收敛或者收敛的速度。

可是的你的计算很正常，只是收敛的慢些，还不至于出错。可能会有另一种求解格式能快一些收敛。但是差分格式有没有问题要看物理场合不合预想，结果合不合实际。

如果你那么在意单时间步内收敛步数。那就再缩短时间步。时间步越短，单时间步收敛的越快。当然计算同样的时间，需要的时间步越多。

[炎冰灬灬]

@bestucan 谢谢，涨知识了。我在考虑要不要最大迭代次数变成50次，看看它是不是要迭代到50次，如果迭代不到50次，那就应该是30次不够收敛，如果迭代到了50次，就有点问题了感觉

[炎冰灬灬]

@bestucan 我设置迭代50次，它还真的就迭代到了50次。。。

[bestucan]

@炎冰灬灬 所以要看残差线，如果随着迭代残差不再下降，迭代到天荒地老也不行啊

[炎冰灬灬]

@bestucan 对的，我正在监测，我看了自带算例的残差，0-2秒残差下降，2秒这一刻突然集体升高，然后2-4秒一直上升，4-6秒上下波动趋于稳定。然后我现在主要就是两秒这个时候在那疯狂迭代，我两秒之前和自带算例一样，残差是下降的