Example 7-2 PSSH Applied to Thermal Cracking of Ethane |POLVER05_3
d(C1)/d(t) = -k1*C1-k2*C1*C2-k4*C1*C6 #
d(C2)/d(t) = 2*k1*C1-k2*C1*C2 #
d(C6)/d(t) = k3*C4-k4*C6*C1 #
d(C4)/d(t) = k2*C1*C2-k3*C4+k4*C6*C1-k5*C4^2 #
d(C7)/d(t) = k4*C1*C6 #
d(C3)/d(t) = k2*C1*C2 #
d(C5)/d(t) = k3*C4 #
d(C8)/d(t) = 0.5*k5*C4^2 #
d(CP5)/d(t) = k3*(2*k1/k5)^0.5*CP1^0.5 #
d(CP1)/d(t) = -k1*CP1-2*k1*CP1-(k3*(2*k1/k5)^0.5)*(CP1^0.5) #
k5 = 3980000000 #
T = 1000 #
k1 = 10*exp((87500/1.987)*(1/1250-1/T)) #
k2 = 8450000*exp((13000/1.987)*(1/1250-1/T)) #
k4 = 2530000000*exp((9700/1.987)*(1/1250-1/T)) #
k3 = 3200000*exp((40000/1.987)*(1/1250-1/T)) #
t(0)=0
C1(0)=0.1
C2(0)=0
C6(0)=0
C4(0)=0
C7(0)=0
C3(0)=0
C5(0)=0
C8(0)=0
CP5(0)=0
CP1(0)=0.1
t(f)=12