PROBLEMAS RESUELTOS

DE
INGENIERÍA DE LA REACCIÓN QUÍMICA

TEMA 6

03.- Se pretende caracterizar el flujo en un reactor de laboratorio mediante la realización de un ensayo

estímulo-respuesta inyectando una pulsación de trazador. Para ello se toma el efluente de salida del
reactor en los intervalos de tiempo indicados y se determina la concentración del trazador (CT) en
cada muestra (en unidades arbitrarias). Los resultados obtenidos del ensayo son los siguientes:
Tiempo (s) 0-5 5-10 10-15 15-20 20-25 25-40 40-45 45-60 60-70 70-80 80-90
CT (ua) 0,000 0,025 0,160 0,200 0,220 0,275 0,285 0,300 0,320 0,330 0,300

Tiempo (s) 90-100 100-150 150-200 200-220 220-300 300-340 340-420 420-480 480-520 520-600
CT (ua) 0,285 0,240 0,190 0,180 0,120 0,100 0,060 0,025 0,015 0,005
Estudie la aplicabilidad de los diferentes modelos de flujo no ideal y proponga un modelo adecuado
para caracterizar el flujo.

SOLUCIÓN

En primer lugar, a la hora de procesar los datos para el análisis de la DTR, se debe proceder siempre
del modo más coherente posible, y teniendo en cuenta la forma en la que los datos se han tabulado.
Así, en este problema se debe aprovechar que se dispone de datos de concentración media de trazador
para determinados intervalos de tiempo. Por ello, se puede aplicar directamente el método de integra-
ción numérica de los rectángulos, ajustándolo a dichos intervalos; ya que resultaría mucho más labo-
rioso aplicarlo a otros posibles intervalos que hubiera que definir ad hoc.

Las tablas de las distintas funciones estadísticas calculadas se incluyen en el archivo anexo al final. En
las columnas A y B se presentan los datos originales primarios y, en la columna C, se presenta el valor
de tiempo (t) que se ha asignado a cada valor de concentración (C’) para establecer las curvas. Como
se puede observar, su valor corresponde al instante final de cada intervalo tabulado. Hay que tener en
cuenta que el trazador contenido en cada muestra corresponde al que ha salido en cada intervalo de
tiempo y esos son precisamente los intervalos que se usan para calcular los intervalos de área por el
método de los rectángulos o de los trapecios.

Por otra parte, se puede observar también que se han incluido arbitrariamente dos datos: el primero de
ellos, C = 0 para t = 0, que debe existir siempre en todos los cálculos; y el último, C = 0 para t = 650.
Este último valor ha tenido que ser estimado, ya que la señal de trazador original no termina en el valor
C = 0, como debería corresponder a un experimento en el que hubiera salido todo el trazador del sis-
tema. Para dicha estimación se ha tenido en cuenta que la pendiente que presentan los dos últimos pun-
tos tabulados es de -1,3·10-4 ua/s, y que, extrapolando dicha pendiente, la señal llegaría al valor 0 justo
en el instante de 640 s. Puesto que en realidad la señal debe ir reduciendo su pendiente progresiva-
mente, se ha tomado como valor estimado para el final de la señal de trazador el instante t = 650 s.

Una vez que se dispone de la tabla de valores de C’ versus t debidamente corregida, se puede proceder
a realizar los cálculos numéricos. En las columnas E, G y J se presentan respectivamente las funciones
E(t), F(t) e I(t). En la casilla H29 se presenta el valor del tiempo medio de residencia; en la casilla I29,
el valor de la varianza de residencia; y, en la I37, se calcula la anchura típica de residencia. Asimismo,

1
en la casilla K29 se presenta el valor de edad media en el interior del sistema; en la casilla L29, el va-
lor de la varianza de edad; y, en la L37, se calcula la anchura típica de edad. Por otra parte, en la co-
lumna N se presenta el valor calculado del tiempo reducido (adimensional). En columnas O, Q y T se
presentan las funciones E(θ), F(θ) e I(θ) respectivamente. Finalmente, en las casillas S33 y T33 se pre-
senta el valor de la varianza de residencia reducida, aplicando las dos ecuaciones de cálculo que se han
definido en las clases y que, obviamente, deben conducir al mismo resultado.

