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

[宝丁]

@东岳 我在这里找到了一个您的邮箱，给您发了一个这个网格量一般的网格，您能收到吗？
振荡是从一开始就振荡。时间步长为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 采用速度入口算出来的结果，我感觉挺好的

[宝丁]

@东岳 这个是用您的对流项格用比较密的网格算出来的结果，稳定了，但是求解的流场还有一些不太理想，我再好好调一下，感谢东岳老师。

[宝丁]

@东岳 东岳老师，感觉问题并没有解决，我用upwind的算法求解稀疏一些的网格算的久一些就会出现发散的情况，网格量大的却能收敛，收敛后的流场如上所示也会有一些不好的地方，现在尝试提高离散项精度，但是很难找到合适的离散项，有些不知所措了，产生这个的根本原因是什么。

[李东岳]

cascade.zip

这是我最后得到的结果，不知道准不准，从下面那个图看起来差不多，你可以顺着我的这个算例，逐个调整参数直到准确位置

不要运行blockMesh，忘记删了

[宝丁]

@东岳 好的，谢谢老师

[宝丁]

@东岳 老师，这个算例运行的时候残差很高，然后到0.01s 左右就会发散了，不知道是不是网格的问题

[李东岳]

@宝丁 降低时间步长试试

[宝丁]

@东岳 好的

[宝丁]

@东岳 @东岳 好像也没有作用，我也用我的网格试了一下，只改网格，出口那个地方也是出现了不稳定现象，但是还能算，感觉收敛不了

[宝丁]

@东岳 可以尝试把速度 压力 温度的离散格式全改成二阶迎风，目前看效果挺好的