REACTOR ISOTÉRMICO IRREVERSIBLE

El reactor es un reactor mezcla completa agitado (CSTR, continued Stirred-tank Reactor). El

componente A reacciona irreversiblemente a una razón específica k para formar el
componente B

Asumase que la concentración de los componentes en el flujo de entrada es como se muestra,

y que la reacción es de primer orden

Asumimos que el tanque esta perfectamente mezclado y que la temperatura del sistema es
constante. Además asumamos que la velocidad de reacción es directamente proporcional a la
concentración instantánea dentro del tanque y que no hay gradientes de concentración.

Suposiciones
• el tanque esta perfectamente mezclado
• la temperatura del sistema es constante
• el volumen del liquido contenido es constante
• La velocidad de reaccion es proporcional a la concentracion instantanea
• no hay gradientes de concentracion.
 Flujo de A a la entrada: Fo⋅ CAo
flujo de A a la salida: F ⋅ CA

velocidad de formación de A: −V ⋅ k ⋅ CA

Moles de A en el reactor V ⋅ CA

Fo⋅ CAo − F ⋅ CA − V ⋅ k ⋅ CA =
d
dt
(V ⋅ CA )
1 ecuacion, 4 incognitas

Flujo de B a la entrada: Fo⋅ CBo

Flujo de B a la salida F ⋅ CB
Velocidad de formacion de B V ⋅ k ⋅ CA
moles de B en el reactor V ⋅ CB

Fo⋅ CBo − F ⋅ CB + V ⋅ k ⋅ CA =
d
dt
(
V ⋅ CB )
2 ecuaciones, 6 incognitas

tomemos como incognitas solo a Ca y a Cb, las otras son variables

libres (variables de diseño, variables de control o perturbaciones)
 d
Fo⋅ CAo − F ⋅ CA − V ⋅ k ⋅ CA = V CA
dt
linealizando los terminos

( ) (
Fo⋅ CAo = Fo ( 0) CAo ( 0) + Fo ( 0) CAo ( t ) − CAo ( 0) + CAo ( 0) Fo ( t ) − Fo ( 0) )
( )
F ⋅ CA = F ( 0) CA ( 0) + F ( 0) CA ( t ) − CA ( 0) + CA ( 0) ( F ( t ) − F ( 0) )

la ecuacion del estado estacionario

Fo ( 0) ⋅ CAo ( 0) − F ( 0) ⋅ CA ( 0) − V ⋅ k ⋅ CA ( 0) = 0
las variables de desviacion seran
Γo ( t) = Fo ( t) − Fo ( 0) DAo ( t ) = CAo ( t ) − CAo ( 0)
Γ ( t ) = F ( t) − F ( 0) DA ( t ) = CA ( t ) − CA ( 0)

sustituyendo en la ecuacion diferencial

d
Fo ( 0) ⋅ DAo ( t ) + CAo ( 0) ⋅ Γo ( t ) − F ( 0) ⋅ DA ( t ) − CA ( 0) ⋅ Γ ( t ) − V ⋅ k ⋅ DA ( t ) = V DA ( t )
dt
de igual manera, para la segunda ecuacion

( ) (
Fo⋅ CBo = Fo ( 0) CBo ( 0) + Fo ( 0) CBo ( t ) − CBo ( 0) + CBo ( 0) Fo ( t ) − Fo ( 0) )
( )
F ⋅ CB = F ( 0) CB ( 0) + F ( 0) CB ( t ) − CB ( 0) + CB ( 0) ( F ( t ) − F ( 0) )

la ecuacion del estado

estacionario
Fo ( 0) ⋅ CBo ( 0) − F ( 0) ⋅ CB ( 0) + V ⋅ k ⋅ CA ( 0) = 0
las variables de desviacion seran
DBo ( t ) = CBo ( t ) − CBo ( 0)
DB ( t ) = CB ( t ) − CB ( 0)