Se puede comprobar que todas las funciones cumplen las propiedades necesarias. Por ejemplo, que las
funciones de distribución F(t) y F(θ) comienzan en el valor 0 (marcado en rojo e introducido arbitraria-
mente) y terminan en el valor 1. Por el contrario, la función I(θ) comienza en 1 y termina en 0.

NOTA 1

Se puede comprobar que, si se cambia el valor del último dato de la señal t = 650 s (que ha sido esti-
mado), cambian también los valores de las funciones estadísticas. Incluso puede que cambie el valor
del último dato de la función I(t), de forma que es posible que no sea igual a 0. Por ejemplo, si el úl-
timo dato de la señal se fija en t < 650, el último valor de I(t) no es exactamente 0; mientras que sí vale
0 si se introduce en la tabla cualquier otro valor de t igual o mayor. Por eso se ha tomado precisamente
ese valor de t en la tabla.

FIN DE LA NOTA 1

Una vez que se han realizado todos los cálculos pertinentes, para la aplicación del modelo de disper-
sión axial o del de tanques en serie, interesa conocer el valor correspondiente a la varianza de residen-
cia reducida (casillas S33 ó T33). En este caso, su valor resulta ser de 0,604049.

Modelo de dispersión axial

Según este modelo, debemos suponer inicialmente que el módulo de dispersión axial del sistema es
bajo (D/uL < 0,01); por lo tanto, la ecuación a aplicar para su cálculo sería la siguiente:

𝐷 𝜎𝜃2 0,604049
= = = 0,302034
𝑢𝐿 2 2

Como se puede observar, el valor resultante no es menor de 0,01 como se había supuesto, por lo que
no sería aplicable dicha ecuación. No obstante, en las casillas AC3, AD3 y AE3, se realiza el cálculo
de la conversión para una supuesta reacción de primer orden y densidad constante con kA = 0,01 s-1. El
valor resultante es de 0,729424.

Puesto que el módulo de dispersión ha resultado alto (D/uL > 0,01), pero no demasiado (D/uL < 1), se
pueden aplicar las ecuaciones correspondientes a dicha situación.

𝐷 2 𝐷
𝑎𝑏𝑖𝑒𝑟𝑡𝑜 − 𝑎𝑏𝑖𝑒𝑟𝑡𝑜 8 ( ) + 2 ( ) − 𝜎𝜃2 = 0
𝑢𝐿 𝑢𝐿

𝐷 2 𝐷
𝑎𝑏𝑖𝑒𝑟𝑡𝑜 − 𝑐𝑒𝑟𝑟𝑎𝑑𝑜 𝑜 𝑐𝑒𝑟𝑟𝑎𝑑𝑜 − 𝑎𝑏𝑖𝑒𝑟𝑡𝑜 3 ( ) + 2 ( ) − 𝜎𝜃2 = 0
𝑢𝐿 𝑢𝐿

𝐷 2 𝐷 𝐷
𝑐𝑒𝑟𝑟𝑎𝑑𝑜 − 𝑐𝑒𝑟𝑟𝑎𝑑𝑜 2 ( ) (1 − 𝑒 −1/(𝑢𝐿) ) + 2 ( ) − 𝜎𝜃2 = 0
𝑢𝐿 𝑢𝐿

2
Los módulos resultantes se muestran en las casillas AA6, AA9 y AA17 respectivamente. En el último
caso es necesario realizar un cálculo iterativo, que se ha comenzado utilizando el valor del módulo ob-
tenido anteriormente cuando se supuso su valor bajo.

𝐷 𝐷 𝐷
| = 0,176879 | = 0,225649 | = 0,311004
𝑢𝐿 𝑎/𝑎 𝑢𝐿 𝑎/𝑐−𝑎/𝑜 𝑢𝐿 𝑐/𝑐

