TEMA 1: BALANCES MOLARES

1.1 VELOCIDAD DE REACCIÓN

La velocidad de reacción indica con que rapidez se consume una cierta cantidad de moles de
una especie química para formar otra especie química.

Como ya sabéis, en una reacción química la masa total no varía, tan sólo se transforman unas
especies en otras. No obstante, si hacemos balances de materia por componentes, podemos
hablar de velocidad de formación o de consumo de un componente concreto.

Consideremos la reacción de producción del insecticida DDT (diclorodifenilo tricloroetano) a

partir de clorobenceno y cloral:

CCl3CHO + 2 C6H5Cl —> (C6H4Cl)2CHCCl3 + H2O

Que podemos representar como:

A + 2 B —> C + D

La velocidad de desaparición del reactivo A se expresa como –rA y siempre es una número
positivo (por ejemplo, –rA = 4 moles A / L·s).

Así pues, podemos definir la velocidad de reacción –rA como el número de moles de A que
reaccionan por unidad de tiempo y por unidad de volumen.

¿Cómo se define matemáticamente la velocidad de reacción química?

¿Qué os parecería definir matemáticamente la desaparición de A de la siguiente forma?

dc A
−r A = (1)
dt

Donde cA es la concentración de A y t es el tiempo.

Parece una buena opción pero en realidad es una definición válida sólo para unas
determinadas condiciones: un reactor discontinuo con volumen de reacción constante.

Eso significa que si aplicamos un balance de materia a cada tipo de reactor, obtendremos,
para cada caso, una expresión matemática diferente que relacione –r A con las condiciones
específicas de ese reactor.

Pero, por otro lado, siempre existe una expresión de -rA que es independiente del tipo de
reactor. Es la denominada ecuación de velocidad, por ejemplo:

−r A =k· c A o −r A =k· c A
2
Ambas expresiones, la obtenida con el balance de materia y la ecuación de velocidad son,
obviamente, correctas y se utilizará una u otra en función de lo que se requiera.

1.2 EQUACIÓN GENERAL DEL BALANCE DE MATERIA

Como recordareis, para un sistema de límites definidos, el balance de materia (en moles) para
la especie j, en cualquier instante de tiempo t, es:

Entrada−Salida+ Generación=Acumulación

[ ][ ][ ][ ]
Velocidad de
Caudal molar Caudal molar Velocidad de generación
acumulación
de jque entra de jque sale de j por reacción química
de j dentro
al sistema − del sistema + dentro del sistema =
del sistema
(moles
tiempo )moles
tiempo ( ) moles
tiempo ( )moles
tiempo ( )
Que se puede expresar como: Entrada−Salida+ Generación=Acumulación

dN j
F jE −F jS +G j= (2)
dt

donde, FjE es el caudal molar de j en la entrada, FjS es el caudal molar de j en la salida, Gj es la

velocidad de generación de j y Nj es el número de moles de j dentro del sistema a tiempo t.

Si todas las variables del sistema (temperatura, concentración de la especie, etc.) son
uniformes en todo el volumen del sistema, la velocidad de generación Gj es:

G j=r j ∙ V (3)

donde, V es el volumen útil del sistema, por ejemplo, en un reactor sería el volumen ocupado
por la mezcla reactiva.

Pero la velocidad rj puede variar con la posición dentro del volumen del sistema:

Entonces, la velocidad de generación en el volumen ∆V1 sería:

∆ G j 1=r j1 ∙ ∆ V 1 (4)
y en el volumen ∆V2 sería:

∆ G j 2=r j 2 ∙ ∆ V 2 (5)

Así, la velocidad de generación total sería:

M M
G j=∑ ∆ G ji =∑ r ji ∙ ∆ V i (6)
i=1 i=1

y considerando que M —> ∞ y ∆V —> 0, se transforma en una integral:

V
G j=∫ r j ∙ dV (7)
0

Finalmente, el balance de materia para la especie j se puede escribir como:

V
dN j
F jE −F jS +∫ r j ∙ dV = (8)
0 dt

A partir de este balance se pueden desarrollar las ecuaciones de diseño de los diferentes
reactores continuos y discontinuos. Con estas ecuaciones se puede determinar el volumen del
reactor o el tiempo necesario para conseguir convertir una determinada cantidad de reactivos
en productos.

Entonces, la velocidad de generación en el volumen ∆V1 sería:

∆ G j 1=r j1 ∙ ∆ V 1

y en el volumen ∆V2 sería:

∆ G j 2=r j 2 ∙ ∆ V 2

1.3 REACTORES DISCONTINUOS (BATCH REACTORS)

¿En qué casos se utilizan estos reactores?

ü En operaciones a pequeña escala, como por ejemplo cuando se prueban nuevos

procesos para evaluar su viabilidad a escala industrial.

ü Para fabricar productos costosos.

ü Para procesos que son difíciles de desarrollar en continuo, como procesos con cambio
continuado del producto producido.

¿Qué ventajas aporta la operación en discontinuo?

ü La principal es que se puede conseguir una gran conversión dejando que los reactivos
reaccionen durante mucho tiempo.
¿Qué desventajas presenta la operación en discontinuo?

ü Requiere un elevado coste en mano de obra (carga, descarga, limpieza).

ü Puede presentar variabilidad en la calidad del producto de un lote a otro.

ü Es muy difícil producir a gran escala trabajando en discontinuo.

Un reactor discontinuo no tiene flujo de entrada de reactivos ni flujo de salida de productos

mientras se está produciendo la reacción:

F jE =F jS =0 (1)

Y, por lo tanto, el balance de materia (en moles) de la especie j es:

dN j V
=∫ r j ∙ dV (2)
dt 0

Si el reactor discontinuo está perfectamente agitado, la velocidad de reacción es igual en

cualquier punto del volumen y la anterior ecuación se transforma en:

dN j
=r j ∙V (3)
dt

Esta es la ecuación de diseño de uno de los reactores ideales, el REACTOR DISCONTINUO DE

TANQUE AGITADO (RDTA). Este reactor es ideal porque se supone que está perfectamente
agitado.

El RDTA se presenta con un esquema como el siguiente:

Consideremos la isomerización de la especie A en un RDTA:

A —> B

¿Qué tiempo se necesitaría para reducir el número inicial de moles de A, NA0, hasta un valor
deseado NA?

Partimos de la ecuación de diseño obtenida a partir del balance de materia:

dN A
=r A ∙V (4)
dt

reordenando,

dN A
dt= (5)
rA∙V

Integrando entre los siguientes límites:

• t = 0; NA = NA0
 • t = tf; NA = NAf
tf N Af
dN A
∫ dt = ∫ r (6)
0 N A0 A∙V

Resolvemos la integral de la izquierda y cambiamos el orden de los límites de integración de la

integral de la derecha:
NA0
dN A
tf =∫ (7)
N Af
−r A ∙V

1.4 REACTORES DE FLUJO CONTINUO

Los reactores de flujo continuo operan, casi siempre, en estado estacionario.

Existen dos tipos diferentes:

q El reactor continuo de mezcla completa

q El reactor tubular

Estos reactores son de uso habitual en procesos industriales, el de mezcla completa con
mezclas reactivas líquidas y el tubular con mezclas reactivas gaseosas.

En la fotografía podéis observar un reactor continuo de mezcla completa industrial:

