Challenge

Week 5.1 - Compact Notation Derivation for a simple Mechanism

Compact Notation Derivation for a simple Mechanism

Aim:

Derive the reaction rate ODEs of a simple reaction mechanism using compact notation.

Description:

Chemical kinetics: Chemical kinetics is also known as reaction kinetics. It is the branch of physical chemistry that is concerned with understanding the rates of chemical reactions. it concerns how rates of chemical reactions are measured and predicted, and how reaction-rate data can be used to deduce probable reaction mechanisms. It is to be contrasted with thermodynamics, which deals with the direction in which a process occurs but in itself tells nothing about its rate. Chemical kinetics includes investigations of how experimental conditions influence the speed of a chemical reaction and yield information about the reaction mechanism and transition states, as well as the construction of mathematical models that also can describe the characteristics of a chemical reaction.

Reaction rate:

The reaction rate or rate of reaction is the speed at which a chemical reaction takes place. Reaction rates can vary dramatically. For most reactions, the rate decreases as the reaction proceeds. A reaction's rate can be determined by measuring the changes in concentration over time.

Reaction Mechanism:

A reaction mechanism is the step-by-step sequence of elementary reactions by which overall chemical change occurs.

`CO + O_2 ⇌ CO_2 + O`

`O + H_2O ⇌ OH + OH`

`CO + OH ⇌ CO_2 + H`

`H + O_2 ⇌ OH + O`

Compact notation is a useful tool to describe the properties of each element of a set, rather than writing out all elements of a set.
describe any size mechanism to computer
Reactions:- Reactions are denoted by index letter `i`
  1. Reaction 1 :- `CO + O_2 ⇌ CO_2 + O`
  2. Reaction 2 :- `O + H_2O ⇌ OH + OH`
  3. Reaction 3 :- `CO + OH ⇌ CO_2 + H`
  4. Reaction 4 :- `H + O_2 ⇌ OH + O`
Species participating in reactions(N):- Species are denoted by index letter `j`
  1. `CO`
  2. `O_2`
  3. `CO_2`
  4. `O`
  5. `H_2O`
  6. `OH`
  7. `H`

There is a total of 7 species participating in the reaction.

Reactants stoichiometric coefficient matrix (`nu_(ij)^'`):-

`nu_(ij)^' = [(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)]`

Products stoichiometric coefficient matrix (`nu_(ij)^('')`):

`nu_{ij}^('') = [(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)]`
Net stoichiometric coefficient matrix (`nu_(ij) = nu_(ij)^('')-nu_(ij)^'`):-
`nu_(ij) = [(-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)]`
Net reaction:-
`q_i = k_(fi)prod_(j=1)^N [X_j]^(nu'_(ij)) - k_(ri)prod_(j=1)^N [X_j]^(nu''_(ij))`
`q_1 = k_{f1}[CO]^1[O_2]^1[CO_2]^0[O]^0[H_2O]^0[OH]^0[H]^0 - k_{r1}[CO]^0[O_2]^0[CO_2]^1[O]^1[H_2O]^0[OH]^0[H]^0`
`q_2 = k_{f1}[CO]^0[O_2]^0[CO_2]^0[O]^1[H_2O]^1[OH]^0[H]^0 - k_{r1}[CO]^0[O_2]^0[CO_2]^0[O]^0[H_2O]^0[OH]^2[H]^0`
`q_3 = k_{f1}[CO]^1[O_2]^0[CO_2]^0[O]^0[H_2O]^0[OH]^1[H]^0 - k_{r1}[CO]^0[O_2]^0[CO_2]^1[O]^0[H_2O]^0[OH]^0[H]^1`
`q_4 = k_{f1}[CO]^0[O_2]^1[CO_2]^0[O]^0[H_2O]^0[OH]^0[H]^1 - k_{r1}[CO]^0[O_2]^0[CO_2]^0[O]^1[H_2O]^0[OH]^1[H]^0`
`q_1 = k_{f1}[CO][O_2] - k_{r1}[CO_2][O]`
`q_2 = k_{f2}[O][H_2O] - k_{r2}[OH]^2`
`q_3 = k_{f3}[CO][OH] - k_{r3}[CO_2][H]`
`q_4 = k_{f4}[O_2][H] - k_{r4}[O][OH]`
Reaction rate:-
`(d[X_j])/(dt)=omega_j = sum_(i=1)^L nu_(ij) q_i`
L = number of reactions `L=4`
`[(d[CO]"/"dt),(d[O_2]"/"dt),(d[CO_2]"/"dt),(d[O]"/"dt),(d[H_2O]"/"dt),(d[OH]"/"dt),(d[H]"/"dt)]=((w1),(w2),(w3),(w4),(w5),(w6),(w7)) =[[-1,0,-1,0],[-1,0,0,-1],[1,-1,0,1],[0,2,-1,1],[0,0,1,-1],[1,0,1,0],[0,-1,0,0]]*((q1),(q2),(q3),(q4))`
  1. `omega_1 = -q_1-q_3`
  2. `omega_2 = -q_1-q_4`
  3. `omega_3 = q_1+q_3`
  4. `omega_4 = q_1-q_2+q_4`
  5. `omega_5 = -q_2`
  6. `omega_6 = 2q_2-q_3+q_4`
  7. `omega_7 = q_3-q_4`

ODE'S

  1. `{d[CO]}/{dt} = - k_{f1}[CO][O_2] + k_{r1}[CO_2][O] - k_{f3}[CO][OH] + k_{r3}[CO_2][H]`
  2. `{d[O_2]}/{dt} =- k_{f1}[CO][O_2] + k_{r1}[CO_2][O] - k_{f4}[O_2][H] + k_{r4}[O][OH]`
  3. `{d[CO_2]}/{dt} = k_{f1}[CO][O_2] - k_{r1}[CO_2][O] + k_{f3}[CO][OH] - k_{r3}[CO_2][H]`
  4. `{d[O]}/{dt} = k_{f1}[CO][O_2]-k_{r1}[CO_2][O] - k_{f2}[O][H_2O] + k_{r2}[OH]^2+ k_{f4}[O_2][H] - k_{r4}[O][OH]`
  5. `{d[H_2O]}/{dt} = - k_{f2}[O][H_2O] + k_{r2}[OH]^2`
  6. `{d[OH]}/{dt} = 2k_{f2}[O][H_2O] - 2k_{r2}[OH]^2 - k_{f3}[CO][OH] + k_{r3}[CO_2][H] + k_{f4}[O_2][H] - k_{r4}[O][OH]`
  7. `{d[H]}/{dt} = k_{f3}[CO][OH] - k_{r3}[CO_2][H] - k_{f4}[O_2][H] + k_{r4}[O][OH]`
These ODE's can easily fed to computed and easily solved by using numerical techniques.
Your portfolio, next

Ship work like this
in your own course.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.