请教各位同学和老师，我在尝试用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]
边界条件如下表
感觉应该是颗粒太大或者是网格划的不太好，有没有大佬能帮我看看怎么调整更准确一些