En la fotografía podéis observar un reactor tubular industrial:

El reactor continuo de mezcla completa o reactor continuo de tanque agitado (RCTA) es el
segundo tipo de reactor ideal. Se asume que su agitación produce una mezcla perfecta y que
este tipo de reactor opera en estado estacionario.

Se representa con un esquema como el siguiente:

En un RCTA, la temperatura, las concentraciones de reactivos y productos y la velocidad de

reacción no dependen del tiempo ni de la posición espacial.

Eso significa que todas las variables del sistema (T, c j, rj) son iguales en todos los puntos
interiores del reactor y, por lo tanto, también son iguales al valor en la salida del reactor.

Todo esto cambiaría si la mezcla no fuese perfecta, algo que puede ocurrir en determinados
reactores industriales.

Vamos ahora a obtener la ecuación de diseño del RCTA partiendo del balance de materia:
V
dN j
F jE −F jS +∫ r j ∙ dV = (8)
0 dt

Como consideramos que opera en estado estacionario (EE):

dN j
=0 (9)
dt

También consideremos que la mezcla es completa, por lo tanto:

V

∫ r j ∙dV =r j ∙ V (10)
0

Por consiguiente, combinando las anteriores ecuaciones, obtenemos:

F jE−F jS
V= (11)
−r j

Esta es la ecuación de diseño de un RCTA.

Cuando el caudal volumétrico es constante, como ocurre con las reacciones en fase líquida y
en determinados casos en fase gas, los caudales molares pueden expresarse como:

F j=c j ∙ Q (12)

Donde Fj es el caudal molar del componente j, cj es la concentración del componente j y Q es el

caudal volumétrico.

Por consiguiente, la ecuación de diseño del RCTA puede reescribirse como:

c jE ∙ Q E −c jS ∙Q S Q ∙ ( c jE −c jS )
V= = (13)
−r j −r j
1.4.2. Reactor tubular

Un reactor tubular es un tubo cilíndrico en el que los reactivos se consumen de forma continua
a medida que circulan a lo largo del reactor. Como el anterior reactor continuo, el reactor
tubular también opera normalmente en EE.

Se representa con un esquema como el siguiente:

En un reactor tubular las concentraciones de las diferentes especies pueden variar a lo largo
del eje axial del reactor y también pueden variar radialmente.

No obstante, cuando la circulación del fluido es turbulenta, la variación de las concentraciones

puede ocurrir sólo a lo largo del eje axial. En ese caso, el reactor se comporta de forma ideal y
se denomina reactor continuo de flujo pistón (RCFP).

En un RCFP, las concentraciones de reactivos y productos varían a lo largo del eje axial y, por
lo tanto, también varía la velocidad de reacción.

Para deducir la ecuación de diseño de un RCFP vamos a hacer un balance de materia para la
especie j en un volumen muy pequeño del reactor.

∆V

F F
j ∆G j
j
V V+∆V

En un volumen tan pequeño, ∆V, podemos suponer que la variación de la velocidad de

reacción a lo largo del eje axial será despreciable. Por lo tanto:
V
∆ G j=∫ r j ∙ dV =r j ∙ ∆ V (14)
0

Si hacemos un balance molar en ∆V tendremos:

Entrada−Salida+ Generación=Acumulación

][ ][ ]
Velocidad de Velocidad de
[ ][
Caudal molar − Caudal molar +
de la especie j de la especie j
generación
de j dentro de ∆ V
= acumulación
de jdentro de ∆ V

F j|V − F j|V +∆V + r j ∙ ∆ V =0 (15)

Reordenando,

F j|V −F j|V +∆V

r j= (16)
∆V

Si consideramos que ∆V —> 0 y aplicamos límites:

F j|V −F j|V +∆V dF j

r j= lim =¿ ¿ (17)
∆V →0 ∆V dV

Obtenemos:

dF j
r j= (18)
dV

que es la ecuación de diseño de un RCFP.

¿Cómo podemos calcular el volumen de RCFP necesario para reducir el caudal molar del
reactivo A desde FAE hasta FAS?

Reordenamos la ecuación de diseño:

dF A
dV = (19)
rA

e integramos entre los siguientes límites:

• V = 0; FA = FAE

• V = V; FA = FAS
F AS F AE

( )
dF A −1
V =∫ =¿ ∫ ¿ dF A (20)
F AE
rA F
rA
AS
TEMA 2: DISEÑO DE REACTORES ISOTERMOS

2.1 DEFINICIÓN DE CONVERSIÓN

Para definir una conversión debemos escoger como base de cálculo una de las especies que
reaccionan y, después, relacionar el resto de especies con esa base de cálculo.

Consideremos una reacción general:

ν A ∙ A+ ν B ∙ B — > ν C ∙ C+ ν D ∙ D

Si tomamos el reactivo A como base de cálculo y dividimos por su coeficiente estequiométrico:

νB νC νD
A+ ∙ B — > ∙C + ∙ D
νA νA νA

Cuántos moles de C se formarán por cada mol de A que se consume?

Una manera de responder a esa pregunta es definiendo una CONVERSIÓN. La conversión XA es

el número de moles de A que reaccionan por mol de A alimentado al sistema

moles de A que reaccionan

X A= (1)
moles de A alimentados

Dado que siempre definiremos la conversión respecto a la misma base de cálculo,

eliminaremos el subíndice A para simplificar: X =X A
 Ø En reacciones irreversibles, la conversión máxima que se puede lograr es 1.

Ø En reacciones reversibles, la conversión máxima que se puede lograr es la conversión

de equilibrio.

¿Cuántos moles de C se formarán por cada mol de A que se consume?

Una manera de responder a esa pregunta es definiendo una CONVERSIÓN. La conversión XA es

el número de moles de A que reaccionan por mol de A alimentado al sistema

moles de A que reaccionan

X A= (1)
moles de A alimentados

Dado que siempre definiremos la conversión respecto a la misma base de cálculo,

eliminaremos el subíndice A para simplificar: X =X A

Ø En reacciones irreversibles, la conversión máxima que se puede lograr es 1.

Ø En reacciones reversibles, la conversión máxima que se puede lograr es la conversión

de equilibrio.

¿Cómo se puede calcular la conversión máxima de equilibrio?

La variación de la constante de equilibrio con la temperatura viene dada por la ecuación de

Van’t Hoff:
T2 o
K ∆ HT
ln 2 =∫ ∙ dT (2)
K1 T R ∙ T2
1

o
Donde ∆ H T es la entalpía estándar de la reacción a la temperatura de operación, K 1 y K 2 son
las constantes de equilibrio a las temperaturas T 1 y T 2, respectivamente y R es las constante
de los gases ideales.

En general:
T
∆ H T =∆ H 0+ ∫ ∑ ( ν j ∙ C pj ) ∙ dT
0 0
(3)
298 j

0
Donde ∆ H 0 es la entalpía estándar de la reacción a una temperatura de referencia y C pj es la
capacidad calorífica del compuesto j.
0 0
En numerosas situaciones podemos asumir que ∆ H T ≈ ∆ H 0 o que no son muy diferentes y se
puede calcular un valor medio entre las temperaturas 1 y 2:
0 0
0∆ H1+ ∆ H 2
∆H =T
(4)
2

Así, la ecuación de Van’t Hoff queda:

 o

( )
K2 ∆ HT 1 1
ln = ∙ − (5)
K1 R T1 T2

Vamos a ver un ejemplo de calculo:

Queremos conocer la conversión de equilibrio entre 0 y 100ºC de la reacción A <=> B. La

reacción se realiza en fase líquida y las capacidades caloríficas de ambos componentes son
iguales y constantes con la temperatura.

Datos
0
∆ G298 = -3375 cal/mol
0
∆ H 298 = -18000 cal/mol

Como las capacidades caloríficas son iguales:

T T
+ ∫ ∑ ( ν j ∙ C pj ) ∙ dT =∆ H + ∫ ( C pA−C pB ) ∙ dT =∆ H 298
0 0 0 0
∆ H =∆ H
T 298 298 (6)
298 j 298

Por lo tanto, la variación de la constante de equilibrio con la temperatura será:

o
ln
K
K 298
=
∆ HT
R
∙
1
−
298 T
1
( ) (7)

0
El valor de K298 se obtiene a partir de ∆ G298 :

0
∆ G298 =−R ∙ T ∙ ln K 298 (8)

0
El valor de K298 se obtiene a partir de ∆ G298 :

0
∆ G298 =−R ∙ T ∙ ln K 298

−3375=−1.98 ∙298 ∙ ln K 298

( )
cal
K 298 =exp ∙ (
3375
1.98 ∙298
=300
mol
cal
∙K
)
mol ∙ K

Y de la ecuación de Van’t Hoff:

ln
K −18000
300
=
1.99
∙
1
−
298 T
1
( ) (9)

Y de la ecuación de Van’t Hoff:

ln
K −18000
300
=
1.99
∙
1
−
298 T
1
( )
18000
lnK = −24.7 (10)
1.99∙ T

Por otro lado, la ecuación de equilibrio relaciona la constante de equilibrio con las
concentraciones de reactivos y productos en el equilibrio:

C Beq
K= (11)
C Aeq

Y estas concentraciones se pueden expresar en función de la conversión:

c B =c BE −c AE ∙ X ∙
( )
νB
νA
(12)

c A =c AE ∙ ( 1− X ) (13)

Así,

K=
C Beq
=
c BE−c AE ∙ X eq ∙ ( −11 ) = X eq (14)
C Aeq c AE ∙ ( 1−X eq ) 1−X eq

Si damos valores a la temperatura entre 0 y 100ºC, obtenemos la Keq de la ecuación (9) y

substituyendo ese valor en la ecuación (14) obtenemos la Xeq (ver valores en la siguiente tabla):

La representación gráfica nos muestra la dependencia de la Xeq con la temperatura. La reacción

que estamos estudiando es exotérmica y podemos observar que, a medida que aumenta la
temperatura, disminuye la conversión de equilibrio. Esta circunstancia ocurre siempre con
reacciones reversibles exotérmicas.

X eq
Si consideramos que la reacción es endotérmica, con el mismo valor del calor de reacción pero
en positivo, podemos observar que, a medida que aumenta la temperatura, aumenta la
conversión de equilibrio. Esta circunstancia ocurre siempre con reacciones reversibles
endotérmicas.

X eq

T (℃ )

2.2 ECUACIONES DE DISEÑO PARA REACTORES DISCONTINUOS

En los reactores discontinuos, a medida que transcurre el tiempo, se forman más productos.
Por tanto, en estos reactores, la conversión es función del tiempo.

Si N A 0 es el número de moles iniciales en el reactor, N A 0 ∙ X es el número de moles que

habrán reaccionado cuando haya transcurrido un tiempo t:

reaccionan ] [ alimentados ] [ moles de A alimentados ]

