solids4foam运行算例，流固耦合交界面报错

[洱聿]

在solids4foam中运行验证算例（孤立波下斜坡上的直墙变形和受力），但是当波浪到固体交界面时，就会出现浮点报错问题，修改了网格和固体交界面的问题，还是没有改善这个问题。请问有大佬们知道怎么解决这个问题吗？
报错时间步alpha场：（每次到这个地方就出现报错）

报错部分log文件：

Time = 8.933, iteration: 9
Current fsi under-relaxation factor (Aitken) of interface: 3.24861e-06
Maximal accumulated displacement of interface 0: 2.02709e-06
GAMG:  Solving for cellMotionUx, Initial residual = 0.27053, Final residual = 4.88471e-05, No Iterations 2
GAMG:  Solving for cellMotionUz, Initial residual = 0.276424, Final residual = 7.60589e-05, No Iterations 2
GAMG:  Solving for cellMotionUx, Initial residual = 5.19396e-05, Final residual = 4.62925e-07, No Iterations 4
GAMG:  Solving for cellMotionUz, Initial residual = 8.68504e-05, Final residual = 8.18649e-07, No Iterations 3
Evolving fluid model: wave2solids
PIMPLE: iteration 1
MULES: Solving for alpha1
alpha1 BC on patch inlet
alpha1 BC on patch inlet
Liquid phase volume fraction = 0.259489  Min(alpha1) = -2191.92  Max(alpha1) = 2360.17
MULES: Solving for alpha1
alpha1 BC on patch inlet
alpha1 BC on patch inlet
Liquid phase volume fraction = 0.259489  Min(alpha1) = -158062  Max(alpha1) = 163488
Velocity BC on patch inlet
3D_2D Absorption BC on patch outlet
"Correction Levels" 1( 0 )

DICPCG:  Solving for pd, Initial residual = 0.10778, Final residual = 0.00369688, No Iterations 8
DICPCG:  Solving for pd, Initial residual = 6.32049e-06, Final residual = 2.94799e-07, No Iterations 644
DICPCG:  Solving for pd, Initial residual = 4.65247e-05, Final residual = 2.1618e-06, No Iterations 637
DICPCG:  Solving for pd, Initial residual = 2.0431e-07, Final residual = 9.74494e-08, No Iterations 110
Velocity BC on patch inlet
3D_2D Absorption BC on patch outlet
"Correction Levels" 1( 0 )

DICPCG:  Solving for pd, Initial residual = 4.64318e-05, Final residual = 2.21746e-06, No Iterations 638
DICPCG:  Solving for pd, Initial residual = 2.16117e-07, Final residual = 8.27667e-08, No Iterations 132
Velocity BC on patch inlet
3D_2D Absorption BC on patch outlet
"Correction Levels" 1( 0 )

Setting traction on solid interfaces
Mapping face values using directMap
Total force on fluid interface 0: (4.3321e+10 -2.82466e-19 -0.673933)
Total force on solid interface 0: (-4.3321e+10 2.82466e-19 0.673933)

Evolving solid solver
Solving the updated Lagrangian form of the momentum equation for DD
--------------------------------------------------------------------------
mpirun noticed that process rank 6 with PID 168600 on node node03 exited on signal 8 (Floating point exception).