[sahm]

Dear FOAMers
Recently I've been trying to solve an axi-symmetric case of a diffusion flame. The case is based on this Paper by M.D. Smooke.

The reaction is 2 step-methane combustion defined in chemkin format.
Here I`ve uploaded the case files. Just have to blockMesh and run reactingFoam.
When I run this case, after some time, the solver stops giving this error:

#0 Foam::error::PrintStack(Foam::Ostream&) in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 __exp1 at interp.c:0
#4 __ieee754_pow at interp.c:0
#5 pow in "/lib64/libm.so.6"
#6 Foam::ODEChemistryModel<Foam::PsiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::PerfectGas> > > >::Omega(Foam::Reaction<Foam::sutherlandTranspor t< Foam::specieThermo<Foam::janafThermo<Foam::Perfect Gas> > > > const&, Foam::Field<double> const&, double, double, double&, double&, int&, double&, double&, int&) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#7 Foam::ODEChemistryModel<Foam::PsiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::PerfectGas> > > >::Omega(Foam::Field<double> const&, double, double) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#8 Foam::ODEChemistryModel<Foam::PsiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::PerfectGas> > > >::derivatives(double, Foam::Field<double> const&, Foam::Field<double>&) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#9 Foam::SIBS::SIMPR(Foam::ODE const&, double, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::SquareMatrix<double> const&, double, int, Foam::Field<double>&) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libODE.so"
#10 Foam::SIBS::solve(Foam::ODE const&, double&, Foam::Field<double>&, Foam::Field<double>&, double, Foam::Field<double> const&, double, double&, double&) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libODE.so"
#11 Foam::ODESolver::solve(Foam::ODE const&, double, double, Foam::Field<double>&, double, double&) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libODE.so"
#12 Foam::Ode<Foam::PsiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::PerfectGas> > > >::solve(Foam::Field<double>&, double, double, double, double) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#13 Foam::ODEChemistryModel<Foam::PsiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam::PerfectGas> > > >::solve(double, double) in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#14 main in "/home/sahm/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"
#15 __libc_start_main in "/lib64/libc.so.6"
#16 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"
Floating point exception

This error looks similar to what I have seen from other people in this forum, but this is in chemistry. I know something goes wrong in the chemistry solver, but I have no idea what's the problem? Can anybody help me with this?

[markusrehm]

Hello,

you can try to use another ODE solver first, reduce the time step and switch off chemistry. If that doesn't help try to run it in debug mode to find out where it died.

Regards, Markus.

[sahm]

Ok, The problem with chemistry was my Chemkin File, I had to include only the forward reaction for CO, using <=> means that it is a reversible reaction.
Now I have this new error:

Solving chemistry


attempt to use janafThermo<equationOfState> out of temperature range 200 -> 5000; T = 5074.26#0 Foam::error:rintStack(Foam::Ostream&) in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> >::H(double) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#3 Foam::ODEChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >::solve(double, double) in "/home/sahm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#4 main in "/home/sahm/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"
#5 __libc_start_main in "/lib64/libc.so.6"
#6 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/sahm/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"


From function janafThermo<equationOfState>::checkT(const scalar T) const
in file /home/sahm/OpenFOAM/OpenFOAM-1.6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 64.

FOAM aborting

which means my temperature is going out of range. I have seen this error before, but I don't know how to solve it.
Any help is appreciated.

[markusrehm]

Hello Sahm,

there is no easy solution and the reasons for the temperature going out of range can be VERY different (might be the solver step, the chemical mechanism, the ODE sover or just a wrong flow profile).

Please follow the hints posted by me before here:

http://www.cfd-online.com/Forums/ope...and-error.html

Generally speaking deactivate all models that make problems and successively activate them.

Furthermore you should run the case in debug - mode and observe all relevant quantities.

Regards, Markus.

[sahm]

Markus
I have checked many things and finaly I came to this point that I think the solver of reactingFoam is wrong in versions 1.5 and 1.6, but it is changed and works properly in 1.6.x and 1.7. I have described my problems and tests in these links:
Trying to figure them out.
I guess its diffusion.
Also there is one test case that I can send to you in mail.
Alberto helped a lot, and Finally I can run my cases in 1.7, But I have a new problem.
My Lab-mate and I are trying to use your cantera interface in 1.5, but since the hEqn is not right( based on what I have discussed with other people) I have to correct it. At this point, I need to have access to hi ( enthalpy of each specie) Can you help me with it? Is the species and mixture thermophysical models different in cantera interface than OpenFoam interface?