[ quemoles de A = moles de A ∙ moles de A que reaccionan =N ∙X
A0

También podemos expresar el número de moles de A que quedan en el reactor transcurrido

un tiempo t, N A , en función de N A 0 y X :

[ ][
moles de A
que quedan en el =
reactor a tiempo t
alimentados ][
moles de A − moles de A que
han reaccionado ]

N A =N A 0 −N A 0 ∙ X=N A 0 ∙ ( 1− X ) (1)
La anterior expresión la podemos generalizar para cualquier compuesto j (reactivo o producto)
implicado en la reacción:

νj
N j=N j0 −N A 0 ∙ ∙X (2)
νA

Cuando el reactor discontinuo tiene una mezcla perfecta, hablamos de un RDTA y, en ese caso,
la velocidad de reacción es la misma en cualquier punto del reactor a un determinado tiempo,
aunque dicha velocidad va variando a medida que transcurre el tiempo.

Sabemos que la ecuación de diseño de un RDTA es:

dN A
=r A ∙V (3)
dt

Si multiplicamos por (-1), podemos expresar la ecuación de diseño en función de la velocidad

de desaparición de A:

dN A
- =( −r A ) ∙ V (4)
dt

Si queremos expresar la ecuación de diseño en función de la conversión, podemos substituir

N A por la expresión de la ecuación 1:

d ( N A 0 ∙ ( 1−X ) )
- =( −r A ) ∙ V (5)
dt

Derivando, obtenemos:

dX
N A 0∙ =(−r A ) ∙ V (6)
dt

Esta es la ecuación de diseño de un RDTA en función de la conversión

Si el volumen de reacción es constante, V =V 0, podemos escribir la ecuación de diseño de la

siguiente manera:

−dN A 1
∙ =
−d
( )
NA
V0
=
−dc A V
=( −r A ) ∙ =( −r A )
(7)
dt V0 dt dt V0

Es decir,

−dc A
=(−r A ) (8)
dt

Si lo que nos interesa es conocer el tiempo necesario para alcanzar una determinada
conversión, lo podemos obtener integrando la ecuación de diseño:
 dX
dt=N A 0 ∙
−r A ∙V
X
dX
t=N A 0 ∙∫ (9)
0 −r A ∙V

En un reactor continuo, la conversión aumenta a medida que aumenta el tiempo que los
reactivos permanecen dentro del reactor, pero en este caso, ese tiempo es función del
volumen del reactor.

A este tiempo se le denomina TIEMPO DE RESIDENCIA y se define como:

V
τ= (10)
Q

Donde τ es el tiempo de residencia, V es el volumen útil del reactor y Q es el caudal

volumétrico que pasa a través del reactor.

Como mencionamos en una anterior ocasión, consideremos que los reactores continuos
operarán en EE por lo que, en este caso, no tiene sentido hablar de condiciones iniciales y
finales si no de condiciones de entrada y de salida.

Tampoco tiene sentido trabajar en número de moles Nj si no en caudales molares Fj

Así, podemos expresar el caudal molar en función de la conversión como:

F A=F AE ∙ ( 1− X ) (11)

νj
F j=F jE −F AE ∙ ∙ X (12)
νA

El caudal molar de entrada se puede calcular como:

F A E =Q ∙ c A E (13)

2.3.1. Reactor continuo de tanque agitado (RCTA)

Recordemos que, en un RCTA, los parámetros del sistema (temperatura, concentraciones,

velocidad de reacción) no varían, ni en el espacio (por la mezcla perfecta), ni en el tiempo
(operación en EE).

Así, la ecuación de diseño de este reactor es:

F AE−F A
V= (14)
−r A

Podemos expresar esta ecuación en función de la conversión substituyendo la ecuación 11:

F AE−( F AE ∙ ( 1− X ) ) F AE ∙ X
V= = (15)
−r A −r A
¿En qué condiciones debe evaluarse -rA en un RCTA?

Como es un reactor de mezcla perfecta, las condiciones en el interior del reactor son las
mismas que en la salida. Por lo tanto, (−r A ) es (−r A )Salida

2.3.2. Reactor continuo de flujo pistón (RCFP)

Recordemos que, en un RCFP, los parámetros del sistema (temperatura, concentraciones,

velocidad de reacción) varían a lo largo del eje axial del reactor pero no así en el tiempo
(operación en EE).

Así, la ecuación de diseño de este reactor es:

dF A
=r A (16)
dV

Podemos expresar esta ecuación en función de la conversión, multiplicando en primer lugar

por (-1) para tener la ecuación en función de (−r A ):

dF A
- =(−r A ) (17)
dV

Substituyendo ahora la ecuación 11:

d ( F AE ∙ (1−X ) )
- =( −r A ) (18)
dV

Derivando nos queda:

dX
F AE ∙ =(−r A ) (19)
dV

Y si el objetivo final es conocer el volumen del reactor, integramos la ecuación de diseño:

X
dX
V =F AE ∙∫ (20)
0 (−r A )
Para poder resolver cualquiera de las ecuaciones de diseño de los 3 tipos de reactores ideales
es necesario conocer como varía (−r A ) con la concentración, es decir, como varía (−r A ) con la
conversión.

2.4. APLICACIÓN DE LAS ECUACIONES DE DISEÑO PARA EL DISEÑO DE REACTORES

CONTINUOS

En este apartado veremos como dimensionar un RCTA y un RCFP conociendo la velocidad de

reacción en función de la conversión.
Escogemos un caso sencillo pero muy común, una ecuación de velocidad de primer orden:

(−r A )=k ∙c A =k ∙ c AE ∙ ( 1−X ) (1)

Donde k es la constante cinética, que sólo depende de la temperatura y c AE es la

concentración de A en la entrada del reactor.

Si recordamos las ecuaciones de diseño de un RCTA y un RCFP estudiadas en anteriores

apartados, podemos ver que el volumen del reactor es función, en ambos casos, de la inversa
1
de la velocidad de reacción
(−r A )
Esta inversa podemos expresarla como:

1 1 1
= ∙ (2)
(−r A ) k ∙ c AE ( 1−X )

En numerosas ocasiones, los datos cinéticos disponibles para el diseño de un reactor son datos
experimentales obtenidos en un laboratorio, como los del siguiente ejemplo:

Para la reacción: A—>B, se dispone de los siguientes datos cinéticos obtenidos a T = 300 K y P =
830 kPa. La mezcla reactiva utilizada era A puro y FAE era de 0.4 mol/s.

X -rA (mol/m3·s)
0 0.45
0.1 0.37
0.2 0.30
0.4 0.195
0.6 0.113
0.7 0.079
0.8 0.05

Si procesamos esos datos podemos obtener:

(FAE/-rA) (m3)
0.89
1.08
1.33
2.05
3.54
5.06
8

Estos datos se pueden representar

gráficamente y se obtiene la denominada
GRÁFICA DE LEVENSPIEL
Si analizamos la gráfica de Levenspiel podemos observar que, si la reacción se realiza a
temperatura constante, es decir, de forma isotérmica, la máxima velocidad de reacción se
obtiene con la mayor concentración de reactivo, es decir, cuando la conversión es casi 0. En
esas circunstancias, la inversa de la velocidad de reacción tendrá un valor pequeño.

Al contrario, si la concentración de reactivo es baja, la conversión será elevada y la velocidad

de reacción baja. En esas circunstancias, la inversa de la velocidad de reacción será elevada.

1
Para reacciones irreversibles, a medida que X—>1, (−r A )—>0 y, por lo tanto, —>∞
(−r A )
En consecuencia, y según las ecuaciones de diseño, V —> ∞

Esto implica que, aunque teóricamente es posible alcanzar una conversión del 100% (X=1), se
requeriría un volumen infinito para alcanzar ese resultado.

Para reacciones reversibles, la máxima conversión posible es la de equilibrio X MÁX =X eq

1
A medida que X—>Xeq, (−r A )—>0 y, por lo tanto, —>∞
(−r A )
En consecuencia, y según las ecuaciones de diseño, V —> ∞

Esto implica que, aunque teóricamente es posible alcanzar la conversión de equilibrio, se

requeriría un volumen infinito para alcanzar ese resultado.

2.4.1. Determinación del volumen de un RCTA

Utilizando los datos cinéticos que ya hemos visto anteriormente, se nos propone dimensionar
un RCTA para realizar la reacción A—>B de forma isotérmica a 300 K, tratando 0.4 mol/s de A
puro y alcanzando una conversión del 80%.

X -rA (mol/m3·s) (1/-rA) (m3·s/mol)

0 0.45 2.22
0.1 0.37 2.70
 0.2 0.30 3.33
0.4 0.195 5.13
0.6 0.113 8.85
0.7 0.079 12.7
0.8 0.05 20

El volumen del RCTA se puede calcular a partir de su ecuación de diseño:

F AE ∙ X
V= (3)
−r A

En un RCTA isotérmico, la temperatura y la conversión de la salida son idénticas a las del

interior del reactor. Por tanto, necesitamos el valor de la inversa de la velocidad a la
conversión de salida.

( )
3
1 m ∙s
=20
(−r A ) X =0.8
mol

Por tanto,

V=
F AE ∙ X
−r A (
= 0.4
mol
s ) (
∙ ( 0.8 ) ∙ 20
m3 ∙ s
mol )
=6.4 m3

Utilizando la gráfica de Levenspiel también podríamos calcular el volumen de forma gráfica. En

la siguiente figura, el volumen del RCTA corresponde al rectángulo gris:

2.4.2. Determinación del volumen de un RCFP

Vamos a hacer ahora los cálculos para realizar la misma operación en un RCFP.

Reacción A—>B de forma isotérmica a 300 K, tratando 0.4 mol/s de A puro y alcanzando una
conversión del 80%.

X -rA (mol/m3·s) (1/-rA) (m3·s/mol)

0 0.45 2.22
0.1 0.37 2.70
 0.2 0.30 3.33
0.4 0.195 5.13
0.6 0.113 8.85
0.7 0.079 12.7
0.8 0.05 20

La ecuación de diseño de un RCFP es:

X
dX
V =F AE ∙∫ (4)
0 (−r A )
Si no tenemos una expresión analítica de−r A y sólo disponemos de datos puntuales, debemos
resolver la integral de forma numérica, por ejemplo mediante el método de Simpson.

Considerando un intervalo de h = 0.2, el resultado será:

0.2
∙ [ 0.89+ 4 ∙ ( 1.33+3.54 ) +2 ∙ ( 2.05 ) +8.0 ] =2.17 m
3
V= (5)
3

También podríamos obtenir el volumen del RCFP utilizando la gráfica de Levenspiel. En la

siguiente figura, el volumen del RCFP corresponde al área en negro, ya que esa área
corresponde a la integral de la función:

2.4.3. Comparación entre un RCTA i un RCFP isotérmicos

Los resultados del cálculo del volumen de ambos reactores con las gráficas de Levenspiel son
extrapolables a cualquier pareja de RCTA y RCFP isotérmicos que realicen la misma operación.

En la siguiente figura se superponen las dos áreas que representan ambos volúmenes:
En una gráfica de Levenspiel, el volumen de un RCTA isotérmico siempre corresponderá a un
rectángulo mientras que el volumen de un RCFP isotérmico será el área bajo la función.

Si ambos reactores realizan la misma operación para alcanzar la misma conversión de salida, el
volumen del RCFP isotérmico siempre será menor al volumen del RCTA isotérmico.

2.5. REACTORES EN SERIE

Muchas veces, los reactores están conectados en serie, de manera que la salida de un reactor
es la entrada del siguiente. En esas condiciones es mejor definir la conversión en términos de
la ubicación a lo largo de la serie antes que definir una conversión individual para cada reactor.

Así, para reactores en serie como los de la figura:

Se puede definir la conversión en un punto i como:

moles totales de A que han reaccionado hasta el punto i

Xi=
moles de A alimentados en el primer reactor

Así,

F Ai =F AE −F AE ∙ X i=F AE ∙ ( 1− X i ) (1)

5.2.1 RCTA en serie

Consideremos 2 RCTA colocados en serie:

En el primer reactor, la velocidad de reacción es −r A 1 y la conversión en la salida es X 1

Si hacemos un balance de materia en ese reactor tendremos:

F AE−F A 1+ r A 1 ∙V 1=0 (2)

El caudal molar de A en la salida del primer reactor es:

F A 1=F AE−F AE ∙ X 1=F AE ∙ ( 1−X 1 ) (3)

Y combinando ambas expresiones, obtenemos el volumen del primer reactor:

V 1=F AE ∙ ( )
1
−r A 1
∙X1 (4)

En el segundo reactor, la velocidad de reacción es −r A 2 y la conversión de salida es X 2

Si hacemos un balance de materia en ese reactor tendremos:

F A 1−F A 2 +r A 2 ∙V 2=0 (5)

El caudal molar de A en la salida del segundo reactor es:

F A 2=F AE−F AE ∙ X 2=F AE ∙ ( 1−X 2 ) (6)

Hay que tener en cuenta que la ecuación 6 es válida porque hemos definido X 2 como una
conversión global del sistema de dos reactores y no como la conversión individual del segundo
reactor.

Combinando las ecuaciones 5 y 6, obtenemos el volumen del segundo reactor:

F A 1−F A 2 F AE ∙ ( 1− X 1 )−F AE ∙ ( 1−X 2 ) F AE

V 2= = = ∙ ( X 2 −X 1 ) (7)
−r A 2 −r A 2 −r A 2

Comparación de volúmenes para RCTA en serie

Si en el anterior esquema de dos reactores en serie se consigue un 40% de conversión en el

primer reactor, ¿qué volumen debe tener cada reactor para que la conversión global sea del
80%?

Datos

FAE = 0.4 mol/s

 X -rA (mol/m3·s)
0 0.45
0.1 0.37
0.2 0.30
0.4 0.195
0.6 0.113
0.7 0.079
0.8 0.05

F AE
En primer lugar, utilizamos los datos cinéticos para obtener
−r A

(FAE/-rA) (m3)
0.89
1.08
1.33
2.05
3.54
5.06
8

Primer reactor

Sabemos que X 1 = 0.4 y con la anterior tabla obtenemos: ( )

F AE
−r A X1
= 2.05 m3

Entonces, V 1= ( )
F AE
−r A 1 X1
∙ X 1=( 2.05 ) ∙ ( 0.4 )=0.82 m
3

Segundo reactor

Sabemos que X 2 = 0.8 y con la anterior tabla obtenemos:

( )
F AE
−r A X2
= 8.0 m3

Entonces, V 2=
( )
F AE
−r A 2 X2
∙ ( X 2 −X 1 )=( 8.0 ) ∙ ( 0.8−0.4 )=3.3 m
3

El volumen total de la serie de reactores es:

3
V =V 1+V 2=0.82+3.0=4.02 m

En la figura podemos ver la resolución gráfica de la serie de RCTA:

¿Qué pasaría si sólo hubiese un RCTA?

En ese caso, X S= 0.8 y con la anterior tabla obtenemos:

( )
F AE
−r A XS
= 8.0 m3

Entonces, V RCTA=
( )
F AE
−r A XS
∙ X S= ( 8.0 ) ∙ ( 0.8 )=6.4 m
3

Acabamos de comprobar que, para conseguir la misma conversión, el volumen global de una
serie de 2 RCTA es menor que un único RCTA.

¿Qué pasa si dividimos el volumen global en más RCTA?

En la figura podemos ver el resultado con una

serie de 5 RCTA

Si comparamos el resultado de una serie de 2 y otra de 5 reactores:

Podemos observar que cuanto mayor es el número de RCTA en serie, menor es el área global,
es decir, menor es el volumen global del sistema.

El límite sería una serie infinita de RCTA colocados en serie. En ese caso el volumen global sería
el área bajo la función cinética, es decir, la serie infinita de RCTA se comportaría como un
RCFP.
2.5.2. RCFP en serie

Ya hemos visto que colocar RCTA en serie disminuye el volumen total necesario para alcanzar
una determinada conversión.

¿Pasa lo mismo en el caso del RCFP?

Analicemos una serie con dos RCFP y el caso de un solo RCFP

Comparación de volúmenes para RCFP en serie

Primer reactor

El volumen del primer RCFP sería:

X1
dX
V 1=∫ F AE ∙ (8)
0
−r A

Segundo reactor

El volumen del segundo RCFP sería:

X2
dX
V 2=∫ F AE ∙ (9)
X1
−r A

Por lo tanto, el volumen total de la serie de reactores es:

 X1 X2
dX dX
V =∫ F AE ∙ +∫ F AE ∙ (10)
0
−r A X 1
−r A

Un único reactor

En el caso de un único reactor, el volumen del RCFP sería:

X2
dX
V RCFP=∫ F AE ∙ (11)
0
−r A

Se puede comprobar fácilmente que el volumen de la ecuación 10 y el de la ecuación 11 son

iguales. Es decir, en el caso de un RCFP, no se obtiene ninguna ventaja competitiva colocando
varios reactores en serie.

TEMA 3: REACTORES NO ISOTERMOS EN ESTADO ESTACIONARIO

3.1. BALANCE DE ENERGÍA EN UN REACTOR CONTINUO

Para poder diseñar reactores que operen en condiciones no-isotermas se necesita, además de la
ecuación de diseño, un balance de energía.

La razón es que, en condiciones no-isotermas, no se conoce la temperatura o temperaturas a la que

operará el reactor y, por lo tanto, tenemos dos incógnitas (V y T), por lo que necesitamos dos
ecuaciones.

El primer apartado de este tema es deducir el balance de energía para un reactor continuo. En
apartados posteriores, utilizaremos este balance de energía para diseñar reactores continuos no-
isotermos.

Vamos a deducir el balance de energía para un reactor continuo:

q w

FiE FiS

HiE HiS

Dades:

Fi: caudal molar del componente i

Hi: entalpia del componente i

q: flujo de calor

w: flujo de trabajo

¿
El calor de reacción a una determinada temperatura T, se puede expresar como:

°
∆ H R=∆ H R + ∆ C p ·(T −T R )

Y la expresión final del balance de energía en E.E: para un reactor en continuo es:

ws se puede menospreciar

n
F iE
q−w s−F AE · ∑ ·C pi · ( T −T E )−[ ∆ H °R + ∆ C p · ( T −T R ) · F AE · X ]=0
i=1 F AE

En el sumatorio se tiene que incluir todo lo que entra al reactor (reactivos i inertes), pero en el cálculo
de ∆ C psolo los reactivos:

∆ C p=∑ υ j · C pj
j

3.2. OPERACIÓ ADIABÀTICA

Muchas reacciones a escala industrial se realizan en condiciones adiabáticas, esto es, aislando
el reactor de su entorno, de forma que no se produce ningún intercambio de calor entre el
interior del reactor y su entorno (q=0).

La operación adiabática es un caso particular de la operación no-isoterma pero, como es muy

habitual a escala industrial, la estudiaremos por separado.

3.2.1 Balance de energía adiabático

Partiremos del balance de energía global para un reactor continuo desarrollado en un

apartado anterior:
n
F iE
q−w S −F AE ∙ ∑ ∙ C pi ∙ ( T S −T E ) −[ ∆ H 0R +∆ C p ∙ ( T −T R ) ] ∙ F AE ∙ X =0 (1)
i=1 F AE

En un reactor en el que se pueda despreciar el trabajo (w S ≈ 0) y que opere de forma

adiabática (𝑞=0), el balance quedaría como:
n
F
∑ F iE ∙C pi ∙ ( T −T E )
i=1 AE (2)
X=
−[ ∆ H +∆ C p ∙ ( T −T R ) ]
0
R

La ecuación 2 recibe el nombre de ADIABÁTICA DE CONVERSIÓN y nos permite relacionar

temperatura y conversión en un reactor adiabático.
0
Muchas veces se cumple que ∆ H R ≫ ∆ C p ∙ ( T −T R ), por lo que la adiabática de conversión se
transforma en:
n
F
∑ F iE ∙C pi ∙ ( T −T E ) (3)
X = i =1 AE

−[ ∆ H 0
R ]
Y la relación entre T y X en el reactor adiabático se vuelve lineal.

3.2.2. RCFP adiabático

Podemos reescribir la adiabática de conversión de la siguiente manera:

n
F iE
X ∙ [−∆ H 0R ] +∑ ∙ C ∙T + X ∙ ∆ C p ∙T R
i=1 F AE pi E
T= n
(4)
F iE
∑ ∙C + X ∙ ∆ C p
F AE pi
i=1
Para obtener los perfiles de temperatura, conversión y concentraciones a los largo del reactor,
necesitamos acoplar la adiabática de conversión con la ecuación de diseño del RCFP:

dX
F AE ∙ =−r A (5)
dV

Veamos cuál es el algoritmo de resolución:

Consideremos la siguiente reacción elemental y reversible en fase gas, con un alimento que
sólo contiene el reactivo A. Las Cp de ambos componentes son iguales:

A≤¿ B

Ecuación de diseño

dX −r A
= (6)
dV F AE

Ecuación de velocidad

(
−r A =k ∙ c A −
cB
K ) (7)

Ecuación de Arrhenius

k =k R ∙exp
[ (E 1 1
∙ −
R TR T )] (8)

Ecuación de Van’t Hoff

[ )]
0

