Diseño de Procesos Químicos Capítulo 3

ESTRUCTURAS DE RECIRCULACIÓN DE DIAGRAMAS DE FLUJO

DE PROCESOS
Si el potencial económico de la Estructura de Entradas-Salidas de un proceso químico
simple es positivo, entonces se puede avanzar a la etapa de diseño de la Estructura de
Recirculación. En este nivel se consideran explícitamente los flujos de las corrientes internas
del proceso, el diseño del reactor químico y del compresor, si es el caso, de cada corriente
gaseosa de recirculación. Se siguen suponiendo separaciones perfectas y el sistema de
separación se representa por un solo bloque de procesamiento, al que entra la corriente de
salida del reactor y salen diversas corrientes acordes a los códigos de destino asignados a los
compuestos.
En relación al reactor químico se deben tomar varias decisiones de diseño en esta
etapa, en particular lo referente a su operación adiabática o isotérmica, así como el porcentaje
del exceso de uno de los reactivos y el tipo de reactor. También debe diseñarse en esta etapa
el compresor asociado a cada una de las recirculaciones gaseosas, debido a que éste es uno
de los equipos más caros de un proceso químico.
La evaluación económica de la Estructura de Recirculación de un proceso adiciona a
los egresos los costos totales del reactor químico y, si es el caso, de los compresores de las
corrientes de recirculación. Para ello es necesario diseñar el reactor químico y el compresor
de una recirculación gaseosa tomando en consideración los flujos molares o másicos de las
corrientes internas de un proceso, así como las condiciones de operación de tales unidades
de proceso.
En estas notas se muestra cómo desarrollar y calcular el potencial económico de la
Estructura de Recirculación de la Opción o Destino 3 del proceso HDA visualizada en la
Figura 2.9; de igual manera, se presenta la Estructura de Recirculación del Proceso de
Producción de Ciclohexano a partir de la hidrogenación de Benceno, con el propósito de
ilustrar cómo plantear y resolver los balances de materia de sistemas de recirculación con
reacciones químicas gobernadas por el equilibrio termodinámico.

3.1 BALANCES DE MATERIA DE ESTRUCTURAS DE RECIRCULACIÓN

La formulación de los balances de materia de la Estructura de Recirculación de un
proceso se realiza en esta sección para el proceso HDA mostrado en la Figura 3.1. Siempre
se inicia con el balance de materia del reactivo limitante, que en este proceso es el Tolueno.
El balance del Tolueno en el mezclador previo al reactor se expresa como
FT  FFT  RL (1)
donde la corriente de recirculación en fase líquida está constituida totalmente por el Tolueno
que no reacciona en el reactor, debido a que se supone que todos los compuestos valiosos
experimentan una separación perfecta y que el reactivo limitante se recircula completamente:
RL  FT 1  x  (2)

Combinando estas dos ecuaciones se obtiene el Tolueno alimentado al reactor

FFT
FT  (3)
x

Dr. Medardo Serna González – FIQ - UMSNH 1

Diseño de Procesos Químicos Capítulo 3

Figura 3.1 – Estructura de Recirculación de la Opción 3 de la Figura 2.9 del Proceso HDA

Una vez que se ha obtenido el flujo del reactivo limitante a la entrada del reactor se
procede a realizar el balance de los demás compuestos, teniendo presente que el exceso de
un reactivo, MR, se puede expresar como la relación del flujo molar del reactivo en exceso
entre el flujo molar del reactivo limitante a la entrada del reactor. Por lo tanto, el flujo del
Hidrógeno a la entrada del reactor es dado por las siguientes expresiones
FH ,in  yFH FG  yPH RG (4)

FH ,in   MR   FT ,in  (5)

Combinando las Ecs. (3) a (5) se llega a

1  MR PB 
RG    y FH FG  (6)
y PH  x S 
Sustituyendo la Ec. (43b) del capítulo 2 de las notas en la Ec. (4) se obtiene
PB  MR y FH 1  1  y PH 1  S  2
RG     (7)
y PH S  x y FH  y PH 
El flujo molar del metano que se alimenta al reactor es dado por
FM ,in  1  yFH  FG  1  yPH  RG (8)

Los flujos molares de los compuestos a la salida del reactor se calculan de las
ecuaciones de balance de materia de los compuestos en esta unidad de proceso, esto es:
Hidrógeno: FH ,out  FH ,in   R1   R 2
(9)
Metano: FM ,out  FM ,in   R1 (10)

Dr. Medardo Serna González – FIQ - UMSNH 2

Diseño de Procesos Químicos Capítulo 3

Benceno: FB ,out  PB   R1  2 R 2 (11)

Tolueno: FT ,out  FT ,in   R1 (12)

Difenilo: FD ,out   R 2 (13)

Una vez especificado un valor para el exceso del Hidrógeno a la entrada del reactor,
MR, es posible calcular los flujos molares de los compuestos en las diversas corrientes
internas del proceso en las que están presentes para valores dados de PB, yPH, yFH y x. Se infiere
correctamente, por ende, que el modelo matemático que incluye los balances globales de
materia y los balances de materia de la estructura de recirculación de esta Estructura de
Recirculación del proceso HDA tiene 5 grados de libertad.

3.2 RELACIÓN ENTRE yPH Y OTRAS VARIABLES INDEPENDIENTES

Los balances globales de materia y los balances de materia de la Estructura de
Recirculación de la opción 3 del proceso HDA se han formulado con base en la fracción
molar del Hidrógeno en la corriente de purga, yPH. Otra opción, que es la que comúnmente se
emplea en simulación de procesos, modela el divisor mostrado explícitamente en la Figura
3.2 en términos de la fracción de división  DIV para representar la división del flujo de la
corriente gaseosa proveniente del sistema de separación en dos corrientes con la misma
composición: la purga y la recirculación gaseosa.
Una de las ventajas de la formulación basada en la fracción de división es que con
certeza se sabe que los valores factibles de este parámetro se encuentran en la región

0   DIV  1.0 (14)

independientemente del número de compuestos que estén presentes en la purga y de la

participación de tales compuestos como reactivos y/o productos en el subsistema
reaccionante del proceso. Su desventaja es que usualmente requiere cálculos iterativos para
obtener los flujos molares de las corrientes internas del proceso.
En contraposición, la formulación basada en la fracción molar del Hidrógeno en la
purga simplifica los cálculos de los balances a consta de no saber a priori el rango preciso de
los valores factibles de yPH. Por lo tanto, se debe generar esta información antes de simular la
Estructura de Recirculación del proceso HDA para evitar cálculos innecesarios y/o errores
de cálculo.
Considerando la restricción impuesta sobre la sumatoria de las fracciones molares de
los compuestos en la purga, que toma la forma:

yPH  y PM  1 (15)

se determina correctamente que para cualquier valor factible de  DIV la fracción molar de cada
uno de los compuestos en la purga siempre será menor que 1.0 y mayor que 0. Sin embargo,
este análisis es insuficiente para definir exactamente la región de los valores factibles de
yPH y, por consiguiente, de yPM en el rango general [0 1], ya que no toma en cuenta los
siguientes aspectos:

Dr. Medardo Serna González – FIQ - UMSNH 3

 Diseño de Procesos Químicos Capítulo 3

 El divisor de flujo solamente es una parte del proceso HDA, que es un sistema
constituido por varias unidades interconectadas con interacciones entre ellas. En otras
palabras, el funcionamiento del divisor depende de los parámetros y valores de las
variables independientes de las demás unidades de proceso, así como los
funcionamientos de tales unidades dependen de la fracción de divisor.
 En la corriente de purga hay dos compuestos: el Hidrógeno, que es uno de los
reactivos en la reacción principal, y el Metano que, a la vez, es una impureza y un
subproducto. El papel que desempeñan en el reactor estos compuestos es también un
factor a considerar, lo que se hace al formular sus correspondientes balances de
materia.
Por lo tanto, con el propósito de conocer con exactitud la región de valores factibles
de yPH, esta variable se debe expresar explícitamente como una función de la fracción del
divisor y de las variables independientes de la Estructura de Recirculación del proceso
HDA mostrada en la Figura 3.2.
 y P F
y R
PH G
DIV PH G H , out

DIV yPH PG
FH ,in
  MR  F
T ,in F
H ,out
FH ,in
 
R1 R2

yFH FG
SISTEMA DE PB
REACTOR
FFT SEPARACIÓN PD

Balances del Hidrógeno

Figura 3.2 – Fracción de división de la Opción 3 de la Figura 2.9 del Proceso HDA

Como se muestra en la Figura 3.2, la fracción de división en términos del Hidrógeno

purgado y suponiendo separaciones perfectas se expresa como

yPH PG
 DIV  (16)
FH ,out

Los balances globales de Hidrógeno y Metano son dados por las Ecs. (17) y (18),
respectivamente,

y PH PG  yFH FG   R1   R 2 (17)

(1  y PH ) PG  (1  yFH ) FG   R1 (18)

Combinando las Ecs. (16) y (17) se llega a

FG   DIV FH ,out   R1   R 2  y FH (19)

Dr. Medardo Serna González – FIQ - UMSNH 4

Diseño de Procesos Químicos Capítulo 3

Sustituyendo esta ecuación en la Ec. (18) se obtiene

(1  yFH )
(1  yPH ) PG 
yFH
 DIV FH ,out   R1   R 2    R1 (20)

