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

[宝丁]

@东岳 好的，谢谢老师

[宝丁]

@东岳 老师，这个算例运行的时候残差很高，然后到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

[宝丁]

@东岳 还有，老师，您这些都是在哪里看的，感觉从来没见过的用法

[李东岳]

老铁，Minmod是一个特别出名的TVD格式啊。这些都是在书上看到的

[宝丁]

@东岳 好的，谢谢老师，做这个的时间太短了，后面一定加强总结

[宝丁]

@东岳 真心求教这是什么书