(
∆ HR 1 1
K=K R ∙ exp ∙ − (9)
R TR T

Estequiometría (fase gas, ε =0, P=P E )

TE
c A =c AE ∙ ( 1− X ) ∙ (10)
T

TE
c B =c AE ∙ X ∙ (11)
T

Combinando las ecuaciones 7, 10 y 11:

(
−r A =k ∙ c AE ( 1− X )− )
X TE
∙
K T
(12)

Adiabática de conversión
 n
F iE
X ∙ [−∆ H 0R ] +∑ ∙ C ∙T + X ∙ ∆ C p ∙T R
i=1 F AE pi E
T= n
(13)
F iE
∑ ∙C + X ∙ ∆ C p
F AE pi
i=1

Como A es puro y ∆ C p=0:

X ∙ [−∆ H 0R ]
T =T E + (14)
C pA

Resolución numérica

En primer lugar, debemos recordar que lo que tenemos que hacer para encontrar el volumen
del RCFP es integrar su ecuación de diseño:
X1
F AE
V =∫ ∙ dX (15)
0
−r A

Como −r A =f ( X ,T ), esta integral sólo puede solucionarse de forma numérica:

a) Consideramos una X=0

b) Calculamos T con la ecuación 14

c) Calculamos k con la ecuación 8

