OpenFOAM中的方程类为fvMatrix，该类为一个基于有限容积的稀疏矩阵类，更确切的说它是一个类模板。对于定义标量方程 fvMatrix<scalar>，也可以写成fvScalarMatrix，两者一样，因为openfoam运用了 typedef进行了类型别名定义
typedef fvMatrix<scalar> fvScalarMatrix ; //标量
typedef fvMatrix<vector> fvVectorMatrix ; //矢量
fvMatrix中有两个基本变量A 和 b初始化为0，用以存储系数矩阵A和右边值b。

@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方程。