CFD Online Logo CFD Online URL
Home > Forums > FLUENT

rate exponent for combustion modeling..

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

LinkBack Thread Tools Display Modes
Old   September 12, 2007, 07:28
Default rate exponent for combustion modeling..
palani kumar
Posts: n/a
hi.. i am working on boundary layer combustion inside a cylindrical chamber..oxidiser (O2 - gas) comes thru the port and pyrolysed fuel (butadiene - gas) comes from the wall..a boundary layer develops on the wall along the axis.. fuel and oxidiser mix here and combust..the reactions tht i use are

C4H6 + 3.5 O2 -----> 4 CO + 3 H2O

CO + 0.5 O2 <------> CO2

can i use laminar flamelet model under "species transport and finite rate chemistry"..??

if yes, how to set the rate exponents..?? i dint find much abt its setting in fluent's user guide..Stephen R Turns in "An intro to combustion" says that they are determined experimentally..Kenneth K Kuo in his book "Principles of combusition" uses examples having stoichiometric coefficients as rate exponents..

if no, then what model can i use..??with "eddy dissipation" model i got unbelievable wrong results..

i let them combust in fuel rich condition (equivalence ratio 2.33) can i use the actual number of moles reacting as rate exponents..?? For ex. i have 1 mole of C4H6 reacting with 1.5 moles of O2..So can i set the stoichiometric coefficients as 1 and 3.5 and the rate exponents as 1 and 1.5 for C4H6 and O2 respectively..???

thanks in advance.. Palani Kumar.
  Reply With Quote

Old   September 12, 2007, 08:52
Default Re: rate exponent for combustion modeling..
Posts: n/a
Hi Palani,

If you are using the pressure based solver the best approach to model hybrid rocket combustion is to create a prePDF look up table for use with the non-premixed combustion model(non adiabatic beta shape PDF). To my knowledge it is the most accurate model for this type of problem and makes things alot easier when setting the model up as you can define the equivalence ratio through the relative mole and species fractions. You are able to generate either single or multiple flamelets although I think a single stretched laminar flamelet is most suitable for modelling the diffusion flame. Although you can use equilibrium chemistry instead which does not require flamelet generation. If your modelling compressible flow I'd go for the finite rate/eddy-dissipation or EDC (use the ISAT tables as direct integration of chemistry is CPU intensive). I'd recommend looking at the following AIAA papers:

Investigation of Metallized and Nonmetallized Hydroxyl Terminated Polybutadiene/Hydrogen Peroxide Hybrid Rockets

Simulation of the Cold-Wall Swirl Driven Combustion Chamber

Of course this is only my view on the subject but the papers mentioned should help you develop an understanding how to model combustion. Hope this helps you and would be interested to hear how your simulations go.

  Reply With Quote

Old   September 13, 2007, 05:56
Default Re: rate exponent for combustion modeling..
palani kumar
Posts: n/a
yes..i am modelling compressible i use density based solver..

even in finite-rate/eddy dissipation model i need to give stoichiometric coefficients and rate exponents for the reactions that i am huv do i decide upon the rate exponents..???

i will look at usage of ISAT of now, i am not yet into PDF models..


  Reply With Quote

Old   September 13, 2007, 11:08
Default Re: rate exponent for combustion modeling..
Posts: n/a
This is something I'm also currently having problems with but I'd suggest typing 'chemical equilibrium' into wikipedia and follow the various links from the article. From what I understand it is not a particularly straight forward process and requires estimates for chamber pressure and combustion temperature. As well as the pre-exponential factor, temperature exponent and activation energy to find the rate constant of the reaction. For the CO + 0.5O2 <---> CO2 reaction I calculated a forward rate constant of around 1.6e+13 for for a chamber pressure of 20.41 bar and assuming a reaction temperature of 3000K. But am unsure if this is correct as I have nothing else to compare this value to or whether the I used the correct units in the calculation. Hopefully someone else might be able to provide more specific details or figures but hope this helps a bit for now.

  Reply With Quote

Old   April 26, 2015, 16:18
sanjeet Limbu
Join Date: Mar 2015
Posts: 91
Rep Power: 4
sanjeetlimbu is on a distinguished road
I am trying to do auto ignition for nheptane/Ar/O2/N2 mixture... with adiabatic wall
Can anyone suggest how to do ..

1. I tried the laminer - species transport: It shows no rise in temp / pressure (I did the cold flow calculation + patch temp 1400K and species reactant), as there is no combustion even after patch temp... just increase in temperature as in attached snapshot

2. I tried the flame let model in Partially premixed but that show error in the unsteady initialisation that the mixture fraction limit is grater than the largest domain fraction even as I entered 0 in it
Attached Images
File Type: jpg model+op_condition.jpg (64.0 KB, 51 views)
File Type: jpg patch.jpg (62.2 KB, 44 views)
File Type: jpg patch _using 1400, 9000_No combustion.jpg (27.7 KB, 41 views)
sanjeetlimbu is offline   Reply With Quote

Old   May 17, 2016, 03:27
New Member
Join Date: Apr 2013
Location: Chennai
Posts: 1
Rep Power: 0
LathaShan is on a distinguished road
Hello Palani,
Did you get the exact rate exponent value? how did you calculate those? please explain, Im also working on that
LathaShan 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
High rate thickener CFD modeling musfiqur Main CFD Forum 0 October 3, 2011 20:38
UDF changing the rate exponent of a reaction Stefan H Fluent UDF and Scheme Programming 0 September 16, 2009 13:20
rate exponent changing Vitali FLUENT 0 July 27, 2007 11:40
Modeling reduced flow rate in 2D Tom Smith FLUENT 2 April 27, 2007 09:04
chemical reaction - decompostition La S. Hyuck CFX 1 May 23, 2001 00:07

All times are GMT -4. The time now is 08:16.