Puesto que no se ha indicado si el experimento se ha realizado en un sistema de tipo a/a, a/c, c/a o c/c,
en las casillas AE6, AE9 y AE17 se realiza el cálculo de la conversión para la supuesta reacción que se
ha indicado antes. Como se puede observar, en todos los casos se alcanza una conversión que se en-
cuentra entre 0,72 y 0,76.

Modelo de tanque en serie

Según este modelo, el valor del número de tanques se calcula según la siguiente ecuación:

1 1
𝑁= 2 = 0,604049 = 1,655495 ≈ 2
𝜎𝜃

El resultado no es lo suficientemente alto como para aplicar el modelo con precisión (N > 50); ni si-
quiera es lo suficientemente alto como para aplicarlo con reservas (N > 30). No obstante, en la casilla
AE27 se realiza el cálculo de la conversión que resultaría para un sistema de 2 RFMC en serie, con la
supuesta reacción antes indicada.

𝐴 1 1
𝑥𝐴 = 1 − = 1− =1− = 0,711465
𝐴𝑜 𝑘𝑡̅
𝑁
0,01 · 172,3324 2
1 + (𝑁) 1+( )
2

Como se puede observar, dado que el número de tanque no ha resultado lo suficientemente alto, la
conversión estimada resulta inferior a la calculada antes.

Integral combinada

Puesto que se está suponiendo una reacción de orden 1 (lineal), se puede estimar también la conversión
mediante la integral combinada. Así, en la columna X se calcula la conversión alcanzada para cada
tiempo medido. Igualmente, en la columna Y se calcula el valor del incremento de área para cada in-
tervalo de la integral combinada. Finalmente, en la casilla Y29 se calcula la conversión media resul-
tante (x̅A = 0,662901).

Como se puede observar, la conversión resultante es inferior a la estimada con los dos modelos ante-
riores (066 < 0,71); lo que indica que el sistema no se simula bien suponiendo que se trata de un reac-
tor tubular cercano al RFP. En realidad, al ser la conversión de la integral combinada menor que la de
dichos moldeos, se debe esperar que el sistema se parezca más a un reactor de tanque agitado, más cer-
cano a un RFMC.

Modelo de Cholette-Cloutier

Para este modelo, se debe representar ln I(θ) frente a θ. Como resultado, se debería obtener una orde-
nada en el origen (b) negativa, relacionada con el grado de cortocircuito; y una pendiente (a) también
negativa, relacionada con el grado de volumen muerto, además de con el parámetro anterior. De este
modo:
3
 𝑉𝑆𝐼𝑆 𝑄𝑀𝐶 𝜏𝑆𝐼𝑆 𝑄𝐴
ln 𝐼(𝜃) = 𝑎 𝜃 + 𝑏 𝑎=− · =− 𝑏 = ln ( )
𝑉𝑀𝐶 𝑄𝑆𝐼𝑆 𝜏𝑀𝐶 𝑄𝑆𝐼𝑆

𝑉𝑆𝐼𝑆 = 𝑉𝑀𝐶 + 𝑉𝑉𝑀 𝑄𝑆𝐼𝑆 = 𝑄𝑀𝐶 + 𝑄𝐶𝐶

La tabla de datos correspondiente se muestra en las columnas AG, AH y AI. La representación de di-
chos datos se muestra en la figura que se incluye al final (Figura 1). Como se puede observar, aunque
los datos no se alinean de forma muy ajustada a una recta, en todo caso, no se obtiene de ningún modo
una ordenada en el origen negativa. En consecuencia, no resulta de aplicación este modelo.

No obstante, si se hubiera obtenido un resultado adecuado, la estimación de la conversión debería rea-

lizarse conforme a las siguientes ecuaciones de diseño:

t̅
𝜏𝑆𝐼𝑆 t̅ 𝑘 𝜏𝑀𝐶 −𝑘 𝑎 𝑄𝐴
𝜏𝑀𝐶 =− =− 𝑥𝐴1 = = 𝑥𝐴2 = · 𝑥 = 𝑒 𝑏 𝑥𝐴1
𝑎 𝑎 1 + 𝑘 𝜏𝑀𝐶 t̅ 𝑄𝑆𝐼𝑆 𝐴1
1−𝑘𝑎

t̅
−𝑘 𝑎
𝑥𝐴2 = 𝑒 𝑏
t̅
1−𝑘𝑎

El resultado que se hubiera obtenido aplicando el ajuste se presenta en las casillas AL41 y AL42, re-
sultando una conversión final de xA2 = 0,734356.

Modelo de Hovorka-Adler

Para este modelo, se debe representar también ln I(θ) frente a θ. Como resultado se debería obtener una
ordenada en el origen (b) positiva, relacionada con el tamaño de los dos reactores ideales; y una pen-
diente (a) negativa, relacionada con el tamaño del RFMC. En este modelo, el valor absoluto de la pen-
diente (a) debe ser superior a 1, puesto que el volumen del RFMC no puede ser superior al volumen de
todo el sistema. Así:

𝑉𝐹𝑃 𝑉𝑆𝐼𝑆
ln 𝐼(𝜃) = 𝑎 𝜃 + 𝑏 𝑏= 𝑎=−
𝑉𝑀𝐶 𝑉𝑀𝐶

𝑉𝑆𝐼𝑆
𝑉𝑀𝐶 = − 𝑉𝐹𝑃 = 𝑏 · 𝑉𝑀𝐶 𝑉𝑆𝐼𝑆 = 𝑉𝑀𝐶 + 𝑉𝐹𝑃 + 𝑉𝑉𝑀
𝑎

La representación de los datos es la misma que en el modelo anterior (Figura 1). Como se puede obser-
var, sería posible ajustar los datos a este modelo con suficiente precisión (r2 > 0,99), sólo si se descar-
tan algunos de los datos de los extremos de la curva; es decir, del principio y del final de la tabla. En
este caso, se han descartado los datos marcados en rojo, que corresponden a los valores de θ < 0,2 y de
θ > 2,0. Por lo tanto se han conservado sólo los datos de θ entre 0,2 y 2,0. El ajuste de dichos datos
(0,2 < θ < 2,0) ofrece el siguiente resultado:

𝑉𝑆𝐼𝑆
𝑟 2 = 0,992325 𝑏 = 0,232817 𝑎=− = −1,07972
𝑉𝑀𝐶

4
 𝑉𝐹𝑃
𝑉𝑀𝐶 1 1 𝑏 𝑉𝑀𝐶 𝑉𝐹𝑃 0,232817
= = = 0,926163 = = = = 0,215626
𝑉𝑆𝐼𝑆 −𝑎 1,07972 −𝑎 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 1,07972
𝑉𝑀𝐶

𝑉𝑉𝑀 𝑉𝑆𝐼𝑆 − 𝑉𝑀𝐶 − 𝑉𝐹𝑃 𝑉𝑀𝐶 𝑉𝐹𝑃

𝑉𝑉𝑀 = 𝑉𝑆𝐼𝑆 − 𝑉𝑀𝐶 − 𝑉𝐹𝑃 = =1− −
𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆

𝑉𝑉𝑀
= 1 − 0,926163 − 0,215626 = −0,14179
𝑉𝑆𝐼𝑆

Se observa que, a pesar de la bondad del ajuste, no se obtienen unos valores coherentes de los diferen-
tes volúmenes del modelo, ya que se obtiene un volumen muerto negativo. Sin embargo, si se hubiera
obtenido un resultado adecuado, la estimación de la conversión debería realizarse conforme a las si-
guientes ecuaciones de diseño:

𝑉𝑀𝐶
𝜏𝑀𝐶 𝜏𝑀𝐶 𝑄𝑀𝐶 𝑉𝑀𝐶 1 𝜏𝑀𝐶 1 t̅
= = = = = 𝜏𝑀𝐶 = = 159,6079 𝑠
𝜏𝑆𝐼𝑆 t̅ 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 −𝑎 t̅ −𝑎 −𝑎
𝑄𝑆𝐼𝑆

𝑉𝐹𝑃
𝜏𝐹𝑃 𝜏𝐹𝑃 𝑄 𝑉𝐹𝑃 𝑏 𝜏𝐹𝑃 𝑏 t̅ 𝑏
= = 𝐹𝑃 = = = 𝜏𝐹𝑃 = = 37,1594 𝑠
𝜏𝑆𝐼𝑆 t̅ 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 −𝑎 t̅ −𝑎 −𝑎
𝑄𝑆𝐼𝑆

− 𝑘 𝜏𝐹𝑃
𝑥𝐴1 + 𝑘 𝜏𝑀𝐶 (1 − 𝑒 − 𝑘 𝜏𝐹𝑃 ) + 𝑘 𝜏𝑀𝐶
𝑥𝐴1 = 1 − 𝑒 𝑥𝐴2 = 𝑥𝐴2 =
1 + 𝑘 𝜏𝑀𝐶 1 + 𝑘 𝜏𝑀𝐶

0,310366 + (0,01 · 159,6079)

𝑥𝐴1 = 1 − 𝑒 − 0,01·37,1594 = 0,310366 𝑥𝐴2 = = 0,734356
1 + (0,01 · 159,6079)

Se observa que la conversión estimada según este modelo (xA = 0,73) resulta mayor que la estimada
mediante la integral combinada; por lo que se confirma que este modelo tampoco sería completamente
asumible.

NOTA 2

Si se hubiera aplicado el ajuste de este último modelo al conjunto completo de datos experimentales
sin descartar ninguno (Figura 2), los resultados hubieran sido los siguientes:

𝑉𝑆𝐼𝑆
𝑟 2 = 0,932313 𝑏 = 0,394268 𝑎=− = −1,55495
𝑉𝑀𝐶

𝑉𝐹𝑃
𝑉𝑀𝐶 1 1 𝑏 𝑉𝑀𝐶 𝑉𝐹𝑃 0,232817
= = = 0,643108 = = = = 0,253557
𝑉𝑆𝐼𝑆 −𝑎 1,55495 −𝑎 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 1,55495
𝑉𝑀𝐶

𝑉𝑉𝑀 𝑉𝑆𝐼𝑆 − 𝑉𝑀𝐶 − 𝑉𝐹𝑃 𝑉𝑀𝐶 𝑉𝐹𝑃

𝑉𝑉𝑀 = 𝑉𝑆𝐼𝑆 − 𝑉𝑀𝐶 − 𝑉𝐹𝑃 = =1− −
𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆 𝑉𝑆𝐼𝑆

5
 𝑉𝑉𝑀
= 1 − 0,643108 − 0,253557 = 0,103335
𝑉𝑆𝐼𝑆

Se puede observar que ahora sí se obtienen valores coherentes de los diferentes volúmenes del modelo,
aunque la regresión resulta muy baja (r2 = 0,93). En este caso, si se realiza la estimación de la conver-
sión con dichos coeficientes de ajuste, el resultado es el siguiente:

t̅ 172,3324 t̅ 𝑏 172,3324 · 0,394268

𝜏𝑀𝐶 = = = 110,8283 𝑠 𝜏𝐹𝑃 = = = 43,6961 𝑠
−𝑎 1,55495 −𝑎 1,55495

(1 − 𝑒 − 𝑘 𝜏𝐹𝑃 ) + 𝑘 𝜏𝑀𝐶 (1 − 𝑒 − 0,01·43,6961 ) + (0,01 · 110,8283)

𝑥𝐴2 = = = 0,693591
1 + 𝑘 𝜏𝑀𝐶 1 + (0,01 · 110,8283)

