CRE Notes AMU
CRE Notes AMU
Course Objective Application of basic Chemical Engineering Principles to understand the chemical
kinetics and reactor design.
Reference Books:
1. Smith,J.M., Chemical Engineering Kinetics, 5th Edition (2002),McGraw-Hill,New York.
2. Schmidt, Lanny D., the Engineering of Chemical Reactions, 2nd Edition (2005), Oxford University Press, Oxford,
New York.
3. Froment,G. F. and K.B. Bischoff,Chemical Reactor Analysis and Design,John Wiley & Sons,New York (1990).
4. Hill, Charles, G. Jr., An Introduction to Chemical Engineering Kinetics and Reactor Design, John Wiley & Sons, New
York (1977).
5. Holland, Charles D., and R.G. Anthony, Fundamentals of Chemical Reaction Engineering, 2nd Edition (1989), Prentice-
Hall Inc, New Jersey.
UNIT I:
OVERVIEW OF CHEMICAL REACTION ENGINEERING
• CRE is concerned with chemical treatment of the above shown typical chemical process.
• This chemical step is mainly consisted of inconsequential unit, or simply a mixing tank (or reactor) that makes or
breaks the process economically.
• Reactor design pre-requisites – chemical kinetics, thermodynamics, fluid mechanics, heat and mass transfers, and
economics.
• Chemical reaction engineering is the synthesis of all these factors with the aim of properly designing a chemical
reactor.
OVERVIEW OF CHEMICAL REACTION ENGINEERING
To know what a reactor can do, one must know the kinetics, the contacting patterns, and the performance
equation,as shown by figure 1.2.
-------- (1)
Output = f [input, kinetics,contacting pattern]
The above equation is known as performance
equation.
It relates input and output of the system
(reactor).
It can be used to
compare different designs and conditions,
Figure 1.2: Information needed to know what a Reactor can do. Find which is best, and then scale up to larger units.
OVERVIEW OF CHEMICAL REACTION ENGINEERING
Homogeneous and Heterogeneous Reactions:
In CRE, the classification of reactions can be made on the basis of types/number of phases involved in a reaction.
A homogeneous reaction is the reaction in which the reaction takes place in the same phase i.e., either solid, or
liquid, or gas.
A reaction is said to be heterogeneous if the reaction takes place with at least two phases involved.
Table 1.1: Classification of Chemical Reactions useful in Reactor Design
OVERVIEW OF CHEMICAL REACTION ENGINEERING
Rate of Reaction:
Selecting one reaction component “i” and defining the rate of reaction in terms of i.
Number of ways in which rate of reaction can be defined are as follows –
1. Rate of Reaction based on Volume of the Reacting Fluid:
-------- (2)
3. Based on Unit Interfacial Surface in two-fluids system or on Unit Surface of solid in solid-fluid system:
-------- (4)
OVERVIEW OF CHEMICAL REACTION ENGINEERING
4. Based on Unit Volume of Solid in Gas-Solid System:
-------- (5)
5. Based on Unit Volume of Reactor (if rate based on volume of reacting fluid is different):
-------- (6)
For homogeneous systems, volume of fluid is often identical to volume of reactor therefore, in that case equations
(2) and (6) can be used interchangeably.
In heterogeneous systems,all of the above equation can be encountered.
Equations (2) to (6) are interconnected to each other in the following manner –
OVERVIEW OF CHEMICAL REACTION ENGINEERING
or,
-------- (7)
Simple Reactors:
Also, the rates of reaction of all the elements are related to each other as –
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
2. Elementary and Non-Elementary Reactions: Reactions, in which rate equation corresponds to stoichiometric
equation of the system, are known as Elementary Reactions. For example, for a given single reaction, the
stoichiometric equation can be given as –
When there is no direct correspondence between stoichiometry and rate, then the reactions are known as non-
elementary reactions. The classical example of a non-elementary reaction is that between hydrogen and bromine,
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
The rate expression of given stoichiometric equation is given as –
The non-elementary reactions can be explained as a sequence of single reactions or elementary reactions.
3. Molecularity and Order of Reaction:
Molecularity of an elementary reaction refers to the number of molecules taking part in the reaction.
It can have a value of 1, 2, 3, … which means molecularity will have only integral value.
Molecularity is defined only for elementary reactions.
For a reaction,rate is given as –
The dimensions of the rate constant,k, for nth order reaction can be given as
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
For first order reaction, the rate constant has the following dimensions
Writing the rate equation in terms of equivalent measures does not have any affect on its order but will affect the
rate constant,k.
Elementary reactions are represented by equations that contain molecularity as well as rate constant in it.
Therefore,
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
6. Kinetic Models for Non-Elementary Reactions:
In non-elementary reactions, the stoichiometry does not match with the kinetics of the reaction therefore, a
multistep reaction model is developed to explain the kinetics.
A non-elementary reaction can be simplified as series of single elementary reactions to identify the kinetics easily.
For example –
b) Ions and Polar Substances: These are electrically charged atoms, molecules, or fragments of molecules known as
ions.They may acts as active intermediates in reactions. For example –
R is said to be the active intermediate but if its activity is increased then its quantity becomes very small or
unobservable as it reacts to form B very fast therefore in such situations, R is referred to as reactive intermediate.
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
Postulated reaction schemes involving these three kinds of intermediates can be of two types.
Non-chain Reactions: In the non-chain reaction the intermediate is formed in the first reaction and then
disappears as it reacts further to give the product.Thus,
Chain Reactions: In chain reactions the intermediate is formed in a first reaction, called the chain initiation step. It
then combines with reactant to form product and more intermediate in the chain propagation step. Occasionally the
intermediate is destroyed in the chain termination step.Thus,
The essential feature of the chain reaction is the propagation step. In this step the intermediate is not consumed but
acts simply as a catalyst for the conversion of material. Thus, each molecule of intermediate can catalyze a long chain
of reactions, even thousands,before being finally destroyed.
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
Examples pertaining to the above mentioned mechanisms are:
a. Free Radicals,Chain Reaction Mechanism:The reaction
and can be explained by the following scheme which introduces and involves the intermediates H 2. and Br2.,
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
b. Molecular intermediates,Non-chain mechanism:The general class of enzyme- catalyzed fermentation reactions
In such reactions the concentration of intermediate may become more than negligible, in which case a special
analysis,first proposed by Michaelis and Menten (1913),is required.
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
7. Testing Kinetic Models:
There are 2 major complex problems associated with the kinetics mechanism postulated earlier.
i. The reaction may proceeds with more than one reaction mechanisms.
ii. More than one mechanism can be consistent with the kinetic data.
While considering the elementary reactions, two types of hypothesis regarding types of intermediates are made.
Type 1: An unseen and unmeasured intermediate X usually present at such small concentration that its rate of
change in the mixture can be taken to be zero.Thus,we assume
where,
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
Temperature DependentTerm:
1. Temperature dependency from Arrhenius’ Law: For many reactions, and particularly elementary reactions, the rate
expression can be written as a product of a temperature-dependent term and a composition- dependent term,
or –
ri = f1(temperature) . f2(composition)
ri = k . f2(composition)
where,k in most of the cases,as predicted experimentally, is given by Arrhenius’ equation –
Where, k0 is the frequency or pre exponential factor, E is the activation energy of reaction and, R is the universal
gas constant.
KINETICS/REACTION MECHANISM OF HOMOGENEOUS REACTIONS
At the same concentration but different temperatures and at constant
E,Arrhenius’ law indicates that –
At a first sight this could very well be an elementary biomolecular reaction with two molecules of ammonia
combining to yield directly the four product molecules. From this principle, however, the reverse reaction would
then also have to be an elementary reaction involving the direct combination of three molecules of hydrogen with
one of nitrogen. Because such a process is rejected as improbable, the bimolecular forward mechanism must also
be rejected.
6. For multiple reactions a change in the observed activation energy with temperature indicates a shift in the
controlling mechanism of reaction. Thus, for an increase in temperature Eobs rises for reactions or steps in
parallel, Eobs falls for reactions or steps in series. Conversely, for a decrease in temperature Eobs falls for reactions
in parallel,Eobs rises for reactions in series.These findings are illustrated in Fig. 2.3.
KINETICS/REACTION
MECHANISM OF
HOMOGENEOUS
REACTIONS
Figure 2.3 A change in activation energy indicates a shift in controlling mechanism of reaction.
CHAPTER 3:
INTERPRETATION OF BATCH REACTOR DATA
The determination of the rate equation is usually a two step procedure:
1. The concentration dependent term is found at fixed temperature then,
2. The temperature dependent term is found yielding the complete rate equation.
Equipment by which empirical information is obtained can also be divided into two categories: the batch and flow
reactors.
Batch reactor is a simple container which holds the contents while they react with each other.
To get the empirical information, extent of reaction at various times is to be calculated. This can be done in the
following ways –
1. By following the concentration of a given component.
2. By following change in some physical property of the fluid (example refractive index, electrical conductivity etc.).
3. For a constant-volume system, a change in total pressure can be followed.
4. For a constant-pressure system,a change in volume can be followed.
INTERPRETATION OF BATCH REACTOR DATA
Experimentally, the batch reactor operates at isothermal and constant volume process conditions as it is easy to
interpret homogeneous kinetic data for batch reactor at these conditions.
It is a simple device adaptable to small-scale lab set-ups with very little auxiliary equipment or instrumentation.
There are two methods for interpreting kinetic data of a batch reactor – the integral and the differential.
In the integral method of analysis we guess a particular form of rate equation and, after appropriate integration
and mathematical manipulation, predict that the plot of a certain concentration function versus time should yield a
straight line. The data are plotted, and if a reasonably good straight line is obtained, then the rate equation is said
to satisfactorily fit the data.
In the differential method of analysis we test the fit of the rate expression to the data directly and without any
integration. However, since the rate expression is a differential equation, we must first find (1/V)(dN/dt) from the
data before attempting the fitting procedure.
The integral method is easy to use and is recommended when testing specific mechanisms, or relatively simple
rate expressions, or when the data are so scattered that we cannot reliably find the derivatives needed in the
differential method.
INTERPRETATION OF BATCH REACTOR DATA
The differential method is useful in more complicated situations but requires more accurate or larger amounts of
data.
The integral method can only test this or that particular mechanism or rate form; the differential method can be
used to develop or build up a rate equation to fit the data.
In general, integral analysis is to be attempted first, and, if not successful, then the differential method can be tried.
-------- (3.2)
It is evident from the above equations that rate of reaction of any component is given by the rate of change of its
concentration or partial pressure (only for constant volume systems).
For gas reactions with changing numbers of moles, rate of reaction can be obtained by following the change in
total pressure of the system.
1. Analysis of Total Pressure Data Obtained in a Constant-Volume System: For isothermal gas reactions, writing
the general stoichiometric equation, and under each term indicating the number of moles of that component:
INTERPRETATION OF BATCH REACTOR DATA
Total number of moles present in the system initially –
where,
Assuming Ideal gas law and writing for reactant A in terms of system volume V –
-------- (3.4)
Combining eqns. 3.3 and 3.4, we get –
or -------- (3.5)
INTERPRETATION OF BATCH REACTOR DATA
Equation 3.5 gives the concentration or partial pressure of reactant A as a function of the total pressure π at time t,
initial partial pressure of A, pA0, and initial total pressure of the system, π0.
Similarly, for any product R we can find –
-------- (3.6)
Equations 3.5 and 3.6 are the desired relationships between total pressure of the system and the partial pressure of
reacting materials.
If in case, the precise stoichiometry is not known, or if more than one stoichiometric equation is needed to
represent the reaction, then the "total pressure" procedure cannot be used.
2. The Conversion: Also known as fractional conversion, or the fraction of any reactant, say A, converted to
product, or the fraction of A reacted away.
In simple words, it can be called out as, conversion of A, with symbol XA.
INTERPRETATION OF BATCH REACTOR DATA
Suppose that NA0 is the initial amount of A in the reactor at time t = 0, and that NA is the amount present at time t.
Then the conversion of A in the constant volume system is given by:
-------- (3.7)
and,
-------- (3.8)
-------- (3.9)
Let us consider the below mentioned first-order rate equation and test it for integral analysis –
-------- (3.10)
or,
-------- (3.11)
or,
-------- (3.12)
-------- (3.13a)
with corresponding rate equation,
-------- (3.13b)
The amounts of A and B that have reacted at any time t are equal and given by CA0XA, we may write Eqs. 3.13a and
3.13b in terms of XA as –
INTERPRETATION OF BATCH REACTOR DATA
Letting M = CB0/CA0 be the initial molar ratio of reactants, we obtain –
After breakdown into partial fractions, integration, and rearrangement, the final result in a number of different
forms is
-------- (3.14)
Figure 3.2 shows two equivalent ways of obtaining a linear plot between the concentration function and time for
this second-order rate law.
INTERPRETATION OF BATCH REACTOR DATA
Figure 3.2: Test for Bimolecular Mechanism with CA0 does not equals to CB0, or for the second-order reaction
INTERPRETATION OF BATCH REACTOR DATA
If CB0 is much larger than CA0, CB remains approximately constant at all times, and Eq. 3.14 approaches Eq. 3.11
or 3.12 for the first-order reaction. Thus, the second- order reaction becomes a pseudo first-order reaction.
Case 1: for the second-order reaction with equal initial concentrations of A and B, or for the reaction
-------- (3.15a)
Figure 3.3: Test for Bimolecular Mechanism with CA0 equals to CB0, or for the second-order reaction
INTERPRETATION OF BATCH REACTOR DATA
Case 1I: If the reaction is first order with respect to both A and B, hence second order overall, or
-------- (3.17a)
-------- (3.17b)
-------- (3.19)
These two cases apply to all reaction types. Thus, special forms for the integrated expressions appear whenever
reactants are used in stoichiometric ratios, or when the reaction is not elementary.
INTERPRETATION OF BATCH REACTOR DATA
3. Irreversible Trimolecular-Type Third-Order Reactions: Consider the trimolecular reaction –
-------- (3.20a)
with corresponding rate equation,
-------- (3.20b)
in terms of XA as –
On separation of variables, breakdown into partial fractions, and integration, we obtain after manipulation
-------- (3.21)
INTERPRETATION OF BATCH REACTOR DATA
Now if CD0 is much larger than both CA0 and CB0 the reaction becomes second order and Eq. 3.21 reduces to Eq.
3.14.
All trimolecular reactions found so far are of the form of Eq. 3.22 or 3.25. Thus,
With,
-------- (3.22)
-------- (3.23)
-------- (3.25)
The order “n” cannot be found explicitly from Eq. 3.29, so a trial-and-error solution must be made.
Just select a value for n and calculate k. The value of n which minimizes the variation in k is the desired value of n.
Reactions with order n > 1 can never go to completion in finite time.
For orders n < 1 this rate form predicts that the reactant concentration will fall to zero and then become negative at
some finite time.
INTERPRETATION OF BATCH REACTOR DATA
Since the real concentration cannot fall below zero we should not carry out the integration beyond this time for n
< 1. Also, as a consequence of this feature, in real systems the observed fractional order will shift upward to unity
as reactant is depleted.
5. Zero Order Reactions: A reaction is said to be of zero order when the rate of reaction is independent of the
concentration of materials, that is –
-------- (3.30)
Integrating the above equation in the light that CA can never becomes negative we will get –
-------- (3.31)
INTERPRETATION OF BATCH REACTOR DATA
The above equation can be plotted as follows:
6. Overall Order of Irreversible Reactions from the Half-Life t1/2: Consider an irreversible reaction as follows –
If the reactants are present in their stoichiometric ratios, they will remain at that ratio throughout the reaction.
INTERPRETATION OF BATCH REACTOR DATA
Thus, for reactants A and B at any time CB/CA = β/α, and we may write –
Or,
-------- (3.32)
Integrating for n ≠ 1 gives,
Defining the half-life of the reaction, t1/2, as the time needed for the concentration of reactants to drop to one-half
the original value, we obtain –
-------- (3.33a)
INTERPRETATION OF BATCH REACTOR DATA
This expression shows that a plot of log t,, vs. log CA0 gives a straight line of slope 1 - n, as shown in Fig. 3.5.
From the above figure, it can be concluded that the fractional conversion in a given time:
rises with increased concentration for orders greater than one,
drops with increased concentration for orders less than one, and
is independent of initial concentration for reactions of first order.
INTERPRETATION OF BATCH REACTOR DATA
6. Fractional Life Method tF:
The half-life method can be extended to any fractional life method in which the concentration of reactant drops to
any fractional value F = CA/CA0 in time tF .
-------- (3.33b)
7. Irreversible Reactions in Parallel: Consider the simplest case, A decomposing by two competing paths, both
elementary reactions:
-------- (3.34)
INTERPRETATION OF BATCH REACTOR DATA
-------- (3.35)
-------- (3.36)
The k values are found using all three differential rate equations. First of all, Eq. 34, which is of simple first order,
is integrated to give –
-------- (3.37)
The above eqn. when plotted yields, (k1 +k2) as slope, as shown in figure 3.6.
Now taking eqns. 3.35 and 3.36 and dividing them will results in the following ratio:
-------- (3.38)
The slope of a plot (figure 3.6) of CR versus CS gives the ratio k1/k2. Knowing k1/k2 as well as k1 + k2, gives k1 and
k2. Typical concentration-time curves of the three components in a batch reactor for the case where CR0 = CS0 = 0
and k1 > k2 are shown in Fig. 3.7.
Figure 3.6: Evaluation of the Rate Constants for two competing Figure 3.7: Typical Concentration-Time curves for Competing
Elementary First Order Reactions Reactions
INTERPRETATION OF BATCH REACTOR DATA
8. Homogeneous Catalyzed Reactions: Consider the reaction rate for a homogeneous catalyzed system as the sum of
rates of both the uncatalyzed and catalyzed reactions –
-------- (3.39)
INTERPRETATION OF BATCH REACTOR DATA
On integrating eqn. 3.39 with taking the fact into consideration that the catalyst concentration remains unchanged
throughout the reaction, yields –
-------- (3.40)
Making a series of runs with different catalyst concentrations allows us to find k1 and k2.
Plot the observed k value against the catalyst concentrations as shown in Fig. 3.8. The slope of such a plot is k2
and the intercept k1.
If there is no change in the number of moles of A and R, we may write that at any time:
In terms of initial reactant ratio M = CR0/CA0 and fractional conversion of A, this can be written as –
-------- (3.43)
If the reaction is to proceed, then some product R must be present in the reactor.
As the formation of R increases, the rate of reaction starts to increase.
When A is just used up or get exhausted in reaction, the rate of reaction starts to decrease and must drop to zero as
there is no more reactant present in the reactor for further conversion.
The rate of reaction is at its highest when A and R are present equally.
This characteristic of an autocatalytic reaction is shown in Figure 3.9
INTERPRETATION OF BATCH REACTOR DATA
-------- (3.44)
-------- (3.45)
-------- (3.46)
Starting with a concentration CA0 of A, no R or S present, and by integration of eqn. 3.44 the concentration of A is
found to be:
-------- (3.47)
INTERPRETATION OF BATCH REACTOR DATA
To get the changing concentration of R, insert the value of CA from eqn. 3.47 in to eqn. 3.45 and rearrange the eqn.
-------- (3.48)
By multiplying through with the integrating factor, the solution can be obtained as follows:
Applying this general procedure to the integration of Eqn. 3.48, the integrating factor is found to be ek2t.
The constant of integration is found to be –k1CA0/(k2 – k1) from the initial conditions CR0 = 0 at t = 0, and the final
expression for the changing concentration of R is given as –
INTERPRETATION OF BATCH REACTOR DATA
-------- (3.49)
Thus eqns. 3.47, 3.49, and 3.50 gives the variation of concentrations of components A, R, and S with time t.
Case I: If k2 much larger then k1, then eqn. 3.50 reduces to –
Or, the rate is determined by k1 or the first step of the two-step reaction.
INTERPRETATION OF BATCH REACTOR DATA
Case I: If k1 much larger then k2, then eqn. 3.50 reduces to –
Or, the rate is determined by k2 or the second step of the two-step reaction, which is the slowest step.
In general, for any number of reactions in series it is the slowest step that has the greatest influence on the overall
reaction rate.
The maximum concentration of R may be found out by differentiating eqn. 3.49 and setting dCR/dt = 0. Therefore,
the time at which the maximum concentration of R occurs is gives by –
-------- (3.51)
Therefore, the maximum concentration of R is found by combining eqns. 3.49 & 3.51 as –
-------- (3.52)
INTERPRETATION OF BATCH REACTOR DATA
Figure 3.11 shows the general characteristics of the
concentration-time curves for the three components.
A decreases exponentially, R rises to a maximum and then
falls, and S rises continuously, the greatest rate of increase of
S occurring where R is a maximum.
In particular, this figure shows that one can evaluate k1 and k2
by noting the maximum concentration of intermediate and the
time when this maximum is reached.
-------- (3.53a)
Starting with a concentration ratio M = CB0/CA0 the corresponding rate equation can be given as –
-------- (3.53b)
Now at equilibrium dC A/dt = 0. Hence from Eqn. 3.53b we find the fractional conversion of A at equilibrium
conditions to be –
With conversions measured in terms of XAe, this may be looked on as a pseudo first-
order irreversible reaction which on integration gives
-------- (3.54)
With restrictions that CA0 = CB0 and CR0 = CS0 = 0, the integrated form of
the final eqns. Can be given as –
-------- (3.56)
Figure 3.14: Test for the bimolecular type
reversible reactions of Eq. 3.55.
A plot between conversion and time is shown in fig. 3.14.
INTERPRETATION OF BATCH REACTOR DATA
13. Reversible Reactions in General:. For orders other than one or two, integration of the rate equation becomes
cumbersome. Then search for an adequate rate equation is best done by the differential method.
14. Reactions of Shifting Orders:. In searching for a kinetic equation it may be found that the data are well fitted by
one reaction order at high concentrations but by another order at low concentrations. Consider the reaction
-------- (3.57)
From the above rate eqn. one can observe that (fig. 3.15)
I. At high CA (or k2CA is much larger than 1) – the reaction is of zero order with rate constant k1/k2.
II. At low CA (or k2CA is much lower than 1) – the reaction is of first order with rate constant k1.
To apply the integral method, separate variables and integrate eqn. 3.57 which gives –
-------- (3.58a)
INTERPRETATION OF BATCH REACTOR DATA
On rearranging, eqn. 3.58a can be re-written as –
-------- (3.58b)
or,
-------- (3.58c)
INTERPRETATION OF BATCH REACTOR DATA
-------- (3.59)
INTERPRETATION OF BATCH REACTOR DATA
Equation 3.59 shifts from order m – n at high concentration to order m at low concentration, with the transition
taking place at k2CnA ≈ 1. this type of equation can then be used to fit data of any two orders.
Another form of equation which can account for this shift is –
-------- (3.60)
or,
(b) testing an nth-order form -rA = kCnA by taking logarithms of the rate equation (Fig. 3.18).
INTERPRETATION OF BATCH REACTOR DATA
Differential Method of Analysis of Data
Figure 3.17: Test for the particular rate form –rA = Figure 3.18: Test for an nth – order rate form by the
kf(CA) by the differential method. differential method.
INTERPRETATION OF BATCH REACTOR DATA
For graphical testing in differential analysis, some mathematical manipulations are required to get the straight line
curve.
For example, consider a set of CA vs t data to which we have to fit M-M equation:
-------- (3.57)
Eqn. 3.57 has already been treated by integral analysis. However to perform graphical testing the above eqn. can
be manipulated by taking reciprocals as –
-------- (3.61)
-------- (3.62)
INTERPRETATION OF BATCH REACTOR DATA
A plot of -rA vs. (-rA)/CA is linear, as shown in Fig. 3.19.
The progress of the reaction is followed by noting the movement of the bead with time, a much simpler procedure
than trying to measure the composition of the mixture, especially for microreactors. Thus,
V0 = initial volume of the reactor
V = the volume at time t.
INTERPRETATION OF BATCH REACTOR DATA
This kind of reactor can be used for isothermal constant pressure operations, of reactions having a single
stoichiometry. For such systems the volume is linearly related to the conversion, or –
-------- (3.63a)
Or,
-------- (3.63b)
Where, ƐA is the fractional change in volume of the system between no conversion and complete conversion of
reactant A. Therefore,
-------- (3.64)
But with 50% inerts into consideration will change the fractional change in volume as two volumes of reactant
mixture yield, on complete conversion, five volumes of product mixture –
ƐA accounts for both the reaction stoichiometry and the presence of inerts.
Now,
-------- (3.65)
INTERPRETATION OF BATCH REACTOR DATA
Therefore, -------- (3.66)
Eqn. 3.66 is the relationship between conversion and concentration for isothermal varying-volume (or varying-
density) systems satisfying the linearity assumption of eqn. 3.63.
Now, the rate of reaction in general is given by –
Replacing V from eqn. 3.63a and NA from eqn. 3.65 we will get final rate equation as –
-------- (3.68)
On integration, it gives,
-------- (3.70)
INTERPRETATION OF BATCH REACTOR DATA
As shown in Fig. 3.21, the logarithm of the fractional change in volume versus time yields a straight line of slope
kƐA/CA0 .
INTERPRETATION OF BATCH REACTOR DATA
b) First order Reactions: For a unimolecular type first-order
reaction the rate of change of any reactant A is given as:
-------- (3.71)
Or,
nth-Order: For all rate forms other than zero-, first-, and second-order the integral method of analysis is not useful.
INTERPRETATION OF BATCH REACTOR DATA
Temperature and Reaction Rate:
So far we have studied the effect of concentration of reactants and products on the rate of reaction, all at a given
temperature.
To obtain the complete rate equation, we also need to know the role of temperature on reaction rate.
The concentration dependent terms f(C) usually remain unchanged at different temperature.
Chemical theory predicts that the rate constant should be temperature-dependent in the following manner:
However, since the exponential term is much more temperature-sensitive than the power term, therefore we can
reasonably consider the rate constants to vary approximately as e-E/RT.
Thus we can examine for the variation of the rate constant with temperature by an Arrhenius-type relationship –
Summary:
In chemistry, reaction rates and activation energies are developed in
terms of concentration.
For homogeneous reactions, E and –rA normally depends on
concentrations as the rate constants unit are often written as s-1,
liter/mole.sec, etc., without pressure appearing in the units.
It is better to convert p values into C values while making runs at
different temperatures.
INTERPRETATION OF BATCH REACTOR DATA
The Search for a Rate Equation:
All the rate equations are manipulated mathematically into a linearized
form because of the simplicity of family of straight lines that allows it to
be tested and rejected.
Three methods are commonly used to test for the linearity of a set of
points. These are as follows:
Calculation of k from Individual Data Points: For a given rate equation, the
rate constant can be determined for each experimental point by either
integral or differential method. If there is no appreciable trend in k values
then the equation is said to be satisfactory and the k values are averaged.
k values are calculated by taking the slopes of the lines connecting
individual data point to the origin. So for the same magnitude, k values
calculated near the origin (low conversion) will vary widely as compared
to the k values calculated farther from the origin. This will create confusion
about being the k value, a constant.
4. INTRODUCTION TO REACTOR DESIGN
In reactor design, generally the
study perform is to know what size
and type of reactor and method of
operation are best for a given job.
There are many factors which must
be accounted for in predicting the
performance of a reactor.
Equipment in which homogeneous
reactions are carried out can be one
of the three general types; the batch,
the steady-state flow, and the
unsteady-state flow or semibatch
reactor. The semibatch reactor
includes all reactors that do not fall
into the first two categories.
4. INTRODUCTION TO REACTOR DESIGN
(rate of reactant flow into element of volume) = (rate of reactant flow out of element of volume) + (rate of loss due to
chemical reaction within the element of volume) + (rate of
accumulation of reactant in element of volume) ---- (4.1)
4. INTRODUCTION TO REACTOR DESIGN
For uniform reactor composition (independent of position), the balance can be made over the whole reactor.
For non-uniform reactor composition, the balance must be made over a differential element of volume and then
integrated across the whole reactor for the appropriate flow and concentration conditions.
For various reactors, eqn. 4.1 ca be simplified and integrated so that the resultant expression gives the basic
performance equation of that particular type of reactor.
In batch reactor, the first two terms of eqn. 4.1 equals to zero as there is no inflow and outflow of the materials.
Similarly for steady state flow reactors, the last term of accumulation disappears and, for semi-batch reactor, all of
the four terms of the equation have to be considered.
For non-isothermal operations/processes an energy balances are used in addition to the material balances.
Therefore –
(rate of heat flow into element of volume) = (rate of heat flow out of element of volume) + (rate of disappearance of
heat by reaction within the element of volume) + (rate of accumulation of
heat within element of volume) ---- (4.2)
4. INTRODUCTION TO REACTOR DESIGN
A satisfactory design for a given problem is achieved only when –
a) We can predict the response of the reacting system to changes in operating conditions (how rates and equilibrium conversion
change with temperature and pressure),
b) We are able to compare yields for alternative designs (adiabatic vs. isothermal operations, single vs. multiple reactors, flow
vs. batch system), and
c) We can estimate the economics of these different alternatives.
There are two related measures of the extent of reaction, the concentration CA and the conversion XA.
The relationship between these two measures depends on a number of factors. Consider the three special cases in
order to understand this relationship.
4. INTRODUCTION TO REACTOR DESIGN
4. INTRODUCTION TO REACTOR DESIGN
Special Case 1: Constant Density Batch and Flow Systems
Includes most liquid reactions and also gas reactions (which runs at constant temperature and density).
CA and XA are related to each other as:
-------- (4.3)
-------- (4.4)
4. INTRODUCTION TO REACTOR DESIGN
Special Case 2: Batch and Flow Systems of Gases of Changing Density but with T and π Constant
In this case, the density changes because of the change in number of moles during reaction. In addition, the
volume of a fluid element changes linearly with conversion as:
-------- (4.5)
-------- (4.6)
4. INTRODUCTION TO REACTOR DESIGN
Special Case 3: Batch and Flow Systems of Gases in General (varying density, T, and π)
In this case, consider the following reaction which react according to –
Pick one reactant (called as key reactant) as the basis for determining the conversion. Let A be the key reactant,
then for the ideal gas behavior,
-------- (5.1)
-------- (5.3)
This is the general equation showing the time required to achieve a conversion XA for either isothermal or non-
isothermal operation. The volume of reacting fluid and the reaction rate remain under the integral sign, for in
general they both change as reaction proceeds.
This equation may be simplified for a number of situations. If the density of the fluid remains constant, we obtain:
-------- (5.4)
For all reactions in which the volume of reacting mixture changes proportionately with conversion, such as in
single gas-phase reactions with significant density changes, eqn. 5.3 becomes –
-------- (5.5)
5. IDEAL REACTORS FOR A SINGLE REACTION
Eqns. 5.2 – 5.5 are applicable to both isothermal and non-isothermal operations.
For the non-isothermal operations, the variation of rate with temperature, and the variation of temperature with
conversion, must be known before solution is possible.
Figure 5.2 is a graphical representation of two of these equations.
5. IDEAL REACTORS FOR A SINGLE REACTION
Space time and Space velocity:
Just as the reaction time t is the natural performance measure for a batch reactor, so are the space-time and space-
velocity the proper performance measures of flow reactors. These terms are defined as follows:
Space time:
-------- (5.6)
Space velocity:
-------- (5.7)
Thus, a space-velocity of 5 hr-l means that five reactor volumes of feed at specified conditions are being fed into
the reactor per hour. A space-time of 2 min means that every 2 min one reactor volume of feed at specified
conditions is being treated by the reactor.
5. IDEAL REACTORS FOR A SINGLE REACTION
Space time and Space velocity:
We may arbitrarily select the temperature, pressure, and state (gas, liquid, or solid) at which we choose to measure
the volume of material being fed to the reactor.
If they are of the stream entering the reactor, the relation between s and r and the other pertinent variables is
-------- (5.8)
It may be more convenient to measure the volumetric feed rate at some standard state, especially when the reactor
is to operate at a number of temperatures. For example, the material is gaseous when fed to the reactor at high
temperature but is liquid at the standard state, care must be taken to specify precisely what state has been chosen.
5. IDEAL REACTORS FOR A SINGLE REACTION
The relation between the space-velocity and space-time for actual feed conditions (unprimed symbols) and at
standard conditions (designated by primes) is given by
-------- (5.9)
-------- (5.10)
Now, evaluating each term of eqn. 5.10, if FA0 = v0CA0 is the molar feed rate of component A to the reactor, then
considering the reactor as a whole we have –
Input of A, moles/time = FA0(1 – XA) = FA0
Output of A, moles/time = FA0(1 – XA)
5. IDEAL REACTORS FOR A SINGLE REACTION
Disappearance of A by reaction, moles/time = (–rA)V
Introducing the above evaluated terms in eqn. 3.10 we get:
More generally, if the feed on which conversion is based, subscript 0, enters the reactor partially converted,
subscript i, and leaves at conditions given by subscript f, we have
or -------- (5.12)
5. IDEAL REACTORS FOR A SINGLE REACTION
For the special case of constant-density systems XA = 1 – CA/CA0, in which case the performance equation for mixed
reactors can also be written in terms of concentrations or,
These expressions relate in a simple way the four terms XA, -rA, V, FA0 ; thus, knowing any three allows the fourth to be
found directly.
In design, then, the size of reactor needed for a given duty or the extent of conversion in a reactor of given size is found
directly.
In kinetic studies each steady-state run gives, without integration, the reaction rate for the conditions within the reactor.
The ease of interpretation of data from a mixed flow reactor makes its use very attractive in kinetic studies, in particular
with messy reactions (e.g., multiple reactions and solid catalyzed reactions). Figure 5.4 is a graphical representation of
these mixed flow performance equations. For any specific kinetic form the equations can be written out directly.
5. IDEAL REACTORS FOR A SINGLE REACTION
Figure 5.4 is a graphical representation of these mixed flow performance equations. For any specific kinetic form
the equations can be written out directly.
5. IDEAL REACTORS FOR A SINGLE REACTION
As an example, for constant density systems CA/CA0 = 1 – XA, thus the performance expression for first-order
reaction becomes –
-------- (5.14a)
Therefore, first order reaction with variable volume condition, the performance 5.11 can be rewritten as:
-------- (5.14b)
For second-order reaction,A → products, -rA, = kC2A, ƐA = 0, the performance equation of eqn. 5.11 becomes
-------- (5.15)
5. IDEAL REACTORS FOR A SINGLE REACTION
Similar expressions can be written for any other form of rate equation. These expressions can be written either in
terms of concentrations or conversions.
Using conversions is simpler for systems of changing density, while either form can be used for systems of
constant density.
Example 5.1
Example 5.2
Example 5.3
5. IDEAL REACTORS FOR A SINGLE REACTION
III. STEADY STATE PLUG FLOW REACTOR:
In a plug flow reactor the composition of the fluid varies from point to point along a flow path; consequently, the
material balance for a reaction component must be made for a differential element of volume dV. Thus for reactant
A, eqn. 4.1 becomes
-------- (5.10)
5. IDEAL REACTORS FOR A SINGLE REACTION
From figure 5.5, for a elemental volume dV we can write
Input of A, moles/time = FA
Output of A, moles/time = FA + dFA
Disappearance of A by reaction, moles/time = (–rA)V
Introducing the above evaluated terms in eqn. 3.10 we get:
Also,
For the reactor as a whole the expression must be integrated. Therefore, we obtain
5. IDEAL REACTORS FOR A SINGLE REACTION
Thus,
Equation 5.17 allows the determination of reactor size for a given feed rate and required conversion.
Compare Eqns. 5.11 and 5.17. The difference is that in plug flow rA varies, whereas in mixed flow rA is constant.
As a more general expression for plug flow reactors, if the feed on which conversion is based, subscript 0, enters
the reactor partially converted, subscript i, and leaves at a conversion designated by subscript f, we have
or -------- (5.18)
5. IDEAL REACTORS FOR A SINGLE REACTION
For the special case of constant-density systems XA = 1 – CA/CA0 and dXA = - dCA/CA0 in which case the
performance equation can also be written in terms of concentrations or,
These performance equations, Eqns. 5.17 to 5.19, can be written either in terms of concentrations or conversions.
For systems of changing density it is more convenient to use conversions; however, there is no particular
preference for constant density systems.
Whatever its form, the performance equations interrelate the rate of reaction, the extent of reaction, the reactor
volume, and the feed rate, and if any one of these quantities is unknown it can be found from the other three.
5. IDEAL REACTORS FOR A SINGLE REACTION
Figure 5.6 displays these performance equations and shows that the spacetime needed for any particular duty can
always be found by numerical or graphical integration. However, for certain simple kinetic forms analytic
integration is possible-and convenient. To do this, insert the kinetic expression for r, in Eq. 17 and integrate.
5. IDEAL REACTORS FOR A SINGLE REACTION
Some of the simpler integrated forms for plug flow are as follows:
a) Zero-order homogeneous reaction, any constant ƐA,
-------- (5.20)
-------- (5.21)
c) First-order reversible reaction, CR0/CA0 = M, kinetics approximated or fitted by –rA = k1CA – k2CR with an
observed equilibrium conversion XAe, any constant ƐA,
-------- (5.22)
d) Second-order irreversible reaction, A + B → products with equimolar feed or 2A → products, any constant ƐA,
-------- (5.23)
5. IDEAL REACTORS FOR A SINGLE REACTION
Where the density is constant, put ƐA = 0 to obtain the simplified performance equation.
By comparing the batch expressions of derived earlier with these plug flow expressions we find:
1. For systems of constant density (constant-volume batch and constant-density plug flow) the performance equations are
identical, τ for plug flow is equivalent to t for the batch reactor, and the equations can be used interchangeably.
2. For systems of changing density there is no direct correspondence between the batch and the plug flow equations and the
correct equation must be used for each particular situation. In this case the performance equations cannot be used
interchangeably.
Example 5.4
Example 5.5
Example 5.6
5. IDEAL REACTORS FOR A SINGLE REACTION
Holding time and Space Time for Flow Reactors:
They are defined as:
Space time:
Holding Time:
-------- (5.24)
For constant density systems (all liquids and constant density gases)
For changing density systems τ ≠ t and t ≠ V/v0 in which case it becomes difficult to find how these terms are
related.
5. IDEAL REACTORS FOR A SINGLE REACTION
As a simple illustration of the difference between these 2, consider two cases of the steady-flow popcorn popper of
which takes in 1 liter/min of raw corn and produces 28 liters/min of product popcorn.
Consider three cases, called X, Y, and Z, which are shown in Fig. 5.7.
5. IDEAL REACTORS FOR A SINGLE REACTION
In the first case (case X) all the popping occurs at the back end of the reactor.
In the second case (case Y) all the popping occurs at the front end of the reactor.
In the third case (case Z) the popping occurs somewhere between entrance and exit.
In all three cases, irrespective f where the popping occurs –
However, we see that the residence time in the three cases is very different, or
5. IDEAL REACTORS FOR A SINGLE REACTION
This shows that, the value of t depends on what happens in the reactor, while the value of τ is independent of what
happens in the reactor.
This example shows that t and τ are not, in general, identical.
For batch systems, the proper performance measure is the time of reaction.
However, for flow systems, τ or V/FA0 is the proper performance measure.
The above simple example shows that in the special case of constant fluid density the space-time is equivalent to
the holding time; hence, these terms can be used interchangeably.
This special case includes practically all liquid phase reactions.
However, for fluids of changing density, e.g., non-isothermal gas reactions or gas reactions with changing number
of moles, a distinction should be made between t and τ and the correct measure should be used in each situation.
5. IDEAL REACTORS FOR A SINGLE REACTION
5. IDEAL REACTORS FOR A SINGLE REACTION
6. DESIGN FOR SINGLE REACTIONS
With the wide choice of systems available and with the many factors to be considered, no neat formula can be
expected to give the optimum setup.
Experience, engineering judgment, and a sound knowledge of the characteristics of the various reactor systems are
all needed in selecting a reasonably good design and, the choice in the last analysis will be dictated by the
economics of the overall process.
Single reactions are reactions whose progress can be described and followed adequately by using one and only one
rate expression coupled with the necessary stoichiometric and equilibrium expressions.
For such reactions product distribution is fixed; hence, the important factor in comparing designs is the reactor
size.
A. Size Comparison of Single Reactors:
1. Batch Reactor:
The batch reactor has the advantage of small instrumentation cost and flexibility of operation (may be shut down
easily and quickly).
6. DESIGN FOR SINGLE REACTIONS
1. Batch Reactor:
It has the disadvantage of high labor and handling cost, often considerable shutdown time to empty, clean out, and
refill, and poorer quality control of the product.
Hence the batch reactor is well suited to produce small amounts of material and to produce many different
products from one piece of equipment.
Regarding reactor sizes, on comparing eqns. 5.4 – 5.19 for a given duty and for ε = 0 shows that an element of
fluid reacts for the same length of time in the batch and in the plug flow reactor.
Thus, the same volume of these reactors is needed to do a given job.
On a long-term production basis one must correct the size requirement estimate to account for the shutdown time
between batches.
Still, it is easy to relate the performance capabilities of the batch reactor with the plug flow reactor.
6. DESIGN FOR SINGLE REACTIONS
2. Mixed Versus Plug Flow Reactors, First and Second Order Reactions:
For a given duty the ratio of sizes of mixed and plug flow reactors will depend on the extent of reaction, the
stoichiometry, and the form of the rate equation.
For the general case, a comparison of eqns. 5.11 and 5.17 will give this size ratio. Consider this comparison for the
large class of reactions approximated by the simple nth-order rate law:
Where n varies anywhere from 0 to 3. For mixed flow, eqn. 5.11 gives
-------- (6.1)
Equations 1 and 2 are displayed in graphical form in Fig. 6.1 to provide a quick comparison of the performance of
plug flow with mixed flow reactors.
5. IDEAL REACTORS FOR A SINGLE REACTION
6. DESIGN FOR SINGLE REACTIONS
For
The ordinate becomes the volume ratio Vm,/Vp or space-time ratio τm/τp if the same quantities of identical feed are
used.
For identical feed composition CA0 and flow rate FA0, the ordinate of this figure gives directly the volume ratio
required for any specified conversion. Figure 6.1 shows the following.
1. For any particular duty and for all positive reaction orders the mixed reactor is always larger than the plug flow reactor. The
ratio of volumes increases with reaction order.
2. When conversion is small, the reactor performance is only slightly affected by flow type. The performance ratio increases
very rapidly at high conversion; consequently, a proper representation of the flow becomes very important in this range of
conversion.
3. Density variation during reaction affects design; however, it is normally of secondary importance compared to the difference
in flow type.
6. DESIGN FOR SINGLE REACTIONS
3. Variation of Reactant Ratio for Second Order Reactions:
Second order reactions of two components and of the type:
-------- (3.13)
Can behave as second order reactions of one component when the reactant ratio is unity. Therefore,
-------- (6.3)
For the second case, when a large excess of reactant B is used then its concentration does not change appreciably
(CB ≈ CB0) and the reaction approaches first-order behavior with respect to the limiting component A, or
-------- (6.4)
Thus in Fig. 6.1, and in terms of the limiting component A, the size ratio of mixed to plug flow reactors is
represented by the region between the first-order and the second-order curves.
6. DESIGN FOR SINGLE REACTIONS
4. General Graphical Comparison:
For reactions with arbitrary but known rate the
performance capabilities of mixed and plug flow
reactors are best illustrated in Fig. 6.2.
The ratio of shaded and of hatched areas gives the
ratio of space-times needed in these two reactors.
From Fig. 6.2, for all nth-order reactions, n > 0, it can
be seen that mixed flow always needs a larger volume
than does plug flow for any given duty.