FSFI solver 流固流耦合报错

xye111 · 2019 年6 月 11 日 13:32


Time = 0.00025, iteration: 19
Maximal accumulated displacement of interface points: 0.0088426663
Total force (fluid2) = (-27.346619 0.0053126885 -0.00011606031)
Total force (solid) = (27.359856 -0.0060967625 0.00011609978)
Evolving stress model
GAMG:  Solving for cellMotionUx, Initial residual = 0.38829468, Final residual = 5.5106967e-06, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0.30829545, Final residual = 0.00013109314, No Iterations 1
GAMG:  Solving for cellMotionUz, Initial residual = 0.31272232, Final residual = 5.6901553e-05, No Iterations 1
GAMG:  Solving for cellMotionUx, Initial residual = 8.6729416e-06, Final residual = 9.9984975e-07, No Iterations 2
GAMG:  Solving for cellMotionUy, Initial residual = 0.00018559714, Final residual = 3.5385504e-07, No Iterations 4
GAMG:  Solving for cellMotionUz, Initial residual = 8.2761148e-05, Final residual = 8.3616596e-07, No Iterations 3
Evolving flow model: consistentIcoFlow
Evolving flow model
Courant Number mean: 0.0011329055 max: 0.013267262 velocity magnitude: 0.12300909
DILUPBiCG:  Solving for Ux, Initial residual = 8.6472069e-05, Final residual = 3.1415702e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00087550742, Final residual = 3.3169054e-09, No Iterations 1
GAMG:  Solving for cellMotionUx, Initial residual = 0.88361913, Final residual = 0.00034671998, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00043769359, Final residual = 2.0235658e-09, No Iterations 1
Consistent backwardDdtPhiCorr
Courant Number mean: 0.12846723 max: 2.1358503 velocity magnitude: 5.8791266
GAMG:  Solving for cellMotionUy, Initial residual = 0.9855137, Final residual = 0.00024652486, No Iterations 4
GAMG:  Solving for cellMotionUz, Initial residual = 0.99597287, Final residual = 0.00092237877, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.004098423, Final residual = 9.6863841e-07, No Iterations 68
GAMG:  Solving for cellMotionUx, Initial residual = 0.88361913, Final residual = 0.00034671998, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0017857355, Final residual = 8.2634475e-07, No Iterations 8
time step continuity errors : sum local = 3.4824554e-08, global = -5.2581913e-09
Consistent backwardDdtPhiCorr
GAMG:  Solving for p, Initial residual = 0.00021537424, Final residual = 8.0490705e-07, No Iterations 4
GAMG:  Solving for p, Initial residual = 2.1936024e-05, Final residual = 8.4014851e-07, No Iterations 3
time step continuity errors : sum local = 3.5346556e-08, global = 7.5182808e-09
Consistent backwardDdtPhiCorr
GAMG:  Solving for p, Initial residual = 3.3760269e-06, Final residual = 6.9865106e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 9.3289046e-07, Final residual = 7.4048799e-07, No Iterations 1
time step continuity errors : sum local = 3.1153171e-08, global = 7.5586552e-09
Setting traction on solid patch
Total force (fluid2) = (-27.346619 0.0053126885 -0.00011606031)
Total force (solid) = (27.359856 -0.0060967625 0.00011609978)
Evolving stress model
GAMG:  Solving for cellMotionUy, Initial residual = 0.9855137, Final residual = 0.00024652486, No Iterations 4
--------------------------------------------------------------------------
mpirun noticed that process rank 4 with PID 21678 on node localhost.localdomain exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

使用了8个线程进行了并行运算，第一个时间步长就算不过去，最后就浮点例外

xye111 · 2019 年6 月 12 日 07:27

流固流耦合报错以后，FSI并行运算也发生了同样的错误。情况紧急，如下所示：

