Reactors

Title: Reactors

Authors:

Sean Cabaniss, David Park, Maxim Slivinsky, and Julianne Wagoner ^{[Winter 2014]}

Neil Dalvie, KC Anderson ^{[Winter 2016]}

Steward: David Chen, Fengqi You

Date Presented: February 4, 2014

Introduction

The center of any chemical process is the reactor, where chemical reactions are carried out to transform feeds into products. Reactor design is a vital step in the overall design of a process. It is important to ensure that the equipment specified will be capable of achieving the desired yields and selectivity.

Ideal Reactors

Batch Reactors

In a batch reactor, the reagents are added together and allowed to react for a given amount of time. The compositions change with time, but there is no flow through the process. Additional reagents may be added as the reaction proceeds, and changes in temperature may also be made. Products are removed from the reactor after the reaction has proceeded to completion.

Batch processes are suitable for small-scale production (less than 1,000,000 lb/yr) and for processes where several different products or grades are to be produced in the same equipment (Douglas, 1988). When production volumes are relatively small and/or the chemistry is relatively complex, batch processing provides an important means of quality control.

Plug Flow Reactor (PFR)

A PFR with tubular geometry has perfect radial mixing but no axial mixing. All materials hav the same residence time, τ, and experience the same temperature and concentration profiles along the reactor. Equation for PFR is given by:

$dM=\Re dV$

where M = molar flow rate, dV is the incremental volume, and $\Re$ is the rate of reaction per unit volume.

This equation can be integrated along the length of the reactor to yield relationships between reactor resident time and concentration or conversion.

Continuously Stirred Tank Reactor (CSTR)

The stirred tank reactor models a large scale conventional laboratory flask and can be considered to be the basic chemical reactor. In a CSTR, shown in Figure 1, there is no spatial variation- the entire vessel contents are at the same temperature, pressure, and concentration. Therefore the fluid leaving the reactor is at the same temperature and concentration as the fluid inside the reactor.

The material balance across the CSTR is given by:

$M_{\text{in}}-M_{\text{out}}=\Re V$

Some of the material the enters the reactor can leave immediately, while some leaves much later, so there is a broad distribution in residence time as shown in Figure 1.

Figure 1. Continuously Stirred Tank Reactor (Towler and Sinnott, 2013)

More information on stirred tanks can be found in the Mixing section.

Fluidized Bed Reactor

The fluidized bed reactor is a reactor in which the fluid from the bottom of the reactor keeps the solid catalysts suspended inside the reactor. This reactor has an advantage over traditional packed bed reactor because it has better heat and mass transfer. In industry, the fluidized bed reactor is most prominently used to produce of aniline, through the hydrogenation of nitrobenzene.

General Reactor Design

The design of the reactor should not be carried out separately from the overall process design due to the significant impact on capital and operating costs on other parts of the process (Towler and Sinnott, 2013).

Step 1: Collect Required Data

Out of all process equipment, reactor design requires the most process input data: reaction enthalpies, phase-equilibrium constants, heat and mass transfer coefficients, as well as reaction rate constants. All of the aforementioned parameters can be estimated using simulation models or literature correlations except for reaction rate constant constants, which need to be determined experimentally (Towler and Sinnott, 2013).

Enthalpy of Reaction

The heat given out in a chemical reaction is based on the enthalpies of the component chemical reactions, which are given for standard temperature and pressure (1 atm, 25 C). Values for standard heats of reaction can be found tabulated in literature, or can be calculated from heats of formation or combustion. Care must be taken to quote the basis for the heat of reaction and the states of reactants and products.

The following equation is used to convert enthalpies from standard conditions to the process conditions:

$\Delta H_{\text{r,P,T}}=\Delta H_{\text{r}}^{\circ }+\int _{1}^{P}\left[\left({\frac {\partial H_{prod.}}{\partial P}}\right)_{T}-\left({\frac {\partial H_{react.}}{\partial P}}\right)_{T}\right]dP+\int _{298}^{T}\left[\left({\frac {\partial H_{prod.}}{\partial T}}\right)_{P}-\left({\frac {\partial H_{react.}}{\partial T}}\right)_{P}\right]dT$

If the effect from pressure is not significant and only Temperature needs to be accounted for, the following equation should be used:

$\Delta H_{\text{r,T}}=\Delta H_{\text{r}}^{\circ }+\Delta H_{\text{prod.}}+\Delta H_{\text{react.}}$

Equilibrium Constant and Gibbs Free Energy

$\Delta G=-\mathbf {R} T\ln K$

Where $\Delta G$ is the change in Gibbs free energy from the reaction at temperature $T$, $\mathbf {R}$ is the ideal gas constant, and $K$ is the reaction equilibrium constant, given by:

$K=\prod _{i=1}^{n}{a_{i}}^{\alpha _{i}}$

where $a_{i}$ is the activity of component i, $\alpha _{i}$ is the stoichiometric coefficient of component $i$, and $n$ is the total number of components.

Equilibrium constants can be found in the literature and are useful for evaluating the rates of forward and reverse reactions. Care must be taken to the experimental design used for the literature equilibrium constants to make sure they are consistent with the conditions of the actual process reactor. For more complicated reactions consisting of several sequential or simultaneous reactions, the equilibrium is found by minimizing the Gibbs free energy (Towler and Sinnott, 2013). Commercial process simulation programs use the Gibbs reactor model in this way.

Reaction Mechanisms, Rate Equations, and Rate Constants

In most cases the main process reaction rate equations and rate constants cannot be predicted from first principles and must be approximated (Towler and Sinnott, 2013). This is due to the following:

• Use of heterogeneous catalysis or enzymes which lead to Langmuir-Hinshelwood-Hougen-Watson or Michaelis-Menten kinetics
• Mass transfer between vapor and liquid or two liquid phases
• Multistep mechanisms whose rate expressions do not follow overall reaction stoichiometry
• Competing side reactions

As a result the main process reaction is usually approximated as first- or second-order over a narrow range of process conditions (temperature, pressure, species concentrations) to estimate the residence time required for a target conversion. Rate equations are always a fit for experimental data and should thus be used for interpolation within the data. It is important to collect more data when extrapolating, especially for exothermic reactions which have the potential for runaway (Towler and Sinnott, 2013).

Heat and Mass Transfer Properties

Heat Transfer

The design of internal heating or cooling devices can be found in Heat Transfer Equipment. Correlations for tube-side heat-transfer coefficients for catalyst-packed tubes of a heat exchanger are given below:

