OpenFOAM计算复杂外形报错
-
求助，目前在使用hy2Foam求解器计算类似乘波体外形，已经计算过钝头算例可以得到收敛解。但目前计算外形相对复杂，二维计算，计算域为1*0.55m，网格数量42000。计算过程中经常报错，求问有没有什么调试建议。下面是报错信息。
Mean and max Courant Numbers = 5.29417e-05 0.0491796 deltaT = 3.99074e-11 Time = 2.017990284e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000161695, Final residual = 3.0494e-13, No Iterations 1 DILUPBiCG: Solving for O2, Initial residual = 0.000148597, Final residual = 1.8029e-13, No Iterations 1 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 1.07385e-06, Final residual = 1.49309e-06, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 8.83461e-07, Final residual = 1.19425e-06, No Iterations 1000 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 2.94288e-05, Final residual = 3.15326e+222, No Iterations 400 [1;31mAttempt to use rho2ReactionThermo out of temperature range 11152 times during this iteration. -1.64159e+236 < Tlow: 200[0m Phase no 1.0 ExecutionTime = 3564.14 s ClockTime = 3564 s Iteration no 4641 (0.79 s) Mean and max Courant Numbers = 5.41311e-05 0.0489379 deltaT = 4.07806e-11 Time = 2.018031064e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000162296, Final residual = 1.67236e-16, No Iterations 2 DILUPBiCG: Solving for O2, Initial residual = 0.000148383, Final residual = 6.70388e-18, No Iterations 2 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 #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::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:? #9 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #10 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #11 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"