CFD中文网

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

@李东岳 谢谢李老师！

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 *** 求大神指点一二！

@henry-zhou 多给你一个选择，如果算les，直接改simulationType就可以

@henry-zhou Class method那面可以给定log粒径分布，平均直径满足文献值即可

在文献中看到一个颗粒相的粘度模型，随着颗粒相的体积分数变化。形式如下：7c540afb-dad6-43b0-a57b-38472090e828-image.png e87e0cf3-cd2a-4ced-a653-b008b44bcfd0-image.png 。请教下各位老师，这个在哪里设置？或者在哪里植入。现在完全没有思路:136: ，请各位老师不吝赐教！

第31页（页码13）最上面那段是不是你要的

https://www.jianguoyun.com/p/DYNwnmQQsLv2BRiS3fED

请问这个问题解决了吗？实验中颗粒会粘附在壁面上，首先第一层粘附颗粒肯定跟壁面之间有粘附力的作用。力的来源可能属于分子粘附力。请问这个力的方程是什么，有植入吗？
我最近在做液滴跟壁面附着的情况下受力如何变形以及多大的力使其滑动。。就是因为有附着力的存在，所以想植入，但是没经验。。附着力只作用在了接触面，算是表面力么。。感觉问题跟你的很像，不解决附着力的问题，就不太好算。。

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

@李东岳 东岳老师是想诱惑我再报一次您的辅导班:141:

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

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

tahts difficult to answer, I would suggest to get to know the algorithm first, if you want to know the settings, the knowledge should come to the first place. So, get some books to have a look at the algorithm, afterwards the settings shouldn't be a problem.