fsiFoam 并行运算报错的问题

mpirun noticed that process rank 1 with PID 18430 on node localhost.localdomain exited on signal 11 (Segmentation fault).
GAMG:  Solving for p, Initial residual = 0.0058860803, Final residual = 7.5212105e-07, No Iterations 9
GAMG:  Solving for p, Initial residual = 0.00093355983, Final residual = 8.4947072e-07, No Iterations 6
time step continuity errors : sum local = 4.1825478e-11, global = -6.9814272e-12
GAMG:  Solving for p, Initial residual = 0.00029293937, Final residual = 4.0761422e-07, No Iterations 7
GAMG:  Solving for p, Initial residual = 5.566022e-05, Final residual = 3.5578998e-07, No Iterations 4
time step continuity errors : sum local = 1.7516981e-11, global = -3.2252254e-12
DILUPBiCG:  Solving for epsilon, Initial residual = 6.9644483e-05, Final residual = 1.7462709e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00019727494, Final residual = 6.6749395e-07, No Iterations 1
Setting traction on solid patch
Total force (fluid) = (-0.36750964 -0.0078857045 -2.285684e-07)
Total force (solid) = (0.36733414 0.0079405659 2.2894511e-07)
Evolving stress model
[[email protected] fluid]$ --------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 41651 on node localhost.localdomain exited on signal 11 (Segmentation fault).
--------------------------------------------------------------------------

每次计算完十个步长以后就报错,它这样我一个约化速度都算不出

Time = 0.0025, iteration: 3
Maximal accumulated displacement of interface points: 0.00032631935
GAMG:  Solving for cellMotionUx, Initial residual = 0.4925908, Final residual = 0.00018380911, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0.72466457, Final residual = 0.00070954375, No Iterations 3
GAMG:  Solving for cellMotionUz, Initial residual = 0.49491501, Final residual = 0.00031470224, No Iterations 3
GAMG:  Solving for cellMotionUx, Initial residual = 0.00021824665, Final residual = 7.3600322e-07, No Iterations 4
GAMG:  Solving for cellMotionUy, Initial residual = 0.00047565705, Final residual = 4.9458908e-07, No Iterations 4
GAMG:  Solving for cellMotionUz, Initial residual = 0.00028304055, Final residual = 7.082873e-07, No Iterations 4
Evolving flow model
Courant Number mean: 0.0012308981 max: 0.024856613 velocity magnitude: 0.06863784
DILUPBiCG:  Solving for Ux, Initial residual = 1.2973902e-05, Final residual = 3.2127674e-09, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1.5359233e-05, Final residual = 5.161974e-09, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 9.3678323e-05, Final residual = 2.6852861e-08, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.030644822, Final residual = 8.9685037e-07, No Iterations 13
GAMG:  Solving for p, Initial residual = 0.018629844, Final residual = 3.9145968e-07, No Iterations 10
time step continuity errors : sum local = 1.69041e-12, global = 5.6256055e-13
GAMG:  Solving for p, Initial residual = 0.0038889137, Final residual = 9.9233935e-07, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.00063893118, Final residual = 6.9922725e-07, No Iterations 6
time step continuity errors : sum local = 3.0060493e-12, global = 6.5550613e-14
GAMG:  Solving for p, Initial residual = 0.00012569669, Final residual = 8.5840571e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 2.3154034e-05, Final residual = 4.1048321e-07, No Iterations 3
time step continuity errors : sum local = 1.7647374e-12, global = -2.5456673e-13
DILUPBiCG:  Solving for epsilon, Initial residual = 6.0710465e-07, Final residual = 6.0710465e-07, No Iterations 0
DILUPBiCG:  Solving for k, Initial residual = 0.00026085173, Final residual = 5.3623471e-07, No Iterations 1
Setting traction on solid patch
Total force (fluid) = (-0.040457893 -8.8885722e-05 2.2119027e-09)
Total force (solid) = (0.040437818 8.9133019e-05 -2.2073432e-09)
Evolving stress model
DICPCG: Solving for D, Initial residula = 0.00056605288, Final residual = 7.280206e-05, No outer iterations = 200
Max relative residual = 0.0022285338, Relative residual = 3.592933e-05, enforceLinear = 0
Current fsi relative residual norm: 0.13985263
DICPCG: Solving for D, Initial residula = 0.00056605288, Final residual = 7.280206e-05, No outer iterations = 200
Max relative residual = 0.0022285338, Relative residual = 3.592933e-05, enforceLinear = 0
Current fsi relative residual norm: 0.13985263

这是log文件中最后的内容,不光是并行运算报错,直接fsi也会报错,都在相同的时刻

不知道并行的网格是怎么分的?

z方向分成4块,不使用并行也会报错,问题可能不在并行上。报错内容为

segmentation fault (core dumped) 

这个问题比较紧急,不然没法计算下去

第一循环算完没问题,出在保存上,可能跟数据交换有关系。把 VTK 文件夹删除,重做 createZones ,再进行运算,就没问题了。

并行计算暂时还不能。

大概什么时候能并行运算呢

我看了下程序,程序本身是可以的。并行网格的划分上可能有些技巧。

并行的网格划分可以试一下

decomposeParFsi

我试了下decomposeparfsi,但是并不行