Reactores con Flujos no Ideales

I.Q. Johsueé de Jesús Franco León

Dr.Alejandro González Álvarez

1
UNIVERSIDAD DE GUADALAJARA
Posgrado de Ingenierı́a Quı́mica

2
Distribución de Tiempos de
Residencia (RTD)

1 Función de distribución de tiempos de residencia E(t)

La distribución de tiempos de residencia (RTD) de un reactor quı́mi-

co en un fluido, es una distribución de probabilidad que describe la
cantidad de tiempo que un elemento diferencial de un fluido puede pa-
sar dentro de un reactor. Los ingenieros quı́micos usan el RTD para
caracterizar el mezclado y flujo dentro de un reactor y ası́ comparar
su conducta con sus modelos teóricos ideales correspondientes. Esto es
útil, no solo en la solución de problemas y detección de corto-circuitos
en reactores existentes, también es útil para estimar la eficiencias de
las reacciones involucradas y en el diseño de futuros reactores.

La función de distribución de tiempos de residencia es representa-

da por una distribución de edad de salida,E(t). La función E(t) está
normalizada y se define matemáticamente:

Z ∞
E(t) dt = 1
0

1.1 Determinación experimental de la RTD

La RTD es medida introduciendo un trazador no reactivo dentro del

sistema a la entrada, la concentración del trazador cambia acorde a una
función conocida y su respuesta se encuentra mediante la medición de
la concentración del trazador a la salida. El trazador electo no debe
modificar las propiedades fı́sicas del fluido (densidad y viscosidad) y
la introducción del trazador no debe modificar las condiciones hidro-
dinámicas. La expresión para el calculo de la RTD cuando se emplean
datos obtenidos con los experimento mencionados, es la siguiente:
3
 C(t)
E(t) = R ∞ =1
0 C(t) dt

1.2 Curva de distribución acumulativa

La curva F es otra función que se ha definido como respuesta nor-

malizada a una alimentación especı́fica. En ocasiones, la curva F se
emplea del mismo modo que la RTD para elaborar un modelo de los
reactores quı́micos.

1.3 RTD en CSTR perfectamente mezclado

Primero consideramos un tanque agitado perfectamente mezclado

con un fluido de densidad constante que fluye hacia adentro y hacia
afuera a caudales iguales. El RTD se puede analizar mediante una in-
yección de pulso o de reducción. Supongamos, en primera instancia,
4
que la corriente de entrada al recipiente contiene un trazador a una
concentración de C0, y que el sistema está en estado estable. La con-
centración del tanque y del efluente del marcador es, por lo tanto, igual
a C0. En el momento cero, la concentración del marcador en la alimen-
tación se reduce repentinamente a cero. Debido a que no se inyecta
más marcador, la concentración en el tanque comienza a disminuir. Lo
cual lleva a tener la siguiente expresión que define la concentración en
función del tiempo de residencia medio.

 
t
C(t) = C0exp −
tm

Si esta ecuación la sustituimos en la expresión de RTD para el calculo

experimental con un trazador, obtendrı́amos una expresión que descri-
be la RTD en un CSTR perfectamente mezclado.

 
C(t) C0exp − tmt exp(−(t/tm))
E(t) = R ∞ =R   =
∞ tm
0 C(t) dt
t
0 C0 exp − tm dt

5
Problema 1

Se inyectó una muestra del trazador hytane a 320 K en forma de

pulso a un reactor y se midió su concentración en el efluente en función
del tiempo, obteniéndose los datos que se muestran a continuación

t = (0,1,2,3,4,5,6,7,8,9,10,12,14) (min)

C = (0,1,5,8,10,8,6,4,3,2.2,1.5,0.6,0) (g/m3)

Las mediciones representan las concentraciones exactas en los tiem-

pos listados y no son valores promedio entre las diversas pruebas de
muestreo.

Solución

(a) Construya figuras mostrando C(t) y E(t) en función del tiempo.

6
Z ∞
C(t) dt = 50,65 (min/m3)
0

7
 C(t)
E(t) = R ∞ =1
0 C(t) dt

Esto indica que en el rango de 0 a 14 min sale el 100 % del material,

E(t) nos indica la edad de las partı́culas en el reactor.

(b) Determine tanto la fracción de material que sale del reactor y

que ha pasado entre 3 y 6 minutos dentro de él.

C(t)
E(t) = R 6 = 0,49
3 C(t) dt

Esto indica que en el rango de 3 a 6 min sale el 49 % del material.

8
 (c) Determine la fracción de material que sale del reactor y ha pasado
tres minutos o menos dentro del mismo.

C(t)
E(t) = R 3 = 0,19
0 C(t) dt

Por tanto en el rango de 0 a 3 min sale el 19 % del material.

A continuación presentamos la interface gráfica diseñada con MATLAB

para resolver el problema.

9
10
Problema 2

Este es un proceso critico y no puede durar el material más de 10

min en el reactor, indique que porcentaje de material ha salido en este
tiempo.

Los datos que se muestran a continuación son la respuesta de con-

centración a una entrada de pulso a un recipiente. Grafique función
F(t).

t = (0, 0,1, 0,2, 1, 2, 5, 10, 30) (min)

C = (0,50, 0,40, 0,34, 0,30, 0,25, 0,14, 0,04, 0,002) (mol/L)

Solución

La función F (t), esta dada por la siguiente expresión:

Z t
F (t) = E(t)dt
0

Representa la fracción de efluente que ha estado en el reactor un

tiempo menor a t, con ayuda de MATLAB, graficamos la función y
obtenemos lo siguiente:

11
A los 10 min del proceso, sale el 80 % del material.

A continuación se resuelve con la interface gráfica

12
Problema 3

Distintos experimentos llevados a cabo en un reactor continuo de

tanque agitado hacen sospechar que el comportamiento del reactor no
es ideal. La respuesta ante un ensayo de trazador en impulso viene
dada por la ecuación C(t) = 5texp(−2,5t) (con t en min). Parece ser
que el reactor se puede modelar considerando que los elementos de
fluido de distintas edades no se mezclan entre sı́, aparentando como si
pequeños reactores discontinuos estuvieran operando en el interior del
reactor en continuo. Cómo dato adicional, se sabe qué el trazador duro
aproximadamente 4 min en el reactor.

a) ¿Qué fracción de efluente salio después del primer minuto tras la

inyección del trazador?
Solución

Primero calculamos, la integral de la función de concentración vs

tiempo

13
 Z 4
C(t) dt = 0,7954
0

Con esto calculamos la edad del efluente, gráficamente E(t) vs t

C(t)
E(t) = R 1 = 0,71
0 C(t) dt

