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

Reaction Mechanism Issue (Therm.dat)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 3, 2016, 09:21
Default Reaction Mechanism Issue (Therm.dat)
  #1
Senior Member
 
Tobias
Join Date: May 2016
Location: Germany
Posts: 264
Rep Power: 10
MFGT is on a distinguished road
Hello,

I recently got the Mechanism from Andrae et al. for Toluene reference fuels, C&F 2008, containing 137 species and 633 reactions. I wanted to compare the results with the Jia TRF Mechanism, included in the multi-component fuel example case of the sector engine.
However, there were a few things preventing me from doing so:

The therm.dat file contained species information, that were missing the exponential format:
CH2CHCHCHC.H2
CH2CHCHCO
-> I added the exponent E within the file

Furthermore CONVERGE told me that there was missing species information for HONO, which was caused by a duplicate HNO species:
-> I deleted the duplicate information. The duplicate also had 5 lines instead of 4.

Also a few species were missing the gas phase indicator G. But this was only a warning.
-> I added the G manually.

Finally, CONVERGE told me that there is no gas species of my parcel species C7H8. This is because toluene is named as C6H5CH3 in that mechanism file set:
-> I replaced C6H5CH3 with C7H8 in both mech.dat and therm.dat.

Since I am no expert in the chemistry, I wanted to ask if someone else has run into similar issues? Was my procedure to repair the files correct?

At least CONVERGE is not reporting any errors right now, but i couldnt test the case as i dont have any capacities at the moment.

Thanks for help or further information.

edit:
The message log in Studio still reports an issue:
ERROR: Therm.dat malformed species C7H14O2H. Error high temp
MFGT is offline   Reply With Quote

Old   August 3, 2016, 13:04
Default
  #2
New Member
 
SauravMitra's Avatar
 
Saurav
Join Date: Dec 2015
Location: Convergent Science, Madison, USA
Posts: 27
Rep Power: 10
SauravMitra is on a distinguished road
Hi MFGT,
I do have to clean therm.dat in some instances the way you have done. Here are some tips to ensure that you have a clean therm.dat

a) Case Setup Module in Studio is a great help, as it fixes many errors and check for consistencies.

b) Try out the Chemistry Module in Studio as well and let it parse the therm.dat. And see if you get any new warnings or messages or not.

c) Finally, use the utility "cleantherm" to have it consistent with mech.dat. To run this utility, go to your working directory and type the following command :

converge cleantherm

It generates a new file, "therm_cleaned.dat" which needs to be renamed as therm.dat. We talk about this in our Manual.



We also talk about the numbers and characters that we expect to be there with each species, that define species name, elemental composition, low temperature, common temperature, high temperature, end of a line etc.


We expect the following in in line number 1 for species C7H14O2H
C7H14O2H 7/23/98 C 7H 15O 2 G 300.000 5000.000

By the error log from Studio it looks like, there is a non-numeric character instead of the high temperature value (bold-faced above). Change it to a numeric value, say 5000.

Have a look at our tutorial cases in the downloads website, and go through a typical therm.dat file.

Anyhow, keep posting your varied experiences and let the Forum know how it goes for you as far this issue is concerned.

Thanks
__________________
Saurav Mitra
Senior Research Engineer
CONVERGECFD
SauravMitra is offline   Reply With Quote

Old   August 4, 2016, 04:40
Default
  #3
Senior Member
 
Tobias
Join Date: May 2016
Location: Germany
Posts: 264
Rep Power: 10
MFGT is on a distinguished road
Quote:
Originally Posted by SauravMitra View Post
Hi MFGT,
I do have to clean therm.dat in some instances the way you have done. Here are some tips to ensure that you have a clean therm.dat

a) Case Setup Module in Studio is a great help, as it fixes many errors and check for consistencies.

b) Try out the Chemistry Module in Studio as well and let it parse the therm.dat. And see if you get any new warnings or messages or not.

c) Finally, use the utility "cleantherm" to have it consistent with mech.dat. To run this utility, go to your working directory and type the following command :

converge cleantherm

It generates a new file, "therm_cleaned.dat" which needs to be renamed as therm.dat. We talk about this in our Manual.



We also talk about the numbers and characters that we expect to be there with each species, that define species name, elemental composition, low temperature, common temperature, high temperature, end of a line etc.


We expect the following in in line number 1 for species C7H14O2H
C7H14O2H 7/23/98 C 7H 15O 2 G 300.000 5000.000

By the error log from Studio it looks like, there is a non-numeric character instead of the high temperature value (bold-faced above). Change it to a numeric value, say 5000.

Have a look at our tutorial cases in the downloads website, and go through a typical therm.dat file.

Anyhow, keep posting your varied experiences and let the Forum know how it goes for you as far this issue is concerned.

Thanks
Hi Saurav,

the species looked like that:
Code:
C7H14O2H DUMMIE         C   7H  15O   2    0G    298.000  5000.000 1392.00     1
 0.25262202E+02 0.34665205E-01-0.11881259E-04 0.18468732E-08-0.10723417E-12    2
-0.29337518E+05-0.99492901E+02 0.15137817E+01 0.88557275E-01-0.59245715E-04    3
 0.21180186E-07-0.32074172E-11-0.20814251E+05 0.28808225E+02                   4
If i replace the first line with yours, its not reporting the error anymore. I dont know if its because of the common temperature of 1392.00 which is added after the high themperature value?
But basically every other species has a common temperature added in that row.

Furthermore, i think there should be only 14 exponential numbers, but some species have 15?
MFGT is offline   Reply With Quote

Old   August 4, 2016, 12:55
Default
  #4
New Member
 
SauravMitra's Avatar
 
Saurav
Join Date: Dec 2015
Location: Convergent Science, Madison, USA
Posts: 27
Rep Power: 10
SauravMitra is on a distinguished road
Hi MFGT,

The spaces in between the numbers matter. Refer to the therm.dat in our tutorial cases. If the numbers/characters are in different columns, then there would be a parsing issue.

Secondly, by default the common temperature for any species will be the one that you specify at the very begining of therm.dat file. If you wanted to specify a different one for a particular species, as you did (1392), just make sure the spaces in between the numbers/characters are correct.

We read 14 exponential numbers, more specifically 4 numbers in line 4.

Thanks
__________________
Saurav Mitra
Senior Research Engineer
CONVERGECFD
SauravMitra is offline   Reply With Quote

Old   August 18, 2016, 05:54
Default
  #5
Senior Member
 
Tobias
Join Date: May 2016
Location: Germany
Posts: 264
Rep Power: 10
MFGT is on a distinguished road
All the hard work for nothing, as it turned out replacing the Jia mechanism with the Andrae mechnaism increased my computational time by a factor of 2.4

At the meoment, i can not simply have these large computational times, so i am forced to stick to Jia mechanism (or others of similar size).
MFGT is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about reaction mechanism in ReactingFoam Dan1788 OpenFOAM Running, Solving & CFD 6 December 13, 2016 06:09
Having problem importing surface reaction mechanism from CHEMKIN Quntin FLUENT 0 March 10, 2016 02:25
Jet-A/air reaction mechanism in FLUENT Jenny_W FLUENT 1 November 29, 2013 02:12
Segmentation fault in running alternateSteadyReactingFoam,why? NewKid OpenFOAM 18 January 20, 2011 17:55
chemical reaction - decompostition La S. Hyuck CFX 1 May 23, 2001 01:07


All times are GMT -4. The time now is 12:40.