Time = 1e-05, iteration: 5
Maximal accumulated displacement of interface points: 0.046676295
GAMG:  Solving for p, Initial residual = 0.0032315607, Final residual = 8.9507388e-07, No Iterations 33
GAMG:  Solving for cellMotionUx, Initial residual = 0.00010647769, Final residual = 5.9148222e-07, No Iterations 3
time step continuity errors : sum local = 1.2845725e-08, global = -2.9601403e-09
GAMG:  Solving for cellMotionUy, Initial residual = 9.9913063e-05, Final residual = 5.5415777e-07, No Iterations 3
GAMG:  Solving for cellMotionUz, Initial residual = 8.2015774e-05, Final residual = 3.9931276e-07, No Iterations 3
Evolving flow model
Courant Number mean: 0.0096259259 max: 0.27094059 velocity magnitude: 16.789765
DILUPBiCG:  Solving for Ux, Initial residual = 0.023487621, Final residual = 3.0854568e-08, No Iterations 3
GAMG:  Solving for cellMotionUx, Initial residual = 0.52537202, Final residual = 9.7056478e-05, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0.51824086, Final residual = 9.1221982e-05, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.026018225, Final residual = 3.4290884e-08, No Iterations 3
GAMG:  Solving for cellMotionUz, Initial residual = 0.46546013, Final residual = 8.4091306e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0049264123, Final residual = 9.1385776e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0023274793, Final residual = 9.4479932e-07, No Iterations 18
DILUPBiCG:  Solving for Ux, Initial residual = 0.023487621, Final residual = 3.0854568e-08, No Iterations 3
Evolving flow model
DILUPBiCG:  Solving for Uy, Initial residual = 0.026018225, Final residual = 3.4290884e-08, No Iterations 3
GAMG:  Solving for cellMotionUx, Initial residual = 0.00010647769, Final residual = 5.9148222e-07, No Iterations 3
Courant Number mean: 0.0096259259 max: 0.27094059 velocity magnitude: 16.789765
GAMG:  Solving for cellMotionUy, Initial residual = 9.9913063e-05, Final residual = 5.5415777e-07, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.0049264123, Final residual = 9.1385776e-07, No Iterations 2
GAMG:  Solving for cellMotionUz, Initial residual = 8.2015774e-05, Final residual = 3.9931276e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00047997957, Final residual = 6.3947772e-07, No Iterations 8
time step continuity errors : sum local = 9.1772512e-09, global = -1.5690527e-09
DILUPBiCG:  Solving for Ux, Initial residual = 0.023487621, Final residual = 3.0854568e-08, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.026018225, Final residual = 3.4290884e-08, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.0049264123, Final residual = 9.1385776e-07, No Iterations 2
Evolving flow model
GAMG:  Solving for p, Initial residual = 0.00036105454, Final residual = 8.287778e-07, No Iterations 7
Courant Number mean: 0.0096259259 max: 0.27094059 velocity magnitude: 16.789765
GAMG:  Solving for p, Initial residual = 0.00010991433, Final residual = 4.8045169e-07, No Iterations 6
time step continuity errors : sum local = 6.8954827e-09, global = 1.1207701e-09
DILUPBiCG:  Solving for epsilon, Initial residual = 0.20830993, Final residual = 3.7015182e-07, No Iterations 3
DILUPBiCG:  Solving for Ux, Initial residual = 0.023487621, Final residual = 3.0854568e-08, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.026018225, Final residual = 3.4290884e-08, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.27867101, Final residual = 1.6392718e-08, No Iterations 4
bounding k, min: -0.02642606 max: 0.85658885 average: 0.038971731
Setting traction on solid patch
DILUPBiCG:  Solving for Uz, Initial residual = 0.0049264123, Final residual = 9.1385776e-07, No Iterations 2
Total force (fluid) = (-0.012684776 0.00097461781 -1.2742929)
Total force (solid) = (0.021626368 0.0031878273 1.273759)
Evolving stress model
GAMG:  Solving for p, Initial residual = 0.021515055, Final residual = 9.8993892e-07, No Iterations 85
GAMG:  Solving for p, Initial residual = 0.021515055, Final residual = 9.8993892e-07, No Iterations 85
GAMG:  Solving for p, Initial residual = 0.021515055, Final residual = 9.8993892e-07, No Iterations 85
GAMG:  Solving for p, Initial residual = 0.0032315607, Final residual = 8.9507388e-07, No Iterations 33
time step continuity errors : sum local = 1.2845725e-08, global = -2.9601403e-09
GAMG:  Solving for p, Initial residual = 0.021515055, Final residual = 9.8993892e-07, No Iterations 85
GAMG:  Solving for p, Initial residual = 0.0032315607, Final residual = 8.9507388e-07, No Iterations 33
time step continuity errors : sum local = 1.2845725e-08, global = -2.9601403e-09
GAMG:  Solving for p, Initial residual = 0.0032315607, Final residual = 8.9507388e-07, No Iterations 33
time step continuity errors : sum local = 1.2845725e-08, global = -2.9601403e-09
GAMG:  Solving for p, Initial residual = 0.0023274793, Final residual = 9.4479932e-07, No Iterations 18
GAMG:  Solving for p, Initial residual = 0.0023274793, Final residual = 9.4479932e-07, No Iterations 18
GAMG:  Solving for p, Initial residual = 0.0032315607, Final residual = 8.9507388e-07, No Iterations 33
time step continuity errors : sum local = 1.2845725e-08, global = -2.9601403e-09
GAMG:  Solving for p, Initial residual = 0.021515055, Final residual = 9.8993892e-07, No Iterations 85
GAMG:  Solving for p, Initial residual = 0.00047997957, Final residual = 6.3947772e-07, No Iterations 8
time step continuity errors : sum local = 9.1772512e-09, global = -1.5690527e-09
GAMG:  Solving for p, Initial residual = 0.00047997957, Final residual = 6.3947772e-07, No Iterations 8
time step continuity errors : sum local = 9.1772512e-09, global = -1.5690527e-09
GAMG:  Solving for p, Initial residual = 0.00036105454, Final residual = 8.287778e-07, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.0023274793, Final residual = 9.4479932e-07, No Iterations 18
GAMG:  Solving for p, Initial residual = 0.00036105454, Final residual = 8.287778e-07, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.00010991433, Final residual = 4.8045169e-07, No Iterations 6
time step continuity errors : sum local = 6.8954827e-09, global = 1.1207701e-09
GAMG:  Solving for p, Initial residual = 0.00047997957, Final residual = 6.3947772e-07, No Iterations 8
time step continuity errors : sum local = 9.1772512e-09, global = -1.5690527e-09
GAMG:  Solving for p, Initial residual = 0.0023274793, Final residual = 9.4479932e-07, No Iterations 18
GAMG:  Solving for p, Initial residual = 0.00010991433, Final residual = 4.8045169e-07, No Iterations 6
time step continuity errors : sum local = 6.8954827e-09, global = 1.1207701e-09
DILUPBiCG:  Solving for epsilon, Initial residual = 0.20830993, Final residual = 3.7015182e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0032315607, Final residual = 8.9507388e-07, No Iterations 33
time step continuity errors : sum local = 1.2845725e-08, global = -2.9601403e-09
GAMG:  Solving for p, Initial residual = 0.00047997957, Final residual = 6.3947772e-07, No Iterations 8
time step continuity errors : sum local = 9.1772512e-09, global = -1.5690527e-09
DILUPBiCG:  Solving for k, Initial residual = 0.27867101, Final residual = 1.6392718e-08, No Iterations 4
bounding k, min: -0.02642606 max: 0.85658885 average: 0.038971731
Setting traction on solid patch
DILUPBiCG:  Solving for epsilon, Initial residual = 0.20830993, Final residual = 3.7015182e-07, No Iterations 3
Total force (fluid) = (-0.012684776 0.00097461781 -1.2742929)
Total force (solid) = (0.021626368 0.0031878273 1.273759)
Evolving stress model
GAMG:  Solving for p, Initial residual = 0.00036105454, Final residual = 8.287778e-07, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.00010991433, Final residual = 4.8045169e-07, No Iterations 6
time step continuity errors : sum local = 6.8954827e-09, global = 1.1207701e-09
DILUPBiCG:  Solving for k, Initial residual = 0.27867101, Final residual = 1.6392718e-08, No Iterations 4
bounding k, min: -0.02642606 max: 0.85658885 average: 0.038971731
GAMG:  Solving for p, Initial residual = 0.00036105454, Final residual = 8.287778e-07, No Iterations 7
Setting traction on solid patch
Total force (fluid) = (-0.012684776 0.00097461781 -1.2742929)
Total force (solid) = (0.021626368 0.0031878273 1.273759)
Evolving stress model
GAMG:  Solving for p, Initial residual = 0.0023274793, Final residual = 9.4479932e-07, No Iterations 18
DILUPBiCG:  Solving for epsilon, Initial residual = 0.20830993, Final residual = 3.7015182e-07, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.27867101, Final residual = 1.6392718e-08, No Iterations 4
bounding k, min: -0.02642606 max: 0.85658885 average: 0.038971731
GAMG:  Solving for p, Initial residual = 0.00010991433, Final residual = 4.8045169e-07, No Iterations 6
time step continuity errors : sum local = 6.8954827e-09, global = 1.1207701e-09
Setting traction on solid patch
Total force (fluid) = (-0.012684776 0.00097461781 -1.2742929)
Total force (solid) = (0.021626368 0.0031878273 1.273759)
Evolving stress model
GAMG:  Solving for p, Initial residual = 0.00047997957, Final residual = 6.3947772e-07, No Iterations 8
time step continuity errors : sum local = 9.1772512e-09, global = -1.5690527e-09
DILUPBiCG:  Solving for epsilon, Initial residual = 0.20830993, Final residual = 3.7015182e-07, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.27867101, Final residual = 1.6392718e-08, No Iterations 4
bounding k, min: -0.02642606 max: 0.85658885 average: 0.038971731
Setting traction on solid patch
Total force (fluid) = (-0.012684776 0.00097461781 -1.2742929)
Total force (solid) = (0.021626368 0.0031878273 1.273759)
Evolving stress model
GAMG:  Solving for p, Initial residual = 0.00036105454, Final residual = 8.287778e-07, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.00010991433, Final residual = 4.8045169e-07, No Iterations 6
time step continuity errors : sum local = 6.8954827e-09, global = 1.1207701e-09
DILUPBiCG:  Solving for epsilon, Initial residual = 0.20830993, Final residual = 3.7015182e-07, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.27867101, Final residual = 1.6392718e-08, No Iterations 4
bounding k, min: -0.02642606 max: 0.85658885 average: 0.038971731
Setting traction on solid patch
Total force (fluid) = (-0.012684776 0.00097461781 -1.2742929)
Total force (solid) = (0.021626368 0.0031878273 1.273759)
Evolving stress model
^C
[xye@localhost fluid]$ --------------------------------------------------------------------------
mpirun noticed that process rank 5 with PID 6065 on node localhost.localdomain exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

