18&19 Maria Gonzales

dQ g dQ e
> (35)
dT dT
En cambio, los A y C tienen la característica inversa. La (35) es condición de inestabilidad. En
cambio, excepto para el reactor adiabático, no es condición suficiente de estabilidad, coma se verá
enseguida.
8.7. ANALISIS DEL ESTADO NO ESTACIONARIO
PLANTEAMIENTO MATEMATICO
Analizaremos aquí el estado no estacionario o transitorio. Para ello usaremos las ecuaciones (10) y
(21) del capítulo 6:
dX (6-10)
θ =−x +θ r (T , X )
dt
dX (6-21)
θ =T 0−T + θ Jr ( T , X ) −Q(T )
dt
Simultáneamente escribiremos las ecuaciones de estado estacionario, para un determinado
estado estacionario caracterizado por T ee y X ee:
θ=−X ee +θ r (T ee , X ee) (36)
θ=T 0 −T ee +θ Jr ( T ee , X ee ) −Q(T ee) (37)
Restemos de las ecuaciones (6-10) y (6-21) las (36) y (37), respectivamente:
d (X −X ee ) (38)
θ =−(X −X ee )+θ {r ( T , X ) −r ( T ee , X ee ) }
dt
d ( X −X ee ) (39)
θ =−( T −T ee) +θ J { r ( T , X )−r ( T ee , X ee ) }−{Q (T )−Q ( T ee ) }
dt
Obsérvese que T ee y X ee no son funciones del tiempo y, por lo tanto, sus derivadas con respecto al
tiempo son nulas. Las (18) y (39) muestran cómo varían las perturbaciones.
x=X −X ee ; y=T −T ee (40)
con el tiempo. Estas ecuaciones valen para cualquier valor que tomen las perturbaciones. Pero
supongamos que las perturbaciones toman valores absolutos pequeños. Entonces, las partes no
lineales (entre llaves en las (38) y (39)) pueden ser linealizadas alrededor del punto de estado
estacionario, T ee, X ee, desarrollando por series de Taylor:
( dTdr ) ∗y +( dXdr ) ∗x
r ( T , X ) −r ( T ee , X ee ) =
ee ee
Q ( T ) −Q ( T ) =(
dT )
dQ
ee ∗y
ee
en donde el subíndice ee implica que las derivadas están evaluadas en el estado estacionario
considerado. Reemplazando en las (38) y (39) se tiene:
θ
dx
dt
=−Lx+
N
J
y ( ) (41)
dy
θ =J ( 1−L ) x−( M −N ) y (42)
dt
En donde
L=1−θ ( dTdr ) , M =1+( dQ

ee dT ) , N =Jθ ( )
dr
dT ee ee
(43)
Las (41) y (42) forman un sistema de dos ecuaciones lineales diferenciales ordinarias. En los libros
de matemáticas se demuestra que la solución de este tipo de sistema de ecuaciones diferenciales
tiene la siguiente forma:
t t
m1 m2
x ( t )=cte1 e θ
+ cte 2 e θ (44 a)
t t (44 b)
m1 m2
θ θ
x ( t )=cte3 e +cte4 e
en donde cte i , i = 1, 2, 3 y 4 son constantes de integración y m1 y m2 son las raíces de la ecuación

cuadrática:
m 2 + ( L+ M −N ) m+ ( LM −N )=0 (45)
Si cualquiera de las dos raíces, m1 o m2, tiene parte real positiva, alguna (o todas) de las
exponenciales de las (44), crecerá sin límite a medida que pasa el tiempo. En otras palabras, las
perturbaciones x e y se harán más y más grandes a medida que aumenta el tiempo. También
puede decirse que el reactor se alejará del estado estacionario, cada vez más. En consecuencia, el
estado estacionario será inestable.
Imponiendo a las raíces características, m1 y m2, la condición de que tengan su parte real
negativa, se obtiene
L+ M > N y> LM > N . (46)
Estas son las condiciones de estabilidad del estado estacionario, puesto que las exponenciales,
todas negativas, tenderán a cero y, con ellas, las perturbaciones.
La forma en que las perturbaciones varían con el tiempo muestra en la figura 16. Si m1 y m2 son
reales, se tienen evoluciones del tipo mostrado en la figura 16 (a). Las curvas A, B y C son casos
estables en que las perturbaciones tienden a cero con el transcurrir del tiempo. En cambio, las D y
E son casos en que el estado estacionario es inestable y las perturbaciones no cesan de crecer.
Este crecimiento en realidad se detiene cuando el reactor pasa a trabajar a otro estado
estacionario estable. Cuando las raíces son complejas, las perturbaciones son oscilatorias, como se
muestra en la figura 16 (b}. El caso A es uno estable y las oscilaciones se amortiguan hasta
anularse. Por el contrario, el caso B es inestable y las oscilaciones crecen en amplitud con el
tiempo. Si las raíces son imaginarias puras, el sistema oscilar permanentemente con amplitud
constante, como se muestra por trazos en la figura 16(b). En este caso la definici6n de estabilidad
o inestabilidad es ambigua pe-
FIGURA 16. Evolución de las perturbaciones en el tiempo

18&19 Maria Gonzales

Cargado por

Copyright:

Formatos disponibles

18&19 Maria Gonzales

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

18&19 Maria Gonzales

Cargado por

Copyright:

Formatos disponibles

dQ g dQ e

8.7. ANALISIS DEL ESTADO NO ESTACIONARIO

θ=−X ee +θ r (T ee , X ee) (36)

θ=T 0 −T ee +θ Jr ( T ee , X ee ) −Q(T ee) (37)

Restemos de las ecuaciones (6-10) y (6-21) las (36) y (37), respectivamente:

x=X −X ee ; y=T −T ee (40)

L=1−θ ( dTdr ) , M =1+( dQ

en donde cte i , i = 1, 2, 3 y 4 son constantes de integración y m1 y m2 son las raíces de la ecuación

L+ M > N y> LM > N . (46)

FIGURA 16. Evolución de las perturbaciones en el tiempo

También podría gustarte