For heating: ${{h_{i}d_{t}} \over \lambda _{f}}=.813{\left({\frac {\rho _{f}ud_{p}}{\mu }}\right)}^{.9}e^{-6d_{p}/d_{t}}$

and for cooling: ${{h_{i}d_{t}} \over \lambda _{f}}=3.50{\left({\frac {\rho _{f}ud_{p}}{\mu }}\right)}^{.7}e^{-4.6d_{p}/d_{t}}$

where $h_{i}$ is the tube-side heat transfer coefficient for a packed tube, $d_{t}$ is the tube diameter, $\lambda _{f}$ is the fluid thermal conductivity, $\rho _{f}$ is the fluid density, $u$ is the superficial velocity, $d_{p}$ is the effective particle diameter, and $\mu$ is the fluid viscosity.

Diffusion Coefficients

Diffusion coefficients are necessary when mass transfer can limit the rate of reaction, such as in catalytic reactions or reactions involving mass transfer processes such as gas absorption, distillation, and liquid-liquid extraction.

The diffusivity for gases can be estimated by the following correlation (Fuller, Schettler, Giddings):

$D_{v}={\frac {1.013\times 10^{-7}T^{1.75}{\left({\frac {1}{M_{a}}}+{\frac {1}{M_{b}}}\right)}^{1/2}}{P{\left[{\left(\sum _{a}v_{i}\right)}^{1/3}+{\left(\sum _{b}v_{i}\right)}^{1/3}\right]}^{2}}}$

where $D_{v}$ is the diffusivity, $T$ is temperature, $M_{a},M_{b}$ are the molecular masses of components $a$ and $b$, $P$ is the total pressure, and $\sum _{a}v_{i},\sum _{b}v_{i}$ are the summation of special diffusion volume coefficients for components $a$ and $b$, given in the table below:

(volume coefficient table from towler)

Wilke and Chang developed a correlation for estimating the diffusivity of components in the liquid phase:

$D_{L}={\frac {1.173\times 10^{-13}{(\phi M_{w})}^{1/2}T}{\mu V_{m}^{.6}}}$

where $D_{L}$ is the liquid diffusivity, $\phi$ is an association factor for the solvent, $M_{w}$ is the molecular mass of the solvent, $\mu$ is the solvent viscosity, $T$ is the temperature, and $V_{m}$ is the molar volume of the solute at its boiling point. This correlation holds for organic compounds in water but not for water in organic solvents.

Mass Transfer

For multiphase reactors it is necessary to estimate the mass transfer coefficient.

The equation of Gupta and Thodos predicts the mass transfer coefficient for a packed bed of particles:

${\frac {kd_{p}}{D}}=2.06{\frac {1}{\epsilon }}{Re}^{.425}{Sc}^{.33}$

where $k$ is the mass transfer coefficient, $d_{p}$ is the particle diameter, $D$ is the diffusivity, $Re$ is the Reynolds number calculated using the superficial velocity through the bed, $Sc$ is the Schmidt number, and $\epsilon$ is the bed void fraction.

Mass transfer between vapor and liquid in an agitated vessel can be described by the Van't Riet equations:

For air-water: $k_{L}a=0.026{\left({\frac {P_{a}}{V_{liq}}}\right)}^{.4}Q^{1/2}$

and for air-water-electrolyte: $k_{L}a=0.002{\left({\frac {P_{a}}{V_{liq}}}\right)}^{.7}Q^{.2}$

where $k_{L}$ is the mass transfer coefficient, $a$ is the interfacial area per unit volume, $Q$ is the gas volumetric flow rate, $V_{liq}$ is the liquid volume, and $P_{a}$ is the agitator power input.

Fair's method for calculating the mass transfer coefficient for low viscosity systems is given by:

${\frac {{(k_{L}a)}_{system}}{{(k_{L}a)}_{air-water}}}={\left({\frac {D_{L,system}}{D_{L,air-water}}}\right)}^{1/2}$

where $D_{L}$ is the liquid phase diffusivity.

Mass transfer correlations for vapor-liquid systems should be used with caution when there are surfactants (Towler and Sinnott, 2013).

Step 2: Select Reaction Conditions

A major determining factor in reactor type selection is the choice of operating conditions. Optimal process operation usually involves optimizing process yield and not necessarily reactor yield. Based on the preliminary economical analysis a target range of yields and selectivities can be chosen. The final reaction conditions must be verified experimentally to ensure target yields and selectitivities are realized (Towler and Sinnott, 2013).

Chemical or Biochemical Reaction

If the desired product is to be produced by a biochemical reaction the chosen conditions must maintain the viability of the biological agent (e.g. microorganisms or enzymes). Proteins denature outside of their specific temperature and pH ranges, while living organisms require specific concentrations of oxygen and other solutes to survive and cannot withstand high shear rates. See bioreactors for further information on their design.

Catalyst

A catalyst is used to increase the reaction rate by lowering the activation energy without being consumed in the reaction. The use of catalyst imposes operating condition constraints as the catalyst must maintain activity for a period of time between catalyst regenerations. Catalyst deactivation can be accelerated by high temperatures as well as contaminants in the feed or recycle streams.

Temperature

Increasing the reaction temperature will increase the reaction rate, diffusivities, and mass-transfer rates. Temperature also affects the equilibrium constant: higher temperature increases equilibrium constant for endothermic reactions and decreases it for exothermic reaction- see the figure below.

Figure 2. Effect of temperature on equilibrium constant (Towler and Sinnott, 2013)

Increased reaction temperature will reduce the cost of reactor design except for the following scenarios/considerations (Towler and Sinnott, 2013):

• Biochemical reactions where living organisms could die at high temperatures
• Presence of organic compounds that undergo thermal degradation
• Unwanted side reactions that accelerate with higher temperature, such as polymerization or auto-oxidation
• Oxidation reactions where selectivity decreases at higher temperatures as product oxidation tends to increase
• Exothermic reactions as it is more difficult to control the temperature and there is risk of reaction run away
• Construction cost of the reactor can become prohibitive at extremely high temperatures

Pressure

The main consideration when choosing the reactor pressure is to maintain the reaction at the desired phase for the selected temperature. The pressure can also be chosen to allow for vaporization of a component, making separation of a product easier, shifting the reaction equilibrium, or removing heat from the reactor. Increasing pressure for reactions that take place in the gas phase increases reactant activity and thus the reaction rate. Reactor yields follow Le Chatelier's principle: for reactions that increase number of moles lower pressure will increase equilibrium conversion, for reactions that decrease number of moles lower pressure will decrease equilibrium conversion. Increasing the pressure in gas-liquid reactions increases the solubility of the gas in the liquid which increases the reaction rate.

