CAPITULO

Comportamiento de procesos
• •
trntis üomoe

13.1 INTRODUCCION

Por lo general, se ha considerado que los modelos para los procesos que
se desarrollan en este texto dependen del tiempo. Es muy importante
comprender la forma en que dichos procesos dependen del tiempo, para
efectuar el diseño y el análisis de los experimentos y para saber cómo
operará un proceso en relación a las especificaciones de diseño para el
estado estable. Generalmente, se ha hecho énfasis en lo primero, como por
ejemplo, en los capítulos S y 8. Aquí se hace un estudio de un reactor tipo
tanque con agitación, en la operación real de un proceso que ha sido
diseñado para proporcionar un comportamiento específico. Particularrnen-
te, se considera el efecto de ,las variaciones en la alimentación sobre la
conversión del reactor y también el diseño de un sistema de control simple
para mantener una producción constante, aun cuando fluctúen las condi-
ciones de alimentación. Al realizar 10 anterior, se hace hincapié en los
aspectos de predicción y diseño del proceso del análisis, y se dan algunas
técnicas con el fin de determinar el comportamiento del modelo para
sistemas de segundo orden.

13.2 ECUACIONES DEL REACTOR

Las ecuaciones que describen la reacción de pseudoprimer orden
A --+ D en un reactor tipo tanque con agitación y con camisa, se dedujeron
en la sección 12.9 como

