为什么计算过程中动量方程U不收敛？

[Jimmy Yao]

各位大佬好，我又来提问个问题了
我在计算一个小薄壁件的流固耦合问题，之前主要是流体出问题，通过不断调整网格（更换了N多软件，最终确定了用starccm+的切割体网格，属实真爱了）

现在质量应该是非常好了，checkMesh也通过了
但是从输出的残差情况来看，压力泊松方程计算貌似还不错，但动量方程好像是不收敛，请问这种情况下我应该怎么做算例的修改？
下面附上log文件、边界条件以及计算设置，麻烦各位大佬帮忙出出主意，谢谢啦！！

3d540928-8b55-4848-8aed-39e35aa518dd-log(1)
045f87f1-bf52-43a9-92e4-3710cc8b0b9b-controlDict
f2a44212-72f0-4c08-bf85-d260e5910237-fvSchemes bc9f46a9-5381-4120-ab33-be2f6fded262-fvSolution
c4200769-ceed-46d3-a950-9dcfd717ef68-0.rar

[李东岳]

为什么设置这么多的pimple循环？

[Jimmy Yao]

@李东岳 李老师，我是参考的别人的流固耦合算例做的设置，其实这么多确实没啥用，算的时候懒得改了

[李东岳]

首先你应该试试用piso同时减小时间步长

[Jimmy Yao]

@李东岳 好的李老师，我现在试一下

[Amadeus]

流固耦合问题你的固体是有限元法计算吗? 用什么软件和OpenFOAM配合呢? 是不是考虑流固界面上插值计算导致的速度场难以收敛?

[Jimmy Yao]

@Amadeus 大佬您好，我是用的calculix做耦合，之前不收敛是流场网格的问题，在画完网格之后想试一下流场行不行，现在是在单纯算流场
刚刚根据李老师建议改成PISO，暂时看能往下算，等我再算一段时间看看效果

[Jimmy Yao]

@李东岳 感谢李老师！PISO+小时间步长确实是可行的，但是我对这个还有一些疑问，希望您那有时间能帮忙看看。为什么输出的收敛信息并没有达到fvsolution的要求，pisoFoam就判定收敛了？

同时我后面又试了下小时间步长下的PIMPLE，也可以算，但是不知道为啥pimpleFoam必须算到最大循环次数。而且我之前计算时候都会显示“PIMPLE：converged in X ierations”，但是这个算例没有，从残差上看好像还收敛了

然后计算结果也很怪，我照常给了1m/s的速度入口、0压力出口，但是在第一个时刻整个流体域内都变成了1m/s，如下图9所示，是第一个时间步的速度U计算结果。通过查阅料（https://www.zhihu.com/question/405674832），我发现是不可压算法的问题的，但是之前的算例都没有发现该问题，这种情况我该如何进行算例修改？是修改边界条件吗？

[Jimmy Yao]

下面是我之前的边界条件设置

上面的知乎链接好像不太好打开，我放在下面
https://www.zhihu.com/question/405674832

[李东岳]

为什么输出的收敛信息并没有达到fvsolution的要求，pisoFoam就判定收敛了？

为什么说没有达到？我看已经收敛了

[Jimmy Yao]

@李东岳 哦哦，李老师我想错了，瞬态问题是应该看final residual是吧，没琢磨明白

还有就是这个气体突然充满的问题，我又看了下模型，一开始我以为模型尺寸的问题，但用paraview看了下好像导入openfoam没变化，就不知道这个怎么解决了

[李东岳]

所有问题都是initial residual

还有就是这个气体突然充满的问题，我又看了下模型，一开始我以为模型尺寸的问题，但用paraview看了下好像导入openfoam没变化，就不知道这个怎么解决了

这个我不知道你说的是什么

另外从你的问题历史记录来看，你应该刚上手openfoam就开始了你的课题的计算，感觉基本功不太牢固。你最好把基本的openfoam应用的面铺开，多跑跑基础算例，看看都是怎么设置的，各种设置对结果有什么影响，有个初步了解之后再上手你复杂的各种耦合算例。

[Jimmy Yao]

@李东岳 好的，谢谢李老师！

我刚说的问题其实就是这个，我设置的入口是1m/s，按道理应该是速度一点点进入流域，但是第一个时刻就充满整个计算域了，这个不知道咋搞
下面就是第一个时刻的计算结果，步长是1e-5

[李东岳]

https://cfd-china.com/topic/7872 你跟这个是一个问题吧？这是正常的

[Jimmy Yao]

@李东岳 好的，谢谢，李老师！