Tarea 1

UNCP
FIQ
ESCUELA PROFESIONAL: Ingeniería Química
CATEDRA: Ingeniería de las reacciones Químicas II
ALUMNO: Ronaldo José Orihuela Orihuela
Semestre: VIII
Tarea 1
INGENIERÍA DE LA CINÉTICA QUÍMICA – SMITH (EJEMPLO 5.1)

5.1. En un estudio de la producción de aceites secantes por medio de la descomposición de aceite de ricino
acetilado, Grummitt y Fleming correlacionaron los datos de descomposición en base a una reacción de primer
orden representado como:
Aceite de ricino acetilado (l) → CH, COOH (g) + aceite secante (f)
r=kC
Donde r es la velocidad de descomposición, en gramos de ácido acético producida por minuto por mililitro, y C
es la concentración de ácido acético, en gramos por mililitro, equivalente al aceite de ricino acetilado. Los
datos obtenidos en el intervalo de temperatura de 295 a 340 ℃ indicaban una energía de activación de 44
500 cal/mol g, en concordancia con la siguiente expresión para la constante de velocidad específica de la
reacción, k:
−44500
ln k = +35.2
RT
Donde T está en grados Kelvin.
Si un reactor por lotes contiene inicialmente 227 kg de aceite de ricino acetilado a 340 °C (densidad 0.90) y la
operación es adiabática, trace las curvas de conversión (fracción del aceite acetilado que se descompone) y
de temperatura en función del tiempo. Se estima que el efecto calorífico endotérmico de esta reacción es 62
760 J/mol de vapor de ácido acético. El aceite acetilado que se carga al reactor contiene 0.156 kg del
equivalente de ácido acético por kg de aceite, esto es, la descomposición total de 1 kg del aceite produciría
0.156 kg de ácido acético. Suponga que el calor específico de la mezcla reaccionante líquida es constante e
igual a 2.51 x 103 J/ (kg) (K). Suponga también que el vapor de ácido acético que se forma sale del reactor a
la temperatura de la mezcla reaccionante.
SOLUCION: El balance de masa, suponiendo que no hay cambios de volumen durante el curso de la
reacción, es la ecuación:
dx 1 dt C o
= r = ………(A)
dt C 0 dx r
Donde C 0 y r se miden en términos de ácido acético equivalente. La relación concentración-conversión es:
c=C 0 (1−x)
Si reemplazamos r con la ecuación de velocidad expresada en términos de la conversión y la temperatura, la
ecuación (A) se transforma en:
dt Co
=
dx (k /60)C0 (1−x)
dt 60
=f ( T , x )=
dx
[ (
exp 35.2−
44500
RT )]
(1−x )
…………… (B)
Puesto que el reactor opera adiabáticamente, se puede aplicar la ecuación (5-11). Si tomamos 1 kg de aceite
como base, las moles iniciales de ácido acético equivalente son V0CA0 = 0.156/60. Sustituyendo números en la
ecuación (5-11) se obtiene:
( )
0.156
−62760 60 ……………. (C )
T −T f = 3 −3
(x−0)
2.51∗10 10
O bien T =T O−65 x=( 340+ 273 )−65 x ………………(D)
La parte derecha de la Ec.(B) puede ya expresarse en términos de x exclusivamente, sustituyendo la

temperatura por la ecuación (D). De esta forma, la ecuación (B) puede resolverse en forma numérica usando,
por ejemplo, el procedimiento de Runge-Kutta. Las ecuaciones de trabajo para este problema de una variable
dependiente son:
K o =( ∆ X ) f ( X n , T X n
)
Ax
K 1= ( ∆ X ) f ( X n + , T ∆x )
2 (x + 2 )
n
Ax
K 2=( ∆ X ) f ( X n + , T ∆x )
2 (x + 2 )
n
K 3= ( ∆ X ) f ( X n +∆ x , T ( x + ∆ x ))
n
1
t n+1=t n + (K o +2 K 1+ 2 K 2 + K 3 )
6
Donde las ecuaciones (B) y (D) nos dan f(T, x) y n representan el número de incrementos de Ax.Para el primer
incremento se selecciona la Ecua (D) da T=T o=613 K . Entonces de la Ecua (B).
60
f ( X o , T xo )=f ( 0.613 )=
°
=259 s
{ [
exp 35.2−
44500
]}
1.98 ( 613 )
(1−0)
Los valores de k son:
( k , )=0.1 f ( 0.613° )=0.1 ( 259 )=25.9 s
0.1
k 1=0.1 f (0+ ,T 0.1 )
2 0+
2
0 8
¿ 0.1 f ( 0.05 , 610 ) =0.1 ( 327 )=32.7
(
k 2=0.1 f 0+
0.1
2 )
, Ta, , , ,, =32.7 s
k 3=0.1 f ( 0+ 0.1, T 0.1 )=f ( 0.1 ,606.5 0 )=42.7 s
1
t 1=t 0+ ¿
6
¿ 0+33=33 s
Los resultados para los siguientes incrementos se muestran en la tabla 5-1 y la fig. 5-4.Las curvas de
Temperatura y conversión en función del tiempo muestran la necesidad de suministrar energía (en forma de
calor) esta reacción altamente endotérmica cuando se desean conversiones altas. En este caso, sin un
sumistro de energía, la temperatura disminuye tan rápidamente que la reacción de hecho se suspende al
llegar a una conversión del 50 % . Si en lugar de operarse adiabáticamente, se añadiera el “reactor un flujo
constante de energía Q=52700 J/s”, el balance de energía, de acuerdo en las ecuación (5.9) y (5.7) seria.
dT d X A dQ
mcp =−∆ H R ( V 0 C AO ) +
dt dt dt
O numéricamente,
( )
0.156
3 dT 60 dx
227 ( 2.51∗10 ) =−62760 227 +Q ´
dt 10
−3
dt
Integrando desde x=0 , T =T 0 cuando t=0 , esta expresión se transforma en