(j dCA k (j -E/RT
- = CAl - CA - o e cA. (13.1)
dt

(j dT
&= T -
lOA
T + Jk (je-E/RTc - u[l + K][T - Te] (13.2)

449
450 Sistemas no isotérmicos

(13.3)

Aquí se ha despreciado el trabajo a la flecha, y se ha definido () = Vlq,

J = -t1y R/P.E.p, K = ha/Pe%.E.pe, u = Pc%.E.pcK/[1 + K]PQ.E.1P en donde el
subíndice e se refiere a la camisa. Es conveniente efectuar una simplifica-
ción para disminuir al mínimo las manipulaciones, sin afectar serrarnente
la estructura básica del problema. Se supone que el tiempo de residencia ()c
de la camisa. es pequeño en comparación al del reactor (volumen pequeño o
alta velocidad de flujo), y formalmente se establece que ()c -- O. (En el
problema 13.1, se bosqueja una justificación rigurosa de este aspecto.) En
este caso, la ecuación 13.3 se simplifica. dando como resultado la siguiente
ecuación algebraica

(13.4)
Con ella se puede calcular T~ substituir en la ecuación 13.2. La ecuación
resultante para la temperatura del reactor es:

(13.5)

Las ecuaciones 13.1 y 13.5 unidas a las condiciones en que el tiempo t = 0,

determinan el comportamiento transitorio. Ya se ha estudiado la solución
de estas ecuaciones en la operación bajo condiciones constantes.
La diferencia entre las ecuaciones lineales y no lineales se estudia en las
secciones 17.1 Y 17.2. Las ecuaciones 1~.1 Y 13.5 incluyen el exponencial
e-F;/ RI', así como el prod ucto CAe-E/ IlT, y, por lo tanto, son no lineales.
Debido a lo anterior, es muy poco probable que puedan obtenerse las
soluciones en términos de funciones simples. Se puede emplear la integra-
ción numérica del tipo utilizado en la sección 12.11 para obtener respuestas
especificas para un conjunto dado de parámetros, pero usualmente se
requerirán expresiones analíticas y de validez general para construcciones
más amplias. Al final de la Seco 2.4 se pueden encontrar los indicios para
efectuar las aproximaciones necesarias que conducirán a un conjunto de
ecuaciones significativas que tengan soluciones analíticas simples. Ahí se
observó que las pequeñas desviaciones de un valor designado se podrían
aproximar bastante. mediante una función lineal. Este concepto tiene
aplicaciones muy útiles.
Es necesario hacer una concesión fisica, con objeto de obtener un
problema matemático de fácil manipulación y, ésta consiste en considerar
solamente las temperaturas y concentraciones que estén cercanas a las
condiciones constantes de diseño. Posteriormente, se tratará el enunciado
 Comportamiento de procesos transitorios 4S1

matemático preciso de esta suposición y su necesidad. En el estado estable

las ecuaciones 13.1 y 13.5 se convierten en:

O = CAfs - c..1s - k o 8 e -E/RT.. CAs

(13 . 6)
-E/Rl',
O = Tf,; - T, + Jko8e CAs - u.jTs - TcfJ (13.7)

Para comprender el fenómeno físico, simplificando cuanto sea posible el

tratamiento matemático, es conveniente considerar que las condiciones
co ns tan tes de operación son el punto de referencia, y medir las cantidades
en términos de su diferencia a partir de dicho valor. De esta forma,
haciendo que x(t) sea la diferencia entre la concentración de A en cualquier
tiempo y el valor bajo condiciones constantes:
(13.8a)
en forma similar, si y(t) es la diferencia de temperatura:

y(t) = T(t) - T, (13.8b)

Es preciso considerar el caso en que puedan variar la concentración de la

alimentación y la temperatura, de tal manera, que las fluctuaciones se
definan como xf(t) y Yf(t), respectivamente,
•
(I3.8c)

(l3.8d)

Para no complicar el desarrollo, se supondrá que las velocidades de flujo q

y qc permanecen constantes, y también la temperatura de alimentación del
refrigerante, TCf' En otra sección se utilizarán estas suposiciones.
Las ecuaciones 13.1 Y 13.5 pueden escribirse en términos de las variables
de tluctuación x y y como:

O
d[ct<
.'
+ x] = 0-
d»
dt dt
= e•..t.fs + xf - c.I" - X - koOe-E/R[Ts+!}J[C_1S + x]
(13.9)

o d[T< + ?J] D dy
dt el!
- T Ls + '1)
• f
- T .<. - 7,/ + J k.o De-J.·'/ R,[Ts+yJ[c, .t» + x]

(13.10)
452 Sistemas no isotérmicos

Se sabe que c_l.~y T', son constantes, de tal manera, que dcAJdt =dTjdt = O.
Estas ecuaciones comprenden un número de términos de la operación bajo
condiciones constantes, y pueden ordenarse de nuevo y simplificarse al
substituir las ecu acrones 13.6 y 13.7 para el estado estable. Esto da:

(13.11)

(13.12)

En este aspecto se utilizará la idea de estar muy cercanos a la operación

bajo condiciones constantes y que, por lo tanto, x y y son de pequeña
magnitud. El exponencial puede expresarse en términos de una serie de
Taylor (ecuación 1S.20), alrededor de y = O como:

Entonces la ecuación 13.11 es

ko ec As
Ee-E1 RT.
RT2 s Y

(13.13)

Por supuesto, la ecuación 13.13 es equivalente completamente a la ecuación

or igin al de balance de masa, 13.1. que, sin embargo, está en una forma más
 Comportamiento de procesos transitorios 453

apropiada para efectuar un cálculo aproximado. Se supone

que x y y son
sujicienremente pequeñas en magnitud. de tal manera que y2, xy, etc; son
despreciables en comparación con x y y. Se pueden obtener resultados muy
útiles mediante esta aproximación, pero ¡hay que notar que implica una
gran restricción física! No obstante, ahora la ecuación 13.13 se convierte en:

dx k Bc Ee- EIRT•
a _
u
= XI - X -
In-El
l(oUe
RT.
X -
-.::.o---,.~.J...::..s
R T 2
_
Y (13.14a)
dt s

De manera semejante, la ecuación de temperatura 13.12 se convierte en:

d a E -EIRT.
n ~
u
__. YI _ y + Jk e -EIRT Sx
o e + Jk OuCAs 2
e y - u!J
dt RT s

(13.15a)

o reagrupando:

dx [k Bc Ee-EIRTS]
B- = - [1 +k ee-EIRTs]x _ o As y + X
dt o RT2 I
s

(13.14b)

dy
e - = + [Jko8e-EIRTs]x -
[Jk
1 +u _ o
ec
As
Ee-EIRTS]
y +y
dt RT2s I

(13.15b)

Estas ecu aciones son ahora lineales, y los coeficientes de x y y son

constantes conocidas. cuyos valores dependen. del estado constante muy
cerca del cual van a efectuarse los cálculos. Limitándose a las concentra-
ciones y temperaturas cercanas al estado estable, se han obtenido las
ecuaciones que pueden tratarse analíticamente mediante los métodos
bosquejados en el capítulo 17.
Para su aplicación, es conveniente convertir las dos ecuaciones lineales
diferenciales de primer orden a una ecuación única de segundo orden, ya
sea para la concentración o la temperatura. Las ecuaciones 13.14 y 13.15
son de la forma de las ecuaciones 17.13. En la sección 17.10 se demuestra
que estas ecuaciones se pueden reemplazar con cualquiera de las siguientes
ecuaciones equivalentes de segundo orden:
454 Sistemas 110 isotérmicos

(13.16)

(13.17)

(13.18a)

(13.18b)

j;() __ 1 1+ _ Jk e
OCAs E e -EIRTS]
1 t - u XI
e2 [ Ry2s

1 dx / k C Ee-EIRTs
+ __ ---'-0 --=..-'=1-' _
(13.19a)
() dt - 8Ry2 y/
s

{',,(t) =
J'.. e2 [1 +
-.L k
o
ee-EIRT~]1j
. r
+ 1e dYI'
dt
+ 2e J k o e-E1RT"x I

(13.19b)

Para los parárnetros del reactor utilizados en los ejemplos 16.1 y 16.4, los
valores de al' ao, 11(1) Y 12(t) de cada uno de los tres estados estables son
como sigue:

CAs = 0.499, Ts=311.4: a1=6.57 a¿ = 10.80

dx , -4
fl(r) = 10.78x/ + 3.33 - - 9.92 x 10 YI
dt
ely f
f"l(r) = 11.13!J/ + 3.33 -
tlr
+ 7.01 x¡
e.ls = 0.257, T, = 393.4: al = -6.71 00 = -33.46
dx,
flet) = -44.00x¡ + 3.33 -
elt
- 0.16321/

d u¡ '~
/2(t) = 21.6521/ + 3.33 -' + 3.56 x 10 x¡
elr

C.Js = 0.008, 'T, = 477.5: a1=188.3 00 = 616.6

 Comportamiento de procesos transitorios 455

dx,
fl(t) = - 64.11 xf + 3.33 - - 0.222y f
dt

13.3 OBTENCION DEL ESTADO ESTABLE

En la sección 12.10 se comenzó a investigar la cuestión de si sería posible

alcanzar las condiciones de operación en estado estable calculado previa-
mente. Se desarrolló un criterio gráfico para la imposibilidad de alcanzar-
las. y ahora se está en posición de tratar totalmente con el problema. Se
supone que las perturbaciones en la alimentación han desplazado la
concentración y la temperatura en una cantidad Xo y Yo, respectivamente,

a partir de las condiciones -de diseño para estado estable. Para t > O se
puede suponer, para simplificar, que no hay perturbaciones posteriores
(Xf = Yf = O). Entonces la ecuación 13.16 que describe la respuesta del
reactor se convierte en:

(13.20)

Si al pasar el tiempo x(t) disminuye- hasta cero con el tiempo, el diseño bajo
condiciones constantes es estable y puede mantenerse. Si x{t) no regresa a
cero con el tiempo. no puede realizarse la operación en condiciones
constantes.
La ecuación 13.20 es una ecuación lineal de segundo orden con
coeficientes constantes. El procedimiento de resolución se bosqueja en la
sección 17.3. Se supone una solución de la forma:

x = ce='

Substituyendo en la ecuación 13.20, se tiene:

Ce [m2 lHt
+ a.m + aol = O
y. dado que el exponencial no desaparece, se obtiene la ecuación ca-
racter is tica:
111.
2
+ al171 + ao = O (13.21)
Esta tiene las raíces m¡ Y 1712 dadas por:

111, ~ - ~' + [a~'- ao)''' (13.22a)

al'!. J1/2
/11'2 = [ 4 - a)
( (13.22b)
2
456 Sistemas no isotérmicos

Note que:

Inl + 1112 = -al (13.23)

111lm2 = +ao (13.24)
La solución general es:
x(t) = Cle1lt1t + C2em2t (13.25)
El comportamiento de la concentración y la temperatura después del
tiempo cero es regido completamente por el carácter de mI Y m2.
Inicialmente, es conveniente distinguir entre los casos cuando mI Y m ¿ son
reales. ' y cuando éstos son complejos, dependiendo del signo del discri-
minante. a12 - 4ao.
Caso I~ a12 - 4ao > O, mI Y m2 son reales.
En este caso, existen tres posibilidades: dos raíces positivas, dos raíces
negativas, una raíz positiva y una negativa. El término e= t tenderá a cero si
m es negativo, y crecerá sin límite si m es positivo. De esta forma, tanto mI
como 171'2 son negativos, x(t) irá hasta cero, como se muestra en la figura
13.1. Si 1111 o m2 son positivos, la solución crecerá hasta ± ex), según se
indica en la figura. La operación en condiciones constantes es estable y
puede mantenerse a pesar de las pequeñas perturbaciones, si y sólo si 1111 y
1112 son negativos. Es obvio que, cuando una o ambas raíces son positivas, la

conversión en realidad no irá hasta infinito. Después de un corto periodo de

tiempo, la suposición de que x2 es pequeño en comparación con x, no será
válida, y la ecuación 13.20 no podrá describir la .situación física.
Las ecuaciories 13.23 y 13.24 proporcionan un criterio sencillo para
determinar cuándo son negativas ambas raíces. La suma de dos números
negativos es negativa, de tal manera, que a partir de la ecuación 13.23, -al
deberá ser negativo, o al > O. El producto de dos números negativos es
positivo, de tal manera, que a partir de la ecuación 13.24 se tiene a¿ > O. A
partir de la ecuación 13.18 éstos pueden expresarse en términos de los
parámetros del reactor como:

_ Jc A,~E] >O (13.26a)

RT2s

_ JC .., .• EJ >0 (13.26b)

2
RT s

Note que para el reactor adiabático. u = 0, la segunda condición lleva

implícita la primera, dado que ea1 = (Pao + l. En general, esto no se
aplica.
 Comportamiento de procesos transitorios 457

Una positiva.
una negativa

X (t)

t
FIGURA 13.1 Respuesta transitoria de un sistema lineal homogéneo de segundo
orden con raíces características reales a12 - 4a~ > o. x(t) tiende a cero, si y sólo si
ambas raíces son negativas.

Caso 11: a12 - 4ao < 0, mI Y m2 son complejos.

Este caso puede ocurrir solamente, si a¿ > O.. m¡ Y 1112 son ahora

J
2
al. al
1112 = - 2- 1 ao - 4
donde i2 = -l. Entonces, la ecuación 13.25 se convierte en:

Siguiendo la sección 17.4 se puede utilizar la igualdad:

y escribir otra vez x(t) en la forma:

x(t) =
-a t/2
e 1 [Kl cos
J ao - 4 t + K sen
al
2
2
J ao - 4t
al
2 J (13.27)

Aquí KI = CI + C2, K2 = i[CI - C2]. Entonces, si al > 0, la respuesta es

amortiguada, según se demuestra en la figura 13.2(a). Si al < la fluctua- °
ción de la concentración crece a partir de cero en forma osci1atoria, según
se muestra en la figura 13.2(b).
 458 Sistemas no isotérmicos

.x (t) }---+----I-----\---.,---- ."C (t ) t----\----+------\----r---

t t

w w
He; u HA I 3.2 Respuesta tra ns itoria de un sistema lineal homogéneo de segundo
2
orden con \'a'íccs ca ructet-ixticas complejas. al - 4ao < O. x(t) tiende a cero en
forma oxc iln toriu. si y sólo si. ambas raíces tienen partes negativas reales
L(a) al > O], Y a u menta en forma oscilatoria si existe una parte positiva real
[(b), al < O].

De esta forma, el diseño en estado estable se alcanza, si y sólo si,al > O.

Dado que el comportamiento oscilatorio (raíces complejas) requiere que
ao > 0, conduce a las mismas condiciones que el caso T,al > 0, ao > 0, que
se expresa en términos de los parámetros del reactor en las ecuaciones 13.26.

13.4 CONDICION DE LA PENDIENTE •

En la sección 12.10 se utilizaron argumentos especiales para concluir que

el estado estable no se puede alcanzar en la práctica, siempre que la
pend iente de la línea I de la figura 12.7 sea más negativa que la de la línea
11. Esta condición de pendiente puede relacionarse con el criterio seguido
en la sección anterior, según se indica:
La línea I es una gráfica de c.L respecto a T definida por la ecuación 12.55

C.lr
C..:LI = 1 + (12.55)
kof)e-l~'/RT

Se ha utilizado el subíndice 1 para expresar "a partir de la línea 1". La

pendiente se forma a partir de la diferenciación con respecto a T

"
dCA 1 k o 8e-E/RTc ..:ir E

dT

En estado estable CAl = CAs' T _.- Ts,'y la ecuación 12.55 puede substituirse
para obtener la pendiente en la intersección de las líneas 1 y 11:
 Comportamiento de procesos transitorios 459

klJE CAse-EIRT.
Estado estable: (13.28)
dT RTs21 + kOee-EIRT6

La línea 11 es una línea recta definida por la ecuación 12.54,

CAII = C.A.I + .! [T, + uTct] _ 1 +u T (12.54)

J J

Entonces la pendiente es:

1+ u
(13.29}
dT J

Al comparar las ecuaciones 13.28 y 13.29 con la ecuacion 13.26b, se

demuestra que la condición a¿ > O es idéntica a la condición de que
dCAI/dT> dCAIT/dT en estado estable. De esta forma, las condiciones
constantes medias deben tener a¿ < O y, nunca pueden alcanzarse en la
práctica. Las condiciones de operación constante a alta y baja temperatura
siempre tienen a¿ > O. A menos, que el reactor sea adiabático (u = O)
siempre es necesario comprobar que al > O para asegurar que la operación
bajo condiciones constantes es estable. Es posible construir un reactor con
condiciones Go > O, l?l < O. Un ejemplo de esto, es la reacción entre cloruro
de acetilo yagua para formar ácido acético y ácido clorhídrico en un
reactor de tlujo continuo tomado de una publicación de Baccaro, Gaitonde,
y Douglas en el AIChE Journal. 16, 249 (1970.).

13.5 CONDICION DE LA RESPUESTA

Las ecuaciones lineales simples del modelo que describen el comporta-

miento del reactor cuando las concentraciones y las temperaturas son
cercanas al estado estable, permiten calcular rápidamente los efectos de los
cambios en las condiciones de operación. Suponga, por ejemplo, que el
proceso está operando en estado estable cuando la temperatura de alimen-
tación cambia en una cantidad constante Yt. Entonces puede calcularse el
cambio en la concentración del efluente.
A partir de las ecuaciones 13.16 y 13.19, el cambio en la concentración
.x(t) = cA - cAs se describe mediante la ecuación

k C Ee-E1RT•
o VI = constante
As (13.30)
eRT'2e
Se ha tijado XI = O en la función h (1). La solución a esta ecuación no
homogénea se obtiene, utilizando los métodos bosquejados en la sección
 460 Sistemas no isotérmicos

17.7. La solución general es una suma de una solución homogénea y una

particular

X(t) = xh(t) + xp(t)

En la sección 13.3 se estudió la solución homogénea y es ésta:

xh(t) = Clem1t + C2e m2t

donde mI Y m2 están dadas por la ecuación 13.22. La solución particular

puede encontrarse del método de coeficientes no determinados en el
ejemplo 17.12 y es

kOC.:1sEe-EIRT.y I
xp = aof1RTs
2

De esta forma, la solución general es la siguiente:

k o eA.~Ee-E1RT.y I
(13.31)
aofJRTs2

Las constantes el y C2 se determinan a partir de las condiciones

iniciales. Cuando t = O el reactor está operando en estado estable, de tal
manera que x(O) = y(O) = O. A partir de la ecuación 13.14 es evidente que •
dxjdt = O a t = O cuando x(O), y(O), y x, igual a cero. De esta forma:

k o eAs Ee-E1RT'1J .':JI = O

x(O) = Cl + C2 -
aof1RTs 2

Después de calcular Cl y C2 se obtiene:

k o eAs Ee-E1RT.y I

aofJRTs2

1111 e11l2t] (
13.32)
m1 - m¿

El cambio de concentración con respecto al tiempo, para mI Y m2 reales se

muestra en la figura 13.3(a), o en la figura 13.3(b) cuando mI Y m2 son
complejos.
 Comportamiento de procesos transitorios 461

-EIRTs
- k o cAS E e Yf
- - - - - - - :=-::;:--.------
ao 8 RT~

x tt )

t
(a)

-EIRTs
-ko CA S Ee Yf

ao 8 RT;

xtt )

t
(b)

FIGUI{/\ 13.3 Cambio en la concentración de un RTTAFC calculada a partir de

las ccuacioncs linealizadas para un cambio constante en la temperatura de
alimentación. (a) Raíces características reales .•(b) Raíces características complejas.

13.6 ANALOGIAS

En numerosas ocasiones se ha indicado que se pueden describir

diferentes fenómenos fisicos mediante la misma ecuación matemática,
siempre y cuando, se asigne el significado adecuado a los' coeficientes de
dicha ecuación. Por ejemplo, los modelos de transferencia de masa y trans-
ferencia de calor entre fases tienen la misma forma matemática. A menudo,
es posible aprovechar dichas semejanzas en una estructura matemática,
462 Sistemas no isotérmicos

utilizando el conocimiento de un fenómeno fis ico para inferir el comporta-

miento de otro, y estas analogías juegan un papel muy importante en el
diseño de ingeniería.
J{x
/\1\1\1\1\
F (t)
/1/

<Ix
fJ.-
dt

FI(; U RA IJA Sis tc m a mu su-resortc-cá rnuru de a rre. Fuerza del resorte = - Kx,
lúer/.~I ,'iscos~1 de vucia m icruo = -,Il d.cl dt .

Para esto, no es necesario restringir el estudio solamente a fenómenos que

ocurren en situaciones de ingeniería química. Por ejemplo, considere el
sistema mecánico que se muestra en la figura 13.4. A partir de los
principios estudiados en los cursos básicos de física. es evidente que el
desplazamiento de la masa desde el equilibrio, x, se determina mediante la
ecuación:
d?» f.1., dx K 1
dt?
+--+
111 d t
x - F(t)' (13.33)
111 m
K es la constante del resorte, !l es el coeficiente de viscosidad de
amortiguación, In la masa, y Ftt) la fuerza externa aplicada a la masa. La
ecuación 13.33 es idéntica en forma a la ecuación }O3.16para la concentra-
ción del reactor. K / l1l corresponde a ao. Note que la constante del resorte
debe ser positiva para que el sistema pueda regresar al equilibrio. La suma
de energía potencial y energía cinética de la masa se puede expresar como:

E =--
111 [ti .1:J "!. 1Kx2
+"¡)
2. d t -
Entonces

1 .u: d » cI"2:~:
---+--:t:=
K dx _ !:!...[clX].!
J71 dt d t dt? /11 dt 1// dt

La energía de la masa debe disminuir hasta cero a medida que se aproxima

al equilibrio. de tal manera, que la condición de estabilidad requiere que
d Eldt < 0, p/m > O. De esta forma. para el sistema mecánico se tiene
ó

al > O. Por supuesto. la conservación de la energía requiere que el

amortiguador mecánico se caliente por el trabajo friccional. "
La ecuación para la caída del voltaje x a través del ca pacitor para el
 Comportamiento de procesos transitorios 463

circuito resistencia-inductancia-capacitancia, que se muestra en la figura

13.5 es la siguiente:
d2x 1 dx 1 1 dg
-+--+-x=-- (13.34)
dt2 RC dt Le RC dt
R es la resistencia. L la inductancia, e la capacitancia, y <ff' la fuente de l
voltaje. También esto es idéntico a la ecuación del reactor. Las analogías

+ + e
L x{

f'IGURA 13.5 Circuito R-L-C

son particularmente convenientes, puesto que construyendo un circuito

eléctrico con una ecuación correspondiente a otra situación fisica, es
posible seguir la variación del voltaje del circuito con el tiempo y
determinar el comportamiento del proceso en cuestión. Este es el principio
básico por el cual operan las com pu tadoras analógicas. Es necesario
relacionar los parámetros de circuito eléctrico con los parámetros análogos
del proceso real. También, se conocen computadoras analógicas mecánicas
y de fluido. aunque se emplean con menor frecuencia.

13.7 DISEÑO DEL SISTEMA DE CONTROL

En la sección 13.5 se calculó el cambio en la concentración del efluente

cuando ocurre un cambio en la temperatura de alimentación. En la mayoría
de los procesos, es deseable mantener la corriente de salida tan constante
como sea posible a pesar de las fluctuaciones en las entradas. Para esto, es
necesario diseñar sistemas de control que cambian ciertas condiciones de
operación para mantener la salida deseada. En la sección 4.3 se examinó el
problema de diseñar un controlador de nivel para un tanque. Aquí se
considerará el diseño del sistema de control para un reactor de tanque
agitado.
Las ecuaciones básicas del reactor son 13.1 y 13.5, Yse repiten aquí:
 464 Sistemas no isotérmicos

dCA
() -- k () -E/RT
- c_'u - eA. - o e cA (13.1) .
dt

dT -E/RT
()- = T, - T + Jko()e cA - u[T - Te,] ( 13.5)
dt

Se supone que la concentración o la temperatura de la corriente de

alimentación puede variar, pero se continuará suponiendo para simplificar,
que Te, es constante. El primer problema consiste en saber qué variable
debe medirse. Es obvio, que se desea mantener CA constante, pero es difícil
medir las concentraciones rápidamente y con exactitud. Por consiguiente, en
la mayoría de los casos habrá que medir la temperatura del reactor, y de
ella deducir las desviaciones en la concentración.
A continuación, se necesita establecer la cantidad que va a manipularse.
En este caso, probablemente se variará la velocidad de flujo de enfriamiento
qe o. en forma equivalente, u. Cuando la concentración CA sea demasiado
baja, se aumentará u, 10 que aumenta la velocidad de enfriamiento y por
consiguiente, disminuye la temperatura del. reactor y la velocidad de
reacción. Esto tenderá a compensar el efecto aumentando nuevamente cA
En forma similar, si CA se hace muy alta, se disminuirá u. En la figura 13.6
se muestra esquemáticamente el sistema de control.

Control de flujo
Fe t--------,

Po.sicionodor
L..-....---J neum6tico

V61vula de Sensor de
controi temperatura

FIGURA 13.6 Sistema de control de retroalimentación a un RTTAFC. La

velocidad de flujo de enfriamiento se ajusta midiendo la temperatura del reactor.

Las ecuaciones del diseño se derivan por Iinealización, exactamente en la

misma forma que se utilizó para obtener las ecuaciones 13.14 y 13.15. Sea
x = CA - CAs' x, =CA, - CA/s' Y = T - Ts' Y¡ = TI - TIs'y defina una
nueva variable l~ = U - us,donde us se calcula en condiciones constantes,
 Comportamiento de procesos transitorios 465

Cuando estas variables se substituyen en las ecuaciones 13.1 y 13.5 se

obtiene, después de lineálizar,·
dx 1
()_ = _ [" + k o ()e .
-El RT S] x - [k(}C
o As
EcmR~J + x,
2 Y
dt RTs
(13.35)
1 Jko(}C.AsEe-EIRTj
[ + u, - RT2
s
Y

- [Ts - Tc,]w + y, (13.36)

La ecuación 13.35 es idéntica a la 13.1, mientras que la ecuación 13.36
difiere de la 13.5, solamente por el término adicional - [Ts - Te,]w.
La estrategia de control consiste en aumentar u cuando se desea
disminuir la temperatura. Por consiguiente, w y y deben tener el mismo
signo algebraico. El esquema más simple es el control proporcional, en
donde el cambio de u se considera proporcional al cambio en T.
w=Ky (13.37)

De esta forma, cuando la temperatura varia considerablemente, el

esfuerzo para controlar es grande, y para pequeños cambios en temperatu-
ra, es pequeño. Note que la acción de control no se basa en la observación
directa de la fluctuación de la concentración, que es la cantidad en que
realmente se está interesado. Al substituir la ecuación 13.37 en la 13.36 se
tiene la ecuación final de trabajo para la temperatura:
o dy = [Jko(}e-EIRTa]x
dt
Jk (}Ee-EIRT,,}
- {1 + Us + K[Ts - Te,]·- O
RTs
2 Y + y,
(13.38)
Combinando las ecuaciones 13.35 y 13.38 en una sola ecuación de
segundo orden para la concentración, se obtiene

(13.39)
 466 Sistemas no isotérmicos

al' ao, y fl(l) son definidas por las ecuaciones 13.18 y 13.19.
El efecto principal del sistema de control, es aumentar los coeficientes de
la ecuación característica y, por consiguiente, hacer las raíces más
negativas. De esta forma, se acelera la respuesta transitoria. Reconside-
rando el problema relativo a la respuesta parcial de la sección 13.5, puede
obtenerse una medida cuantitativa del efecto del control. En ese caso, para
xt = 0, Yt = constante, la solución análoga a la ecuación 13.32 es:
- k e Ee-E/RT$y
x(t) = o As t

{ao + ;2 K[Ts - Tct][l + ko8e-E/RTS]}RTs2

x [1 + M MI - M2
2
eMl.t

MIY M2 son raíces de la ecuación característica:

(13.41)

x(t)

FIGURA 13.7 Cambio en la concentración de un H.TTAFC con raíces característi-

cas reales para un cambio gradual constante en la temperatura de alimentación a
diversos niveles de control proporcional de retroalimentación.

En la figura 13.7 se muestran las respuestas transitorias para diversos

valores de la ganan cia de retroalimentación K. La selección de K puede
hacerse sobre la base del valor máximo esperado de IYtl y el valor mínimo
permitido de Ix(t)I, en forma similar a como se realizó en el capitulo 4.
Note que cuando k es suficientemente grande, el estado estable final x( (0)
calculado de la ecuación 13.41 puede hacerse tan cercano a cero como se
desee.
 Comportamiento de procesos transitorios 467

También es interesante comparar la respuesta de la concentración del

etluente con un cambio en la temperatura de alimentación, calculada
anteriormente, con la respuesta a un cambio brusco en la concentración de
alimentación. En este caso, Yf = O, xf = constante y:
1{ J k Bc Ee-E.'/ RTS}
FI(t) = --; 1 + u + K[Ts - Te!] - O As 2 Xf
e- RT s

La solución homogénea es la misma, C1eM1t + C2eM2t, y la solución

particular se determina corno antes. En el tiempo cero se tiene x(O) =
y(O) = O. La condición inicial para dxjdt se obtiene con la ecuación 13.35.

dx 1
t = o: - (O) = - x
dt f e
Utilizando estas condiciones iniciales para evaluar las constantes, se
encuentra que la solución es:

x ( t) {1 + II + K[T. - Te!] - [JkofJcAsEe-E/RTs/RT/]}Xf

= -'--------------=..:'-'----=-----='---..:...:..::----~-..::......::~~
Q oe
2
+ K[Ts - Tef][I + koee-E/RT.]

X {l + M2
MI - M2
e"lf1t _.
MI -
MI
M2
eM2t}

(13.42)

En este caso, K aparece tanto en el numerador corno en el denominador, de

manera quex(t)nunca puede dirigirse a cero para la constante Xf diferente
de cero. Para K grande se obtiene:

Por consiguiente, no hay ninguna ventaja en hacer la ganancia excesiva-

mente grande. Esta imposibilidad de eliminar completamente una fluctua-
ción en la concentración de alimentación, es una consecuencia de la estruc-
tura que se ha adoptado para el sistema de control. Más aún, toda la
estrategia está basada en las temperaturas, a pesar de que la cantidad que
realmente interesa es la concentración.
Si se considera que w tiene una dependencia más general de y que la de
la ecuación 13.37, de todas formas se conserva la estructura lineal de las
ecuaciones, y la capacidad de diseñar dispositivos eléctricos y mecánicos
simples para efectuar el control. La mayoría de los sistemas de control
utilizados en las aplicaciones de proceso entran en la categoría general de
controles de tres form as
 468 Sistemas no isotérmicos

W = K lY
dt
i
+ K'!. -dy + K 3
o
t
Y( T) d-r (13.43)

En los problemas se exponen las propiedades que introducen las formas de

derivada e integral. Debe aclararse que el ejemplo del reactor que se utiliza
aquí, es único en un aspecto importante. Aunque puede no ser útil
económicamente, lograr que la ganancia proporcional K sea excesivamente
grande. no hay nada perjudicial en ello. Para situaciones distintas a la que
se considera aquí, es posible hacer un sistema inestable considerando una
ganancia demasiado grande. Por consiguiente, habrá un límite superior en
los parámetros de control. Esto, también, se deja a los problemas.

13.8 CONCLUSIONES

En este capítulo no se introducen principios físicos nuevos. Sin embargo,

el capítulo es una prueba para la capacidad del lector para leer matemáti-
cas. Frecuentemente, se hace referencia al capítulo 17 sobre ecuaciones di-
ferenciales. y debe consultarse, o bien consultar un texto sobre ecuaciones
diferenciales si hay dificultad. Fuera del uso de las ecuaciones lineales dife-
renciales de segundo orden para expresar el comportamiento del modelo, el
nuevo concepto es el de linealizar las ecuaciones del modelo para operar lo
más cerca posible al estado estable.
Este capítulo contiene una introducción a lo que se conoce generalmente
como dinámica del proceso y control. La mayoría de los textos en esta área
suponen que el lector es capaz de desarrollar un modelo, y enfatizan las
consecuencias de la manipulación matemática de las ecuaciones del
modelo. Algunos de los textos básicos escritos para ingenieros químicos
son:
13.1 P. S. Buckley. Tech niqu es of Process Control. Wiley. New York,
1964.
13.2 D. R. Coughanowr and L. B. Koppel, Process Systems Arialysis and
Con t rol. McGtaw-Hill, New York. 1965.
13.3 J. M. Douglas, Process Dynamics and Control. Prentice-Hall, En-
glewood Clitfs, 1972.
13.4 P. Harriott. Process Control. McGraw-Hill, New York, 1964.
13.5 D. D. Per lmuttcr, Lntroduction to Process Control. Wiley, New
York. 1965.
En la sección 22 de Ch ernical Engineers' Handbook hay un estudio bastante
extenso sobre instrumentación del control.
13.6 J. H. Perry, Ch emical Engineers' Handbook, 4th ed., McGraw-Hill,
New York, 1963.
Las analogías mencionadas en la sección 13.6 son la base del estudio en la
obra:
 Comportamiento de procesos transitorios 469

13.7 J. L. Shearer, A. T. Murphy, and H. H. Richardson, Introduction to

System Dynamics. Addison- Wesley, Reading, Mass., 1967.
En la obra siguiente se trata la aplicación a los problemas fisiológicos:
13.8 J. H. Milsum, Biological Control Systems A'nalysis. McGraw-Hill,
New York, 1969.
En los primeros capítulos de la obra siguiente, se encuentra una introduc-
ción a la teoría del control óptimo:
13.9 M. M. Denn, Optimization by Variational Methods, McGraw-Hill,
New York, 1969.
Las implicaciones de la retroalimentación se tratan en términos cualitativos
y algunas veces filosóficos, en una colección de artículos de la revista
Scientitic American que se ha publicado como libro:
13.10 Automatic Control, Simon And Schuster, New York, 1955.

13.9 PROBLEMAS

13.1 Al pasar de la ecuación 13.3 a la 13.4, se indicó que cuando ()c es

pequeña en comparación con (), debe despreciarse el término
0c[dTc/dt] en la ecuación para T'; Este problema conduce a los
pasos necesarios para demostrar esto en forma rigurosa.
(a) El "tiempo natural" para la respuesta de temperatura y concen-
tración es (). Sea T = t/(J, e = ()c/()[l + K]. Demuestre que la
•
solución a la ecuación 13.3 es

+ T[t - +
Ti T)'= Tcoe-Tlc.
i O e
e-Ale
JKT(T
1
-
+
A)
K
T el dA

TeO es el valor de Te en el tiempo cero. Es posible. que desee

referirse a la sección 15.10.
(b) Considere e un número pequeño, por ejemplo = 0.05. Grafi- 1;":

que [l/e]e-A//: respecto a A. Note que la función es diferente de

cero. solamente para A entre cero y Se. En forma similar. para T
mayor que un valor Se, e-TIC. ~ O. .
(e) Ahora puede escribirse:

5C.[1- e-A/c. ] +
Te( T) ,-.....J
1
O e
KT(T
1
-
+
A)
K
T el dA

Demuestre la expresión T(T - }.) ,-.....J T(T) para O< }.< Se. Por
consiguiente, escriba:
470 Sistemas no isotérmicos

Tc"""" KT
l+K
+ -r o
-1 e-)./e dA
e .

(d)Lleve a cabo la integración, y demuestre que conduce a la

ecuación 13.4.
13.2 El crecimiento intermitente de poblaciones microbianas puede ser
descrito, aproximadamente, mediante la ecuación:
dx
- = k1x(1 - k2x)
dt
donde x es la masa total del bio-organismo.
(a) Demuestre que la ecuación puede ponerse en la forma dyjdT
= y[1 - y] y calcule y( T). Grafique la solución para las
condiciones iniciales:
y(O) = 0, 0.05, 0.50, 0.95, 1.00
(b) ¿Cuáles son los estados estables del sistema? ¿Son estables o
inestables? Dé una explicación desde el punto de vista físico.
(e) Linealice la ecuación para y en cada uno de los estados estables
y resuélvala. Calcule cada solución linealizada junto con la solu-
ción exacta para las tres condiciones iniciales intermedias en la
parte (b). ¿Qué conclusiones se pueden obtener?
13.3 Considere un reactor tipo tanque agitado de flujo continuo adiabá-
rico descrito por las ecuaciones 13.1 y 13.2 con u = O.
(a) Demuestre que cuando las condiciones de (= O satisfacen
la relación especial
JCA(O) + T(O) = JCAf + r,
se obtiene una sola ecuación diferencial de primer orden para T.
(b) Manipule la ecuación para T, transformándola a:
rT(t) dT
t = JT(O) F(T)
Utilizando los parámetros del reactor del ejemplo 16.1, aplique
la regla del trapezoide (sección 15.6) para graficar T con respecto
a t en las siguientes condiciones iniciales:

(i) T(O) =
3000K CACO) = 0.533 moles/ piee
(ii) T(O) = 350 e..! (O) = 0.385
(iii) T(O) = 390 CACO) = 0.267
(iv) T(O) = 395 <. (O) = 0.252
(v) T(O) = 425 CACO) = 0.163
(vi) T(O) = 480.2 cA(O) = O
 Comportamiento de procesos transitorios 471

(c) Para los casos (i), (iii), (iv), y (vi) de la parte tb), que están cerca
del estado estable, utilice las ecuaciones 13.17 linealizadas y los
parámetros dados al final de la sección 13.2 para obtener T( r ).
¿ Cómo se compara la solución lineal con la solución a la
ecuación no lineal?
(d) Para los casos (ii) y (v) en la parte (b) compare la solución no
lineal con la solución de la ecuación linealizada. Utilice el estado
estable. con temperatura superior e inferior a la inicial para los
parámetros de la ecuación lineal izada.
13.4 Desarrolle métodos numéricos análogos a las ecuaciones 17.18 Y
17.19 para las ecuaciones diferenciales simultáneas.

dx
-;¡¡ = F(x, y) x(O) = Xo

d1f
d~ = G(x, y) y(O) = Yo

13.5 (a) Utilice el resultado del problema 13.4 o una rutina de biblioteca
que pueda conseguir de su centro de cálculo, para repetir
algunos de los cálculos de la parte (e) del problema 13.3, resol-
viendo nu méricarnente las ecuaciones linealizadas del reactor,
•
13.14 y 13.15. La comparación de las soluciones numéricas y las
soluciones analíticas a las ecuaciones linealizadas proporciona
una verificación del método numérico. Puede utilizarse una
calculadora, ya que una regla de cálculo no proporciona la
exactitud suficiente.
(b) Obtenga soluciones numéricas para los reactores adiabáticos no
lineales. Ecuaciones 13.1 y 13.2 donde u = O para las condicio-
nes iniciales dadas en los casos (ii) y (v) del problema 13.3, parte
(b). Verifique sus cálculos, comparando con la integración
numérica en la parte (b) del problema 13.3.
(c) Si está utilizando un programa de computadora, seleccione otras
condiciones iniciales que no satisfagan JeA (O) + T(O) = Je A, +
T, e integre las ecuaciones no lineales. Grafique los resultados
en el "plano de fase" corno se hizo en la figura 12.1 S, y compare
con los cálculos, utilizando el método de isocIinas (ecuación
12.61). Es posible, que su centro de cálculo disponga de un pro-
grama para la graficación en el plano de fase.
13.6 La reacción irreversible A -- productos, tiene una velocidad de
reacción rA- = koe-E/RT en _ Deduzca las ecuaciones transitorias
linealizadas para un reactor agitado tipo tanque de flujo continuo
adiabático que opera cerca del estado estable. Verifique su resulta-
 472 Sistemas no isotérmicos

do haciendo n = 1, comparando con las ecuaciones 13.14 y 13.15

para u == O.
13.7 Considere el problema de control en la secci6n 13.7. Demuestre que
para un control proporcional-más-integral:

w = K1y + K3itY(T) d-r

la t1uctuaci6n de temperatura y(t) irá a un valor de cero en el

estado estable para x,
y Y, constantes con K1 finita.
13.8 Para el problema de control en la secci6n 13.7, suponga que se
desea utilizar un control proporcional-más-derivada

Deduzca el efecto de la constante K2 para la forma de derivada.

13.9 Considere el problema de control de la secci6n 13.7. Demuestre que
una fluctuaci6n en la concentraci6n resultante de un cambio brusco
en la concentraci6n de alimentaci6n, no puede eliminarse mediante
un sistema de control de retroalimentaci6n que solamente ajuste la
velocidad de flujo del fluido de enfriamiento, basándose en la des-
viaci6n de temperatura.
13.10 Para el problema de control descrito en la secci6n 13.7, suponga
que u se mantiene constante, pero se ajusta q la velocidad de flujo a
través del reactor. Por consiguiente, tome el tiempo de residencia ()
como variable del control, y sea

T = V - ()s
q
Obtenga las ecuaciones linealizadas apropiadas, y determine el
efecto de un controlador proporcional en donde T es proporcional a
y.
13.11 La siguiente ecuaci6n describe un proceso controlado:
3 2
-d x + 6 -d x + 2
dx
11 - + 6x = u(t)
dt" dt dt
u(t) = -Kx

Demuestre que cuando K> 60 el sistema se inestabiliza. (Nota: la

ecuaci6n característica puede resolverse simplemente en términos de
K. Utilice la f6rmula para una ecuaci6n cúbica qu.e aparece en
Chemical Engineers' Handbook o en un libro de álgebra. Si 10
desea, suponga K»
1 para simplificar las relaciones.