All Courses
All Courses
Courses by Software
Courses by Semester
Courses by Domain
Tool-focused Courses
Machine learning
POPULAR COURSES
Success Stories
AIM: Derive the compact Notation of simple Mechanism. Reaction Mechanism: In chemical kinetics, we use measurement of the macroscopic properties like,rate of change in the concentration of reactants or products with time, to discover the sequence of events that occur at the molecular level during a reaction. This…
Vishavjeet Singh Yadav
updated on 21 Jan 2022
AIM: Derive the compact Notation of simple Mechanism.
Reaction Mechanism:
In chemical kinetics, we use measurement of the macroscopic properties like,rate of change in the concentration of reactants or products with time, to discover the sequence of events that occur at the molecular level during a reaction. This molecular level describes how individual atoms, ions, molecules interect to form products and all these stepwise changes are collectively called reaction mechanism.
Species Notations:
j |
Species |
1 | CO |
2 | O2 |
3 | CO2 |
4 | O |
5 | H2O |
6 | OH |
7 | H |
Reaction Notation:
i |
Reactions |
R1 | CO+O2⇔CO2+O |
R2 | O+H2O⇔OH+OH |
R3 | CO+OH⇔CO2+H |
R4 | H+O2⇔OH+O |
Creating matrix for reactant side
vji'=⎡⎢ ⎢ ⎢ ⎢⎣1100000000110010000100100001⎤⎥ ⎥ ⎥ ⎥⎦
matrix for product side
vji''=⎡⎢ ⎢ ⎢ ⎢⎣0011000000002000100010001010⎤⎥ ⎥ ⎥ ⎥⎦
Thus,
vji=vji''−vji'
vji=⎡⎢ ⎢ ⎢ ⎢⎣−1−111000000−1−120−10100−110−10101−1⎤⎥ ⎥ ⎥ ⎥⎦
j represents species from 1 to 7 and values in matrix represents the moles of each species in reactions. Also row represents the number of reactions.
Therefor, net reaction rate can be calculated as
qi=kfiN∏J=1[Xj]vji'−kfiN∏J=1[Xj]vji''
Now, net reaction rate for r1 is
q1=kf1[CO][O2]−kr1[CO2][O]
Similarly net reaction rate
q2=kf2[O][H2O]−kr2[OH][OH]
q3=kf3[CO][OH]−kr3[CO2][H]
q4=kf4[H][O2]−kr4[OH][O]
ODE's for each species
dCOdt=kr1[CO2][O]−kf1[CO][O2]+kr3[CO2][H]−kf3[CO][OH]
dO2dt=kr1[CO2][O]−kf1[CO][O2]+kr4[OH][O]−kf4[H][O2]
dCO2dt=kr1[CO2][O]−kf1[CO][O2]+kr3[CO2][H]−kf3[CO][OH]
dOdt=kr1[CO2][O]−kf1[CO][O2]+kr2[OH][OH]−kf2[O][H2O]+kr4[O][OH]−kf4[H][O2]
dH2Odt=kr2[OH][OH]−kf2[O][H2O]
dOHdt=kr2[OH][OH]−kf2[O][H2O]+kr3[CO2][H]−kf3[CO][OH]
dHdt=kr3[CO2][H]−kf3[CO][OH]+kr4[OH][O]−kf4[H][O2]
Product rate can be defined as:
wj=L∑i=1vjiqi
Thus the product rate for the reactions can be written as,
For species CO
w1=v11q1+v12q2+v13q3+v14q4
w1=(−1)q1+(0)q2+(−1)q3+(0)q4
w1=−kf1[CO][O2]+kr1[CO2][O]−kf3[CO][OH]+kr3[CO2][H]
ODE can be written as
dw1dt=dCOdt
We have seen here that the equation is same as we extracted above.
Similarly product rate can be write for other species
w2=v21q1+v2q2+v23q3+v24q4
.
.
.
Thus the equations can write in matrix form....
⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣w1w2w3w4w5w6w7⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦=⎡⎢ ⎢ ⎢ ⎢⎣−1−111000000−1−120−10100−110−10101−1⎤⎥ ⎥ ⎥ ⎥⎦TX⎡⎢ ⎢ ⎢ ⎢⎣q1q2q3q4⎤⎥ ⎥ ⎥ ⎥⎦
X represents the each species
Leave a comment
Thanks for choosing to leave a comment. Please keep in mind that all the comments are moderated as per our comment policy, and your email will not be published for privacy reasons. Please leave a personal & meaningful conversation.
Other comments...
Week 5.2 - Literature review: ODE Stability
AIM: To understand the stability of an ODE and its effect on engineering simulation. Ordinary Differential Equation: It contains derivatives of one or more functions of an Independent variable. We need to determine the function or set of functions that satisfy the equation to solve an ODE. When does a ODE become…
29 Jan 2022 12:30 PM IST
Week 5.1 - Compact Notation Derivation for a simple Mechanism
AIM: Derive the compact Notation of simple Mechanism. Reaction Mechanism: In chemical kinetics, we use measurement of the macroscopic properties like,rate of change in the concentration of reactants or products with time, to discover the sequence of events that occur at the molecular level during a reaction. This…
21 Jan 2022 11:06 AM IST
Week 4.2 - Combustion Efficiency Calculation after Preheating
AIM To calculate the combustion efficiency after Preheating. OBJECTIVES Find the effect of the range of inlet air preheating from 298k to 600k on the adiabatic flame temperature. Find the effect of pre-heating temperature on combustion efficiency. Theory A Recuperator is often used in power engineering…
17 Jan 2022 11:14 AM IST
Week 4.1- Handling Mixtures with Cantera
AIM To work on the quantity class of the Cantera to create various mixtures. Objective Perform calculation of the adiabatic flame temperature of a gas mixture using quantity class. 1) Use the "moles" method/function of the A object(air) and explain how it was calculated. Now consider any hydrocarbon…
20 Apr 2021 07:08 PM IST
Related Courses
0 Hours of Content
Skill-Lync offers industry relevant advanced engineering courses for engineering students by partnering with industry experts.
© 2025 Skill-Lync Inc. All Rights Reserved.