# Residence time distribution

The residence time distribution (RTD) of a chemical reactor is a probability distribution function that describes the amount of time a fluid element could spend inside the reactor. Chemical engineers use the RTD to characterize the mixing and flow within reactors and to compare the behavior of real reactors to their ideal models. This is useful, not only for troubleshooting existing reactors, but in estimating the yield of a given reaction and designing future reactors.

The concept was first proposed by MacMullin and Weber in 1935, but was not used extensively until P.V. Danckwerts analyzed a number of important RTDs in 1953.

## TheoryEdit

An RTD curve for a reasonably well-mixed reactor

The theory of residence time distributions generally begins with three assumptions:

1. the reactor is at steady-state,
2. transports at the inlet and the outlet takes place only by advection, and
3. the flow is incompressible.

The incompressibility assumption is not required, but compressible flows are more difficult to work with and less common in chemical processes. A further level of complexity is required for multi-phase reactors, where a separate RTD will describe the flow of each phase, for example bubbling air through a liquid.

The distribution of residence times is represented by an exit age distribution, ${\displaystyle E(t)}$ . The function ${\displaystyle E(t)}$  has the units of time−1 and is defined such that

${\displaystyle \int _{0}^{\infty }E(t)\,dt=1}$ .

The fraction of the fluid that spends a given duration, ${\displaystyle t}$  inside the reactor is given by the value of ${\displaystyle E(t)dt}$ .

The fraction of the fluid that leaves the reactor with an age less than ${\displaystyle t_{1}}$  is

${\displaystyle \int _{0}^{t_{1}}E(t)\,dt}$ .

The fraction of the fluid that leaves the reactor with an age greater than ${\displaystyle t_{1}}$  is

${\displaystyle \int _{t_{1}}^{\infty }E(t)\,dt=1-\int _{0}^{t_{1}}E(t)\,dt}$ .

The average residence time is given by the first moment of the age distribution:

${\displaystyle {\bar {t}}=\int _{0}^{\infty }t\cdot E(t)\,dt}$ .

If there are no dead, or stagnant, zones within the reactor then ${\displaystyle {\bar {t}}}$  will be equal to ${\displaystyle \tau }$ , the residence time calculated from the total reactor volume and the volumetric flow rate of the fluid:

${\displaystyle \tau ={\frac {V}{v}}}$ .

The higher order central moments can provide significant information about the behavior of the function ${\displaystyle E(t)}$ . For example, the second central moment indicates the variance (${\displaystyle \sigma ^{2}}$ ), the degree of dispersion around the mean.

${\displaystyle \sigma ^{2}=\int _{0}^{\infty }(t-{\bar {t}})^{2}\cdot E(t)\,dt}$

The third central moment indicates the skewness of the RTD and the fourth central moment indicates the kurtosis (the "peakedness").

One can also define an internal age distribution ${\displaystyle I(t)}$  that describes the reactor contents. This function has a similar definition as ${\displaystyle E(t)}$ : the fraction of fluid within the reactor with an age of ${\displaystyle t}$  is ${\displaystyle I(t)dt}$ . As shown by Danckwerts, the relation between ${\displaystyle E(t)}$  and ${\displaystyle I(t)}$  can be found from the mass balance:

${\displaystyle I(t)={\frac {1}{\tau }}\left(1-\int _{0}^{t}E(t)\ dt\right)\qquad E(t)=-\tau {\frac {dI(t)}{dt}}}$

## Determining the RTD experimentallyEdit

Residence time distributions are measured by introducing a non-reactive tracer into the system at the inlet. The concentration of the tracer is changed according to a known function and the response is found by measuring the concentration of the tracer at the outlet. The selected tracer should not modify the physical characteristics of the fluid (equal density, equal viscosity) and the introduction of the tracer should not modify the hydrodynamic conditions. The tracer chemical should have the following properties:

1. Non-reactive
2. Easily detectable
3. Properties similar to the reaction mixture
4. Completely soluble

In general, the change in tracer concentration will either be a pulse or a step. Other functions are possible, but they require more calculations to deconvolute the RTD curve, ${\displaystyle E(t)}$ .

### Pulse experimentsEdit

This method required the introduction of a very small volume of concentrated tracer at the inlet of the reactor, such that it approaches the Dirac delta function. Although an infinitely short injection cannot be produced, it can be made much smaller than the mean residence time of the vessel. If a mass of tracer, ${\displaystyle M}$ , is introduced into a vessel of volume ${\displaystyle V}$  and an expected residence time of ${\displaystyle \tau }$ , the resulting curve of ${\displaystyle C(t)}$  can be transformed into a dimensionless residence time distribution curve by the following relation:

${\displaystyle E(t)={\frac {C(t)}{\int _{0}^{\infty }C(t)\ dt}}}$

### Step experimentsEdit

The concentration of tracer in a step experiment at the reactor inlet changes abruptly from 0 to ${\displaystyle C_{0}}$ . The concentration of tracer at the outlet is measured and normalized to the concentration ${\displaystyle C_{0}}$  to obtain the non-dimensional curve ${\displaystyle F(t)}$  which goes from 0 to 1:

${\displaystyle F(t)={\frac {C(t)}{C_{0}}}}$ .

The step- and pulse-responses of a reactor are related by the following:

${\displaystyle F(t)=\int _{0}^{t}E(t')\,dt'\qquad E(t)={\frac {dF(t)}{dt}}}$

The value of the mean residence time and the variance can also be deduced from the function ${\displaystyle F(t)}$ :

${\displaystyle {\overline {t}}=\int _{0}^{\infty }t\cdot E(t)\,dt=\int _{0}^{1}t\,dF(t)=-\int _{1}^{0}t\,d(1-F(t))=\int _{0}^{\infty }(1-F(t))\,dt}$
${\displaystyle \sigma ^{2}=\int _{0}^{\infty }(t-{\overline {t}})^{2}\cdot E(t)\,dt=\int _{0}^{1}(t-{\overline {t}})^{2}\,dF(t)=\int _{0}^{1}t^{2}\,dF(t)-{\overline {t}}^{2}=2\int _{0}^{\infty }t(1-F(t))\,dt-{\overline {t}}^{2}}$

A step experiment is often easier to perform than a pulse experiment, but it tends to smooth over some of the details that a pulse response could show. It is easy to numerically integrate an experimental pulse response to obtain a very high-quality estimate of the step response, but the reverse is not the case because any noise in the concentration measurement will be amplified by numeric differentiation.

## RTDs of ideal and real reactorsEdit

The residence time distribution of a reactor can be used to compare its behaviour to that of two ideal reactor models: the plug-flow reactor and the continuous stirred-tank reactor (CSTR), or mixed-flow reactor. This characteristic is important in order to calculate the performance of a reaction with known kinetics.

### Plug flow reactorsEdit

In an ideal plug-flow reactor there is no axial mixing and the fluid elements leave in the same order they arrived. Therefore, fluid entering the reactor at time ${\displaystyle t}$  will exit the reactor at time ${\displaystyle t+\tau }$ , where ${\displaystyle \tau }$  is the residence time of the reactor. The residence time distribution function is therefore a Dirac delta function at ${\displaystyle \tau }$ .

${\displaystyle E(t)=\delta (t-\tau )\,}$

The variance of an ideal plug-flow reactor is zero.

The RTD of a real reactor deviate from that of an ideal reactor, depending on the hydrodynamics within the vessel. A non-zero variance indicates that there is some dispersion along the path of the fluid, which may be attributed to turbulence, a non-uniform velocity profile, or diffusion. If the mean of the ${\displaystyle E(t)}$  curve arrives earlier than the expected time ${\displaystyle \tau }$  it indicates that there is stagnant fluid within the vessel. If the RTD curve shows more than one main peak it may indicate channeling, parallel paths to the exit, or strong internal circulation.

### Continuous stirred-tank reactorsEdit

An ideal continuous stirred-tank reactor is based on the assumption that the flow at the inlet is completely and instantly mixed into the bulk of the reactor. The reactor and the outlet fluid have identical, homogeneous compositions at all times. An ideal CSTR has an exponential residence time distribution:

${\displaystyle E(t)={\frac {1}{\tau }}e^{-t/\tau }\,}$

In reality, it is impossible to obtain such rapid mixing, especially on industrial scales where reactor vessels may range between 1 and thousands of cubic meters, and hence the RTD of a real reactor will deviate from the ideal exponential decay. For example, there will be some finite delay before ${\displaystyle E(t)}$  reaches its maximum value and the length of the delay will reflect the rate of mass transfer within the reactor. Just as was noted for a plug-flow reactor, an early mean will indicate some stagnant fluid within the vessel, while the presence of multiple peaks could indicate channeling, parallel paths to the exit, or strong internal circulation. Short-circuiting fluid within the reactor would appear in an RTD curve as a small pulse of concentrated tracer that reaches the outlet shortly after injection.

### OceanographicEdit

In chemical oceanography, residence time (t) of an element is defined as the amount of an element in the ocean at steady state divided by the rate at which that element is added to the ocean:

${\displaystyle t=({\text{Mean concentration in ocean}})\times ({\text{Ocean volume}})/({\text{Input per year}})}$

where the ocean volume is (1.37×1021 L). The input sums all inputs to the ocean. For many elements, the major input is from rivers and the input per year is the mean river concentration times the continental runoff rate. If the concentration of an element is not changing, then the input and output of an element must be equal (steady state). The residence time can then be calculated using the estimated output, if that is known.

### BiochemicalEdit

Residence time is also used when studying bacteria, in which context it is given the symbol Г. It is the inverse of the eigenvalue derived from the mass balance method.

## ReferencesEdit

1. ^ H. SCOTT FOGLER, Elements of Chemical Reaction Engineering Fourth Edition.
1. 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.
2. 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.
3. Levenspiel, Octave (1999). Chemical Reaction Engineering (3rd ed.). John Wiley & Sons. ISBN 0-471-25424-X.
4. Nauman, E. Bruce (2004). "Residence Time Distributions". Handbook of Industrial Mixing: Science and Practice. Wiley Interscience. pp. 1–17. ISBN 0-471-26919-0.