MPPIC中，尝试复杂形状网格计算一直不收敛

你的小祖宗

请教各位同学和老师，我在尝试用MPPIC在一个存在斜面的三维复杂网格中模拟固液两相流时，计算一直没办法收敛。尝试不同的边界条件和不同的k和epsilon值，都没有任何改善。以下是stack trace：

Lookup solver for epsilon.waterFinal
[16] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[16] #1  Foam::sigFpe::sigHandler(int) at ??:?
[16] #2  ? in "/lib64/libc.so.6"
[16] #3  Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[16] #4  Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:?
[16] #5  Foam::RASModels::kEpsilon<Foam::PhaseIncompressibleTurbulenceModel<Foam::singlePhaseTransportModel> >::correct() at ??:?
[16] #6  ? at ??:?
[16] #7  __libc_start_main in "/lib64/libc.so.6"
[16] #8  ? at ??:?
[cas019:410749:0:410749] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
==== backtrace ====
    0  /lustre/opt/cascadelake/linux-centos7-skylake_avx512/gcc-8.3.0/ucx-1.6.0-dacxhp7gepme2j7jmvpcmbm2qrptc6k7/lib/libucs.so.0(+0x1bf44) [0x2b968ef64f44]
    1  /lustre/opt/cascadelake/linux-centos7-skylake_avx512/gcc-8.3.0/ucx-1.6.0-dacxhp7gepme2j7jmvpcmbm2qrptc6k7/lib/libucs.so.0(+0x1c2ca) [0x2b968ef652ca]
    2  /lib64/libc.so.6(gsignal+0x37) [0x2b9689945337]
    3  /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0xd5) [0x2b9688a0dca5]
    4  /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4FoammlINS_7volMeshEEENS_3tmpINS_16DimensionedFieldIdT_EEEERKS6_RKS5_+0x164) [0x2b9683ba5d74]
    5  /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libDPMTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_34PhaseIncompressibleTurbulenceModelINS_25singlePhaseTransportModelEEEE7correctEv+0x547) [0x2b9687558727]
    6  MPPICFoam() [0x433d76]
    7  /lib64/libc.so.6(__libc_start_main+0xf5) [0x2b9689931505]
    8  MPPICFoam() [0x4352e0]

边界条件如下表

感觉应该是颗粒太大或者是网格划的不太好，有没有大佬能帮我看看怎么调整更准确一些