scalarTransportFoam计算某一相的停留时间分布

[青山]

首先使用multiphaseEulerFoam模拟了一个气液固三相连续搅拌釜，想获得固体相在釜中的停留时间分布，所以想将固体的速度场和通量场给到scalarTransportFoam中，从入口进入示踪剂以获得固体停留时间分布。这时遇到一个问题，整个釜的空时大概4个小时，但是计算大概1000s之后整个釜内的示踪剂浓度分布就不变化了，出口的浓度也达到稳定并且极低，请问这是什么问题？附上我的设置

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,T)     Gauss linearUpwind grad(T);
}

laplacianSchemes
{
    default         none;
    laplacian(DT,T) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

solvers
{
    T
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 1;
}
relaxationFactors
{

    equations
    {
        ".*"         1;
    }
}

[青山]

并且算着算着迭代次数逐渐变为0了

[李东岳]

T的进出口边界怎么给的

[青山]

@李东岳 进口固定值1，其他wall和出口都是零梯度

[李东岳]

你的扩散系数给的多少

[青山]

@李东岳 老师，因为我想固体颗粒在液体中应该没有扩散，所以我只考虑了对流，扩散系数设了0，是这个问题吗

[李东岳]

不是问题。把你的云图贴一个看看

[青山]

[李东岳]

你算的是合理的

[青山]

@李东岳 老师，我继续计算了4个小时也还是没有变化，出口一直没有浓度，而且T的迭代次数已经都是0了，就是说已经稳态了，可是这样如何得到停留时间分布？

[青山]

@李东岳 或者说是不是多相流中的单相速度场并不能直接拿到scalarTransportFoam中使用，因为单相速度场没有相含率信息，例如一个网格的固含率为0，但是其还是会有一个固相速度，然而scalarTransportFoam仅读取速度场无法知道这一点，会不会导致质量不守恒的问题。

[李东岳]

我继续计算了4个小时也还是没有变化，出口一直没有浓度，

你这个的原因，是因为流体的流场，并没有向出口走。你计算流场的时候，是存在进口速度的么？

[青山]

@李东岳 老师，问题解决了，要把固体的通量场也作为初值给scalarTransportFoam，并且要把phi里面得fixedValue类型改为calculated就可以了，谢谢老师

[李东岳]

scalarTransportFoam为什么会存在固体通量场呢，另外scalarTransportFoam不需要，也不能给phi的边界条件

[青山]

@李东岳 我的意思是把multiphaseEulerFoam算出来的U.particles和phi.particles一起复制过来作为scalarTransportFoam的初值，然后把原先phi.particles中的边界条件修改一下就可以了

[青山]

@李东岳 另外老师还想请教个问题，scalarTransportFoam中库朗数还需要小于1吗，还是说时间步长可以大一些让库朗数大一些

[李东岳]

不需要小于1，不过小于1会比较精准，另外大于1，要注意越界行为。如果没啥问题，可以大于1

[青山]

@李东岳 好的，谢谢老师