多相流射流求解

yfclark

各位大佬好，我最近正在用compressibleInterFoam计算射流问题，求解出来的相体积分数分布如下：

感觉分布不是很正常。
我的fvScheme设置如下：

对流项和扩散项：
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss linear;

    div(rhoPhi,U)   Gauss linearUpwindV grad(U);
    div(rhoPhi,T)   Gauss linearUpwind grad(T);

    div(phi,p)      bounded Gauss Gamma 0.2;
    div(rhoPhi,K)   bounded Gauss Gamma 0.2;

    div(phi,k)         bounded Gauss Gamma 0.2;
    div(phid1,p_rgh)   bounded Gauss Gamma 0.2;
    div(phid2,p_rgh)  bounded Gauss Gamma 0.2;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
    //turbulence model
    div(rhoPhi,k)                      bounded Gauss Gamma 0.2;
    div(rhoPhi,epsilon)                bounded Gauss Gamma 0.2;
    //other turbulence model
    div(rhoPhi,R)                      bounded Gauss Gamma 0.2;
    div(rhoPhi,omega)                  bounded Gauss Gamma 0.2;
    div((rho*R))                    Gauss linear;
    div(R)                          Gauss linear;
    div(U)                          Gauss linear;
    div((muEff*dev(T(grad(U)))))    Gauss linear;
div((muEff*dev2(T(grad(U)))))   Gauss linear;

边界条件：
入口速度入口压力零梯度
出口为压力无反射边界条件

yfclark

速度云图

ldymig

时间算的应该还不够久，还没有充分的发展

ldymig

此处的压力无反射边界条件是哪一个？