What is the plan? In the following pages we manipulate Equation (11-9) in order to apply it to each of the reactor types we have been discussing: batch, PFR, PBR, and CSTR. The end result of the application of the energy balance to each type of reactor is shown in Table 11-1. These equations can be used in Step 5 of the algorithm discussed in Example E11-1. The equations in Table 11-1 relate temperature to conversion and to molar flow rates and to the system parameters, such as the overall heat-transfer coefficient and area, Ua, with the corresponding ambient temperature, Ta, and the heat of reaction, ΔHRx.
Example 11-1. Energy Balances of Common Reactors
For an exothermic reaction (–ΔHRx) > 0
In general most of the PFR and PBR energy balances can be written as
3A. PFR in terms of conversion
3B. PBR in terms of conversion
3C. PBR in terms of molar flow rates
3D. PFR in terms of molar flow rates
The equations in Table 11-1 are the ones we will use to solve reaction engineering problems with heat effects.
[Nomenclature: U = overall heat-transfer coefficient, (J/m2 · s · K); A = CSTR heat-exchange area, (m2); a = PFR heat-exchange area per volume of reactor, (m2/m3); CPi = mean heat capacity of species i, (J/mol/K); CPc = the heat capacity of the coolant, (J/kg/K), = coolant flow rate, (kg/s); ΔHRx = heat of reaction, (J/mol);
ΔHRxij = heat of reaction wrt species j in reaction i, (J/mol);
= heat added to the reactor, (J/s); and
All other symbols are as defined in Chapters 5 and 6.]
Examples of How to Use Table 11-1. We now couple the energy balance equations in Table 11-1 with the appropriate reactor mole balance, rate law, and stoichiometry algorithm to solve reaction engineering problems with heat effects. For example, recall the rate law for a first-order reaction, Equation (E11-1.5) in Example 11-1.
If the reaction is carried out adiabatically, then we use Equation (T11-1.B) for the reaction in Example 11-1 to obtain
Consequently, we can now obtain –rA as a function of X alone by first choosing X, then calculating T from Equation (T11-1.B), then calculating k from Equation (E11-1.3), and then finally calculating (–rA) from Equation (E11-1.5).
We can use this sequence to prepare a table of (FA0/–rA) as a function of X. We can then proceed to size PFRs and CSTRs. In the absolute worst case scenario, we could use the techniques in Chapter 2 (e.g., Levenspiel plots or the quadrature formulas in Appendix A). However, instead of using a Levenspiel plot, we will most likely use Polymath to solve our coupled differential energy and mole balance equations.
If there is cooling along the length of a PFR, we could then apply Equation (T11-1.D) to this reaction to arrive at two coupled differential equations.
which are easily solved using an ODE solver such as Polymath.
Similarly, for the case of the reaction A → B carried out in a CSTR, we could use Polymath or MATLAB to solve two nonlinear algebraic equations in X and T. These two equations are the combined mole balance
and the application of Equation (T11-1.C), which is rearranged in the form
From these three cases, (1) adiabatic PFR and CSTR, (2) PFR and PBR with heat effects, and (3) CSTR with heat effects, one can see how one couples the energy balances and mole balances. In principle, one could simply use Table 11-1 to apply to different reactors and reaction systems without further discussion. However, understanding the derivation of these equations will greatly facilitate their proper application and evaluation to various reactors and reaction systems. Consequently, we now derive the equations given in Table 11-1.
Why bother to derive the equations in Table 11-1? Because I have found that students can apply these equations much more accurately to solve reaction engineering problems with heat effects if they have gone through the derivation to understand the assumptions and manipulations used in arriving at the equations in Table 11.1.
We will now dissect the molar flow rates and enthalpy terms in Equation (11-9) to arrive at a set of equations we can readily apply to a number of reactor situations.
18.191.211.66