CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

ReactingFoam Error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 3, 2010, 20:21
Default ReactingFoam Error
  #1
Member
 
sahm's Avatar
 
S. Ali H.M.
Join Date: Nov 2009
Location: Chicago
Posts: 69
Rep Power: 7
sahm is on a distinguished road
Send a message via Yahoo to 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?
Attached Files
File Type: zip Bunsen-3-Reaction.zip (85.6 KB, 68 views)
__________________
SAHM
sahm is offline   Reply With Quote

Old   August 4, 2010, 04:55
Default
  #2
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 176
Rep Power: 8
markusrehm is on a distinguished road
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.
markusrehm is offline   Reply With Quote

Old   August 5, 2010, 12:32
Default New problem
  #3
Member
 
sahm's Avatar
 
S. Ali H.M.
Join Date: Nov 2009
Location: Chicago
Posts: 69
Rep Power: 7
sahm is on a distinguished road
Send a message via Yahoo to 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.
__________________
SAHM
sahm is offline   Reply With Quote

Old   August 6, 2010, 03:41
Default
  #4
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 176
Rep Power: 8
markusrehm is on a distinguished road
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:

Cantera: Help me understand this error

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.
markusrehm is offline   Reply With Quote

Old   August 20, 2010, 15:22
Default New problem
  #5
Member
 
sahm's Avatar
 
S. Ali H.M.
Join Date: Nov 2009
Location: Chicago
Posts: 69
Rep Power: 7
sahm is on a distinguished road
Send a message via Yahoo to 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?
__________________
SAHM
sahm is offline   Reply With Quote

Reply

Tags
chemkin, methane, openfoam, reactingfoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM 300 October 29, 2014 19:00
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08
Compiling Netgen on Fedora Core is driving me crazy jango Open Source Meshers: Gmsh, Netgen, CGNS, ... 3 November 9, 2007 14:29
user defined function cfduser CFX 0 April 29, 2006 10:58


All times are GMT -4. The time now is 07:53.