@bestucan 谢谢老师，这个问题我再去认真了解一下并尝试尝试。

@caijinjin 老师，这是我的离散格式和求解器设置

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    limited         cellLimited Gauss linear 1;
}

divSchemes
{
    default                         none;

    "div\(phi,alpha.*\)"            Gauss vanLeer;
    "div\(phir,alpha.*\)"           Gauss vanLeer;

    "div\(alphaRhoPhi.*,U.*\)"      Gauss linearUpwind limited;
    "div\(phi.*,U.*\)"              Gauss linearUpwind limited;

    //"div\(alphaRhoPhi.*,Yi\)"       Gauss linearUpwind limited;
    "div\(alphaRhoPhi.*,(h|e).*\)"  Gauss linearUpwind limited;
    "div\(alphaPhi.*,f.*\)"         Gauss linearUpwind limited;
    "div\(alphaRhoPhi.*,K.*\)"      Gauss linearUpwind limited;
    "div\(alphaRhoPhi.*,\(p\|thermo:rho.*\)\)" Gauss linearUpwind limited;
    "div\(alphaRhoPhi.*,(k|epsilon|omega).*\)"  Gauss linearUpwind limited;

    "div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
}

wallDist
{
    method          meshWave;
    nRequired       yes;
}

求解器设置

solvers
{
    "alpha.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 1;
    }

    droplets
    {
        nCorr                   1;
        tolerance               1e-4;
        scale                   true;
        solveOnFinalIterOnly    true;
        sourceUpdateInterval    20;
    }

    p_rgh
    {
        solver          GAMG;
        smoother        DIC;
        tolerance       1e-8;
        relTol          0.01;
        maxIter         20;
        minIter         2;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "U.*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-7;
        relTol          0;
        minIter         1;
    }

    "(e|h).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1;
        relTol          0;
        minIter         0;
        maxIter         0;
    }

    "(k|epsilon|omega).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-7;
        relTol          0;
        minIter         1;
    }

    "f.*"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0;
    }
}

PIMPLE
{
    nOuterCorrectors    3;
    nCorrectors         1;
    nNonOrthogonalCorrectors 0;
    pRefPoint (0.0276 0 0);
    pRefValue 101325;
}

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

@李东岳 老师，我这个网格是三维的，图上展示了二维，用的结构网格，of9计算。目前发现我只要把流速设的大一点，就不会有这种问题出现，但是流速变小就不行了。另外，网格壁面不加密且数量少的时候它也不报错，自己也还在找原因中。

目前本人在用multiphaseEulerFoam求解器求解孔板中的油水两相流动，孔板架构如图

模拟过程中水是分散相且体积分数设为了0.005，运算时发现有时候运算会报错，计算过程显示Mules算法求解相分数时，不管连续相还是分散相，最小相分数都为0，偶尔为负。请问各位老师和同学们有没有了解这一问题的？不吝赐教。

我自己暂时的理解是，程序中fjfj.phase()/fj.x()和fkfk.phase()/fk.x()指代的是j和k的颗粒数目，也就是右边公式中的Nj和Nk。但是其中的fj.phase和fk.phase是什么含义就有点不太明白了。还有就是聚并速率的量纲是m3/s的，程序里面为什么还要再乘以i颗粒体积fi.x()？这样子的话量纲就不是每秒了。

在离散法求解群体平衡方程中，程序文件中颗粒聚并生成项（birthByCoalescence）和文献中的给的式子形式不太一样，请问哪位大佬可以帮忙解释以下程序中为什么还要多乘以颗粒的直径？

方程里面η在程序里面应该对应gamma，指的是颗粒的分配体积，c对应程序里的coalescenceRate

@李东岳 另外请问分块方式对计算速度有影响吗，我都是用simple方法进行分块的

@李东岳 我台式是AMD 3600 6核， 工作站AMD EPYC 7502 32核，平时用都快不少，就这个算例基本没有加速效果，可能和网格有关系吗？

最近画了个釜式反应器，然后用的是结构化网格，一共50多万个六面体，使用自己台式机分四块并行计算和工作站分24块计算作比较，发现计算速度是一样的，请问这是什么原因？按理说速度最起码要提升一点啊，一点没提升是什么问题？网格是这样比较规整的

@tens 谢谢，这种我也考虑过，但是在进口凸起部分应该会出现涡流，流场结果可能会有偏差。我在尝试用mappeed边界条件来实现。