Fo⋅ CBo − F ⋅ CB + V ⋅ k ⋅ CA =
d
dt
(
V ⋅ CB )
al sustituir se tiene que
d
Fo ( 0) ⋅ DBo ( t ) + CBo ( 0) ⋅ Γo ( t ) − F ( 0) ⋅ DB ( t ) − CB ( 0) ⋅ Γ ( t ) + V ⋅ k ⋅ DA ( t ) = V DB ( t )
dt
de la primera ecuacion, despejando se tiene que

Fo ( 0) CAo ( 0) F ( 0) CA ( 0) d
⋅ DAo ( t ) + ⋅ Γo ( t ) − ⋅ DA ( t ) − ⋅ Γ ( t) − k ⋅ DA ( t) = DA ( t)
V V V V dt
al aplicar la transformada de laplace

Fo ( 0) CAo ( 0) F ( 0) CA ( 0)
⋅ DAo ( s) + ⋅ Γo ( s) − ⋅ DA ( s) − ⋅ Γ ( s) − k ⋅ DA ( s) = s⋅ DA ( s)
V V V V

al hacer lo mismo con la segunda ecuacion

Fo ( 0) CBo ( 0) F ( 0) CB ( 0)
⋅ DBo ( s) + ⋅ Γo ( s) − ⋅ DB ( s) − ⋅ Γ ( s) + k ⋅ DA ( s) = s⋅ DB ( s)
V V V V

F o ( 0) CAo ( 0) CA ( 0)
V V V
DA ( s) = ⋅ DAo ( s) + ⋅ Γo ( s) − ⋅ Γ ( s)
F ( 0) F ( 0) F ( 0)
s+k+ s+k+ s+k +
V V V

llamemos CA ( 0)
F o ( 0) CAo ( 0)
1
τ1 = V V K3 =
V
F ( 0) K1 = K2 =
k+ F ( 0) F ( 0) k+
F ( 0)
V k+ k+ V
V V
 K1 K2 K3
DA ( s) = ⋅ DAo ( s) + ⋅ Γo ( s) − ⋅ Γ ( s)
τ1⋅ s + 1 τ1⋅ s + 1 τ1⋅ s + 1

de igual manera, de la segunda ecuacion

Fo ( 0) CBo ( 0) F ( 0) CB ( 0)
⋅ DBo ( s) + ⋅ Γo ( s) − ⋅ DB ( s) − ⋅ Γ ( s) + k ⋅ DA ( s) = s⋅ DB ( s)
V V V V

Fo ( 0) CBo ( 0) CB ( 0)
V V V k
DB ( s) = ⋅ DBo ( s) + ⋅ Γo ( s) − ⋅ Γ ( s) + ⋅ D ( s)
F ( 0) F ( 0) F ( 0) F ( 0) A
s+ s+ s+ s+
V V V V
llamemos

CBo ( 0)
F o ( 0) CB ( 0) k
F ( 0) V K7 =
τ2 = V K5 = V F ( 0)
V K4 = F ( 0) K6 = V
F ( 0) F ( 0)
V
V V

K4 K5 K6 K7
DB ( s) = ⋅ DBo ( s) + ⋅ Γo ( s) − ⋅ Γ ( s) + ⋅ DA ( s)
τ2⋅ s + 1 τ2⋅ s + 1 τ2⋅ s + 1 τ2⋅ s + 1
el modelo es, entonces, en resumen

K1 K2 K3
DA ( s) = ⋅ DAo ( s) + ⋅ Γo ( s) − ⋅ Γ ( s)
τ1⋅ s + 1 τ1⋅ s + 1 τ1⋅ s + 1
K4 K5 K6 K7
DB ( s) = ⋅ DBo ( s) + ⋅ Γo ( s) − ⋅ Γ ( s) + ⋅ DA ( s)
τ2⋅ s + 1 τ2⋅ s + 1 τ2⋅ s + 1 τ2⋅ s + 1

Γo ( t) = Fo ( t) − Fo ( 0) Γ ( t ) = F ( t) − F ( 0)
DAo ( t ) = CAo ( t ) − CAo ( 0) DA ( t ) = CA ( t ) − CA ( 0)
DBo ( t ) = CBo ( t ) − CBo ( 0) DB ( t ) = CB ( t ) − CB ( 0)

