REACTORES

HOMOGENEOS
IDEALES
UNIDAD 1

S
! Definir la etapa de reacción química como base del diseño de reactores,
destacando la importancia de la cinética química, tanto en sistemas homogéneos
como heterogéneos, catalíticos o no.
! Formular la ecuación cinética en función de los factores que afectan a la
velocidad de una reacción química: concentración, temperatura y medio.
! Desarrollar los principales parámetros para utilizar la ecuación cinética, que ha
de basarse siempre en datos experimentales bajo condiciones controladas.
! Establecer los fundamentos del diseño de reactores, analizando la información
que es necesario poseer para proceder al mismo.
! Definir los tres modelos básicos para el estudio de los reactores químicos,
desarrollando en cada caso sus respectivas ecuaciones de diseño.
! Esquematizar los equipos utilizados para llevar a cabo las reacciones químicas,
tanto homogéneas como heterogéneas.
Ingeniería de la Reacción Química: Disciplina que sintetiza la información,
los conocimientos y la experiencia en una gran variedad de campos para
diseñar el equipo necesario en el que se puedan desarrollar las reacciones
químicas.

Reactor químico: Unidad de proceso diseñada para llevar a cabo una o

varias reacciones química o procesos en los que unos “reactivos” se
convierten en unos “productos”.

Cambios energéticos que se producen en las reacciones químicas

(Termodinámica Química):

! Efecto térmico durante la reacción (calor de reacción).

! Grado de extensión de la reacción (equilibrio químico).
Reacción en un tiempo aceptable y con un rendimiento adecuado (Cinética
Química):

! Velocidad de reacción.
! Variables que influyen sobre la velocidad de reacción.

Si interviene más de una fase (la velocidad global puede ser la velocidad
de transporte):

! Etapas físicas de transporte en las fases.

! Etapas físicas de transporte entre las fases.
Según el número y tipos de fases implicadas, las reacciones se consideran
como:
! Reacciones homogéneas: Transcurren solamente en una fase
(normalmente entre gases o entre líquidos miscibles entre sí).
! Reacciones heterogéneas: Requieren, al menos, la presencia de más
de una fase para transcurrir del modo en que lo hacen (normalmente,
bifásicas, fluido-fluido o fluido-sólido).

En ambos tipos puede influir la presencia de catalizadores.

Catalizador: Sustancia que cambia la velocidad de reacción (acelerándola

o retardándola) y que no es ni un reactivo ni un producto, ni necesita estar
presente en grandes cantidades; generalmente no sufre variación.
Ecuación cinética: Expresión de la velocidad de la reacción que contempla
la influencia de los factores que influyen sobre ella.

Velocidad de reacción: Cantidad (moles) de reactivo o producto transfor-

mado por unidad de volumen y unidad de tiempo:

La velocidad de reacción depende (comprobado experimentalmente) de:

! Composición.
! Temperatura.
! Otros factores del medio (reacciones heterogéneas).
El modelo más sencillo y útil de expresar la dependencia de la velocidad
con la concentración es:

! Exponentes: Órdenes de reacción (individuales o global).

! Parámetro: Coeficiente cinético o velocidad reacción

específica de (incluye la dependencia con
la temperatura).
Se ha encontrado experimentalmente una dependencia del coeficiente
cinético con la temperatura del tipo (Arrhenius, 1885):

! ko: “Factor de frecuencia”.

! E: “Energía de activación”.

Los parámetros de esta ecuación se determinan a partir de datos

experimentales a varias temperaturas.
En las reacciones heterogéneas se producen fenómenos de transferencia
de materia y transmisión de calor simultáneos a la propia reacción y que
pueden controlar la velocidad global.

Estos fenómenos siguen modelos diferentes a los de las reacciones y han

de obtenerse específicamente para cada caso concreto:

! Reacciones fluido-fluido.
! Reacciones fluido-sólido no catalíticas.
! Reacciones fluido-sólido catalíticas.

IngQui-5 [9]
Definiciones de uso generalizado en la Ingeniería de la Reacción Química:

! Reactivo limitante: Reactivo que se agota en primer lugar si la

reacción es completa; normalmente es el reactivo que interviene en
menor proporción con respecto al valor estequiométrico y es el que se
utiliza generalmente para aplicar la ecuación cinética.

! Conversión: También llamada “conversión fraccional”, es la fracción

de un reactivo (generalmente, el reactivo limitante) que se ha converti-
do en productos:
 A volumen
constante:

y la relación entre concentración y conversión será:

La ecuación cinética expresa la velocidad en función de concentraciones
de las sustancias o de factores relacionados con el transporte.

Objetivo del diseño: Determinar una ecuación cinética aplicable para

diseñar el reactor.

Formulación de la ecuación cinética: Conocimiento de los valores

numéricos de los parámetros (coeficiente cinético y órdenes de reacción)
determinados experimentalmente por:

! Ajuste a ecuaciones cinéticas enunciadas empíricamente.

! Ajuste a ecuaciones cinéticas emanadas de modelos de reacción o de
transporte-reacción.
Obtención de datos experimentales en condiciones adecuadas:

! Sistema de laboratorio de operación sencilla y fácil control.

! Seguimiento del avance de la reacción mediante la concentración de
alguna especie o propiedad proporcional a ella.

Ajuste de los datos experimentales a:

! Ecuaciones de velocidad basadas en mecanismos hipotéticos.

! Ecuaciones empíricas potenciales.

El análisis de los datos experimentales se hace mediante dos métodos

(considerando sistemas isotermos y a volumen constante):
El método integral es fácil de aplicar y está recomendado cuando:

" Se ensayan ecuaciones cinéticas sencillas.

" Los datos son escasos o están dispersos.

! Método diferencial: Diferenciar directamente los datos experimentales de

concentración frente al tiempo (obtención de la velocidad) y compararlos con la
ecuación cinética propuesta.

El método diferencial es muy útil para el análisis de comportamientos

complejos pero:

" Requiere mayor cantidad de datos.

Ejemplo
 Solución
DATOS METODO DIFERENCIAL A==>P
EXPERIMENTALES
0
t(s) Ca(mol/L) LN CA LN R 0 0.5 1 1.5 2 2.5
2.07944154 -2.302585
20 8 2 093 -1
1.79175946 -2.99573227
40 6 9 4
1.60943791 -3.401197
-2
60 5 2 382
y = 1.7215x - 6.0202
1.09861228 -4.094344 R² = 0.98696
120 3 9 562
LN R

0.69314718 -4.787491 -3
180 2 1 743

-4

-5
k=411.66 n=1.7215

-6
LN CA
 Solución
METODO INTEGRAL
350
y = 328.34x - 10.888
integral R² = 0.9692
300
t(s) Ca(mol/L) 1/CA t(s)
0 10 0.1 0
250
20 8 0.125 20
0.16666666
40 6 7 40 200
1/CA
60 5 0.2 60
0.33333333 150
120 3 3 120
180 2 0.5 180 100
300 1 1 300
50

0
k=328.34 n=2 0 0.2 0.4 0.6 0.8 1 1.2
t(s)
 EJERCICIO 2
124 Ingeniería de la cinética química

Proponga una ecuación de velocidad que concuerde con estos datos.

2-10. La velocidad de esterificación de ácido acético y alcohol se puede aumentar con HCl
como catalizador homogéneo. A 100 OC, la velocidad de la reacción directa es

r, = k2 C,C,, moles g/(L)(min)

k, = 4.76 x 10-’ L/(mol g)(min)

y la velocidad de la reacción inversa de

r; = k; CEC, moles g/(L)(min)

k; = 1.63 x 10-’ L/(mol g)(min)

donde C, = concentración del ácido acético

CoH = concentración del alcohol
CE = concentración del éster
C, = concentración del agua
Una mezcla inicial consta de masas iguales de 90% en peso de solución acuosa de ácido y 95%
en peso de solución de etanol. Calcule la conversión de ácido a éster para diversos tiempos de
reacción, en condiciones de volumen constante. Suponiendo una miscibilidad completa, esti-
me la conversión de equilibrio.
 2.0 PLANTEAMIENTO

1) PARA CALCULAR LA CONVERSIÓN A

DIFERENTES TIEMPOS DE REACCIÓN, ES
NECESARIO PLANTEAR Y RESOLVER LA
ECUACIÓN DE VELOCIDAD DE REACCIÓN
POR EL MÉTODO INTEGRAL

2) PARA CALCULAR LA CONVERSIÓN DE

EQUILIBRIO, HAY QUE RESOLVER EN LA
CONDICIÓN DE EQUILIBRIO
Diseño: Determinación del tamaño (volumen del reactor) para obtener una
determinada cantidad de producto bajo unas condiciones establecidas.

La información necesaria sobre el sistema de reacción se muestra en la

Figura siguiente:
: Información necesaria para el diseño de reactores químicos.
 EJERCICIO 2
124 Ingeniería de la cinética química

Proponga una ecuación de velocidad que concuerde con estos datos.

2-10. La velocidad de esterificación de ácido acético y alcohol se puede aumentar con HCl
como catalizador homogéneo. A 100 OC, la velocidad de la reacción directa es

r, = k2 C,C,, moles g/(L)(min)

k, = 4.76 x 10-’ L/(mol g)(min)

y la velocidad de la reacción inversa de

r; = k; CEC, moles g/(L)(min)

k; = 1.63 x 10-’ L/(mol g)(min)

donde C, = concentración del ácido acético

CoH = concentración del alcohol
CE = concentración del éster
C, = concentración del agua
Una mezcla inicial consta de masas iguales de 90% en peso de solución acuosa de ácido y 95%
en peso de solución de etanol. Calcule la conversión de ácido a éster para diversos tiempos de
reacción, en condiciones de volumen constante. Suponiendo una miscibilidad completa, esti-
me la conversión de equilibrio.
 1.0. TRADUCCION

CH3-COOH + CH3-CH2 – OH CH3 – COO – CH2 – CH3 + H2O

A+B C+D
A: ACIDO ACETICO Datos adicionales
B: ETANOL PMA= 60 g/gmol
C:ACETATO DE ETILO ρA=1.05 g/ml
D: AGUA
PMB= 46 g/gmol
R= k1CACB – k2CCCD gmol/L-min
k1=4.76 x 10-4 L/(gmol-min) ρB=0.0789 g/ml
k2=1,63 x 10-4 L/(gmol-min)
SOL. A= 90% EN PESO DE A
SOL. B = 95 % EN PESO DE B
MASA MEZCLA= MSOL A+MSOL B
MSOL A= MSOL B

