关于multiphaseEulerFoam运算分散相体积分数问题

Caijinjin

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

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

李东岳

openfoam9这面现在multphaseEULERfoam比之前稳定多了。你看看是不是你设置的有问题。另外你这个是结构网格吧？

Caijinjin

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

Caijinjin

@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;
    }
}

bestucan

你看你图中圈着的上一个Phase-sum volume fractionmin 是正的，圈着的地方变成负的了，当然max那上下都负的。max先数值溢出，min跟着溢出。

说明这是计算本身的不稳定。

速度快了，能够不发散，是把不稳定多出来的那部分赶紧吹出去了；慢了，吹不出去就原地爆炸了。

我看了格式，里面有不少迎风格式，好像迎风格式不能考虑回流？
这种管子冲快了，刚出狭管上下两个角落应该会有回流涡吧？
可能因为这个

Caijinjin

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