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

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