1) CONVERSIÓN DEL ESTER A DIFERENTES TIEMPOS?

2) CONVERSIÓN DE EQUILIBRIO?
 2.0 PLANTEAMIENTO

1) PARA CALCULAR LA CONVERSIÓN A

DIFERENTES TIEMPOS DE REACCIÓN, ES
NECESARIO PLANTEAR Y RESOLVER LA
ECUACIÓN DE VELOCIDAD DE REACCIÓN
POR EL MÉTODO INTEGRAL

2) PARA CALCULAR LA CONVERSIÓN DE

EQUILIBRIO, HAY QUE RESOLVER EN LA
CONDICIÓN DE EQUILIBRIO
 3.0. CALCULOS
3.1.1. Expresión cinética
−dCa
= k1CaCb − k2CcCd
dt

3.1.2. Transformación a conversión

A + B = C + D

t=0 CA0 CB0 0 0

T=t in XCA0 XCA0 XCA0 XCA0

T=t fin CA=CA0- CB=CB0-XCA0 CC=XCA0 CD=XCA0

XCA0

dx
= k1(1− x)(Cb0 − xCa0) − k2Ca0x 2
dt
 3.1.2. Integrando
dx
t=∫ dx
k1(1− x)(Cb0 − xCa0) − k2Ca0x 2 t=∫
a(1− x)(b − cx) − ex 2
Sean:
k1=a
b=CB0 (2ac − 2e)x − ( a(ac 2 − 2abc + ab2 + 4eb) − a(c + b))
c=CA0 Ln
e=k2CA0
(2ac − 2e)x + ( a(ac 2 − 2abc + ab2 + 4eb) − a(c + b))
t=
a(ac 2 − 2abc + ab2 + 4eb)

(2k1Cao − 2k2Cao)x − ( k1(k1Cao 2 − 2k1CboCao + k1Cbo 2 + 4k2CaoCbo) − k1(Cao + Cbo))

Ln
(2k1Cao − 2k2Cao)x + ( k1(k1Cao 2 − 2k1CboCao + k1Cbo 2 + 4k2CaoCbo) − k1(Cao + Cbo))
t=
k1(k1Cao 2 − 2k1CboCao + k1Cbo 2 + 4k2CaoCbo)
 3.2. Conversión de equilibrio

3.2.1. Condición de equilibrio

R = k1CACB- K2CCCD=0
R = k1(1− x)(Cb0 − xCa0) − k2Ca0x 2

La conversión de equilibrio estará dada por:

(Cao + Cb0) ± (Cao + Cbo) − 4CaoCbo( k1− k2 )
xeq = k1
2Cao( k1− k2 )
k1
Reactores homogéneos ideales
 IDEAL??
Es una representación mental
que se relaciona con algo real

La expresión ideal pude indicar a un

prototipo, modelo o ejemplar de algo
 MODELO?
Aquello que se toma como referencia para
tratar de producir algo igual
Los modelos matemáticos son
sistemas de ecuaciones y
UNA DE TANTAS CLASIFICACIONES proposiciones lógicas que
MENTALES: Representación mental de estructuras y intentan representar las
procesos relaciones entre variables
SEMÁNTICOS: Formalización linguística de un modelo
(propiedades mensurables del
mental
FÍSICOS: Modelos a escala reducida construidos
sistema cuyas magnitudes
generalmente en el laboratorio con materiales de distinta varían en el tiempo) y
naturaleza; se asume una similitud dinámica entre el modelo parámetros (cantidades
y el sistema real temporalmente invariables que
ANALÓGICOS: Utilizan sistemas mecánicos o eléctricos caracterizan al sistema)
análogos a los sistemas bajo consideración, como por
ejemplo, el flujo de electricidad puede ser capaz de simular
el flujo de agua.
NUMÉRICOS: Representación formal en base al uso de
matemáticas.
PV=nRT
 Reactor tipo tanque
Aspectos mecánicos o
detalles de construcción
Reactores tubulares

Reactores continuos
Intercambio de materiales con Reactores discontinuos
su ambiente.
Reactores semicontinuos

Reactores homogéneos
Fases presentes
Reactores heterogéneos

CLASIFICACION DE LOS Reactores cataliticos

REACTORES QUIMICOS.
Presencia de catalizadores Reactores autocatalíticos

Reactores no catalíticos

Reactores isotérmicos

Reactores pseudoisotérmicos

Temperatura de operación. Reactores no isotérmicos

Reactores adiabáticos

Reactores especiales Reactores no adiabáticos

 TIPOS DE REACTORES
IDEALES
S Se definen tres tipos de reactores, combinando el modo de operación

y los modelos ideales de flujo:

!Reactor discontinuo
!Reactor tipo tanque agitado
!Reactor tubular

SSe analizan mediante ecuaciones de conservación y se obtienen sus ecuaciones de

diseño para reacciones homogéneas (para reacciones heterogéneas se combinan o
adaptan estos tres tipos básicos).
Por que reactores ideales?

Se llama mezcla ideal, a aquella en la que no

hay interacción entre sus partes y por tanto
cualquier magnitud conservativa y aditiva
será igual a la media ponderada de los valores
que correspondan a cada especie en su estado
de referencia(Propiedades molares parciales)

