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),它也会发散
综上原因,就是当离心力很大的时候,难道计算就无法稳定了吗???