用twoPhaseEulerFoam求解水砂输送的问题

[李东岳]

液相、固相注入速度都是（0.5 0 0），但计算过程中底部达到4m/s+，

这是用层流模拟的么？好像速度确实有点大。不过也不好判断。没准是重力的作用加速了？

[Zhy2022]

@东岳 RAS，但是重力也不该水平方向这么大的速度哇，主要是没堆积

[Zhy2022]

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

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

[李东岳]

把你fvScheme发一下

[Zhy2022]

@东岳 老师，您看：

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

[Zhy2022]

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

[Zhy2022]

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

[Zhy2022]

@李东岳 东岳老师，在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哇？

[Zhy2022]

同样工况下，fluent模拟结果和OF模拟结果存在差异：
fluent的砂堤快速堆积，而OF的砂堤会先跑满整个空间，再堆积回缩，堆积也达不到最大堆积程度0.63，
而且OF中固相Y方向速度有局部异常向下的特点，不知道是不是固相计算有什么问题，不知东岳老师@李东岳 或者其他前辈能否给一些建议？

[李东岳]

你这个是否有参考的sci或者参考数据可以对比？有的话可以把openfoam算例清理干净上传给我我给你调一下，如果不涉密的话

[Zhy2022]

@李东岳 谢谢东岳老师！暂无可参考的sci，目前仅能通过实验结果进行粗略对比。
最大问题还是固相分数回流、沉降达不到最大堆积分数0.62和固相Y速度不太对劲。
case如何上传给您？实验视频有10个G不太方便发送，我以截图的方式一并发给您

[Zhy2022]

@李东岳 东岳老师，有一篇Mfix的模拟结果，出口设置和我的模型略有不同，影响不大，一并发给您

[veen]

@Zhy2022
之前在openfoam里边算过空气中沙物质沉降，算出来最大体积堆积分数和fluent是可以对上的，会不会是alphaMax 这个参数在设置的问题呢？

[Zhy2022]

@veen 谢谢前辈，我打的最大相分数给的0.62，但是堆积只有0.59左右，这也是个很奇怪的问题

[Zhy2022]

@veen 前辈，您的设置能否分享一下（不用网格）？

[李东岳]

@Zhy2022 不好意思算例太大了，你把算例清理好只保留必要文件发给我就好，不需要结果文件

[veen]

@Zhy2022 刚刚看了我的alphamax也没改，链接在这：https://github.com/Veenxz/Veen_OpenFOAM_Cases/tree/main/2D_Plane_EE_Multiphase

[veen]

@Zhy2022 你直接blockmesh改一下，然后改一下流体参数应该就好了，曳力函数也需要改一下。

[Zhy2022]

@李东岳 好的东岳老师，其实里面没有结果文件，只是我实验做的一点结果和文章 我再给您发一个只有算例的

[Zhy2022]

@veen 谢谢前辈！ 目前打算先用gidaspow算算 后面会修正曳力模型