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 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	`O_2`
3	`CO_2`
4	`O`
5	`H_2O`
6	`OH`
7	`H`

Reaction Notation:

i	Reactions
`R_1`	`CO + O_2 iff CO_2 +O`
`R_2`	`O + H_2O iff OH +OH`
`R_3`	`CO + OH iff CO_2 + H`
`R_4`	`H + O_2 iff OH +O`

Creating matrix for reactant side

`v_(ji)' = [[1,1,0,0,0,0,0],[0,0,0,1,1,0,0],[1,0,0,0,0,1,0],[0,1,0,0,0,0,1]]`

matrix for product side

`v_(ji)'' = [[0,0,1,1,0,0,0],[0,0,0,0,0,2,0],[0,0,1,0,0,0,1],[0,0,0,1,0,1,0]]`

Thus,

`v_(ji) = v_(ji)'' - v_(ji)'`

`v_(ji) = [[-1-,1,1,1,0,0,0],[0,0,0,-1,-1,2,0],[-1,0,1,0,0,-1,1],[0,-1,0,1,0,1,-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

`q_i = k_(fi) prod_(J=1)^N [X_j] ^(v_(ji)' )-k_(fi) prod_(J=1)^N [X_j] ^(v_(ji)'' )`

Now, net reaction rate for `r_1` is

`q_1 = k_(f1) [CO][O_2]-k_(r1) [CO_2][O]`

Similarly net reaction rate

`q_2 = k_(f2) [O][H_2O]-k_(r2) [OH][OH]`

`q_3 = k_(f3) [CO][OH]-k_(r3) [CO_2][H]`

`q_4 = k_(f4) [H][O_2]-k_(r4) [OH][O]`

ODE's for each species

`d[CO]/dt = k_(r1)[CO_2][O]-k_(f1)[CO][O_2] +k_(r3)[CO_2][H]- k_(f3)[CO][OH]`

`d[O_2]/dt = k_(r1)[CO_2][O]-k_(f1)[CO][O_2]+k_(r4)[OH][O]- k_(f4)[H][O_2]`

`d[CO_2]/dt = k_(r1)[CO_2][O]-k_(f1)[CO][O_2]+k_(r3)[CO_2][H]-k_(f3)[CO][OH]`

`d[O]/dt = k_(r1)[CO_2][O]-k_(f1)[CO][O_2]+k_(r2)[OH][OH]- k_(f2)[O][H_2O]+k_(r4)[O][OH]- k_(f4)[H][O_2]`

`d[H_2O]/dt = k_(r2)[OH][OH]- k_(f2)[O][H_2O]`

`d[OH]/dt =k_(r2)[OH][OH]- k_(f2)[O][H_2O] +k_(r3)[CO_2][H]- k_(f3)[CO][OH]`

`d[H]/dt =k_(r3)[CO_2][H]- k_(f3)[CO][OH]+k_(r4)[OH][O]- k_(f4)[H][O_2]`

Product rate can be defined as:

`w_j = sum_(i=1)^L v_(ji)q_i`

Thus the product rate for the reactions can be written as,

For species CO

`w_1 =v_11q_1+v_12q_2+v_13q_3+v_14q_4`

`w_1 =(-1)q_1+(0)q_2+(-1)q_3+(0)q_4`

`w_1 = -k_(f1)[CO][O_2]+k_(r1)[CO_2][O] -k_(f3)[CO][OH]+k_(r3)[CO_2][H]`

ODE can be written as

`d(w_1)/dt =d[CO]/dt`

We have seen here that the equation is same as we extracted above.

Similarly product rate can be write for other species

`w_2 =v_21q_1+v_2q_2+v_23q_3+v_24q_4`

.

.

.

Thus the equations can write in matrix form....

`[[w_1],[w_2],[w_3],[w_4],[w_5],[w_6],[w_7]] = [[-1,-1,1,1,0,0,0],[0,0,0,-1,-1,2,0],[-1,0,1,0,0,-1,1],[0,-1,0,1,0,1,-1]]^T X[[q_1],[q_2],[q_3],[q_4]]`

X represents the each species