OpenFOAM SRF求解器有大问题呀？

[东方白杨]

最近一直在尝试用SRFpimpleFoam求解器算一个旋转槽道流，如图1，可是发现计算一直发散，经过再三尝试，目前我测出来的结论如下：

计算参数，进口主流0.1335m/s，旋转轴为Z轴（0 0 1），旋转中心为原点（0 0 0），转速rpm=120，槽道几何比例为4:2:2，我尝试了不同scale下旋转槽道的计算，发现有问题。

首先，计算域为0.04m-0.02m-0.02m时，旋转状态下是可以算的，监测通道进口、中心点相对速度，都可以稳定。但是当计算域为4m-2m-2m时，通道进口、中心点相对速度就无限增大，然后发散了。

经过各种测试，排除原因，我得到一个结论，就是这种情况是离心力导致的，SRF求解的动量方程如下：

离心力已经写成势函数形式，我发现上面的计算发散就是由于离心力太大导致的，经过测试，有三种情况会导致发散，
1转速omga很大
2旋转半径很大，就如前面的0.02m变为了2m的旋转半径，导致离心力很大
3旋转中心里通道进口很远，如果我把旋转中心从0点调为(-5 0 0)，它也会发散
综上原因，就是当离心力很大的时候，难道计算就无法稳定了吗？？？

[aurora16]

网格？边界层？时间步长？

[东方白杨]

@aurora16 这些因素已经排除过了，况且你说的这几个点也不会导致我上面提到的问题

[aurora16]

你进出口设的周期吗

[Amadeus]

源项太大了又是显式离散可能导致对时间步长(库朗数)的要求很严格?

[东方白杨]

@Amadeus @aurora16 感谢各位回复，最后发现是由于对流-扩散方程中，离心力项太大了(高达几百量级)，而其它项还在10的量级，离散求解方程时，量级差异太大了，导致发散。计算还是要保证方程中各项量级差异不能太大