Reaction Phase

Reactions are usually carried out in liquid or gas phases as fluids are easier to handle, heat and cool, and transport than solids. For reagents or products in the solid phase a suspension in liquid or gas is usually used. The phase of the reaction is usually determined by reactor temperature and pressure. Liquid-phase operation is usually preferred due to the highest concentrations and greatest compactness. However, at temperatures above the critical temperature there cannot be a liquid phase. The pressure can sometimes be adjusted to keep all reagents in the liquid phase, however when this is not possible a multiphase reactor will be necessary. If mass transfer limitations become too significant it can be beneficial to reduce the pressure such that the reaction temperature is above the dew point and the reaction is carried out in the vapor phase.

Solvent

Solvents are used for liquid-phase reactions and can be used for the following:

• Dilution of feed to improve selectivity
• Increasing solubility of gas-phase components
• Dissolving solids in the reacting phase
• Increasing thermal mass which lowers temperature change per unit volume from reaction
• Improving miscibility of mutually insoluble components

Solvents should be inert in the main reaction and should not react with products or feed contaminants. Solvents should also be inexpensive and easily separated from the reaction products. Some widely used process solvents and their properties are given in the table below:

Commonly used process solvents (Towler and Sinnott, 2013)

Concentrations

Higher concentrations of feed can lead to higher reaction rate, however for exothermic reactions high feed concentrations should be avoided. Feed compounds are usually not supplied in stoichiometric ratio as using a higher concentration of one feed can lead to increased selectivity towards the desired product.

Understanding the effect of feed contaminants and by-products is essential to reactor design; they can play significant roles in reactor selectivity and performance. When recycling attention must be paid to by-products; those formed through reversible reactions can be recycled leading to improved overall selectivity. Feed contaminants generally pose a greater issue than by-products due to their ability to poison catalysts or kill biological organisms. If a feed contaminant is particularly detrimental to the reactor performance it should be removed upstream of the reactor.

Inert compounds will usually increase reactor cost due to the larger volume required, as well as increase downstream separation costs; they can still be advantageous for the following circumstances:

• Inerts in gas-phase reactions reduce partial pressure of reagents which can increase equilibrium conversion in reactions that lead to an increase in number of moles
• Feed compound reacting with itself or products can be reduced by dilution using inerts
• Inerts can allow operation outside of the flammability envelope
• Reaction solutions can be buffered to control pH

Step 3: Determine Materials of Construction

A preliminary analysis of the materials of construction for the reactor can be conducted after the reaction conditions have been specified. Particularly important in this analysis are the temperatures and pressures the process will run at. At extreme conditions, costly alloys may need to be used. In addition, the designer must ensure that process streams will not react with materials used in process equipment.

Step 4: Determine Rate-Limiting Step and Critical Sizing Parameters

The key parameters that determine the extent of reaction must be identified by carrying out an experiment plan with a broad range of conditions. In general, the rate of reaction is usually limited by the following fundamental processes. The first three have been discussed in previous sections. Mixing will be developed in more detail in its own section.

• Intrinsic kinetics: There will usually be one slowest step that governs the overall rate.
• Mass-transfer rate: In multiphase reactions and processes that use porous heterogeneous catalysis, mass transfer can be particularly important. Often, careful experimentation will be needed to separate the effects of mass transfer and the rate of reaction to determine which is the rate-limiting step.
• Heat-transfer rate: The rate of heat addition can become the governing parameter for endothermic reactions. Heat-transfer devices such as heat exchangers or fired heaters may need to be used.
• Mixing: The time taken to mix the reagents can be the limiting step for very fast reactions.

Once rate data have been collected, the designer can fit a suitable model of reaction kinetics. Next, a critical sizing parameter can be specified for the reactor. This will usually be one of the parameters given in Figure 1.

Figure 1. Reactor Sizing Parameters (Towler and Sinnott, 2013)

Step 5: Preliminary Sizing, Layout, and Costing of Reactor

The designer can estimate the reactor and catalyst volume from the sizing parameter. This calculation will yield a value for the active reacting volume necessary. Clearly, the actual reactor will need additional space. The geometry of the reactor will depend on the desired flow pattern and mixing requirements (Towler and Sinnott, 2013). The cost of most reactors can be estimated by determining the cost of a pressure vessel with the same dimensions and adding in the cost of the internals (Towler and Sinnott, 2013).

Step 6: Estimate Reactor Performance

At this point in the design process, it is important to verify that the proposed reactor will achieve the target conversions and selectivities. A combination of experimental methods, such as pilot plants, and computer simulations can be used to predict the full-scale reactor performance.

Step 7: Optimize the Design

The reactor is typically a relatively small fraction of the total capital cost (Towler and Sinnott, 2013), so minimal time should be devoted to optimization to reduce the reactor cost. However, if the target conversion, yields, and selectivities are not met, the process economics could be significantly impacted. Therefore, steps 2 to 6 should be repeated at least until the minimum specifications are met (Towler and Sinnott, 2013).

How to Design a Fluidized Bed Reactor

In order to manually design a fluidized bed reactor which would yield a useful information, myriads of parameters such as porosity at minimum fluidization, minimum fluidization velocity, bubble size, minimum fluidization velocity, velocity of bubble rise, fraction of the total bed occupied by bubbles, fraction of the bed consisting of wakes, volume of the catalyst in the bubbles, clouds and emulsion, mass transfer coefficient between bubble and cloud and additional information must be obtained experimentally. For most purposes, such manual calculation is not only tedious but also has a great chance of yielding an incorrect result. Instead of manual calculation, it is efficient to take advantage of fluidized bed reactor modeling capability offered in Aspen Plus.

Mixing in Industrial Reactors

Mixing plays an important role in many processing stages, including reactor performance. It is critical to select the appropriate method of mixing in order to ensure the process produces the desired process yields, product purity, and cost effectiveness.

Correlations such as the Reynolds number can be used to determine the extent of mixing and correlate power consumption and heat transfer to the reactor shell (Towler, 2012). In some cases, simple correlations may not be adequate:

• If dead zones cannot be tolerated for reasons of product purity, safety, etc.
• If reactor internals are complex
• If reaction selectivity is very sensitive to mixing

In these cases, it is usually necessary to carry out a more sophisticated analysis of mixing:

• Use computational fluid dynamics to model the reactor
• Use physical modeling (“cold flow”) experiments
• Use tomography methods to look at performance of real reactor