F o ( 0) CAo ( 0) CA ( 0)
1 F ( 0)
τ1 = τ2 = V V V
F ( 0) V K1 = K2 = K3 =
k+ F ( 0) F ( 0) F ( 0)
V k+ k+ k+
V V V

F o ( 0) CBo ( 0) CB ( 0)
V V V k
K4 = K5 = K6 = K7 =
F ( 0) F ( 0) F ( 0) F ( 0)
V V V V
Parámetros:
Fo ( t ) := 1 F ( t ) := 1 V := 20 CAo ( t ) := 1 CBo ( t ) := 0.5

CAi := 1 CBi := 0.2 k := 0.1 tend := 40

Dado
Fo ( t ) F ( t) d
⋅ CAo ( t ) − ⋅ CA ( t ) − k ⋅ CA ( t ) = CA ( t ) CA ( 0) = CAi
V V dt
Fo ( t ) F ( t) d CB ( 0) = CBi
⋅ CBo ( t ) − ⋅ CB ( t ) + k ⋅ CA ( t ) = CB ( t )
V V dt

 CA   C   MODELO NO LINEAL
  := Odesolve A  , t , tend
 CB   CB  
    
1.5

CA ( t) 1

CB ( t)
0.5

0
0 10 20 30 40
t
 CAeq := CA ( tend ) CBeq := CB ( tend )

Γo ( t) := Fo ( t) − Fo ( 0) = 0 DAo ( t ) := CAo ( t ) − CAo ( 0) = 0

Γ ( t ) := F ( t) − F ( 0) = 0 DBo ( t ) := CBo ( t ) − CBo ( 0) = 0.0

Dado
Fo ( 0) CAo ( 0) CA ( 0)
⋅ DAo ( t ) + ⋅ Γo ( t ) −
F ( 0)
⋅ DA ( t ) − ⋅ Γ ( t) − k ⋅ DA ( t) =
d
DA ( t) DA ( 0) = CAi − CAeq
V V V V dt
Fo ( 0) CBo ( 0) F ( 0) CB ( 0) d
⋅ DBo ( t ) + ⋅ Γo ( t ) − ⋅ DB ( t ) − ⋅ Γ ( t) + k ⋅ DA ( t) = DB ( t) DB ( 0) = CBi − CBeq
V V V V dt
 DA   D  
  := Odesolve A  , t , tend MODELO LINEAL
 DB   DB  
    
 1

0.5
DA ( t)
0
DB ( t)

− 0.5

−1
0 10 20 30 40
t

Comparación modelo lineal y no lineal

1.2

1
CA ( t)
0.8
DA ( t) +CAeq

CB ( t) 0.6

DB ( t) +CBeq 0.4

0.2

0
0 10 20 30 40
t
REACTOR NO ISOTÉRMICO IRREVERSIBLE
Considérese el reactor con una chaqueta de enfriamiento que puede remover el calor de reacción
exotérmico λ (joule/mol de A consumido). Asumimos que λ es negativo para una reacción
exotérmica y positivo para una reacción endotermica. La velocidad de generación de calor por
unidad de tiempo viene dada por la velocidad del reactivo consumido multiplicado por λ.
El reactor es un reactor mezcla completa agitado (CSTR, continued Stirred-tank Reactor). El
componente A reacciona irreversiblemente a una razón específica k para formar el
componente B
Asumase que la concentración de los componentes en el flujo de entrada es como se muestra,
y que la reacción es de primer orden
Suposiciones
• el tanque esta perfectamente mezclado
• el volumen del liquido contenido es constante
• La velocidad de reaccion es proporcional a la concentracion
instantánea
• la entalpia es solo funcion de la temperatura
• Cp no cambia con la temperatura
• Cv y Cp son esencialmente iguales para liquidos
• La densidad no cambia con la temperatura, y la densidad de
todas las sustancias es esencialmente la misma