Y comprobamos el resultado con la curva F(t), la cual se presenta a

continuación.

14
 Lo anterior se puede resolver de la siguiente manera con la interface
gráfica.

15
2 Momentos de la RTD

Es muy común para comparar RTD emplear momentos, en lugar de

intentar emplear toda la distribución . Para tal finalidad, se emplean
normalmente tres momentos.

El primero es el tiempo medio de residencia. El segundo momento,

que se emplea de manera común, se toma en torno a la media y se
llama varianza o cuadrado de la desviación estándar.

2.1 Tiempo de residencia

El primer momento proporciona el tiempo promedio que las molécu-

las del efluente pasan dentro del reactor. Cuya ecuación es la siguiente

Z ∞
tm = tE(t) dt
0

El uso más importante de la teorı́a del tiempo de residencia es su

aplicación a equipos que ya están construidos y en funcionamiento.

Si no hay zonas muertas o estancamientos dentro del reactor, en-

tonces tm será igual al tiempo de residencia teórico τ , este se define
como la relación entre el volumen del reactor y el flujo volumétrico de
entrada del fluido:

V
τ=
ν

2.1.1 Volumen muerto y volumen activo en reactores

CSTR

Varios factores influyen en el tener zonas muertas en los reactores

continuos, por ejemplo:
16
 -Diseño geométrico del reactor.
-Mezclado y agitación.
-Viscosidad del fluido.
-Distribución de relleno en lechos empacados

El volumen muerto esta dado por la siguiente expresión:

Vm = V − (tmυ)

Y el volumen activo por :

Va = V − Vm

2.2 Varianza

Ecuaciones de momentos de grado superior pueden proveer informa-

ción sobre la conducta del la función E(t). Por ejemplo, el segundo
momento central indica la varianza (σ 2). El grado de dispersión com-
parado con la media.

Z ∞
2
σ = (t − tm)2E(t) dt
0

17
 Problema 4
Calcule el tiempo medio de residencia y la varianza para el reactor
caracterizado en el ejemplo 1 por la RTD obtenida con una alimentación
por pulso a 320 K.

Solución

Para calcular el tiempo de residencia partimos de la siguiente expre-

sión

Z ∞
tm = tE(t) dt
0

Esto significa que el tiempo de residencia es el área bajo la curva,

para la gráfica de t vs E(t) que es la gráfica c) de la figura.

Z 14
tm = tE(t) dt = 5,123min
0

Por otra parte, para el cálculo de la varianza, es necesario que cal-

culemos primero la desviación estándar. Con la siguiente expresión:

Z ∞
2
σ = (t − tm)2E(t) dt
0

El valor de la desviación estándar es el área bajo la curva, para la

gráfica de E(t)(t − tm)2 vs E(t) que es la gráfica d) de la figura.

18
 Z 14
2
σ = (t − tm)2E(t) dt = 5,951min2
0

por tanto la varianza es:

σ = 2,44min

A continuación se muestra la imagen donde el problema se resuelto

con la interface gráfica.

19
20
Problema 5

Al estudiar las caracteristicas de mezclado de los reactores quimicos,

se inyectan un pulso de un indicador que no reacciona en el reactor, el
indicador se inyecta a tiempo igual a cero. La concentración de material
en el efluente del reactor es medido en función del tiempo C(t). En este
estudio en el que se estudieron las caracteristicas de mezclado en un
reactor de flujo continuo, se obtuvieron los siguientesdatos

t = (0, 1, 2, 3, 4, 5, 6, 7, 8) (s)

C = (0, 2, 4, 7, 6, 5, 2, 1, 0) (mg/L)

Calcular la función RTD, el tiempo de residencia y la varianza de la

función de RTD

Solución

Calculando el tiempo medio de residencia:

Z 8
tm = tE(t) dt = 3,66min
0

La desviación estándar es la siguiente:

Z 8
2
σ = (t − tm)2E(t) dt = 2,22
0

21
 A continuación se presentan las gráficas que demuestran el compor-
tamiento de la RTD, tm y σ 2

Resolviendo el problema con la interface gráfica.

22
Problema 6

En un reactor de volumen V, el cual es igual a 15 L se desea llevar

a cabo la reacción A → B en fase lı́quida a temperatura constante. El
alimento se introduce con un caudal volumétrico v de 4.0 L/min.Para
caracterizar previamente el comportamiento del reactor, y sospechando
que éste se comporta como un reactor tipo tanque con volumen muerto,
se realiza un ensayo con trazador en impulso utilizando el v indicado
anteriormente y se mide la variación de la concentración de trazador
con el tiempo a la salida del reactor. Los resultados obtenidos son los
que se muestran a continuación:

t = (0,0.8,2.0,4.0,6.8,8.0,12.8,16.8,24.0) min

C = (15.37,12.15,9.43,5.24,2.08,1.62,0.36,0.12,0.001) mg/L

Calcular el tiempo medio de residencia y los volúmenes muerto y

activo del reactor.

Solución

Primero con los datos experimentales de la RTD, calculamos el tiem-

po medio de residencia y posteriormente, calculamos el volumen muerto
partimos de la siguiente expresión

tm = 3,3022

Vm = V − (tmυ) = 15,0 − (3,3 ∗ 4,0) = 1,79L

Y el volumen activo seria de :

23
 Va = V − Vm = 15,00 − 1,79 = 13,209L

Comprobando con la interface gráfica, tenemos lo siguiente

24
2. MODELOS DE REACTORES NO
IDEALES
No todos los reactores de tanque están perfectamente mezclados ni
todos los reactores tubulares presentan comportamiento de flujo tapón.
En tales casos, es necesario aplicar métodos para considerar desviacio-
nes del comportamiento ideal.

Los modelos son útiles para representar el flujo en recipientes reales,

para el escalamiento y para identificar el flujo defectuoso. Existen dife-
rentes tipos de modelos, lo que depende de si el flujo es cercano a tanque
agitado, flujo pistón o intermedio entre estos. La siguiente ecuación no
ayudara a entender los paso necesarios par a hacer una adecuada pre-
dicción en reactores no ideales.

Datos de RTD + Cinética + Modelo = Predicción

2.1 Modelo de Dispersión

El modelo de dispersión considera un fluido en flujo pistón que tiene

sobrepuesto un cierto grado de retromezclado, cuya magnitud es inde-
pendiente de la posición dentro del recipiente, esto quiere implica que
no existen corto circuitos o zonas muertas.

Este fenómeno se rige por la ley de fick, en su forma normalizada se

representa de la siguiente manera:

∂ 2C ∂C
 
