|
|
| Line 1: |
Line 1: |
| 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 engineering|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.
| | Hi there. Let me start by introducing the writer, her title is Myrtle Cleary. For years he's been residing in North Dakota and his family enjoys it. It's not a typical factor but what she likes performing is foundation leaping and now she is trying to make cash with it. Bookkeeping is what I do.<br><br>My web-site: [http://Chorokdeul.Co.kr/index.php?document_srl=454418&mid=customer21 at home std testing] |
| | |
| The concept was first proposed by MacMullin and Weber in 1935, but was not used extensively until [[Peter Victor Danckwerts|P.V. Danckwerts]] analyzed a number of important RTDs in 1953.
| |
| | |
| ==Theory==
| |
| [[Image:Typical CSTR RTD curve.svg|thumb|right|An RTD curve for a reasonably well-mixed reactor]]
| |
| | |
| The theory of residence time distributions generally begins with three assumptions:
| |
| # the reactor is at steady-state,
| |
| # transports at the inlet and the outlet takes place only by [[advection]], and
| |
| # the fluid 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 time]]s is represented by an exit age distribution, <math>E(t)</math>. The function <math>E(t)</math> has the units of time<sup>-1</sup> and is defined such that
| |
| : <math>\int_{0}^\infty E(t)\, dt = 1</math>.
| |
| | |
| The fraction of the fluid that spends a given duration, <math>t</math> inside the reactor is given by the value of <math>E(t) dt</math>.
| |
| | |
| The fraction of the fluid that leaves the reactor with an age less than <math>t_1</math> is
| |
| : <math>\int_{0}^{t_1} E(t)\, dt</math>.
| |
| | |
| The fraction of the fluid that leaves the reactor with an age greater than <math>t_1</math> is
| |
| : <math>\int_{t_1}^\infty E(t)\, dt = 1-\int_{0}^{t_1} E(t)\, dt</math>.
| |
| | |
| The average residence time is given by the first [[moment (mathematics)|moment]] of the age distribution:
| |
| : <math> \bar t = \int_{0}^\infty t \cdot E(t)\, dt</math>.
| |
| | |
| If there are no dead, or stagnant, zones within the reactor then <math>\bar t</math> will be equal to <math>\tau</math>, the [[residence time]] calculated from the total reactor volume and the volumetric flow rate of the fluid:
| |
| | |
| : <math>\tau = \frac {V}{v}</math>.
| |
| | |
| The higher order [[central moment]]s can provide significant information about the behavior of the function <math>E(t)</math>. For example, the second central moment indicates the variance (<math>\sigma^2</math>), the degree of dispersion around the mean.
| |
| : <math>\sigma^2 = \int_{0}^\infty (t- \bar t)^2 \cdot E(t)\, dt</math>
| |
| | |
| 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 <math>I(t)</math> that describes the reactor contents. This function has a similar definition as <math>E(t)</math>: the fraction of fluid within the reactor with an age of <math>t</math> is <math>I(t) dt</math>. As shown by Danckwerts, the relation between <math>E(t)</math> and <math>I(t)</math> can be found from the mass balance:
| |
| : <math>I(t) = \frac{1}{\tau}\left(1-\int_0^t E(t)\ dt\right) \qquad E(t) = -\tau \frac{dI(t)}{dt}</math>
| |
| | |
| ==Determining the RTD Experimentally==
| |
| 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.
| |
| | |
| 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 [[deconvolution|deconvolute]] the RTD curve, <math>E(t)</math>.
| |
| | |
| ===Pulse Experiments===
| |
| 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, <math>M</math>, is introduced into a vessel of volume <math>V</math> and an expected residence
| |
| time of <math>\tau</math>, the resulting curve of <math>C(t)</math> can be transformed into a dimensionless residence time distribution curve by the following relation:
| |
| | |
| : <math>E(t) = \frac{C(t)}{\int_0^\infty C(t)\ dt}</math>
| |
| | |
| ===Step Experiments===
| |
| In a step experiment the concentration of tracer at the reactor inlet changes abruptly from 0 to <math>C_0</math>. The concentration of tracer at the outlet is measured and normalized to the concentration <math>C_0</math> to obtain the non-dimensional curve <math>F(t)</math> which goes from 0 to 1:
| |
| : <math>F(t) = \frac {C(t)}{C_{0}}</math>.
| |
| | |
| The step- and pulse-responses of a reactor are related by the following:
| |
| : <math>F(t) = \int_{0}^t E(t)\, dt \qquad E(t) = \frac {dF(t)}{dt}</math>
| |
| | |
| The value of the mean residence time and the variance can also be deduced from the function <math>F(t)</math>:
| |
| : <math>\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</math>
| |
| : <math>\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</math>
| |
| | |
| 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 reactors==
| |
| The residence time distribution of a reactor can be used to compare its behavior to that of two ideal reactor models: the [[plug flow reactor model|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 Reactors===
| |
| In an ideal [[Plug flow reactor model|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 <math>t</math> will exit the reactor at time <math>t+\tau</math>, where <math>\tau</math> is the residence time of the reactor. The residence time distribution function is therefore a [[dirac delta function]] at <math>\tau</math>.
| |
| | |
| : <math>E(t) = \delta(t-\tau)\,</math> | |
| | |
| 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 <math>E(t)</math> curve arrives earlier than the expected time <math>\tau</math> 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.
| |
| | |
| ===Mixed Flow Reactors===
| |
| 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:
| |
| : <math>E(t) = \frac{1}{\tau} e^{-t/\tau}\,</math> | |
| | |
| In reality, it is impossible to obtain such rapid mixing, especially on industrial scales where reactor vessels may range between 1 and several tens 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 <math>E(t)</math> 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.
| |
| | |
| ===Oceanographic===
| |
| | |
| 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:
| |
| | |
| t = (Mean Concentration in Ocean) × (Ocean Volume) / (Input per year)
| |
| where the ocean volume is (1.37×10^21 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 × 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.
| |
| | |
| ===Biochemical===
| |
| | |
| Residence time not only relates to hydraulic residence time but bacterial residence time as well. It has a symbol Г (tow). It is the inverse of the eigenvalue derived form the mass balance method.
| |
| | |
| ==References==
| |
| #{{cite journal
| |
| | author = R.B. MacMullin and M. Weber
| |
| | title = The theory of short-circuiting in continuous-flow mixing vessels in series and kinetics of chemical reactions in such systems
| |
| | journal = Transactions of American Institute of Chemical Engineers
| |
| | year = 1935
| |
| | volume = 31
| |
| | issue = 2
| |
| | pages = 409–458 }}
| |
| #{{cite journal
| |
| | author = P.V. Danckwerts
| |
| | title = Continuous flow systems. Distribution of residence times.
| |
| | journal = Chemical Engineering Science
| |
| | year = 1953
| |
| | volume = 2
| |
| | pages = 1–13
| |
| | doi = 10.1016/0009-2509(53)80001-1 }}
| |
| #{{cite book
| |
| | last = Levenspiel
| |
| | first = Octave
| |
| | title = Chemical Reaction Engineering
| |
| | edition = 3rd
| |
| | publisher = John Wiley & Sons
| |
| | year = 1999
| |
| | isbn = 0-471-25424-X }}
| |
| #{{cite book
| |
| | last = Nauman
| |
| | first = E. Bruce
| |
| | chapter = Residence Time Distributions
| |
| | title = Handbook of Industrial Mixing: Science and Practice
| |
| | publisher = Wiley Interscience
| |
| | year = 2004
| |
| | pages = 1–17
| |
| | isbn = 0-471-26919-0}}
| |
| | |
| [[Category:Chemical engineering]]
| |