Reactores Qumicos II.

Reactores Homogneos

II. Reactores Homogneos

Tema 2. Reactor discontinuo de tanque agitado Tema 3. Reactor continuo de tanque agitado Tema 4. Reactor tubular continuo Tema 5. Comparacin de reactores

Reactor Continuo de tanque agitado

Caractersticas Diseo del reactor: operacin isotrmica y no isotrmica Secuencia de reactores continuos de tanque agitado Estabilidad del reactor en estado estacionario. Anlisis general de estabilidad. Planos de fases: periodo de puesta en marcha

Reactor Continuo de tanque agitado

APLICACIONES:
PRODUCCIONES GRANDES REACCIONES MUY SENSIBLES A LA TEMPERATURA
(PELIGRO DE EXPLOSIN)

TIEMPOS DE RESIDENCIA ELEVADOS PREFERENTEMENTE REACCIONES EN FASE LQUIDA

FLUJO DE MEZCLA PERFECTA

HOMOGENEIDAD EN TODO EL REACTOR COMPOSICIN, PROPIEDADES FSICAS, TEMPERATURA... IGUALES EN EL INTERIOR DEL REACTOR Y EN LA CORRIENTE DE SALIDA

Tema 2. Reactor Continuo de tanque agitado

Diseo del reactor:

Tiempo de residencia Intercambio de calor Balance de Materia Balance de Energa

Sistema de agitacin

Tipo, Potencia

Volumen del Reactor

Tiempo de residencia

Reactor Continuo de Tanque Agitado

Balance de materia (componente j)
FA0 cA0 FA cA FA cA
velocidad de flujo de entrada flujo neto de = - salida de j + [conversin de j / tiempo] acumulacin de j de j

dN j dt

= FjE FjS + r j V

Balance de energa calorfica

velocidad de flujo de entrada flujo de flujo de calor acumulacin de calor = de calor - salida de calor + [calor de reaccin/tiempo ] + intercambiado

MTcp

dT = dt

F (H H ) + V (H ) r + q
jE jE jS ri i

Reactor Continuo de Tanque Agitado

balance de materia:
flujo de flujo de - = FjE - FjS entrada de j salida de j
- generacin = rj V + desaparici nREACCIN

[Acumulacin] =

dN j dt

balance de energa:
cantidad de calor aportada o retirada por intercambio de calor: calor generado en la reaccin:

[calor intercambi ado] = US(Tm - T ) = q

q = 0 (adiabtic o)

[calor

reaccin] = (- Hr ) r V = (- Hr ) rj j V

calor acumulado en el reactor:

calor acumulado dT = MT Cp dt en el reactor

Flujo neto convectivo de calor:

flujo de entrada flujo de - = de calor salida de calor

F
j

jE

