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

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

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

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

就是一个这样的结构，望各位大佬指点一二！

现在想模拟一个罐内多相流的循环，就是流体从底部出口流出后从上面的入口重新进入，来模拟泵的循环，目前进出口面积相等的案例已经使用cyclic边界加fvoptions模拟出来了，结果没啥大问题，质量都守恒。

但是如果我想让上面的进口大一点，下面保持不变，那cyclic边界就不能使用了。是否有其他边界边界可以代替？cyclicAMI边界是否可以实现？使用周期边界的话，进出口大小不一样是否会导致进多出少，质量的不守恒？

@李东岳 好的，谢谢东岳老师

@李东岳 用的of7，那请问我这个网格的非正交性算大还是正常，如果网格质量不好那使用reactingTwoPhaseEulerFoam来计算是否会网格的敏感性有所改善？还是用最新的of9里面的multiphaseEulerFoam会好一点？

另外附上我的离散格式，出现负值是否因为离散格式设置错误？未包含有界限制？

ddtSchemes
{
    default     Euler;
}

gradSchemes
{
    default     Gauss linear;
}

divSchemes
{
    default                         none;

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

    "div\(alphaRhoPhi.*,U.*\)"      Gauss limitedLinearV 1;
    "div\(phi.*,U.*\)"              Gauss limitedLinearV 1;

    "div\(alphaRhoPhi.*,(h|e).*\)"  Gauss limitedLinear 1;
    "div\(alphaRhoPhi.*,K.*\)"      Gauss limitedLinear 1;
    "div\(alphaPhi.*,p\)"           Gauss limitedLinear 1;

    div(alphaRhoPhi.particles,Theta.particles) Gauss limitedLinear 1;

    "div\(alphaRhoPhi.*,(k|epsilon).*\)"  Gauss limitedLinear 1;

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

    div((((thermo:rho.particles*nut.particles)*dev2(T(grad(U.particles))))+(((thermo:rho.particles*lambda.particles)*div(phi.particles))*I)))  Gauss linear;
}

laplacianSchemes
{
    default     Gauss linear uncorrected;
    bounded     Gauss linear uncorrected;
}

interpolationSchemes
{
    default     linear;
}

snGradSchemes
{
    default     uncorrected;
    bounded     uncorrected;
}

如图，是一个搅拌釜内的液固两相流，用的twoPhaseEulerFoam,计算过程中逐渐发散，出现了很多负值，请问这个可能是什么原因导致的？对网格进行了质量检查，网格为非结构网格，这个非正交性的值是偏大还是正常范围？

fvSolution设置为：（为了不求解温度，我将h/e的迭代数设为0，是否对计算有影响？）

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

        smoothLimiter   0.1;

        implicitPhasePressure yes;
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-9;
        relTol          0;
        minIter         1;
    }

    p_rgh
    {
        solver          GAMG;
        smoother        DIC;
        tolerance       1e-8;
        relTol          0;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

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

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

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

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

PIMPLE
{
    nOuterCorrectors 3;
    nCorrectors      1;
    nNonOrthogonalCorrectors 0;
}

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

希望各位大佬帮忙解答，感谢！