UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA QUÍMICA
ESCUELA ACADÉMICA PROFESIONAL DE INGENIERÍA QUÍMICA
“INGENIERÍA DE LAS REACCIONES QUÍMICAS II 082D”
TEMA: PFR NO ISOTÉRMICOS ADIABÁTICOS
CATEDRÁTICO: Dr. ORE VIDALON, SALVADOR TEODULO
ESTUDIANTE: COSSIO SIERRA KRUGGER
PFR NO ISOTÉRMICOS ADIABÁTICOS
REFERENCIA:
Fogler S. Elementos de Ingeniería de las Reacciones Químicas. Tercera edición.
Ejemplo 8-7 página 462. Producción de anhídrido acético.
_____________________________________________________________________________
La pirólisis en fase vapor de acetona para dar ceteno y metano de
acuerdo a la siguiente reacción:
CH 3 COC H 3 → CH 2 CO+C H 4
es de primer orden respecto a la acetona y la velocidad de reacción
específica se puede expresar como:
34 222
lnk =34.34−
T
Donde k se da en s-1 y T se da en Kelvin. Se desea alimentar 8 000 kg de
acetona por hora a un PFR. El reactor consiste en un banco de 1000 tubos
cédula 40 de 1 pulgada.
La temperatura y la presión de alimentación es de 1035 K y 162 kPa (1.6
atm), respectivamente. Para un reactor que opera adiabáticamente
grafique la conversión y la temperatura en función del volumen del
reactor.
Datos Adicionales:
H oR ( T R ) Cp, J/mol·K
kJ /mol
Aceton -216.67 26.63 + 0.1830T - 45.86 x 10-6T2
a
Ceteno -61.09 20.04 + 0.0945T – 30.95 x 10-6T2
Metano -74.81 13.39 + 0.0770T – 18.71 x 10-6T2
H oR ( T R ) = Calor estándar de formación a la temperatura de referencia
Cp = Capacidad calorífica
SOLUCIÓN
BALANCE DE MATERIA: BALANCE DE ENERGÍA:
dV (−r A )=F Ao d x A dH =Ft CpdT + ∆ H R (−r ¿¿ A )dV ¿
dQ=U o ( T s−T ) d A h=0
d x A (−r A ) dT −∆ H R (−r A )
= =
dV F Ao dV F t Cp
Para una reacción de primer orden:
(−r A )=k C A
34 222
(
k =exp 34.34−
T )
−34 222
k =e 34.34 · e T
−34 222
14 T
k =8.197 ×10 · e
Para un reactor de flujo:
moles de A
F tiempo [ ] moles de A
C A= A [ ¿ ] ¿
v litros litro
tiempo
Para calcular el caudal volumétrico usaremos la ecuación general de los gases
ideales:
PV =nRT
Pv=FT RT
F T RT
v=
P
FA
C A=
F T RT
P
P FA
C A= ( )
RT FT
ESTEQUIOMETRÍA:
Si tenemos una reacción general de la forma:
aA +bB → dD+eE
Tomando A como base de cálculo:
b d e
A+ B → D+ E
a a a
El flujo molar de cada uno de los componentes quedará establecido como:
F A=F Ao−x A F Ao
b
F B=F Bo − x A F Ao
a
d
F D =F Do + x A F Ao
a
e
F E =F Eo+ x A F Ao
a
F I =F Io
____________________________________________________
F T =(F ¿ ¿ Ao+ F Bo + F Do + F Eo + F Io )¿
e d b
(
+ + − −1 x A F Ao
a a a )
F T =F ¿ + δ x A F Ao
F ¿=(F ¿ ¿ Ao+ F Bo + F Do + F Eo + F Io )¿
e d b
(
δ = + − −1
a a a )
Para nuestro problema donde la reacción es de la forma:
CH 3 COC H 3 → CH 2 CO+C H 4
A → D + E
F T =F ¿ + δ x A F Ao
F ¿=(F ¿ ¿ Ao+ F Do + F Eo + F Io )¿
F ¿=(F¿¿ Ao+0+ 0+0)¿
e d b
(
δ = + − −1
a a a )
δ =1+1−1=1
F T =F Ao + x A F Ao=F Ao ( 1+ x A )
tenemos:
F −x F
C A= Ao A Ao
F T RT
P
F Ao ( 1−x A )
C A=
F Ao ( 1+ x A ) RT
P
P 1−x A
C A= ( )(
RT 1+ x A )
Para las condiciones iniciales:
V=0 , xA = 0
P
C Ao=
RTo( )
Para un proceso a presión constante:
NOTA IMPORTANTE:
Si consideramos significativa la variación de la temperatura, se tiene el
siguiente artificio:
P T o 1−x A
C A= ( )( )(
RT T o 1+ x A )
Trasladamos T o al primer paréntesis, y resulta:
P T o 1−x A
C A=
( )( )(
R T o T 1+ x A )
Reemplazando la concentración inicial tenemos:
T o 1−x A
C A=C Ao ( )(
T 1+ x A )
Entonces, la velocidad de reacción será:
T o 1−x A
( )(
(−r A )=k C Ao T 1+ x
A
)
Por tanto, el balance de masa quedará expresado como:
d xA
=
k C Ao ( TT )( 1−x
o
1+ x )
A
dV F Ao
BALANCE DE
ENERGÍA:
dT −∆ H R (−r A )
=
dV F t Cp
T
dT
=
[ o
−(−r A ) ∆ H + ∫ ∆ CpdT
R
298
]
dV F t Cp
Cpi=α i + β i T + γ i T 2
∆ Cp=∆ α + ∆ βT + ∆ γ T 2
T
dT
=
[ o
−(−r A ) ∆ H R + ∫ ( ∆ α + ∆ βT +∆ γ T ) dT
298
2
]
dV F t Cp
Integrando:
∆β 2 ∆γ 3
dT
=
[
−(−r A ) ∆ H oR + ∆ α ( T −T R ) +
2
( T −T 2R ) +
3
( T −T 3R ) ]
dV F t Cp
F t Cp=∑ F i Cpi=F A Cp A + F B CpB + F D CpD + F E CpE
Recordemos que:
F A=F Ao−x A F Ao
b
F B=F Bo − x A F Ao
a
d
F D =F Do + x A F Ao
a
e
F E =F Eo+ x A F Ao
a
F I =F Io
F t Cp=( F Ao −x A F Ao ) Cp A +¿
b
(F Bo )
− x A F Ao CpB +¿
a
d
(F Do + x A F Ao
a )Cp +¿
D
e
(F + x A F Ao
Eo
a E) Cp
_____________________________________________________________
F t Cp=F Ao Cp A + F Bo Cp B + F Do Cp D + F Eo Cp E + ¿
x A F Ao ( ae Cp + da Cp − ba Cp −Cp )
E D B A
e e
Cp E= ( α E + β E T + γ E T 2 )
a a
d d
CpD = ( α D + β D T +γ D T 2 )
a a
−b −d
Cp B= ( α B+ βB T + γ BT 2)
a a
−a −a
a
Cp A =
a
( α A+ β A T + γ A T 2)
_______________________________________
∆ Cp=∆ α + ∆ βT + ∆ γ T 2
F t Cp=∑ F i Cp i=¿ ( F Ao Cp A + F Bo Cp B + F Do Cp D + F Eo Cp E ) +¿ ¿
x A F Ao ( ae Cp + da Cp − ba Cp −Cp )
E D B A
F t Cp=∑ F i Cp i=( ∑ Fio Cpi ) + x A F Ao ∆ Cp
Para el caso del presente problema:
CH 3 COC H 3 → CH 2 CO+C H 4
A → D + E
F t Cp=∑ F i Cp i=¿ ( F Ao Cp A + F Do Cp D + F Eo Cp E ) +¿ ¿
x A F Ao ( ∆ α +∆ βT +∆ γ T 2 )
Si consideramos lo establecido por las condiciones:
F Do=F Eo =F Io=0
F t Cp=∑ F i Cpi=¿ ( F Ao Cp A ) + x A F Ao ( ∆ α + ∆ βT + ∆ γ T 2 ) ¿
F t Cp=∑ F i Cp i=¿ F Ao [ Cp A + x A ( ∆ α + ∆ βT + ∆ γ T 2 ) ] ¿
F t Cp=∑ F i Cp i=¿ F Ao [ (α A + β A T + γ A T 2)+ x A ( ∆ α + ∆ βT +∆ γ T 2 ) ] ¿
Teniendo la ecuación:
∆β 2 ∆γ 3
dT
=
[
−(−r A ) ∆ H oR + ∆ α ( T −T R ) +
2
( T −T 2R ) +
3
( T −T 3R ) ]
dV F t Cp
Reemplazando:
∆β 2 ∆γ 3
dT
=
[
−(−r A ) ∆ H oR + ∆ α ( T −T R ) +
2
( T −T 2R ) +
3
( T −T 3R ) ]
dV F Ao [ (α A + β A T + γ A T 2)+ x A ( ∆ α + ∆ βT +∆ γ T 2) ]
Calor estándar de formación:
ACETONA -216.67 kJ/mol
CETENO -61.09
METANO -74.81
∆ H oR 80.77 kJ/mol
Capacidad calorífica:
ACETONA 26.63 + 0.183 T + -4.59E-05 T^2
CETENO 20.04 + 0.0945 T + -3.10E-05 T^2
METANO 13.39 + 0.077 T + -1.87E-05 T^2
6.8 + -0.0115 T + -3.80E-06 T^2
J
∆ α=6.8
mol·K
∆ β=−0.0115=−11.5· 10−3
∆β
=−5.75 ·10−3
2
∆ γ =−3.8 ·10−6
∆γ
=−1.27 ·10−6
3
RESOLUCION EN MATLAB:
clc,clf,clear all
disp('*********INGENIERIA DE LAS REACCIONES QUIMICAS II*****')
disp('*******PFR NO ISOTERMICO ADIABATICO*****')
%Produccion de anhidrido acetico
%CH3COCH3?CH2CO+CH4
%ejemplo8.7 pagina 462 Fogler S. Elementos rx qcas 3era ed.
% Ecuacion del balance de masa
FA0=38.3;
T0=1035; %Temperatura del enunciado en K
R=8.31;
P=162; %presion del enunciado en kpa
%Ecuacion del balance de Energía
TR=298; %K
DHR0=80770; %unidad convertida a J/mol
%CAPACIDAD CALORIFICA 6.8+(-0.0115)*T+-3.8*(10^(-6))*(T^2)
AA1=6.8; AB1=-0.0115; AY1=-3.8*(10^(-6));
AB2=-5.75*(10^(-3));
AY3=-1.27*(10^(-6));
%capacidad calorifica Acetona
CcAc1=26.63; CcAc2=0.183; CcAc3=-4.59*(10^(-5));
%M=Ft*Cp
dXA=@(V,XA,T) (((exp(34.34-34222/T))*(((P/(R*T0))*1000)*(T0/T)*((1-
XA)/(1+XA))))/FA0); %primera ecuacion EDO
dT=@(V,XA,T) ((-(exp(34.34-34222/T))*(((P/(R*T0))*1000)*(T0/T)*((1-
XA)/(1+XA))))*((DHR0+AA1*(T-TR)+(AB2)*(T^2-(TR)^2)+(AY3)*(T^3-
(TR)^3))))/(FA0*((CcAc1+ CcAc2*T+CcAc3*(T^2))
+XA*(AA1+AB1*T+AY1*(T^2))));
V=[0;0];
XA=[0;0];
T=[1035,0];
i=0;
n=110;
h=(4-0)/n;
disp('----------------------------')
disp(' i XA T(K) V (L)')
disp('----------------------------')
for i=1:15
k1=dXA(V(i),XA(i),T(i));
c1=dT(V(i),XA(i),T(i));
k2=dXA(V(i)+h/2,XA(i)+h*k1/2,T(i)+h*c1/2);
c2=dT(V(i)+h/2,XA(i)+h*k1/2,T(i)+h*c1/2);
k3=dXA(V(i)+h/2,XA(i)+h*k2/2,T(i)+h*c2/2);
c3=dT(V(i)+h/2,XA(i)+h*k2/2,T(i)+h*c2/2);
k4=dXA(V(i)+h,XA(i)+h*k3,T(i)+h*c3);
c4=dT(V(i)+h,XA(i)+h*k3,T(i)+h*c3);
XA(i+1)=double(XA(i)+(1/6)*h*(k1+2*k2+2*k3+k4));
T(i+1)=double(T(i)+(1/6)*h*(c1+2*c2+2*c3+c4));
V(i+1)=V(i)+h;
fprintf('%2d %1.4f %3.2f %3.2f\n',i,XA(i),T(i),V(i));
end
subplot(1,2,1),plot(V,XA,'m')
subplot(1,1,1),plotyy(V,XA,V,T),grid on, title('Concentracion vs
volumen vs Temperatura°K)'),
xlabel('volumen'),ylabel('Concentracion(x)')
CORRIENDO EL PROGRAMA:
----------------------------
i XA T(K) V (L)
----------------------------
1 0.0000 1035.00 0.00
2 0.0435 1014.32 0.04
3 0.0679 1002.56 0.07
4 0.0848 994.35 0.11
5 0.0977 988.06 0.15
6 0.1081 982.98 0.18
7 0.1168 978.71 0.22
8 0.1242 975.04 0.25
9 0.1307 971.82 0.29
10 0.1365 968.96 0.33
11 0.1416 966.38 0.36
12 0.1463 964.04 0.40
13 0.1506 961.90 0.44
14 0.1546 959.92 0.47
15 0.1583 958.09 0.51
Concentracion vs volumen vs Temperatura°K)
0.2 1050
Concentracion(x)
0.1 1000
0 950
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
volumen