[H jE

- H jS = V cp (TE - T )

RCTA - ISOTRMICO
Diseo del Reactor:
B.M. para la especie j =A (reactivo limitante):
FAo - FA = ( rA ) V F - FA V = Ao ( rA )
Fase lquida (V = cte) FA = c A V FAo = c Ao V
. .

F X V = Ao A ( rA )

tM =

V C A0 - C A = & ( rA ) V

RCTA - ISOTRMICO
resolucin de las ecuaciones de diseo mediante:
mtodos grficos
representando 1/(-rA) vs XA: mtodo de las reas representando (-rA) vs XA: mtodo de las pendientes

RCTA - ISOTRMICO
resolucin de las ecuaciones de diseo mediante:
mtodos grficos mtodos analticos: para una reaccin A B ( r = k cAn)

si n = 0:

r =k

tM =

c A0 - c A k
c A0 - c A X A (1 + A ) = k cA k (1 - X A ) XA k (1 - X A )

si n = 1:

C ( 1- XA ) r = k A0 (1 + A X A )

tM =

Si A = 0 tM =
tM =

si n = 2:

c2 ( 1 - X A ) r = k A0 (1 + A X A )2
2

c A0 - c A X A (1 + A ) (1 + A X A ) = 2 k c2 k c A0 (1 - X A ) A k c A0 (1 - X A ) XA
2

Si A = 0 tM =

RCTA - ISOTRMICO
balance de energa:
simplificaciones:
reactor isotrmico (T = cte)

& V Cp (TE - T ) + (- Hr ) r V +
Obtencin de:

q
U S (T )

=0

cantidad de calor que se debe intercambiar (q) superficie de intercambio (S) temperatura del fluido de intercambio (Tm)

Reactor continuo de tanque agitado

Ejemplo: En
un reactor continuo de tanque agitado se produce la

reaccin en fase lquida A + B C, para la cual k = 1011exp (7000/T) L/molh. Si el reactor se alimenta con 100 L/h de una mezcla de ambos reactivos (CA0 = 2,5 mol/L; CB0 = 3,5 mol/L), opera en condiciones isotrmicas a 60 C y se desea alcanzar una conversin del 95%, determine el volumen necesario. las concentraciones de reactivos y productos en el reactor y en la corriente de salida son: C = C - (1 - X ) = 2,5 (1 - 0,95 ) = 0,125 mol A/L
A A0 A

CB = CB0 - C A0 X A = 3,5 - 2,5 0,95 = 1,125 mol B/L CC = C A0 X A = 2,5 0,95 = 2,375 mol C/L

r viene dada por: r = k C A CB = 10 11 exp y, el volumen: t = V = C A0 - C A = & (- rA ) V

M

(2,5 - 0,125 ) mol

10,44 mol L h L

mol mol 7000 L mol 0,125 1,125 = 10,44 L L L h 273 + 60 mol h

= 0,227 h

L & V = tM V = 0,227 h 100 23 L h

Secuencia de reactores continuos de tanque agitado (Isotrmico)

con un solo RCTA la conversin a la salida es baja si se trabaja con una batera de reactores en serie, es posible incrementar la XA para el mismo V del sistema el diseo se realiza a partir de los balances de materia:
CA0 CA1 CA2 CAn-2 CAn-1 CAn

n-1

C A1 = C A0 - (- rA1 ) t M1 C A2 = C A1 - (- rA2 ) tM2 para un tanque i : C An = C An-1 - (- rAn ) tMn

C Ai = C Ai-1 - (- rAi ) tMi

Secuencia de reactores continuos de tanque agitado (Isotrmico)

CA1 CAn-1

CA0

CA2

CAn-2

CAn

n-1

puesto que: es decir,

C Ai = C A0 (1 X Ai ) X Ai - X Ai-1 t Mi = C AO ( rAi ) C Ai-1 = C A0 (1 X Ai1 )

X An - X An-1 = Vn (- rAn ) FA0

para un tanque cualquiera:

X An = X A0 +

1 FA0

V (- r
n n=1

An

el volumen total del sistema: VT =

V
n =1

Secuencia de reactores continuos de tanque agitado (Isotrmico)

la resolucin puede realizarse analtica o grficamente:
resolucin analtica:
las ecuaciones anteriores son fcilmente resolubles cuando se lleva a cabo una reaccin de primer orden:
A P N =1 N=2 N=3
M M

r = kC A C A0 = C A1 + k 1C A1tM1 = C A1 (1 + k1t M1 )
-1 C A1 = C A2 + k 2 C A2 tM2 = C A2 (1 + k 2 tM2 ) = C A0 (1 + k1t M1 ) -1 1 C A2 = C A3 + k 3 C A3 tM3 = C A3 (1 + k 3 t M3 ) = C A0 (1 + k 1t M1 ) (1 + k 2 t M2 )

C A1 = C A0 (1 + k1t M1 )

-1

-1 -1 C A2 = C A0 (1 + k 1t M1 ) (1 + k 2 t M2 ) -1 -1 -1 C A3 = C A0 (1 + k1t M1 ) (1 + k 2 t M2 ) (1 + k 3 t M3 )

N=N

C An-1 = C An + k n C An tMn = C An (1 + k n t n ) = C A0 (1 + k 1t M1 ) (1 + k 2 t M2 ) (1 + k 3 t M3 ) K (1 + k n t Mn-1 )

-1 1 1

es decir,

C An = C A0

(1 + k t
i i =1

Mi

Secuencia de reactores continuos de tanque agitado (Isotrmico)

la resolucin puede realizarse analtica o grficamente:
resolucin analtica:
si todos los tanques operan a la misma temperatura y presentan el mismo tamao: k 1t M1 = k 2 t M2 = k 3 t M3 = KK = k n t Mn = k t M

C An =

(1 + kt M )n

C A0

tM =

1 C A0 n - 1 k C An

el lmite cuando N , proporciona una ecuacin anloga a la de un reactor tubular en el que tiene lugar una reaccin de primer orden:

tP =

1 C 1 C A0 - 1 = ln A0 n lm k C An k CA n

Secuencia de reactores continuos de tanque agitado (Isotrmico)

la resolucin puede realizarse analtica o grficamente:
resolucin analtica: resolucin grfica:
mtodo de las pendientes:

C An-1 = C An + (- rAn ) t Mn

(- rAn ) = -

1 1 C An + C An-1 t Mn t Mn

Secuencia de reactores continuos de tanque agitado (Isotrmico)

la resolucin puede realizarse analtica o grficamente:
resolucin analtica: resolucin grfica:
mtodo de las pendientes: mtodo de las reas: de forma anloga, representando 1/(-rA) vs XA

Secuencia de reactores continuos de tanque agitado (Isotrmico)

la resolucin puede realizarse analtica o grficamente:
resolucin analtica: resolucin grfica:
mtodo de las pendientes mtodo de las reas mtodo nmrico de Eldridge y Piret

C An =C An-1 - t M (- rAn )

Secuencia de reactores continuos de tanque agitado (Isotrmico)

Ejemplo:
Un reactivo A se descompone en fase acuosa en un sistema de

tres reactores asociados en serie. El primer tanque, de 60 L, resulta ser de doble volumen que el segundo y seis veces mayor que el tercero. Cuando se alimentan al sistema 1000 L/h de una disolucin acuosa de A con una concentracin inicial CA0 = 2 mol/L, la conversin de A resulta ser del 50, 66.5 y 71.27 % respectivamente. Determine la ecuacin cintica que describe la reaccin de descomposicin de A. la concentracin de A a la salida de cada reactor:
C A1 = C A0 (1 - X A1 ) = 2 (1 - 0,50 ) = 1 mol A/L C A2 = C A0 (1 - X A2 ) = 2 (1 - 0,665 ) = 0,67 mol A/L

C A3 = C A0 (1 - X A3 ) = 2 (1 - 0,712 ) = 0,576 mol A/L

a partir de la ecuacin de diseo, se determina (-rA) para cada reactor:

10

Secuencia de reactores continuos de tanque agitado (Isotrmico)

a partir de la ecuacin de diseo, se determina (-rA) para cada reactor:
t M1 = c -c V1 c A0 - c A1 2-1 mol = (- rA )1 = A0 A1 = = 16,66 & V1 (- rA )1 60 1000 L h V1 & V
1

t M2 =

V2 c A1 - c A2 c -c 1 - 0,67 mol = (- rA )2 = A1 A2 = = 11,00 & V2 (- rA )2 30 1000 L h V2 & V

2

t M3 =

V3 c A2 - c A3 c -c mol 0,67 - 0,576 = (- rA )3 = A2 A3 = = 9,40 & V3 (- rA )3 L h 10 1000 V3 & V3

si n = 1,

(- rA )1 (- rA )2 (- rA )3
C A1 = C A2 = C A3

= k = 16,44 h-1

Secuencia de reactores continuos de tanque agitado (Isotrmico)

optimizacin de reactores en serie:
determinar el volumen del sistema (XA,mx) si XA es conocida, minimizar el volumen del sistema

optimizacin de dos reactores conectados en serie:

CA0 CA1 CA2

B.M.a los dos reactores:

t M1 =

slo puede actuarse sobre t y CA1 t M2 CA1 depende de la cintica

C A0 - C A1 (- rA )1 t M = t M1 + t M2 C A1 - C A2 = (- rA )2

siendo t M1 t M2

11

Secuencia de reactores continuos de tanque agitado (Isotrmico)

Ejemplo: la reaccin irreversible A P se lleva a cabo en dos reactores conectados en serie de distinto volumen y que operan a la misma temperatura. La concentracin inicial de reactivo A es de 2 mol/L y se desea obtener una conversin final del 80 %. Determine la disposicin ms adecuada para los reactores si n = 1, n = 2 y n = 0,5. Datos: k = 0,8 h-1; FA0 = 100 mol/h
el flujo volumtrico de entrada: la concentracin de salida: n=
100 mol/h & F V = A0 = = 50 L/h C A0 2 mol/L

C A2 = C A0 (1 - X A ) = 2 (1 - 0,8) = 0,4 mol/L

(C A0 + C A1 )2
4 C3 A1

1 = C A1 = 0,774 mol/L C A2

t M1 = t M2

C A0 - C A1 2 - 0,774 L = = 4,6 h V1 = 50 4,6 h = 229 L 1/2 h k 1 C1/2 0,8 (0,774 ) A1 C - C A2 0,774 - 0,4 L = A1 1/2 = = 0,74 h V2 = 50 0,74 h = 37 L 1/2 h k 2 C A2 0,8 (0,4 )

V1 > V2

Secuencia de reactores continuos de tanque agitado (Isotrmico)

Ejemplo (cont.)
n=1
C A1 = C A0 C A2 = 2 0,4 = 0,89 mol/L
t M1 = t M2 C A0 - C A1 2 - 0,89 L = = 1,55 h V1 = 50 1,55 h = 77,5 L h k 1 C A1 0,8 (0,89 ) C A1 - C A2 0,89 - 0,4 L = = = 1,50 h V2 = 50 1,50 h = 77 L h k 2 C1/2 0,8 (0,4 ) A2 V1 V2

n=2
2 C A0 - C A1 1 = 2 y C A1 = 1,04 mol/L C3 C A2 A1
t M1 = t M2 C A0 - C A1 L 2 - 1,04 = = 3,42 h V1 = 50 3,42 h = 170 L 2 h k 1 C1/2 0,8 (1,04 ) A1 C - C A2 1,04 - 0,4 L = A1 1/2 = = 5 h V2 = 50 5 h = 250 L 2 h k 2 C A2 0,8 (0,4 ) V1 < V2

12

Reactor Continuo de Tanque Agitado

diseo de un RCTA no isotrmico:
las variaciones de temperatura se deben a la entalpa de reaccin y van a afectar a la cintica y/o el equilibrio de la reaccin

La forma de operar puede ser:

adiabtica: sin intercambio de energa con los alrededores (q = 0) no isotrmica no adiabtica: intercambio con un fluido (q)

Clculo: resolucin simultnea de los balances de materia y energa

(estado estacionario)
FAo - FA = ( rA ) V c Ao - c A = ( rA ) t M

& V Cp (TE - T ) + (- Hr ) r V + U S (Tm - T ) = 0

[ C.I.:

t=0

To , cAo ]

RCTA: anlisis de estabilidad

operacin autotrmica del reactor es aquella en la que la regin estable se alcanza sin entrada de energa al sistema: la operacin es adiabtica el calor necesario procede del sistema puede haber ms de un conjunto de valores de las variables de operacin que conduzcan al estado estacionario en ocasiones puede ser necesario el encendido del reactor o el sistema puede presentar inercia al apagado del mismo, producindose un ciclo de histresis en el encendidoapagado del reactor y existiendo valores crticos para el encendido y apagado del sistema

13

RCTA: anlisis de estabilidad

el balance de energa al sistema est formado por dos trminos:
UG =

(- H r ) r V
& V Cp

= tM J r

calor generado en la reaccin calor eliminado/aportado en la reaccin

U E = T (1 +m )- (TE + m Tm ) Funcin de Generacin de Calor UG

Funcin de Eliminacin/Aporte de Calor UE

En estado estacionario debe cumplirse que UG = UE

para una reaccin de primer orden:

UG = J C A0

k tM 1 + k tM

UE = T (1 + m) - (TE + m Tm )

Se representan las funciones UG / UE vs. T (K)

RCTA NO ISOTRMICO
Ejemplo: La reaccin A P se lleva a cabo en un reactor continuo de tanque agitado de 2000 L de volumen. El reactor se alimenta con 300 L/min de una disolucin de A puro conteniendo 4 mol/L. Determine la conversin y temperatura de estado estacionario si el reactor opera adiabticamente y: T0 = 290 K, T0 = 298 K y T0 = 305 K Datos:
Cp = 3,5 J/gK

se plantean los balances de materia y energa:

B.M. FAo - FA = ( rA ) V

12000 = 1,15 g/cm3 Hr = - 50 kJ/mol k = 2,4 1015 exp min-1 T

c Ao - c A = ( rA ) t M

B.E.

& V Cp (TE - T ) + (- Hr ) r V + U S (Tm - T ) = 0

La resolucin del sistema de ecuaciones proporcionar las condiciones de operacin

14

RCTA NO ISOTRMICO
Ejemplo (cont.)

Para T0 = 290 K, hay una sola solucin: XA = 0,02 T = 291 K

UAA G X E,U X

1,0 1,0 0,9 0,9 0,8 0,8 0,7 0,7 0,6 0,6 0,5 0,5 0,4 0,4 0,3 0,3 0,2 0,2 0,1 0,1

Para T0 = 298 K, hay tres soluciones: XA = 0,08 T = 302 K XA = 0,42 T = 319 K XA = 0,91 T = 343 K

T0 = 305 K T0 = 290 K

Para T0 = 305 K, hay una sola solucin: XA = 0,96 T = 353 K

T0 = 298 K
290 290 300 300 310 310 320 320 330 330 340 340 350 350 360 360

0,0 0,0 280 280

T (K) T (K)

RCTA NO ISOTRMICO
a la vista de los resultados del ejemplo anterior cabe hacerse las siguientes preguntas:
en qu intervalo de temperaturas de alimentacin existen mltiples estados estacionarios? son todos ellos estables? cul de los posibles estados estacionarios tiene lugar en la operacin en estado estacionario?

15

RCTA: anlisis de estabilidad

en qu intervalo de temperatura de alimentacin existen mltiples estados estacionarios?
XA UE,UG XA
1,0 1,0 0,9 0,9 0,8 0,8 0,7 0,7 0,6 0,6 0,5 0,5 0,4 0,4 0,3 0,3 0,2 0,2 0,1 0,1
0 T0 0,0 0,0 280 290 280 290

B B

osi T0 est comprendida entre T0 y T0,se darn mltiples estados estacionarios osi T0 > T0 , se dar un nico estado estacionario en la zona de encendido

A A T T0 T0
300 310 300 TA 310 T
A

320 320

330 TB 340 330 TB 340

350 350

360 360

TT(K) (K)

osi T0 < T0 , se dar un nico estado estacionario en la zona de apagado

RCTA: anlisis de estabilidad

son todos los estados de operacin estables?
50

2 c 3

45 40 35 30 25 20 15 10 5 0 280

1
290 300

a
310 320 330 340 350 360

T (K)

16

RCTA: anlisis de estabilidad

son todos ellos estables?
situacin 1: la reaccin transcurre a T baja y con XA tambin pequea situacin 3: la reaccin transcurre a T alta y con XA elevada situacin 2: tres posibles estados estacionarios: a, b y c
si la temperatura aumenta alrededor del punto c, tanto UG como UE aumentan, por los que el sistema se enfriar hasta volver al punto c. Si la temperatura disminuye, UE < UG y la temperatura del sistema aumenta hasta volver de nuevo al punto c en el punto a la situacin es similar: si la temperatura aumenta, T disminuye hasta alcanzar de nuevo el estado estacionario y viceversa su T disminuye en el punto b, si la T aumenta, UG aumenta ms deprisa de lo que aumenta UE, de forma que la T seguir aumentando hasta el punto c. Si T disminuye, la T disminuir hasta el punto a

RCTA: anlisis de estabilidad

cul de los posibles estados estacionarios tiene lugar en la operacin en estado estacionario?
para cualquier valor de la temperatura de alimentacin entre T0 y T0 no es posible predecir cul de los dos estados estacionarios prevalece
El estado estacionario en la zona de encendido es preferible al estado estacionario en la zona de apagado, si se desean conversiones finales elevadas

para evitar la posibilidad de mltiples estados estacionarios, debe seleccionarse la temperatura de alimentacin para evitar stos y tambin deber ser mayor que la temperatura de ignicin (T0 ).

17

RCTA: Anlisis de estabilidad

Sensibilidad paramtrica:
Influencia de la temperatura de alimentacin, TE, en el estado o estados estacionarios del sistema de reaccin Influencia del flujo volumtrico de alimentacin en el estado o estados estacionarios del sistema de reaccin

RCTA: anlisis de estabilidad

Temperatura de alimentacin
Partiendo de las ecuaciones de las curvas UG y UE:
UG,, UE UG UE
50 50 7 7
8

40 40
6

30 30

20 20

10 10
1 2

24

0 0 275 275

3 3 300 300 325 325 350 350 375 375

T (K) T (K)

18

RCTA: anlisis de estabilidad

Temperatura de alimentacin
Partiendo de las ecuaciones de las curvas UG y UE: la representacin T vs T0 constituye la curva de ignicinapagado del reactor

RCTA: anlisis de estabilidad

Tiempo espacial
de forma anloga a la anterior es posible determinar los valores crticos para el tiempo espacial

UG =

(- Hr ) r V = t
& V Cp

Jr

UE = T (1 + m) - (TE + m Tm )

& & Si V = 0 UG V UG 0 & & Si V m V m

19

RCTA: anlisis de estabilidad

Tiempo espacial
de forma anloga a la anterior es posible determinar los valores crticos para el tiempo espacial

UG =

(- Hr ) r V = t
& V Cp

Jr

UE = T (1 + m) - (TE + m Tm )

& & Si V = 0 UG V UG 0 & & Si V m V m

Secuencia de reactores continuos de tanque agitado (No Isotrmico)

La forma de operar puede ser:

adiabtica: sin intercambio de energa con los alrededores (q = 0) no isotrmica no adiabtica: intercambio con un fluido (q)

Clculo: Para cada uno de los tanques (l) que conforman la batera se lleva a cabo la resolucin simultnea de los balances de materia y energa, proporcionando las condiciones de salida de de dicho tanque (cAl, Tl) que correspondern con las de entrada del tanque siguiente.

20

RCTA: Anlisis de Estabilidad (local)

Estabilidad de los estados estacionarios
Anlisis de bifurcacin
Condicin necesaria de unicidad

mximo de

d UG (T) <A dT

f (T ) = U E (T ) - U G (T )

f (T) = A T - B - UG (T) = 0 dUG (T) f =A=0 T T dT T

Condicin necesaria de estabilidad

dUE (T) dUG (T) > dT dT

RCTA: Anlisis de Estabilidad (local)

Comportamiento ante perturbaciones
Balances en estado no estacionario
B.M. dNA dt = FAo FA + rA V

dT & & B. E. V cp = V cp (TE T) + ( Hr ) r V + U S (Tm T) dt

Balances en estado estacionario

B.M. 0 = FAo FA + rA V B.E. & 0 = V cp (TE T) + ( Hr ) r V + U S (Tm T)

Restando miembro a miembro: dc A B.M. tM = (c A c As ) + A t M (r rs ) dt dT B.E. tM = (1 + m)(T Ts ) + t M J(r rs ) dt

21

RCTA: Anlisis de Estabilidad (local)

Comportamiento ante perturbaciones
siendo: x = C A - C AS y = T - TS es decir, t M dx = x + A t M (r rs )

dt dy tM = (1 + m) y + t MJ(r rs ) dt

linealizando: r = rS + rCAS x + rTS y (Taylor)

x 1 = ( ArcAs ) x + ( ArTs ) y t tM y (1 + m) = + JrTs y + (JrcAs ) x t tM

RCTA: Anlisis de Estabilidad (local)

Comportamiento ante perturbaciones
aplicando transformadas de Laplace y resolviendo: x = A 11 e s1 t + A 12 e s 2 t

y = A 21 e s1 t + A 22 e s 2 t
s 2 - (a + b ) s + (ab - c ) = 0

siendo s1 y s2 las races del polinomio:

s=

1 (a + b ) 2

(a - b )2 + 4 c

con:

a = rC

AS

1 tM

1 b = J rTs - (1 + m) tM c = J rTs rC
AS

r rC AS = C A

T = TS

r rT = T C A = C AS

las races pueden ser reales o imaginarias

22

RCTA: Anlisis de Estabilidad (local)

Comportamiento ante perturbaciones
races reales
s1 y s2 < 0: s1 s2 = 0: s1 y s2 > 0: desaparece la perturbacin con el tiempo las perturbaciones permanecen en un valor final cte. las perturbaciones crecen exponencialmente

RCTA: Anlisis de Estabilidad (local)

Comportamiento ante perturbaciones
races reales races imaginarias
parte real < 0: parte real = 0: parte real > 0: oscilacin amortiguada (perturbacin desaparece) oscilacin sostenida amplitud de las perturbaciones crece exponencialmente

23

RCTA: Anlisis de Estabilidad (local)

Comportamiento ante perturbaciones
races reales
ambas <o la perturbacin desaparece con el tiempo

races imaginarias
parte real < 0: oscilacin amortiguada (perturbacin desaparece)

debe cumplirse:

s1 + s 2 < 0

s1 s 2 > 0

a >b

ab > c

RCTA: Anlisis de Estabilidad (global)

Plano de fases
estudio de la estabilidad global del reactor resolucin de las ecuaciones T-t y CA-t en estado no estacionario para diferentes condiciones iniciales Ti y CAi (mtodo de Euler)

B.M. B.E.

dc A = (c A c As ) + A t M (r rs ) dt dT tM = (1 + m)(T Ts ) + t M J(r rs ) dt tM

24

RCTA: Anlisis de Estabilidad (global)

Plano de fases
estudio de la estabilidad global del reactor resolucin de las ecuaciones T-t y CA-t en estado no estacionario trayectorias en planos de fases: lnea separatriz
cuando se parte de diferentes condiciones de partida, se observa una aproximacin al estado estacionario TS en forma espiral, con diferente trayectoria para cada condicin inicial

RCTA: Anlisis de Estabilidad (global)

Plano de fases
estudio de la estabilidad global del reactor resolucin de las ecuaciones T-t y cA(XA)-t en estado no estacionario trayectorias en planos de fases: lnea separatriz
cuando se parte de diferentes condiciones de partida, se observa una aproximacin al estado estacionario TS en forma espiral, con diferente trayectoria para cada condicin inicial

25

RCTA: Anlisis de Estabilidad (global)

Plano de fases
estudio de la estabilidad global del reactor resolucin de las ecuaciones T-t y -t en estado no estacionario trayectorias en planos de fases: lnea separatriz
cuando se parte de diferentes condiciones de partida, se observa una aproximacin al estado estacionario TS en forma espiral, con diferente trayectoria para cada condicin inicial para un caso tpico con tres estados estacionarios A, B y C

26