d) Calculamos K con la ecuación 9

e) Calculamos −r A con la ecuación 12

F AE
f) Calculamos
−r A

g) Escogemos otra X i +1=X i + ∆ X y volvemos al paso b)

F AE
h) Preparamos una tabla de valores X vs
−r A

i) Resolvemos numéricamente la ecuación 15

3.2.3. Isomerización en fase líquida del n-butano

Queremos isomerizar n-butano a isobutano en un RCFP adiabático:

n−C 4 H 10≤¿i−¿ C 4 H 10
El isobutano es un producto que se utiliza en la producción de aditivos para la gasolina. La
reacción se realiza en fase líquida y a alta presión, utilizando trazas de un catalizador líquido.
La cinética es de primer orden en ambos sentidos.

Se pretende tratar un caudal molar de 163 kmol/h de una mezcla que contiene un 90% de n-
butano y un 10% de isopentano. Dicha mezcla entra al reactor a una temperatura de 330K. La
conversión que se desea obtener es del 70%.

¿Qué volumen de reactor es necesario? ¿Cuál será la temperatura de salida del reactor?

Datos:

0 J kmol
∆ H R=−6900 n−butano c AE=9.3 3
mol m

J kJ
C p n−butano=141 ∙K E=65.7
mol mol

J
C p isobutano=141 ∙K K=3.03 ( a 333 K )
mol

