|
[Sponsors] |
![]() |
![]() |
#1 |
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 |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Franco Marra
Join Date: Mar 2009
Location: Napoli - Italy
Posts: 69
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 |
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 |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
units in OpenFOAM | megacrout | OpenFOAM | 4 | July 19, 2011 08:04 |
units of geometrical parameters lost | SoobY | ANSYS Meshing & Geometry | 2 | January 10, 2011 03:58 |
Mass transfer units in fluent | nite_surfer | FLUENT | 0 | July 12, 2007 12:51 |
Question on CHEMKIN INTERFACE | Rui Yang | Phoenics | 7 | January 30, 2002 04:04 |
To dimensionalize or not to dimensionalize !! | tolani owosina | Main CFD Forum | 8 | November 13, 1999 13:45 |