∂C D
= −
∂Θ uL ∂z 2 ∂z

Donde D/uL es el módulo de dispersión del recipiente.

25
D/uL → 0 la dispersión es insignificante, por lo tanto hay flujo pistón

D/uL → ∞ gran dispersión, por lo tanto hay tanque agitado.

Se puede calcular este módulo a partir de diferentes curvas de RTD,

por ejemplo las curvas de RTD normalizadas en las que se tiene un
flujo piston para D/uL menores a 0.01 y si es mayor se tiene flujo de
tanque agitado.

P er se conoce como número de Peclet del reactor y emplea la lon-

gitud del reactor, L, como longitud caracterı́stica, de modo que P er =
uL/D
26
 El número de Peclet se puede encontrar experimentalmente usando
una inyección de pulso en la entrada del tubo y encontrando la función
RTD como se describe en la primera sección. El número de Peclet
depende de la varianza de la función RTD normalizada y se calcula de
la siguiente manera:

2
   2
σ D D
σΘ2 = 2 = 2 −2 [1 − e−uL/D ]
tm uL uL

Para pequeños D/uL se puede utilizar la siguiente ecuación, ya que

el segundo termino de la expresión anterior, es insignificante.

2
 
σ D
σΘ2 = 2 = 2
tm uL

2.1.1 Conversión en Modelo de Dispersión

Cundo la reacción es de primer orden, la conversión se calcula con

la siguiente expresión

4qexp(P er /2)
X =1−
(1 + q)2exp(P er q/2) − (1 − q)2exp(−P er q/2)

p
q= 1 + (4Da/P er )

A Da se le conoce como el número de Damkohler para la convección

27
 Da = ktm

Si la dispersión fuera cero, implicarı́a que tendrı́amos el compor-

tamiento de un PFR ideal, y el modelo anterior se convertirı́a en el
siguiente.

X = 1 − exp(−tmk)

Cuando se tiene una reacción de segundo orden, se puede linealizar

la ley de velocidad de la siguiente manera:

CA0
−rA = kCA2 = k CA = k 0 CA
2

Da = k 0tm

Para un PFR ideal se segundo orden, la expresión seria la siguiente

ktmCA0
X=
1 + ktmCA0

28
Problema 7
Se realizaron estudios de RTD en un reactor tubular, se inyecto
rápidamente a la entrada del reactor una solución de NaCl y se tomaron
muestras a la salida y se obtuvieron los siguientes resultados

t = (10,22.5,27.5,32.5,37.5,42.5,47.5,52.5,62.5) h

C =(0,60,210,170,75,35,10,5,0) mol/m3

Calcule el modulo de dispersión

Solución

Para calcular el modulo de dispersión es necesario calcular el tiempo

medio de residencia y la desviación estándar

Z 57,5
tm = tE(t) dt = 30,7449
10

Z 57,5
2
σ = (t − tm)2E(t) dt = 39,9808
00

Con estos datos procedemos a calcular el modulo de dispersión

   2
D D
0,04058 = 2 −2 [1 − e−uL/D ]
uL uL

El modulo de dispersión para este experimento es el siguiente

29
 D
= 0,02114
uL
Utilizando la interface gráfica para resolver el problema, tenemos lo
siguiente

30
Problema 8
Se caracterizo un reactor no ideal con un experimento de RTD, en
el cual se determino un tiempo de residencia medio de 7.2 min y un
modulo de dispersión igual a 0.0988. La concentración inicial del ex-
perimeto fue de CA0 = 0.3 mol/min y una constante de velocidad k =
0.2.

Calcular la conversión para el modelo de dispersión y el modelo de

PFR ideal, hacerlo para reacciones de primer y segundo orden.

Solución
a) Modelo de dispersión de primer orden

Calculamos el numero de Pe, haciendo la inversa del modulo de

dispersión P e = (D/uL)−1 = 0,0988−1 = 10,1214, ahora sustituimos
en la ecuación

Da = 1,44

p
q = 1 + (4(1,44)/10,1214) = 1,2526

4qexp(P er /2)
X =1−
(1 + q)2exp(P er q/2) − (1 − q)2exp(−P er q/2)

X = 0,7250

b) Modelo de PFR ideal de primer orden

31
 X = 1 − exp(−1,44)

X = 0,7640

Utilizando la interface gráfica, resolvemos los do primero incisos y

comprobamos las respuestas.

c) Modelo de dispersión para segundo orden

En este caso el numero Da cambia porque este ya esta en función de

la concentración inicial, esto después de linealizar la ley de velocidad.
Por otro lado el modulo de dispersión sera el mismo y por tanto el
numero de Pe también.

Da = 0,432

32
 p
q= 1 + (4(0,432)/10,1214) = 1,0820

4qexp(P er /2)
X =1−
(1 + q)2exp(P er q/2) − (1 − q)2exp(−P er q/2)

X = 0,3406

d) Modelo de PFR ideal para segundo orden

ktmCA0
X=
1 + ktmCA0

0,432
X=
1 + 0,432

X = 0,3016

Utilizando la interface gráfica, tenemos lo siguiente

33
34
2.2 Modelo de tanques en serie

El modelo de tanques en serie es de un solo parámetro. Este mo-

delo analiza la RTD para determinar el número de tanques ideales, n,
colocados en serie que darán aproximadamente la misma RTD que el
reactor no ideal.

Para hacer la relación de la RTD a un sistema de tanques en serie,

se parte de la siguiente expresión:

tn−1
E(t) = n e−t/τi
(n − 1)!τi

Pero para que esta función sea más manejable en la mayorı́a de la

literatura se utiliza de forma normalizada y las expresiones son las
siguientes.

35
 Las ecuaciones anteriores se pueden ver gráficamente de la siguiente
manera.

36
2.2.1 Conversión en Modelo de Tanques es serie

Para una reacción de primer orden, cuando modelamos un tanque

CSTR ideal se representa bajo la siguiente expresión

1
X=
(1 + tmk)

De la expresión anterior se parte para modelar el modelo de tanque

en serie para una reacción de primer orden, el cual es el siguiente.

1
X =1−
(1 + τik)n

NOTA: El modelo de tanques en serie acepta un numero y fracción

para los tanques cuando la reacción es de orden uno, pero esto no
aplica cuando la reacción es de un orden mayor a uno, para estos caso
el modelo de tanques en serie únicamente acepta numero de tanques
enteros.

Cuando la reacción es de segundo orden, los siguientes modelos se

aproximan a la solución, tanto para el modelo de CSTR ideal como el
modelo de tanques en serie.

√
−1 + 1 + 4tmkCA0
X =1−
2tmkCA0

