在fsifoam中运行算例报错
-
目前基于两相流求解器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) -
-
fields是场松弛，equations是方程松弛。在simpleFoam解析里面有提及 http://dyfluid.com/simpleFoam.html