−6
T =T 0+ 1 ,75∗10 Q ´ ( t−0 )−65 ( x−0 ) (D´)
Por Matlab:
PROGRAMA PARA LA RESOLUCION DEL EJERCICIO

METODO RUNGE-KUTTA
clc ,clear all

disp('---------------------------------------------------')
disp(' METODO Runge-Kutta ')
disp('INGENIERIA DE LAS REACCIONES QUIMICAS II - FIQ')
disp('---------------------------------------------------')
syms x y z
f1=-65;%funcion Numero 1
f2=60/(exp(35.2-44500/(1.987*(y)))*(1-x));%f(T,x)...ECUACION (B)
n0=0;xf=0.6; n=6;x0=0;Y0=613;Z0=0;
h=(xf-n0)/n;
a=zeros(1,n);
for i=1:n
a(i)=n0+h;
n0=a(i);
end
disp('i Conversión(x) k1 c1 k2 c2 k3 c3 k4 c4 T
t')
T=zeros(1,n);t=zeros(1,n);
for i=1:n
x=x0; y=Y0; z=Z0;
k=zeros(1,4);
c=zeros(1,4);
for j=1
k(j)=f1;%K1=(?X)f(Xn+Ax/2,T((xn+?x/2) ))
c(j)=eval(f2);
end
for l=1:2
x=x0+h/2; y=Y0+h*k(l)/2; z=Z0+h*c(l)/2;%K2=(?X)f (Xn+Ax/2,T((xn+?x/2))
k(l+1)=f1;
c(l+1)=eval(f2);
end
for m=3
x=x0+h; y=Y0+h*k(m); z=Z0+h*c(m); %K3=(?X)f(Xn+?x,T((xn+?x)))
k(m+1)=f1;
c(m+1)=eval(f2);
end
y1=Y0+(h/6)*(k(1)+2*k(2)+2*k(3)+k(4));%t(n+1)= tn+1/6(Ko+2K1+2K2+K3)
z1=Z0+(h/6)*(c(1)+2*c(2)+2*c(3)+c(4));
T(i)=y1;
t(i)=z1;
fprintf('%d %10.2f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %10.2f %10.2f\
n',i,a(i),k(1),c(1),k(2),c(2),k(3),c(3),k(4),c(4),y1,z1)
Y0=y1;
Z0=z1;
x0=a(i);
end
subplot(1,2,1),plot(t,a,'m'), title('Conversión en función del tiempo'),
xlabel('Tiempo'),ylabel('Conversión(x)')
subplot(1,2,2),plot(t,T,'c'), title('Temperatura en función del tiempo'),
xlabel('Tiempo'),ylabel('T (K)')
subplot(1,3,1),plotyy(t/60,a,t/60,T),grid on
subplot(1,3,2),plot(t,a,'g'), title('Temperatura en función del tiempo'),
xlabel('Tiempo'),ylabel('T (K)'),grid on
subplot(1,3,3),plot(t,T,'b'), title('Temperatura en función del tiempo'),
xlabel('Tiempo'),ylabel('T (K)'),grid on
disp(' Fin del Programa')
SOLUCION:
---------------------------------------------------
METODO Runge-Kutta
INGENIERIA DE LAS REACCIONES QUIMICAS II - FIQ
---------------------------------------------------
i Conversión(x) k1 c1 k2 c2 k3 c3 k4 c4 T t
1 0.10 -65.000 227.857 -65.000 291.413 -65.000 291.413 -65.000 374.513 606.50 29.47
2 0.20 -65.000 374.513 -65.000 483.825 -65.000 483.825 -65.000 628.565 600.00 78.44
3 0.30 -65.000 628.565 -65.000 821.607 -65.000 821.607 -65.000 1081.133 593.50 161.71
4 0.40 -65.000 1081.133 -65.000 1433.172 -65.000 1433.172 -65.000 1915.553 587.00 307.20
5 0.50 -65.000 1915.553 -65.000 2584.238 -65.000 2584.238 -65.000 3523.771 580.50 570.14
6 0.60 -65.000 3523.771 -65.000 4865.172 -65.000 4865.172 -65.000 6817.921 574.00 1066.84
Fin del Programa

Tarea 1

Cargado por

Copyright:

Formatos disponibles

Tarea 1

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tarea 1

Cargado por

Copyright:

Formatos disponibles

UNCP

INGENIERÍA DE LA CINÉTICA QUÍMICA – SMITH (EJEMPLO 5.1)

Donde T está en grados Kelvin.

Donde C 0 y r se miden en términos de ácido acético equivalente. La relación concentración-conversión es:

O bien T =T O−65 x=( 340+ 273 )−65 x ………………(D)

La parte derecha de la Ec.(B) puede ya expresarse en términos de x exclusivamente, sustituyendo la

( k , )=0.1 f ( 0.613° )=0.1 ( 259 )=25.9 s

k 3=0.1 f ( 0+ 0.1, T 0.1 )=f ( 0.1 ,606.5 0 )=42.7 s

Integrando desde x=0 , T =T 0 cuando t=0 , esta expresión se transforma en

PROGRAMA PARA LA RESOLUCION DEL EJERCICIO

clc ,clear all

También podría gustarte