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

[宝丁]

@东岳 不对，是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格式啊。这些都是在书上看到的

[宝丁]

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

[宝丁]

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

[李东岳]

Finite Volume Methods for Hyperbolic Problems

• 对于刚学CFD的，我不推荐这本书，这本书单纯讲求解$\frac{\p \phi}{\p t}+\nabla\cdot(\bfU\phi)=0$，并且略针对可压缩空气动力学

• 学CFD几年的可以看看，这本书讲解的很细致思路很独特，扩展性很强

[宝丁]

@东岳 收到

[宝丁]

@东岳 老师，国内有没有比较好的讲这种算法的书，我的英语水平如果要看这样的书花费的时间太长了

[李东岳]

国内有没有比较好的讲这种算法的书

@宝丁 那本书即使英语水平好，看起来也费尽，哈哈 目前中文版的我只看过《数值传热学》，但《数值传热学》没有介绍这部分的内容。目前没有其他中文书推荐了。博士的话得加强英语啊。现在好多学生写SCI，我要是改的话基本等于重写。

莫非是硕士？硕士毕业就无所谓了，怎么快怎么来

[宝丁]

@东岳 我是硕士，入门不久，但是也想多学点东西。现在对老师真是钦佩的不行。

[李东岳]

@宝丁 要读博士的话，趁硕士期间多看CFD书，等到博士中后期，好多时间都要用在发SCI写论文上了。现在打基础完全来得及。我当时硕士期间研究方向没定，但比较喜欢推方程，自己硬看了2年CFD教材和C++。CFD这面东西太多了，算法太多，感觉不看个几十本大厚书搞不出个全貌..硕士时间太紧了

不读博士的话，硕士这样直接算毕业完全也没问题。要多学点东西的话，看《数值传热学》，CFD的基本流程也能有点感觉

[宝丁]

@东岳 好的，谢谢老师

[宝丁]

计算虽然可以收敛，但是出口总压却一直不能固定，而且进口马赫数一旦达到0.9，整个计算就会崩掉，这难道是求解器的问题？

[宝丁]

@宝丁 出口静压不能固定

[宝丁]

这应该与我的出口边界条件有关

        /*type           fixedValue;
        //meanValue      $pOut;
        value          uniform 1e5;*/

        type            waveTransmissive;
        field           p;
        psi             thermo:psi;
        gamma           1.4;
        fieldInf        120000;
        lInf            0.2;
        value           uniform 110000;

原先想用fixedMean这个边界条件，但是一用出口边界的压力值是固定了，但是出口附近的压力会异常的高，好像出口把流道给堵住了一样；
但是用 waveTransmissive这个边界条件，边界的压力值就会一直下降，进口的压力也下降，速度升高，Ma数达到0.9左右整个计算就会崩掉。

[宝丁]

目前来看这个问题是不是已经要修改ＯＦ的一些底层代码
不知道这个难度如何，应该要放弃这个算例了