Gas Mixing

Gases mix easily because of their low viscosities. The mixing given by turbulent flow in a length of pipe is usually sufficient for most purposes (Towler and Sinnott, 2013). Orifices, vanes, and baffles can be used to increase turbulence.

Liquid Mixing

• Inline Mixing Inline mixers can be used for the continuous mixing of low-viscosity fluids. One inexpensive method involves the use of static devices that promote turbulent mixing in pipelines. Some typical designs are shown in Figures 2(a), (b), and (c).

Figure 2. Inline mixers: (a) tee; (b) injection; (c) annular (Towler and Sinnott, 2013)

When mixing low viscosity fluids (<50 mNs/m²) with similar densities and flow rates, a simple mixing tee, Figure 2(a), followed by a length of pipe equal to 10 to 20 pipe diameters, is suitable (Towler and Sinnott, 2013).

When one flow is much lower than the other, an injection mixer, Figure 2(b&c), should be used. A satisfactory blend will be achieved in about 80 pipe diameters (Towler and Sinnott, 2013). Baffles or other flow restrictions can be used to reduce the mixing length required. These mixers work by introducing one fluid into the flowing stream of the other through a concentric pipe or an annular array of jets (Towler and Sinnott, 2013).

• Stirred Tanks Stirred tanks were discussed in the Ideal Reactors section. Mixing is conducted by an impeller mounted on a shaft driven by a motor. The reactor usually contains baffles or other internals to induce turbulence and prevent the contents from swirling and creating a vortex. Typically, baffles are 1/10 of diameter and located 1/20 of diameter from wall (Towler, 2012). A typical arrangement of agitator and baffles in a stirred tank, and the flow pattern generated, is shown in Figure 3. Mixing occurs through the bulk flow of the liquid and by the motion of the turbulent eddies created by the agitator. Bulk flow is the predominant mixing mechanism required for the blending of miscible liquids and for solids suspension. Turbulent mixing is important in operations involving mass and heat transfer, which can be considered as shear-controlled processes (Towler and Sinnott, 2013).

Figure 3. Agitator arrangements and flow patterns (Towler and Sinnott, 2013)

At high Reynolds numbers (low viscosity), one of the three basic types of impeller shown in Figure 4 should be used. For processes controlled by turbulent mixing, the flat-bladed (Rushton) turbines are appropriate. For bulk mixing, the propeller and pitched-bladed turbines are appropriate (Towler and Sinnott, 2013).

Figure 4. Basic impeller types (Towler and Sinnott, 2013)

For more viscous fluids, paddle, anchor, and helical ribbon agitators (Figures 5(a), (b), and (c)), are used (Towler and Sinnott, 2013). The selection chart given in Figure 6 can be used to make a preliminary selection of the agitator type, based on the liquid viscosity and tank volume (Towler and Sinnott, 2013).

Figure 5. Low-speed agitators (Towler and Sinnott, 2013)

Figure 6. Agitator selection guide (Towler and Sinnott, 2013)

Gas-Liquid Mixing

Gases can be mixed into liquids using the inline mixing or stirred tank methods discussed previously. A special type of gas injector, called a sparger (shown in Figure 7) can also be used. This is a long injection tube with multiple holes drilled in it.

Figure 7. Gas sparger (Towler and Sinnott, 2013)

A small flow of liquid can be dispersed into a gas stream using a spray nozzle (Figure 8).

Figure 8. Liquid injection into gas (Towler and Sinnott, 2013)

Solid-Liquid Mixing

Solids are usually added to a liquid in a stirred tank at atmospheric pressure. In order to allow more accurate control of dissolved solid concentration, mixing of solids and liquids is often carried out as a batch operation (Towler and Sinnott, 2013).

Types of Reactors

Most reactors used in industry approximate the ideal batch reactor, PFR, or CSTR. In fact, real reactors can be modeled as networks or combinations of multiple plug-flow and stirred-tank reactors (Towler and Sinnott, 2013). Examples of real reactors that approximate the flow pattern of ideal reactors are shown in Figure 10. These reactors will be discussed in more detail in the following sections.

Figure 10. Ideal reactors and some real reactors that approximate the same flow pattern (Towler and Sinnott, 2013)

Vapor-Liquid Reactors

Vapor-liquid reactions are important in many chemical processes. For example, oxygenation and hydrogenation reactions are usually carried out with the organic component in the liquid phase (Towler and Sinnott, 2013). A summary of common goals for vapor-liquid reactors and the reactors used to achieve those goals is shown in Table 1.

Goal	Types of Vapor-Liquid Reactors	Examples
Maintain low concentration of gas component in liquid	Sparged stirred tank reactor Sparged tubular reactor	Liquid phase oxidations using air Fermenters
Contact gas and liquid over catalyst	Trickle bed reactor Slurry phase reactor	Catalytic hydrogenation
React a component out of the gas phase to high conversion	Multi-stage V/L contactor (reactive absorption column) Venturi scrubber	Chemisorption Acid gas scrubbing

Table 1. Summary of Vapor-Liquid Reactors (Towler, 2012)

If the residence time requirements are short enough, vapor-liquid contacting columns are preferred because of the high area for mass transfer. Trayed or packed columns can be used to contact vapor and liquid for reaction. The column packing may be catalytically active or could be inert packing (Towler, 2012). Please see the separation processes section of this website for more information on the types of processes used for the third goal listed.

Stirred tanks or tubular reactors are used when long residence time is needed for the liquid phase (Towler and Sinnott, 2013). These types of reactors and more will be discussed in the catalytic processes section of this page.

The reactors listed under the first goal in the table are unique to vapor-liquid processes. The basic concept of a sparger was discussed in the mixing section. Sparged reactors are shown in Figure 11.

Figure 11. Sparged stirred tank and tubular reactors (Towler, 2012)

The gas is bubbled up through the liquid in a sparged reactor. For smaller bubbles, a porous pipe diffuser can be used instead (Towler, 2012). The designer must allow some disengaging space at the top of the reactor, or entrainment will be excessive. If the gas flow rate is large then the gas flow can be used as the primary means of agitation. Perry's Handbook suggests the following air rates (ft³/ft²min) for agitating an open tank full of water at 1 atm:

Degree of agitation	Liquid depth 9 ft	Liquid depth 3 ft
Moderate	0.65	1.3
Complete	1.3	2.6
Violent	3.1	6.2

Table 2. Summary of suggested flow rates for gas flow as agitation (Towler, 2012)

Catalytic Processes