√ 1/n
−1 + (1 + 2(1 + 1 + 4tmkCA0)
X =1−
2tmkCA0

37
Problema 9
Calcular los tanques en serie que darán aproximadamente la misma
RTD que el reactor no ideal en el cual se obtuvieron los siguientes datos
en el experimento con el trazador.

t = (0,1,2,3,4,5,6,7,8,9,10,12,14) min

C = (0,1,5,8,10,8,6,4,3,2.2,1.5,0.6,0) g/m3

Solución

Iniciamos por calcular el tiempo medio de residencia y la varianza,

esto ya fue calculado en el problema 4 para estos datos experimentales.

tm = 5,123min

σ 2 = 5,951min2

Calculamos con estos datos la varianza normalizada.

σ2
σΘ2 = 2 = 0,2267
tm

Finalmente calculamos el numero de tanques en serie

1
N= = 4,41tanques
σΘ2
38
 Utilizando la interface gráfica, observamos que se comprueban los
resultados.

39
Problema 10
Se tiene instalado un tanque CSTR en una planta farmacéutica,
por cuestiones de un nuevo producto es necesario eliminar dispersión
y tener un perfil de flujo pistón.

A través de un pulso de rastreador en el CSTR se obtienen los si-

guientes datos.

t = (0,5,10,15,20,30,40,50,70,100,125,150,175,200) min

C = (112,96,82.2,70.6,61,45,34,26,15.7,7.7,5.1,2.5,1.7,0.9) mol/m3

partir de la funciones RTD, calcule los tanques necesarios para tener

un flujo pistón y acercarlos a un PFR ideal.

Solución

Para tener un flujo tapón (Reactor PFR) en una serie de reactores

de tanque agitados es necesario que la dispersión del flujo dentro del
reactor sea mı́nima, esto significa tener un D/uL < 0,01. Donde el
numero de P e esta dado por la siguiente ecuación:

σΘ2
D/uL =
2
Las funciones RTD normalizadas para tanques en serie nos dan la
respuesta, cuando utilizamos 51 tanques en serie tenemos lo siguiente.

1 1
σΘ2 = = = 0,01961
N 51

0,01961
D/uL = = 0,00903
2
40
 Al tener un D/uL = 0,00903 < 0,01 podemos decir que tenemos un
flujo tapón, si graficamos es la funcion RTD normalizada deberiamos
tener una campana Gaussina mas estrecha y alargada a medida que
los reactores en serie aumentan y el D/uL disminuye. A continuacion
se muestra la grafica en la cual se describe claramente la prediccion.

41
Problema 11
Consideremos la reacción de primer orden A → B que se lleva a cabo
en un reactor CSTR no ideal en fase lı́quida. La constante cinética es
0.25 min−1. En el reactor se realiza una experimento con trazador en
impulso y se obtienen los siguientes resultados:

t = (0,1,2,3,4,5,6,7,8,9,10,12,14) min

C = (0,1,5,8,10,8,6,4,3,2.2,1.5,0.6,0) g/m3

Calcular la conversión que se obtendrá en el sistema utilizando: a)

modelo de CSTR ideal; b) modelo de tanques en serie. c) ¿ cuales son
las conversiones si la reacción es de segundo orden, la concentración
inicial CA0 = 2mol/m3

Solución

a) Conversión en modelo de CSTR ideal primer orden

Primero es necesario calcular el tiempo de residencia medio (tm),

para esto utilizamos la interface gráfica:

tm = 5,12734 min

1 1
X= =
(1 + tmk) 1 + 1,281835

X = 0,5617

b) Conversión en modelo de tanques en serie primer

orden
42
 Los tanques en serie que representan esta RTD fueron calculados
en el problema 9 y fueron N = 4.41 Tanques. Con esto partimos para
calcular la conversión.

1
X =1−
(1 + τik)n

Donde τi = 5,12734/4,41 = 1,1626

X = 0,6754

Aplicando la interface gráfica, comprobamos el resultado obtenido

para los dos primeros casos.

43
2.3 Modelo de Segregación

En un CSTR ”perfectamente mezclado”; se asume que el fluido de

entrada se distribuye de inmediato y equitativamente en toda la mezcla
de reacción. Se asume que dicho mezclado tiene lugar inclusive a mi-
croescala, es decir, todos los elementos de distintas edades se mezclan
a la perfección para formar un fluido en una micro mezcla. Cuando los
elementos del fluido de distintas edades no se mezclan en absoluto y
éstos permanecen segregados unos de otros, se dice que el fluido está
totalmente segregado. Los extremos del micromezclado y la segregación
totales son los lı́mites del micro mezclado de una mezcla de reacción.

2.3.1 Conversión en Modelo de Segregación

Para determinar la conversión media en la corriente del efluente, es

necesario sacar el promedio de las conversiones de los diversos glóbulos
en la corriente de salida:

dX̄
= X(t)E(t)
dt
Haciendo la sumatoria de todos los glóbulos, la conversión media es

Z ∞
X̄ = = X(t)E(t)dt
0

Para una reacción de primer orden, la X(t) puede calcularse con la

siguiente manera

X(t) = 1 − e−ktm

NOTA: Algo muy importante y que ya fue estudiado, es el echo

que para reacciones de primer orden, la conversión tanto en modelos
44
de tanques en serie, modelo de segregación y el modelo de mezclado
máximo es el la misma, pero esto solo aplica para reacciones de primer
orden.

45
Problema 12
Calcular la conversión para una reacción de primer orden, con los
datos experimentales del problema 11, comprobar que los resultados
de la conversión con los mismos para el modelo de tanques en serie y
el modelo de segregación.

Solución

Para la RTD del problema anterior la constante de velocidad fue de

k = 0,25min−1 y el tiempo medio de residencia tm = 5,12734min.
Sustituyendo en la integral, tenemos:

Z 14
X̄ = = (1 − e−1,218)E(t)dt
0

Lo anterior lo integramos en el intervalo de los datos obtenido en el

experimento los cuales fueron:

t = (0,1,2,3,4,5,6,7,8,9,10,12,14) min

C = (0,1,5,8,10,8,6,4,3,2.2,1.5,0.6,0) g/m3

Con ayuda de un programa matemático como MATLAB calculamos

la integral definida y la conversión obtenida es la siguiente:

X̄ = 0,6746

46
 Efectivamente para los dos modelo, tanto el modelo de tanques en
serie y el modelo de segregación para una reacción de primer orden, la
conversión es la misma.

Aplicando esto a la inteface gráfica, comprobamos los resultados.

47
2.4 Modelo de Mezclado Máximo

