只修改了粘度出现一下报错想问问该如何解决，粘度为2e-6m2/s
[8] #0 Foam::error::printStack(Foam::Ostream&)[11] #0 [14] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) at ??:?
[11] #1 Foam::sigFpe::sigHandler(int) at ??:?
[8] #1 Foam::sigFpe::sigHandler(int) at ??:?
[14] #1 Foam::sigFpe::sigHandler(int) at ??:?
[11] #2 ? at ??:?
[8] #2 ? at ??:?
[14] #2 ? in "/lib64/libpthread.so.0"
[8] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libpthread.so.0"
[11] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libpthread.so.0"
[14] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
at ??:?
[8] #4 void Foam::multiply<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&)[11] #4 void Foam::multiply<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[14] #4 void Foam::multiply<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[8] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[11] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[14] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[8] #6 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() at ??:?
[11] #6 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() at ??:?
[14] #6 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() at ??:?
[8] #7 at ??:?
[11] #7 at ??:?
[14] #7 ??? at ??:?
[8] #8 __libc_start_main at ??:?
[11] #8 __libc_start_main at ??:?
[14] #8 __libc_start_main in "/lib64/libc.so.6"
[8] #9 in "/lib64/libc.so.6"
[11] #9 in "/lib64/libc.so.6"
[14] #9 ??? at ??:?
[ibnode77:65243:0:65243] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
at ??:?
[ibnode77:65246:0:65246] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
at ??:?
[ibnode77:65249:0:65249] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
==== backtrace (tid: 65243) ====
0 0x000000000000f4fb raise() ???:0
1 0x000000000069a44d Foam::multiply() ???:0
2 0x00000000000d4418 Foam::multiply<Foam::fvPatchField, Foam::volMesh>() ???:0
3 0x000000000015f8b6 Foam::operator*<Foam::fvPatchField, Foam::volMesh>() ???:0
4 0x00000000001ab676 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() ???:0
5 0x000000000042a05e main() ???:0
6 0x0000000000022555 __libc_start_main() ???:0
7 0x000000000042e172 _start() ???:0

[ibnode77:65243] *** Process received signal ***
[ibnode77:65243] Signal: Floating point exception (8)
[ibnode77:65243] Signal code: (-6)
[ibnode77:65243] Failing at address: 0x4080000fedb
[ibnode77:65243] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b424c49c630]
[ibnode77:65243] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b424c49c4fb]
[ibnode77:65243] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b424c49c630]
[ibnode77:65243] [ 3] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6+0x11d)[0x2b424b51044d]
[ibnode77:65243] [ 4] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4Foam8multiplyINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9+0xa8)[0x2b4247910418]
[ibnode77:65243] [ 5] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4FoammlINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_RKS8+0x156)[0x2b424799b8b6]
[ibnode77:65243] [ 6] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels11RNGkEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x386)[0x2b42479e7676]
[ibnode77:65243] [ 7] cfdemSolverPiso[0x42a05e]
[ibnode77:65243] [ 8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b424c6cb555]
[ibnode77:65243] [ 9] cfdemSolverPiso[0x42e172]
[ibnode77:65243] *** End of error message ***
==== backtrace (tid: 65246) ====
0 0x000000000000f4fb raise() ???:0
1 0x000000000069a44d Foam::multiply() ???:0
2 0x00000000000d4418 Foam::multiply<Foam::fvPatchField, Foam::volMesh>() ???:0
3 0x000000000015f8b6 Foam::operator*<Foam::fvPatchField, Foam::volMesh>() ???:0
4 0x00000000001ab676 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() ???:0
5 0x000000000042a05e main() ???:0
6 0x0000000000022555 __libc_start_main() ???:0
7 0x000000000042e172 _start() ???:0

[ibnode77:65246] *** Process received signal ***
[ibnode77:65246] Signal: Floating point exception (8)
[ibnode77:65246] Signal code: (-6)
[ibnode77:65246] Failing at address: 0x4080000fede
[ibnode77:65246] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b3d893ab630]
[ibnode77:65246] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b3d893ab4fb]
[ibnode77:65246] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b3d893ab630]
[ibnode77:65246] [ 3] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6+0x11d)[0x2b3d8841f44d]
[ibnode77:65246] [ 4] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4Foam8multiplyINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9+0xa8)[0x2b3d8481f418]
[ibnode77:65246] [ 5] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4FoammlINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_RKS8+0x156)[0x2b3d848aa8b6]
[ibnode77:65246] [ 6] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels11RNGkEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x386)[0x2b3d848f6676]
[ibnode77:65246] [ 7] cfdemSolverPiso[0x42a05e]
[ibnode77:65246] [ 8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b3d895da555]
[ibnode77:65246] [ 9] cfdemSolverPiso[0x42e172]
[ibnode77:65246] *** End of error message ***
==== backtrace (tid: 65249) ====
0 0x000000000000f4fb raise() ???:0
1 0x000000000069a44d Foam::multiply() ???:0
2 0x00000000000d4418 Foam::multiply<Foam::fvPatchField, Foam::volMesh>() ???:0
3 0x000000000015f8b6 Foam::operator*<Foam::fvPatchField, Foam::volMesh>() ???:0
4 0x00000000001ab676 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() ???:0
5 0x000000000042a05e main() ???:0
6 0x0000000000022555 __libc_start_main() ???:0
7 0x000000000042e172 _start() ???:0

[ibnode77:65249] *** Process received signal ***
[ibnode77:65249] Signal: Floating point exception (8)
[ibnode77:65249] Signal code: (-6)
[ibnode77:65249] Failing at address: 0x4080000fee1
[ibnode77:65249] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b14b9db5630]
[ibnode77:65249] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b14b9db54fb]
[ibnode77:65249] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b14b9db5630]
[ibnode77:65249] [ 3] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6+0x11d)[0x2b14b8e2944d]
[ibnode77:65249] [ 4] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4Foam8multiplyINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9+0xa8)[0x2b14b5229418]
[ibnode77:65249] [ 5] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4FoammlINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_RKS8+0x156)[0x2b14b52b48b6]
[ibnode77:65249] [ 6] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels11RNGkEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x386)[0x2b14b5300676]
[ibnode77:65249] [ 7] cfdemSolverPiso[0x42a05e]
[ibnode77:65249] [ 8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b14b9fe4555]
[ibnode77:65249] [ 9] cfdemSolverPiso[0x42e172]
[ibnode77:65249] *** End of error message ***

Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.

mpirun noticed that process rank 14 with PID 0 on node ibnode77 exited on signal 8 (Floating point exception).

@小刘lyw 哈哈哈我主要还是对颗粒数有疑惑，实在不明白应该具体算多少，有好几种公式，可以探讨一下，

@李东岳 文章中没有具体的流场结果，都是流场和粒子耦合在一起的结果，但是改变流体速度值域是可以对的上的

你好我也是做相关方面的方便交流一下吗

各位好，关于复现文章出现了一些问题，和文章设置一样的参数，同时开启了DiFeliceDrag、gradPForce、viscForce、 virtualMassForce，但结果是粒子堆积现象严重，没有悬浮起来，（感觉拖曳力的影响较小）想问问可能造成的原因有哪些？
麻烦各位帮忙看看！非常感谢