A catalyst increases the rate of a chemical reaction without itself becoming permanently changed by the reaction. Catalysts allow reactions to be run in smaller reactors and operated at lower temperatures and improve selectivity. Therefore, catalysts will almost always lead to a more economically attractive process than a noncatalytic route (Towler and Sinnott, 2013). Catalysts are normally selected based on performance rather than price since increases catalysts selectivity will almost always quickly pay back any price premium expected by the manufacturer. It is important to test the catalysts under conditions that are representative of process conditions (Towler and Sinnott, 2013).

Catalyst activity often deteriorates over time (Towler, 2012). Common causes of deactivation include:

• Poisoning by components in feed (e.g. base destroys acid catalyst)
• Blockage of pores or active sites by byproducts such as coke
• Thermal or hydrothermal modification of catalyst structure

Slow activity loss can be compensated by:

• Putting in more catalyst (lower space velocity)
• Slowly raising reactor temperature

Rapid activity loss may require moving the catalyst to a continuous regeneration zone (Towler, 2012).

Catalytic reactions can be either homogenous (catalyst is in the same phase as the reagents) or heterogeneous (catalyst is not in the same phase as the reagents).

Homogeneous Catalysis

Homogeneous catalysis can be conducted in the basic batch reactors, PFRs, or CSTRs that have already been discussed. However, when the catalyst is in the same phase as the reagent, recovering this catalyst after the reaction can be difficult and expensive, particularly if the catalyst is sensitive to high temperatures (Towler, 2012). Providing adequate interfacial area is also a challenge of homogeneous catalysis. A reaction often only occurs at the interface or in the boundary layer between the catalyst and the reagents. Increased mixing can increase the rate and selectivity of the reaction, but this can require detailed and expensive mixing equipment (Towler, 2012). For these reasons, reactions requiring homogenous catalysts are not usually used unless an easy separation can be found to recover the catalyst.

Heterogeneous Catalysis

Catalyst recovery in processes involving heterogeneous catalysis is much easier. However, the rate of reaction is limited by the available inter-phase surface area and the mass transfer of reagents and products to and from the interface (Towler, 2012). Therefore, reactors for these processes are design to reduce these limitations.

• Fixed Bed Reactors

In a fixed-bed reactor, the reagent flows over a stationary bed of packed catalyst (Towler and Sinnott, 2013). This is the most common type of reactor used for heterogeneous catalysis as long as the catalyst does not require continuous regeneration and the reaction mixture does not require high agitation (Towler, 2012). The amount of catalyst necessary can be found using the following equations:

The ratio of the bed height (L) to the diameter (D) determines the distribution of reagents and the pressure drop across the bed. An increased L/D ratio creates a more even distribution and less change of localized deactivation or "hot spots." However, increasing the L/D ratio increases the pressure drop, requiring higher compression and pumping costs (Towler, 2012). The Ergun equation can be used to calculate the pressure drop in packed beds.

Where V is the superficial velocity (volume flowrate divided by cross-sectional area), μ is the viscosity, D_p is the particle diameter and ε is the porosity of the packed bed (Towler, 2012). Given these trade-offs, it may make sense to split the catalyst over several beds (Towler, 2012).

• Radial Flow Reactors

When there is very little pressure drop available, the L/D ratio must be much less that one (Towler, 2012). A common solution to this is to use a radial flow reactor with the catalyst contained in an annulus between vertical perforated or slotted screens. The fluid flows radially through the bed and the direction of flow can be either inwards or outwards (Towler and Sinnott, 2013). An example of a radial flow reactor is shown in Figure 12.

Figure 12. Radial flow reactor (Towler, 2012)

• Moving Bed Reactors

A moving bed reactor is similar to a radial flow reactor, but the catalyst is moved through the annular space (Towler, 2012).

• Fluidized Bed Reactors

If the fluid flow is up through the catalyst bed then the bed can become fluidized if the pressure drop is high enough to support the weight of the catalyst. Fluidized beds usually have a lower pressure drop than down flow at high flow rates (Towler, 2012). In addition, fluidizing the catalyst eases the transition from one reaction zone to another.

The catalyst bed is fluidized using a distributor to inject fluidization fluid, which is not necessarily the feed. Fluidization occurs when the bed pressure drop balances the weight of the particles, or

Where ∆P is the pressure drop, ρ_p and ρ_g are the densities of the particle and gas respectively, ε_m is the porosity at minimum fluidization, and L is the height of the bed (Towler, 2012). Fluidization can only be used with relatively small sized particles (<300 micrometers with gases). The solid material must be strong enough to withstand attrition in the fluidized bed and cheap enough to allow for make-up to replace attrition losses (Towler and Sinnott, 2013). A fluidized-bed reactors must also make allowance for separating the fluid-phase product from entrained solids so that solids are not carried out of the reactor (Towler and Sinnott, 2013).

• Trickle Bed Reactors

Trickle bed reactors are used when all three phases are involved in the reaction. They must ensure good distribution of both the vapor and the liquid, without channeling of either phase (Towler, 2012). In a trickle bed reactor, the liquid flows down over the surface of a stationary bed of solids. The gas phase usually also flows downwards with the liquid, but countercurrent flow is feasible as long as flooding conditions are avoided (Towler and Sinnott, 2013). This requires a more sophisticated distributor like those used for packed distillation columns (Towler, 2012). An example of a trickle bed reactor is shown in Figure 13.

Figure 13. Example of trickle bed reactor (Towler, 2012)

• Slurry Reactors

Liquid is mixed up in the liquid in slurry phase reactions. Slurry reactors are prone to attrition of the solids, caused by pumping or agitation of the liquid (Towler and Sinnott, 2013). Slurry-phase operation is usually not preferred for processes that use heterogeneous catalysts because the catalyst tends to become eroded and can be difficult to recover from the liquid (Towler and Sinnott, 2013).

Bioreactors

Bioreactors have requirements that add complexity compared to simpler chemical reactors. These reactions often are three-phase (cells, water, and air), need sterile operation, and require heat removal (Towler, 2012). However, biological systems have the following advantages:

• Some products can only be made by biological routes
• Large molecules such as proteins can be made
• Selectivity for desired product can be very high
• Products are often very valuable

Enzyme Catalysis

Enzymes are the biological equivalent of catalysts. They can sometimes be isolated from host cells. They are usually proteins and, therefore, most are thermally unstable above ~60 degrees Celsius and active only in water at a restricted pH (Towler, 2012). Enzymes can sometimes be absorbed onto a solid or encapsulated in a gel without losing their structure. In this case, they can be used in a conventional fixed bed reactor. Typically, homogenous reactions are carried out in batch reactors.

