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

Reaction Model

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

Like Tree1Likes
  • 1 Post By megacrout

Reply
 
LinkBack Thread Tools Display Modes
Old   August 17, 2011, 00:43
Default Reaction Model
  #1
New Member
 
Yopi
Join Date: Jul 2011
Posts: 7
Rep Power: 5
zovie is on a distinguished road
Hi everyone,

in openfoam201 i found different reaction model between chemFoam and reactingFoam (both combustion). chemFoam using chemkin and reactingFoam using something else.

how to change reactingfoam case, so that can read chemkin like chemFoam ?

tengkyu ...
zovie is offline   Reply With Quote

Old   August 17, 2011, 11:55
Default
  #2
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 5
megacrout is on a distinguished road
Hi Yopi,

The two applications for reading and using physical (Sutherland and Janaf coefficients) and chemical (reactions equations, preexponential factor, T-dependency of this latter one and activation energy) are chemkinReader and foamChemistryReader.

Moving from one model to the other is pretty tricky as they use different "syntax" and do not require the same units. This is mentioned and further explained here for example: reactingFoam constant/reactions file. So this gets messy if you have many reactions (which is necessary if you want to get realistic results) and no automatic way to changing these data.

So, answer to your question: change the file containing the reaction equations (usually chem.inp, always in constant/) as well as the name of the chemistry reader (reason why I wrote that first paragraph that seems at first to miss the point of your question) and the reaction file to be read if necessary. Those last two modifications have to be done in constant/thermophysicalProperties.

Im not an expert, but I think this is all you have to do.
Good luck.
zovie likes this.
megacrout is offline   Reply With Quote

Old   August 17, 2011, 23:06
Default
  #3
New Member
 
Yopi
Join Date: Jul 2011
Posts: 7
Rep Power: 5
zovie is on a distinguished road
thanks Tibo ... i am using chemkinToFoam and it works

at first i am confused how to compute highcpcoeff lowcpcoeff (in thermo.compressibleGas - reactingFoam). i think there are a table or equation for every particel or molecule. but now it solved.

thanks.
zovie is offline   Reply With Quote

Old   August 24, 2011, 04:18
Default
  #4
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 5
megacrout is on a distinguished road
Big correction: There actually is a nice utility called chemkinToFoam to convert your chemkin files to the foamChemistry format.

It requires the following syntax:
chemkinToFoam <CHEMKINChemistryFile> <CHEMKINThermodynamicsFile> <FOAMChemistryFile> <FOAMThermodynamicsFile>

The corresponding foamToChemkin utility does not seem to exist though.

enjoy.
megacrout is offline   Reply With Quote

Old   April 6, 2012, 16:25
Thumbs down
  #5
New Member
 
Mostafa Moghaddami
Join Date: Oct 2009
Posts: 13
Rep Power: 6
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:

------------------------------------------------------------------------------------
--> FOAM FATAL ERROR:
Inconsistent thermo package selected:

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

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

3
(
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

Reply

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
Eulerian multiphase model for species trnasport: UDF for reaction rate harsh Main CFD Forum 2 February 20, 2013 09:02
Wrong calculation of nut in the kOmegaSST turbulence model FelixL OpenFOAM Bugs 27 March 27, 2012 09:02
Turbulence Model and limitation to Reynolds number qascapri FLUENT 0 January 24, 2011 10:48
How to use the competing reaction model in CFX 11 NewCoalman CFX 0 May 5, 2008 03:15
discrete phase model with reaction Art FLUENT 0 September 7, 2004 10:51


All times are GMT -4. The time now is 03:52.