Se puede observar que este resultado se aproxima bastante al que se obtuvo mediante la integral com-
binada (x̅A = 0,662901). Por lo tanto, se puede decir que a pesar de que la regresión obtenida es muy
baja (r2 < 0,95)), se puede asumir el modelo si se usa el conjunto completo de datos.

FIN DE LA NOTA 2

6
Figura 1

Figura 2

7
 2
intervalo C' t Δa[C'] C(t) = E(t) Δa[E(t)] γE(t) = F(t) t·E(t)·Δt t ·E(t)·Δt I(t) t·I(t)·Δt t2·I(t)·Δt
s ua s (ua)·s 1/s - - s s2 1/s - -
0 0 0 0 0,005803
0-5 0 5 0 0 0 0 0 0 0,005803 0,145068 0,725342
5-10 0,025 10 0,0625 0,000316 0,000791 0,000791 0,015813 0,158128 0,005798 0,289908 2,899076
10-15 0,16 15 0,4625 0,002024 0,005851 0,006641 0,151803 2,27704 0,005764 0,432315 6,484727
15-20 0,2 20 0,9 0,00253 0,011385 0,018027 0,253004 5,060089 0,005698 0,569814 11,39627
20-25 0,22 25 1,05 0,002783 0,013283 0,031309 0,347881 8,697027 0,005621 0,702633 17,56581
25-40 0,275 40 3,7125 0,003479 0,046964 0,078273 2,087287 83,49146 0,005349 3,209124 128,365
40-45 0,285 45 1,4 0,003605 0,01771 0,095984 0,811195 36,5038 0,005246 1,180299 53,11344
45-60 0,3 60 4,3875 0,003795 0,055503 0,151486 3,41556 204,9336 0,004924 4,431333 265,88
60-70 0,32 70 3,1 0,004048 0,039216 0,190702 2,83365 198,3555 0,004696 3,287302 230,1111
70-80 0,33 80 3,25 0,004175 0,041113 0,231815 3,339658 267,1727 0,004458 3,566061 285,2849
80-90 0,3 90 3,15 0,003795 0,039848 0,271664 3,41556 307,4004 0,004226 3,803713 342,3341
90-100 0,285 100 2,925 0,003605 0,037002 0,308665 3,605313 360,5313 0,004012 4,011635 401,1635
100-150 0,24 150 13,125 0,003036 0,166034 0,4747 22,7704 3415,56 0,003048 22,86136 3429,205
150-200 0,19 200 10,75 0,002404 0,13599 0,610689 24,03542 4807,084 0,002259 22,59068 4518,136
200-220 0,18 220 3,7 0,002277 0,046806 0,657495 10,01898 2204,175 0,001987 8,74485 1923,867
220-300 0,12 300 12 0,001518 0,151803 0,809298 36,43264 10929,79 0,001107 26,55827 7967,481
300-340 0,1 340 4,4 0,001265 0,055661 0,864959 17,2043 5849,462 0,000784 10,65707 3623,405
340-420 0,06 420 6,4 0,000759 0,080961 0,94592 25,50285 10711,2 0,000314 10,54403 4428,493
420-480 0,025 480 2,55 0,000316 0,032258 0,978178 9,108159 4371,917 0,000127 3,646807 1750,468
480-520 0,015 520 0,8 0,00019 0,01012 0,988299 3,946869 2052,372 6,79E-05 1,41233 734,4118
520-600 0,005 600 0,8 6,33E-05 0,01012 0,998419 3,036053 1821,632 9,18E-06 0,440436 264,2614
0 650 0,125 0 0,001581 1 0 0 0 0 0

A A A = t̅ A A = t̅ A
2
(ua)·s (ua)·s s s s s2
79,05 1 172,3324 47637,77 133,085 30385,05
residencia edad
σ2 σ2
0,0045
C(t) = E(t) s2 s2
0,004 17939,32 12673,42

0,0035
±2s ±2s
0,003 s s
267,8755 225,1526
0,0025

0,002

0,0015

0,001