Microorganism Design

As an alternative to an enzyme catalyst, engineered microorganisms can be used to produce a chemical of interest. This product could be complex biological compounds, therapeutic proteins, or commodity plastics and fuels (Westfall 2011). Host cells as a platform for modification have so far included bacteria, yeast, and mammalian cells (Schmidt 2005). Cells must be engineered to produce a heterologous product through recombinant DNA. For a therapeutic protein, this includes identifying the DNA sequence coding for the protein, and expressing that DNA in a cell host (Seider 2004). For a commodity molecule, enzymes that catalyze the synthesis of that molecule must be identified and expressed in the host cell.

Several features of an engineered cell line must be considered and potentially altered. These considerations are derived from the fact that a process must produce the product at high quantities, must enable growth at reasonable conditions, must be able to purify the product, and must be economically viable. Therefore, cells must grow to high concentrations at a reasonable temperature (usually between ambient and 37 C), and must grow on relatively cheap feedstocks such as glucose or methanol (Towler 2012). To achieve high product concentrations in batch or fed batch processes, it can be advantageous to only induce production of the product after cells have grown to a high concentration. This requires the heterologous DNA to be expressed with an inducible promoter. For example, production of a product could be induced when the feedstock is switched to methanol (Yurimoto 2000).

In addition to producing high amounts of product, it is important to design an organism to facilitate simple extraction of the purified product. For example, an organism can be designed to secrete a product into the surround media. This allows continuous fermentation, where media is fed to maintain cell growth as the product is separated off (Huang 2008). Another case involves the lysis of cells to extract a product that is not secreted. This process would be executed in batches, which can be timed optimally to maximize use of fermentation and separation equipment (Biegler 1997). With both of these options, cells need to tolerate high levels of product either intracellularly or in the surrounding environment.

The engineering of microorganisms presents a number of formidable challenges. Many companies avoid this issue by screening known microorganisms for strains that naturally produce high titers of product, or close precursors. Expressing heterologous genes in cells causes high stress, and disrupts natural metabolic balance. Metabolic engineering studies ways to relieve these stresses by “rewiring” synthesis networks within cells (Stephanopoulos 1998). It is also important to balance the feedstock to provide enough nutrients for growth and product synthesis. Finally, microorganisms often lack the appropriate machinery to create biologics from more complex organisms (like a human antibody). For example, glycosylation of proteins is a mechanism only recently achieved in bacteria (Nothaft, 2010).

Cell Growth

Cell growth goes through several phases during a batch, shown in Figure 15.

Figure 15. Cell growth and product formation in batch fermentation (Towler and Sinnott, 2013)

• I: Innoculation: slow growth while cells adapt to new environment
• II: Exponential growth: growth rate proportional to cell mass
• III: Slow growth as substrate or other factors begin to limit rate
• IV: Stationary phase: cell growth rate and death rate are equal
• V: Decline phase: cells die or sporulate, often caused by product build-up

Intracellular product accumulation is slow at first because there are a limited number of cells (Towler, 2012). However, it is important to note that product accumulation continue even after the live cell count falls, since dead cells still contain product.

The growth rate of cells can be limited by factors such as:

• The availability of the primary subtrate
  • Typically glucose, fructose, sucrose, or other carbohydrate
• The availability of other metabolites
  • Vitamins, minerals, hormones, or enzyme cofactors
• The availability of oxygen
• Mass transfer properties of the reaction system
• Inhibition or poisoning by products or byproducts
• High temperature caused by inadequate heat removal

All of these factors are exacerbated at higher cell concentrations (Towler, 2012). Clearly, biological reactions must be carefully controlled. An addition complication in dealing with biological reactions is that the product formation is often not closedly tied to the rate of consumption of the substrate (Towler, 2012). This is because of the fact that the product may be made by the cells at a relatively low concentration and the fact that some cell metabolic processes may not be involved in formation of the desired product (Towler, 2012).

Batch or Continuous

When designing a process involving biological organisms, it is important to select the right reactor method. The following equations represent the extremes:

Batch culture:

$r_{\text{b}}$ = rate of cell mass production in 1 batch cycle

Types of Bioreactors

• Stirred Tank Fermenter

The stirred tank fermenter is the most common reactor used for biological reactions (Towler, 2012) and is similar to the stirred tanks discussed previously. It can be used in both batch and continuous mode. Figure 14 shows a stirred tank fermenter.

Figure 14. Fermentation reactor (Towler and Sinnott, 2013)

• Shaftless Bioreactors

Shaftless bioreactors are used when the pump shaft seal is considered a non-permissible source of contamination. These reactors use gas flow to provide agitation of the liquid. The design requires careful attention to hydraulics (Towler, 2012). Examples of shaftless bioreactors are shown in Figure 15.

Figure 15. Examples of shaftless bioreactors (Towler, 2012)

• WAVE Bioreactors

WAVE bioreactors represent an alternative to standard stainless steel bioreactors. These reactors are flexible and single-use, cutting down time between batches and allow for a more sterile environment (citation). These disposable reactors are mostly used in mammalian cell culture. Three layers of plastic are the minimum necessary for construction. The first is a structural layer, followed by a barrier layer that allows for permeability. The last layer, the fluid contact layer, is designed to take into account inertness and maintain a good seal (citation).

Figure 16. Example of a WAVE bioreactor (citation).

• Packed Bed Bioreactors

Packed Bed Bioreactors are structured so that the cells are immobilized and placed on large particles. Although they are relatively simple to construct, they can have blockage issues or poor oxygen transfer. There are three types of flow: downward flow, upward flow and the recycling method. In industry, upward flow is preferred, especially when there is gas production during the fermentation (citation).

Figure 17. Example of a packed bed bioreactor (citation).

Preventing Contamination

Since cells are easily affected by both unwanted chemicals and other species in the reactor, bioreactors must be designed in order to avoid contamination (Towler, 2013). Bacterial spores are the most demanding sterilization challenge in a bioreactor. Bacterial spores are dormant and non-reproductive structures produced by a small number of bacteria. Since they are meant to ensure survival of bacteria in times of environmental stress, they are heat resistant. In order to ensure that all spores within the medium are killed, the medium must be sterilized (Shuler, 2002).

• Sterilization