vid(T , p, x) = ∑ xivi(T , p)
 BALANCE DE MATERIA EN
UN ELEMENTO DE
VOLUMEN REACCIONANTE

q Q 0 : Flujo volumétrico que entra al

elemento de volumen[L3/tiempo]
CA0:Concentración inicial del reactivo
limitante(entrada del elemento de
QS, volumen)[M/L3]
Q0, CA, T0: Temperatura a la entrada del elemento
CA0, TS de volumen.[ºC,ºK]
T0 Q: Flujo volumétrico que sale del elemento
V: Volumen del elemento de de volumen[L3/tiempo]
reacción [L3] CA: Concentración final del reactivo
CA: Concentración del RL en el
limitante(salida del elemento de volumen)
interior del elemento[M/L3]
q: Calor transferido entre el [M/L3]
e l e m e n t o d e vo l u m e n y l o s T: Temperatura a la salida del elemento de
alrededores [Energía/tiempo] volumen.[ºC,ºK]
 Balance con reacción
para el reactivo limitante
ACUMULACION
Nº de moles de reactivo limitante
por unidad de tiempo que entran TRANSPORTE DE MASA
al elemento de volumen

Nº de moles de reactivo limitante

por unidad de tiempo que se
acumulan en elemento de
Nº de moles de reactivo volumen
limitante por unidad de
tiempo que salen del elemento
de volumen

Nº de moles de reactivo REACCION

limitante por unidad de
tiempo que se transforman por
reacción
 Balance con reacción
para el reactivo limitante

QCA0 QSCA VR d(VCA)/dt

Q0: Flujo volumétrico que entra al elemento de volumen[L3/tiempo]

CA0:Concentración inicial del reactivo limitante(entrada del elemento de volumen)[M/L3]
T0: Temperatura a la entrada del elemento de volumen.[ºC,ºK]
Q: Flujo volumétrico que sale del elemento de volumen[L3/tiempo]
CA: Concentración final del reactivo limitante(salida del elemento de volumen)[M/L3]
T: Temperatura a la salida del elemento de volumen.[ºC,ºK]
V: Volumen del elemento de reacción [L3]
CA: Concentración del RL en el interior del elemento[M/L3]
q: Calor transferido entre el elemento de volumen y los alrededores [Energía/tiempo]
R: Ecuación cinética, rapidez de reacción = f(T, C)
 BALANCE DE ENERGÍA EN
UN ELEMENTO DE
VOLUMEN REACCIONANTE

q Q 0 : Flujo volumétrico que entra al

elemento de volumen[L3/tiempo]
CA0:Concentración inicial del reactivo
limitante(entrada del elemento de
QS, volumen)[M/L3]
Q0, CA, T0: Temperatura a la entrada del elemento
CA0, TS de volumen.[ºC,ºK]
T0 Q: Flujo volumétrico que sale del elemento
V: Volumen del elemento de de volumen[L3/tiempo]
reacción [L3] CA: Concentración final del reactivo
CA: Concentración del RL en el
limitante(salida del elemento de volumen)
interior del elemento[M/L3]
q: Calor transferido entre el [M/L3]
e l e m e n t o d e vo l u m e n y l o s T: Temperatura a la salida del elemento de
alrededores [Energía/tiempo] volumen.[ºC,ºK]
 BALANCE DE ENERGÍA EN
UN ELEMENTO DE VOLUMEN
REACCIONANTE

ENERGÍA ENERGÍA

Capacidad de realizar trabajo,

de producir movimiento, de
generar cambio

POTENCIAL
CINÉTICA MOVIMIENTO
MECÁNICA

ELECTRICA
ENERGÍA
NUCLEAR
MAGNETICA
TÉRMICA
SONORA
SOLAR
 CONSERVACION DE
LA ENERGÍA
Para un sistema cerrado en un proceso
cíclico, el calor transferido es igual al trabajo
realizado por el sistema

!∫ dq = !∫ dw !∫ dq − !∫ dw = 0
DIFERENCIALES
INEXACTAS =
PROCESO: CONJUNTO PROPIEDADES DE
DE CAMBIOS DE ESTADO TRAYECTORIA
ESTADO: PUNTO EN EL QUE UN
SISTEMA TIENE SUS PROPIEDADES
DEFINIDAS(PVT)
 BALANCE DE ENERGÍA EN
UN ELEMENTO DE VOLUMEN
REACCIONANTE

m 2 2
W = (v2 − v1 )
2

! ! W = mg(h2 − h1 )
W= ∫ F • dS
1
W = k(x22 − x12 )
2

W= ∫ P dV
 CONSERVACION DE
LA ENERGÍA

Para un sistema cerrado en un proceso

cíclico, el calor transferido es igual al trabajo
realizado por el sistema

!∫ dq = !∫ dw !∫ dq − !∫ dw = 0
∫ dq − ∫ dw = ΔU PROPIEDAD DE
ESTADO
PROPIEDAD DE
TRAYECTORIA

q = ΔU
q− ∫ p dv = U 2
−U1 q = (U 2 + PV2 ) − (U1 + PV1 ) = ΔH
q − p(V2 −V1 ) = U 2 −U1
 CONSERVACION DE
LA ENERGÍA

q ENERGÍA QUE ENTRA AL ELEMENTO

• ENTALPIA DE LOS REACTIVOS (T0)
• ENERGÍA(CALOR)DE REACCIÓN

