Supercritical flow
 

This is a supercritical flow （O2 、H2）, but it always diverges in the calculation process：疯狂的：：疯狂的：. I hope you can give me some suggestions. Thank you very much.


：
Courant Number mean: 2.51125 max: 18.699
Time = 0.0003855

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
min/max(U)) = 0, 149.098
update D with binary diffusion model...
species transprotation equation
DILUPBiCG: Solving for O2, Initial residual = 0.00116254, Final residual = 2.02936e-06, No Iterations 3
DILUPBiCG: Solving for H2, Initial residual = 0.000823762, Final residual = 5.40367e-07, No Iterations 3
Calculating inert species N2
Solving energy equation
DILUPBiCG: Solving for h, Initial residual = 0.00355332, Final residual = 4.344e-09, No Iterations 6
Starting correcting themo
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 in "/lib/x86_64-linux-gnu/libm.so.6"
#4 pow in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::PengRobinson<Foam::realGasSpecie>::Z(double, double) const at ??:?
#6 Foam::realGasJanafThermo<Foam::PengRobinson<Foam:: realGasSpecie> >::cp(double, double) const at ??:?
#7 Foam::realGasHePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::calculate() at ??:?
#8 Foam::realGasHePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::correct() at ??:?
#9
at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
at ??:?

> Courant Number mean: 2.51125 max: 18.699

You should not have Courant numbers beyond 1. I recommend even lower values.

Use the automatic time step control which ensures that.

Thank you very much. I have made changes based on your suggestion. Unfortunately, the result is still divergent.

Courant Number mean: 0.0499972 max: 0.11485
Time = 0.00164182

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
min/max(U)) = 0, 43.6587
update D with binary diffusion model...
species transprotation equation
DILUPBiCG: Solving for O2, Initial residual = 4.01292e-06, Final residual = 6.74252e-19, No Iterations 3
DILUPBiCG: Solving for H2, Initial residual = 3.93064e-06, Final residual = 2.12152e-18, No Iterations 3
Calculating inert species N2
Solving energy equation
DILUPBiCG: Solving for h, Initial residual = 2.2811e-05, Final residual = 1.54394e-10, No Iterations 3
Starting correcting themo
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 in "/lib/x86_64-linux-gnu/libm.so.6"
#4 pow in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::PengRobinson<Foam::realGasSpecie>::Z(double, double) const at ??:?
#6 Foam::realGasJanafThermo<Foam::PengRobinson<Foam:: realGasSpecie> >::cp(double, double) const at ??:?
#7 Foam::realGasHePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::calculate() at ??:?
#8 Foam::realGasHePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::correct() at ??:?
#9
at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
at ??:?
浮点数例外 (核心已转储)