求助: SonicFoam求解二维叶栅时无法收敛(求解文件可下载)
-
源文件里面的松弛因子给的不合理,修改以后依然没什么效果
修改后relaxationFactors { fields { p 0.7; rho 0.8; } equations { p 0.7; U 0.3; "(e|h|k|epsilon|omega)" 0.2; }
-
现在我能确定的是我的整个边界条件在低速下是可以用的,用pisofoam,湍流模型选用kepsilon收敛的很好,现在正在尝试kOmega。但是在高速求解器如我现在要用的sonicfoam 和原先尝试的rhoCentralFoam都会出现问题。在sonicfoam 里面的问题就是会出现内部流场十分紊乱,所以我觉得是不是我求解和算法控制那块是不是有问题
-
我也尝试过,确实是可以的,但是我的工况要求我定义进口的方向,尤其是这种跨声速流场,对方向很敏感
-
-
@东岳 好的 文件不大,就直接放到这里了
R67_70.zip -
速度角度和网格不成直线,对流项格式不好设置。不太清楚你这个结果应该是什么样,最后我这个稳定到这样了。因为我重新画了网格,你用你的网格试试下面的fvScheme,我感觉可能是格式的问题,因为看起来是震荡。
还得进一步研究一下,不能用这么迎风的格式,你用稳态求解器试过没
ddtSchemes { default Euler; } gradSchemes { default cellLimited leastSquares 01; } divSchemes { default none; div(phi,U) Gauss upwind; div(phi,e) Gauss upwind; div(phid,p) Gauss upwind; div(phi,K) Gauss upwind; div(phiv,p) Gauss upwind; //div(phi,U) Gauss limitedLinearV 1; //div(phi,e) Gauss limitedLinear 1; //div(phid,p) Gauss limitedLinear 1; //div(phi,K) Gauss limitedLinear 1; //div(phiv,p) Gauss limitedLinear 1; div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.5; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; }
2019年4月12日 11:59
11/72
2019年5月20日 03:01