rhoPimpleFoam在大网格下的发散问题

[mathea]

各位大佬好，我这边用openfoam10的rhoPimpleFoam算了一个亚音速马赫数的叶珊，通过调整测试，我下面的这些设置分别在15万网格的2D算例和200万网格的测试算例上面都能够稳定运行，但是在大网格1000万的网格级别就会发散，1000万网格计算的初始条件是200万网格的mapFields后开始计算，计算了大概9万步后发散。下面会有对比图，我觉得可能在涡的附近算的有问题，但是怎么修改这个问题，想请各位前辈大佬们指点。
case.zip

[李东岳]

你是否确保时间不长足够小

[mathea]

@李东岳 在 rhoPimpleFoam在大网格下的发散问题 中说：

你是否确保时间不长足够小

李老师好，时间步长我试过Co=0.2或者Co=0.1.我也试过给固定的deltaT=1e-7，都会形成同样的结果。

[李东岳]

你看一下你的湍流粘度是否正常
另外这个算例你是不是也看过 http://dyfluid.com/openfoam.html#id8

[mathea]

湍流黏度看起来不正常的原因是要保证跟实验和仿真的雷诺数相同。这个算例我是参照过的。

[李东岳]

湍流黏度看起来不正常

你湍流变量的边界条件调整一下，把这个搞正常看看

[mathea]

@李东岳 在 rhoPimpleFoam在大网格下的发散问题 中说：

湍流黏度看起来不正常

你湍流变量的边界条件调整一下，把这个搞正常看看

调整后跑了十来天，现在能够稳定计算，特来感谢李老师！李老师妙手回春！

[mathea]

李老师，上面算的算例因为网格问题要重新计算了，我参照您发我的算例在ICEM中画网格。如图所示A的严格对称点应该是B点，您画在了C点，所以要在openfoam里面用cyclicAMI来周期上下两个面。出现了面不匹配的老问题，我尝试在ICEM中设定周期，但是好像这种叶珊上下面应该没法设置节点的周期把？随后我仔细观察了您的算例，发现在boundary文件中的cyclicAMI的容差设置为1，随后我在createPatchDict采用了相同的容差。这一步算是可以通过了，但是在decomposePar中采用scotch方法进行分块时，还是在cyclicAMI这里报错了，麻烦李老师给看看我上述的操作是漏了什么嘛？
另附一下 网格处理命令：
fluentMeshToFoam ICEM.mesh
autoPatch 45 -overwrite
createPatch -overwrite