基于interCondensatingEvaporatingFoam求解器，我加入了多孔介质模型。
第一个是stefanProblem，两个case的差别是孔隙率和渗透率的取值不同。可以看到当孔隙率取1（没有多孔介质），它的最大液相分数是可以到1的（不知道这个是不是可以保证我的求解器是对的），但是第二个case的最大液相分数是在0.6左右。


第三个case结果中液相百分数甚至只有0.02~

下面是我fvsolution中关于alpha的约束，请问我应该怎么查找原因。感谢各位老师~

http://blog.sina.com.cn/s/blog_5fdfa7e601010rkx.html 苏老师的博客

基于interCondensatingEvaporatingFoam求解器，编写了一个可求解porous区域的求解器。在我设置孔隙率eps为1时（即没有多孔区域，用以和标准的stefanProblem结果做验证），求解的结果和stefanProblem结果相吻合。在我设置孔隙率eps为0.8/0.6/0.5/0.2时，求解结果均发散。
我试图修改了fvSolution中的参数，将relaxationFactors从1逐渐降低到0.1，发散的时间点也从1.6s推迟到了4s，这和relaxationFactors发作用相符合。将p_rgh的solver从PCG改为GAMG，并没有推迟发散时间，反而求解速度更慢
请问我还应该注意哪几个方面用来减小发散的可能性。感谢各位老师~

solvers
{
    alpha.liquid
    {
        nAlphaCorr      2;
        nAlphaSubCycles 2;
        cAlpha          1;

        MULESCorr       no;
        nLimiterIter    5;
    }

    rho
    {
        solver           diagonal;
        tolerance       1e-7;
        relTol          0.1;
    }

    rhoFinal
    {
        $rho;
        tolerance       1e-7;
        relTol          0;
    }

    p_rgh
    {
        solver           GAMG;
        tolerance        1e-9;
        //solver          PCG;
        preconditioner  DIC;
        relTol           0.001;
        smoother         DIC;
    }

    mDotSmearFinal
    {
        solver           PCG;
        tolerance        1e-6;
        preconditioner   DIC;
        relTol           0.00;
        smoother         DIC;
    }

    p_rghFinal
    {
        $p_rgh;
        tolerance        1e-9;
        relTol           0;
        minIter 10;
    }

    pcorrFinal
    {
        $p_rgh;
        tolerance        1e-9;
        relTol           0;
    }


    "(U|h|T.*|k|epsilon|R)"
    {
        solver           smoothSolver;
        smoother        symGaussSeidel;
        tolerance        1e-7;
        relTol           0.;
        minIter          15;
        maxIter          50;
    }

    "(U|h|T.*|k|epsilon|R)Final"
    {
        $U;
        tolerance        1e-7;
        relTol           0;
        maxIter          50;
    }
}

PIMPLE
{
    momentumPredictor   yes;
    nCorrectors         5;
    nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
    equations
    {
        ".*"           0.1;
    }
}

@李东岳 感谢老师 发现是我Co计算公式写错辽~已经修改并且成功运行了 谢谢~

补充~
算例几何为规则长方体（前后面为empty），10mm✖0.2mm✖0.02mm，划分网格数为400✖2乘以1。
用相同的model使用不同的求解器，就不会出现这样的情况，库朗数是在逐渐增大，Δt也越来越大，和想象的情况一样。

请问我在跑算例的时候，存在一个界面库朗数的突然增大，之后又稳定在一个值。导致我的deltat变得非常小（负十五次方），请问这可能是什么原因导致的。运行的log文件我附在下方~感谢各位老师~

log.doc

@bestucan 终于解决啦 谢谢老师 感谢分享~

已经编译成功。错误的原因是，有其他的求解器中有相同的库，把其他的求解器编译的删除掉（因为那些求解器只是过渡 所以就选择删除了）或者命名为不同的库 可以减少这类错误的发生。这也是为啥放在of源代码的位置上就没有所述错误的原因。
不过编译结束之后，会出现Duplicate entry错误，还没有解决这个问题，希望能得到各位老师的指导。谢谢~

上面的办法并不可行~会出现Duplicate entry constantAlphaContactAngle in runtime selection table fvPatchField 等类似的错误。还是应该放在USER下

@hongjiewang 已解决~不过我是把所有的库都编译在了FOAM_LIBBIN和FOAM_APPBIN下，所以应该是我的编译路径有问题~