CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   Reaction Model (

zovie August 17, 2011 00:43

Reaction Model
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 ...

megacrout August 17, 2011 11:55

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: 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 August 17, 2011 23:06

thanks Tibo ... i am using chemkinToFoam and it works :D

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.


megacrout August 24, 2011 04:18

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.


Mostafa April 6, 2012 16:25

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.

All times are GMT -4. The time now is 00:41.