CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

simulating comb with reactingFoam

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

LinkBack Thread Tools Display Modes
Old   January 17, 2011, 13:04
Default simulating comb with reactingFoam
Josť Rodrigues
Join Date: Jun 2010
Location: IN+/IST Lisbon
Posts: 53
Rep Power: 9
jose_rodrig is on a distinguished road

Im trying to simulate combustion inside a CAN Type Gas Turbine Combustion Chamber with methane (CH4).

In case you know, the input conditions and the mesh are similar to the 3-step EDU tutorial from the straccm+ manual. If thats not the case, shortly, the setup consists of a 60 degree wedge with 2 swirling air inlets to the main chamber, 1 swirling fuel(CH4) inlet and one 3rd air inlet running in an outer channel for cooling. Both the 2 air inlets and the fuel inlet are characterized with cylindricalInletVelocity BC, supplied with OF-1.7.1.

The Can is 1m long and .25m radius and the domain only represents a 60 degree wedge section of the Can. The mesh as 1800 cells (coarse mesh). The wedge sides have symmetryPlane BC.

I tried to run this simulation with reactingFoam solver modified to account ignition (as purposed in the OFwiki).
Plus, most of the settings of the several dicts were taken from the counterFlowFlame2D tutorial (such as initialChemicalTimeStep and other coeff). RANS was used for turbulence modeling.

At 1st, I used the ode chemistrySolver in the chemistryProperties dict but i got the janaf out of range error. The temp exceeded the 5000k limit.

Next, I used the sequential chemistySolver and then I got a solution. However, the residuals were very unstable with a courant no. of 0.5 and the MAXtemp would oscillate between 2600k and 3100k.

1. What are the basis of either the sequential and the ode chemistrySolvers? Can you give some references to this models?

2. Which Courant number should I set?

3. Can I use different ddt schemes such as CrankNicholson with reactingFoam? What about steadyState?

4. Which settings would you recommend to tune so I can achieve convergence. I m avoiding refining the mesh to keep the setup as close as possible to the tut in starccm+.

5. Is it possible to use dieselFoam without injectors or spray?

That is all =P
I sincerely hope this thread will be a contribution to the OFOAM community!
Thank you.
jose_rodrig is offline   Reply With Quote

Old   January 20, 2011, 05:47
Senior Member
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 9
nakul is on a distinguished road

This JANAF error is quite a common problem with reactingFoam. It occurs mainly due to wrong case setup and there could be many reasons for setup to go wrong,

As for your questions:

1) Its advisable that you check out the source code in "src/thermophysical/chemistry" to find out the details of these chemistry solvers. If however you come across some reference do post it here. I would also be interested in goin through them.

2) Based on my experience I have found that for reactingFoam, max Co=0.1 is best suited.

3) There could be some stability issues if you change the ddt scheme. But you can give it a try. Also for steady state a separate solver called steadyReactingFoam is available.

4) You may do a "checkMesh" to see if everything is alright with your mesh. If not then you would have to change it. Apart from this to achieve convergence have a look at your BC and IC. They are the primary source of "JANAF" error.

5) Yes I think there shouldn't be any problem in doing that. Its code is exactly similar to that of reactingFoam. You may give it a try too.
nakul is offline   Reply With Quote

Old   February 7, 2011, 12:16
Post Ignition method
New Member
Join Date: Oct 2010
Posts: 1
Rep Power: 0
Tang is on a distinguished road
Hi Josť,

I'm also trying to simulate a one step methane combustion inside a cylindrical burner with reactingFoam and I'd like to modify this solver to account for ignition. I modified it as mentioned on the OFwiki and run the case it runs without error but not ignition occurs. Apparently the combustionProperties file is not read at all since when I start the run it get the following message, and (correct me if I'm wrong) it should be read together with "chemistryProperties" "g" and "thermophysicalProperties":

/************************************************** ************************************************** ***************/
Create time

Create mesh for time = 0

Reading chemistry properties

Reading g

Reading thermophysicalProperties
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysic s>>
Selecting chemistryReader foamChemistryReader
Selecting chemistrySolver ode
Selecting ODE solver SIBS
ODEChemistryModel: Number of species = 5 and reactions = 1
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.
/************************************************** ************************************************** ***************/
Did you encounter the same problem? Is the procedure also valid for OF-1.7.7? You should know that I'm quite unfamiliar with the use of OF and C++.

Thanks a lot for any kind of help
Tang is offline   Reply With Quote

Old   February 7, 2011, 12:36
Josť Rodrigues
Join Date: Jun 2010
Location: IN+/IST Lisbon
Posts: 53
Rep Power: 9
jose_rodrig is on a distinguished road
Hi Tang,

You re right: your solver is not reading the combustionProperties files. The beguining of the solver output should be something like:

Create time

Create mesh for time = 0

Reading chemistry properties

Reading combustion properties

Found ignition cells:

Ignition on

Reading g


I recommend that you recheck all the steps from the OFwiki tutorial. Make sure you use a different name for the solver like reactingFoamIgnition (as i did).

jose_rodrig is offline   Reply With Quote

Old   April 6, 2012, 16:31
Thumbs down
New Member
Mostafa Moghaddami
Join Date: Oct 2009
Posts: 13
Rep Power: 9
Mostafa is on a distinguished road
Dear All

I am going to run a case includes some liquid species that react with each other. I am going to use the reactingFoam as a solver.
All species have constant proprties and I want to use the following thermo type:
thermoType hsPsiThermo<reactingMixture<constTransport<specieT hermo<eConstThermo<icoPolynomial>>>>>>

but when I set it in the constat/thermophysicalProperties file and run the case, I get this error:

Inconsistent thermo package selected:

hsPsiThermo<reactingMixture<constTransport<specieT hermo<eConstThermo<icoPolynomial>>>>>>

Please select a thermo package based on gasThermoPhysics. Valid options include:

hsPsiMixtureThermo<singleStepReactingMixture<gasTh ermoPhysics>>
hsPsiMixtureThermo<multiComponentMixture<gasThermo Physics>>
hsPsiMixtureThermo<reactingMixture<gasThermoPhysic s>>

From function autoPtr<hsCombustionThermo> hsCombustionThermo::NewType(const fvMesh&, const word&)
in file combustionThermo/hsCombustionThermo/hsCombustionThermoNew.C at line 116.

FOAM exiting

What should I do to use the reactingFoam for constant propertie species?

Thanks in advance for your help.
Mostafa is offline   Reply With Quote


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
reactingFoam wedge handling wrong U dhondupant OpenFOAM Bugs 1 December 9, 2010 08:34
reactingFoam in openfoam1.7.1 Farshad_Noravesh OpenFOAM 10 November 30, 2010 03:38
reactingFoam floating point exception pajofego OpenFOAM 0 November 6, 2010 18:29
Boundary condition setting for non-premixed combustion using reactingFoam skyopener OpenFOAM 0 May 23, 2010 22:55
reactingFoam - turbulent reacting flow hamburgFoam OpenFOAM 0 December 7, 2009 13:57

All times are GMT -4. The time now is 20:25.