POLYMATH Results
No Title 03-01-2017, Rev5.1.233
Calculated values of the DEQ variables
Variable initial value minimal value maximal value final value
v 0 0 10 10
fa 0.00344 0.0024962 0.00344 0.0024962
fb 0 0 8.972E-04 8.972E-04
fc 0 0 8.613E-04 8.613E-04
fd 0 0 1.078E-05 1.078E-05
fe 0 0 1.078E-05 1.078E-05
ff 0 0 3.588E-05 3.588E-05
fg 0 0 3.588E-05 3.588E-05
T 800 765.24421 800 765.24421
HIa 1.18E+05 1.18E+05 1.18E+05 1.18E+05
H2a 1.052E+05 1.052E+05 1.052E+05 1.052E+05
H3a -5.39E+04 -5.39E+04 -5.39E+04 -5.39E+04
phi 0.4 0.4 0.4 0.4
KI 0.0458535 0.019636 0.0458535 0.019636
sr 14.5 14.5 14.5 14.5
fi 0.04988 0.04988 0.04988 0.04988
ft 0.05332 0.05332 0.054228 0.054228
Pa 0.1548387 0.1104739 0.1548387 0.1104739
Pb 0 0 0.0397074 0.0397074
Pc 0 0 0.0381196 0.0381196
p 2137 2137 2137 2137
r2b 2.991E-06 5.162E-07 2.991E-06 5.162E-07
rd 2.991E-06 5.162E-07 2.991E-06 5.162E-07
re 2.991E-06 5.162E-07 2.991E-06 5.162E-07
r3t 0 0 4.196E-06 4.151E-06
rf 0 0 4.196E-06 4.151E-06
rg 0 0 4.196E-06 4.151E-06
rIs 2.138E-04 2.479E-05 2.138E-04 2.479E-05
rb 2.138E-04 2.479E-05 2.138E-04 2.479E-05
rc 2.138E-04 2.064E-05 2.138E-04 2.064E-05
ra -2.167E-04 -2.167E-04 -2.946E-05 -2.946E-05
ODE Report (RKF45)
Differential equations as entered by the user
[1] d(fa)/d(v) = ra
[2] d(fb)/d(v) = rb
[3] d(fc)/d(v) = rc
[4] d(fd)/d(v) = rd
[5] d(fe)/d(v) = re
[6] d(ff)/d(v) = rf
[7] d(fg)/d(v) = rg
[8] d(T)/d(v) = -(rIs*HIa+r2b*H2a+r3t*H3a)/(fa*299+fb*273+fc*30+fd*201+fe*90+ff*68+fi*40)
Explicit equations as entered by the user
[1] HIa = 118000
[2] H2a = 105200
[3] H3a = -53900
[4] phi = .4
[5] KI = exp(-17.34-1.302e4/T+5.051*ln(T)+((-2.314e-10*T+1.30e-6)*T-0.004931)*T)
[6] sr = 14.5
[7] fi = sr*.00344
[8] ft = fa+fb+fc+fd+fe+ff+fg+fi
[9] Pa = fa/ft*2.4
[10] Pb = fb/ft*2.4
[11] Pc = fc/ft*2.4
[12] p = 2137
[13] r2b = p*(1-phi)*exp(13.2392-25000/T)*Pa
[14] rd = r2b
[15] re = r2b
[16] r3t = p*(1-phi)*exp(.2961-11000/T)*Pa*Pc
[17] rf = r3t
[18] rg = r3t
[19] rIs = p*(1-phi)*exp(-0.08539-10925/T)*(Pa-Pb*Pc/KI)
[20] rb = rIs
[21] rc = rIs-r3t
[22] ra = -rIs-r2b-r3t
Comments
[1] d(fa)/d(v) = ra
fa0=Kmol/s
[8] d(T)/d(v) = -(rIs*HIa+r2b*H2a+r3t*H3a)/(fa*299+fb*273+fc*30+fd*201+fe*90+ff*68+fi*40)
Energy Balance
[13] KI = exp(-17.34-1.302e4/T+5.051*ln(T)+((-2.314e-10*T+1.30e-6)*T-0.004931)*T)
Kp1
[14] sr = 14.5
vapor
[16] ft = fa+fb+fc+fd+fe+ff+fg+fi
flujo total
[20] r2b = p*(1-phi)*exp(13.2392-25000/T)*Pa
Ley de velocidad para el benceno
[23] r3t = p*(1-phi)*exp(.2961-11000/T)*Pa*Pc
Ley de velocidad para el tolueno
[26] rIs = p*(1-phi)*exp(-0.08539-10925/T)*(Pa-Pb*Pc/KI)
Ley de velocidad para el estireno
[29] ra = -rIs-r2b-r3t
Velocidad de reaccion total
Independent variable
variable name : v
initial value : 0
final value : 10
Precision
Step size guess. h = 0.000001
Truncation error tolerance. eps = 0.000001
General
number of differential equations: 8
number of explicit equations: 22
Data file: C:\Users\ULISES\Documents\ulises\ulises lap\reactores\polymath\tarea simulacion.pol