0,0005

0
0 100 200 300 400 500 600 700
 θ E(θ) Δa[E(θ)] γE(θ) = F(θ) θ·E(θ)·Δθ θ2·E(θ)·Δθ 2
I(θ) θ·I(θ)·Δθ θ ·I(θ)·Δθ k xA(t) xA·E(t)Δt
- - - - - - - - - 1/s - -
0 0 0 1 0,01 0
0,029014 0 0 0 0 0 1 0,000842 2,44E-05 0,048771 0
0,058027 0,054501 0,000791 0,000791 9,18E-05 5,32E-06 0,999209 0,001682 9,76E-05 0,095163 0,00015
0,087041 0,348807 0,005851 0,006641 0,000881 7,67E-05 0,993359 0,002509 0,000218 0,139292 0,00141
0,116055 0,436009 0,011385 0,018027 0,001468 0,00017 0,981973 0,003306 0,000384 0,181269 0,002293
0,145068 0,479609 0,013283 0,031309 0,002019 0,000293 0,968691 0,004077 0,000591 0,221199 0,003078
0,23211 0,599512 0,046964 0,078273 0,012112 0,002811 0,921727 0,018622 0,004322 0,32968 0,017203
0,261123 0,621312 0,01771 0,095984 0,004707 0,001229 0,904016 0,006849 0,001788 0,362372 0,006532
0,348164 0,654013 0,055503 0,151486 0,01982 0,0069 0,848514 0,025714 0,008953 0,451188 0,025684
0,406192 0,697614 0,039216 0,190702 0,016443 0,006679 0,809298 0,019075 0,007748 0,503415 0,020379
0,464219 0,719414 0,041113 0,231815 0,019379 0,008996 0,768185 0,020693 0,009606 0,550671 0,022988
0,522247 0,654013 0,039848 0,271664 0,01982 0,010351 0,728336 0,022072 0,011527 0,59343 0,022521
0,580274 0,621312 0,037002 0,308665 0,020921 0,01214 0,691335 0,023278 0,013508 0,632121 0,02279
0,870411 0,52321 0,166034 0,4747 0,132131 0,115008 0,5253 0,132659 0,115467 0,77687 0,117931
1,160548 0,414208 0,13599 0,610689 0,139471 0,161863 0,389311 0,131088 0,152134 0,864665 0,103913
1,276603 0,392408 0,046806 0,657495 0,058138 0,074219 0,342505 0,050744 0,06478 0,889197 0,040495
1,740822 0,261605 0,151803 0,809298 0,211409 0,368026 0,190702 0,154111 0,268279 0,950213 0,115396
1,972932 0,218004 0,055661 0,864959 0,099832 0,196962 0,135041 0,06184 0,122007 0,966627 0,048912
2,437151 0,130803 0,080961 0,94592 0,147986 0,360665 0,05408 0,061184 0,149115 0,985004 0,059811
2,785315 0,054501 0,032258 0,978178 0,052852 0,14721 0,021822 0,021161 0,058941 0,99177 0,018819
3,017425 0,032701 0,01012 0,988299 0,022903 0,069107 0,011701 0,008195 0,024729 0,994483 0,007548
3,481644 0,0109 0,01012 0,998419 0,017617 0,061338 0,001581 0,002556 0,008898 0,997521 0,005048
3,771781 0 0,001581 1 0 0 0 0 0 0,998497 0

A A=θ A A=θ A A = x̅A

- - - - - -
1 1 1,604049 0,772258 1,023119 0,662901
residencia edad
1,2 σ2 σ2
I(θ) - -
0,604049 0,604049 0,426737 0,426737
1

±2s ±2s
0,8 - -
1,554412 1,306502
0,6

0,4

0,2

0
0 1 2 3 4
MODELO DE DISPERSIÓN AXIAL θ I(θ) ln I(θ)
D/uL pequeño a CA/CAo xA - - - 0
0 1 2 3 4
0,302024 >0,01 1,755547 0,270576 0,729424 0 1 0 -0,5
0,029014 1 0 y = -1,0797x + 0,2328
D/uL grande aa a CA/CAo xA 0,058027 0,999209 -0,00079 -1 R² = 0,9923

0,176879 1,489725 0,243152 0,756848 0,087041 0,993359 -0,00666

-1,5
0,116055 0,981973 -0,01819
D/uL grande ca-ac a CA/CAo xA 0,145068 0,968691 -0,03181 -2
0,225649 1,598581 0,254949 0,745051 0,23211 0,921727 -0,08151
-2,5
0,261123 0,904016 -0,10091
D/uL grande cc a CA/CAo xA 0,348164 0,848514 -0,16427 -3
0,302024 itera 0 1,755547 0,270576 0,729424 0,406192 0,809298 -0,21159
-3,5
0,310498 itera 1 1,772104 0,272136 0,727864 0,464219 0,768185 -0,26373
0,310976 itera 2 1,773033 0,272223 0,727777 0,522247 0,728336 -0,31699 -4
0,311003 itera 3 1,773086 0,272228 0,727772 0,580274 0,691335 -0,36913
0,311004 itera 4 1,773089 0,272228 0,727772 0,870411 0,5253 -0,64378 -4,5
0,311004 itera 5 1,773089 0,272228 0,727772 1,160548 0,389311 -0,94338 1
1,276603 0,342505 -1,07147
módulo de dispersión axial xA 1,740822 0,190702 -1,65704 0
0,176879 min min 0,727772 1,972932 0,135041 -2,00218 0 1 2 3 4

0,311004 max max 0,756848 2,437151 0,05408 -2,9173 -1

y = -1,5549x + 0,3943
0,243942 med med 0,74231 2,785315 0,021822 -3,82485 R² = 0,9323
-2
3,017425 0,011701 -4,44804
MODELO DE TANQUES EN SERIE 3,481644 0,001581 -6,44952
-3
N xA 3,771781 0
1,655495 <30 <<50 0,306949 0,693051 -4
2 0,288535 0,711465 a b
-5
-1,07972 0,232817
0,030028 0,029891
-6
r2 > 0,99 0,992325 0,058756
1292,9 10 -7
4,463498 0,034523

MODELO DE HOVORKA-ADLER MODELO DE CHOLETTE-CLOUTIER

VSIS (L) 100 QSIS (L/s) 100
VMC (L) 92,61633 a = -VSIS/VMC QMC (L/s) 126,215 b = ln(QMC/QSIS)
VFP (L) 21,56263 b = VFP/VMC QCC (L/s) -26,215
VVM (L) -14,179 VMC (L) 116,8957 a = -(VSIS/VMC)(QMC/QSIS) = τSIS/τMC
τFP (s) 37,1594 t ̅ b/a VVM (L) -16,8957
τMC (s) 159,6079 t/̅ a τMC (s) 159,6079 τSIS/a = t/̅ -a
xA1 0,310366 xA1 0,614804
xA2 0,734356 xA2 0,775975

a b
-1,55495 0,394268
0,093686 0,136732
r2 < 0,95 0,932313 0,467269
275,476 20
60,14748 4,366803

MODELO DE HOVORKA-ADLER MODELO DE CHOLETTE-CLOUTIER

VSIS (L) 100 QSIS (L/s) 100
VMC (L) 64,31079 a = -VSIS/VMC QMC (L/s) 148,3298 b = ln(QMC/QSIS)
VFP (L) 25,35569 b = VFP/VMC QCC (L/s) -48,3298
VVM (L) 10,33353 VMC (L) 95,39206 a = -(VSIS/VMC)(QMC/QSIS) = τSIS/τMC
τFP (s) 43,69606 t ̅ b/a VVM (L) 4,607936
τMC (s) 110,8283 t/̅ a τMC (s) 110,8283 τSIS/a = t/̅ -a
xA1 0,354003 xA1 0,52568
xA2 0,693591 xA2 0,779741