xye111 · 2019 年6 月 12 日 15:11

Time = 0.00025, iteration: 9
Maximal accumulated displacement of interface points: 0.0034953781
GAMG:  Solving for cellMotionUx, Initial residual = 0.73894454, Final residual = 3.7578043e-06, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0.94227586, Final residual = 0.00041286968, No Iterations 1
GAMG:  Solving for cellMotionUz, Initial residual = 0.78251055, Final residual = 0.00014535672, No Iterations 1
GAMG:  Solving for cellMotionUx, Initial residual = 2.6103666e-05, Final residual = 9.3200422e-07, No Iterations 3
GAMG:  Solving for cellMotionUy, Initial residual = 0.0024125592, Final residual = 1.2830121e-06, No Iterations 5
GAMG:  Solving for cellMotionUz, Initial residual = 0.00031814594, Final residual = 9.5479077e-07, No Iterations 4
Evolving flow model: consistentIcoFlow
Courant Number mean: 0.00061364502 max: 0.0045949843 velocity magnitude: 0.045928833
DILUPBiCG:  Solving for Ux, Initial residual = 0.00019214644, Final residual = 3.3203918e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.001843396, Final residual = 3.8025855e-09, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00090397517, Final residual = 1.7356225e-09, No Iterations 1
Consistent backwardDdtPhiCorr
GAMG:  Solving for p, Initial residual = 0.0067869395, Final residual = 9.9594813e-07, No Iterations 88
GAMG:  Solving for p, Initial residual = 0.0029821264, Final residual = 9.511378e-07, No Iterations 10
time step continuity errors : sum local = 1.3800486e-08, global = -5.6654649e-10
Consistent backwardDdtPhiCorr
GAMG:  Solving for p, Initial residual = 0.00035889845, Final residual = 7.113431e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 3.6504751e-05, Final residual = 7.5432329e-07, No Iterations 3
time step continuity errors : sum local = 1.0912172e-08, global = -6.9188946e-10
Consistent backwardDdtPhiCorr
GAMG:  Solving for p, Initial residual = 5.2808361e-06, Final residual = 8.6651752e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 1.3443669e-06, Final residual = 5.0129706e-07, No Iterations 1
time step continuity errors : sum local = 7.2516322e-09, global = -5.8100331e-10
Setting traction on solid patch
Total force (fluid2) = (-9.5286294 0.00083656471 -0.00029799967)
Total force (solid) = (9.5340998 -0.002029606 0.0002982776)
Evolving stress model
DICPCG: Solving for D, Initial residula = 1, Final residual = 0.077432708, No outer iterations = 11
Max relative residual = 1, Relative residual = 0.0002614466, enforceLinear = 0
Current fsi relative residual norm: 1

