分散相体积分数增加，报错！
-
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 *** 求大神指点一二！