5ta Práctica Calificada

Cinética Química y Diseño de

Reactores
A. Ugarte, 19991351H
1. Piden: Entonces las concentraciones serán:
a. Conversión de ácido acético en función del
𝐹𝑂 (𝑤𝐻 )𝑂 𝑡(1 − 𝑥)
tiempo, desde 0 min hasta 120 min. [𝐻] =
𝑀𝐻 𝑉
Datos (1/𝑀𝑂𝐻 )(𝑚𝑂𝐻 )𝑖 − 𝑥 𝐹𝑂 (𝑤𝐻 )𝑂 𝑡(1/𝑀𝐻 )
[𝑂𝐻] =
b. Reactor semicontinuo de tanque agitado a 𝑉
temperatura constante (1/𝑀𝑊 ) 𝐹𝑂 (𝑤𝑊 )𝑂 𝑡 + 𝑥 𝐹𝑂 (𝑤𝐻 )𝑂 𝑡(1/𝑀𝐻 )
[𝑊] =
c. Reacción de esterificación del ácido acético, 𝑉
 H= acido, OH= alcohol etílico
𝑥 𝐹𝑂 (𝑤𝐻 )𝑂 𝑡(1/𝑀𝐻 )
 W=agua, E= éster [𝐸] =
𝑉
d. Condiciones:
i. Masa de alcohol en el reactor al inicio: Reaccion elemental de segundo orden:
 mOH=400Lb
𝑟 = 𝑘1 [𝐻][𝑂𝐻] − 𝑘2 [𝑊][𝐸]𝑀𝐻
ii. Flujo másico de la solución de ácido acético:
 Fo=3.92lb/min El volumen de la reacción varía según:
iii. Fracción en peso de ácido acético:
𝐹𝑂 𝑚𝑜
 (wH)o=0.426 𝑉 = 𝑉𝑜 + 𝜌
𝑡 , 𝑉𝑜 = 𝜌
 (ww)o=1-0.426=0.574
iv. Densidad de la mezcla reaccionante: Luego la ec. De diseño se puede exresar como sigue:
 d=59.8 lb/pie3
𝑑(𝑥𝑡) 𝑀𝐻 𝑟𝑉
e. Cinética =
𝑑𝑡 𝐹𝑂 (𝑤𝐻 )𝑂
i. k1=0.00763 pie3/mol-lb min
ii. k2=0.00262 pie3/mol-lb min Se reemplaza la velocidad de reacción y se tiene una
ecuación diferencial ordinaria de primer orden: y’ = f(x,y)
Asunciones:
Se resuelve por el método de Runge-Kutta de cuarto
f. Densidad constante e igual al agua.
orden: y’ = f(x,y) , y(xo)=yo
Resolución:
El método de RK4 indica que:
La ecuación de diseño será: Yi+1 = yi + h/6(k1+2k2+2k3+k4)
Dónde:
Se define x, a cualquier t, como la relación de la cantidad k1=f(xi,yi)
del reactante A convertido con respecto a la cantidad k2=f(xi+h/2,yi+hk1/2)
total de A añadida en el tiempo dado, entonces, la k3=f(xi+h/2,yi+hk2/2)
concentración de A en el reactor estará relacionada con k4=f(xi+h,yi+hk3)
x como sigue: En este caso:
y=xt, x=t
[A]1=[(mA)i+Fo(wA)ot](1-x)/MAV
 Conversion vs tiempo EM
0.35

En Matlab con el Método de Runge-Kutta de 4°orden: 0.3

%Los valores de entrada sera fo=2.2 y

Conversion del ácido acético

