CFD中文网

@yfclark 非常感谢。

请问paraview能够统计网格数吗？比如液相体积分数在0.5-0.6的网格数有多少。或者其它软件或者程序能够计算出来的方法，有了解的大神指点下，谢谢。

@玉玺 有可能是湍流模型选择的问题，你可以试试其它湍流模型，另外，这是用于多组分的，不但要设定相体积分数，还要设定组分的质量分数场。

@玉玺 因为我不做反应，后来就没有深究，不好意思，我也给你不了其他建议了。

我记得这是求解器设定的限制，好久没用那个了，具体忘了。建议你具体看一下模型，我当时的设置就是跟所选模型不符。

是的，求解器自带的是蒸发，我想计算空化，二者的计算原理大不相同，我更改了质量源相，导致质量分数越界。

@东岳 我看MULES代码，感觉他就是为求解相方程定制的，没有扩散相的处理，所以我不知道这里的扩散相放在哪个位置，怎么处理。

@东岳 那么下面带拉普拉斯项的组分输运方程如何采用MULES求解呢？谢谢。

        fvm::ddt(alpha, rho, Yi)
      + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")

      - fvm::laplacian
        (
            fvc::interpolate(alpha)
           *fvc::interpolate(this->muEff()/Sct_),
            Yi
        )
     ==
        alpha*this->R(Yi)
       + Source
      + fvc::ddt(residualAlpha_*rho, Yi)
      - fvm::ddt(residualAlpha_*rho, Yi)

@yfclark 非常感谢你的建议。源相我确实是直接添加的，没有再做处理，这个Sp会有什么作用（虽然有时我也会模仿着这样处理）？越界的主要原因是什么？

@东岳 对，就是自带的Y方程。

@东岳 非常感谢东岳老师的回答，我有个疑问，这个求解器是可以算蒸发的，而且用的就是界面组分，他们是怎么处理的呢？就避免了越界。

reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界，这是什么原因，为什么没添加源相之前有界呢？如何处理才能避免质量分数远大于1？

@yfclark 源相是需要加到压力泊松方程里面，当然由于源相的产生导致的动量源相和质量源相也要加到对应方程中的。不设置任何东西应该是不考虑曳力模型的。不清楚你要怎么加拉格朗日模型，为什么不考虑相间曳力？

@东岳 那跟商软的层流模型有区别吗？商软也是这样?谢谢！

@Hungryandfool 刚才发错了，我用的社会上的超算，0.1元每核时，我们学校也是这样的。价格基本都一样。

@zwl 刚才发错了，我用的社会上的超算，0.1元每核时，我们学校也是这样的。价格基本都一样。

@队长别开枪 非常感谢，我看了redistributePar的help，看到这个命令选项里确实有并行，我研究下怎么用。

@队长别开枪 非常感谢，我看了redistributePar的help，看到这个命令选项里确实有并行，我研究下怎么用。

@dzw05 不知道你说的直接处理是在processor路径下直接paraFoam吗？这样只能看到分得这一块网格数据，而我要看整体结果，怎么操作？有相关资料吗？求分享，谢谢。

请问下东岳老师，openFoam里面的求解器如果选成层流，是不是都是直接求解动量方程？可以作为直接数值模拟求解器？