Scribd
Cargar
582 vistas7 páginas

Problemas de Matlab

El documento presenta la solución a un examen parcial de termodinámica de procesos químicos. El problema involucra calcular la temperatura a la cual debe operarse un reactor para convertir el 30% de 1-buteno a 1,3-butadieno a una presión de 2 atm. Se presentan los cálculos termodinámicos y cinéticos para determinar que la temperatura requerida es de aproximadamente 770°C.

Cargado por

Irene Hinostroza
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
582 vistas7 páginas

Problemas de Matlab

El documento presenta la solución a un examen parcial de termodinámica de procesos químicos. El problema involucra calcular la temperatura a la cual debe operarse un reactor para convertir el 30% de 1-buteno a 1,3-butadieno a una presión de 2 atm. Se presentan los cálculos termodinámicos y cinéticos para determinar que la temperatura requerida es de aproximadamente 770°C.

Cargado por

Irene Hinostroza
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 7

(Garcia)

Un ingeniero químico una empresa

EXAMEN PARCIAL III


SOLUCIONARIO FILA 1

CÁTEDRA: TERMODINÁMICA DE LOS PROCESOS

QUÍMICOS II

CATEDRÁTICO: Ing. SALVADOR ORE VIDALON

ESTUDIANTES: Hinostroza Meza Ayda Irene

SEMESTRE: V I

HUANCAYO – PERÚ
2019
1. El 1,3 butadieno (C4H6) se puede preparar por la deshidrogenación catalítica en
fase gaseosa del 1-buteno (C4H8):
C4H8 = C4H6 + H2
Con el fin de suprimir las reacciones secundarias, el etilbenceno se diluye con
vapor de agua antes de que pase en el reactor.
Estimar la temperatura a la cual el reactor debe ser operado con el fin de
convertir el 30% del 1-buteno a 1,3-butadieno a una presión de 2 atm a partir de
una alimentación que consiste de 12 en mol de vapor de agua y 1 en mol de 1-
buteno.

NOTA:
𝑚𝑜𝑙𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠 − 𝑚𝑜𝑙𝑒𝑠 𝑓𝑖𝑛𝑎𝑙𝑒𝑠
𝑐𝑜𝑛𝑣𝑒𝑟𝑠𝑖ó𝑛 =
𝑚𝑜𝑙𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠

SOLUCIÓN:
C4H8 = C4H6 + H2
Con avance de una reacción:

𝑑𝑛𝐶4𝐻8 𝑑𝑛𝐶4𝐻6 𝑑𝑛𝐻2


= = = 𝑑ɛ
−1 1 1
Hallando ɛ :
𝑚𝑜𝑙𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠 − 𝑚𝑜𝑙𝑒𝑠 𝑓𝑖𝑛𝑎𝑙𝑒𝑠
𝑐𝑜𝑛𝑣𝑒𝑟𝑠𝑖ó𝑛 =
𝑚𝑜𝑙𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠

1−1+ ɛ
0.3 =
1
ɛ = 0.3
Integrando:

 Para C4H8:
𝐶4𝐻8 ɛ
∫ 𝑑𝑛𝐶4𝐻8 = ∫ −ɛ
𝐶4𝐻8° 0

°
𝑛𝐶4𝐻8 − 𝑛𝐶4𝐻8 = −ɛ

𝑛𝐶6𝐻8 = 1 − ɛ
𝑛𝐶4𝐻8 = 0.7
 Para C4H6:
𝐶4ℎ6 ɛ
∫ 𝑑𝑛𝐶4𝐻6 = ∫ ɛ
𝐶4𝐻6° 0
°
𝑛𝐶4𝐻6 − 𝑛𝐶4𝐻6 = ɛ

𝑛𝐶4𝐻6 = ɛ
𝑛𝐶4𝐻6 = 0.3

 Para H2:
𝐻2 ɛ
∫ 𝑑𝑛𝐻2 = ∫ ɛ
𝐻2° 0
°
𝑛𝐻2 − 𝑛𝐻2 = ɛ

𝑛𝐻2 = ɛ
𝑛𝐻2 = 0.3
Numero de moles de H20:

𝑛𝐻2𝑂 = 12

Numero de moles totales:

𝑛𝑡 = 1 − ɛ + ɛ + ɛ + 12
𝑛𝑡 = 13 + ɛ
𝑛𝑡 = 13.3

Hallando la fracción molar de cada componente:

 Para C4H8:

𝑛𝐶4𝐻8
𝑦𝐶4𝐻8 =
𝑛𝑇
1− ɛ
𝑦𝐶4𝐻8 =
13 + ɛ

𝑦𝐶4𝐻8 = 0.05263

 Para C4H6:
𝑛𝐶4𝐻6
𝑦𝐶4𝐻6 =
𝑛𝑇
ɛ
𝑦𝐶4𝐻6 =
13 + ɛ

𝑦𝐶4𝐻6 = 0.0225

 Para H2:
𝑛𝐻2
𝑦𝐻2 =
𝑛𝑇
ɛ
𝑦𝐻2 =
13 + ɛ

𝑦𝐻2 = 0.0225

Si:

𝐾 = 𝐾𝑃 𝐾𝑦

Donde:

𝐾𝑃 = 𝑃1+1−1

