CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Error calculating thermodynamic properties in reactingFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 16, 2009, 05:43
Default Error calculating thermodynamic properties in reactingFoam
  #1
Member
 
Michele Vascellari
Join Date: Mar 2009
Posts: 70
Rep Power: 17
mighelone is on a distinguished road
Dear all,

i'm trying to solve a common test-case for gas combustion (the CH4 bluff-body flame of Sandia).

I'm considering a skeletal mechanism with 16 species and 41 reactions (but I've obtained the same errors using a simple one-step mechanism) for CH4 combustion. I've solved the cold flow, and I've patched the initial temperature field, using setField creating an ignition zone near the bluff body with a temperature of 1000 K.

After some iterations I've obtained the following error message:

Code:
DILUPBiCG:  Solving for h, Initial residual = 0.993298, Final residual = 4.85856e-07, No Iterations 58
[0] 
[0] 
[0] attempt to use janafThermo<equationOfState> out of temperature range 200 -> 3000;  T = 170.194#0  Foam::error::printStack(Foam::Ostream&) in "/home/michele/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::error::abort() in "/home/michele/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> >::H(double) const in "/home/michele/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#3  Foam::hPsiMixtureThermo<Foam::reactingMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > > >::calculate() in "/home/michele/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#4  Foam::hPsiMixtureThermo<Foam::reactingMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > > >::correct() in "/home/michele/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#5  main in "/home/michele/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"
It seems an error calculating the temperature field from enthalpy.

Any idea how to solve this problem?

I'm starting from the reactingFoam tutorials to define chemistry dictionary.

Thanks in advance for the attention

Michele
mighelone is offline   Reply With Quote

Old   October 19, 2009, 10:47
Default
  #2
Member
 
Michele Vascellari
Join Date: Mar 2009
Posts: 70
Rep Power: 17
mighelone is on a distinguished road
No ideas how to solve this problem?

I'm trying to modify some simulation parameters with no success, like for example the time step and others.

Thank you for the attention

Michele
mighelone is offline   Reply With Quote

Old   October 21, 2009, 05:40
Default
  #3
New Member
 
Roger Jove
Join Date: Mar 2009
Location: Barcelona
Posts: 4
Rep Power: 17
roger.jove is on a distinguished road
Hi Michele,

I am sorry but going out of temperature range is a quite common problem. I have found this problem with different solvers that use janafThermo.

Sometimes I have solved it using a finer mesh. I have tried other options that I found in the link below but they didn't worked for me (I am still a newbie!). However, take a look, it will be provably interesting for you.


Roger


http://www.cfd-online.com/Forums/ope...naf-range.html
roger.jove is offline   Reply With Quote

Old   October 21, 2009, 06:56
Default
  #4
Member
 
Michele Vascellari
Join Date: Mar 2009
Posts: 70
Rep Power: 17
mighelone is on a distinguished road
Hi Roger,

thank you for your suggestion, in the topic I've found some interesting suggestions to solve the problem.

I guess the CHEMKIN chemistry solver is not enough robust to work at low temperature range, that could be found during the ignition of combustion problem.

Probably considering another chemistry solver (http://openfoamwiki.net/index.php/Co...ateReactinFoam) the problem can be solved.

Michele
mighelone is offline   Reply With Quote

Old   October 21, 2009, 10:25
Default
  #5
Member
 
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 41
Rep Power: 17
dohnie is on a distinguished road
Mighelone,
can you locate where in your field this temperature occurs?
A non-suitable boundary condition (reaction close to a BC interfers with BC) is a frequent error source.
Moreover, a too high Courant Number in reactingFoam can also give problems with the reaction.
dohnie is offline   Reply With Quote

Old   October 21, 2009, 17:47
Default
  #6
Member
 
Michele Vascellari
Join Date: Mar 2009
Posts: 70
Rep Power: 17
mighelone is on a distinguished road
I'm simulating a bluff body flame. The flame is developed downstream of the bluff body. I have modeled a little part of the duct upstream, therefore the flame is not attached to the inlet boundaries.

I've tried to modifiy the Courant number, reducing it until 0.01, but without success. On the contrary the Janaf error seems to appear before reducing the Courant number.

Tomorrow I'm trying to modify the mesh, reducing the size and improving the quality.

Michele
mighelone is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
Laminar diffusion flames and reactingFoam Transport properties jgaricano OpenFOAM Running, Solving & CFD 0 June 4, 2008 17:58
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 13:02
Thermodynamic Properties of fluids at high temp hagupta CFX 3 March 15, 2006 13:55
thermodynamic properties Jorge Main CFD Forum 1 November 1, 2004 09:55


All times are GMT -4. The time now is 06:27.