Distribución de tiempos de residencia Índice Teoría Determinación Experimental del RDT Referencias Menú de navegación10.1016/0009-2509(53)80001-1
Ingeniería química
reactor químicofluidoingenieros químicosMacMullinWeber1935fluido compresibleasimetríacurtosisdeconvolucionarDirac delta
La distribución de tiempos de residencia (RTD) de un reactor químico en un fluido, es una distribución de probabilidad que describe la cantidad de tiempo que un elemento diferencial de un fluido puede pasar 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.
El concepto fue propuesto por MacMullin y Weber en 1935, pero no fue usado intensivamente hasta que P.V. Danckwerts analizó un número importante de RTDs en 1953.
Índice
1 Teoría
2 Determinación Experimental del RDT
2.1 Experimento en pulso
2.2 Experimento en escalón
3 Referencias
Teoría
Está basada en 3 suposiciones principales:
- Sistema en estado estacionario.
- Transporte en la entrada y salida tiene lugar sólo por advección.
- El fluido es incompresible.
La última suposición no es requerida pero un fluido compresible es más difícil de modelar y manipular, además de ser menos común en procesos industriales. Se requiere un nivel de complejidad superior para modelar reactores heterogéneos donde el RTD describe el flujo de cada fase, por ejemplo en los sistemas de borboteo.
La función de distribución de tiempos de residencia es representada por una distribución de edad de salida, E(t)displaystyle E(t). La función E(t)displaystyle E(t) está normalizada y se define matemáticamente:
∫0∞E(t)dt=1.displaystyle int _0^infty E(t),dt=1.
La fracción del fluido que tiene una duración de tiempo tdisplaystyle t dentro del reactor, está dada por el valor de
E(t)dtdisplaystyle E(t)dt.
La fracción del fluido que deja el reactor con una edad inferior a t1displaystyle t_1 es
∫0t1E(t)dt.displaystyle int _0^t_1E(t),dt.
La fracción del fluido que deja el reactor con una edad mayor que t1displaystyle t_1 es
∫t1∞E(t)dt=1−∫0t1E(t)dt.displaystyle int _t_1^infty E(t),dt=1-int _0^t_1E(t),dt.
El tiempo de residencia promedio se calcula con una ecuación de momentos:
t¯=∫0∞t⋅E(t)dt.displaystyle bar t=int _0^infty tcdot E(t),dt.
Si no hay zonas muertas o estancamientos dentro del reactor, entonces t¯displaystyle bar t será igual al tiempo de residencia teórico τdisplaystyle tau , éste se define como la relación entre el volumen del reactor y el flujo volumétrico de entrada del fluido:
τ=Vv.displaystyle tau =frac Vv.
Ecuaciones de momentos de grado superior pueden proveer información sobre la conducta del la función E(t)displaystyle E(t). Por ejemplo, el segundo momento central indica la varianza (σ2)displaystyle (sigma ^2), El grado de dispersión comparado con la media.
σ2=∫0∞(t−t¯)2⋅E(t)dtdisplaystyle sigma ^2=int _0^infty (t-bar t)^2cdot E(t),dt
El tercer momento central indica el grado de asimetría de la distribución E(t)displaystyle E(t) y el cuarto momento central indica el grado de curtosis presente.
Se puede definir una distribución integral de edad I(t)displaystyle I(t) que describe el contenido del reactor. Su definición es un tanto similar a E(t)displaystyle E(t): la fracción de fluido dentro del reactor con una edad tdisplaystyle t es I(t)dtdisplaystyle I(t)dt. Demostrado por Danckwerts, la relación entre E(t)displaystyle E(t) y I(t)displaystyle I(t) puede ser encontrada mediante un balance de masa:
I(t)=1τ(1−∫0tE(t) dt)E(t)=−τdI(t)dtdisplaystyle I(t)=frac 1tau left(1-int _0^tE(t) dtright)qquad E(t)=-tau frac dI(t)dt
Determinación Experimental del RDT
El RTD es medido 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, viscosidad) y la introducción del trazador no debe modificar las condiciones hidrodinámicas.
En general, el cambio en la concentración del trazador será entre una función pulso o una función escalón. Es posible usar otras funciones pero se requiere mas cálculos para deconvolucionar la curva RTD, E(t).
Experimento en pulso
El método requiere la introducción de un volumen muy pequeño del trazador muy concentrado al interior del reactor, tal que su función se aproxime a la función Dirac delta. Debido a que una función infinitamente corta no puede ser producida, esta suele ser producida de tal modo que sea mucho más pequeña que el tiempo de residencia del recipiente. Si una masa del trazador, M, es introducida dentro del recipiente de volumen V y un tiempo de residencia esperado τdisplaystyle tau , la curva resultante de C(t)displaystyle C(t) puede ser transformada en una curva de tiempo de de residencia adimensional mediante la siguiente relación:
E(t)=C(t)∫0∞C(t) dtdisplaystyle E(t)=frac C(t)int _0^infty C(t) dt
Experimento en escalón
En un experimento en escalón la concentración del trazador a la entrada cambia abrupta-mente en un tiempo dado t desde 0 a C0displaystyle C_0. La concentración del trazador a la salida se mide y normaliza a la concentración C0displaystyle C_0 para obtener la curva no-dimensional F(t) cuyos valores rondan entre 0 a 1:
F(t)=C(t)C0.displaystyle F(t)=frac C(t)C_0.
Las respuestas de salida de la función escalón y función pulso de un reactor se relacionan matemáticamente:
F(t)=∫0tE(t)dtE(t)=dF(t)dtdisplaystyle F(t)=int _0^tE(t),dtqquad E(t)=frac dF(t)dt
El valor del tiempos de residencia medio y la varianza puede ser deducidos de la función F(t)displaystyle F(t):
t¯=∫0∞t⋅E(t)dt=∫01tdF(t)=−∫01td(1−F(t))=∫0∞(1−F(t))dtdisplaystyle overline t=int _0^infty tcdot E(t),dt=int _0^1t,dF(t)=-int _0^1t,d(1-F(t))=int _0^infty (1-F(t)),dt
σ2=∫0∞(t−t¯)2⋅E(t)dt=∫01(t−t¯)2dF(t)=∫01t2dF(t)−t¯2=2∫0∞t(1−F(t))dt−t¯2displaystyle sigma ^2=int _0^infty (t-overline t)^2cdot E(t),dt=int _0^1(t-overline t)^2,dF(t)=int _0^1t^2,dF(t)-overline t^2=2int _0^infty t(1-F(t)),dt-overline t^2
Un experimento en escalón frecuentemente es más sencillo de realizar que un experimento en pulso, pero tiende a difuminar algunos de los detalles que una respuesta en pulso suele mostrar. Es fácil integrar numéricamente un experimento en pulso y suele obtenerse un estimado de alta calidad, lo contrario ocurre si se quiere obtener una función inversa mediante diferenciación debido a que cualquier ruido en la medición de la concentración se amplifica por la diferenciación numérica.
Referencias
R.B. MacMullin and M. Weber (1935). «The theory of short-circuiting in continuous-flow mixing vessels in series and kinetics of chemical reactions in such systems». Transactions of American Institute of Chemical Engineers 31 (2): 409–458.
P.V. Danckwerts (1953). «Continuous flow systems. Distribution of residence times.». Chemical Engineering Science 2: 1–13. doi:10.1016/0009-2509(53)80001-1.
Levenspiel, Octave (1999). Chemical Reaction Engineering (3rd edición). John Wiley & Sons. ISBN 0-471-25424-X.
Nauman, E. Bruce (2004). «Residence Time Distributions». Handbook of Industrial Mixing: Science and Practice. Wiley Interscience. pp. 1–17. ISBN 0-471-26919-0.