En un reactor con fluido segregado no ocurre mezclado entre partı́cu-

las de fluido hasta que dicho fluido sale del reactor. Por supuesto, la
salida del reactor es el punto más tardı́o posible, en el cual puede ocurrir
el mezclado.Esto quiere decir que en el mezclado máximo, el mezcla-
do ocurre en el punto más tardı́o posible, caso contrario al modelo de
segregación en el que el mezclado ocurre en el punto más temprano
posible.

Al desplazamos por el reactor de izquierda a derecha, λ disminuye

y llega a cero en la salida. En el extremo izquierdo del reactor, λ se
aproxima al infinito o el tiempo de residencia máximo, si es distinto de
infinito.

2.4.1 Conversión en Modelo de Mezclado Máximo

Diferentes estudios demuestran que para una reacción de primer

orden, la conversión sera la misma tanto en el modelo de mezclado
máximo y el modelo de segregación y por tanto en el modelo de tanques
en serie.

Cuando tenemos una reacción de segundo orden, la ecuación que

describe el modelo de mezclao máximo es la siguiente.

48
 dX θN −1N N /(N − 1)!
= PN −1 X − kCA0τ (1 − X)
dθ i i
i=0 θ N /i!

Si queremos obtener la solución, la ecuación se integra hacia atrás

numéricamente comenzando con un valor muy grande de θ y terminan-
do con la conversión final en θ = 0. Con una RTD dada y órdenes de
reacción mayores de uno, el modelo de mezclado máximo proporciona
el lı́mite inferior de la conversión.

Se puden utilizar paquetes de solución matemática para resolver es-

tos problemas pero es necesario hacer varias modificaciones, otra forma
es resolverlos con métodos numéricos, como el método de integración
de euler.

49
Problema 13
Calcular la conversión para el modelo de mezclado máximo en una
reacción de segundo orden,con un numero de convección ktmCA0 = 10.
Considere que los tanques en serie que describen el modelo son 3.

Solución

La conversión se calcula con la siguiente ecuación y que aplica para

N tanques

dX θN −1N N /(N − 1)!

= PN −1 X − kCA0τ (1 − X)
dθ i i
i=0 θ N /i!

Aplicando el modelo para tres tanques y con la ayuda de MATLAB

calculamos la conversión, esto se debe hacer integrando hacia atrás.

50
La conversión con este modelo seria la siguiente:

XAE = 0,8341

51
 PROBLEMAS EXTRAS

Problema 14
º 
La reacción de dimerización de 2 orden 2A B (rA = -kCA 2) tiene
·
lugar en fase lı́quida, siendo k = 0.01 dm3 /(mol min) a la temperatura
de reacción. El alimento es A puro con una concentración de 8 mol/L.
El volumen teórico del reactor es de 1000 L y el caudal de alimentación
para la dimerización es de 25 L/min.

Para poder conocer el comportamiento del reactor, el cual no se

comporta idealmente, se ha desarrollado un ensayo con trazador en
impulso con un caudal de 25 L/min, cuyos resultados se muestran en
la siguiente tabla:

Se desea conocer los lı́mites entre los que puede variar la conversión
según el grado de micromezcla.

Solución

Primero calculamos la conversión para el modelo de segregación com-

pleta con ayuda de la interface gráfica

52
 La conversión para el modelo de segregación completa es:

X = 0,6058

A hora procedemos a calcular el la conversión para el modelo de

mezclado , con ayuda de la interface gráfica primero calculamos el
tiempo medio de residencia el cual es tm = 37,2514 min.Resolviendo
con la interface gráfica obtenemos la siguiente gráfica. La cual es la
solución de la EDO que describe la conversión en el modelo de mezclado
máximo.

53
 La conversión para el modelo de mezclado máximo es:

X = 0,5647

Por tanto los lı́mites entre los que puede variar la conversión según
el grado de micromezcla son X : 0,5647 − 0,6058

54
Problema 15
La respuesta a una entrada de pulso a un reactor de lecho empaque-
tado se da en la tabla a continuación. Calcule el número de Peclet para
el sistema y el número de tanques en serie que se usarı́an para modelar
este reactor.

Solución

Resolvemos el problema aplicando la interface gráfica:

55
 Los tanques en serie que describen la RTD son N = 12,9998 = 13

Para el calculo del numero de Pe, la interface grafica nos da el valor

de D/uL pero sabemos que el numero de Pe es el reciproco del modulo
de dispersión.

2
   2
σ D D
σΘ2 = 2 = 2 −2 [1 − e−uL/D ]
tm uL uL

D/uL = 0,04039 P e = 24,758

56
Problema 16
El flujo que se alimenta a un reactor es de 10 dm3/min. Una prueba
por pulso dio las siguientes medidas de concentración a la salida:

A) Obtenga las gráficas de E(t) y F(t)

B) ¿Qué fracción del material pasa de 2 a 4 minutos en el reactor?

C) ¿Qué fracción del material pasa menos de 6 minutos en el reactor?

D) ¿Qué fracción del material pasa menos de 3 minutos en el reactor?

E) ¿Qué valor tienen el tiempo medio de residencia, tm, y la varianza,

σ2?

F) ¿ Qué conversión se lograrı́a en un PFR ideal y un CSTR ideal

para una reacción de segundo orden con kCA0 = 0,1min−1yCAO =
1mol/dm3?

57
Solución

A) Gráficas de E(t) y F(t)

58
 B) Fracción del material pasa de 2 a 4 minutos en el
reactor

59
 Z 4
E(t)dt = 0,162 o el 16,2
2

C) Fracción del material pasa más de 6 minutos en el

reactor

Z 6
0E(t)dt = 0,588 o el 58,8
6

D) Fracción del material pasa menos de 3 minutos en

el reactor

60
 Z 3
E(t)dt = 0,188 o el 18,8
0

E) Valor del tiempo medio de residencia, tm, y la va-

rianza, σ 2

Utilizamos la interface gráfica para encontrar estos dos valores.

El tm = 9,8633 min y σ 2 = 74,55 min2

61
 F) Conversión para PFR y CSTR ideal

Igualmente utilizamos la interface gráfica el valor de ktmCA0 =

0,98633

La conversión para el CSTR ideal es de X = 0,3796 y la conversión

para el PFR ideal es de X = 0,4965
62
Problema 17
Ingenieros quı́micos, desean realizar una gráfica que muestre la con-
versión en función del numero de tanques en serie y número de Damkoh-
ler para la convección, esto en el modelo de mezclado máximo para una
reacción de segundo orden. El rango del los parametros de estudio son
de hasta 10 tanques en serie y un Da = [1,2,3,4,5,6,7,8,9,10,15,20,25,30]

Solución

Con ayuda de MATLAB e integrando la ecuación diferencial hacia

atras, calculamos los valores con la siguiente ecuación.

dX θN −1N N /(N − 1)!

= PN −1 X − kCA0τ (1 − X)
dθ i i
i=0 θ N /i!

A continuación se presentan los resultados y la gráfica obtenida.

63
64
65
 El perfil muestra que a medida que aumentan los tanques en serie y
el numero de convección, la conversión sera mayor.

66
Problema 18
Los datos que se muestran a continuación son la respuesta de con-
centración a una entrada de pulso a un recipiente. Calcule la función
RTD y compare sus valores con un CSTR que tenga el mismo tiempo
medio de residencia.

Solución

Iniciamos por calcular el tm, el cual puede ser calculado con la in-
terface gráfica.

Con el tm = 4,6927 podemos calcular la RTD del CSTR perfecta-

mente mezclado con la siguiente expresión:

67
 exp(−(t/tm))
E(t) =
tm

Posteriormente calculamos los valores de la RTD para los datos del

experimento del pulso.

C(t)
E(t) = R 30
0 C(t) dt

Y donde:

Z 30
C(t) dt = 2,068
0

Si graficamos las dos funciones obtenemos los siguientes perfiles.

68
Existe una similitud entre las dos funciones pero no son perfectas.

69
REFERENCIAS
FFogler, H.S., .Elements of Chemical Reaction Engineering”, Prentice
Hall, N.J., 2008.4ta. Ed.

Levenspiel, O., Ïngenierı́a de las Reacciones Quı́micas”, Reverte, Bs.

As., 1981, 3rd. Ed

R.E. Hayes, J.P. Mmbaga., Ïntroduction to Chemical Reactor Analy-

sis”,2013, 2da Ed. Escrito por R.E. Hayes, J.P. Mmbaga

70
APÉNDICE

Código MATLAB para funciones RTD

clc
f u n c i o n e s = menu ( ’ Funciones RTD’ , ’ Datos Vectores ’ . . .
, ’ Funcion ’ )
switch funciones
case 1
t = i n p u t ( ’ Datos de tiempo= ’ ) ;
C = i n p u t ( ’ Datos de c o n c e n t r a c i o n= ’ ) ;
case 2
xn = 0 ;
x l = i n p u t ( ’ Tiempo f i n a l a l que s a l e e l t r a z a d o r \n ’ ) ;
t = [ xn : . 1 : x l ] ;
* *
% ejemplo de f u n c i o n ( 5 . t ) . ( exp ( −2.5. t ) ) ;*
C = input ( ’ I n t r o d u c i r funcion = ’ ) ;
end
f1 = figure ;
p l o t ( t , C, ’ k ’ , ’ LineWidth ’ , 1 . 5 )
Area = t r a p z ( t ,C ) ;
p = a r e a ( t , C, ’ LineWidth ’ , 1 . 5 , ’ FaceAlpha ’ , 0 . 4 ) ;
s e t ( p , ’ FaceColor ’ , ’ k ’ ) ;
g r i d on
$ $
x l a b e l ( ’ t , ( min ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , . . .
’ fontsize ’ ,16)
$ $
y l a b e l ( ’ C , ( g/mˆ3) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ . . .
, ’ fontsize ’ ,16)
E = C. / Area ;
f = trapz ( t ,E) ;
F = cumtrapz ( t , E ) ;
x0 = i n p u t ( ’ Tiempo i n i c i a l en e l que s e . . .
71
 d e s e a c o n o c e r l a s a l i d a d e l m a t e r i a l=\n ’ ) ;
x f = i n p u t ( ’ Tiempo f i n a l en e l que s e d e s e a . . .
c o n o c e r l a s a l i d a d e l m a t e r i a l=\n ’ ) ;
ind = ( t > x0 ) & ( t < x f ) ;
xw = t ( ind ) ;
yw = E( ind ) ;
ya = i n t e r p 1 ( t , E, x0 ) ;
yb = i n t e r p 1 ( t , E, x f ) ;
xw = [ x0 , xw , x f ] ;
yw = [ ya , yw , yb ] ;
A = t r a p z (xw , yw ) ;
*
A1 = A 1 0 0 ;
f2 = figure ;
p l o t ( t , E, ’ k ’ , ’ LineWidth ’ , 1 . 5 )
$ $
x l a b e l ( ’ t , ( min ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , . . .
’ fontsize ’ ,16)
$ $
y l a b e l ( ’ E( t ) , ( 1 / min ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , . . .
’ fontsize ’ ,16)
g r i d on
hold on
h=a r e a (xw , yw , ’ FaceAlpha ’ , 0 . 4 ) ;
s e t ( h , ’ FaceColor ’ , ’m’ ) ;
f p r i n t f ( ’ P o r c e n t a j e de m a t e r i a l que pasa en . . .
e l r e a c t o r en e s e lampso de tiempo : %10.1 f \n\n ’ , A1)
% Grafica F
f 3=f i g u r e ;
x =[ xf , x f ] ;
y =[0 ,A ] ;
x1 =[0 , x f ] ;
y1 =[A,A ] ;
p l o t ( t , F , ’ k−o ’ , x , y , ’m: ’ , x1 , y1 , ’m: ’ , ’ LineWidth ’ , 1 . 5 )
x l a b e l ( ’ t , ( min ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , ’ f o n t s i z e ’ , 1 6 )
y l a b e l ( ’F( t ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , ’ f o n t s i z e ’ , 1 6 )

72
g r i d on

Código MATLAB para módulos de RTD

clc
f u n c i o n e s = menu ( ’ Modulos DTR’ , . . .
’ Datos Vectores ’ , ’ Funcion ’ )
switch funciones
case 1
t = i n p u t ( ’ Datos de tiempo= ’ ) ;
C = i n p u t ( ’ Datos de c o n c e n t r a c i o n= ’ ) ;
case 2
xn = 0 ;
x l = i n p u t ( ’ Tiempo f i n a l a l que s a l e e l t r a z a d o r \n ’ ) ;
t = [ xn : . 1 : x l ] ;
C = input ( ’ I n t r o d u c i r funcion = ’ ) ;
end
subplot (2 ,2 ,1);
p l o t ( t , C, ’ k ’ , ’ LineWidth ’ , 1 . 5 )
Area = t r a p z ( t ,C ) ;
p = patch ( t , C, ’m’ , ’ LineWidth ’ , 1 . 5 , ’ FaceAlpha ’ , 0 . 3 ) ;
g r i d on
$ $
x l a b e l ( ’ t , ( min ) , a ) ’ , ’ i n t e r p r e t e r ’ , . . .
’ latex ’ , ’ fontsize ’ ,16)
$$
ylabel ( ’ C ’ , ’ interpreter ’ , ’ latex ’ , ’ fontsize ’ ,16)
E = C. / Area ;
x0 = i n p u t ( ’ Tiempo i n i c i a l de i n t e g r a c i o n \n ’ ) ;
x f = i n p u t ( ’ Tiempo f i n a l de i n t e g r a c i o n \n ’ ) ;
ind = ( t > x0 ) & ( t < x f ) ;
xw = t ( ind ) ;
yw = E( ind ) ;
ya = i n t e r p 1 ( t , E, x0 ) ;
yb = i n t e r p 1 ( t , E, x f ) ;
73
xw = [ x0 , xw , x f ] ;
yw = [ ya , yw , yb ] ;
A = t r a p z (xw , yw ) ;
*
A = A 100;
subplot (2 ,2 ,2);
p l o t ( t , E, ’ k ’ , ’ LineWidth ’ , 1 . 5 )
$ $
x l a b e l ( ’ t , ( min ) , b ) ’ , ’ i n t e r p r e t e r ’ , . . .
’ latex ’ , ’ fontsize ’ ,16)
$ $
y l a b e l ( ’ E( t ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , . . .
’ fontsize ’ ,16)
g r i d on
hold on
a r e a (xw , yw , ’ FaceAlpha ’ , 0 . 3 )
f p r i n t f ( ’ P o r c e n t a j e de m a t e r i a l que pasa en . . .
e l r e a c t o r en e s e lampso de tiempo : %10.1 f \n\n ’ ,A)
%TIEMPO DE RESIDENCIA
tE = t . E ;*
subplot (2 ,2 ,3);
p l o t ( t , tE , ’ k ’ , ’ LineWidth ’ , 1 . 5 )
$ $
x l a b e l ( ’ t , ( min ) , c ) ’ , ’ i n t e r p r e t e r ’ , . . .
’ latex ’ , ’ fontsize ’ ,16)
$ $
y l a b e l ( ’ tE ( t ) ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , . .
’ fontsize ’ ,16)
g r i d on
hold on
tm = t r a p z ( t , tE ) ;
p = patch ( t , tE , ’ y ’ , ’ LineWidth ’ , 1 . 5 , ’ FaceAlpha ’ , 0 . 3 ) ;
f p r i n t f ( ’ Tiempo de R e s i d e n c i a Medio : %10.3 f \n\n ’ , tm)
%VARIANZA
tmE=(t−tm ) . ˆ 2 ;
v =(tmE . E ) ; *
subplot (2 ,2 ,4);
p l o t ( t , v , ’ k ’ , ’ LineWidth ’ , 1 . 5 )

74
 $ $
x l a b e l ( ’ t , ( min ) , d ) ’ , ’ i n t e r p r e t e r ’ , . . .
’ latex ’ , ’ fontsize ’ ,16)
$ $
y l a b e l ( ’ E( t ) ( t−tm)ˆ2 ’ , ’ i n t e r p r e t e r ’ , . . .
’ latex ’ , ’ fontsize ’ ,16)
g r i d on
hold on
vari = trapz ( t , v ) ;
a r e a ( t , v , ’ FaceAlpha ’ , 0 . 3 ) ;
f p r i n t f ( ’ D e s v i a c i o n Estandar : %10.3 f \n\n ’ , v a r i )
S = sqrt ( vari );
f p r i n t f ( ’ Varianza : %10.3 f \n\n ’ , S )

Código MATLAB para volumen muerto y activo de un

reactor CSTR

% C a l c u l o volumen muerto y a c t i v o
clc
close all
V = i n p u t ( ’ Volumen d e l r e a c t o r= ’ ) ;
v = i n p u t ( ’ F l u j o molar= ’ ) ;
tm = i n p u t ( ’ tiempo medio de r e s i d e n c i a= ’ ) ;
Vm = V−(tm v ) ; *
f p r i n t f ( ’ El volumen muerto d e l r e a c t o r e s : . . .
%10.2 f \n\n ’ ,Vm)
Va = V−Vm;
f p r i n t f ( ’ El volumen a c t i v o d e l r e a c t o r e s : . . .
%10.2 f \n\n ’ , Va)

Código MATLAB para tanque en serie y modulo de

dispersión

75
 clc
close all
f u n c i o n e s = menu ( ’ Tanques en s e r i e y D/uL ’ , . . .
’ Datos Vectores ’ , ’ Funcion ’ )
switch funciones
case 1
t = i n p u t ( ’ Datos de tiempo= ’ ) ;
C = i n p u t ( ’ Datos de c o n c e n t r a c i o n= ’ ) ;
case 2
xn = 0 ;
x l = i n p u t ( ’ Tiempo f i n a l a l que s a l e e l t r a z a d o r \n ’ ) ;
t = [ xn : . 1 : x l ] ;
C = input ( ’ I n t r o d u c i r funcion = ’ ) ;
end
A = t r a p z ( t ,C ) ;
E = C. /A;
*
tE = t . E ;
f = t r a p z ( t , tE ) ;
tmE =(t−f ) . ˆ 2 ;
F =(tmE . E ) ;*
v = trapz ( t ,F) ;
v a r i =(v / ( ( f ) ˆ 2 ) ) ;
N = 1/ v a r i ;
f p r i n t f ( ’ Numero de tanques en s e r i e : . . .
%10.3 f \n\n ’ ,N)
MD = v a r i / 2 ;
f p r i n t f ( ’ Modulo de d i s p e r s i o n : %10.3 f \n\n ’ ,MD)

Código MATLAB para conversión de los diferentes

modelos para reacciones de primer orden

%CONVERSION MODELOS NO IDEAL PRIMER ORDEN

primerorden = menu ( ’CONVERSION’ , . . .
76
 ’ Moedelo tanque i d e a l ’ , ’ Modelo tanque en s e r i e ’ , . . .
’ Modelo de s e g r e g a c i o n ’ , ’ Modelo de d i s p e r s i o n ’ , . . .
’ Modelo de PFR i d e a l ’ )
s w i t c h primerorden
case 1
d i s p ( ’MODELO CSTR IDEAL ’ )
k = i n p u t ( ’ Constante de v e l o c i d a d= ’ ) ;
tm = i n p u t ( ’ Tiempo de r e s i d e n c i a= ’ ) ;
*
X = 1 −(1/(1+(k tm ) ) ) ;
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
case 2
d i s p ( ’MODELO TANQUE EN SERIE ’ )
k = i n p u t ( ’ Constante de v e l o c i d a d= ’ ) ;
tm = i n p u t ( ’ Tiempo de r e s i d e n c i a= ’ ) ;
N = i n p u t ( ’ Numero de tanques= ’ ) ;
t i = tm/N;
*
X = 1 −(1/(1+(k t i ) ) ˆN) ;
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
case 3
d i s p ( ’MODELO DE SEGREGACION’ )
k=i n p u t ( ’ Constante de v e l o c i d a d= ’ )
t=i n p u t ( ’ Datos de tiempo= ’ )
C=i n p u t ( ’ Datos de c o n c e n t r a c i o n= ’ )
Area = t r a p z ( t ,C ) ;
E = C. / Area ;
e = exp(−k . t ) *
Ee = e . E *
tEe = t r a p z ( t , Ee )
X = 1−tEe
case 4
d i s p ( ’MODELO DE DISPERSION ’ )
k = i n p u t ( ’ Constante de v e l o c i d a d= ’ ) ;
tm = i n p u t ( ’ Tiempo de r e s i d e n c i a= ’ ) ;

77
Pe = i n p u t ( ’ I n t r o d u c i r numero de Pe= ’ )
Da = tm k *
*
q =(1+((4 Da)/ Pe ) ) ˆ . 5
** * *
X = 1 −((4 q exp ( Pe / 2 ) ) / ( ( ( ( 1 + q ) ˆ 2 ) exp ( Pe ( q / 2 ) ) . . .
* *
−(((1−q ) ˆ 2 ) exp((−Pe ( q / 2 ) ) ) ) ) ) )
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
case 5
d i s p ( ’MODELO DE PFR i d e a l ’ )
k = i n p u t ( ’ Constante de v e l o c i d a d= ’ ) ;
tm = i n p u t ( ’ Tiempo de r e s i d e n c i a= ’ ) ;
Da = tm k *
X = 1−exp(−Da)
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
end

Código MATLAB para conversión de los diferentes

modelos para reacciones de segundo orden

%CONVERSION MODELOS NO IDEALES 2DO ORDEN

segunorden = menu ( ’CONVERSION’ , . . .
’ Moedelo tanque i d e a l ’ , ’ Modelo tanque en s e r i e ’ , . . .
’ Modelo de s e g r e g a c i o n ’ , ’ Modelo de d i s p e r s i o n ’ . . .
, ’ Modelo de PFR i d e a l ’ )
s w i t c h segunorden
case 1
d i s p ( ’MODELO CSTR IDEAL ’ )
Da = i n p u t ( ’ Numero de c o n v e c c i o n= ’ ) ;
* *
X = 1 −((((1+(4 Da) ) ˆ . 5 ) − 1 ) / ( 2 Da ) ) ;
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
case 2
d i s p ( ’MODELO TANQUE EN SERIE ’ )
Da = i n p u t ( ’ Numero de c o n v e c c i o n= ’ ) ;
N = i n p u t ( ’ Numero de tanques= ’ ) ;
78
KTCi = Da/N
* *
X = 1 −( ( ( (( 2 ( ( ( 1 +( 4 KTCi) ) ˆ . 5 ) + 1 ) ) + 1 ) ˆ ( 1 /N) ) − 1 ) . . .
*
/(2 KTCi ) ) ;
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
case 3
d i s p ( ’MODELO DE SEGREGACION’ )
Da = i n p u t ( ’ Numero de c o n v e c c i o n= ’ ) ;
C = i n p u t ( ’ Datos de c o n c e n t r a c i o n= ’ )
Area = t r a p z ( t ,C ) ;
E = C. / Area ;
e = Da/(1+Da)
Ee = e . E *
tEe = t r a p z ( t , Ee )
X = 1−tEe
case 4
d i s p ( ’MODELO DE DISPERSION ’ )
Da = i n p u t ( ’ Numero de c o n v e c c i o n= ’ ) ;
Pe = i n p u t ( ’ I n t r o d u c i r numero de Pe= ’ )
Dai = Da/2
*
q = (1+((4 Dai )/ Pe ) ) ˆ . 5
**
X = 1 −((4 q exp ( Pe / 2 ) ) / ( ( ( ( 1 + q ) ˆ 2 ) . . .
* * * *
exp ( Pe ( q/2)) −(((1 −q ) ˆ 2 ) exp((−Pe ( q / 2 ) ) ) ) ) ) )
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
case 5
d i s p ( ’MODELO PFR IDEAL ’ )
Da =i n p u t ( ’ Numero de c o n v e c c i o n= ’ ) ;
X = Da/(1+Da)
f p r i n t f ( ’X: %10.4 f \n\n ’ ,X)
end

Código MATLAB para conversión en Modelo de

mezclado máximo

79
f u n c t i o n ODE
clc , close a l l
t0 = i n p u t ( ’ Tiempo f i n a l= ’ ) ;
t f = i n p u t ( ’ Tiempo i n i c i a l = ’ ) ;
p e r i o d o = [ t0 t f ] ;
X0 = 0 ;
x0 = X0 ;
[ t , x ] = ode45 (@( t , x ) modelo ( t , x ) , p e r i o d o , x0 ) ;
figure (1)

p l o t ( t , x ( : , 1 ) , ’m’ , ’ LineWidth ’ , 1 )
s e t ( gca , ’ f o n t s i z e ’ , 1 2 )
$ $
x l a b e l ( ’ \ theta ’ , ’ i n t e r p r e t e r ’ , ’ l a t e x ’ , . . .
’ fontsize ’ ,16)
$$
ylabel ( ’ X ’ , ’ interpreter ’ , ’ latex ’ , . . .
’ fontsize ’ ,16)
g r i d on
g r i d minor
hold on
end

f u n c t i o n dx = modelo ( t , x )
% V a r i a b l e s de e s t a d o
X = x(1);
% Parametros
N = 5 % Numero de Tanque
D = 10 % Numero de c o n v e c c i o n
B = f a c t o r i a l (N−1)
*
l = ( ( t ˆ(N−1)) (NˆN) ) /B
i f N == 1
sm = 1
e l s e i f N˜=1
sumy = 1

80
f o r p = 1 : 1 : N−1
b = f a c t o r i a l (p)
*
o = ( ( t ˆp . Nˆp ) . / b )
sumy = sumy + o
end
sm = sumy
end

% Ecuaciones d i f e r e n c i a l e s
* *
dX =(( l . / sm ) . X)−(D. (1 −X) ˆ 2 )
% Vector de r e s u l t a d o s
dx = dX ;
end

81