𝐾𝑃 = 21+1−1
𝐾𝑃 = 2
Y:
1 1
𝑦𝐶8𝐻8 ∗ 𝑦𝐻2
𝐾𝑦 = 1
𝑦𝐶8𝐻10
ɛ ɛ
13 + ɛ ∗ 13 + ɛ
𝐾𝑦 =
1− ɛ
13 + ɛ
𝐾𝑦 = 0.00966

Reemplazando en K:

𝐾 = 𝐾𝑃 𝐾𝑦
ɛ ɛ
∗
𝐾 =2∗ 13 + ɛ 13 + ɛ
1− ɛ
13 + ɛ
ɛ ɛ
13 + ɛ ∗ 13 + ɛ
𝐾=2
1− ɛ
13 + ɛ

𝐾 = 0.019334
Utilizando el programa en MATLAB :
EQUILIBRIO EN LAS REACCIONES QUÍMICAS:
clc,clear
format short
disp('EQUILIBRIO EN LAS REACCIONES QUIMICAS')
disp('*****************************')
disp('C6H8 = C4H6 + H2')
disp('*****************************')
%ESCRIBIMOS LAS CONSTANTES DE REACCION
a=1; c=1; d=1;
%ESCRIBIMOS LOS DATOS DEL DATA BANK
C3=[-1.687*10^0 3.419*10^-1 -2.340*10^-4 6.335*10^-8];
C1=[-2.994*10^0 3.532*10^-1 -1.990*10^-4 4.463*10^-8];
C4=[2.714*10 9.274*10^-3 -1.381*10^-5 7.645*10^-9];
%HALLAMOS LAS VARIACIONES
VC1=c*C3(1,1)+d*C4(1,1)-a*C1(1,1);
VC2=c*C3(1,2)+d*C4(1,2)-a*C1(1,2);
VC3=c*C3(1,3)+d*C4(1,3)-a*C1(1,3);
VC4=c*C3(1,4)+d*C4(1,4)-a*C1(1,4);
%CALORES DE FORMACION Y ENERGIA LIBRE DE GIBBS
COMP3=[1.102*10^5 1.508*10^5];
COMP1=[-1.260*10^2 7.134*10^4];
COMP4=[0 0];
%HALLAMOS SUS RESPECTIVAS VARIACIONES
VH298=c*COMP3(1,1)+d*COMP4(1,1)-a*COMP1(1,1);
VG298=c*COMP3(1,2)+d*COMP4(1,2)-a*COMP1(1,2);
%HALLAMOS VARIACION DE HO
VHO=VH298-VC1*298-(VC2*298^2)/2-(VC3*298^3)/3-(VC4*298^4)/4;
%HALLAMOS VHT
syms T
VHT=VHO+VC1*T+(VC2*T^2)/2+(VC3*T^3)/3+(VC4*T^4)/4;
R=8.314;
%ECUACION DE GIBBS HELMHOLTZ
LK=(VHT/R*T^2);
%INTEGRANDO
I=VHO/R*298-VC1*log(298)/R-VC2*298/2*R-(VC3*298^2)/6*R-
(VC4*298^3)/12*R;
lnk=-VHO/R*T+VC1*log(T)/R+VC2*T/2*R+(VC3*T^2)/6*R+(VC4*T^3)/12*R+I;
%HALLAMOS VGT
VGT=lnk*(-R*T);
IR=-VG298/298+VHO/298-VC1*log(298)-VC2*298/2-(VC3*298^2)/6-
(VC4*298^3)/12;
disp('constante K:')
T=[600 650 700 750 800 850 900 920];
J=R.*T;
VHT=(VHO+VC1.*T+(VC2.*T.^2)/2+(VC3.*T.^3)/3+(VC4.*T.^4)/4)/1;
VGT1=(VHO-VC1.*T.*log(T)-(VC2.*T.^2)./2-(VC3.*T.^3)./6-
(VC4.*T.^4)./12-IR.*T)/1;
K=exp(-(VGT1./J))
disp (' K T ')
for e=1:8
disp ([K(e) T(e)])
end
plot(T,K,'-')
title('VARIACIÓN DE LA CONSTANTE DE EQUILIBRIO CON LA TEMPERATURA'),
xlabel('TEMPERATURA, K'),ylabel('CONSTANTE DE EQUILIBRIO')
legend('DATOS HALLADOS')
grid on
RESPUESTA:

EQUILIBRIO EN LAS REACCIONES QUIMICAS

*****************************

C6H8 = C4H6 + H2

****************************

constante K:

K T
0.0001 600.0000
0.0006 650.0000
0.0031 700.0000
0.0118 750.0000
0.0152 760.0000
0.0193 770.0000
0.0245 780.0000
0.0387 800.0000
0.1110 850.0000
0.2837 900.0000
0.4016 920.0000
0.6585 950.0000
1.2162 990.0000
1.3870 999.0000
GRAFICA:

INTERPOLAMOS:

A0.0193---------------B770.0000

C0.019334 -------------DT
F0.0245---------------G780.0000
𝐹−𝐴 𝐺 −𝐵
=
𝐶 −𝐴 𝐷−𝐵
0.0245 − 0.0193 780 − 770
=
0.019334 − 0.0193 𝑇 − 770
780 − 770
𝑇=( ∗ (0.019334 − 0.0193)) + 770
0.0245 − 0.0193

RESULTADO:

T = 770.06

También podría gustarte