ENERGÍA QUE SALE DEL ELEMENTO

QS, • ENTALPIA DE LOS REACTIVOS (T)
CA, • ENERGIA (CALOR) TRANSFERIDO
Q0,
CA0, TS
T0
VARIACION DE ENERGÍA (ENTALPIA)
EN EL ELEMENTO DE VOLUMEN

dh
(haT 0 + hr ) − (haT + q) =
dt
 CONSERVACION DE
LA ENERGÍA
q

ρ: densidad (masa/volumen)
Cp: Capacidad calorífica a presión constante
QS, (Energía/masa grado)
Q0, CA, ΔH: Entalpía (Calor de reacción)
CA0, TS R: Velocidad de reacción
T0 T: Temperatura
T: tiempo

dT
Qρ0Cp0T0 − QρCpT + (−ΔH r )VR − q =VCp
dt
dT
QρCp(T0 − T ) + (−ΔH r )VR − q =VCp
dt
 Reactor discontinuo
por lotes , intermitente o batch
El reactor intermitente de tanque agitado conocido también, como Reactor Intermitente Reactor Batch, se usa
en investigaciones cinéticas dentro del laboratorio; sin embargo, también se le puede encontrar en operaciones
industriales.

Su operación es sencilla. Es más versátil que un reactor continuo.

El costo inicial es menor que el de los reactores continuos.
Desventajas.
El costo de operación es mayor que el de los reactores continuos. Requiere un ciclo
de operación complicado.

Se utilizan generalmente para líquidos.

Cuando se desea una velocidad de producción pequeña.
Para estudios cinéticos de laboratorio.
Cuando se desea obtener un producto muy puro.

El reactor intermitente es un sistema cerrado; por lo tanto la masa total del mismo
es constante.
El tiempo de residencia de todos los elementos del fluido es el mismo.
Este tipo de reactor opera en estado no estacionario; las concentraciones cambian
con el tiempo.
Sin embargo, y debido a la agitación puede asumirse que en un tiempo
determinado, dentro del reactor la composición y la temperatura son homogéneas.
 El reactor batch
dC A dC A dx
Balance de materia: −VR =V ( ) R=− = C A0
dt dt dt

Balance de energía (entalpías)

dT
(−ΔH R )VR − q =VCp
dt

Q0=0
CA0:Concentración inicial del reactivo limitante
T0: Temperatura a la entrada del elemento de volumen.
Q=0
CA: Concentración final del reactivo limitante
T: Temperatura a la salida del elemento de volumen.[ºC,ºK]
V: Volumen del elemento de reacción [L3]
CA: Concentración del RL en el interior del elemento[M/L3]
q: Calor transferido entre el elemento de volumen y los alrededores [Energía/tiempo]
R: Ecuación cinética, rapidez de reacción = f(T, C)
T: Tiempo de reacción
 Cálculo del reactor
batch
ECUACIONES DE DISEÑO
tiempo operación CONVERSION A V=CTE
PRODUCCIÓN:
t1 Carga de reactivos C A0 − C A
VC A0 x x=
t2 precalentamiento P= C A0
tCICLO
t3=tR Reacción química Sea la reacción:
A =èPRODUCTOS
t4 Enfriamiento
dC dx
t5 Descarga de productos R = − A = C A0 = kC An = kC A0 (1− x)
dt dt
t6 Limpieza
E
tCICLO Tiempo del ciclo (Σti)
− Ecuación de Arrhenius
k = Ae RT

dC A dx BALANCE GENERAL DE MATERIA

R=− = C A0 PARA EL REACTOR POR LOTES
dt dt
dT BALANCE GENERAL DE ENERGÍA
(−ΔH R )VR − q =VCp
dt PARA EL REACTOR POR LOTES
 Cálculo del reactor
batch Isotérmico
ECUACIONES DE DISEÑO
CONVERSION A V=CTE
PRODUCCIÓN: Condiciones T= cte
VC A0 x C A0 − C A k=cte.
P= x=
C A0 BALANCE GENERAL DE ENERGÍA
tCICLO PARA EL REACTOR POR LOTES
Sea la reacción:
A =èPRODUCTOS ISOTERMICO
dC dx (−ΔH R )VR − q = 0
R = − A = C A0 = kC An = kC A0 (1− x)
dt dt
BALANCE GENERAL DE MATERIA
PARA EL REACTOR POR LOTES
ISOTERMICO
dC A dx
R=− = C A0
dt dt
Ecuación de Arrhenius E
−
k = Ae RT
 Ejercicio
LA REACCIÓN DE SAPONIFICACIÓN
NaOH (ac) + C2 H 5 (CH 3CO) → Na(CH 3CO)(ac) + C2 H 5OH (ac)
ES DE SEGUNDO ORDEN E IRREVERSIBLE A
CONVERSIONES BAJAS. UN REACTOR DE LABORATORIO
CON BUEN MEZCLADO SE CARGA CON UN A SOLUCIÓN
ACUOSA QUE CONTIENE NaOH Y ACETATO DE ETILO,
AMBOS A CONCENTRACIONES INICIALES IGUALES 0.1
NORMALES. DESPUÉS DE 15 MIN, LA CONVERSIÓN DE
ACETATO DE ETILO 18 %.
PARA UNA CARGA INICIAL QUE CONTIENE NaOH Y
ACETATO DE ETILO CON HIDROXIDO DE SODIO 0.2
NORMALES, CALCULE EL TIEMPO REQUERIDO PARA
OBTENER UINA CONVERSIÓN DEL 30 % EN UN REACTOR
INTERMITENTE DE TANQUE CON AGITACIÓN A ESCALA
COMERCIAL.¿QUE VOLUMEN ES NECESARIO PARA
PRODUCIR 50 Kg DE ACETATO DE SODIO?
 Cálculo del reactor
