使用multiphaseInterFoam时,加入了动网格,出现了不收敛的情况
-
使用multiphaseInterFoam时,用snappyHexMesh设置了一块区域,然后加入了动网格,让这块区域直线运动,结果在0.65s出现了不收敛的情况,这是log文件,能不能看出不收敛的原因是什么,有大神们可以解答一下吗
Courant Number mean: 0.00263715 max: 0.592641 Interface Courant Number mean: 0.000205746 max: 0.196717 deltaT = 1.86717e-105 Time = 0.671233 PIMPLE: iteration 1 GAMG: Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0 GAMG: Solving for cellDisplacementy, Initial residual = 4.57288e-06, Final residual = 4.57288e-06, No Iterations 0 Execution time for mesh.update() = 0.02 s GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 3.45004e-06, No Iterations 5 time step continuity errors : sum local = 1.13181e-13, global = -2.30621e-14, cumulative = -0.0560662 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.42526e-64 1 Phase-sum volume fraction, min, max = 1 1 1 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.4271e-64 1 Phase-sum volume fraction, min, max = 1 1 1 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.42894e-64 1 Phase-sum volume fraction, min, max = 1 1 1 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.43078e-64 1 Phase-sum volume fraction, min, max = 1 1 1 GAMG: Solving for p_rgh, Initial residual = 0.0133035, Final residual = 0.000493626, No Iterations 2 time step continuity errors : sum local = 4.2474e-05, global = 3.00597e-06, cumulative = -0.0560632 GAMGPCG: Solving for p_rgh, Initial residual = 0.00231444, Final residual = 3.46657e-08, No Iterations 6 time step continuity errors : sum local = 1.5629e-09, global = 6.70208e-11, cumulative = -0.0560632 ExecutionTime = 147.84 s ClockTime = 148 s Courant Number mean: 0.0025889 max: 0.575287 Interface Courant Number mean: 0.000207612 max: 0.194023 deltaT = 1.62282e-105 Time = 0.671233 PIMPLE: iteration 1 GAMG: Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0 GAMG: Solving for cellDisplacementy, Initial residual = 4.57288e-06, Final residual = 4.57288e-06, No Iterations 0 Execution time for mesh.update() = 0.02 s GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.87234e-06, No Iterations 7 time step continuity errors : sum local = 4.07967e-13, global = -9.83711e-14, cumulative = -0.0560632 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.4315e-64 1 Phase-sum volume fraction, min, max = 1 1 1 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.43222e-64 1 Phase-sum volume fraction, min, max = 1 1 1 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.43294e-64 1 Phase-sum volume fraction, min, max = 1 1 1 MULES: Solving for alpha.water water volume fraction, min, max = 0 0 0 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 0 MULES: Solving for alpha.salt salt volume fraction, min, max = 0.696714 0 1 MULES: Solving for alpha.air air volume fraction, min, max = 0.303286 3.43366e-64 1 Phase-sum volume fraction, min, max = 1 1 1 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libc.so.6 #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:? #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:? #9 ? in ~/OpenFOAM2006/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam #10 ? in /lib/x86_64-linux-gnu/libc.so.6 #11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #12 ? in ~/OpenFOAM2006/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/multiphaseInterFoam Floating point exception (core dumped)