El flujo molar total de la corriente de purga es igual a la suma los flujos molares del
Hidrógeno y del Metano purgados, esto es,

PG  yPH PG  (1  yPH ) PG   DIV FH ,out  (1  yPH ) PG (21)

Combinando las Ecs. (20) y (21) se genera la relación

(1  yFH )
PG   DIV FH ,out 
yFH
 DIV FH ,out   R1   R 2    R1 (22)

Teniendo presente que la fracción molar del Hidrógeno yFH en la corriente gaseosa de
alimentación y la razón molar MR son variables independientes conocidas, y que ambos
avances de reacción y el flujo molar del Hidrógeno a la salida del reactor dependen solamente
de la conversión por paso del Tolueno, x, se concluye que el flujo molar total de la corriente
de purga dado por la Ec. (20) depende solamente de la fracción de divisor,  DIV , y de x.
Combinando las Ecs. (16) y (22) se llega a la relación
1
 (1  yFH ) 
yPH   DIV FH ,out   DIV FH ,out    DIV FH ,out   R1   R 2    R1  (23)
 yFH 

que muestra cómo la fracción molar del Hidrógeno en la corriente de purga también depende
únicamente de la fracción de divisor,  DIV , y de la conversión, x.
La fracción del divisor puede tomar valores con significado físico solamente en el
rango dado por la Ec. (14). Un valor de  DIV igual a 1.0 implica, por definición, que no existe
la corriente gaseosa de recirculación. En este límite superior, por consiguiente, la totalidad
del Hidrógeno y Metano que sale del reactor también sale del proceso en la corriente de
purga y dado que el Hidrógeno es el reactivo en exceso, la fracción molar del Hidrógeno
en la purga toma su valor máximo, yPH,máx, para cualquier valor de x. Obviamente la
fracción molar del Metano en este límite toma su valor mínimo: yPM,mín = 1 - yPH,máx.
De modo que al sustituir ξ DIV =1.0 en la Ec. (23) se puede calcular el valor máximo
de la fracción molar del Hidrógeno en la purga para cualquier valor factible de la
conversión, x. De esta manera se define exactamente la región factible de valores que puede
tomar yPH antes de simular la Estructura de Recirculación del proceso HDA:

0  yPH  yPH , máx (24)

Como una consecuencia de las Ecs (15) y (24) se tiene que la fracción molar del
Metano en la purga varía en el rango: yPM,mín < yPM < 1.

Dr. Medardo Serna González – FIQ - UMSNH 5

Diseño de Procesos Químicos Capítulo 3

La sustitución en la Ec. (23) de las expresiones explícitas de las variables que contiene
y de un valor de la fracción de división igual a 1.0
FH ,out   MR  PB  Sx    R1
  R 2 ,  R1  PB S ,  R 2  PB 1  S  2 S y  DIV  1.0

genera, después de la manipulación algebraica requerida para tal fin, una relación explícita
para calcular yPH,máx:

2  MR  y FH  1  S  y FH x
yPH ,máx  (25)
2MR  1  S  yFH x

Las variables independientes MR y yFH del proceso HDA que aparecen en la Ec. (25)
son conocidas. Adicionalmente, no aparece en esta expresión el flujo molar del benceno
producido, PB, y dado que la selectividad S sólo es función de la conversión, x, se concluye
que la variable yPH,máx también depende únicamente de x. Por lo tanto, la Ec. (25) establece
claramente que para cada valor de x existe un valor máximo de la fracción molar del
Hidrógeno en la corriente de purga.
Cuando se alimenta Hidrógeno puro al proceso, la Ec. (25) se simplifica a

2MR  1  S  x
yPH ,máx  (26)
2MR  1  S  x

Del análisis matemático de las Ecs. (25) y (26) se concluye que para cualquier valor
de x siempre se cumple la condición: yPH,máx < yFH.
Si se asignan valores a yPH superiores a su valor máximo fijado por la Ec. (25) o la Ec.
(26), según sea el caso, entonces los flujos molares de la corriente gaseosa de
recirculación serán negativos y la fracción del divisor  DIV será mayor que 1.0, lo que
no tiene significado físico. Por lo tanto, antes de simular la Estructura de Recirculación de la
opción 3 del proceso HDA, para cualquier valor de x se debe calcular yPH,máx a fin de definir
a priori la región de valores factibles de yPH.

EJEMPLO 3.1 Cálculo de yPH,máx en Función de la Conversión x

Desarrolle un programa en MATLAB para calcular la fracción molar máxima del
Hidrógeno en la corriente de purga para diferentes valores de la conversión x. La razón molar
MR y la fracción molar yFH son iguales a 5 y 0.95, respectivamente.

Solución.
La función de MATLAB YPHmax.m calcula la fracción molar máxima del Hidrógeno
en la corriente de purga como una función de la conversión.
Con el propósito de verificar los resultados dados por la Ec. (25) se emplea la Ec. (23)
con  DIV  1 . Ambas ecuaciones dan los mismos resultados, los cuales se presentan en la
Tabla 3.1. A medida que se incrementa la conversión x, el valor de yPH,máx disminuye. Para
todos los valores de x, se cumple la condición: yPH,máx < yFH.

Dr. Medardo Serna González – FIQ - UMSNH 6

Diseño de Procesos Químicos Capítulo 3

La función YPHmax.m requiere como datos tanto a x como a PB para calcular las
variables que aparecen en el lado derecho de la Ec. (23), por lo que se podría concluir
erróneamente que el valor máximo de yPH es una función de esas dos variables. Sin embargo,
no es así, ya que yPH,máx depende únicamente de x de acuerdo a la Ec. (25). Esta dependencia
también se puede demostrar numéricamente al especificar diferentes valores para PB en la
función YPHmax.m a fin de calcular yPH,máx de la Ec. (23); los resultados obtenidos indican
que la variación de PB no altera los valores máximos de yPH para un valor dado de x.

Tabla 3.1 - yPH,máx en función de x

Conversión yPH,máx
0.05 0.9405011
0.15 0.9215052
0.25 0.9025130
0.35 0.8835271
0.45 0.8645525
0.55 0.8455996
0.65 0.8266949
0.75 0.8079189
0.85 0.7896442
0.95 0.7768968

Listado de la función YPHmax.m

function YPHmax

% Ejemplo 3.1
% Cálculo de la fracción molar máxima del Hidrógeno en la corriente de
% purga como una función de la conversión por paso del Tolueno
% MSG | FIQ | UMSNH - 19 de Octubre de 2020

clear, clc

% -----------------------------------------------------------------------
% DATOS
% -----------------------------------------------------------------------
% Flujo molar de salida del Benceno - Producto principal
PB = 100; % kmol/h

% Fracción molar del Hidrógeno en la corriente de alimentación

YFH = 0.95;

% Razón molar de Hidrógeno/Tolueno en la entrada del reactor

MR = 5;

% Conversión por paso del Tolueno en el reactor

x=(0.05:0.1:0.95);

% -----------------------------------------------------------------------
% CÁLCULO DE YPHmax = f(x) | Capítulo 3

Dr. Medardo Serna González – FIQ - UMSNH 7

Diseño de Procesos Químicos Capítulo 3

% -----------------------------------------------------------------------
S = 1 - 0.0036./((1 - x).^1.544);
YPHmax1 = (2*MR*YFH - YFH*x.*(1 + S))./(2*MR + YFH*x.*(1 - S)); % Ec.
(25)

% Verificación de los resultados

Eps_R1 = PB./S;
Eps_R2 = PB.*(1 - S)./(2*S);
FHout = (MR*PB)./(S.*x) - Eps_R1 + Eps_R2;
% Ec. (23) del Capítulo 3 haciendo Ediv = 1
YPHmax2 = FHout./(FHout + ((1 - YFH)/YFH).*(FHout + Eps_R1 - Eps_R2) +
Eps_R1);

% -----------------------------------------------------------------------
% VISUALIZACIÓN TABULAR DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------
Results=[x; YPHmax1; YPHmax2];
disp([' RESULTADOS PARA YPHmax']);
disp(' x YPHmax1 YPHmax2')
fprintf('%4.2f %12.7f %12.7f\n',Results)
disp(' ')

end % function YPHmax

EJEMPLO 3.2 Variación de yPH en Función de ξ DIV

Para una conversión de Tolueno en el reactor de 0.75, determine la variación de la
fracción molar del Hidrógeno y, por lo tanto, del Metano en la corriente de purga como una
función de la fracción de división.

Solución.
La Ec. (23) se resuelve usando la función de MATLAB YPH_vs_Ediv.m para
diferentes valores de la fracción de división que se encuentran en el rango

0   DIV  1.0 (14)

La fracción molar del Metano en la purga se calcula de

y PM  1  y PH (15)

Los resultados obtenidos se presentan en la Figura 3.3, donde se observa que la

fracción molar del Hidrógeno en la purga se incrementa al aumentar la fracción de
división. Como se esperaba, el valor máximo de yPH se obtiene cuando  DIV  1 y es menor
que 1.
La variación de la fracción molar del Metano en la purga con respecto a  DIV es
contraria a la del Hidrógeno, esto es, la fracción molar yPM en la purga y en la recirculación
disminuye cuando aumenta la fracción de división. Por consiguiente, la fracción molar del
Metano alcanza su valor mínimo cuando  DIV = 1.

