Web Module: Ozone |POLVER05_0
d(O2)/d(t) =  -r1-r2+r3+r4+r5+r8+r9+r10+2*r11+r12+r13  #PSSH Equations; initial conditions taken from bibliography
d(O)/d(t) =    (2*k1*(O2)+k3*(O3))/(k2*(O2)+k5*(ClO)+k10*(NO2))    #
d(O3)/d(t) =    r2-r3-r4-r9-r11-r12    #
d(Cl)/d(t) =    (k5*(O)*(ClO)+k7*(Cl2O2)+k8*(ClO2)+k14*(BrCl))/(k4*(O3))    #
d(ClO)/d(t) =    ((k7*(Cl2O2)+k8*(ClO2))/(2*k6))^0.5    #
d(Cl2O2)/d(t) =    r6-r7    #
d(ClO2)/d(t) =    r7-r8    #
d(NO)/d(t) =     -r9+r10+r11     #
d(NO2)/d(t) =     r9-r10-r11     #
d(BrCl)/d(t) =    r13-r14    #
d(Br)/d(t) =    k14*(BrCl)/(k12*(O3))    #
d(BrO)/d(t) =    k14*(BrCl)/(k13*(ClO))    #
r1 =     k1*O2     #Rate Laws
r2 =     k2*(O2)*(O)     #
r3 =     k3*(O3)     #
r4 =     k4*(O3)*(Cl)     #
r5 =     k5*(O)*(ClO)     #
r6 =     k6*(ClO)*(ClO)     #
r7 =     k7*(Cl2O2)     #
r8 =     k8*(ClO2)     #
r9 =     k9*(O3)*(NO)     #
r10 =     k10*(NO2)*(O)     #
r11 =     k11*(NO2)*(O3)     #
r12 =    k12*(Br)*(O3)    #
r13 =    k13*(BrO)*(ClO)    #
r14 =    k14*(BrCl)    #
k1 =  1.03E-94  #Rates of reaction (cm3/molecule*s)
k2 =     5.92E-34     #
k3 =     4.38E-26     #
k4 =     1.22E-11     #
k5 =     3.8E-11     #
k6 =     2.16E-32     #
k7 =     8.53E-15     #
k8 =     6.51E-13     #
k9 =     1.78E-14     #
k10 =     1.05E-11     #
k11 =     1E-18     #
k12 =    1.18E-12    #
k13 =    1.02E-12    #
k14 =  1.06E-8  #*Assumed value, pending real value*
t(0)=0
O2(0)=2E17
O(0)=1E-6
O3(0)=5E12
Cl(0)=6.46E4
ClO(0)=7.22E7
Cl2O2(0)=1E-6
ClO2(0)=1E-6
NO(0)=1E9
NO2(0)=1E9
BrCl(0)=1E-6
Br(0)=3.2E5
BrO(0)=6.4E6
t(f)=7E4