The death kinetics involved in sterilization can usually be described by first order kinetics, but since essentially all contaminants need to be removed, it is usually described in probabilistic terms (Towler, 2013). There are many methods of sterilization, including filtering, chemical, thermal, and radiation. Thermal sterilization using steam is the most common method, as it is the most economical method for large scale reactors. Steam sterilizations either occur in the fermentation vessel as a batch sterilization or in a continuous apparatus upstream of the fermentation vessel (Seider, 2004).

Batch sterilization is used for smaller fermenters (Biegler, 1997). It is usually performed at 121˚C. The disadvantages of batch sterilization are thermal lags and incomplete mixing. The heat up and cool down times from 121˚C to 37˚C are usually longer than the time at the sterilization temperature and can damage the vitamins and protein in the media. The short-exposure and high temperature of continuous sterilization is easier to control, does less damage to the media, and reduces fermenter downtime. The disadvantages, however, are dilution of the medium with steam injection and foaming. (Shuler, 2002). A common process for continuous sterilization consists of steam injected into the media in order to heat it, passing the medium through a holding section to achieve the desired residence time, and then flash cooling the media. Flash cooling prevents contamination from cooling water (Towler, 2013). Diagrams and temperature profiles of batch and continuous sterilization is shown in Fig. 16.

Figure 16. Batch vs. Continuous Sterilization (Shuler, 2002)

If the media being sterilized contains heat-sensitive materials, filter sterilization must be used instead of steam. Filtration is also used to sterilize the process air used in the system. Microporous filters are used, so the media must be prefiltered for larger particulates so that the microporous filter does not get clogged. However, filtration is not as reliable as steam sterilization, as any defects in the membrane can lead to contamination, and viruses can often pass through the filter (Shuler, 2002).

• Cleaning

Cleaning the fermentation vessel at the end of the production run is necessary in order to remove residual substrates that could lead to contamination of future batches. Cleaning consists of multiple washes - washing with high-pressure water jets, an alkaline cleaning solution, tap water, acidic cleaning solution, tap water, and deionized water. The vessel is drained in between each of these steps. For this reason, the vessel should have no internal dead spots where material could accumulate. Also, due to the repeated emptying and filling of the vessel, cleaning leads to significant down time between batches (Towler, 2013).

Heating and Cooling of Reacting Systems

Exothermic and endothermic reactions will require reactors with heat control systems to prevent operating conditions from falling out of the desired range. Reactor performance is often limited by the ability to add or remove heat. Insufficient heat removal can cause runaway reactions, particularly dangerous situations in chemical processing (Turton et al., 2012). Before considering the design of a heating or cooling system to couple with a reactor, a few important questions should be asked (Towler and Sinnott, 2013).

1. Can the reaction be carried out adiabatically?

2. Can the feeds provide the required heating or cooling? Staged addition of feed can help alleviate the cost of adding a heat exchange network or heat transfer jacket. Also consider adding an inert diluent or hot/cold shots (Seider et al., 2004).

3. Would it be more cost effective to carry out the heat exchange outside of the reactor?

4. Would it be more effective to carry out the reaction inside of a heat transfer device? If a reaction requires only a small volume or small quantities of catalyst, it may be possible to utilize a heat exchanger as a temperature controller and as a reaction location.

5. Does the proposed design allow the process to be started up and shut down smoothly?

6. Are there safety concerns with heating or cooling the reactor?

After considering these aspects of the design, commercial design software such as HYSYS or UniSim can be utilized to estimate heating/cooling requirements. Once this is done, design of the heat exchange system can begin, with different reactor types and reactions requiring different design approaches (Towler and Sinnott, 2013).

Stirred Tank Reactors

Heating and cooling of a stirred tank reactor is done to ensure a uniform reaction temperature, so that there do not exist hot or cold spots within the reactor that can negatively affect selectivity (Towler and Sinnott, 2013).

For indirect heat transfer, there are three main alternatives: a heat transfer jacket, an internal coil, and an external heat transfer circuit. A jacket is utilized as long as there is sufficient heat transfer area for the heat exchange to take place. If this is not the case, coils are used, although the inclusion of a heating coil will significantly increase reactor volume and utility requirements, leading to a large increase in price for the reactor. External circuits contain a heat exchanger that will heat or cool the product stream as required and recycle this material to the reactor to control temperature. External circuits are useful because they can be designed independently of the reactor; sizing the required pumps and heat exchangers will not fundamentally change the activity of the reactor. For any of these choices, it should be ensure that no corrosion of the involved piping will occur, as utility streams bleeding into the reactor can have a very negative impact on the selectivity of the reaction and on the operation of the reactor on a whole (Towler and Sinnott, 2013).

Some direct heat transfer alternatives also exist, as long the reaction in question is compatible with the addition of extra water. Steam can be pumped into the reactor to maintain temperature, which will eliminate the need to design heat transfer surfaces. However, steam injected into the system cannot be recovered, so this will lead to an increase in annual utility costs. Additionally, vapor will be produced if it did not exist previously, so reactors will need to be redesigned to accommodate a vapor removal system (Towler and Sinnott, 2013).

Catalytic Reactors

Slurry Reactors

Since slurry reactors already use a mix of solid catalyst and liquid reactants, any of the methods described in the Stirred Tank Reactors section can be applied to slurry reactors. It is not recommended to use internal coils in such a design, as reactor slurry will often corrode heat exchange material very easily (Towler and Sinnott, 2013).

Fixed-bed Reactors

Indirect heat transfer is not often utilized to control the temperature in fixed-bed reactors, as it hard to maintain uniform temperature across the radial section of the catalyst bed. In cases where temperature control is required, the reactor will be split into smaller sections. After each bed, there will be an heat transfer stage, where the product stream is heated or cooled as necessary and returned to the next catalytic segment (Towler and Sinnott, 2013).

Fluidized-bed Reactors

Fluidized bed reactors have high heat-transfer coefficients, so indirect heat transfer is highly effective. The heat capacity of the solid catalyst particles can be used as a heat transfer medium themselves; heated catalyst contains a reaction location and the necessary heat to maintain the required temperature. Deactivated catalyst is heated during reactivation and recycle (Towler and Sinnott, 2013).

Heat Exchangers as Reactors

It is sometimes necessary to design a reactor as a heat transfer device, like when it is necessary to operate a reactor isothermally and there is a large heat of reaction. Some common situations include high-temperature endothermic reactions that quickly quench without continuous heat input and low-temperature exothermic reactions that must be kept at constant temperature to maintain selectivity. The most common heat transfer equipment used for reactions are shell and tube heat exchangers and fired heaters (Towler and Sinnott, 2013).

Homogenous Reactions