Dr. Medardo Serna González – FIQ - UMSNH 8

Diseño de Procesos Químicos Capítulo 3

Figura 3.3 – Fracción molar del Hidrógeno y Metano en la purga en función de  DIV para x = 0.75

El análisis cuidadoso de la Figura 3.3 ayuda a predecir el funcionamiento general del

proceso HDA. En particular, dado que al aumentar  DIV se incrementa el valor de yPH para un
valor dado de x, es posible percibir y describir cualitativamente varias interacciones que
existen entre las variables del proceso:
 Al aumentar yPH se incrementa el flujo molar de la purga y la pérdida de Hidrógeno
en la purga, por lo que también se incrementa el flujo molar de la alimentación
gaseosa rica en Hidrógeno. Este comportamiento contribuye a la disminución del
potencial económico del proceso.
 Al aumentar yPH disminuye el flujo de la corriente gaseosa de recirculación, por lo
que también disminuye la potencia al freno del compresor y, por tanto, el costo total
de compresión. Este comportamiento contribuye al incremento del potencial
económico del proceso.
 Al aumentar yPH disminuye el flujo de la alimentación al reactor. Este comportamiento
contribuye a la disminución marginal del volumen del reactor y, a la vez, al
incremento de la temperatura de salida del reactor.
La verificación numérica de estas interacciones se presentará en el Ejemplo 3.3.

Listado de la función YPH_vs_Ediv.m

function YPH_vs_Ediv

% Ejemplo 3.2
% Variación de la fracción molar del Hidrógeno en la corriente de purga
% como una función de la conversión por paso del Tolueno
% MSG | FIQ | UMSNH - 18 de Octubre de 2020

Dr. Medardo Serna González – FIQ - UMSNH 9

Diseño de Procesos Químicos Capítulo 3

clear, clc

% -----------------------------------------------------------------------
% DATOS - Variables Independientes
% -----------------------------------------------------------------------
% Flujo molar de salida del Benceno - Producto principal
PB = 100; % kmol/h

% Fracción molar del Hidrógeno en la corriente de alimentación

YFH = 0.95;

% Razón molar de Hidrógeno/Tolueno en la entrada del reactor

MR = 5;

% Conversión por paso del Tolueno en el reactor

x = 0.75;

% Fracción de división
Ediv = (0:0.01:1);

% -----------------------------------------------------------------------
% CÁLCULO DE YPH = f(Ediv) para un valor dado de x | Capítulo 3
% -----------------------------------------------------------------------
S = 1 - 0.0036/((1 - x)^1.544);
FHout = (MR*PB)/(S*x);
Eps_R1 = PB/S;
Eps_R2 = PB*(1 - S)/(2*S);
% Ec. (23) del Capítulo 3
YPH = FHout*Ediv./(FHout*Ediv + ((1 - YFH)/YFH)*(FHout*Ediv + Eps_R1 -
Eps_R2) + Eps_R1);
% Ec. (15) del Capítulo 3
YPM = 1 - YPH;

% -----------------------------------------------------------------------
% VISUALIZACIÓN TABULAR DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------
Results=[Ediv; YPH; YPM];
disp(['RESULTADOS - YPH y YPM vs Ediv']);
disp(' ')
disp('Ediv YPH YPM')
fprintf('%4.2f %12.6f %12.6f\n',Results)

% -----------------------------------------------------------------------
% VISUALIZACIÓN GRÁFICA DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------
plot (Ediv, YPH, ':ok', Ediv, YPM, '--xr')
legend ('YPH', 'YPM')
title ('Variación de YPH y YPM en función de Ediv | x = 0.75')
xlabel ('Fracción de división')
ylabel ('Fracción molar del Hidrógeno y Metano')
grid on
end % function YPH_vs_Ediv

Dr. Medardo Serna González – FIQ - UMSNH 10

Diseño de Procesos Químicos Capítulo 3

EJEMPLO 3.3 Visualización de las Interacciones entre las Variables del

Proceso HDA
Desarrolle un programa en MATLAB para calcular los balances de materia de la
Estructura de Recirculación de la opción 3 del proceso HDA. La razón molar MR y la fracción
molar yFH son iguales a 5 y 0.95, respectivamente. La producción de Benceno PB es igual a
100 kmol/h.
Para x = [0.15, 0.50, 0.75] y los valores de la fracción molar yPH que se encuentran en
el rango (0.05:0.05:0.80), grafique los siguientes resultados:
(a) La fracción divisor  DIV en función de yPH.
(b) El flujo molar de la corriente de purga PG en función de yPH.
(c) El flujo molar de la corriente de recirculación RG en función de yPH.
(d) El flujo molar de la alimentación total al reactor FRct_in en función de yPH.

Para yPH = [0.15, 0.50, 0.75] y los valores de la conversión x = (0.05:0.05:0.95),

grafique los siguientes resultados:
(e) La fracción divisor  DIV en función de x.
(f) El flujo molar de la corriente de purga PG en función de x.
(g) El flujo molar de la corriente de recirculación RG en función de x.
(h) El flujo molar de la alimentación total al reactor FRct_in en función de x.

Solución.
La función de MATLAB Recycle_Procesos_HDA_1a.m implementa los cálculos de
los balances de materia y grafica los resultados obtenidos. Las Figuras 3.4 a 3.11
corresponden a las gráficas solicitadas en los incisos (a) a (h).

Figura 3.4 – Variación de  DIV con yPH para diferentes valores de x

Dr. Medardo Serna González – FIQ - UMSNH 11

Diseño de Procesos Químicos Capítulo 3

Figura 3.5 – Variación de PG con yPH para diferentes valores de x

Figura 3.6 – Variación de RG con x para diferentes valores de yPH

Figura 3.7 – Variación de FRctIN con yPH para diferentes valores de x

Dr. Medardo Serna González – FIQ - UMSNH 12

Diseño de Procesos Químicos Capítulo 3

Figura 3.8 – Variación de  DIV con x para diferentes valores de yPH

Figura 3.9 – Variación de PG con x para diferentes valores de yPH

Figura 3.10 – Variación de RG con x para diferentes valores de yPH

Dr. Medardo Serna González – FIQ - UMSNH 13

Diseño de Procesos Químicos Capítulo 3

Figura 3.11 – Variación de FRctIN con x para diferentes valores de yPH

Listado de la función Recycle_ProcesoHDA_1a.m

function Recycle_ProcesoHDA_1a

clear, clc

% Estructura de Recirculación del Proceso HDA

% Opción 3 - Figura 3.1 de las NOTAS
% Cálculo de los flujos de las corrientes internas para diferentes
valores
% de la conversión (x) y de la fracción molar del Hidrógeno en la
corriente
% de purga (YPH)
% Cálculo secuencial de las ecuaciones de balance de materia (ver NOTAS)
% MSG - FIQ - UMSNH | 19 de Octubre de 2020

% -----------------------------------------------------------------------
% DATOS - 5 GRADOS DE LIBERTAD
% -----------------------------------------------------------------------
% (1) Flujo molar de salida del Benceno - Producto principal
PB = 100; % kmol/h

% (2) Fracción molar del Hidrógeno en la corriente de alimentación

YFH = 0.95;

% (3) Fracción molar del Hidrógeno en la corriente de purga

YPH = (0.05:0.05:0.8);

% (4) Conversión por paso del Tolueno en el reactor

x = (0.05:0.05:0.95);

% (5) Razón molar H2/Tolueno en la entrada del reactor

MR = 5;

% -----------------------------------------------------------------------
% CÁLCULOS DE LOS BALANCES DE MATERIA

Dr. Medardo Serna González – FIQ - UMSNH 14

Diseño de Procesos Químicos Capítulo 3

% -----------------------------------------------------------------------

for i=1:length(YPH)

for j=1:length(x)
% ---------------------------------------------------------------
% BALANCES GLOBALES DE MATERIA - MODELO EN CAPÍTULO 2
% ---------------------------------------------------------------
S(j) = 1 - 0.0036/((1 - x(j))^1.544); % Ec. (49)
Eps_R1(j) = PB/S(j); % Ec. (48)
Eps_R2(j) = (PB/2)*((1 - S(j))/S(j)); % Ec. (45a)
FFT(j) = PB/S(j); % Ec. (46a)
PD(j) = (PB/2)*((1 - S(j))/S(j)); % Ec. (47a)
FG(i,j) = (PB/S(j))*(1-(1-YPH(i))*(1-S(j))/2)/(YFH-YPH(i)); %
Ec. (43b)
PG(i,j) = FG(i,j)+(PB/2)*(1/S(j)-1); % Ec. (44b)

% ---------------------------------------------------------------
% CÁLCULOS DE LA ESTRUCTURA DE RECIRCULACIÓN-MODELO EN CAPÍTULO 3
% ---------------------------------------------------------------
% Flujos molares de los compuestos alimentados al reactor
% y de la corriente de recirculación gaseosa
FTin(j) = FFT(j)/x(j); % Ec. (3)
FHin(j) = MR*FTin(j); % Ec. (5)
RG(i,j) = (FHin(j) - YFH*FG(i,j))/YPH(i); % Ec. (4)
FMin(i,j) = (1 - YFH)*FG(i,j) + (1 - YPH(i))*RG(i,j); % Ec. (8)
% FRct_in: Flujo molar de la alimentación al reactor, kmol/h
FRct_in(i,j) = FHin(j) + FMin(i,j) + FTin(j);

