Chemkin units

Rickard.Solsjo · September 1, 2011, 03:50

Hi
Im using chemkin for combustion in dieselFoam and need to change the pre-exponential factor and activation energy for the reactions.
However trying to understand the default one is hard since there are no references to the literature for the values, there is also no unit for the pre-exp factor, and then I dont mean the orders 0.25 and 1.
The default value for A is 5E8 in the chem.inp, and the 1E8 for the chem.inp_15.
My reference value in Turns book is 5.1E11, but there is no straightforward scaling.
Does anyone know where these values come from ?
Regards

francescomarra · September 2, 2011, 05:59

Dear Rickard,

this is more a general question, concerning not only OpenFOAM. In order to understand the right way for translate coefficients of chemical reactions from one units to another, I would suggest you to write down the full balance equations in the form of ODE: d\rho Y/dt=omega.
By keeping consistent, while changing units, the dimensions of left and right hand side, you will find the rule for unit transformation of the reaction coefficients appearing in omega.

The following report by Andreas Lundstrom:
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/AndreasLundstrom/reactingFoam.pdf
from the excellent web site of Prof. Nilsson, will give you further information and details.

Best regards,

Franco

Rickard.Solsjo · September 2, 2011, 06:07

Hi Franco, thanks for the reply.

I read the tutorial by Anders some time ago, and also checked up on how the chemkinReader for OF works. The problem for me, is when the default values are just taken out of thin air.
The values of 5E8 for C7H16, I can imagine are found from some reference, like Turns, who expresses the value in (gmol/cm3)^(1-m-n)/s and not (mol/cm3....)
What they probably did, was just to take the value 5.11E11 and divided it by 1000 to get mol/cm3, without taking into aspect the effect of the orders m and n.
For ethanol I just did as you said but considered the orders. It should work okay now. The activation energy value however, I dont know where they found it. It seems a bit low to me anyway, for my fuel its around 30 000.

Rickard

September 1, 2011, 03:50	Chemkin units	#1
Rickard.Solsjo Senior Member Rickard Join Date: May 2010 Location: Lund, Skåne, Sweden Posts: 143 Rep Power: 16	Hi Im using chemkin for combustion in dieselFoam and need to change the pre-exponential factor and activation energy for the reactions. However trying to understand the default one is hard since there are no references to the literature for the values, there is also no unit for the pre-exp factor, and then I dont mean the orders 0.25 and 1. The default value for A is 5E8 in the chem.inp, and the 1E8 for the chem.inp_15. My reference value in Turns book is 5.1E11, but there is no straightforward scaling. Does anyone know where these values come from ? Regards

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
units in OpenFOAM	megacrout	OpenFOAM	4	July 19, 2011 09:04
units of geometrical parameters lost	SoobY	ANSYS Meshing & Geometry	2	January 10, 2011 04:58
Mass transfer units in fluent	nite_surfer	FLUENT	0	July 12, 2007 13:51
Question on CHEMKIN INTERFACE	Rui Yang	Phoenics	7	January 30, 2002 05:04
To dimensionalize or not to dimensionalize !!	tolani owosina	Main CFD Forum	8	November 13, 1999 14:45

September 2, 2011, 05:59		#2
francescomarra Member Franco Marra Join Date: Mar 2009 Location: Napoli - Italy Posts: 70 Rep Power: 17	Dear Rickard, this is more a general question, concerning not only OpenFOAM. In order to understand the right way for translate coefficients of chemical reactions from one units to another, I would suggest you to write down the full balance equations in the form of ODE: d\rho Y/dt=omega. By keeping consistent, while changing units, the dimensions of left and right hand side, you will find the rule for unit transformation of the reaction coefficients appearing in omega. The following report by Andreas Lundstrom: http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/AndreasLundstrom/reactingFoam.pdf from the excellent web site of Prof. Nilsson, will give you further information and details. Best regards, Franco

September 2, 2011, 06:07		#3
Rickard.Solsjo Senior Member Rickard Join Date: May 2010 Location: Lund, Skåne, Sweden Posts: 143 Rep Power: 16	Hi Franco, thanks for the reply. I read the tutorial by Anders some time ago, and also checked up on how the chemkinReader for OF works. The problem for me, is when the default values are just taken out of thin air. The values of 5E8 for C7H16, I can imagine are found from some reference, like Turns, who expresses the value in (gmol/cm3)^(1-m-n)/s and not (mol/cm3....) What they probably did, was just to take the value 5.11E11 and divided it by 1000 to get mol/cm3, without taking into aspect the effect of the orders m and n. For ethanol I just did as you said but considered the orders. It should work okay now. The activation energy value however, I dont know where they found it. It seems a bit low to me anyway, for my fuel its around 30 000. Rickard