极低雷诺数下，PIMPLE迭代初始残差保持为1，降不下来

[李东岳]

今天才看到这个消息。不知道后来debug的怎么样了。

[lwjetmann]

@李东岳 最后还是没搞明白咋回事，怀疑是动网格不能用来算很低雷诺数的情况，外边界条件可能不匹配

假如圆柱不动，多低雷诺数（甚至拿掉对流项）应该都可以算吧

[李东岳]

你上次发帖子的时候我刚办完课，然后去美国了。
这次回帖子的时候我下周要回国办课了。
两次都赶上我要么出国，要么回国。4月中旬从英国回来，或许可以给你看看。

[lwjetmann]

@李东岳 感觉这个问题挺有意义的，肯定可以从原理上说清楚为啥发散。给的是一个极大的单块网格平动（入口速度零，零法向压力梯度；出口压力零，零速度梯度），一直调不出来就放弃了，可能雷诺数低到某个程度，就不能靠放大计算域来使边界条件匹配了。

祝李老师上课愉快！

[CHRISVC]

最近也在搞类似的，debug了很久，不过没有研究过1e-3这么低的雷诺数，提供几个思路：

1. 关闭动网格，只算静止圆柱的看看还发散吗？确认一下是不是dynamicMeshDict设置的问题；

2. 雷诺数降这么低，对y+的影响是怎么考虑的？是否要改壁面函数或不能用壁面函数了？

[AppleKiller]

@lwjetmann 在 极低雷诺数下，PIMPLE迭代初始残差保持为1，降不下来 中说：

@李东岳 感觉这个问题挺有意义的，肯定可以从原理上说清楚为啥发散。给的是一个极大的单块网格平动（入口速度零，零法向压力梯度；出口压力零，零速度梯度），一直调不出来就放弃了，可能雷诺数低到某个程度，就不能靠放大计算域来使边界条件匹配了。

祝李老师上课愉快！

老哥，这个问题最后解决了吗

[lwjetmann]

@CHRISVC 感谢回复！

1. 观察到：关闭动网格不发散，减少圆柱数量（降低网格复杂度）发散的雷诺数会降低，怀疑是动网格或者网格非正交修正的问题

2. 层流不用壁面函数

[lwjetmann]

@AppleKiller

没有

做了一个能复现问题的干净算例，感兴趣的同学可以试试（网格量57万，需要算一段时间才发散，我保留了发散后的log.pimpleFoam文件，可以用文件里的plot_residual脚本来查看残差。使用的OpenFOAM版本为v2006或v2406）。

通过网盘分享的文件：2D_21cylinders.zip
链接: https://pan.baidu.com/s/1lHExA2yRnUzXFxYYtAc6rw?pwd=8888 提取码: 8888

[李东岳]

有可能把算例网格降低到5万以下么

[lwjetmann]

@李东岳 很奇怪，又画了个很稀的6万网格量的网格，先前会发散的极低雷诺数，现在却不发散了。看来雷诺数极低，网格就不能太密？

通过网盘分享的文件：2D_21cylinders_Re0.0001.zip
链接: https://pan.baidu.com/s/1NdfFZrV2k_T8mAY7fA7ktg?pwd=8888 提取码: 8888