Unformatted text preview: 229 7.1. ISOTHERMAL, ISOCHORIC KINETICS The standard model for chemical reaction, which will be generalized and discussed in
more detail later, induces the following two ordinary diﬀerential equations for the evolution
of O and O2 molar concentrations:
−E 13
RT dρO
= −2 a13 T β13 exp
dt ρO ρO ρM +2 a14 T β14 exp =k13 (T ) = a13 T β13 exp ρO2 ρM , (7.6) =k14 (T ) =r13 dρO2
dt −E 14
RT −E 13
RT =r14 ρO ρO ρM − a14 T β14 exp =k13 (T ) −E 14
RT ρO2 ρM . (7.7) =k14 (T ) =r13 =r14 Here we use the notation ρi as the molar concentration of species i. Also a common usage for
molar concentration is given by square brackets, e.g. ρO2 = [O2 ]. The symbol M represents
an arbitrary third body and is an inert participant in the reaction. We also use the common
notation of a temperaturedependent portion of the reaction rate for reaction j , kj (T ), where
Ej
.
RT
The reaction rates for reactions 13 and 14 are deﬁned as
kj (T ) = aj T βj exp r13 = k13 ρO ρO ρM ,
r14 = k14 ρO2 ρM . (7.8) (7.9)
(7.10) We will give details of how to generalize this form later. The system Eq. (7.67.7) can be
written simply as
dρO
= −2r13 + 2r14 ,
dt
dρO2
= r13 − r14 .
dt
Even more simply, in vector form, Eqs. (7.117.12) can be written as
dρ
= ν · r.
dt (7.11)
(7.12) (7.13) Here we have taken
ρ= ρO
,
ρ O2 (7.14) ν= −2 2
,
1 −1 (7.15) r= r13
.
r14 (7.16) CC BYNCND. 18 November 2011, J. M. Powers. 230 CHAPTER 7. KINETICS IN SOME MORE DETAIL In general, we will have ρ be a column vector of dimension N × 1, ν will be a rectangular
matrix of dimension N × J of rank R, and r will be a column vector of length J × 1. So
Eqs. (7.117.12) take the form
d ρO
dt ρO2 = −2 2
1 −1 r13
.
r14 (7.17) Note here that the rank R of ν is R = L = 1. Let us also deﬁne a stoichiometric matrix φ of
dimension L × N . The component of φ, φli represents the number of element l in species i.
Generally φ will be full rank, which will vary since we can have L < N , L = N , or L > N .
Here we have L < N and φ is of dimension 1 × 2:
φ= 1 2 . (7.18) Element conservation is guaranteed by insisting that ν be constructed such that
φ · ν = 0. (7.19) So we can say that each of the column vectors of ν lies in the right null space of φ.
For our example, we see that Eq. (7.19) holds:
φ·ν = 1 2 · −2 2
1 −1 =0 0. (7.20) The symbol R is the universal gas constant, where
R = 8.31441 J
mole K 107 erg
J = 8.31441 × 107 erg
.
mole K (7.21) Let us take as initial conditions
ρO (t = 0) = ρO , ρO2 (t = 0) = ρO2 . (7.22) Now M represents an arbitrary third body, so here
ρM = ρO2 + ρO . (7.23) Thus, the ordinary diﬀerential equations of the reaction dynamics reduce to
−E 13
ρO ρO ρO2 + ρO
RT
−E 14
+2a14 T β14 exp
ρO2 ρO2 + ρO ,
RT
−E 13
= a13 T β13 exp
ρO ρO ρO2 + ρO
RT
−E 14
ρO2 ρO2 + ρO .
−a14 T β14 exp
RT dρO
= −2a13 T β13 exp
dt dρO2
dt CC BYNCND. 18 November 2011, J. M. Powers. (7.24) (7.25) ...
View
Full
Document
 Fall '11
 ParkSou
 Dynamics

Click to edit the document details