If the reaction does not required a catalyst, than the heat transfer design is the same as a conventional heat transfer device, with some important changes in the thermal design. The usual heat exchanger equations will not apply to the design of a heat exchanger reactor due to the nonlinear behavior of the reaction rate with regards to temperature. In these cases, the usual practice of conservative temperature estimations will not aid in heat transfer design, as greater detail will be required to ensure the proper operation of the reactor. Detailed kinetic models should be developed before designing the internals of the heat transfer device (Towler and Sinnott, 2013).

Heterogenous Reactions

The problems of designing for homogenous reactions still hold for heterogenous ones, with the added complication of solid catalyst beds. Catalyst can be loaded into the tubes of a shell and tube exchanger if the exchanger is mounted vertically and a suitable retaining screen is included at either end of the design. In this instance, hot catalyst can be reliably recycled and heat treated to reactivate the catalysts and reduce the presence of reactor hot spots. High-temperature endothermic reactions will be even more difficult to design for, as their heat requirements often exceed the amount provided by a heated catalyst. In these cases, a "tube in tube" design is utilized, where feed and catalyst are heated simultaneously by an external fired heater. This can be done as long as thermal expansion does not cause damage to the tubes, or else significant catalyst poisoning can occur. The same concerns as detailed in homogenous reactions will still apply for any design utilized for heterogenous ones, so it is again recommended to develop a detailed kinetic model before determining the amount of heat transfer required to maintain proper selectivity (Towler and Sinnott, 2013).

Safety Considerations in Reactor Design

Reactors require much attention to safety details in the design process due to the hazards they impose. They are often the highest temperature point in the process, heat of reaction may be released, and residence times can be long leading to a large inventory of chemicals. Guidelines exist for inherently safer design principles which seek to remove or reduce process hazards, limiting the impact of unforeseen events. These design methods should be applied throughout the design process as part of good engineering practice; they cannot be retroactively added by a process safety specialist. Some examples are given in the table below:

Some Applications of Inherently Safer Design Approaches in Reactor Design (Towler and Sinnott, 2013)

Exothermic reactions require special consideration due to their potential to runaway (temperature rises from heat of reaction being released, increasing reaction rate, releasing more heat, and so on). The reactor must be designed such that temperature can be precisely controlled and the reaction shut down if temperature control is lost. The use of solvents or inert species also allows for temperature control by adjusting heat capacity flow rate relative to rate of heat release from the reaction. An additional safety feature would allow the reactor to be flooded with cold solvent or diluent.

If there is a cooling system it should be designed to return the process to desired temperature if the maximum temperature is reached.

Venting and relief of reactors is complicated by the potential to keep reacting if containment is lost or material is discharged into the pressure relief system. The relief system should be designed according to guidelines outlined in the Design Institute for Emergency Relief Systems (DIERS) methodology. The reactor design team must understand the reaction mechanism and kinetics, including the role of any compounds which may accelerate the reaction. Details may be found on the AIChE website, here.

Capital Cost of Reactors

Reactors are classified as pressure vessels, and as such the pressure vessel design methods can be used to estimate wall thickness and thus determine capital cost. Additional costs come from reactor internals or other equipment. Jacketed stirred-tank reactors require more in depth analysis than that provided by pressure vessel design. The wall of the reaction vessel may be in compression due to the jacket. For preliminary cost estimating a correlation for jacketed stirred tank reactors operating at pressures below 20 bar can be used:

$C_{e}=a+bS^{n}$

where $C_{e}$ is the purchased equipment cost on a U.S. Gulf Coast Basis, $a,b$ are cost constants, $S$ is the size parameter, and $n$ is the exponent for that type of equipment. Values for $a,b,S,n$ are given in the table below:

Purchased Equipment Cost Factors (Towler and Sinnott, 2013)

Conclusions

The conversion of feed to products is the essence of a chemical process and, thus, the reactor is the heart of a chemical plant. When designing a reactor, an engineer must first collect data about the chemical reaction and then select appropriate reaction conditions, which will help determine suitable materials of construction. Next, the designer should determine the rate-limiting step and, from this, the critical sizing parameter. Next, preliminary sizing, layout, and costing can be conducted for the reactor. At this point, simulations and experiments can be conducted to verify that the proposed reactor will meet the desired specifications. The design is optimized until these targets are met. Throughout the design process, it is important for the engineer to consider the most appropriate type of reactor to use, any mixing or heat transfer equipment that must be added, and safety considerations.

References

Douglas JM. Conceptual Design of Chemical Processes. New York: McGraw-Hill; 1988.

Seider WD, Seader JD, Lewin DR. Process Design Principles: Synthesis, Analysis, and Evaluation. 3rd ed. New York: Wiley; 2004.

Towler G. Chemical Engineering Design. PowerPoint presentation; 2012.

Towler G, Sinnott R. Chemical Engineering Design: Principles, Practice and Economics of Plant and Process Design. 2nd ed. Boston: Elsevier; 2013.

Turton R, Bailie RC, Whiting WB, Shaewitz JA, Bhattacharyya D. Analysis, Synthesis, and Design of Chemical Processes. 4th ed. Upper Saddle River: Prentice-Hall; 2012.

Westfall, Patrick J., and Timothy S. Gardner. "Industrial fermentation of renewable diesel fuels." Current opinion in biotechnology 22.3 (2011): 344-350.

Schmidt, F. R. "Optimization and scale up of industrial fermentation processes." Applied microbiology and biotechnology 68.4 (2005): 425-435.

Yurimoto, Hiroya, et al. "Regulation and evaluation of five methanol-inducible promoters in the methylotrophic yeast Candida boidinii." Biochimica et Biophysica Acta (BBA)-Gene Structure and Expression 1493.1 (2000): 56-63.

Huang, Hua-Jiang, et al. "A review of separation technologies in current and future biorefineries." Separation and Purification Technology 62.1 (2008): 1-21.

Stephanopoulos, George, Aristos A. Aristidou, and Jens Nielsen. Metabolic engineering: principles and methodologies. Academic press, 1998.

Nothaft, Harald, and Christine M. Szymanski. "Protein glycosylation in bacteria: sweeter than ever." Nature Reviews Microbiology 8.11 (2010): 765-778.

Shuler, ML, Kargi, F. Bioprocess Engineering Basic Concepts. 2nd ed. Upper Saddle River: Prentice-Hall. 2002: 315-320

L.T. Biegler, I.E. Grossmann, A.W. Westerberg, Systematic Methods of Chemical ProcessDesign, Prentice-Hall: Upper Saddle River, 1997.