batch Adiabático
ECUACIONES DE DISEÑO
CONVERSION A V=CTE Condiciones q= 0
PRODUCCIÓN: k=Ae-E/RT
VC A0 x C A0 − C A
x = BALANCE GENERAL DE ENERGÍA
P= C A0
tCICLO PARA EL REACTOR POR LOTES
Sea la reacción: ADIABÁTICO dT
A =èPRODUCTOS (−ΔH R )VR =V ρCp
dC dt
dx
R = − A = C A0 = kC An = kC A0 (1− x)
dt dt Temperatura de (−ΔH r )C A0 x
T = T0 +
reacción adiabática ρCp

BALANCE GENERAL DE MATERIA

PARA EL REACTOR POR LOTES dC A dx
ADIABÁTICO R=− = C A0
dt dt
Ecuación de Arrhenius −
E
k = Ae RT
 3.0. CÁLCULOS
NO HAY REACCIÓN DURANTE EL
PRECALENTAMIENTO

S3.1. SECUENCIA 1
S3.1.1. PRECALENTAMIENTO A 55 ºC Ley de Fourier
S3.1.1.1. BALANCE DE MATERIA
dC A dx
R=− = C A0 = kC An = 0
dt dt
S3.1.1.2. BALANCE DE ENERGÍA

dT R=0 q = UAΔT = UA(Ts − T )

(−ΔH R )VR − q =VCp
dt dT
−q = ρVCp U: Coef. De transferencia de calor
dt A: Área de transferencia

3.1.1.3. TIEMPO DE PRECALENTAMIENTO

DEL BALANCE DE ENERGIA CON R=0 Y LA EC. DE TRANSFERENCIA DE
CALOR:
dT ρVCp T dT ρVCp (Ts − T0 )
−UA(Ts − Ts ) = ρVCp t =− ∫ t= ln
dt UA T0
(Ts − T ) UA (Ts − T )
 3.0. CÁLCULOS
S3.1. SECUENCIA 1
S3.1.2.REACCION ADIABATICA
S3.1.2.1. BALANCE DE MATERIA
dC A dx
R=− = C A0 = kC An −
E
dt dt k = Ae RT
S3.1.2.2. BALANCE DE ENERGÍA

dT (−ΔH r )C A0 x
(−ΔH R )VR =VCp T = T0 +
dt ρCp

3.1.2.3. TIEMPO DE REACCIÓN ADIABATICA

SE REQUIERE SOLUCIÓN NUMÉRICA SIMULTANEA ENTRE LOS BALANCES
DE MATERIA Y ENERGÍA
 a) Análisis cualitativo: Analice ambas políticas en un X vs T y prediga en que caso se obtendrá una mayor producción.
Indique que política insumirá un mayor gasto desde el punto de vista del intercambio calórico.
b) Halle el tiempo de reacción para cada una de las políticas propuestas.
Datos:
conversión final: 90%

3.0. CÁLCULOS
6
r = k.C , con k=4.10 .exp(-7900/T)(l/seg)
A
ΔH=-40.000cal/mol
δ=1Kg/l
cp=1cal/(g.°C)
CAo = 4 mol/lt
3

V=5m
S3.1. SECUENCIA 1
BALANCE DE ENERGIA PÁRA UNA UNIDAD DE VOLUMEN REACCIONANTE
S3.1.2.REACCION ADIABATICA
S3.1.2.4. RUNGE- KUTTA
QCp(T0DE
– T) +4º
V ∑ORDEN
(-ΔHi) Ri = VCpdT/dt – q

RUNGE – KUTTA
SEA LA ECUACION dx/dz = f(x)
Ecuaciones de trabajo
K0=ΔZ f(xn) K3= ΔZ f(xn + K2)
K1=ΔZ f(xn + K0/2)
K2= ΔZ f(xn + K1/2) (x)n+1= xn + (1/6) (K0 + 2K1 + 2K2 + K3 )

SDel balance de materia y la ecuación de Arrhenius :

E
dx n
−
C A0 = kC A0 (1− X ) n = Ae RT C A0
n
(1− X ) n
dt

dx (−ΔH r )C A0 x
dt = C A0 E T = T0 +
Ae
−
RT n
C A0 (1− X ) n ρCp
 3.0. CÁLCULOS
Δx = Incremento de la conversión
S3.1.2.4. RUNGE- KUTTA DE 4º ORDEN
Sea: Establecer un incremento
1 Algoritmo:
f (x) = C A0
−
E
n
1. Con el valor x0 = 0 calcular T0 usando el
Ae RT
C A0 (1− X ) n balance de energía (−ΔH r )C A0 x
T = T0 +
k0 = Δx ∗ f (xn ,Txn ) ρCp
2. con x0 y T0 calcular f(x0)
k1 = Δx ∗ f (x Δx
,Tx ) 3. Calcular k0
n+ f (x ,Tx )
2 n+
Δx
2 4. Calcular x Δx = x Δx Tx =T
0+
Δx n+
Δx
Δx
n+ 0+ n+
Δx
2 2 n+
2 2 2 2