% Flujos molares de los compuestos a la salida del reactor

FHout(j) = FHin(j) - Eps_R1(j) + Eps_R2(j); % Ec. (9)
FMout(i,j) = FMin(i,j) + Eps_R1(j); % Ec. (10)
FBout = PB; % Ec. (11)
FTout(j) = FTin(j) - Eps_R1(j); % Ec. (12)
FDout(j) = Eps_R2(j); % Ec. (13)

% Fracción de división - Divisor

Ediv(i,j) = PG(i,j)/(PG(i,j) + RG(i,j));
end
end

% -----------------------------------------------------------------------
% VISUALIZACIÓN TABULAR DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------

for i=1:length(YPH)
Results1=[x; S; Eps_R1; Eps_R2; FFT; FG(i,:); PD; PG(i,:); RG(i,:);
FRct_in(i,:); Ediv(i,:)];
disp([' RESULTADOS PARA YPH =
',num2str(YPH(i))]);
disp(' x S Eps_R1 Eps_R2 FFT FG PD PG
RG FRin Ediv')
fprintf('%4.2f %7.4f %8.3f %6.3f %8.3f %8.3f %6.3f %8.3f %10.3f
%11.3f %8.4f\n',Results1)
disp(' ')
end

Dr. Medardo Serna González – FIQ - UMSNH 15

Diseño de Procesos Químicos Capítulo 3

% -----------------------------------------------------------------------
% VISUALIZACIÓN GRÁFICA DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------

%
plot (YPH, Ediv(:,3), ':ok', YPH, Ediv(:,10), '--xr', YPH, Ediv(:,15), '-
b')
legend ('x=0.15', 'x=0.5','x=0.75')
title ('Figura 3.4 - Variación de Ediv con YPH para diferentes valores de
x')
xlabel ('Fracción molar del Hidrógeno en la purga')
ylabel ('Fracción de división - Divisor')
grid on
%
pause
clf
%
plot (YPH, PG(:,3), ':ok', YPH, PG(:,10), '--xr', YPH, PG(:,15), '-b')
legend ('x=0.15', 'x=0.5','x=0.75')
title ('Figura 3.5 - Variación de PG con YPH para diferentes valores de
x')
xlabel ('Fracción molar del Hidrógeno en la purga')
ylabel ('Flujo molar de la purga, kmol/h')
grid on
%
pause
clf
%
plot (YPH, RG(:,3), ':ok', YPH, RG(:,10), '--xr', YPH, RG(:,15), '-b')
legend ('x=0.15', 'x=0.5','x=0.75')
title ('Figura 3.6 - Variación de RG con YPH para diferentes valores de
x')
xlabel ('Fracción molar del Hidrógeno en la purga')
ylabel ('Flujo molar de la corriente de recirculación, kmol/h')
grid on
%
pause
clf

plot (YPH, FRct_in(:,3), ':ok', YPH, FRct_in(:,10), '--xr', YPH,

FRct_in(:,15), '-b')
legend ('x=0.15', 'x=0.5','x=0.75')
title ('Figura 3.7 - Variación de FRct-In con YPH para diferentes valores
de x')
xlabel ('Fracción molar del Hidrógeno en la purga')
ylabel ('Flujo molar de la alimentación al reactor, kmol/h')
grid on
%
pause
clf
%
plot (x, Ediv(3,:), ':ok', x, Ediv(10,:), '--xr', x, Ediv(15,:), '-b')
legend ('YPH=0.15', 'YPH=0.5','YPH=0.75')
title ('Figura 3.8 - Variación de Ediv con x para diferentes valores de
YPH')
xlabel ('Conversión por paso del Tolueno en el reactor')

Dr. Medardo Serna González – FIQ - UMSNH 16

Diseño de Procesos Químicos Capítulo 3

ylabel ('Fracción de división - Divisor')

grid on
%
pause
clf
%
plot (x, PG(3,:), ':ok', x, PG(10,:), '--xr', x, PG(15,:), '-b')
legend ('YPH=0.15', 'YPH=0.5','YPH=0.75')
title ('Figura 3.9 - Variación de PG con x para diferentes valores de
YPH')
xlabel ('Conversión por paso del Tolueno en el reactor')
ylabel ('Flujo molar de la purga, kmol/h')
grid on
%
pause
clf
%
plot (x, RG(3,:), ':ok', x, RG(10,:), '--xr', x, RG(15,:), '-b')
legend ('YPH=0.15', 'YPH=0.5','YPH=0.75')
title ('Figura 3.10 - Variación de RG con x para diferentes valores de
YPH')
xlabel ('Conversión por paso del Tolueno en el reactor')
ylabel ('Flujo molar de la recirculación, kmol/h')
grid on
%
pause
clf
%
plot (x, FRct_in(3,:), ':ok', x, FRct_in(10,:), '--xr', x, FRct_in(15,:),
'-b')
legend ('YPH=0.15', 'YPH=0.5','YPH=0.75')
title ('Figura 3.11 - Variación de FRct-In con x para diferentes valores
de YPH')
xlabel ('Conversión por paso del Tolueno en el reactor')
ylabel ('Flujo molar de la alimentación al reactor, kmol/h')
grid on
%
pause

clf

end % function Recycle_ProcesoHDA_1a

3.3 EFECTOS TÉRMICOS DEL REACTOR

Los balances de materia y de energía de la Estructura de Recirculación bajo estudio
están desacoplados. Las variables independientes seleccionadas, esto es, PB, yPH, yFH, x y MR
causan el particionado de las ecuaciones del modelo, lo que permite calcular primeramente
los flujos molares de las corrientes de entrada, de salida e internas del proceso y, en segundo
lugar, determinar los efectos térmicos del reactor al resolver el balance de energía
correspondiente.
Cuando se ignoran los cambios de energía cinética y potencial, el balance de energía
del reactor es dado por

Dr. Medardo Serna González – FIQ - UMSNH 17

Diseño de Procesos Químicos Capítulo 3

QRCT  F1,out H 1,out  F2,out H 2,out  F3,out H 3,out  F4,out H 4,out  F5,out H 5,out
 F1,in H 1,in  F2,in H 2, in  F4,in H 4,in (27)
 R1H Ro1   R 2 H Ro 2
En esta ecuación los compuestos han sido representados por subíndices numéricos
en vez de por sus fórmulas químicas, esto es:
H2 CH 4 C6 H 6 C 7 H8 C12 H10
1 2 3 4 5
Después de resolver los balances de materia del proceso HDA, los flujos molares y
los avances de reacción que aparecen en la Ec. (27) son variables conocidas. Por lo tanto,
cuando el sistema reaccionante se comporta como un gas ideal (esto es, si las entalpías
específicas de los compuestos dependen únicamente de la temperatura), el balance de energía
incorpora al modelo matemático solamente 1 ecuación independiente y 3 nuevas variables:
(a) El calor transferido entre el reactor y el medio ambiente, QRCT.
(b) La temperatura de entrada del reactor, Tin.
(c) La temperatura de salida del reactor, Tout.
Ambas temperaturas son variables implícitas en la Ec. (27), pero aparecen
explícitamente como variables en las ecuaciones constitutivas que definen las entalpías
específicas de los compuestos, como se muestra en el Ejemplo 3.4.
El modelo matemático resultante tiene ahora una ecuación y tres variables más, por
lo que el sistema tiene 7 variables independientes o grados de libertad. Esto significa que
es necesario especificar los valores de dos de las tres variables del balance de energía para
poder calcular la variable restante usando la Ec. (27).

EJEMPLO 3.4 Temperatura de Salida del Reactor del Proceso HDA

Cuando el reactor opera adiabáticamente (QRCT = 0) y es conocida la temperatura de
entrada, Tin, la temperatura de salida es la única variable desconocida.
Desarrolle un programa en MATLAB para calcular la temperatura de salida del
reactor de la opción 3 del proceso HDA suponiendo que el sistema reaccionante es un gas
ideal para las siguientes especificaciones 1:

PB , kmol/h y FH yPH x MR Tin , K QRCT , J/s

100 0.95 0.4 0.75 5 894.2611 0

Solución
Se han especificado los valores de 7 variables independientes del modelo matemático
que permite resolver primero los balances de materia y, al final, el balance de energía.

1
La fuente de las propiedades termodinámicas de los compuestos es el Apéndice A del libro: Bruce E. Poling, John M.
Prausnitz and John P. O´Connell. The Properties of Gases and Liquids. Fith edition. McGraw-Hill (2001).

Dr. Medardo Serna González – FIQ - UMSNH 18

Diseño de Procesos Químicos Capítulo 3

Paso 1. Resolver los balances de materia del proceso para determinar los flujos
molares de los compuestos presentes en las corrientes de entrada al y de salida del reactor,
así como ambos avances de reacción. Los resultados obtenidos se reportan en las segunda y
cuarta columnas, así como en las últimas dos filas de la siguiente tabla:

Fi ,in Hi ,in Fi , out H i ,out

Compuesto kmol/h kmol/h
kJ/kmol kJ/kmol
Hidrógeno - 1 687.719 586.140
Metano - 2 776.052 879.210
Benceno - 3 ------- ------- 100.000
Tolueno - 4 137.544 34.386
Difenilo - 5 ------- ------- 1.579
Avance de la reacción 1,  R1 = 103.158 kmol/h
Avance de la reacción 2,  R 2 = 1.579 kmol/h

Paso 2. Calcular las entalpías específicas de los compuestos en la entrada del reactor.
Tin Cpi
H i ,in  R  dT
298.15
R
 a1i 2 a 2i 3 3 
 a 0i Tin  298.15   2 Tin  298.15   3 Tin  298.15   (28)
2

 R
 a3i T 4  298.154  a 4i T 5  298.155 
 4
 in  5  in  

donde R = 8.3143 kJ kmol-1 K-1 y Tin = 894.2611 K. Queda claro que las entalpías calculadas
tienen como estado de referencia los reactivos como gases ideales a 25°C y 1 atm.
Los valores obtenidos de las entalpías específicas de los compuestos alimentados al
reactor se reportan en la tercera columna de la tabla de resultados:

Fi ,in Hi ,in Fi , out H i ,out

Compuesto kmol/h kmol/h
kJ/kmol kJ/kmol
Hidrógeno - 1 687.719 17500.059 586.140
Metano - 2 776.052 31020.788 879.210
Benceno - 3 ------- ------- 100.000
Tolueno - 4 137.544 112437.033 34.386
Difenilo - 5 ------- ------- 1.579
Avance de la reacción 1,  R1 = 103.158 kmol/h
Avance de la reacción 2,  R 2 = 1.579 kmol/h
Estados de referencia: Reactivos y productos como gases ideales a 25°C y 1 atm.

Paso 3. Calcular los calores estándar de reacción a una temperatura de 25°C y una
presión total de 1 atm considerando como gases ideales a los reactivos y productos.

Dr. Medardo Serna González – FIQ - UMSNH 19

Diseño de Procesos Químicos Capítulo 3

H Ro1  H of 2  H of 3  H of 1  H of 4 (29)

H Ro 2  H of 5  H of 1  2H of 3 (30)

Paso 4. Calcular la temperatura de salida del reactor.

Para cada uno de los cinco compuestos en la corriente de salida la entalpía específica
es dada por

Cpi
Tout
H i ,out  R  dT
298.15
R
 a1i 2 a 2i 3 3 
 a 0i Tout  298.15   2 Tout  298.15   3 Tout  298.15   (31)
2

 R
 a3i T 4  298.154  a 4i T 5  298.155 
 4  out  5  out  

En las ecuaciones constitutivas de las entalpías específicas de salida la única incógnita

es la temperatura de salida, Tout.
Las expresiones resultantes para las entalpías específicas de los compuestos a la salida
del reactor, así como los valores calculados de los flujos molares de los compuestos, las
entalpías específicas a la entrada del reactor, los calores de reacción y los avances de ambas
reacciones se sustituyen en la Ec. (27). Así se llega a una ecuación no lineal que se resuelve
numéricamente para determinar la temperatura de salida, Tout; una vez conocida esta
temperatura se calculan las entalpías específicas de salida para cada uno de los cinco
compuestos usando la Ec. (31). Las entalpías de salida calculadas se reportan en la última
columna de la tabla de resultados:

Fi ,in Hi ,in Fi , out H i ,out

Compuesto kmol/h kmol/h
kJ/kmol kJ/kmol
Hidrógeno - 1 687.719 17500.059 586.140 18901.507
Metano - 2 776.052 31020.788 879.210 34261.172
Benceno - 3 ------- ------- 100.000 100845.154
Tolueno - 4 137.544 112437.033 34.386 124104.483
Difenilo - 5 ------- ------- 1.579 194408.354
Avance de la reacción 1,  R1 = 103.158 kmol/h
Avance de la reacción 2,  R 2 = 1.579 kmol/h
Calor estándar de la reacción principal (R1) = -41810 kJ/kmol
Calor estándar de la reacción secundaria (R2) = 16660 kJ/kmol
Temperatura de entrada al reactor = 894.2611 K
Temperatura de salida del reactor = 941.1211 K
Estados de referencia: Reactivos y productos como gases ideales a 25°C y 1 atm.

Dr. Medardo Serna González – FIQ - UMSNH 20

Diseño de Procesos Químicos Capítulo 3

Listado de la función Reactor_Tout_HDA.m

function Reactor_Tout_HDA

% Temperatura de salida del reactor adiabático del proceso HDA - Opción 3

% El fluido reaccionante se considera como un gas ideal
% Propiedades termodinámicas tomadas del Apéndice A de Poling y col.
(2001)
% MSG | FIQ | UMSNH - 28 de Septiembre de 2020

clear, clc

% -----------------------------------------------------------------------
% DATOS - 7 Variables Independientes y Propiedades Termodinámicas
% -----------------------------------------------------------------------
PB = 100; % Flujo de Benceno producido en kmol/h
YFH = 0.95; % Fracción molar del H2 en la alimentación
gaseosa
x = 0.75; % Conversión por paso del Tolueno en el
reactor
YPH = 0.4; % Fracción molar del Hidrógeno en la purga
MR = 5; % Razón molar H2/Tolueno en la entrada del
reactor
Tin = (1150 - 32)/1.8 + 273.15; % Temperatura de entrada en K (1150 °F)
Q_Rct = 0; % Reactor adiabático

% Indices asignados a los compuestos del proceso

% Hidrógeno - 1
% Metano - 2
% Benceno - 3
% Tolueno - 4
% Difenilo - 5

% Coeficientes estequiométricos
Coef_R1 = [-1 1 1 -1 0];
Coef_R2 = [1 0 -2 0 1];

% Calores de formación
DHf = [0 -74.52 82.88 50.17 182.42]; % kJ/mol

% Capacidades caloríficas Cp/R

CpG = [2.883 3.681e-3 -0.772e-5 0.692e-8 -0.213e-11;
4.568 -8.975e-3 3.631e-5 -3.407e-8 1.091e-11;
3.551 -6.184e-3 14.365e-5 -19.807e-8 8.234e-11;
3.866 3.558e-3 13.356e-5 -18.659e-8 7.690e-11;
-0.843 61.392e-3 6.352e-5 -13.754e-8 6.169e-11];

% -----------------------------------------------------------------------
% Flujos de Entrada y Salida del Reactor y Avances de Reacción Calculados
% -----------------------------------------------------------------------
FH = 687.719; % Flujo del Hidrógeno alimentado al reactor, kmol/h
FM = 776.052; % Flujo del Metano alimentado al reactor, kmol/h
FT = 137.544; % Flujo del Tolueno alimentado al reactor, kmol/h

FHout = 586.140; % Flujo del Hidrógeno a la salida del reactor, kmol/h

Dr. Medardo Serna González – FIQ - UMSNH 21

Diseño de Procesos Químicos Capítulo 3

FMout = 879.210; % Flujo del Metano a la salida del reactor, kmol/h

FBout = 100.000; % Flujo del Benceno a la salida del reactor, kmol/h
FTout = 34.386; % Flujo del Tolueno a la salida del reactor, kmol/h
FDout = 1.579; % Flujo del Difenilo a la salida del reactor, kmol/h

Eps_R1 = 103.158; % Avance de la reacción 1, kmol/h

Eps_R2 = 1.579; % Avance de la reacción 2, kmol/h

% -----------------------------------------------------------------------
% CALCULO DE LA TEMPERATURA DE SALIDA DEL REACTOR
% -----------------------------------------------------------------------
Fin = [FH FM 0 FT 0]; % Fi,in - kmol/h
Fout = [FHout FMout FBout FTout FDout]; % Fi,out - kmol/h
Tref = 298.15; % Temperatura de referencia,
K
Hin = [Hcalc(CpG(1,:),Tref,Tin) Hcalc(CpG(2,:),Tref,Tin) 0 ...
Hcalc(CpG(4,:),Tref,Tin) 0]; % Hi,in en kJ/kmol, Ec. (15)
Hrxn_R1 = 1000*Coef_R1*DHf'; % Ec. (16) - kJ/kmol
Hrxn_R2 = 1000*Coef_R2*DHf'; % Ec. (17) - kJ/kmol

Tout_s = Tin + 50; % Estimado de la temperatura de salida (Tout >

Tin)

[Tout] = fzero(@EB_Function,Tout_s); % Temperatura de salida del reactor,

K

% -----------------------------------------------------------------------
% VISUALIZACIÓN DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------
disp('Valores Especificados para las Variables Independientes')
fprintf(' Benceno producido = %5.2f kmol/h\n',PB)
fprintf(' Fracción molar del Hidrógeno en la alimentación =
%.2f\n',YFH)
fprintf(' Conversión por paso del Tolueno = %.2f\n',x)
fprintf(' Fracción molar del Hidrógeno en la purga = %.2f\n',YPH)
fprintf(' Razón molar H2/Tolueno en la alimentación al reactor =
%4.2f\n',MR)
fprintf(' Temperatura de entrada - Reactor = %7.4f K\n',Tin)
fprintf(' Flujo de calor transferido entre el reactor y el medio
ambiente = %4.2f\n',Q_Rct)
disp(' ')
disp('Valores Calculados de las Variables Dependientes')
fprintf(' Hidrógeno alimentado al reactor = %7.3f kmol/h \n',FH)
fprintf(' Metano alimentado al reactor = %7.3f kmol/h \n',FM)
fprintf(' Tolueno alimentado al reactor = %7.3f kmol/h \n',FT)
disp(' ')
fprintf(' Avance de la reacción principal (R1) = %7.3f kmol/h
\n',Eps_R1)
fprintf(' Avance de la reacción secundaria (R2) = %7.3f kmol/h
\n',Eps_R2)
disp(' ')
fprintf(' Entalpía del Hidrógeno alimentado al reactor = %7.3f kJ/kmol
\n',Hin(1))
fprintf(' Entalpía del Metano alimentado al reactor = %7.3f kJ/kmol
\n',Hin(2))

Dr. Medardo Serna González – FIQ - UMSNH 22

Diseño de Procesos Químicos Capítulo 3

fprintf(' Entalpía del Tolueno alimentado al reactor = %7.3f kJ/kmol

\n',Hin(4))
disp(' ')
fprintf(' Calor estándar de la reacción principal (R1) = %7.3f kJ/kmol
\n',Hrxn_R1)
fprintf(' Calor estándar de la reacción secundaria (R2) = %7.3f
kJ/kmol \n',Hrxn_R2)
disp(' ')
fprintf(' Hidrógeno a la salida del reactor = %7.3f kmol/h \n',FHout)
fprintf(' Metano a la salida del reactor = %7.3f kmol/h \n',FMout)
fprintf(' Benceno a la salida del reactor = %7.3f kmol/h \n',FBout)
fprintf(' Tolueno a la salida del reactor = %7.3f kmol/h \n',FTout)
fprintf(' Difenilo a la salida del reactor = %7.3f kmol/h \n',FDout)
disp(' ')
fprintf(' Temperatura de salida del reactor = %8.4f K\n',Tout)
disp(' ')
fprintf(' Entalpía del Hidrógeno a la salida del reactor = %7.3f
kJ/kmol \n',Hout(1))
fprintf(' Entalpía del Metano a la salida del reactor = %7.3f kJ/kmol
\n',Hout(2))
fprintf(' Entalpía del Benceno a la salida del reactor = %7.3f kJ/kmol
\n',Hout(3))
fprintf(' Entalpía del Tolueno a la salida del reactor = %7.3f kJ/kmol
\n',Hout(4))
fprintf(' Entalpía del Difenilo a la salida del reactor = %7.3f
kJ/kmol \n',Hout(5))

%
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
function [fTout] = EB_Function(Tout)

Hout = [Hcalc(CpG(1,:),Tref,Tout) Hcalc(CpG(2,:),Tref,Tout)

Hcalc(CpG(3,:),Tref,Tout) ...
Hcalc(CpG(4,:),Tref,Tout) Hcalc(CpG(5,:),Tref,Tout)]; % Hi,out
en kJ/kmol, Ec. (18)

fTout = Fout*Hout' - Fin*Hin' + Eps_R1*Hrxn_R1 + Eps_R2*Hrxn_R2 -

Q_Rct; % Ec. (14)

end % EB_Function

%
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
function [H] = Hcalc(Cp,T1,T2)
% This function calculates enthalpy change given Tl, T2 and the Cp
% constants. The units of R determine the units of H.
% The form of the polynomial as given by Poling et al (2001).
% Cp/R = Cp(1) + Cp(2)*T + Cp(3)*T^2 + Cp(4)*T^3 + Cp(5)*T^4;
% Temperature in K
R = 8.3143; % kJ/kmol K
H = Cp(1)*(T2 - T1) + Cp(2)/2*(T2^2-T1^2) + Cp(3)/3*(T2^3-T1^3) ...
+ Cp(4)/4*(T2^4-T1^4) + Cp(5)/5*(T2^5-T1^5);
H = R*H; % kJ/kmol
end % function Hcalc

end % function Reactor_Tout_HDA

Dr. Medardo Serna González – FIQ - UMSNH 23

Diseño de Procesos Químicos Capítulo 3

3.4 DISEÑO DEL REACTOR

En esta sección se ilustra el procedimiento de diseño de un reactor tubular para el
proceso HDA.
Debido a que la cantidad de difenilo producido es muy pequeña, el diseño del reactor
considerará solamente la reacción principal del proceso HDA.
La cinética de la reacción principal es presentada por Silsby y col. 2 y por Zimmerman
3
y York :

r  k T C H 2  CTol 
1/ 2
(32)

donde la constante de reacción es dada por

 52, 000 cal / ( gmol ) 

k  6.3 x1010  gmol / L   s  exp 
1/ 2 1
 (33)
 RTR 

Para estimar el volumen del reactor se supone que opera isotérmicamente, a una
temperatura igual al promedio aritmético de las temperaturas de entrada y de salida, esto es:
TR = (Tin + Tout)/2. Esta aproximación requiere resolver solamente el balance de materia
diferencial representado por una ecuación diferencial para diseñar el reactor tubular:

dV F
  4,in (34)
dx r

EJEMPLO 3.5 Volumen del Reactor Tubular del Proceso HDA

Desarrolle un programa en MATLAB para determinar el volumen del reactor del
proceso HDA para las siguientes especificaciones de las variables independientes:

PB , kmol/h y FH yPH x MR Tin , K Tout , K P, atm

100 0.95 0.4 0.75 5 894.2611 941.1211 34.0229

Listado de la función Reactor_Design_HDA.m

function Reactor_Design_HDA

% Diseño del reactor del proceso HDA - Opción 3

% El fluido reaccionante se considera como un gas ideal
% Para simplificar los cálculos se supone que el reactor opera
isotérmicamente,
% a una temperatura T = (Tin + Tout)/2

2
R. I. Silsby y E. W. Sawyer; J. Appl. Chem., 6: 347 (August 1956); W.D. Betts, F. Popper y R. I. Silsby, J. Appl. Chem.,
7: 497 (September 1957).
3
C. C. Zimmerman y R. York, I&EC Proc. Des. Dev., 3: 254 (July 1964).

Dr. Medardo Serna González – FIQ - UMSNH 24

Diseño de Procesos Químicos Capítulo 3

% Propiedades termodinámicas tomadas del Apéndice A de Poling y col.

(2001)
% MSG | FIQ | UMSNH - 28 de Septiembre de 2020

clear, clc

% -----------------------------------------------------------------------
% DATOS - 8 Variables Independientes del proceso y los flujos molares de
% los compuestos a la entrada del reactor
% -----------------------------------------------------------------------
PB = 100; % Flujo de Benceno producido en kmol/h
YFH = 0.95; % Fracción molar del H2 en la alimentación gaseosa
x = 0.75; % Conversión por paso del Tolueno en el reactor
YPH = 0.4; % Fracción molar del Hidrógeno en la purga
MR = 5; % Razón molar H2/Tolueno en la entrada del reactor
Tin = (1150 - 32)/1.8 + 273.15; % Temperatura de entrada en K (1150 °F)
Tout = 941.121; % Temperatura de salida en K
P_Rct = 500/14.696; % Presión de operación del reactor en atm

FH = 687.719; % Flujo del Hidrógeno alimentado al reactor,

kmol/h
FM = 776.052; % Flujo del Metano alimentado al reactor, kmol/h
FT = 137.544; % Flujo del Tolueno alimentado al reactor, kmol/h

% -----------------------------------------------------------------------
% CÁLCULO DEL VOLUMEN DEL REACTOR
% -----------------------------------------------------------------------
% Indices asignados a los compuestos
% Hidrógeno - 1
% Metano - 2
% Benceno - 3
% Tolueno - 4

T_Rct = (Tin + Tout)/2; % Temperatura del reactor, K

V0 = 0; % Condición inicial del volumen

xspan = [0 x]; % Límites de la conversión
[x_Rct V] = ode45(@Rct_Function,xspan,V0); % Ec. (21)
V_Rct = V(end,1)*35.3145; % Volumen del reactor en ft3

% -----------------------------------------------------------------------
% VISUALIZACIÓN DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------
disp('Valores Especificados para las Variables Independientes')
fprintf(' Benceno producido = %5.2f kmol/h\n',PB)
fprintf(' Fracción molar del Hidrógeno en la alimentación =
%.2f\n',YFH)
fprintf(' Conversión por paso del Tolueno = %.2f\n',x)
fprintf(' Fracción molar del Hidrógeno en la purga = %.2f\n',YPH)
fprintf(' Razón molar H2/Tolueno en la alimentación al reactor =
%4.2f\n',MR)
fprintf(' Presión de operación del reactor = %7.4f atm\n',P_Rct)
fprintf(' Temperatura de entrada del reactor = %7.4f K\n',Tin)
fprintf(' Temperatura de salida del reactor = %7.2f\n',Tout)
disp(' ')
disp('Valores Calculados de las Variables Dependientes')

Dr. Medardo Serna González – FIQ - UMSNH 25

Diseño de Procesos Químicos Capítulo 3

fprintf(' Hidrógeno alimentado al reactor = %7.3f kmol/h \n',FH)

fprintf(' Metano alimentado al reactor = %7.3f kmol/h \n',FM)
fprintf(' Tolueno alimentado al reactor = %7.3f kmol/h \n',FT)
disp(' ')
fprintf(' Temperatura de operación del reactor = %7.3f K\n',T_Rct)
fprintf(' Volumen del reactor = %7.3f ft3\n',V_Rct)

%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
function dVdx = Rct_Function(xf,V)
% Flujos molares de los compuestos en función de la conversión xf
F1=FH - FT*xf;
F2=FM + FT*xf;
F3=FT*xf;
F4=FT*(1 - xf);

% Fracciones molares del H2 y Tolueno

Y1 = F1/(F1 + F2 + F3 + F4);
Y4 = F4/(F1 + F2 + F3 + F4);

CTot = P_Rct/(0.08206*T_Rct); % Concentración total, gmol/L

% R = 0.08206 L atm/(gmol K)

C1 = Y1*CTot; % Concentración del Hidrógeno, gmol/L

C4 = Y4*CTot; % Concentración del Tolueno, gmol/L

k = 6.3e10*exp(-52000/(1.987*T_Rct)); % (gmol/L)^(-1/2) (s)^(-1) -

Ec. (20)
% R = 1.987 cal/(gmol K)

Rate = -k*sqrt(C1)*C4; % gmol/(L s) - Ec. (20)

Rate = Rate*3600; % kmol/(m3 h)

dVdx = 0 - FT/Rate;
end % Rct_Function

end % function Reactor_Design_HDA

EJEMPLO 3.6 Modelo y Diseño de un Compresor

Formule un modelo matemático y desarrolle un programa en MATLAB para realizar
el diseño termodinámico del compresor de la corriente de recirculación gaseosa del proceso
HDA, suponiendo que esa corriente se comporta como un gas ideal. Dado que se supone que
el sistema de separación es perfecto, el Hidrógeno (compuesto 1) y el Metano (compuesto 2)
son los únicos compuestos presentes en esa corriente de recirculación.
Los valores especificados de las variables independientes del proceso son:

Pin , Pout ,
PB , kmol/h y FH yPH x MR Tin , °F C
psia psia
100 0.95 0.4 0.75 5 100 465 550 0.8

Dr. Medardo Serna González – FIQ - UMSNH 26

Diseño de Procesos Químicos Capítulo 3

Solución
Primero se formula el modelo matemático siguiendo un procedimiento sistemático
ilustrado paso a paso.

Paso 1. Definición del problema.

Los compresores son equipos usados para incrementar la presión de un gas o vapor y
consumen trabajo mecánico al hacerlo, que es proporcionado por una fuente externa a través
de un eje rotatorio. Un compresor tiene solamente una corriente de entrada y una corriente
de salida.
El objetivo es desarrollar un modelo matemático para llevar a cabo el diseño
termodinámico de compresores de gases ideales y, en particular, el diseño termodinámico del
compresor del proceso HDA para las variables independientes especificadas.

Paso 2. Suposiciones del modelo.

 El compresor opera en estado estacionario.
 Los cambios de energía cinética y potencial del gas entrando al y saliendo del
compresor son despreciables.
 No importa la descripción de los fenómenos de transporte internos que ocurren dentro
del compresor.
 Debido a que la compresión es muy rápida en comparación con la velocidad de
transferencia de calor, se supone que el compresor es adiabático.
 La corriente gaseosa se comporta como un gas ideal.

Paso 3. Formulación del modelo matemático.

Ecuaciones de Balance.
El compresor solamente tiene una entrada y una salida, por lo que el balance de
materia en una base molar para cada componente es trivial
F1,out  F1,in (35)

F2,out  F1,in (36)

Con las suposiciones enunciadas, el balance de energía en estado estacionario se

expresa como
We  H out  H in (37)

donde We es el trabajo de eje real requerido por el compresor, H in la entalpía total de la

corriente de entrada y H out es la entalpía total de la corriente de salida del compresor real

H out  F1,in H1,out  F2,in H 2,out (38)

H in  F1,in H1,in  F2,in H 2,in (39)

Dr. Medardo Serna González – FIQ - UMSNH 27

Diseño de Procesos Químicos Capítulo 3

donde Hi,in y Hi ,out son las entalpías específicas del compuesto i = 1, 2 (esto es: Hidrógeno,
Metano) en las corrientes de entrada y salida, respectivamente,
Cuando se asignan valores a nc flujos molares individuales, la temperatura y la presión
de la corriente de entrada4 y a la presión de la corriente de salida, que son las variables
independientes que usualmente se especifican cuando se diseña un compresor adiabático, los
balances de materia y energía no son suficientes para representar su funcionamiento dado
que la Ec. (37) contiene dos variables desconocidas We y T25. Se necesita, por tanto, más
información para fijar el estado final de una compresión adiabática. La información adicional
se obtiene usando el balance de entropía y la eficiencia isoentrópica del compresor.

Ecuaciones Constitutivas.
Si el compresor es reversible y adiabático, entonces toda la energía mecánica que
consume se emplea en incrementar la presión del gas o vapor. La compresión reversible y
adiabática es ideal, por lo que representa un límite teórico para la compresión real de un gas.
El balance de entropía para el compresor ideal del proceso HDA bajo estudio es dado por

rev
Sout  Sin
F1,in S1,revout  F2,in S2,revout  F1,in S1,in  F2,in S2,in (40)
f  Tout
rev
  F1,in S1,revout  F2,in S2,revout  F1,in S1,in  F2,in S2,in  0
rev
donde Sin es la entropía total de la corriente de entrada, Sout es la entropía total reversible
de la corriente de salida y, para cada compuesto i, Si ,in es la entropía específica en la corriente
de entrada y Sirev
, out es la entropía específica reversible o ideal en la corriente de salida.

La Ec. (40) establece que un compresor adiabático y reversible es isoentrópico. Esta

rev
ecuación es no lineal en la temperatura reversible de salida, Tout . Por lo tanto, para determinar
el valor de esta variable desconocida se debe resolver numéricamente la Ec. (40).
La fricción mecánica y la fricción del fluido son las principales fuentes de la
irreversibilidad que caracteriza a los compresores reales. La irreversibilidad está asociada a
pérdidas de energía, lo que ocasiona que los compresores reales consuman más potencia que
los compresores isoentrópicos.
Para comparar el funcionamiento de un compresor irreversible con el de un compresor
reversible y adiabático, ambos realizando la misma tarea de compresión entre el estado inicial
de la corriente de entrada y la presión de salida para un gas dado, se utiliza la eficiencia de
compresión isoentrópica,  C , que es definida como

Wrev
C  (41)
We

4
En este caso nc = 2, por lo que la corriente de entrada al compresor tiene 4 variables independientes.
5
Tener presente que, para gases ideales, las entalpías específicas de los compuestos en la salida del compresor Hi ,out y,

por tanto, la entalpía total de la corriente de salida, H out , dependen solamente de T2.

Dr. Medardo Serna González – FIQ - UMSNH 28

Diseño de Procesos Químicos Capítulo 3

donde el trabajo reversible es dado por

Wrev   H out
rev
- H in  (42)

rev
donde H out es la entalpía total de la corriente de salida para la compresión isoentrópica:

rev
H out  F1,in H1,revout  F2,in H 2,revout (43)

donde Hirev
,out es la entalpía reversible específica del compuesto i en la corriente de salida.

Las entalpías y entropías específicas para cada componente i considerado como un

gas ideal se expresan como

Tin Cpi
H i ,in  R  dT
298.15
R
 a1i 2 a 2i 3 3 
 a 0i Tin  298.15   2 Tin  298.15   3 Tin  298.15   (44)
2

 R
 a3i T 4  298.154  a 4i T 5  298.155 
 4  in  5  in  

Tout Cpi
H i ,out  R  dT
298.15
R
 a1i 2 a 2i 3 3 
 a 0i Tout  298.15   2 Tout  298.15   3 Tout  298.15   (45)
2

 R
 a3i T 4  298.154  a 4i T 5  298.155 
 4  out  5  out  

rev
Tout Cpi
H irev
, out
 R dT
298.15
R
 rev a1i rev 2 2

a 2i
 a 0i Tout  298.15   2 Tout   298.15  3   Trev
out 
3
 
 298.153  (46)
 R

 4

 out  5
 out  
 a3i T rev 4  298.154  a 4i T rev 5  298.155
  


Tin  Cpi  dT  P 
Si ,in  R     R ln  in 
298.15
 R  T P
 ref 
  Tin  a 2i 2 2 
 a 0i ln  298.15   a1i Tin  298.15   2 Tin  298.15   (47)
 
 R
 a3 a 4  P  
 i Tin3  298.153   i Tin4  298.154   ln  in  
 3 4 P  
 ref 

Dr. Medardo Serna González – FIQ - UMSNH 29

Diseño de Procesos Químicos Capítulo 3

 Cpi  dT
rev
Tout  Pout 
Sirev
, out
 R    R ln  
298.15
 R  T  Pref 


a 0
 i  ln
 Toutrev 
 298.15 

a2
2

 a1i Toutrev  298.15   i Toutrev   298.152
2
 (48)
 R
 
 3
 3
 4
 4

 a3i Toutrev   298.153  a 4i Toutrev   298.154  ln  Pout 
P 
 ref 



En estas funciones se observa que las condiciones de entrada y la presión de salida

del gas son las mismas tanto para el compresor real como para el isoentrópico. La temperatura
de la corriente de salida es la que depende de si la ruta de comprensión es adiabática e
irreversible (Tout) o adiabática y reversible ( Toutrev ).

Algoritmo de Cálculo
(a) Calcular los flujos molares de los compuestos en la corriente de salida de las Ecs.
(36) y (37).
(b) Calcular las entropías específicas de los compuestos en la corriente de entrada de
la Ec. (47).
rev
(c) Calcular la temperatura reversible de la corriente de salida, Tout , de la ecuación
resultante de la combinación de las Ecs. (40) y (48).
(d) Calcular las entalpías específicas de los compuestos en la corriente de entrada de
la Ec. (44).
(e) Calcular la entalpía total de la corriente de entrada de la Ec. (39).
(f) Calcular las entalpías específicas reversibles de los compuestos en la corriente de
salida de la Ec. (46).
(g) Calcular la entalpía total de la corriente de salida para la compresión isoentrópica,
rev
H out , de la Ec. (43).
(h) Calcular el trabajo del compresor isoentrópico, Wrev , de la Ec. (42).
(i) Calcular el trabajo del compresor real, Wrev , de la Ec. (41).
(j) Calcular la temperatura real de la corriente de salida, Tout , de la ecuación resultante
de la combinación de las Ecs. (37), (38) y (45).

Comentarios
A medida que se mejora el diseño de un compresor decrece su irreversibilidad y se
incrementa su eficiencia, es decir, el comportamiento del compresor real se aproxima más a
la compresión isoentrópica ideal. Los compresores bien diseñados tienen eficiencias
isoentrópicas en el rango de 0.80 a 0.90.
Las turbinas operan de manera inversa a los compresores, disminuyendo la presión
de un gas mientras generan trabajo, esto es: Pin > Pout. La eficiencia de una turbina es definida
como

Dr. Medardo Serna González – FIQ - UMSNH 30

Diseño de Procesos Químicos Capítulo 3

T 
We

 H out - H in  (49)
Wrev  H out
rev
- H in 

donde  T es la eficiencia isoentrópica de la turbina.

El modelo matemático del compresor se aplica a una turbina si se reemplaza la Ec.
(41) por la Ec. (49).

Listado de La función Compresor_Ideal_HDA.m

function Compressor_Ideal_HDA

% Diseño termodinámico del compresor del proceso HDA

% La corriente de recirculación gaseosa se considera como un gas ideal
% Cps tomados del Apéndice A de Poling y col. (2001)
% MSG | FIQ | UMSNH - 28 de Septiembre de 2020

clear, clc

% -----------------------------------------------------------------------
% DATOS - Variables Independientes y Propiedades Termodinámicas
% -----------------------------------------------------------------------
x = 0.5; % Conversión por paso del Tolueno en el reactor
YPH = 0.4; % Fracción molar del Hidrógeno en la purga
MR = 5; % Razón molar H2/Tolueno en la entrada del reactor
RG = 1.2779e+03; % Flujo de la recirculación gaseosa en kmol/h

Pin_Comp = 465; % Presión de entrada en psia

Pout_Comp = 550; % Presión de salida en psia
Tin_Comp = 310.9278; % Temperatura de entrada en K (100 °F)
Efic_Comp = 0.8; % Eficiencia del compresor

Tref = 298.15; % Tref en K

Pref = 14.696; % Pref en psia

% Indices asignados a los compuestos

% Hidrógeno - 1
% Metano - 2

% Capacidades caloríficas Cp/R - Poling et al (2001)

CpG = [2.883 3.681e-3 -0.772e-5 0.692e-8 -0.213e-11;
4.568 -8.975e-3 3.631e-5 -3.407e-8 1.091e-11];

% -----------------------------------------------------------------------
% CÁLCULO DE LA POTENCIA DEL COMPRESOR PARA UN GAS IDEAL
% -----------------------------------------------------------------------
% Flujos molares de los compuestos en la entrada y salida del compresor
Fin_Comp = [YPH*RG (1 - YPH)*RG]; % Fi,in en kmol/h
Fout_Comp = Fin_Comp; % Fi,out en kmol/h - Ecs. (22) y (23)

% Entropías de los compuestos en la entrada - Ec. (34)

Dr. Medardo Serna González – FIQ - UMSNH 31

Diseño de Procesos Químicos Capítulo 3

Sin = [Scalc(CpG(1,:),Tref,Tin_Comp,Pref,Pin_Comp) ...

Scalc(CpG(2,:),Tref,Tin_Comp,Pref,Pin_Comp)]; % S en kJ/(kmol K)

% Temperatura de salida del compresor reversible - Ec. (27)

Trev_s = Tin_Comp + 20; % Estimado de Tout reversible, K
[Trev] = fzero(@SB_Function,Trev_s); % Trev del compresor, K

% Trabajo reversible y trabajo real

Hin_Comp = [Hcalc(CpG(1,:),Tref,Tin_Comp) ...
Hcalc(CpG(2,:),Tref,Tin_Comp)]; % Hi,in en kJ/kmol - Ec. (31)
Hout_Rev = [Hcalc(CpG(1,:),Tref,Trev) ...
Hcalc(CpG(2,:),Tref,Trev)]; % Hi,out en kJ/kmol - Ec. (33)
Wrev = Fout_Comp*Hout_Rev' - Fin_Comp*Hin_Comp'; % Wrev en kJ/h, Ec.
(29)
Wreal = Wrev/Efic_Comp; % Wreal en kJ/h, Ec. (28)
We_hp = Wreal/3600*1000*1.341e-3; % Wreal en hp

% Temperatura de salida del compresor real

[Tout_real] = fzero(@EB_Compressor,Trev); % Temperatura en K -
Ecs. (24), (25) y (32)

% -----------------------------------------------------------------------
% VISUALIZACIÓN DE LOS RESULTADOS EN LA PANTALLA
% -----------------------------------------------------------------------
disp('Valores Especificados para las Variables Independientes')
fprintf(' Conversión por paso del Tolueno en el reactor = %3.2f\n',x)
fprintf(' Fracción molar del Hidrógeno en la recirculación =
%3.2f\n',YPH)
fprintf(' Razón molar H2/Tolueno en la alimentación al reactor =
%3.2f\n',MR)
fprintf(' Flujo molar de la recirculación gaseosa = %7.3f kmol/h
\n',RG)
disp(' ')
fprintf(' Presión de entrada = %5.2f psia\n',Pin_Comp)
fprintf(' Presión de salida = %5.2f psia\n',Pout_Comp)
fprintf(' Temperatura de entrada = %5.2f K\n',Tin_Comp)
fprintf(' Eficiencia del compresor = %3.2f\n',Efic_Comp)
disp(' ')
disp('Valores Calculados de las Variables Dependientes del Compresor')
fprintf(' Temperatura de salida del compresor reversible = %6.3f
K\n',Trev)
fprintf(' Trabajo de eje reversible = %6.3f kJ/h\n',Wrev)
fprintf(' Trabajo de eje real = %6.3f kJ/h\n',Wreal)
fprintf(' Trabajo de eje real = %6.3f hp\n',We_hp)
fprintf(' Temperatura de salida del compresor real = %6.3f
K\n',Tout_real)

%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
function [fTrev] = SB_Function(Trev)

Sout = [Scalc(CpG(1,:),Tref,Trev,Pref,Pout_Comp) ...

Scalc(CpG(2,:),Tref,Trev,Pref,Pout_Comp)]; % Sout en kJ/(kmol
K) - Ec. (35)

fTrev = Fout_Comp*Sout' - Fin_Comp*Sin'; % Ec. (27)

Dr. Medardo Serna González – FIQ - UMSNH 32

Diseño de Procesos Químicos Capítulo 3

end % SB_Function

%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
function [fTreal] = EB_Compressor(Treal)

Hout_Comp = [Hcalc(CpG(1,:),Tref,Treal) Hcalc(CpG(2,:),Tref,Treal)];

% kJ/kmol - Ec. (32)

fTreal = Fout_Comp*Hout_Comp' - Fin_Comp*Hin_Comp' - Wreal; % Ec.

(24)

end % EB_Compressor

%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
function [H] = Hcalc(Cp,T1,T2)
% This function calculates enthalpy change given Tl, T2 and the Cp
% constants. The units of R determine the units of H.
% The form of the polynomial as given by Poling et al (2001).
% Cp/R = Cp(1) + Cp(2)*T + Cp(3)*T^2 + Cp(4)*T^3 + Cp(5)*T^4;
% Temperature in K
R = 8.3143; % kJ/kmol K
H = Cp(1)*(T2 - T1) + Cp(2)/2*(T2^2-T1^2) + Cp(3)/3*(T2^3-T1^3) ...
+ Cp(4)/4*(T2^4-T1^4) + Cp(5)/5*(T2^5-T1^5);
H = R*H;
end % function Hcalc

%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
function [S] = Scalc(Cp,T1,T2,P1,P2)
% This function calculates entropy change given Tl, T2, P1, P2 and
the Cp
% constants. The units of R determine the units of S.
% The form of the polynomial as given by Poling et al (2001).
% Cp/R = Cp(1) + Cp(2)*T + Cp(3)*T^2 + Cp(4)*T^3 + Cp(5)*T^4;
% Temperature in K
R = 8.3143; % kJ/kmol K
S = Cp(1)*log(T2/T1) + Cp(2)*(T2 - T1) + Cp(3)/2*(T2^2 - T1^2) ...
+ Cp(4)/3*(T2^3 - T1^3) + Cp(5)/4*(T2^4 - T1^4);
S = R*S - R*log(P2/P1); % kJ/kmol
end % function Scalc

end % function Compressor_Ideal_HDA.m

Dr. Medardo Serna González – FIQ - UMSNH 33