Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

SOLUCIÓN DEL EJEMPLO 5.1 Smith

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 3

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ

CÁTEDRA: INGENIERÍA DE LAS REACCIONES QUIMICAS II


ALUMNO: GASPAR DAVIRAN Daniel G.
E.A.P: Ingeniería Química
SOLUCIONARIO DEL EJEMPLO 5.1
En un estudio de la producción de aceites secantes por medio de la descomposición de
aceite de resino acetilado, Grummitt y Fleming correlacionaron los datos de
descomposición en base a una reacción de primer orden representada 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 “C 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
= + 35.2

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×l03 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.

Solución:
Datos:
Base de cálculo: 1kg de aceite acetilado
Ea: 44500 cal/molg
Q: 62760 J/mol
Ce: 2.51×l03 J/(kg.K).

1. Realizando el balance de materia:


= …………(1)
La relación concentración – conversión es:
= (1 − )
Si reemplazamos r con la ecuación de velocidad expresada en términos de la conversión
y la temperatura, la ecuación (1) obtenemos:
=
(60)(1 − )

= ( , )=[ . ]( )
………..(2)

Puesto que el reactor opera adiabáticamente, aplicará la siguiente ecuación. Si tomamos


1 kg de aceite como base, las moles iniciales de ácido acético equivalente son V,C,O =
0.156/60.
.
− =− . ×
( )( − 0)………(3)

= (340 + 273) − 65 ………………….(4)


2. Reemplazando la ecuación (4) en (1) y resolviendo por el método de Rung
Kutta:
*insertando comandos al programa Matlab:
clc, clear all
disp('================================================')
disp('**************MÉTODO DE RUNGE-KUTTA*************')
disp(' INGENIERIA DE LAS REACCIONES QUIMICAS II' )
disp('================================================')
disp('************************************************')
syms x T t
R=1.987;
dTdx=-65+0*x;
dtdx=60/(exp(35.2-44500/(R*T))*(1-x));
x0=0;
xf=0.4;
Ti=340+273;
To(1)=Ti;
ti=0;
n=4;
h=(xf-x0)/n;
i=0;
fprintf(' i x(i) T(i) t(i)\n')
fprintf('************************************************\n')
fprintf(' %2d %10.6f %10.6f %10.6f\n',i,x0,Ti,ti)
while i<n;
x=x0; T=Ti;t=ti;
k1=eval(dTdx); c1=eval(dtdx);
T=Ti+h*k1/2; t=ti+h*c1/2; x=x0+h/2;
k2=eval(dTdx); c2=eval(dtdx);
T=Ti+h*k2/2; t=ti+h*c2/2; x=x0+h/2;
k3=eval(dTdx); c3=eval(dtdx);
T=Ti+h*k3; t=ti+h*c3; x=x0+h;
k4=eval(dTdx); c4=eval(dtdx);
Ti=Ti+h/6*(k1+2*k2+2*k3+k4);
ti=ti+h/6*(c1+2*c2+2*c3+c4);
x0=x0+h;
i=i+1;
xi(i+1)=x0; yi(i+1)=ti; To(1+i)=Ti;
fprintf(' %2d %10.6f %10.6f
%10.6f\n',i,x0,Ti,ti)
end
fprintf('************************************************\n')
plotyy(yi/60,xi,yi/60,To); title('T° Y CONVERSIÓN EN FUNCIÓN DEL
TIEMPO'),grid on
xlabel('Tiempo, min')
ylabel('Conversión, x')

*Ejecutando el programa

================================================
**************MÉTODO DE RUNGE-KUTTA*************
INGENIERIA DE LAS REACCIONES QUIMICAS II
================================================
************************************************
i x(i) T(i) t(i)
************************************************
0 0.000000 613.000000 0.000000
1 0.100000 606.500000 29.467041
2 0.200000 600.000000 78.439990
3 0.300000 593.500000 161.708750
4 0.400000 587.000000 307.198311
************************************************

T° Y CONVERSIÓN EN FUNCIÓN DEL TIEMPO


0.4 620

0.3 610
Conversión, x

0.2 600

0.1 590

0 580
0 1 2 3 4 5 6
Tiempo, min

También podría gustarte