Chemistry:Continuous stirred-tank reactor
The continuous stirred-tank reactor (CSTR), also known as vat- or backmix reactor, mixed flow reactor (MFR), or a continuous-flow stirred-tank reactor (CFSTR), is a common model for a chemical reactor in chemical engineering and environmental engineering. A CSTR often refers to a model used to estimate the key unit operation variables when using a continuous agitated-tank reactor to reach a specified output. The mathematical model works for all fluids: liquids, gases, and slurries.
The behavior of a CSTR is often approximated or modeled by that of an ideal CSTR, which assumes perfect mixing. In a perfectly mixed reactor, reagent is instantaneously and uniformly mixed throughout the reactor upon entry. Consequently, the output composition is identical to composition of the material inside the reactor, which is a function of residence time and reaction rate. The CSTR is the ideal limit of complete mixing in reactor design, which is the complete opposite of a plug flow reactor (PFR). In practice, no reactors behave ideally but instead fall somewhere in between the mixing limits of an ideal CSTR and PFR.
Ideal CSTR
Modeling
A continuous fluid flow containing non-conservative chemical reactant A enters an ideal CSTR of volume V.
Assumptions:
- perfect or ideal mixing
- steady state [math]\displaystyle{ \Bigl(\frac{dN_A}{dt} = 0\Bigr) }[/math], where NA is the number of moles of species A
- closed boundaries
- constant fluid density (valid for most liquids; valid for gases only if there is no net change in the number of moles or drastic temperature change)
- nth-order reaction (r = kCAn), where k is the reaction rate constant, CA is the concentration of species A, and n is the order of the reaction
- isothermal conditions, or constant temperature (k is constant)
- single, irreversible reaction (νA = −1)
- All reactant A is converted to products via chemical reaction
- NA = CA V
Integral mass balance on number of moles NA of species A in a reactor of volume V:[math]\displaystyle{ 1. [\text{Net accumulation of} ~A] = [A~\text{in}] - [A~\text{out}] + [\text{Net generation of} ~A] }[/math]
[math]\displaystyle{ 2. \frac{dN_A}{dt} = F_{Ao} - F_A + V \nu_A r_A }[/math][1]
where,
- FAo is the molar flow rate inlet of species A
- FA is the molar flow rate outlet of species A
- vA is the stoichiometric coefficient
- rA is the reaction rate
Applying the assumptions of steady state and νA = −1, Equation 2 simplifies to:
[math]\displaystyle{ 3. \ 0 = F_{Ao} - F_A - V r_A }[/math]
The molar flow rates of species A can then be rewritten in terms of the concentration of A and the fluid flow rate (Q):
[math]\displaystyle{ 4. \ 0 = QC_{Ao} - QC_A - V r_A }[/math][2]
Equation 4 can then be rearranged to isolate rA and simplified:
[math]\displaystyle{ 5. \ r_A = \frac{Q}{V} (C_{Ao} - C_A) }[/math][2]
[math]\displaystyle{ 6. \ r_A = \frac{1}{\tau} (C_{Ao} - C_A) }[/math]
where,
- [math]\displaystyle{ \tau }[/math] is the theoretical residence time ([math]\displaystyle{ \tau = \tfrac{V}{Q} }[/math])
- CAo is the inlet concentration of species A
- CA is the reactor/outlet concentration of species A
Residence time is the total amount of time a discrete quantity of reagent spends inside the reactor. For an ideal reactor, the theoretical residence time, [math]\displaystyle{ \tau }[/math], is always equal to the reactor volume divided by the fluid flow rate.[2] See the next section for a more in-depth discussion on the residence time distribution of a CSTR.
Depending on the order of the reaction, the reaction rate, rA, is generally dependent on the concentration of species A in the reactor and the rate constant. A key assumption when modeling a CSTR is that any reactant in the fluid is perfectly (i.e. uniformly) mixed in the reactor, implying that the concentration within the reactor is the same in the outlet stream.[3] The rate constant can be determined using a known empirical reaction rate that is adjusted for temperature using the Arrhenius temperature dependence.[2] Generally, as the temperature increases so does the rate at which the reaction occurs.
Equation 6 can be solved by integration after substituting the proper rate expression. The table below summarizes the outlet concentration of species A for an ideal CSTR. The values of the outlet concentration and residence time are major design criteria in the design of CSTRs for industrial applications.
| Reaction Order | CA |
|---|---|
| n=0 | [math]\displaystyle{ C_A = {C_{Ao}}-{k \tau } }[/math] |
| n=1 | [math]\displaystyle{ C_A = \frac {C_{Ao}}{1 + k \tau } }[/math][1] |
| n=2 | [math]\displaystyle{ C_A = \frac {-1 + \sqrt {1 + 4k\tau C_{Ao}}} {2k\tau} }[/math] |
| Other n | Numerical solution required |
Residence time distribution
An ideal CSTR will exhibit well-defined flow behavior that can be characterized by the reactor's residence time distribution, or exit age distribution.[4] Not all fluid particles will spend the same amount of time within the reactor. The exit age distribution (E(t)) defines the probability that a given fluid particle will spend time t in the reactor. Similarly, the cumulative age distribution (F(t)) gives the probability that a given fluid particle has an exit age less than time t.[3] One of the key takeaways from the exit age distribution is that a very small number of fluid particles will never exit the CSTR.[5] Depending on the application of the reactor, this may either be an asset or a drawback.
Non-ideal CSTR
While the ideal CSTR model is useful for predicting the fate of constituents during a chemical or biological process, CSTRs rarely exhibit ideal behavior in reality.[2] More commonly, the reactor hydraulics do not behave ideally or the system conditions do not obey the initial assumptions. Perfect mixing is a theoretical concept that is not achievable in practice.[6] For engineering purposes, however, if the residence time is 5–10 times the mixing time, the perfect mixing assumption generally holds true.
Non-ideal hydraulic behavior is commonly classified by either dead space or short-circuiting. These phenomena occur when some fluid spends less time in the reactor than the theoretical residence time, [math]\displaystyle{ \tau }[/math]. The presence of corners or baffles in a reactor often results in some dead space where the fluid is poorly mixed.[6] Similarly, a jet of fluid in the reactor can cause short-circuiting, in which a portion of the flow exits the reactor much quicker than the bulk fluid. If dead space or short-circuiting occur in a CSTR, the relevant chemical or biological reactions may not finish before the fluid exits the reactor.[2] Any deviation from ideal flow will result in a residence time distribution different from the ideal distribution, as seen at right.
Modeling non-ideal flow
Although ideal flow reactors are seldom found in practice, they are useful tools for modeling non-ideal flow reactors. Any flow regime can be achieved by modeling a reactor as a combination of ideal CSTRs and plug flow reactors (PFRs) either in series or in parallel.[6] For examples, an infinite series of ideal CSTRs is hydraulically equivalent to an ideal PFR.[2] Reactor models combining a number of CSTRs in series are often termed tanks-in-series (TIS) models.[7]
To model systems that do not obey the assumptions of constant temperature and a single reaction, additional dependent variables must be considered. If the system is considered to be in unsteady-state, a differential equation or a system of coupled differential equations must be solved. Deviations of the CSTR behavior can be considered by the dispersion model. CSTRs are known to be one of the systems which exhibit complex behavior such as steady-state multiplicity, limit cycles, and chaos.
Cascades of CSTRs
Cascades of CSTRs, also known as a series of CSTRs, are used to decrease the volume of a system.[8]
Minimizing Volume
As seen in the graph with one CSTR, where the inverse rate is plotted as a function of fractional conversion, the area in the box is equal to [math]\displaystyle{ \frac{V}{F_{Ao}} }[/math] where V is the total reactor volume and [math]\displaystyle{ F_{Ao} }[/math] is the molar flow rate of the feed. When the same process is applied to a cascade of CSTRs as seen in the graph with three CSTRs, the volume of each reactor is calculated from each inlet and outlet fractional conversion, therefore resulting in a decrease in total reactor volume. Optimum size is achieved when the area above the rectangles from the CSTRs in series that was previously covered by a single CSTR is maximized. For a first order reaction with two CSTRs, equal volumes should be used. As the number of ideal CSTRs (n) approaches infinity, the total reactor volume approaches that of an ideal PFR for the same reaction and fractional conversion.
Ideal Cascade of CSTRs
From the design equation of a single CSTR where [math]\displaystyle{ \tau = \frac{C_{Ao} - C_A}{-r_A} }[/math], we can determine that for a single CSTR in series that [math]\displaystyle{ \tau_i = \frac{C_{A(i-1)}-C_{Ai}}{-r_{Ai}} }[/math]
where [math]\displaystyle{ \tau }[/math] is the space time of the reactor, [math]\displaystyle{ C_{Ao} }[/math] is the feed concentration of A, [math]\displaystyle{ C_{A} }[/math] is the outlet concentration of A, and [math]\displaystyle{ -r_{A} }[/math] is the rate of reaction of A.
First order
For an isothermal first order, constant density reaction in a cascade of identical CSTRs operating at steady state
For one CSTR: [math]\displaystyle{ C_{A1} = \frac{C_{Ao}}{1 + k\tau} }[/math], where k is the rate constant and [math]\displaystyle{ C_{A1} }[/math] is the outlet concentration of A from the first CSTR
Two CSTRs: [math]\displaystyle{ C_{A1} = \frac{C_{Ao}}{1 + k\tau} }[/math] and [math]\displaystyle{ C_{A2} = \frac{C_{A1}}{1 + k\tau} }[/math]
Plugging in the first CSTR equation to the second: [math]\displaystyle{ C_{A2} = \frac{C_{Ao}}{(1 + k\tau)^2} }[/math]
Therefore for m identical CSTRs in series: [math]\displaystyle{ C_{Am} = \frac{C_{Ao}}{(1 + k\tau)^m} }[/math]
When the volumes of the individual CSTRs in series vary, the order of the CSTRs does not change the overall conversion for a first order reaction as long as the CSTRs are run at the same temperature.
Zeroth order
At steady state, the general equation for an isothermal zeroth order reaction at in a cascade of CSTRs is given by [math]\displaystyle{ C_{Am} = C_{Ao} - \sum_{i=1}^{m} k_i\tau_i }[/math]
When the cascade of CSTRs is isothermal with identical reactors, the concentration is given by [math]\displaystyle{ C_{Am} = C_{Ao} - mk_i\tau_i }[/math]
Second order
For an isothermal second order reaction at steady state in a cascade of CSTRs, the general design equation is [math]\displaystyle{ C_{Ai} = \frac{-1 + \sqrt{1 + 4k_i\tau_iC_{A(i-1)}}}{2k_i\tau_i} }[/math]
Non-ideal cascade of CSTRs
With non-ideal reactors, residence time distributions can be calculated. At the concentration at the jth reactor in series is given by
[math]\displaystyle{ \frac{C_j}{C_o} = 1 -e^{-\frac{nt}{\bar{t}}}[1 + \frac{nt}{\bar{t}} + \frac{1}{2!}(\frac{nt}{\bar{t}})^2 + ... + \frac{1}{(j-1)!}(\frac{nt}{\bar{t}})^{j-1}] }[/math]
where n is the total number of CSTRs in series, and [math]\displaystyle{ \bar{t} }[/math] is the average residence time of the cascade given by [math]\displaystyle{ \bar{t} = \frac{V}{Q} }[/math] where Q is the volumetric flow rate.
From this, the cumulative residence time distribution (F(t)) can be calculated as
[math]\displaystyle{ F(t) = \frac{C_n}{C_o} = 1 -e^{-\frac{nt}{\bar{t}}}[1 + \frac{nt}{\bar{t}} + \frac{1}{2!}(\frac{nt}{\bar{t}})^2 + ... + \frac{1}{(n-1)!}(\frac{nt}{\bar{t}})^{n-1}] }[/math]
As n → ∞, F(t) approaches the ideal PFR response. The variance associated with F(t) for a pulse stimulus into a cascade of CSTRs is [math]\displaystyle{ \sigma_{t}^2 = \frac{\bar{t}^2}{n} }[/math].
Cost
When determining the cost of a series of CSTRs, capital and operating costs must be taken into account. As seen above, an increase in the number of CSTRs in series will decrease the total reactor volume. Since cost scales with volume, capital costs are lowered by increasing the number of CSTRs. The largest decrease in cost, and therefore volume, occurs between a single CSTR and having two CSTRs in series. When considering operating cost, operating cost scales with the number of pumps and controls, construction, installation, and maintenance that accompany larger cascades. Therefore as the number of CSTRs increases, the operating cost increases. Therefore, there is a minimum cost associated with a cascade of CSTRs.
Zeroth order reactions
From a rearrangement of the equation given for identical isothermal CSTRs running a zeroth order reaction: [math]\displaystyle{ \tau = \frac{C_{Ao} - C_{Am}}{mk} }[/math], the volume of each individual CSTR will scale by [math]\displaystyle{ \frac{1}{m} }[/math]. Therefore the total reactor volume is independent of the number of CSTRs for a zeroth order reaction. Therefore, cost is not a function of the number of reactors for a zeroth order reaction and does not decrease as the number of CSTRs increases.
Selectivity of parallel reactions
When considering parallel reactions, utilizing a cascade of CSTRs can achieve greater selectivity for a desired product.
For a given parallel reaction [math]\ce{ A -> B }[/math] and [math]\ce{ A -> C }[/math] with constants [math]\displaystyle{ k_1 }[/math] and [math]\displaystyle{ k_2 }[/math] and rate equations [math]\displaystyle{ \frac{d[\ce B]}{dt}=k_1[\ce A]^{n_1} }[/math] and [math]\displaystyle{ \frac{d[\ce C]}{dt}=k_2[\ce A]^{n_2} }[/math], respectively, we can obtain a relationship between the two by dividing [math]\displaystyle{ \frac{d[\ce B]}{dt} }[/math] by [math]\displaystyle{ \frac{d[\ce C]}{dt} }[/math]. Therefore [math]\displaystyle{ \frac{d[\ce B]}{d[\ce C]} = \frac{k_1}{k_2}[\ce A]^{n_1 - n_2} }[/math]. In the case where [math]\displaystyle{ n_1 \gt n_2 }[/math] and B is the desired product, the cascade of CSTRs is favored with a fresh secondary feed of [math]\ce{ A }[/math] in order to maximize the concentration of [math]\ce{ A }[/math].
For a parallel reaction with two or more reactants such as [math]\ce{ A + D -> B }[/math] and [math]\ce{ A + D -> C }[/math] with constants [math]\displaystyle{ k_1 }[/math] and [math]\displaystyle{ k_2 }[/math] and rate equations [math]\displaystyle{ \frac{d[\ce B]}{dt}=k_1[\ce A]^{n_1}[\ce D]^{m_1} }[/math] and [math]\displaystyle{ \frac{d[\ce C]}{dt}=k_2[\ce A]^{n_2}[\ce D]^{m_12} }[/math], respectively, we can obtain a relationship between the two by dividing [math]\displaystyle{ \frac{d[\ce B]}{dt} }[/math] by [math]\displaystyle{ \frac{d[\ce C]}{dt} }[/math]. Therefore [math]\displaystyle{ \frac{d[\ce B]}{d[\ce C]} = \frac{k_1}{k_2}[\ce A]^{n_1 - n_2}[\ce D]^{m_1 - m_2} }[/math]. In the case where [math]\displaystyle{ n_1 \gt n_2 }[/math] and [math]\displaystyle{ m_1 \gt m_2 }[/math] and B is the desired product, a cascade of CSTRs with an inlet stream of high [math]\ce{ [A] }[/math] and [math]\ce{ [D] }[/math] is favored. In the case where [math]\displaystyle{ n_1 \gt n_2 }[/math] and [math]\displaystyle{ m_1 \lt m_2 }[/math] and B is the desired product, a cascade of CSTRs with a high concentration of [math]\ce{ A }[/math] in the feed and small secondary streams of [math]\ce{ D }[/math] is favored.[9]
Series reactions such as [math]\ce{ A -> B -> C }[/math] also have selectivity between [math]\ce{ B }[/math] and [math]\ce{ C }[/math] but CSTRs in general are typically not chosen when the desired product is [math]\ce{ B }[/math] as the back mixing from the CSTR favors [math]\ce{ C }[/math]. Typically a batch reactor or PFR is chosen for these reactions.
Employing CFD in CSTR
Computational Fluid Dynamics (CFD) can be beneficial in analyzing and optimizing the performance of a CSTR in several ways [10]. Below are a few examples:
CFD can simulate the fluid flow patterns inside the CSTR, providing insights into how well the reactants are mixed and distributed within the reactor. Efficient mixing is crucial for ensuring that the reactants are thoroughly and uniformly reacted, leading to better conversion rates and higher product quality.
CFD can help in determining the residence time distribution of reactants inside the reactor. Understanding the RTD is essential for assessing the reactor's performance, as it indicates how long the reactants stay inside the reactor and how effectively the reactions occur.
CFD can predict the temperature distribution within the CSTR and analyze heat transfer phenomena. This information is crucial for controlling exothermic or endothermic reactions and ensuring that the reactor operates under safe and optimal conditions.
CFD can be coupled with reaction kinetics models to predict the reaction rates within the reactor. By considering mass transport effects, CFD can provide valuable data on how different operating conditions impact the reaction rate and product formation.
CFD simulations allow engineers to optimize the design and operating parameters of the CSTR. By testing different configurations virtually, they can identify the most efficient setup, reducing the need for costly and time-consuming physical prototypes. Additionally, CFD can aid in scale-up studies, providing insights into how the reactor's performance might change at different scales.
CFD can help in troubleshooting existing CSTR systems. If there are issues like dead zones, inadequate mixing or temperature variations, CFD simulations can pinpoint the problems and suggest modifications to improve reactor performance.
Applications
CSTRs facilitate rapid dilution of reagents through mixing. Therefore, for non-zero-order reactions, the low concentration of reagent in the reactor means a CSTR will be less efficient at removing the reagent compared to a PFR with the same residence time.[3] Therefore, CSTRs are typically larger than PFRs, which may be a challenge in applications where space is limited. However, one of the added benefits of dilution in CSTRs is the ability to neutralize shocks to the system. As opposed to PFRs, the performance of CSTRs is less susceptible to changes in the influent composition, which makes it ideal for a variety of industrial applications:
Environmental engineering
- Activated sludge process for wastewater treatment[2]
- Lagoon treatment systems for natural wastewater treatment[2]
- Anaerobic digesters for the stabilization of wastewater biosolids[11]
- Treatment wetlands for wastewater and storm water runoff[12]
Chemical engineering
- Loop reactor for production of pharmaceuticals[13]
- Fermentation[13]
- Biogas production
See also
Notes
References
- ↑ 1.0 1.1 Schmidt, Lanny D. (1998). The Engineering of Chemical Reactions. New York: Oxford University Press. ISBN 0-19-510588-5.
- ↑ 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 Metcalf & Eddy (2013-09-03). Wastewater engineering : treatment and resource recovery. Tchobanoglous, George,, Stensel, H. David,, Tsuchihashi, Ryujiro,, Burton, Franklin L. (Franklin Louis), 1927-, Abu-Orf, Mohammad,, Bowden, Gregory (Fifth ed.). New York, NY. ISBN 978-0-07-340118-8. OCLC 858915999.
- ↑ 3.0 3.1 3.2 Benjamin, Mark M. (2013-06-13). Water quality engineering : physical/chemical treatment processes. Lawler, Desmond F.. Hoboken, New Jersey. ISBN 978-1-118-63227-7. OCLC 856567226.
- ↑ Bolin, Bert; Rodhe, Henning (January 1973). "A note on the concepts of age distribution and transit time in natural reservoirs". Tellus 25 (1): 58–62. doi:10.3402/tellusa.v25i1.9644. ISSN 0040-2826.
- ↑ Monsen, Nancy E.; Cloern, James E.; Lucas, Lisa V.; Monismith, Stephen G. (September 2002). "A comment on the use of flushing time, residence time, and age as transport time scales" (in en). Limnology and Oceanography 47 (5): 1545–1553. doi:10.4319/lo.2002.47.5.1545. Bibcode: 2002LimOc..47.1545M.
- ↑ 6.0 6.1 6.2 Davis, Mark E. (2003). Fundamentals of chemical reaction engineering. Davis, Robert J. (International ed.). Boston: McGraw-Hill. ISBN 978-1-62870-437-2. OCLC 880604539.
- ↑ Stokes, R. L.; Nauman, E. Bruce (1970). "Residence Time Distribution Functions for Stirred Tanks in Series". Canadian Journal of Chemical Engineering 48 (6): 723-725.
- ↑ Hill, Charles G.; Root, Thatcher W. (2014). Introduction to Chemical Engineering Kinetics and Reactor Design, Second Edition. Hoboken, new Jersey: Wiley. pp. 241-253, 349-358. ISBN 9781118368251.
- ↑ Levenspiel, Octave (1998). Chemical Reaction Engineering, 3rd Edition. Wiley. ISBN 978-0-471-25424-9.
- ↑ "Employing CFD in CSTR". https://www.jimmylea.com/cfd-consultants.html. Retrieved 17 July 2023.
- ↑ Hurtado, F.J.; Kaiser, A.S.; Zamora, B. (March 2015). "Fluid dynamic analysis of a continuous stirred tank reactor for technical optimization of wastewater digestion". Water Research 71: 282–293. doi:10.1016/j.watres.2014.11.053. ISSN 0043-1354. PMID 25635665.
- ↑ Kadlec, Robert H.; Wallace, Scott D. (2009). Treatment Wetlands (second ed.). Boca Raton, FL, USA: CRC Press. p. 181. ISBN 978-1-56670-526-4.
- ↑ 13.0 13.1 "Visual Encyclopedia of Chemical Engineering". http://encyclopedia.che.engin.umich.edu/Pages/Reactors/CSTR/CSTR.html.
 |  |
