CFD中文网

@bestucan 万分感谢！感谢赵老师，感谢您！ 已经向您邮箱发了邮件zhou_hli@163.com

@bestucan 非常感谢您的指点！之前的确没有查到过赵海波的这本书。可是目前新书没卖的了，二手的好贵啊

@李东岳 谢谢李老师！

请教下各位老师，最近想详细学习下群体平衡模型，包括模型的构建过程，模型求解方法，聚并、破碎、成核、生长、溶解核函数模型等内容，有哪些中英文学习资料可以学习？恳请各位老师推荐，谢谢！

最近想详细学习下群体平衡模型，包括模型的构建，模型求解，聚并、破碎、成核、生长、溶解模型等内容，有哪些中英文学习资料可以看，麻烦各位老师推荐一下，谢谢啦！

multiphaseEulerfoam，低离散相分数的时候，没问题，离散相增加到0.35，就报错，尝试修改了边界条件、相间作用力模型、离散格式和网格都无果，差别就是跑的时间长短，最终都发散报错如下：

[8] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[8] #1  Foam::sigFpe::sigHandler(int) at ??:?
[8] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[8] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[8] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/zhl/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[8] #5  Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel>::nu() const at ??:?
[8] #6  Foam::RASModels::kEpsilon<Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> >::correct() at ??:?
[8] #7  Foam::phaseSystem::correctTurbulence() at ??:?
[8] #8  ? in "/home/zhl/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[8] #9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[8] #10  ?[ubuntu:09382] *** Process received signal ***
或者这样：
[8] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[8] #1  Foam::sigFpe::sigHandler(int) at ??:?
[8] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[8] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[8] #4  Foam::operator/(Foam::tmp<Foam::Field<double> > const&, Foam::UList<double> const&) at ??:?
[8] #5  Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel>::nu(int) const at ??:?
[8] #6  Foam::epsilonWallFunctionFvPatchScalarField::calculate(Foam::momentumTransportModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
[8] #7  Foam::epsilonWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::momentumTransportModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
[8] #8  Foam::epsilonWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
[8] #9  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:?
[8] #10  Foam::RASModels::kEpsilon<Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> >::correct() at ??:?
[8] #11  Foam::phaseSystem::correctTurbulence() at ??:?
[8] #12  ? in "/home/zhl/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
[8] #13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[8] #14  ?[ubuntu:07707] *** Process received signal ***
求大神指点一二！

请教下各位老师，在multiphaseEuler中湍流模型设置中，经产看到：

simulationType  RAS;
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
}
LES
{
    LESModel NicenoKEqn;
    turbulence      on;
    printCoeffs     on;
    delta           cubeRootVol;
    cubeRootVolCoeffs   
  {  
  }
}

或者：

simulationType  laminar;
RAS
{
    RASModel continuousGasKEpsilon;
    turbulence      off;
    printCoeffs     on;
}

类似的设置，不太明白已经是RAS或者 laminar了，下边的LES和RAS有何作用？
非常感谢老师们赐教。

@李东岳 文献中实验结果有平均直径数据，入口处给的是各相的体积分数， 大小范围可以参考结果，不知道初始粒径分布怎们给。另外液滴聚并用的CT模型，破碎用的LuoSvendsen，不知道合理不。非常感谢李老师的解答！

各位老师好！我在用multiphaseEulerFoam中pbm模型模拟液滴或者气泡聚并破碎，初始粒径的范围和分布，对平均直径的结果影响有多大？如果没有实验数据，该如何设置初始粒径的范围和分布呢？再有，一般要多少时间流动可以稳定？各位老师有没有好的经验。

在文献中看到一个颗粒相的粘度模型，随着颗粒相的体积分数变化。形式如下： 。请教下各位老师，这个在哪里设置？或者在哪里植入。现在完全没有思路 ，请各位老师不吝赐教！

东岳老师，我最近在做一个类似的模型，颗粒粘度与其体积分数有关的，用的of8的multiphaseEulerFoam，现在完全没有思路，不知道要将其放在哪里 望您指点一二

请教下各位大神，openfoam中是否可以加入自定义粘度，比如两相流动，某一相的粘度与其相分数相关的粘度模型，是否可以加入，应在哪里加入？ thermophysicalProperties中可以定义常数或者与温度相关的粘温模型，但是没有看到与相分数相关。

非常感谢！我之前把连续相的直径调小，效果好点，但是还是不理想。我一会儿试一下，把Sct调小，看看有变化没。

最近想做一个液固输送，颗粒较多，并且固体颗粒会粘附在壁面上生长。尝试了很多方法，但是都不太理想，固体颗粒总是跑掉，无法在壁面附近聚集，体积分数最多1e-5。
已经设置了壁面类型：

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U.particles;
}
dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (0 0 0);
    WALL
    {
        type            JohnsonJacksonParticleSlip;
        restitutionCoefficient 0.8;
        specularityCoefficient   0.00001;
        muF             0.7;
        sigma           2;
        value           uniform (0 0 0);
    }

颗粒间的作用设置如下：

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      turbulenceProperties.particles;
}

simulationType  RAS;

RAS
{
    model    kineticTheory;

    turbulence      on;
    printCoeffs     on;

    kineticTheoryCoeffs
    {
        equilibrium             off;
        continuousPhase         liquid;
        e                       0;
        alphaMax                0.63;
        alphaMinFriction        0.5;
        residualAlpha           1e-4;

        viscosityModel          Gidaspow;
        conductivityModel       Gidaspow;
        granularPressureModel   Lun;/
        frictionalStressModel   JohnsonJacksonSchaeffer;
        radialModel             SinclairJackson;

        JohnsonJacksonSchaefferCoeffs
        {
            Fr                      0.05;
            eta                     2;
            p                       5;
            phi                     28.5;
            alphaDeltaMin           0.05;
        }
    }

在thermophysicalProperties中的颗粒相粘度比较小1e-6量级。
同时在运行中，能看到总的体积分数Phase-sum volume fraction, min, max = 1.00007496 0.988884069 1.16288542
最大值超过1的情况，虽然没有发散，但是超过的有点多，很不合理，不知道什么原因。
先在这里感谢各位老师指教，谢谢！

大家好，我用multiphaseEulerFoam做一个案例，对某一相使用多组分模型multiComponentPhaseModel，但是相中组分分数分层。因为组分会扩散，正常情况是不应该分层的，调了一段时间，只是换了分层的位置。不知道大家遇到过这种情况不？谢谢各位老师的指点！

大家好，我在做一个算例，用multiphaseEuler模拟管道固液两相流动，其中固相会吸附在管壁上，形成沉积。最近一直在想如何实现，看了很多算例和边界条件都没有找到办法。所以想问下各位大神，openfoam是否有相应的速度壁面类型可以使用？

大家好，最近我在用双流体模型做一个沉积案例，其中一相由于重力会沉积到壁面，并且附着在上边，想请教一下各位大神，如何设置壁面类型或者某一相的参数，才能出现附着的效果。谢谢！！！

@浪迹天大 非常感谢！今天仔细看了下，终于搞定了。

@浪迹天大 十分感谢！但是如何获取我需要的分量？初学请指教一下 ，或者哪里可以查到？

@李东岳 是的李老师，of8中的PBE，加一个成核模型，就是Y-Csat_这一项出的问题，编译之后出现红色标记，出现上述错误。