k2 = Δx ∗ f (x Δx
,Tx ) 5. Calcular k1
n+ n+
Δx
2 2 6. Calcular k2
k3 = Δx ∗ f (xn+Δx ,Tn+Δx ) 7. Calcular xn+Δx = x0+Δx Tn+Δx = T0+Δx f (xn+Δx ,Tn+Δx )
8. Calcular k3
k0 + 2k1 + 2k2 + k3 9. Calcular xn+1=x1
xn+1 = x0 + Fin del primer incremento, retornar a 1
6
 3.0. CÁLCULOS
S3.1.2.5. INTEGRACIÓN NUMERICA
SDel balance de materia y energía ( Temperatura de reacción adiabática):
Regla del trapecio
C A0 x dx
t= n ∫ E
AC A0 0 −
(−ΔH R )C A0 x
R(T0 +
e ρCp
(1− x) n Regla del trapecio compuesta

Regla de Simpson
 3.0. CÁLCULOS

S3.1.3. TIEMPO DE ENFRIAMIENTO

SDel balance de energía con R=0

ρVCp (T0 − Te )
t= ln
UA (T − Te )

3.1.4. CICLO DE OPERACIÓN

tc = ∑ ti = tllenado + tvaciado + t prec + treacción + tenfriamiento

i

3.1.5. VAPOR CONSUMIDO

Gλ = ρVCp(ΔT )
 4.0. OPERACIONES

S4.1.. TIEMPO DE PRECALENTAMIENTO

ρVCp (Ts − T0 )
t= ln
UA (Ts − T )

5m3 * 4.2 *106 J / m3 ºC 120 − 20

t= Ln = 2704.873s
1360(W / m2 ºC) *3.5m2 120 − 55

4.2. TIEMPO DE REACCIÓN ADIABATICA

SBalance de energía:

1.67 *106 (J / Kg) *1(Kgmol / m3 ) *100Kg / kgmol

T = 20ºC + 6 3
*x
4.2 *10 (J / m ºC)
 4.0. OPERACIONES

4.2. TIEMPO DE REACCIÓN ADIABATICA

SRUNGE-KUTTA
TIEMPO X TºK
10079.5072 0.1 293.15
7898.18691 0.2 297.1262
6335.70692 0.3 301.1024
5222.31294 0.4 305.0786
4449.08762 0.5 309.0548
3955.82321 0.6 313.031
3736.16096 0.7 317.0072
3890.53277 0.8 320.9834
4962.84058 0.9 324.9596

https://youtu.be/qgcs2X35Thc

SAPONIFICACIÓN
 Ejercicio
Ensayos de laboratorio sobre la hidrólisis de anhídrido
acético con una concentración inicial CA0= 0.215 mol/l,
dieron como resultado que la velocidad específica de la
reacción es k = 5.882 x 107 e-5888/T min-1.
Un reactor discontinuo se carga con 250 Kg de la
disolución de anhídrido acético en agua con la misma
concentración y una densidad de 1050 kg/l a 15 ºC. La
capacidad calorífica puede considerarse constante (0.9
kcal/kg.ºC.) y el calor de reacción es de -50 kcal/mol.
Calcule el tiempo necesario para alcanzar una
conversión del 80 % en las siguientes condiciones:
1. Operación isotérmica a 15 ºC.
2. Operación adiabática con una temperatura inicial de
15 ºC.
REACTOR IDEAL
DE FLUJO PISTON
 REACTOR IDEAL
DE FLUJO PISTON
BALANCE BALANCE
DE DE
MATERIA ENERGÍA

C A0 − C A
x=
C A0
d(VC A ) dC A dV
Q(C A0 − C A ) −VR = =V + CA
dt dt dt
d(VC A )
regimen − permanente =
dt
=0 QρCp(T0 − T ) + (−ΔH r )VR − q = 0
Q(C A0 − C A ) −VR = 0
−QρCpdT + (−ΔH r )RdV = dq
Q(C A0 − C A ) =VR
−QdC A = RdV
−QρCpdT + (−ΔH r )RdV = dq
dV dC
=− A
Q R
V C A dC
(−ΔH r )C A0 x
Θ= =−∫ A T = T0 +
Q C A0
R ρCp
 REACTOR IDEAL
DE TANQUE CONTINUO
AGITADO(CSTR)
BALANCE BALANCE
DE DE
MATERIA ENERGÍA

C A0 − C A
x=
C A0
d(VC A ) dC A dV
Q(C A0 − C A ) −VR = =V + CA
dt dt dt
regimen − permanente =
d(VC A )
=0 QρCp(T0 − T ) + (−ΔH r )VR = 0
dt
Q(C A0 − C A ) −VR = 0 (−ΔH r )VR = q
Q(C A0 − C A ) =VR
(−ΔH r )C A0 x
T = T0 +
Θ=
V (C A0 − C A )
= ρCp
Q R
 TAREA

REACTOR ISOTERMICO ADIABATICO NO ADIABATICO/NO ISOTERMO

BALANCE DE BALANCE BALANCE DE BALANCE DE ENERGÍA BALANCE DE BALANCE DE ENERGÍA
MATERIA DE MATERIA MATERIA
ENERGÍA
BATCH R=−
dC A
= C A0
dx (−ΔH r )VR =Rq= − dC A = C dx (−ΔH r )C A0 x dC A dx dT
dt dt A0 T = T0 + R=− = C A0 (−ΔH R )VR − q = VCp
dt dt dt dt dt
ρCp
FLUJO V C A dC (−ΔH r )VR = q V C A dC (−ΔH r )C A0 x V C A dC −Q ρCpdT + (−ΔH )RdV
PISTON Θ= =−∫ A
Θ= =−∫ A
T = T0 + Θ= =−∫ A r
Q C A0
R Q C A0
R ρCp Q C A0
R

CSTR V (C A0 − C A ) (−ΔH r )VR = q V (C A0 − C A )QρCp(T − T ) + (−ΔH )VR = 0 V (C A0 − C A )QρCp(T0 − T ) + (−ΔH r )VR

Θ= = Θ= = 0 r Θ= =
Q R Q R (−ΔH r )C A0 x Q R
T = T0 +
ρCp
CONVERSI C A0 − C A (−ΔH r )VR = q C A0 − C A C A0 − C A
ON x= x= x=
C A0 C A0 C A0
 CSTR - ADIABATICO
1. Los gases de escape de los automóviles contiene algo de combustible y
monóxido de carbono sin quemar, sin embargo la combustión puede
continuar si se le instala un convertidor catalítico en la línea de escape,
Suponga que el convertidor opera como un reactor de mezcla completa
adiabático con un tiempo de residencia de 9 segundos y que solo se oxida
el CO a CO2 , además, que con un exceso de aire la reacción es de primer
orden irreversible con los siguientes datos:
K(s-1) = 1.5 x 1010 EXP(- 65000/RT) en donde E (cal/ gmol)
3
Cp= 6.95 cal/gmol °K, ρ = 1.20 Kg/m ΔH = - 5500 Cal/gmol
Si los gases de combustión entran al convertidor a 1073 ° K, cuales serian
los valores posibles de las conversiones y temperaturas de salida de los
gases?

q
Q, Ca0,
T0 Ca, T
 DATOS BM Y BE
X T BM T BE
0 1073 #¡DIV/0!
0.005 1076.29736 1051.03688
0.01 1079.59472 1075.31647
0.02 1086.18945 1100.93505
0.03 1092.78417 1116.66046
0.04 1099.3789 1128.21389
0.05 1105.97362 1137.43825
0.06 1112.56835 1145.16959
0.07 1119.16307 1151.86018
0.08 1125.75779 1157.78312
0.09 1132.35252 1163.11631
0.1 1138.94724 1167.98226
0.11 1145.54197 1172.46905
0.12 1152.13669 1176.64221
0.13 1158.73141 1180.5518
0.14 1165.32614 1184.23701
0.15 1171.92086 1187.72911
0.16 1178.51559 1191.05345
0.17 1185.11031 1194.23089
0.18 1191.70504 1197.27883
0.19 1198.29976 1200.21194
0.2 1204.89448 1203.0427
0.21 1211.48921 1205.78185
0.22 1218.08393 1208.43868
0.23 1224.67866 1211.02133
0.24 1231.27338 1213.53694
0.25 1237.86811 1215.99185
0.26 1244.46283 1218.39168
0.27 1251.05755 1220.74149
0.28 1257.65228 1223.04585
0.29 1264.247 1225.30886
 GRAFICA TvsX BM Y BE
REACTOR CSTR ADIABATICO

1250

1200
TEMP °K

1150

1100

1050

1000
0 0.05 0.1 0.15 0.2 0.25
CONVERSION

DOS SOLUCIONES (0.011, 1075.31 K) (0.19, 1200 K)

 PROBLEMA 3
La figura 2 muestra Ca0/R contra X, para una descomposición del reactivo A en fase líquida, no
isotérmica, no elemental de múltiples reacciones
Considere los sistemas que se muestran en la figura 3. En donde un CSTR y un RFP están
conectados en serie. La conversión intermedia es de 0.3 y la conversión final es de 0.45.
¿Cómo deben acomodarse los reactores para que el volumen total sea mínimo? Explique
Si el flujo volumétrico es de 50 l/min, ¿qué volumen tendrán los reactores?
Si utiliza dos CSTR en serie del mismo volumen del inciso anterior, que conversión alcanza?
Si utiliza dos RFP en serie del mismo volumen del inciso anterior, que conversión alcanza?
12

10 X=0.45

0.45
X=0.45
1/R

2 0.45

0
0 0.1 0.2 0.3 0.4 0.5 0.6
x(conversion)
 Si el flujo volumétrico es de 50 l/min, ¿qué volumen tendrán lo
Si utiliza dos CSTR en serie del mismo volumen del inciso ant
Si utiliza dos RFP en serie del mismo volumen del inciso anter
12

10

θCSTR = BxH = 0.3* 2.2 = 0.66

1/R

6 VCSTR
θCSTR =
Q
VCSTR AREA DEL
= θ *Q = 0,66 *50 = 33L TRAPECIO
4
(B + b) * H
B A=
2
2

H
0
0 0.1 0.2 0.3 0.4 0.5 0.6
x(conversion)