求助，目前在使用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"