Balance de materia:

Fo⋅ CAo − F ⋅ CA − V ⋅ k ⋅ CA =
d
dt
(V ⋅ CA )
Energia Entrante: Fo⋅ ρo ⋅ ho

Energia en el flujo de salida F ⋅ρ⋅h

Calor extraido Q
Calor generado por la reaccion λ⋅ V ⋅ k ⋅ CA
energia interna del sistema ρ⋅V ⋅U

Balance de energia:

d
Fo⋅ ρo ⋅ ho − F ⋅ ρ ⋅ h − Q − λ⋅ V ⋅ k ⋅ CA = (ρ ⋅ V ⋅ U )
dt

( ) ( )
Fo⋅ ρo ⋅ Cp⋅ To − Tref − F ⋅ ρ ⋅ Cp⋅ T − Tref − Q − λ⋅ V ⋅ k ⋅ CA =
d
(
ρ ⋅ V ⋅ Cv⋅ T − Tref 
dt 
)
La constante de velocidad depende de la temperatura
−E
R ⋅ T ( t)
k ( t ) = ko⋅ e

Sustituyendo en el balance de materia, se tiene

−
 −E
Fo ( t ) F ( t) R ⋅ T ( t) d
⋅ CAo ( t ) − ⋅ CA ( t ) − ko⋅ e ⋅ CA ( t ) = CA ( t )
V V dt

derivando y despejando en el balance de energía se tiene

Fo ( t ) ⋅ ρo ⋅ Cp F ( t ) ⋅ ρ ⋅ Cp Q ( t) λ⋅ V ⋅ k ( t) ⋅ CA ( t) d
ρ ⋅ V ⋅ Cv
( )
⋅ To ( t ) − Tref −
ρ ⋅ V ⋅ Cv
(
⋅ T ( t ) − Tref − )
ρ ⋅ V ⋅ Cv
−
ρ ⋅ V ⋅ Cv
= T (t)
dt

al sustituir la expresión de la constante de velocidad:

 −E 

λ⋅ ko⋅ e
R ⋅ T ( t) 
Fo ( t )
⋅ ( To ( t ) − Tref ) −
F (t)
⋅ ( T ( t ) − Tref ) −
Q (t)
−   ⋅ CA ( t ) = d T ( t )
V V ρ⋅V ⋅C v ρ⋅C v dt

El balance de materia en estado transitorio es:

−
 −E
Fo ( t ) F ( t) R ⋅ T ( t) d
⋅ CAo ( t ) − ⋅ CA ( t ) − ko⋅ e ⋅ CA ( t ) = CA ( t )
V V dt
y el balance de materia en estado estacionario
−E
Fo ( 0) F ( 0) R ⋅ T ( 0)
⋅ CAo ( 0) − ⋅ CA ( 0) − ko⋅ e ⋅ CA ( 0) = 0
V V

La linealización de los términos

( ) (
Fo ( t ) ⋅ CAo ( t ) = Fo ( 0) ⋅ CAo ( 0) + CAo ( 0) Fo ( t ) − Fo ( 0) + Fo ( 0) ⋅ CAo ( t ) − CAo ( 0) )
(
F ( t ) ⋅ CA ( t ) = F ( 0) ⋅ CA ( 0) + CA ( 0) ( F ( t ) − F ( 0) ) + F ( 0) ⋅ CA ( t ) − CA ( 0) )
−E
−E −E −E
R ⋅ T ( 0)
R ⋅ T ( t) R ⋅ T ( 0) Ee R ⋅ T ( 0)
e ⋅ CA ( t ) = e ⋅ CA ( 0) +
R T ( 0) 2
⋅ CA ( 0) ( T ( t ) − T ( 0) ) + e (CA ( t ) − CA ( 0) )
Las variables de desviación
ϕo ( t) = Fo ( t) − Fo ( 0) ϕ ( t) = F ( t ) − F ( 0)
DAo ( t ) = CAo ( t ) − CAo ( 0) DA ( t ) = CA ( t ) − CA ( 0) Γ ( t ) = T ( t) − T ( 0)
al restar el balance estacionario al estado transitorio y sustituir las variables de desviación
−E
−E
CAo ( 0) Fo ( 0) CA ( 0) R⋅ T ( 0)
F ( 0) E e R ⋅ T ( 0) d
⋅ ϕo ( t ) + ⋅ DAo ( t ) − ⋅ ϕ (t) − ⋅ DA ( t ) − ko⋅ ⋅ ⋅ CA ( 0) ⋅ Γ ( t ) − ko e ⋅ DA ( t) = DA ( t)
V V V V R T ( 0) 2 dt

