CFD Online Logo CFD Online URL
Home > Forums > FLUENT

Modelling reversible chemical reaction.

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

LinkBack Thread Tools Display Modes
Old   August 20, 2010, 09:19
Default Modelling reversible chemical reaction.
New Member
Join Date: Mar 2010
Location: Poland
Posts: 2
Rep Power: 0
kraviec is on a distinguished road

Iím trying to model a liquid-liquid reversible reactions. What is more that is axissymetric problem with swirl. I have some problems, because Iím relatively new in Ansys Fluent.
To start, my reactions are as follows:

1A + 1B <-> 1C + 1D
1C + 1B <-> 1E + 1D

I have all the data for reaction, k constant and activation energy. Iíve calculated pre-exponential factor based on Arrhenius equation. Also I got all required data in materials panel. I know that constant k for reversible reaction is calculated from standard state entropy and enthalpy.
First question is what is rate exponent and how to calculate it? Chemical kinetics is not my best side, so I donít know, if I calculated it properly.

Still it wonít work properly. I use unsteady solver, laminar flow (there is no need to use turbulence model - it is based on known Reynolds number), 2ddp. In species transport, there is used volumetric type of reaction (laminar finite rate).

Another question is about species transport. At the outflow boundary, there is backflow of last (E) compound (I know that the last one is complementary compund). The same situation was, when I was using pressure-outlet type of boundary. Any ideas?

Last but not least question, is about mesh and time step size. Do these problems can result from too small amounts of mesh elements or too large time step, whether to few iteration per time step?

Thanks in advance for all kinds of ideas or suggestions.
kraviec is offline   Reply With Quote

Old   August 20, 2010, 12:08
Allan Walsh
Join Date: Mar 2009
Posts: 56
Rep Power: 10
Allan Walsh is on a distinguished road
Often the rate exponents for the reactants will be the same as the stoichiometric coefficient in the reaction. For example, in the reaction
Ca + 0.5 O2 -> CaO, then
rate = k [Ca] [O2]^0.5.
If you follow the crude computer notation for the equation, then the rate exponent for oxygen is 0.5.
In practice it can be much more difficult. You should research your specific equations to see if there are suggested literature values for a mechanism with rate exponents. Also, the Fluent (Ansys) materials database has quite a few built-in reaction mechanisms and you may get lucky with yours.
In my experience, it is easy to make mistakes on the units for the pre-exponential factors and I recommend doing a hand calculation as a check of the rate coefficients calculated by Fluent.
Allan Walsh is offline   Reply With Quote

Old   August 20, 2010, 13:49
New Member
Join Date: Mar 2010
Location: Poland
Posts: 2
Rep Power: 0
kraviec is on a distinguished road
The reaction is transeterification and accurate data is hard to find, because it happens in spinning disk reactor (we have only data from batch reactor) and that is a new problem for us.

I've done hand calculation for pre-exponential factor. Based on Fluent help, I found that reaction k constant units are m^3/(kmol*s) so pre-exponential factor should have consitnt units. Am I thinking right?

Also may I ask if you could tell me how to change built-in reaction mechanism, because I think I'm using a default one?
kraviec 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
CFX5.7 - chemical reaction problem... H.L. CFX 4 March 25, 2009 11:33
Modelling batch chemical reaction Frans D FLUENT 0 October 30, 2007 23:59
Specifying a chemical reaction hagupta CFX 3 March 22, 2006 02:46
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:56.