@李东岳 东岳老师，在kineticTheoryModel.C中有这样一段代码：

Foam::tmp<Foam::volScalarField>
Foam::RASModels::kineticTheoryModel::k() const
{
    NotImplemented;
    return nut_;
}


Foam::tmp<Foam::volScalarField>
Foam::RASModels::kineticTheoryModel::epsilon() const
{
    NotImplemented;
    return nut_;
}

但是感觉固相湍流开关打开的话，其湍流黏度完全是按照如下公式计算的：
（其中a为固相，b为液相）
也就是固相被看作连续相后不存在真实的k和epsilon哇？

@东岳 东岳老师，有个问题不太想得通，不知道您对此看法如何：
采用TFM模拟固-液两相流，固相密度设置为常数（fluidisedBed案例中也是如此），视为不可压缩，但固相湍流模型（kineticTheory或phasePressure）都是属于可压缩模型，它们是否矛盾？

@东岳 东岳老师，问题已解决！无法编译与文件路径有关

@东岳 东岳老师，没能找出报错原因 ，只能再来麻烦一下您：

INLET1
    {
        type            codedFixedValue;
        value           $internalField;
        name            stepInjection;
        redictType      variedValue;
        code
        #{
             scalar alphapar = this->patch().fluid().phase();

             scalar t = this->db().time().value();
       
             if(t >= 0 && t<=5)
             {
                  alphapar = 0;
             }
             
             else
	     {
	          alphapar = 0.15;
	     }

	     operator==(alphapar)

         #};
        
    }

@东岳 老师，问题解决了：
固相的速度边界和湍流模型必须配套使用才行！

@东岳 老师，您看：

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;

    div((((alpha.particles*thermo:rho.particles)*nuEff.particles)*dev2(T(grad(U.particles))))) Gauss linear;//改层流后多出
}

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

interpolationSchemes
{
    default     linear;
}

@东岳 东岳老师您好，还想请问一下：我按照of自带算例将固相模型改为从kineticTheory改为phasePressure，以实现层流模拟，报错：

我检查了初始场文件发现没有问题，不知这是怎么回事。。还请您看一看这个gs0.particles是什么东西。。

东岳老师您好，请问如果是两相流中相分数（比如alpha.particle）随时间变化，它在codeFixedValue中应该怎么定义呢？

@hoversoar 嗯 谢谢

@wwzhao 您好，我就是遇到了这种情况。。因为平均库朗数太小了，只要不影响收敛，最大库朗数调到30应该是没有问题的吧？