CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   CHEMKIN new reaction type (reactingFoam) (https://www.cfd-online.com/Forums/openfoam-programming-development/121401-chemkin-new-reaction-type-reactingfoam.html)

vbchris July 27, 2013 15:53

CHEMKIN new reaction type (reactingFoam)
 
1 Attachment(s)
Hi all,

Simple problem. I am working to add a new reaction type to reactingFoam.

I am using the chemKin utility to add a new reaction. However, I'm not sure how to get chemKinReader to use my new reactionRate. I have modified the reactionRate and chemKinReader files to accept a new reactionRate type (or so I thought). But on run time it doesn't seem to like the syntax of the chemkin case file.

The attached file has the modified chemKinReader and reactionRate Files.

The reaction rate model is a solid core model.

Here's the error.
Code:


Selecting chemistryReader chemkinReader


--> FOAM FATAL ERROR:
unknown third-body specie SLDC on line 21
Valid species are :

7
(
CuO
N2
Cu
CO2
O2
CH4
H2O
)



    From function chemkinReader::lex()
    in file chemistryReaders/chemkinReader/chemkinLexer.L at line 1179.

FOAM exiting

Here is the case chemKin file (copper oxide rxn)

Code:

! GRI-Mech Version 3.0 7/30/99  CHEMKIN-II format
! See README30 file at anonymous FTP site unix.sri.com, directory gri;
! WorldWideWeb home page http://www.me.berkeley.edu/gri_mech/ or
! through http://www.gri.org , under 'Basic  Research',
! for additional information, contacts, and disclaimer
ELEMENTS
H  C  N  O CU
END
SPECIES
CH4    O2      H2O    CO2    N2    Cu      CuO   
END
REACTIONS JOULES/MOLE
 CH4 + 4CuO = CO2 + 2H2O + 4Cu            4.5e-4  0.0  60
    FORD / CH4 0.4/
    FORD / CuO 4/
    SLDC / 1.0 0.00003958 1.5 0.4 0.0015 /
 4Cu + 2O2  = 4CuO                        4.7e-6  0.0  15
    FORD / Cu  2/
    FORD / O2  1/
END


Any comments are welcome.

Thanks,

Chris


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