UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA QUÍMICA
INGENIERÍA DE LAS REACCIONES QUÍMICAS II
ESTUDIANTE: CORIÑAUPA ZEVALLOS, CLAUDIA ISABEL.
Se desea producir 120 kg de una solución de ácido acético al 25% en
peso en un reactor Batch a partir de anhídrido acético y agua, de
acuerdo a la reacción:
( C H 3 CO )2 O+ H 2 O → 2C H 3 COOH
A + B → C
los reactantes ingresan simultáneamente a la temperatura de 17°C.
Predecir el tiempo de reacción y la temperatura en función de la
conversión si el proceso se realiza adiabáticamente .
Datos:
cal
∆ H R=−9000
mol de anhidrido acético consumido
1.05 gr
Densidad= ( constante a lo largo de lareacción)
cm3
cal
Cp=0.88 (constante a lo largo de la reacción)
gr ° C
2323
log k =5.233−
T
CB
o
A
(
(−r A )=k C2A ( 1−x A ) C −x A o
o
)
La unidades se expresan en:
L
k [ ¿]
mol ∙min
T [ ¿] K
SOLUCIÓN
1° REALIZAMOS EL BALANCE DE MATERIA
dt C
= A 0 Para (−r A ) :
dx A (−r A )
CB
(−r A )=k C2A ( 1−x A ) C −x A
o
A
( o
o
)
Reemplazando el valor de K:
2323
CB
(−r A )=( 10 )C
5.233−
T 2
Ao
(
( 1−x A ) C −x A
A
o
o
)
Reemplazamos:
dt C A0
=
dx A 2323
C
( 10 5.233−
T )C 2
Ao
B
( 1−x A ) C −x A ( o
)
A o
dt 1
= min
dx A 2323
( 10 5.233−
T )C A0 ( 1−x A ) ( C B −C A x A )o o
Para C A 0 y C B 0
( C H 3 CO )2 O+ H 2 O → 2C H 3 COOH
Ni Ci
N A =N A 0 (1−x¿ ¿ A)¿ C A=C A 0 (1−x A ) R.L. ( C H 3 CO )2 O
N B =N B 0−N A 0 x A C B=C B 0−C A 0 x A
N A0 N B0
C A 0= ;C B 0=
V V
Para N A y N B :
mA (120∗103) g(25 %)
N A= = =294.1176 mol
PM 102 g/mol
mB (120∗103) g(75 %)
NB= = =5000 mol
PM 18 g/mol
Para N A 0 y N B 0:
NA
N A 0=
294.1176
(1−x ¿¿ A)= ¿
(1−x¿¿ A)mol A ¿
294.1176
N B 0=N B + N A 0 x A =5000+
(1−x¿¿ A) x A mol B ¿
Para V:
m 120000 g
V= = =114.2857 L
ρ 1.05∗103 gr /L
Reemplazando:
294.1176
N (1−x ¿¿ A)
C A 0= A 0 = ¿
V 114.2857
294.1176
NB 0 (1−x¿¿ A) x A
C B 0= =5000+ ¿
V 114.2857
dt
Reemplazando en la ecuación de :
dx A
dt 1
= min
dx A 2323
( 10 5.233−
T )C A0 ( 1−x A ) ( C B −C A x A )
o o
dt 1
=
dx A 2323
( 10
5.233−
T ) ¿¿
dt 8.8816∗10−3
= 2323
min
dx A 5.233−
T
10
2° REALIZAMOS EL BALANCE DE ENERGÍA:
dT −ΔH R N A 0
=
dx A mT C P
Reemplazando datos
cal 294.1176
dT
=
9000
molA 1−x A (mol A
)
dx A cal
120000 g 0.88 ( g°C )
dT 25.0668
= K
dx A 1−x A
Entonces resulta que se tiene dos ecuaciones diferenciales ordinarias, para el cual usaremos el método de
Runge-Kutta de cuarto orden:
function Reactor_Batch
clear,clc
disp('-----------------------------------------------')
disp(' REACTOR BATCH ')
disp('-----------------------------------------------')
syms t Te x
%DATOS EXTRAIDOS
x0=0; % COMPOSICIÓN FRACCIONADA INICIAL
Xf=1; % COMPOSICIÓN FRACCIONADA FINAL
t0=0; % TIEMPO INICIAL
Te0=17+273; % TEMPERATURA INICIAL (K)
t01(1)=0;
Te01(1)=17+273;
% ECUACIONES DIFERENCIALES ORDINARIAS
dft=8.8816*10^-3/(10^(5.233-2323/Te));
dfTe=25.0668/(1-x);
h=0.05;
X=[x0:h:xf-h];
n=length(X);
% SE PROCEDE A APLICAR EL MÉTODO RUNGE-KUTTA DE CUARTO ORDEN
disp(' X || t (min) || T (K) ')
disp('-----------------------------------------------')
for i=2:n
k1(i)=double(subs(dft,{x,t,Te},{X(i-1),t0,Te0}));
c1(i)=double(subs(dfTe,{x,t,Te},{X(i-1),t0,Te0}));
k2(i)=double(subs(dft,{x,t,Te},{X(i-1)+h/2,t0+h/2*k1(i),Te0+h/2*c1(i)}));
c2(i)=double(subs(dfTe,{x,t,Te},{X(i-1)+h/2,t0+h/2*k1(i),Te0+h/2*c1(i)}));
k3(i)=double(subs(dft,{x,t,Te},{X(i-1)+h/2,t0+h/2*k2(i),Te0+h/2*c2(i)}));
c3(i)=double(subs(dfTe,{x,t,Te},{X(i-1)+h/2,t0+h/2*k2(i),Te0+h/2*c2(i)}));
k4(i)=double(subs(dft,{x,t,Te},{X(i-1)+h,t0+h*k3(i),Te0+h*c3(i)}));
c4(i)=double(subs(dfTe,{x,t,Te},{X(i-1)+h,t0+h*k3(i),Te0+h*c3(i)}));
t01(i)=t0+h/6*(k1(i)+2*k2(i)+2*k3(i)+k4(i));
Te01(i)=Te0+h/6*(c1(i)+2*c2(i)+2*c3(i)+c4(i));
t0=t01(i);
Te0=Te01(i);
end
for p=1:n
fprintf(' %2.2f | %2.4f | %2.4f\n',X(p),t01(p),Te01(p))
disp('-----------------------------------------------')
end
subplot(2,1,1),plot(t01,X,'-*g'),grid
title('TIEMPO Vs CONVERSIÓN'), xlabel('TIEMPO'),ylabel('CONVERSIÓN')
legend('[CONVERSIÓN]')
subplot(2,1,2),plot(t01,Te01,'-*r'),grid
title('TIEMPO Vs TEMPERATURA'), xlabel('TIEMPO'),ylabel('TEMPERATURA')
legend('[TEMPERATURA]')
end
CORRIDA:
