CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Chemical reaction in buoyantPimpleFoam

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

LinkBack Thread Tools Display Modes
Old   December 17, 2012, 05:50
Default Chemical reaction in buoyantPimpleFoam
New Member
Martin Lubej
Join Date: Jul 2012
Posts: 5
Rep Power: 7
mlubej is on a distinguished road
Hello everyone,

I am making a custom solver from buoyantPimpleFoam base in which I would like to add (for now) a simple chemical reaction a + b -> c

The dissipation of species (volScalarFields) a and b and formation of c follow the simple relation: di/dt = k(T)*a*b, where k(T) is a constant following the Archenius law k(T) = k0 * e^(-Ea/R/T), where are all constants except T - temperature.

So my problem is using the temperature field to calculate the k(T) in a solver with energy balance (h - hEqn).

I tried defining T in createFields.H as any other volScalarField is defined, but the problem is that the program reads T only in the beginning and it stays the same for the whole simulation.

Thank you for any reply.
mlubej is offline   Reply With Quote

Old   December 18, 2012, 05:00
New Member
Martin Lubej
Join Date: Jul 2012
Posts: 5
Rep Power: 7
mlubej is on a distinguished road
Hello again, I would just like to post the solution to my problem. The running temperature at the h - based (hEqn) solvers can be accessed with the function thermo.T(). So I just wrote T = thermo.T() somewhere in the loop, and now all functions of T work. Example:

fvScalarMatrix aEqn
+ fvm::div(phi, Ca)
- fvm::laplacian(Da, Ca)
+ A*pow(2.7183,(-Ea/R/T))*Ca*Cb*rho


This is the equation for transport of specie a with dissipation.

Best regards
mlubej is offline   Reply With Quote


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
Chemical reaction using Eddy Dissipation Model in FLUENT hyperflow FLUENT 5 June 11, 2012 04:54
Segmentation fault in running alternateSteadyReactingFoam,why? NewKid OpenFOAM 18 January 20, 2011 17:55
CFX5.7 - chemical reaction problem... H.L. CFX 4 March 25, 2009 11:33
chemical reaction - decompostition La S. Hyuck CFX 1 May 23, 2001 00:07
chemical reaction - howto Dan Sorensen Siemens 1 February 21, 2000 10:20

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