0.25
who=0.38
function [X,T]=prob(fo,who) 0.2
moh=400;
wwo=1-who; 0.15
d=59.8;
k1=0.00763; 0.1
k2=0.00262;
t=1; 0.05
h=5;
T=0:5:120; 0
0 20 40 60 80 100 120
X=0:5:120; tiempo (min)
while T(t)<=115
K1=(k1*fo*who*T(t)*(1-X(t))*(moh/46-
En Matlab con la función ODE45 de 4°orden
X(t)*fo*who*T(t)/60)-
k2*X(t)*who*T(t)*(fo*wwo*T(t)/18+X(t)*fo*w %Los valores de entrada sera fo=2.2 y
ho*T(t)/60))/((moh+fo*T(t))*fo*who/d); who=0.38
xt=X(t)*T(t)+h*K1/2; function [X,T]=Prob1(Fo,wo)
t1=T(t)+h/2; Fun=@(t,xt)(0.00763*Fo*wo*(t-xt)*(400/46-
x1=xt/t1; xt*Fo*wo/60)-0.00262*xt*wo*(Fo*(1-
K2=(k1*fo*who*t1*(1-x1)*(moh/46- wo)*t/18+xt*Fo*wo/60))/((400+Fo*t)*Fo*wo/5
x1*fo*who*t1/60)- 9.8);
k2*x1*who*t1*(fo*wwo*t1/18+x1*fo*who*t1/60 [T, XT]=ode45(Fun,[0,120],0);
))/((moh+fo*t1)*fo*who/d); i=2;
xt=X(t)*T(t)+h*K2/2; X(1)=0;
x1=xt/t1; n=size(T);
K3=(k1*fo*who*t1*(1-x1)*(moh/46- while i~=n(1,1)+1
x1*fo*who*t1/60)- X(i)=XT(i)/T(i);
k2*x1*who*t1*(fo*wwo*t1/18+x1*fo*who*t1/60 i=i+1;
))/((moh+fo*t1)*fo*who/d); end
xt=X(t)*T(t)+h*K3; plot(T,X)
t1=T(t)+h; grid on
x1=xt/t1; title('X vs t ')
K4=(k1*fo*who*t1*(1-x1)*(moh/46- xlabel('tiempo (min)')
x1*fo*who*t1/60)- ylabel('Conversión de acido')
k2*x1*who*t1*(fo*wwo*t1/18+x1*fo*who*t1/60 end
))/((moh+fo*t1)*fo*who/d);
xt=X(t)*T(t)+h*(K1+2*K2+2*K3+K4)/6; X vs t
X(t+1)=xt/t1; 0.35

T(t+1)=t1;
t=t+1; 0.3

end
plot(T,X); 0.25
Conversión de acido

grid on;
title('Conversion vs tiempo EM'); 0.2

xlabel('tiempo (min) ');

0.15
ylabel('Conversion del ácido acético')
end
0.1

0.05

0
0 20 40 60 80 100 120
tiempo (min)

Se observa claramente, que ambas graficas son

bastante similares, esto comprueba que la función
ODE45 utiliza el método de Runge-Kutta.
Se comprueba el programa mediante una solución  Cl2 + C3H6 → CH2CHCH2Cl + HCl
en hoja de cálculo de Excel, donde se obtuvieron los  Cl2 + C3H6 → CH2ClCHClCH3l
siguientes resultados: f. Información general y datos cinéticos:
i. Entalpias de reacción, en
BTU/mol-lb
ti xi 60 0.202639733
1. ΔHr1(273K) = -48000
0 0 65 0.21318983 2. ΔHr1(355K) = -48060
5 0.02395548 70 0.223044234 3. ΔHr2(273K) = -79380
10 0.04630525 75 0.232246803 4. ΔHr2(355K) = -79200
15 0.06716353 80 0.240838233 ii. Capacidades caloríficas, en
20 0.08663521 85 0.24885632 BTU/mol-lb R
25 0.10481689 90 0.256336194 1. Cp(C3H6) = 25.3
30 0.12179759 95 0.263310536 2. Cp(Cl2) = 8.6
3. Cp(HCl) = 7.2
35 0.13765938 100 0.269809772
4. Cp(CH2CHCH2Cl) = 28.0
40 0.15247795 105 0.275862252
5. Cp(CH2ClCHClCH3l) = 30.7
45 0.16632314 110 0.281494413 iii. Ecuaciones de la velocidad,
50 0.1792595 115 0.286730928 velocidades de reacción, en mol-
55 0.19134664 120 0.291594841 lb/h pie3
1. r1=206000e-27.200/RTpC3H6pCl2
0.35 2. r2=11.7e-6.860/RTpC3H6pCl2
Las presiones en atm, T en R.
0.3
Asunciones:
0.25
a. Las capacidades caloríficas se mantienen
0.2 constantes.
0.15
b. Los calores de reacción se suponen constantes,
tomando los valores a 355K:
0.1
Resolución:
0.05
Se define la conversión, en moles de cloro reaccionado
0 por mol-lb de alimento total, como x1 para la reacción
0 20 40 60 80 100 120 1 y x2 para la reacción 2.

