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

李东岳

首先你应该试试用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

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

CFD中文网