求助： SonicFoam求解二维叶栅时无法收敛（求解文件可下载）

李东岳

不收敛可能性太多了，周一给你调试下你的算例，现在手头没有openfoam

既然rhoCentralFoam能收敛，为何不用rhoCentralFoam

宝丁

@东岳真的是太感谢您了，rhoCentralFoam的总温总压边界再配上进口速度方向，会使得计算发散，我现在正在想办法解决。在压气机内流叶栅的计算里面，一般都会给进口总温总压边界条件，和实验对比的比较好

宝丁

源文件里面的松弛因子给的不合理，修改以后依然没什么效果
修改后

relaxationFactors
{
    fields
    {
        p 0.7;
        rho 0.8;
    }
    equations
    {
        p 0.7;
        U 0.3;
        "(e|h|k|epsilon|omega)" 0.2;
    }

李东岳

你的算例我下来之后，发现lowRe这个计算没有问题？

宝丁

@东岳计算是要过一段时间才会显示发散，看云计算流场应该是乱的

宝丁

@东岳这是5e-5时刻的流场云图

宝丁

现在我能确定的是我的整个边界条件在低速下是可以用的，用pisofoam,湍流模型选用kepsilon收敛的很好，现在正在尝试kOmega。但是在高速求解器如我现在要用的sonicfoam 和原先尝试的rhoCentralFoam都会出现问题。在sonicfoam 里面的问题就是会出现内部流场十分紊乱，所以我觉得是不是我求解和算法控制那块是不是有问题

李东岳

@宝丁这个看起来出现了震荡，我在测试，就是算起来太慢了，

宝丁

@宝丁预计进口流速为60m/s，还是总温总压入口

宝丁

@东岳我在这里找到了一个您的邮箱，给您发了一个这个网格量一般的网格，您能收到吗？
振荡是从一开始就振荡。时间步长为1e-8或者5e-9时，计算到0.0001就可以看到不稳定的流场。

李东岳

我把你的算例，把速度进口替换成zeroGradient计算就没问题了，你试试

但你这个是不是应该是cyclic边界？另外，即使用cyclicAMI，是不是应该是rotational的？
另外，网格正交性是不是可以调整下。个人感觉你调试这个算例，不用满足y+，用3000 5000调试就行，要不算起来太慢了

宝丁

我也尝试过，确实是可以的，但是我的工况要求我定义进口的方向，尤其是这种跨声速流场，对方向很敏感

宝丁

我的边界条件是平移周期边界，cyclicAMI和cyclic相比只是不需要交界面两边网格一一对应

李东岳

把你的tin和blk发来看看，我调整下网格，

宝丁

@东岳当时试的时候一开始确实没什么问题，可是算到后面zeroGradient 边界条件也会发散，但是会好很多，不知道是不是几何不合适，就没有往下尝试

宝丁

@东岳好的文件不大，就直接放到这里了
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;
}

宝丁

@东岳好的，非常非常感谢，我马上就去尝试

宝丁

@东岳我用稳态算过，也是发散

宝丁

@东岳这个是用rhoCentralFoam 采用速度入口算出来的结果，我感觉挺好的

CFD中文网

求助： SonicFoam求解二维叶栅时无法收敛（求解文件可下载）