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

[宝丁]

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

[宝丁]

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

[宝丁]

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

[李东岳]

cascade.zip

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

不要运行blockMesh，忘记删了

[宝丁]

@东岳 好的，谢谢老师

[宝丁]

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

[李东岳]

@宝丁 降低时间步长试试

[宝丁]

@东岳 好的

[宝丁]

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

[宝丁]

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

[李东岳]

@宝丁 多谢分享 楼主好人

[宝丁]

@东岳 额 还没有算收敛，这个好难啊

[李东岳]

用rhoPimpleFoam试试？

[宝丁]

@东岳 好的

[宝丁]

@东岳 老师，不过我在用rhoSimpleFoam中遇到了一个很奇怪的现象，就是进口的压力会出现异常，使得整个计算发散，就像这样。这又会出现数值上面的异常。可以用层流模型试一下。这个问题也是困扰了我好久

[宝丁]

@东岳 不对，是rhoCentralFoam

[宝丁]

@东岳 老师，我已经调通了，非常感谢老师的知道。具体做法就是只修改

    default         none;
   Gauss limitedUpwind limited
    div(phi,e)      Gauss limitedUpwind limited;
    div(phid,p)     Gauss limitedUpwind limited;
    div(phi,K)      Gauss limitedUpwind limited;
    div(phiv,p)     Gauss limitedUpwind limited;
    div(phi,k)      Gauss upwind;
     div(phi,omega)  Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
    div((muEff*dev2(grad(U).T())))  Gauss linear; 

注意，现在计算出来的结果并不是很理想，感觉粘性太大了，所以想问一下东岳老师

Gauss limitedUpwind limited

 wallDist  {method Poisson ；}

这两项的含义
附一张计算的结果

[宝丁]

非常非常非常感谢东岳老师的帮助

[李东岳]

Gauss linearUpwind limited:

• Gauss: 高斯积分计算散度

• linearUpwind: 线性迎风格式

• limited: 在线性迎风格式需要指定U的梯度怎么计算，limited用来调用U的梯度计算格式

wallDist  {method Poisson ；}

是用来计算壁面距离的，Poisson是一种计算避免距离的方法，比较适用于复杂几何，meshWave方法适用于简单几何。kOmegaSST模型会用到壁面距离。

[宝丁]

@东岳 哦哦，不知道怎么复制错了，是那个midmod