总是算到一定时间步,计算停止报错,有哪位大神知道这样的报错信息是问题出在哪吗?是发散的问题吗
-
kEqnCoeffs { Ce 1.048; Ck 0.094; } Reading g Reading hRef Calculating field g.h No MRF models present No finite volume options present GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 8.8753743e-06, No Iterations 8 GAMGPCG: Solving for pcorr, Initial residual = 0.024255991, Final residual = 2.1939681e-06, No Iterations 3 time step continuity errors : sum local = 4.7008361e-15, global = -6.386216e-17, cumulative = -6.386216e-17 Courant Number mean: 0.001130128 max: 63.68082 Starting time loop Courant Number mean: 0.001130128 max: 63.68082 Interface Courant Number mean: 7.3912381e-05 max: 9.8596523 Time = 0.0021 PIMPLE: iteration 1 MULES: Solving for alpha.water Phase-1 volume fraction = 0.50002593 Min(alpha.water) = -2.5864875e-05 Max(alpha.water) = 1.0000001 MULES: Solving for alpha.water Phase-1 volume fraction = 0.50002492 Min(alpha.water) = -0.010606339 Max(alpha.water) = 1.0000001 MULES: Solving for alpha.water Phase-1 volume fraction = 0.50002392 Min(alpha.water) = -0.0074242149 Max(alpha.water) = 1.0000001 MULES: Solving for alpha.water Phase-1 volume fraction = 0.50002292 Min(alpha.water) = -0.0078825052 Max(alpha.water) = 1.0000001 #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::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:? #4 Foam::DICSmoother::DICSmoother(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) at ??:? #5 Foam::lduMatrix::smoother::addsymMatrixConstructorToTable<Foam::DICSmoother>::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) at ??:? #6 Foam::lduMatrix::smoother::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::dictionary const&) at ??:? #7 Foam::GAMGSolver::initVcycle(Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::lduMatrix::smoother>&, Foam::Field<double>&, Foam::Field<double>&) const at ??:? #8 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #9 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? #10 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? #11 ? at ??:? #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? at ??:? 浮点数例外 (核心已转储)
求助,跪谢!