Turbulent Mixing and Chemical Reaction in Stirred Tanks: André Bakker Julian B. Fasano
Turbulent Mixing and Chemical Reaction in Stirred Tanks: André Bakker Julian B. Fasano
Turbulent Mixing and Chemical Reaction in Stirred Tanks: André Bakker Julian B. Fasano
in Stirred Tanks
André Bakker
Julian B. Fasano
Blend time and chemical product distribution in turbulent agitated vessels can be predicted with the
aid of Computational Fluid Mixing (CFM) models. The blend time predictions show good agreement
with an experimental correlation. Calculations for turbulent, time dependent mixing of two
chemicals, exhibiting a competitive pair of reactions, are compared with experimental results. The
effects of the position of the inlet feed stream in the turbulent flow field are studied. It is concluded
that process problems with turbulent chemical reactors can be avoided by incorporating the results
of CFM simulations in the design stage.
Keywords: Mixing, Chemical Reaction, Computational Modeling, Stirred Tank, Turbulent Flow.
Published in “The Online CFM Book” at http://www.bakker.org/cfm.
(c) 1998 André Bakker Updated: February 15, 2000
2 THE ONLINE CFM BOOK
INTRODUCTION
Blending of chemical reactants is a common operation in the chemical process industries. Blend time
predictions are usually based on empirical correlations. When a competitive side reaction is present,
the final product distribution is often unknown until the reactor is built. The effects of the position
of the feed stream on the reaction byproducts are usually unknown. Also, the scale up of chemical
reactors is not straightforward. Thus, there is a need for comprehensive, physical models that can be
used to predict important information like blend time and reaction product distribution, especially
as they relate to scale and feed position.
The objective of this paper is to investigate the extent to which Computational Fluid Mixing
(CFM) models can be used as a tool in the design of industrial reactors. The commercially available
program Fluent™ is used to calculate the flow pattern and the transport and reaction of chemical
species in stirred tanks. The blend time predictions are compared with a literature correlation for
blend time. The product distribution for a pair of competing chemical reactions is compared with
experimental data from the literature.
MODEL
The flow pattern is calculated from conservation equations for mass and momentum, in combination
with the Algebraic Stress Model (ASM) for the turbulent Reynolds stresses, using the Fluent V3
solver. These equations can be found in numerous textbooks and will not be reiterated here. Once
the flow pattern is known the mixing and transport of chemical species can be calculated from the
following model equation:
Here Xi is the mass fraction of chemical species I and Ri is the rate of creation or depletion by
chemical reaction. For a single step, first order reaction like A + B -> R the reaction rate is given
by:
Here CA and CB (upper case) denote the mean molar concentrations of reactants A and B while cA
and cB (lower case) denote the local concentration fluctuations that result from turbulence. When the
species are perfectly mixed the second term on the right hand side, containing the correlation of the
concentration fluctuations, will approach zero. Otherwise, if the species are not perfectly mixed, this
term will be negative and will reduce the reaction rate. The estimation of this correlation term is not
straightforward and numerous models are available. An excellent discussion on this subject was
given by Hannon [1].
The model used here is a slightly modified version of the standard Fluent model [2]. Two
possible reaction rates are calculated, the kinetic reaction rate Rki and a second reaction rate Rmi that
MIXING AND CHEMICAL REACTION 3
is controlled by the turbulent mixing. The kinetic reaction rate for species I is calculated as:
The turbulent mixing limited reaction rate for species I is calculated as:
The "minimum" function gives the minimum value of (ρ.Xj/ν.M j j) of all the reactants j taking part
in this reaction. Finally the reaction rate Ri is calculated as the product of the molar stoichiometry
νi of species I and the minimum of Rki and Rmi:
Here Mi is the molecular weight of species I and Amn is an empirically determined model constant
for reaction n. In the reaction system studied here, νi is +1 for reactants and -1 for products. K is the
kinetic rate constant of the reaction.
The idea behind this model is that in regions with high turbulence levels the eddy lifetime
k/ε will be short, mixing fast and as a result the reaction rate is not limited by small scale mixing.
On the other hand, in regions with low turbulence levels, small scale mixing may be slow and limit
the reaction rate.
This is the reaction system used by Bourne et al. [3] and Middleton et al. [4]. The first reaction is
much faster than the second reaction: K1 = 7300 m3.mole -1.s -1 vs. K2 = 3.5 m .mole 3
.s-1 . -1 The
experimental data published by Middleton et al. were used to determine the model constant Amn.
Two reactors were studied, a 30 l reactor equipped with a D/T=1/2 D-6 impeller and a 600 l reactor
with a D/T=1/3 D-6 impeller. A small volume of reactant B was instantaneously added just below
the liquid surface in a tank otherwise containing reactant A. A and B were added on an equimolar
basis. The transport, mixing and reaction of the chemical species were then calculated using the flow
pattern in Figure 1 as a basis. Experimental data were used as impeller boundary conditions. The
product distribution Xs is then calculated as:
In the reaction model used here it was assumed that small scale mixing only affected the first
reaction and that once this reaction had occurred, the species were locally well mixed. As a result,
Figure 6 The local mass fractions in the 2-D reaction simulation after 0, 1, 2, 4, 10
and 20 seconds.
6 THE ONLINE CFM BOOK
Figure 9 The concentration field in the 3-D blending simulation at 0, 4, 10 and 20s.
8 THE ONLINE CFM BOOK
BLEND TIME
The mixing of two non-reacting species in a tank equipped with a high efficiency impeller
(Chemineer HE-3) impeller was calculated using Fluent V3. The tank diameter was T = 1 m.
Furthermore, Z/T = 1; D/T = 0.33; C/T = 0.32 and RPM = 58. The flow pattern in this tank is shown
in Figure 7. Experimental data were used as impeller boundary conditions. Figure 8 shows the
uniformity of the mixture as a function of time. The model predictions are compared with the results
of the experimental blend time correlation of Fasano and Penney [6]. This graph shows that for a
uniformity above 90% there is excellent agreement between the model predictions and the
experimental correlation. Figure 9 shows the concentration field at t = 0, 4, 10 and 20 seconds
respectively. After 80 s the species are homogeneously mixed.
DISCUSSION
The models presented here correctly predicts blend time and reaction product distribution. The
reaction model correctly predicts the effects of scale, impeller speed and feed location. This shows
that such models can provide valuable tools for designing chemical reactors. Process problems may
be avoided by using CFM early in the design stage. When designing an industrial chemical reactor
it is recommended to determine the values of the model constants on a laboratory scale. The reaction
model can then be used to optimize the product conversion on the production scale varying agitator
speed and feed position.
However, the range of validation of the reaction model was limited. Only one impeller type
and one reaction system were studied. Future work has to concentrate on testing the model for a
wider range of geometries and reaction systems, and if necessary modify the model to increase it's
range of validity.
REFERENCES
NOTATION