目前基于两相流求解器fsifoam,BeamInCross的案例可以运行。
在设置新的算例过程中出现问题,检查网格设置没有问题,但在运行到wavefsifoam中的时候总会出现浮点报错,以下是报错的log文件。有大佬知道是什么问题嘛?
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create dynamic mesh for time = 0
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Reading stress mesh
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
PIMPLE: no residual control data found. Calculations will employ 2 corrector loops
Creating traction displacement boundary conditions
Creating fixed displacement boundary condition
Selecting rheology model linearElastic
Creating constitutive model
Selecting coupling scheme IQN-ILS
Starting time loop
Time = 0.0005 (dt = 0.0005)
Create extended GGI zone-to-zone interpolator
Checking fluid-to-solid face interpolator
Fluid-to-solid face interpolation error: 2.66686e-15
Checking solid-to-fluid point interpolator (GGI)
calcMasterPointAddressing() const
Extended GGI, master point distance, max: 0, avg: 0, min: 0
Solid-to-fluid point interpolation error (GGI): 1.78004e-15
Number of uncovered master faces: 0
Number of uncovered slave faces: 0
calcSlavePointAddressing() const
Extended GGI, slave point distance, max: 0, avg: 0, min: 0
Extended GGI, slave point orientation (<0), max: -0.707107, min: -1, nIncorrectPoints: 0/416
Setting traction on solid patch
Total force (fluid) = (9.09051e-13 -10310.3 0)
Total force (solid) = (-9.17932e-13 10310.3 0)
Evolving solid solver: unsTotalLagrangianSolid
DICPCG: Solving for D, Initial residula = 1, Final residual = 0.00139832, No outer iterations = 46
Max relative residual = 1, Relative residual = 8.57078e-05, enforceLinear = 0
Current fsi relative residual norm: 1
Alternative fsi residual: 1
Time = 0.0005, iteration: 1
Modes before clean-up : 0, modes after clean-up : 0
Current fsi under-relaxation factor: 0.4
Creating accumulated fluid interface displacement
Maximal accumulated displacement of interface points: 0.0176225
GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 0.000274813, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 1, Final residual = 3.40719e-06, No Iterations 1
GAMG: Solving for cellMotionUz, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving fluid solver
Courant Number mean: 0 max: 0 velocity magnitude: 0
volume continuity errors : volume = 1.798, max error = 7.41202e-10, sum local = 3.34731e-14, global = 3.74747e-17
time step continuity errors : sum local = 6.76086e-06, global = -4.79613e-09, cumulative = -4.79613e-09
GAMGPCG: Solving for pcorr: solution singularity
GAMGPCG: Solving for pcorr: solution singularity
time step continuity errors : sum local = 0, global = 0, cumulative = -4.79613e-09
Courant Number mean: 0 max: 0 velocity magnitude: 0
PIMPLE: iteration 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.333151 Min(alpha1_) = 0 Max(alpha1_) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.33315 Min(alpha1_) = -5.43647e-12 Max(alpha1_) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.333149 Min(alpha1_) = -1.63119e-11 Max(alpha1_) = 1
Relaxing time: 0.13 s
DILUPBiCG: Solving for Ux, Initial residual = 0.0139054, Final residual = 5.6888e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.90048, Final residual = 7.68278e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 3.12111e-08, No Iterations 1
GAMG: Solving for pd, Initial residual = 0.000490799, Final residual = 3.64458e-07, No Iterations 6
GAMG: Solving for pd, Initial residual = 5.23309e-07, Final residual = 1.71699e-07, No Iterations 2
time step continuity errors : sum local = 6.15505e-08, global = 2.24553e-08, cumulative = 1.76592e-08
PIMPLE: iteration 2
MULES: Solving for alpha1
Liquid phase volume fraction = 0.33315 Min(alpha1_) = -1.87448e-09 Max(alpha1_) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.333149 Min(alpha1_) = -1.87447e-09 Max(alpha1_) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.333148 Min(alpha1_) = -1.87447e-09 Max(alpha1_) = 1.00001
Relaxing time: 0.06 s
DILUPBiCG: Solving for Ux, Initial residual = 0.00171359, Final residual = 7.95819e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00229415, Final residual = 1.1053e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.971192, Final residual = 3.40668e-10, No Iterations 2
GAMG: Solving for pd, Initial residual = 0.000164208, Final residual = 8.28008e-07, No Iterations 1
GAMGPCG: Solving for pd, Initial residual = 8.28209e-07, Final residual = 3.38658e-09, No Iterations 5
time step continuity errors : sum local = 1.20621e-09, global = -9.02852e-12, cumulative = 1.76502e-08
Setting traction on solid patch
Total force (fluid) = (-0.05684 4.2698e+06 0.000283692)
Total force (solid) = (0.0568298 -4.2698e+06 -0.000283692)
Evolving solid solver: unsTotalLagrangianSolid
DICPCG: Solving for D, Initial residula = 0.99665, Final residual = 0.00109644, No outer iterations = 44
Max relative residual = 1.01269, Relative residual = 9.25099e-05, enforceLinear = 0
Current fsi relative residual norm: 1
Alternative fsi residual: 1.00254
Time = 0.0005, iteration: 2
Current fsi under-relaxation factor: 0.4
Maximal accumulated displacement of interface points: 3.28076
GAMG: Solving for cellMotionUx, Initial residual = 0.993157, Final residual = 0.000144068, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 1, Final residual = 8.88797e-06, No Iterations 1
GAMG: Solving for cellMotionUz, Initial residual = 6.05144e-07, Final residual = 3.33872e-07, No Iterations 1
Evolving fluid solver
Courant Number mean: 0.00051952 max: 0.0335853 velocity magnitude: 0.33609
volume continuity errors : volume = 1.79801, max error = 8.29118e-10, sum local = 3.35534e-14, global = 3.66141e-17
time step continuity errors : sum local = 0.00099525, global = 3.75254e-06, cumulative = 3.75254e-06
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.62621e-05, No Iterations 8
GAMGPCG: Solving for pcorr, Initial residual = 2.80297e-05, Final residual = 1.4592e-06, No Iterations 2
time step continuity errors : sum local = 9.21979e-12, global = -8.2965e-14, cumulative = 3.75254e-06
Courant Number mean: 0.000181841 max: 0.0268574 velocity magnitude: 0.26901
PIMPLE: iteration 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.332773 Min(alpha1_) = -4.11526e-33 Max(alpha1_) = 1.00385
MULES: Solving for alpha1
Liquid phase volume fraction = 0.332869 Min(alpha1_) = 0 Max(alpha1_) = 1.00315
MULES: Solving for alpha1
Liquid phase volume fraction = 0.332965 Min(alpha1_) = 0 Max(alpha1_) = 1.00257
Relaxing time: 0.07 s
DILUPBiCG: Solving for Ux, Initial residual = 0.00750013, Final residual = 1.9279e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.476658, Final residual = 6.03497e-09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.363964, Final residual = 1.19961e-09, No Iterations 2
GAMG: Solving for pd, Initial residual = 0.0552688, Final residual = 4.74688e-07, No Iterations 21
GAMG: Solving for pd, Initial residual = 1.30893e-09, Final residual = 1.80425e-09, No Iterations 1
time step continuity errors : sum local = 6.42096e-08, global = -1.01415e-09, cumulative = 3.75153e-06
PIMPLE: iteration 2
MULES: Solving for alpha1
Liquid phase volume fraction = 0.332835 Min(alpha1_) = 0 Max(alpha1_) = 1.00023
MULES: Solving for alpha1
Liquid phase volume fraction = 0.332993 Min(alpha1_) = 0 Max(alpha1_) = 1.22181
MULES: Solving for alpha1
Liquid phase volume fraction = 0.333151 Min(alpha1_) = 0 Max(alpha1_) = 1.40207
Relaxing time: 0.06 s
DILUPBiCG: Solving for Ux, Initial residual = 0.164439, Final residual = 8.32138e-08, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.201356, Final residual = 6.76597e-07, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.834862, Final residual = 8.61714e-07, No Iterations 5
GAMG: Solving for pd, Initial residual = 3.4311e-05, Final residual = 9.63653e-07, No Iterations 8
GAMGPCG: Solving for pd, Initial residual = 1.72276e-06, Final residual = 3.34751e-09, No Iterations 4
time step continuity errors : sum local = 6.58196e-08, global = 1.14823e-10, cumulative = 3.75164e-06
Setting traction on solid patch
Total force (fluid) = (44895.8 -5.16199e+08 -0.04222)
Total force (solid) = (-44895 5.16233e+08 0.0422304)
Evolving solid solver: unsTotalLagrangianSolid
det: -6.35377, 1.44569
DICPCG: Solving for D, Initial residula = 0.981674, Final residual = 0.000998539, No outer iterations = 46
Max relative residual = 1.04638, Relative residual = 0.000100841, enforceLinear = 0
Current fsi relative residual norm: 1
Alternative fsi residual: 1.00299
Time = 0.0005, iteration: 3
Current fsi under-relaxation factor: 0.4
Maximal accumulated displacement of interface points: 468.826
GAMG: Solving for cellMotionUx, Initial residual = 0.999092, Final residual = 0.000110299, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 1, Final residual = 9.9143e-06, No Iterations 1
GAMG: Solving for cellMotionUz, Initial residual = 0.00878954, Final residual = 5.34767e-06, No Iterations 7
Evolving fluid solver
Courant Number mean: 0.177766 max: 809.001 velocity magnitude: 163.989
volume continuity errors : volume = 1.79747, max error = 2.69711e-06, sum local = 3.54005e-14, global = 2.8054e-17
time step continuity errors : sum local = 0.136847, global = -0.00030095, cumulative = -0.00030095
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 8.68671e-05, No Iterations 67
GAMGPCG: Solving for pcorr, Initial residual = 0.478654, Final residual = 0.00915489, No Iterations 100
time step continuity errors : sum local = -1.61429e-05, global = -1.66005e-06, cumulative = -0.00030261
Courant Number mean: 6.3896 max: 548138 velocity magnitude: 12675.3
PIMPLE: iteration 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.376351 Min(alpha1_) = -13.5602 Max(alpha1_) = 7.08985
MULES: Solving for alpha1
Liquid phase volume fraction = 0.361769 Min(alpha1_) = -295.669 Max(alpha1_) = 154.958
MULES: Solving for alpha1
Liquid phase volume fraction = 0.345683 Min(alpha1_) = -6379.38 Max(alpha1_) = 5395.01
Relaxing time: 0.07 s
DILUPBiCG: Solving for Ux, Initial residual = 0.580998, Final residual = 2.92065e-07, No Iterations 17
DILUPBiCG: Solving for Uy, Initial residual = 0.977539, Final residual = 9.94289e-07, No Iterations 14
DILUPBiCG: Solving for Uz, Initial residual = 0.972502, Final residual = 9.01328e-08, No Iterations 6
Floating point exception (core dumped)