J
C p isopentano=161 ∙K −1
k =31.1 h ( a 360 K )
mol

Ecuación de diseño:

dX
F AE ∙ =−r A
dV

Ecuación de velocidad:
'
−r A =k ∙ c A−k ∙ c B

K=
k
'
(
−r A =k ∙ c A −
cB
K )
k

Ecuación de Arrhenius:
 k =k R ∙exp
[ (
E 1 1
∙ −
R TR T )]
Ecuación de Van’t Hoff:

[ )]
0

(
∆ HR 1 1
K=K R ∙ exp ∙ −
R TR T

Estequiometría: (fase líquida, QE = QS)

c A =c AE ∙ ( 1− X )

c B =c AE ∙ X

Adiabática de conversión:

(−∆ H 0R ) ∙ X
T =T E + n
F
∑ F iE ∙ C pi
i=1 AE

Conversión de equilibrio:

cB c AE ∙ X eq K
K= = X eq=
c A c AE ∙ ( 1−X eq ) 1+ K

Siguiendo el procedimiento de resolución numérica explicado anteriormente, podemos

obtener la siguiente tabla, con la que podremos integrar la ecuación de diseño y, en
consecuencia, obtener el volumen del RCFP.

X T(K) K(h-1) K Xeq -rA(kmol/m3·h) FAE/(-rA) (m3)

0 330 4.22 3.1 0.76 39.2 3.74
0.1 334.3 5.75 3.0 0.75 46.3 3.17
0.2 338.7 7.76 2.9 0.74 52.8 2.78
0.3 343.0 10.47 2.82 0.74 57.8 2.54
0.4 347.3 14.02 2.73 0.73 58.6 2.50
0.5 351.7 18.52 2.65 0.726 53.6 2.74
0.6 356.0 24.27 2.57 0.72 37.7 3.88
0.65 358.1 27.74 2.54 0.718 24.5 5.99
0.7 360.3 31.67 2.5 0.715 6.2 23.29
0.6 0.7
F AE F AE
V =∫ ∙ dX + ∫ ∙ dX =¿ ¿
0 −r A 0.6 −r A

0.1
¿ ∙ ( 3.74+3.88+ 4 ∙ ( 3.17+2.54 +2.74 ) +2 ∙ ( 2.78+2.50 ) ) +¿
3
 0.05 3
∙ ( 3.88+23.29+ 4 ∙ ( 5.99 ) )=1.73+ 0.85=2.58 m
3

3.3. RCFP EN ESTAT ESTACIONARI AMB BESCANVI DE CALOR

En esta sección vamos a considerar un RCFP donde el calor se agrega o se retira a través de las
paredes cilíndricas del mismo.

Para poder elaborar un modelo del reactor, asumiremos que no hay gradientes radiales en el
mismo y que el flujo de calor a través de la pared en una pequeña fracción de volumen es
como sigue:

3.3.1. Deducción del balance de energía per un RCFP

Efectuamos un balance de energía en una fracción de volumen ∆V, considerando (w S ≈ 0):

n n
∆ q+ ∑ [ F i ∙ H i ] V −∑ [ F i ∙ H i ]V +∆ V =0 (1)
i=1 i =1

El flujo de calor al reactor, ∆q, se puede calcular con el coeficiente global de transferencia de
calor, U, el área de intercambio, ∆A y la diferencia entre la temperatura del entorno, T a y la
temperatura del reactor, T.

∆ q=U ∙ ∆ A ∙ ( T a−T )=U ∙ a ∙ ∆ V ∙ ( T a−T ) (2)

Donde a es el área de intercambio de calor por unidad de volumen de reactor.

Sustituyendo ∆q en el balance de energía y dividiendo por ∆V, tenemos:

n n

∑ [ F i ∙ H i ]V −∑ [ F i ∙ H i ]V +∆ V (3)
U ∙ a ∙ ( T a−T ) + i=1 i=1
=0
∆V

Si evaluamos ∆V en el límite, cuando ∆V—>0, obtendremos:

n
d ∑ ( Fi∙ Hi)
i=1 (4)
U ∙ a ∙ ( T a−T ) + =0
dV
Expandiendo el sumatorio:
n n
dF i dH i
U ∙ a ∙ ( T a−T ) + ∑ ∙ H i−¿ ∑ ∙ Fi=0 ¿ (5)
i=1 dV i=1 dV

Si hacemos un balance de materia para la especie i, tendremos:

dF i ϑi
=r i= ∙ r A (6)
dV ϑA

Por otro lado, si diferenciamos la ecuación de la entalpía

T
H i=H +∫ C pi ∙ dT
0
i (7)
TR

Tendremos:

dH i dT
=C pi ∙ (8)
dV dV

Sustituyendo las ecuaciones (6) y (8) en la ecuación (5), obtenemos:

n n
rA dT
U ∙ a ∙ ( T a−T ) + ∑ ϑ i ∙ H i ∙ −¿ ∑ F i ∙C pi ∙ =0 ¿ (9)
i=1 ϑA i =1 dV
n
Donde, ∑ ϑ i ∙ H i=∆ H R (10)
i=1

Y reordenando la ecuación (9), obtenemos la expresión:

−r A
∙ ∆ H R−U ∙ a ∙ ( T −T a )
dT ϑA
= n (11)
dV
∑ F i ∙ C pi
i=1

Donde:

−r A
∙ ∆ H R es el CALOR GENERADO
ϑA

U ∙ a ∙ ( T −T a ) es el CALOR RETIRADO

Podemos escribir la ecuación (11) en función de la conversión:

−r A
∙ ∆ H R−U ∙ a ∙ ( T−T a )
dT ϑA
= (12)

(∑ )
dV n
FiE
F AE ∙ ∙ C +∆ Cp ∙ X
i=1 F AE pi
Si resolvemos simultáneamente la ecuación de diseño (balance de materia):

dX −r A
= (13)
dV F AE

y el balance de energía que acabamos de deducir:

−r A
∙ ∆ H R−U ∙ a ∙ ( T−T a )
dT ϑA
= (12)

(∑ )
dV n
FiE
F AE ∙ ∙ C +∆ Cp ∙ X
i=1 F AE pi

podremos obtener los perfiles de conversión y temperatura en un RCFP no-isotermo, siempre

que la temperatura del entorno, Ta, sea constante.

3.3.2. Isomerización en fase líquida del n-butano

Queremos isomerizar n-butano a isobutano en un conjunto de RCFP no-isotermos:

n−C 4 H 10≤¿i−¿ C 4 H 10

El isobutano es un producto que se utiliza en la producción de aditivos para la gasolina. La

reacción se realiza en fase líquida y a alta presión, utilizando trazas de un catalizador líquido.
La cinética es de primer orden en ambos sentidos.

El conjunto de RCFP está formado por 10 reactores colocados en paralelo. El volumen de cada
reactor es de 5 m3. Los reactores se refrigeran por convección natural y la temperatura
ambiente de la zona donde están instalados es de 37ºC.

Se pretende tratar un total de 163 kmol/h de una mezcla que contiene un 90% de n-butano y
un 10% de isopentano. Dicha mezcla entra al reactor a una temperatura de 310K. La
conversión que se desea obtener es del 75%.

Si la temperatura supera los 325K, empiezan a ocurrir reacciones indeseadas. ¿Podrá llevarse a
cabo la reacción en las condiciones descritas?

Datos:

0 J kmol
∆ H R=−6900 n−butano c AE=9.3 3
mol m

J kJ
C p n−butano=141 ∙K E=65.7
mol mol

J
C p isobutano=141 ∙K K=3.03 ( a 333 K )
mol
 J
C p isopentano=161 ∙K −1
k =31.1 h ( a 360 K )
mol
3
U ∙ a=5000 kJ /h ∙ m ∙ K

Ecuación de diseño:

dX
F AE ∙ =−r A
dV

Ecuación de velocidad:
'
−r A =k ∙ c A−k ∙ c B

k
K= '
(
−r A =k ∙ c A −
cB
K )
k

Ecuación de Arrhenius:

k =k R ∙exp
[ (
E 1 1
∙ −
R TR T )]
Ecuación de Van’t Hoff:

[ )]
0