Expresamos las presiones parciales en función de las

conversiones mediante el siguiente balance (se
2. Piden: supones de entrada 5 moles es irrelevante):
a. Conversión a cloruro de alilo (producto deseado) n(C3H6) = 4 - 5x1 - 5x2
en función de la longitud del tubo. n(Cl2) = 1-5x1 - 5x2
Datos: n(HCl) = 5x1
a. Reactor de flujo pistón de 2.0 cm de n(CH2CHCH2Cl) = 5x1
diámetro interno. n(CH2lCHClCH3l) =5x2
b. Reactantes a 4 moles de propileno/mol n(totales)= 5(1-x2)
de cloro, a 0.85 mol-lb/h. Las presiones parciales quedaran expresados como:
c. Condiciones: Tempera de entrada 200°C, pC3H6 = (1- 5x1 - 5x2)/(5(1-x2))*P=(0.2- x1 - x2)/(1-x2)*P
presión constante a 2 atm
d. Adicionalmente: Coeficiente de
pCl2 = (4 - 5x1 - 5x2)/(5(1-x2))*P=(0.8- x1 - x2)/(1-x2)*P
Las ecuaciones de velocidades quedaran:
transferencia de calor en el interior es
r1=206000e-27.200/RT(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2*P2
5.0 BTU/h pie2 °F.
r2=11.7e-6.860/RT(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2*P2
e. Las reacciones que se consideran en el
Las ecuaciones de diseño serán:
estudio cinético1:
r1dV = Ftdx1 y r2dV = Ftdx2 Balance de masa
1
Reportado por Groll y HEERANE, y Fairbairn, Cheney y FtCpdT=U(Ts-T)dAh-(r1ΔH1+r2ΔH2)dV Balance de energía
Cherniavsky.
El área y el volumen se puede expresar como: (x2’)*i+1= A2 e-27.200/R(Ti+1)(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2
dV=πD2/4dz Luego:
dAh= πDdz Corrector
Reemplazando en los balances de masa y energía: (x1)i+1=(x1)i +( (x1’)i + (x1’)*i+1 )h/2
r1πD2/4dz = Ftdx1 y r2πD2/4dz = Ftdx2 (x2)i+1=(x2)i +( (x2’)i + (x2’)*i+1 )h/2
FtCpdT=(U(Ts-T)πD-(r1ΔH1+r2ΔH2)πD2/4 )dz zi+1=zi + h
Nuestras variables, conversión y temperatura, estarán
en función de la longitud con las siguientes ecuaciones: A. Usando Matlab
dx1/ dz = r1πD2/(4Ft)
dx2/ dz = r2πD2/(4Ft ) Método de Euler Mejorado:
dT/ dz =(U(Ts-T)πD-(r1ΔH1+r2 ΔH2)πD2/4)/(FtCp ) %Los valores de entrada son a=4, b=1,
t=300, P=1.5, d=2
Aparentemente las derivadas de la conversión no function [T,X1,X2,L]=prob2EM(a,b,t,P,d)
dependen de la temperatura, no hay que olvidar la d=d/(12*2.54);
expresión de las velocidades, que por la ec. De m=b/(a+b);
n=1-m;
Arrhenius esta depende de la temperatura. Además Ft=0.85;
estas dos variables están interrelacionados según estas U=5;
dH1=-48000;
ecuaciones.
dH2=-79200;
Se resolverá primero por el método de Euler Cp=[8.6 25.3 28 30.7 7.2];
Modificado: Ao1=206000;
Ao2=11.7;
y’ = f(x,y), y(xo)=yo Ea1=27200;
El método de EM indica que: Ea2=6860;
yi+1/2 = yi + f(xi,yi)h/2 To=(t+273)*1.8;
A1=Ao1*pi*d^2*P^2/(4*Ft);
xi+1/2 = xi + h/2 A2=Ao2*pi*d^2*P^2/(4*Ft);
yi+1 = yi + f(xi+1/2,yi+1/2)h E1=Ea1/1.98584;
En este caso, tenemos dos variables, y definimos como: E2=Ea2/1.98584;
B=U*pi*d;
y = (x1, x2), x=z C=Ao1*dH1*P^2*pi*d^2/4;
Aplicando: D=Ao2*dH2*P^2*pi*d^2/4;
L=0:4:40;
x1’ = A1 e-27.200/RT(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2, x1o=0
X1=0:4:40;
x2’ = A2 e-6.860/RT(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2, x20=0 X2=0:4:40;
T’= B(Ts-T)prom - (C eE1/RT ΔH1+ DeE2/RT ΔH2) (0.2- x1 - x2) T= To:4:(To+40);
h=4;
(0.8- x1 - x2)/(1-x2)2, To=200 z=1;
El termino (Ts – T)prom, se refiere al promedio de esta while z < 11
diferencia en la entrada y al final en el intervalo de dx1=A1*exp(-E1/T(z))*(m-X1(z)-
X2(z))*(n-X1(z)-X2(z))/(1-X2(z))^2;
paso: dx2=A2*exp(-E2/T(z))*(m-X1(z)-
Para la temperatura solo aplicamos el método de X2(z))*(n-X1(z)-X2(z))/(1-X2(z))^2;
Euler: x1=X1(z)+h*dx1;
x2=X2(z)+h*dx2;
Ti+1=Ti + T’i h F=[(m-X1(z)-X2(z)) (n-X1(z)-X2(z))
En las conversiones si aplicamos el método de EM: X1(z) X2(z) X1(z)].*Ft;
FtCp=sum(F.*Cp);
Predictor:
x=-30;
(x1)*i+1 =(x1)i + (x1’)i h (1) e=1;
(x2)*i+1 =(x2)i + (x2’)i h (2) while e>1E-6
y=(2*To-2*T(z)-h*((B*x-(C*exp(-
Luego la temperatura en zi+1: E1/T(z))+D*exp(-E2/T(z)))*(m-X1(z)-
Ti+1=Ti + T’i h, asegurando que (Ts – T)prom, cumpla con la X2(z))*(n-X1(z)-X2(z))/(1-
ecuación. X2(z))^2)/FtCp))/2;
e=x-y;
Se hallan las derivadas en i+1 con x1 y x2 en *i+1, x=y;
obtenidas en (1) y (2) y la temperatura Ti+1: end
(x1’)*i+1= A1 e-27.200/R(Ti+1)(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2
 dT=(B*x-(C*exp(-E1/T(z))+D*exp(- y’ = f(x,y) , y(xo)=yo
E2/T(z)))*(m-X1(z)-X2(z))*(n-X1(z)-
X2(z))/(1-X2(z))^2)/FtCp;
El método de RK4 indica que:
T(z+1)=T(z)+dT*h; Yi+1 = yi + h/6(k1+2k2+2k3+k4)
dx11=A1*exp(-E1/T(z+1))*(m-x1-x2)*(n- Dónde:
x1-x2)/(1-x2)^2;
dx22=A2*exp(-E2/T(z+1))*(m-x1-x2)*(n- k1=f(xi,yi)
x1-x2)/(1-x2)^2; k2=f(xi+h/2,yi+hk1/2)
X1(z+1)=X1(z)+(dx1+dx11)*h/2; k3=f(xi+h/2,yi+hk2/2)
X2(z+1)=X2(z)+(dx2+dx22)*h/2;
z=z+1; k4=f(xi+h,yi+hk3)
end
X1=X1.*(a+b); En nuestro caso se tratara de dos variables x1 y x2
X2=X2.*(a+b);
plot(L,X1+X2); como:
grid on y=(x1, x2), x=z
hold on
plot(L,X1,'r');
Aplicando:
hold on x1’ = A1 e-27.200/RT(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2, x1o=0
plot(L,X2,'r'); x2’ = A2 e-6.860/RT(0.2- x1 - x2) (0.8- x1 - x2)/(1-x2)2, x20=0
title('Conversion vs Longitud');
xlabel('Longitud (pies) '); T’= B(Ts-T)prom - (C eE1/RT ΔH1+ DeE2/RT ΔH2) (0.2- x1 - x2)
ylabel('Conversion de Cloro en la (0.8- x1 - x2)/(1-x2)2, To=200
alimentación'); El termino (Ts – T)prom, se refiere al promedio de esta
end
diferencia en la entrada y al final en el intervalo de
Z X1 X2 X T paso:
0 0 0 0 1031,40 Para la temperatura solo aplicamos el método de
4 0,0107939 0,0118630 0,022657 1043,60 Euler:
8 0,0229032 0,0238503 0,046753 1053,77 Ti+1=Ti + T’i h
12 0,0361142 0,0358569 0,071971 1062,22
16 0,0502060 0,0477962 0,098002 1069,29
En Matlab con Método de Runge-Kutta de cuarto
20 0,0649676 0,0595983 0,124566 1075,20
24 orden:
0,0801952 0,0712062 0,151401 1080,07
28 0,0956916 0,0825726 0,178264 1083,98
32 %Los valores de entrada son a=4, b=1,
0,1112727 0,0936602 0,204932 1087,02
t=300, P=1.5, d=2
36 0,1267734 0,1044402 0,231213 1089,27 function [T,X1,X2,L]=prob2RK4(a,b,t,P,d)
40 0,1420505 0,1148921 0,256942 1090,81 d=d/(12*2.54);
m=b/(a+b);
n=1-m;
Conversion vs Longitud EM
Ft=0.85;
0.35 U=5;
dH1=-48000;
0.3 dH2=-79200;
Conversion de Cloro en la alimentación

Cp=[8.6 25.3 28 30.7 7.2];

0.25 Ao1=206000;
Ao2=11.7;
0.2 Ea1=27200;
Ea2=6860;
0.15
To=(t+273)*1.8;
A1=Ao1*pi*d^2*P^2/(4*Ft);
0.1
A2=Ao2*pi*d^2*P^2/(4*Ft);
E1=Ea1/1.98584;
0.05
E2=Ea2/1.98584;
B=U*pi*d;
0
0 5 10 15 20 25 30 35 40 C=Ao1*dH1*P^2*pi*d^2/4;
Longitud (pies) D=Ao2*dH2*P^2*pi*d^2/4;
L=0:4:40;
Se resuelve usando el método de Runge-Kutta de 4° X1=0:4:40;
X2=0:4:40;
orden: T= To:4:(To+40);
h=4;
z=1; grid on
while z < 11 hold on
F=[(m-X1(z)-X2(z)) (n-X1(z)-X2(z)) plot(L,X1,'r');
X1(z) X2(z) X1(z)].*Ft; hold on
FtCp=sum(F.*Cp); plot(L,X2,'r');
k11=A1*exp(-E1/T(z))*(m-X1(z)- title('Conversion vs Longitud');
X2(z))*(n-X1(z)-X2(z))/(1-X2(z))^2; xlabel('Longitud (pies) ');
k12=A2*exp(-E2/T(z))*(m-X1(z)- ylabel('Conversion de Cloro en la
X2(z))*(n-X1(z)-X2(z))/(1-X2(z))^2; alimentación');
x11=X1(z)+h*k11/2; end
x21=X2(z)+h*k12/2;
x=-30;
e=1; Z X1 X2 X T
while e>1E-6 0 0 0 0 1031,4
y=(2*To-2*T(z)-(h/2)*((B*x-(C*exp(-
4 0,010794 0,011864 0,022658 1043,6
E1/T(z))+D*exp(-E2/T(z)))*(m-X1(z)-
X2(z))*(n-X1(z)-X2(z))/(1- 8 0,022912 0,023854 0,046766 1053,7
X2(z))^2)/FtCp))/2; 12 0,036138 0,035864 0,072003 1062,2
e=x-y; 16 0,050246 0,047806 0,098052 1069,2
x=y;
end
20 0,065022 0,059610 0,124632 1075,2
dT=(B*x-(C*exp(-E1/T(z))+D*exp(- 24 0,080261 0,071221 0,151481 1080,0
E2/T(z)))*(m-X1(z)-X2(z))*(n-X1(z)- 28 0,095766 0,082587 0,178354 1083,9
X2(z))/(1-X2(z))^2)/FtCp; 32 0,111353 0,093675 0,205029 1087,0
t=T(z)+dT*h/2;
36 0,126858 0,104456 0,231314 1089,2
k21=A1*exp(-E1/t)*(m-x11-x21)*(n-x11-
x21)/(1-x21)^2; 40 0,142136 0,114907 0,257043 1090,8
k22=A2*exp(-E2/t)*(m-x11-x21)*(n-x11-
x21)/(1-x21)^2;
x12=X1(z)+h*k21/2;
x22=X2(z)+h*k22/2; 0.35
Conversion vs Longitud RK-4

k31=A1*exp(-E1/t)*(m-x12-x22)*(n-x12-
x22)/(1-x22)^2; 0.3
k32=A2*exp(-E2/t)*(m-x12-x22)*(n-x12-
Conversion de Cloro en la alimentación

x22)/(1-x22)^2; 0.25
x13=X1(z)+h*k31;
x23=X2(z)+h*k32; 0.2
x=-30;
e=1; 0.15
while e>1E-6
y=(2*To-2*T(z)-h*((B*x-(C*exp(- 0.1
E1/T(z))+D*exp(-E2/T(z)))*(m-X1(z)-
X2(z))*(n-X1(z)-X2(z))/(1- 0.05
X2(z))^2)/FtCp))/2;
e=x-y; 0
0 5 10 15 20 25 30 35 40
x=y; Longitud (pies)
end
dT=(B*x-(C*exp(-E1/T(z))+D*exp(-
E2/T(z)))*(m-X1(z)-X2(z))*(n-X1(z)-
X2(z))/(1-X2(z))^2)/FtCp;
T(z+1)=T(z)+dT*h;
k41=A1*exp(-E1/T(z+1))*(m-x13-x23)*(n-
x13-x23)/(1-x23)^2;
k42=A2*exp(-E2/T(z+1))*(m-x13-x23)*(n-
x13-x23)/(1-x23)^2; B. Se comprueba el programa creado mediante la
solución que se realizó en una hoja de cálculo de Excel,
X1(z+1)=X1(z)+h*(k11+2*k21+2*k31+k41)/6; se obtuvieron las siguientes conversiones en función
X2(z+1)=X2(z)+h*(k12+2*k22+2*k32+k42)/6; de la longitud:
z=z+1;
end
X1=X1.*(a+b);
X2=X2.*(a+b);
plot(L,X1+X2);
 Método de Euler Mejorado Método de Runge-Kutta 4° orden:

z T X1 X2 X z T x1i x2i xti

0 1031.400 0.00000 0.00000 0.00000 0 1031.4 0 0 0
4 1043.012 0.01075 0.01185 0.02260 4 1043.011577 0.01076503 0.01185688 0.0226219
8 1052.944 0.02275 0.02381 0.04657 8 1052.967794 0.02278562 0.02382347 0.04660908
12 1061.403 0.03584 0.03579 0.07163 12 1061.465959 0.03589129 0.03580825 0.07169954
16 1068.553 0.04980 0.04771 0.09752 16 1068.661758 0.04988942 0.04773106 0.09762048
20 1074.527 0.06445 0.05950 0.12395 20 1074.681353 0.0645756 0.05952279 0.12409839
24 1079.434 0.07958 0.07109 0.15067 24 1079.630852 0.07974311 0.071125 0.15086811
28 1083.373 0.09498 0.08245 0.17743 28 1083.603533 0.09519137 0.08248951 0.17768088
32 1086.431 0.11047 0.09353 0.20400 32 1086.685134 0.11327788 0.09397999 0.20725786
36 1088.692 0.12589 0.10431 0.23020 36 1088.899371 0.13444146 0.10552315 0.23996462
40 1090.238 0.14109 0.11476 0.25585 40 1090.260718 0.1591217 0.11702802 0.27614972

0.30000 0.3
0.25000 0.25
0.20000 0.2
X1 X1
0.15000 0.15
X2 X2
0.10000 0.1
X X
0.05000 0.05
0.00000 0
0 10 20 30 40 0 10 20 30 40