El balance de energía en estado transitorio

 − 
  −E 

Fo ( t )

λ⋅ ko⋅ e
R ⋅ T ( t) 
⋅ CA ( t ) d
⋅ ( To ( t ) − Tref ) −
F (t)
⋅ ( T ( t ) − Tref ) −
Q (t)
−
  = T ( t)
V V ρ⋅V ⋅C v ρ⋅C v dt

El balance en estado estacionario

 −E 

λ⋅ ko⋅ e
R ⋅ T ( 0) 
Fo ( 0)
⋅ ( To ( 0) − Tref ) −
F ( 0)
⋅ ( T ( 0t ) − Tref ) −
Q ( 0)
−   ⋅ CA ( 0) = 0
V V ρ⋅V ⋅C ρ⋅C
v v
La linealización de términos

(
Fo ( t ) ⋅ To ( t ) − Tref ) = Fo (0)⋅ ( To (0) − Tref ) + ( To (0) − Tref ) ⋅ ϕo (t) + Fo (0)⋅ Γo (t)
F ( t ) ⋅ ( T ( t ) − Tref ) = F ( 0) ⋅ ( T ( 0) − Tref ) + ( T ( 0) − Tref ) ⋅ ϕ ( t ) + F ( 0) ⋅ Γ ( t )

sustituyendo se tiene que

 To ( 0) − Tref Fo ( 0) T ( 0) − Tref F ( 0)  d
 ⋅ ϕo ( t ) + ⋅ Γo ( t ) − ⋅ ϕ ( t) − ⋅ Γ ( t ) ...  = Γ ( t)
 V V V V  dt
−E
 −E 
⋅ R ⋅ T ( 0) ⋅
 −1 λ ko E e λ ko R ⋅ T ( 0) 
 + ρ ⋅ V ⋅ C ⋅ q ( t) − ρ ⋅ C ⋅ R ⋅ 2
⋅ CA ( 0) ⋅ Γ ( t ) −
ρ ⋅ C
⋅e ⋅ DA ( t) 
 v v T ( 0) v 
La suposición de que el volumen del tanque es constante y que la densidad es la misma, implica que
el flujo de entrada y el flujo de salida son iguales, esto implica lo siguiente
ϕo ( t) = ϕ( t)
−E
−E
CAo ( 0) − CA ( 0) Fo ( 0) Fo ( 0) R ⋅ T ( 0)
E e R ⋅ T ( 0) d
⋅ ϕo ( t ) + ⋅ DAo ( t ) − ⋅ DA ( t ) − ko⋅ ⋅ ⋅ CA ( 0) ⋅ Γ ( t ) − ko e ⋅ DA ( t ) = DA ( t )
V V V R T ( 0) 2 dt

 To ( 0) − T ( 0) Fo ( 0) Fo ( 0)  d
 ⋅ ϕo ( t ) + ⋅ Γo ( t ) − ⋅ Γ ( t ) ...  = Γ ( t)
 V V V  dt
−E
 −E 
 −1 λ⋅ ko E e R⋅ T ( 0) λ⋅ ko R ⋅ T ( 0) 
 + ρ ⋅ V ⋅ C ⋅ q ( t) − ρ ⋅ C ⋅ R ⋅ 2
⋅ CA ( 0) ⋅ Γ ( t ) −
ρ ⋅ C
⋅e ⋅ DA ( t) 
 v v T ( 0) v 