Posgrado en Ingeniería Química
División de Ciencias Básicas e Ingeniería
Ingeniería de Reactores Catalíticos
Proyecto 1
Beristain Brito Kevin
Profesor:
Dr. Carlos Omar Castillo Araiza
Resumen:
Se estudia el comportamiento termodinámico de dos reacciones:
1 ⇔
R 1 ¿C 2 H 6+ O2 ❑ C2 H 4 + H 2 O
2
⇔
R2) C 2 H 4 + 3O2 ❑ 2CO 2 +2 H 2 O
Para obtener conversiones y asi optimizar la producción del etileno, variando la temperatura
(1-1000)°C y la presión (1-100) bar, mediante método numéricos puestos en MATlab, para la
obtención de un mejor análisis gráfico, el cual nos permitió entender que en la reacción R1 las
presiones altas llevan a una conversión del etileno, mientras que en la reacción R2, la presión
no afecta, como se hubiera pensado.
Objetivos:
Obtener las Contantes de equilibrio (K) para las reacciones.
1 ⇔
1. C 2 H 6 + O2 ❑ C 2 H 4 + H 2 O
2
⇔
2. C 2 H 4 + 3O2 ❑ 2CO 2 +2 H 2 O
para el estudio termodinámico de la reacción y analizar la cinética.
Metodología:
Tenemos las siguientes ecuaciones
1 ⇔
R 1 ¿C 2 H 6+ O2 ❑ C2 H 4 + H 2 O
2
⇔
R 2 ¿C 2 H 4 +3 O2 ❑ 2 CO2 +2 H 2 O
Primero identificamos cuales son nuestros reactivos limitantes:
Para R1: C 2 H 6 reactivo limitante
Para R2: C 2 H 4 reactivo limitante
Para obtener las constantes de equilibrio K 1 y K 2 de las reacciones R1 y R2 utilizamos la
definición de la constante de equilibrio:
c d
( PC ) ( P D )
K i= a b
( P A) ( PB)
Por lo tanto:
1 1
( PC H ) ( PH O )
K 1= 2
1
4 2
1 /2
( P C H ) ( PO )
2 6 2
2 2
( PCO ) ( PH O )
K 2= 1 3
2 2
( PC H ) ( PO ) 2 4 2
Para las presiones parciales sabemos:
ṅi
Pi= Ptotal
ṅtotal
Donde:
ṅi =flujo de moles del compuesto
ṅtotal =flujo de moles totales
Ptotal=Presion total
Para cada reacción R1:
ṅ C H
PC H = P 2 6
2 6
ṅtotal total
ṅO
PO = 2
P total
2
ṅ total
ṅC H
PC H = P 2 4
2 4
ṅtotal total
ṅH O
PH O= 2
P total
2
ṅ total
Para cada reacción R2:
ṅC H
PC H = P 2 4
2 4
ṅtotal total
ṅO
PO = 2
P total
2
ṅ total
ṅ CO
PCO = P 2
2
ṅtotal total
ṅH O
PH O= 2
P
2
ṅ total total
Como queremos obtener cuanto está saliendo de compuesto, ocupamos la conversión del
compuesto limitante, la definición de conversión es:
ṅi n−ṅ fn
X n=
ṅ in
Donde:
ṅ¿ =Flujo molar inicial
ṅ¿ =Flujo molar final
X n=Convercion
Por lo tanto, para cada compuesto tenemos:
Para R1
ṅ f C H = ṅi C H −( 1−X C
2 6 2 6 2 H6 )
1
ṅ f O = ṅi O − ṅ X
2
2 iC H C2
2 6
2 H6
ṅ f C H =ṅi C H + ṅiC
2 4 2 4 2 H6 XC H 2 6
ṅ f H O =ṅi H O + ṅi C H X C H
2 2 2 6 2 6
Para R2
ṅ f C H =ṅi C H −( 1− X C H
2 4 2 4 2 4
)
ṅ f O = ṅi O −3 ṅi C H X C
2 2 2 4 2 H4
ṅ f CO =ṅi CO + 2˙ni C
2 2 2
H4 XC H 2 4
ṅ f H O =ṅi H O +2 ṅi C
2 2 2 H4 XC H 2 4
Sustituyendo los flujos molares en las presiones parciales y luego en K1 y K2 obtenemos
las siguientes ecuaciones:
Para R1
ṅi C 2 H4 =ṅi H O =0porque al inicio todavía no se producían estos compuestos, por lo tanto:
2
( )( )
1 1
ṅi C H X C H6 ṅi C H X C H6
2 6 2
P total 2 6 2
P total
ṅtotal ṅtotal
K 1=
)( )
1/ 2
1
ṅi O − ṅ X
(
1
ṅi C H6 −( 1− X C H ) 2
2 iC H C H 2 6
2 2 6
Ptotal 2 6
Ptotal
ṅ total ṅtotal
Para R2
ṅi C 2 H4 =ṅi H O =0porque al inicio todavía no se producían estos compuestos por lo tanto:
2
( 2˙ni C H X C
)( )
2 2
H4 2 ṅi C H X C H
2 4 2
Ptotal 2 4 2 4
Ptotal
ṅtotal ṅtotal
K 2=
( )( )
1
ṅ iC H4 −( 1−X C H4 ) ṅiO −3 ṅi C H X C H4
3
Ptotal Ptotal
2 2 2 2 4 2
ṅtotal ṅtotal
Sabemos que K1 y K2 tambien puede expresarse como:
K i=exp ( −△ G(T , P)
RT )
Para obtener:
[ ]
1 1
° ( PC H ) ( P H O )
△ G ( T , P ) =△ G + RTIn
2 4 2
1
1
( PC H ) ( PO )
2 6 2
2
°
△ G =△ H ° reaccion +T 0 △ S ° reaccion
Sustituyendo
R1:
( [ ( )( )
])
1 1
ṅi C H X C H ṅi C H X C H
2 6 2 6
Ptotal 2 6 2 6
Ptotal
ṅtotal ṅtotal
¿
)( )
1 /2
1
ṅi O − ṅ X
(
1
ṅi C H −( 1− X C H ) 2
2 iC H C H 2 6
2 6 2 6
Ptotal Ptotal2 6
° ṅtotal ṅtotal
−△ G
K 1=exp −
RT RT
R2:
[
( ( ( )( )
])
2˙ni C H X C H4
2
2 ṅi C H X C H4
2
2 4 2
Ptotal 2 4 2
Ptotal
ṅtotal ṅ total
¿
)(
1
)
ṅi C H −( 1− X C H ) ṅi O −3 ṅi C H X C H4
3
2 4 2 4
Ptotal 2 2 4 2
Ptotal
° ṅ total ṅ total
−△ G
K 2=exp −
RT RT
Metodología método numérico.
Inicio
Definir los parámetros iniciales: n_C2H4_i, n_O2_i, n_total.
Definir constantes: R, rangos de temperatura T, y presión P_total.
Prealocar Espacio para Resultados
Crear una matriz vacía X_C2H4 para almacenar los resultados de la conversión.
Loop sobre Temperatura y Presión
Para cada temperatura T(i) en el rango definido:
o Para cada presión P(j) en el rango definido:
Convertir la presión a Pascales P.
Definir la Función Objetivo K2_func y K1_func
Usar la ecuación de equilibrio proporcionada y parametrizarla
en términos de la conversión X.
Asegurarse de que la ecuación sea logarítmica y esté
debidamente equilibrada.
Optimización de X usando fminsearch
Llamar a fminsearch para encontrar el valor óptimo de X
que minimice K2_func.
Almacenar el valor optimizado de X en la matriz X_C2H4(i,
j).
Fin del Loop
Resultados:
Para obtener resultados necesitamos valores iniciales por lo tanto para la reacción R1:
mol
ṅi C H6 =50
2
h
mol
ṅi O =50 2
h
mol
ṅtotal =100
h
Se obtiene los siguientes resultados:
Para R1
Para R2:
Para la reacción R2 pondremos como flujo de moles iniciales de ṅi C 2 H4 los mismos que se
producen de la reacción R1.
△ H ° reaccion =−282.53 kJ /mol
△ S ° reaccion =−99.96 kj /mol
Por lo tanto:
°
△ G =△ H ° reaccion +T 0 △ S ° reaccion
°
△ G =−252.71 kJ /mol
Como podemos observar en las gráficas de la reacción R1 al aumentar la presión tiende a favorecer
la dirección de la reacción que produce menos moles de gas. En este caso, como el número de
moles de gas disminuye al pasar de reactivos a productos, una presión más alta favorecería la
formación de C2H4 y H2O.
Mientras que en la reacción R2 al observar las gráficas en este caso, el número de moles de gas es
el mismo en ambos lados de la reacción, por lo que la presión no tendría un efecto significativo en
el equilibrio de esta reacción.
Conclusión.
El análisis termodinámico para una reacción, es de bastante importancia debido a que esto nos
permite manipular las variables de control, para asi obtener los productos deseados, esto no
solo beneficia a la parte económica del proceso, también beneficia al impacto ambiental al
reducir productos no deseados. Al hacer un estudio termodinámico en estas dos reacciones, nos
permitió determinar que para obtener mejores resultados del producto deseado el etileno,
necesitamos subir la presión de la primera reacción, mientras que en la segunda reacción es
mejor enfocarse en la temperatura.
Referencias.
Engineering ToolBox. (n.d.). Standard state - enthalpy of formation, Gibbs free energy, entropy,
molar heat capacity. Retrieved August 16, 2024, from
https://www.engineeringtoolbox.com/standard-state-enthalpy-formation-definition-value-Gibbs-
free-energy-entropy-molar-heat-capacity-d_1978.html
Smith, J. M., Van Ness, H. C., & Abbott, M. M. (2005). Introduction to chemical engineering
thermodynamics (7th ed.). McGraw-Hill.