Time = 0.00025, iteration: 10
Maximal accumulated displacement of interface points: 9.6575407e+15
GAMG:  Solving for cellMotionUx, Initial residual = 1, Final residual = 0.00039886385, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 1, Final residual = 0.00027430961, No Iterations 2
GAMG:  Solving for cellMotionUz, Initial residual = 1, Final residual = 0.00077210207, No Iterations 2
GAMG:  Solving for cellMotionUx, Initial residual = 0.00035229159, Final residual = 6.4657305e-07, No Iterations 4
GAMG:  Solving for cellMotionUy, Initial residual = 0.00024006314, Final residual = 6.2628923e-07, No Iterations 4
GAMG:  Solving for cellMotionUz, Initial residual = 0.00038598462, Final residual = 3.2866915e-07, No Iterations 4
Evolving flow model
Courant Number mean: 16.222417 max: 1355476.9 velocity magnitude: 1.1257356e+16
DILUPBiCG:  Solving for Ux, Initial residual = 0.99108731, Final residual = 1.0439739, No Iterations 1000
DILUPBiCG:  Solving for Uy, Initial residual = 0.97611299, Final residual = 5.7019434, No Iterations 1000
DILUPBiCG:  Solving for Uz, Initial residual = 0.97059061, Final residual = 30.788802, No Iterations 1000

