CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

write forward, reverse reaction coefficient, reaction rate in ChemFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 15, 2012, 09:07
Default write forward, reverse reaction coefficient, reaction rate in ChemFoam
  #1
New Member
 
Terrence Nguyen
Join Date: Jan 2012
Posts: 11
Rep Power: 5
hismother is on a distinguished road
Hi Foamers,

I am using ChemFoam to check whether my chemical mechanism is good. Now I want to print out the forward, reverse reaction coefficient (kf, kr) as well as forward, reverse reaction rate for each reaction from my reaction mechanism for every chemical timestep. After that I can plot the kf, kr over the time. From the code I know it kf, kr are calculated in ODEChemistryModel.C in each iteration during calculating omega but I can't manage to print them out and save in a output file. If anyone has any idea about it, would be great to enlight me.

Terrence.
hismother is offline   Reply With Quote

Old   November 6, 2012, 02:48
Default
  #2
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 6
Hanzo is on a distinguished road
Hi Terrence,

I'm also digging in the source code of chemFoam & friends to solve my problem
(Solver Validation chemFoam using analytical Solution)

and I wonder if you have managed to print out the reaction coefficients and rates.
Hanzo is offline   Reply With Quote

Old   November 6, 2012, 04:53
Default
  #3
New Member
 
Terrence Nguyen
Join Date: Jan 2012
Posts: 11
Rep Power: 5
hismother is on a distinguished road
just go to ODEchemistrysolver.C. There you can find the place where kf, kr as well as reaction rates are computed. Just use Info command to print them out and they will be printed out in your log file.
hismother is offline   Reply With Quote

Old   November 7, 2012, 04:33
Default
  #4
Member
 
Join Date: Nov 2010
Location: Tokyo / Japan
Posts: 40
Rep Power: 6
Hanzo is on a distinguished road
Quote:
Originally Posted by hismother View Post
just go to ODEchemistrysolver.C. There you can find the place where kf, kr as well as reaction rates are computed. Just use Info command to print them out and they will be printed out in your log file.
thanks for the info.

Actually, I tried to avoid editing ODEchemistryModel.C (I guess you were reffering to this file, right?) because I wanted to avoid defining my own lib and include this in chemFoam. Since there is no other option to get these values, I gave up and recompiled it into a mychemistryModel library.
Took some time but now it's fine.
Hanzo is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Helyx-OS (GUI for SnappyHexMesh elvis OpenFOAM Native Meshers: snappyHexMesh and Others 176 July 10, 2015 16:50
Reaction rate under multi phases kongl1986 Fluent UDF and Scheme Programming 3 September 7, 2012 08:43
mesh file for flow over a circular cylinder Ardalan Main CFD Forum 6 April 17, 2010 23:40
Who can write the binary diffusion coefficient ?? Hsin-Ping FLUENT 1 April 15, 2002 02:22
reaction rate as UDF variable sivakumar FLUENT 4 February 12, 2002 01:31


All times are GMT -4. The time now is 20:03.