(
∆ HR 1 1
K=K R ∙ exp ∙ −
R TR T

Estequiometría: (fase líquida, QE = QS)

c A =c AE ∙ ( 1− X )
 c B =c AE ∙ X

Balance de energía:

−r A 0
∙ ∆ H R−U ∙ a ∙ ( T −T a )
dT ϑA
=

( )
dV F IE
F AE ∙ C pA + ∙C
F AE pI
 3
Se puede comprobar que a los 2.9 m de RCFP ya se alcanza la conversión deseada

Se puede comprobar que, en ningún punto del RCFP, se superan los 325K
Además, podemos plantearnos una serie de preguntas adicionales:

1. ¿Nos acercamos a la conversión de equilibrio?

2. ¿Cómo cambia la velocidad de reacción a lo largo del reactor?

Xeq

Se puede comprobar que la conversión en el

RCFP va aumentando hasta acercarse
bastante a la conversión de equilibrio a partir
X
de los 4.5 m3 del RCFP.

Representamos rate vs V porque rate = -r

A

Se puede comprobar que la velocidad de

reacción aumenta en los primeros 0.5 m3 del
RCFP. Después va disminuyendo hasta
hacerse casi 0 a partir de los 4.5 m3 del RCFP.

3.3.3. Balance de energía para el fluido refrigerante/calentador de un RCFP con

intercambiador de calor

El fluido que circule por el intercambiador de calor será un refrigerante si la reacción es

exotérmica y un calentador si es endotérmica. Si el caudal másico de líquido refrigerante o
calentador es suficientemente grande, será capaz de absorber o ceder todo el calor generado
en la reacción sin que cambie su temperatura (T a=cte ).

No obstante, si su caudal másico no es suficientemente grande, la temperatura del líquido

refrigerante o calentador variará a lo largo del eje axial del reactor (T a=f (V )).

El líquido refrigerante o calentador puede circular en paralelo o en contracorriente respecto de

la mezcla reactiva.

A continuación vamos a deducir el balance de energía para la circulación en paralelo en el caso

de realizar una reacción exotérmica.

El balance de energía para

el refrigerante en ∆V es:

][ ]
Calor agregado por
[ entra aV ][
Energía que − Energía que +
sale de V + ∆V
conducción a través =0
de la pared

[ mc ∙ H c ]V − [ mc ∙ H c ]V +∆ V +U ∙ a ∙ ∆ V ∙ ( T −T a ) (1)

Dividiendo por ∆V y analizando del límite cuando ∆V —> 0:

dH c
−mc ∙ +U ∙ a ∙ ( T −T a) =0 (2)
dV

La variación de la entalpía se puede calcular como:

dH c dT a
=C pc ∙ (3)
dV dV

El balance de energía para el refrigerante en ∆V es:

][ ]
Calor agregado por
[ entra aV ][
Energía que − Energía que +
sale de V + ∆V
conducción a través =0
de la pared

[ mc ∙ H c ]V − [ mc ∙ H c ]V +∆ V +U ∙ a ∙ ∆ V ∙ ( T −T a ) (1)
Dividiendo por ∆V y analizando del límite cuando ∆V —> 0:

dH c
−mc ∙ +U ∙ a ∙ ( T −T a) =0 (2)
dV

La variación de la entalpía se puede calcular como:

dH c dT a
=C pc ∙ (3)
dV dV

Y finalmente podemos conocer como varía la temperatura del refrigerante a lo largo del
reactor con la ecuación:

dT a U ∙ a ∙ ( T −T a )
= (4)
dV mc ∙C pc

Los perfiles típicos son:

Para el caso de la circulación en contracorriente, tendríamos el siguiente esquema:

Y el balance de energía en este caso quedaría:

dT a U ∙ a ∙ ( T a −T )
= (5)
dV mc ∙C pc

Observamos que la única diferencia entre ambos tipos de circulación del líquido refrigerante
(ecuaciones 4 y 5) está en el término de las temperaturas: ( T −T a ) y ( T a−T ).

No obstante, hay una diferencia importante a la hora de solucionar un problema de diseño de

un RCFP no-isotermo con circulación en paralelo o con circulación en contracorriente.
Cuando la circulación es en paralelo, la condición inicial para cada ecuación diferencial se
localiza en el mismo punto del reactor:

Cuando la circulación es en contracorriente, la condición inicial para cada ecuación diferencial

no se localiza en el mismo punto del reactor:

Por lo tanto en este caso hay que definir todas las condiciones iniciales en el mismo punto del
reactor, lo que significa que una de ellas, la temperatura de salida del refrigerante, no se
conoce y debe suponerse y hacer resolverse el sistema de ecuaciones con un proceso iterativo:

3.4. OPERACIÓ ADIABÁTICA CON UNA REACCIÓN REVERSIBLE

Sabemos que la máxima conversión que se puede alcanzar en una reacción reversible está
limitada por el equilibrio. La conversión de equilibrio es función de la constante de equilibrio y,
en consecuencia, de la temperatura.

Por ejemplo, para una reacción de primer orden en ambos sentidos:

A≤¿ B
'
r =k ∙ c A −k ∙ c B (1)

k
K= ' (2)
k

K
X eq= (3)
1+ K

En función de si la reacción es exotérmica o endotérmica, la relación de la conversión de

equilibrio con la temperatura es como sigue:

Cuando la reacción reversible se realiza de forma adiabática, el máximo valor de conversión

que se puede alcanzar lo marcará la intersección de la adiabática de conversión con la curva de
equilibrio.

Recordemos la ecuación de la adiabática de conversión, considerando que el calor de reacción

es constante en el rango de trabajo:
n
F
∑ F iE ∙C pi ∙ ( T −T E ) (4)
X = i =1 AE

−[ ∆ H 0
R ]
Como podéis comprobar, la conversión alcanzada en el reactor adiabático es función de la
temperatura.

Consideremos una reacción reversible exotérmica, en ese caso tendríamos:

Cuando la reacción reversible es altamente exotérmica, con un solo reactor adiabático
obtendríamos una conversión muy baja.

No obstante, la simplicidad de la operación adiabática hace que sea muy atractiva a escala
industrial por lo que se ha buscado una estrategia para aumentar la conversión obtenida.
Dicha estrategia consiste en colocar varios reactores adiabáticos en serie con enfriamiento de
la corriente de reactivos y productos entre ellos:
Los puntos donde se cruzan la adiabática de conversión y la curva de equilibrio son sólo puntos
teóricos de trabajo ya que en ese punto la velocidad de reacción sería cero. En realidad, la
salida de cada reactor puede encontrarse cercana a ese punto pero nunca exactamente en él.

En reacciones reversibles altamente endotérmicas nos encontraríamos con las mismas

limitaciones y la solución técnica sería la misma, pero los intercambiadores calentarían la
mezcla reactiva en lugar de enfriarla.

Un ejemplo industrial es el aumento de octanaje de la gasolina. Cuanto más compactas son las
moléculas de hidrocarburos, respecto de un mismo número de átomos de carbono, mayor es
el octanaje.

Por lo tanto, la reacción industrial consiste en transformar hidrocarburos de cadena lineal en

compuestos como naftalenos o aromáticos:

El rango de temperaturas en el que se pueden desarrollar esas reacciones es estrecho, por

encima de 530ºC se dan reacciones secundarias no deseadas y por debajo de 430ºC la
velocidad de reacción es demasiado lenta.

El esquema de reactores y temperaturas que se obtienen a escala industrial es el siguiente:

No todos los reactores adiabáticos utilizados son del mismo tamaño pero sus medidas están
entre 10-20 m de altura y 2-5 m de diámetro.