FSFI的第一个时间不长也算不过去，这是什么原因，没有使用并行运算也浮点报错

cfdwired · 2019 年6 月 13 日 00:10

到了 iteration = 10 那里，interface 最大的位移已经到了 9.6575407e+15，要检查一下边界条件和固体的物性常数

xye111 · 2019 年6 月 13 日 05:15

外侧流体	epsilon	k	nut	p	pointMotionU	U
INLET	inletOutlet	inletOutlet	calculated	zeroGradient	fixedValue	fixedValue
OUTLET	inletOutlet	inletOutlet	calculated	fixedValue	fixedValue	zeroGradient
FRONT_AND_BACK	epsilonWallFunction	kqRWallFunction	nutkWallFunction	zeroGradient	fixedValue	fixedValue
UP_AND_DOWN	epsilonWallFunction	kqRWallFunction	nutkWallFunction	zeroGradient	fixedValue	fixedValue
INTERFACE_INTERFACE	epsilonWallFunction	kqRWallFunction	nutkWallFunction	zeroGradient	fixedValue	movingWallVelocity

内侧流体	epsilon	k	nut	p	pointMotionU	U
INLET	inletOutlet	inletOutlet	calculated	zeroGradient	fixedValue	fixedValue
OUTLET	inletOutlet	inletOutlet	calculated	fixedValue	zeroGradient	zeroGradient
INTERFACEWATER	epsilonWallFunction	kqRWallFunction	nutkWallFunction	zeroGradient	fixedValue	movingWallVelocity

固体	D	pointD	密度	弹性模量	泊松比
OUTERWALL	tractionDisplacement	calculated	6500	7.80E+10	0.37
INNERWALL	tractionDisplacement	calculated
UP_AND_DOWN	fixedValue	fixedValue

chnrdu · 2019 年6 月 13 日 05:29

是不是内侧流体 OUTLET 的 pointMotionU 边界？你设置的是 zeroGradient?

xye111 · 2019 年6 月 13 日 05:33

是设置的零梯度。

chnrdu · 2019 年6 月 13 日 05:34

为什么要设置成 zeroGradient?

xye111 · 2019 年6 月 13 日 05:37

不然就用固定值。

chnrdu · 2019 年6 月 13 日 05:44

pointMotionU 是跟移动网格有关的。你这里的 OUTLET 是不用动的，所以还是用回 fixedValue.

xye111 · 2019 年6 月 13 日 05:54

Time = 0.0005, iteration: 1
Modes before clean-up : 8, modes after clean-up : 8
Maximal accumulated displacement of interface points: 1793919.4
GAMG:  Solving for cellMotionUx, Initial residual = 1, Final residual = 0.00060860902, No Iterations 3
GAMG:  Solving for cellMotionUy, Initial residual = 1, Final residual = 0.00017843867, No Iterations 5
GAMG:  Solving for cellMotionUz, Initial residual = 1, Final residual = 0.00032196348, No Iterations 4
GAMG:  Solving for cellMotionUx, Initial residual = 0.00034398031, Final residual = 6.8291781e-07, No Iterations 4
GAMG:  Solving for cellMotionUy, Initial residual = 6.4110158e-05, Final residual = 9.9743934e-07, No Iterations 3
GAMG:  Solving for cellMotionUz, Initial residual = 0.00015442846, Final residual = 4.5127317e-07, No Iterations 4
Evolving flow model
Courant Number mean: 4.3528161 max: 8558.3166 velocity magnitude: 3131165.9
^C
[xye@localhost fluid]$ ^C
[xye@localhost fluid]$

改了以后还是这样