Assignment 2-Solution
Assignment 2-Solution
Assignment 2
Categories of Questions
I. Design equations for batch, CSTR, PFR and other types of reactors
II. Stoichiometric table analysis
III. Combined use of design equations (batch, CSTR and PFR) and stoichiometric table analysis
(most important for this unit)
IV. Space time and space velocity
3
0 F
V V
Question 1 (II)
10 mol% of A in nitrogen (inert) at 100 °C and 2 atm enters an isothermal open tubular reactor operating at 300 °C,
where the following gas phase reaction occurs:
A → 2B
The reactor was designed to convert at least 50% of A entering the reactor. A technician sampling the reactor exit gas
reported ~10 mol% of B in the leaving gas mixture. There is a 10% pressure drop across the reactor. Does the reactor
operate as per the design specification?
Solutions: method 1
Assume F0 mol of gas enters the reactor per unit time.
It includes 0.1F0 mol A and 0.9F0 mol inert gas.
Under steady state, the reactor exit has 0.1F0(1 − xA) mol A, 0.9F0 mol inert gas and 0.2F0xA mol B.
𝐹
In the inlet, A is 10 mol%: = 10%
𝐹
2𝐹 𝑥
In the outlet, B is 10 mol%: = 10%
𝐹 +𝐹 𝑥
𝐹
2 𝑥
𝐹 2(0.1)𝑥
Divide every term by F0: = 10% = 0.1 𝑥 = 0.526
𝐹 1 + 0.1𝑥
1+ 𝑥
𝐹
6
Question 2 (II)
Consider the following reactions in diluted aqueous phase in a batch reactor:
A+ B →R
R+B→S
The feed contains only A and B. Prove the following concentration relations after a finite time into the reactions.
𝑐 𝑐 𝑐
=1− −
𝑐 𝑐 𝑐
𝑐 2(𝑐 −𝑐 ) 𝑐 −𝑐
= +
𝑐 𝑐 𝑐
Are these relations still valid if the reactions are gas phase reactions in a batch reactor?
Solutions:
𝑐 𝑐 𝑐
=1− − 𝑐 =𝑐 −𝑐 −𝑐 𝑐 +𝑐 +𝑐 =𝑐
𝑐 𝑐 𝑐
𝑐 2(𝑐 −𝑐 ) 𝑐 −𝑐
= + 𝑐 = 2(𝑐 −𝑐 )+𝑐 −𝑐 𝑐 + 2𝑐 − 𝑐 = 2𝑐 −𝑐
𝑐 𝑐 𝑐
7
𝑁 + 𝑁 + 𝑁 = (𝑁 − 𝛼) + (𝛼 − 𝛽) + 𝛽 = 𝑁
Division by V0: + + =
Similarly: 𝑁 + 2𝑁 − 𝑁 = (𝛼 − 𝛽) + 2(𝑁 − 𝛼) − (𝑁 − 𝛼 − 𝛽) = 2𝑁 −𝑁
Division by V0 throughout and assuming a constant density reaction system: 𝑐 + 2𝑐 − 𝑐 = 2𝑐 −𝑐
The total number of moles is not conserved, so the concentration relationships will only hold for a constant density
reacting system.
For gas phase reactions in a batch reactor, the density is constant, so the relations still hold.
8
Question 3 (III)
The gas phase reaction A + B → R + 2S is run in a PFR with a feed where B is present in excess of the stoichiometric
requirement (cA0 = 1 mol/L and cB0 = 2 mol/L). What is the maximum concentration of R possible at the reactor outlet?
Question 3 (III)
The gas phase reaction A + B → R + 2S is run in a PFR with a feed where B is present in excess of the stoichiometric
requirement (cA0 = 1 mol/L and cB0 = 2 mol/L). What is the maximum concentration of R possible at the reactor outlet?
𝐹 𝑥 𝑥 1 3 3 1
Molar percentage of R: = = 𝑥 +1
3𝐹 + 𝐹 𝑥 3+𝑥 3 𝑥 𝑥 3
+1 +1
𝑥 𝑥
Maximum molar percentage of R is achieved at 𝑥 = 1. That is, 25%
𝑁 𝑝 𝑚𝑜𝑙 2𝑚𝑜𝑙
𝑝𝑉 = 𝑁𝑅𝑇 𝑐 = = = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡. = 1 + = 3 𝑚𝑜𝑙/𝐿 𝑐 = 25% × 3 = 0.75 mol/L
𝑉 𝑅𝑇 𝐿 𝐿
10
Question 4 (III)
A gas phase reaction A → 2B with k = 0.5 h−1 is carried out in a (batch) piston reactor at constant temperature and
pressure. 1 mole of A and 9 moles of inert are introduced to the reactor at the start of the reaction. What is the
composition of the gas in the reactor after 1 hour?
0 0 0 0 0
( ) 1−𝑥 =𝑒 / 𝑥 = 0.393
0
11
Question 4 (III)
Based on stoichiometric table, gas composition at 1 h:
𝑁 = 1 − 𝑥 = 1 − 0.393 = 0.607 mol
𝑁 = 2𝑥 = 2 × 0.393 = 0.786 mol
𝑁 = 9 mol
1 𝑑𝑁 1 𝑑𝑁 1 𝑑𝑁 𝑁
Design equation for batch reactor: =𝑟 𝑟 = −𝑘𝑐 = −𝑘𝑐 =− 𝑘
𝑉 𝑑𝑡 𝑉 𝑑𝑡 𝑉 𝑑𝑡 𝑉
𝑑𝑁 𝑑𝑁
=− 𝑘𝑁 =− 𝑘𝑑𝑡
𝑑𝑡 𝑁
𝑡
𝑑𝑁 𝑁 𝑁 . ×
Integrate: = −𝑘𝑑𝑡 ln =− 𝑘𝑡 𝑥 =1− =1−𝑒 =𝑒 = 0.393
𝑁 0 𝑁 𝑁
12
Question 5 (III)
At 650 °C, phosphine vapour decomposes as follows:
xA= 0.75
4PH3 → P4(g) + 6H2(g) rPH3 = −kcPH3 FA0 = 10 mol/h
where k = 10 h−1. What size of PFR operating at 649 °C and 11.4 atm is T = 649 ℃
needed for 75% conversion of 10 mol/h of phosphine in a feed with P = 11.4 atm
V=?
phosphine-inert ratio of 2:1?
Analysis:
𝑑𝑥 𝜈
𝜃 −
𝑥
𝑉=𝐹 𝜈
0 −𝑟 𝑐 = 𝑐
1+𝜀 𝑥
First replace 𝑟 by 𝑐 , then replace 𝑐 by 𝑥 .
Solutions:
𝐹 −4 + 1 + 6 2 𝑃 11.4 × 1.013 × 10 2
𝜀 =𝛿 = × = 0.5 𝑐 = = = 150 mol/m 𝑐 = 𝑐 = 100 mol/m
𝐹 4 2+1 𝑅 𝑇 8.314 × (273 + 649) 3
1−𝑥 1−𝑥
𝑐 = 𝑐 = × 100
1+𝜀 𝑥 1 + 0.5𝑥
0.75 0.75
𝑑𝑥 𝑑𝑥 1 + 0.5𝑥
𝑉=𝐹 = 10 = 𝑑𝑥
0 −𝑟 0 10𝑐 0 (1 − 𝑥 ) × 100
Pay attention to the sign
13
0.75 0.75
1 1.5 1 1.5 1 − 0.75 1
= 𝑑𝑥 − 0.5𝑑𝑥 =− ln − (0.5 × 0.75 − 0.5 × 0) = 0.017 m = 17 L
100 0 1−𝑥 100 0 100 1−0 100
14
Question 6 (I)
A stream of aqueous monomer A (1 mol/L, 4 L/min) enters a 2 L CSTR, is radiated therein and polymerises as follows:
+A +A +A
A→ R→ S→ T…
In the exit stream, 𝑐 = 0.01 mol/L and, for a particular reaction product W (a polymer), 𝑐 = 0.0002 mol/L. Find the
rate of reaction of A and the rate of reaction of W. Assume volume change is insignificant.
Solutions:
Applying the CSTR design equation assuming constant density (diluted Fi0 Fi
aqueous solution): V
𝐹 𝑥
Design equation: 𝑉 =
−𝑟
We know V, FA0 and xA so we can find rA:
𝐹 𝑥 (1 × 4) × (1 − 0.01)
𝑟 =− =− =− 1.98 mol L min
𝑉 2
For product, we use mole balance equation: 𝐹 +𝑟 𝑉=𝐹
𝐹 −𝐹 𝑣 (𝑐 − 𝑐 ) 4 × (0.0002 − 0)
𝑟 = = = = 0.0004 mol L min
𝑉 𝑉 2
15
Question 7 (I)
Consider the following variation to PFR where the feed to the reactor
is distributed uniformly along the reactor length (at α mol h −1
m3reactor). Derive the design equation from the general mole balance
equation.
Hint: Apply the mole balance equation to ΔV, a differential element
of the reactor. Integration of the differential form of the mole balance
equation yields the design equation.
Solutions:
Inflow + Rate of generation = Outflow + Rate of accumulation
Consider the mole balance of A (or any species) over a differential element of volume ΔV:
𝐹 | + 𝛼∆ 𝑉 + 𝑟 ∆𝑉 = 𝐹 | ∆
𝑑𝑓
Taylor’s Expansion: 𝑓(𝑥 + ∆𝑥) = 𝑓(𝑥) + ∆𝑥
𝑑𝑥
𝑑 𝑑𝐹
Substitute 𝐹 | ∆ =𝐹 | + ∆𝑉 =𝑟 +𝛼
𝑑𝑉 𝑑𝑉
𝑑𝐹
𝑉=
0 𝑟 +𝛼
16
Question 8 (I)
.
An irreversible aqueous phase reaction A → products with the rate law 𝑟 =− (with −rA in mol L−1 min−1 and
1+0.6
cA in mol/L) is run in a PFR. A is delivered to the reactor at a concentration of 1.2 mol/L and a volumetric flow rate of
0.7 L/min. Plot the concentration profile along the length of reactor and, from which, determine the volume required for
at least 95% conversion of the feed.
Solutions: This equation is best solved with a computational software such as MATLAB or
Polymath. The following is the concentration profile obtained using Polymath:
Use PFR design equation:
1.5
= 𝑟 =− (1)
𝑑𝑉 1+0.6
𝐹 =𝑐 𝑣 (2)
Question 9 (I)
For a liquid phase reaction between A and B in a semi-batch reactor (where A is initially present in the reactor while B
is fed through the inlet), the design equations describing the variations of concentrations of A and B have been derived
during the lectures. Derive the corresponding design equation governing the concentration of product C. Write down
the appropriate initial conditions.
Solutions:
A+B→C
Mole balance for C (no inflow and outflow, well-mixed reactor contents):
𝑑 𝑑( 𝑉) 𝑑𝑉 𝑑
𝑟 𝑉= = =𝑐 +𝑉 (1)
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
𝑑𝑉 v0, cB0
=𝑣 (2)
𝑑𝑡 V, cA, cB, cC
𝑑𝑐
Combining (1) and (2): 𝑟 𝑉=𝑣 𝑐 +𝑉
𝑑𝑡
𝑑𝑐 𝑣 𝑐 𝑣 𝑐
=𝑟 − =𝑟 − (𝑟 =− 𝑟 =− 𝑟 )
𝑑𝑡 𝑉 𝑉 +𝑣 𝑡
Question 10 (II)
Steam reforming of methane (SRM) is a high temperature (700-1000 °C) catalytic reaction for producing hydrogen by
splitting water. The product is a mixture of CO and H2 generally known as the synthesis gas:
CH4 + H2O → CO + 3H2
The water-gas shift (WGS) reaction also occurs under the steam reforming conditions:
CO + H2O → CO2 + H2
With the help of a stoichiometric table, discuss the effect of WGS on hydrogen production for the following two cases:
(a) Methane as the limiting reactant while water is in large excess in the SRM reaction
(b) Water as the limiting reactant in the SRM reaction
With the help of a stoichiometric table, show how the concentrations of all reaction species may be calculated from the
process variables (temperature, pressure and feed composition) in terms of the conversions of methane and water.
19
Question 11 (II)
Identify the correct statement below after carrying out a stoichiometric table analysis of the following methane coupling
reactions in a flow reactor:
1
2CH4 + O2 ⇌ C2H6 + H2O (1)
2
CH4 + 2O2 → CO2 + 2H2O (2)
(a) This is always an expanding reaction system.
(b) This is always a contracting reaction system.
(c) This is an expanding reaction system if the molar extent of reaction (1) is greater than the molar extent of reaction (2).
(d) This is a contracting reaction system if the molar extent of reaction (2) is greater than the molar extent of reaction (1).
Solutions: Species In Out
Setup the stoichiometric table. Let α and β be the net molar A (CH4) FA0 FA0 – 2α – β
extents of reactions (1) and (2) respectively.
B (O2) FB0 FB0 – 0.5α – 2β
C (C2H6) FC0 FC0 + α
Since α > 0, F < F0 always.
D (H2O) FD0 FD0 + α + 2β
E (CO2) FE0 FE0 + β
In this particular case, the stoichiometric table is not needed.
Reaction 1 is contracting while reaction 2 is volume-conserving. Inert FI0 FI0
Overall it must be contracting. Total F0 F0 – 0.5α
22
Question 12 (III)
A 20% ozone-80% air mixture at 1.5 atm and 93 °C is delivered to a PFR at 1 L/s. Under these conditions, ozone
decomposes by the homogeneous reaction:
2O3 → 3O2
𝑅 = 𝑘𝑐 where k = 0.025 L mol−1 s−1
Analysis:
Question 13 (III)
The gas phase reaction A → B + C is run at constant temperature and pressure in a PFR. The feed is pure A at cA0 = 1
mol/L and the concentration of B in the reactor exit gas is cB = 0.45 mol/L. What is the composition of the reactor exit
gas?
Solutions: 𝑐 =
(1− )
(1) Species In Out
1+
Question 14 (III)
Consider an isothermal single-phase reactor operating at steady state and constant pressure. Given the reaction A + 2B
⇌ R and a gaseous feed of cA0 = cB0 = 100, cB = 20. Find xA and xB.
Solutions: Species In ∆ Out
Setup the stoichiometric table on the right.
A FA0 −FA0xA FA0 − FA0xA
(Batch reactor or flow reactor?)
𝜈 B FB0 −2FA0xA FB0 − 2FA0xA
𝜃 − 𝑥
𝜈
Note: 𝑐 = 𝑐 R 0 FA0xA FA0xA
1+𝜀 𝑥
Total F0 = FA0 + FB0 F = FA0 + FB0 − 2FA0xA
𝐹 𝐹 𝜈
𝜃 = =1 𝜃 = =1 =2
𝐹 𝐹 𝜈
𝐹 −1 − 2 + 1 100
𝜀 =𝛿 = × =− 1
𝐹 1 100 + 100
1−𝑥
For A: 𝑐 = 𝑐 =𝑐 = 100
1−𝑥
1 − 2𝑥
For B: 𝑐 = 𝑐 = 20 ⇒ 𝑥 = 0.444
1−𝑥
𝐹 𝐹 − 2𝐹 𝑥 100 − 2 × 100 × 0.444
𝑥 =1− =1− =1− = 0.888
𝐹 𝐹 100
26
Question 15 (IV)
Consider a gas-phase reaction 2A → R + 2S with unknown kinetics. If a space velocity of l min−1 (based on the inlet
flow rate at process conditions) is needed for 90% conversion of A in a PFR, find the corresponding space time and
mean residence time of fluid in the PFR. What if this is a liquid-phase reaction?
Analysis:
𝑉
Space time: 𝜏 = time necessary to process one reactor volume of fluid based on the inlet conditions
𝑣
𝑣 1
Space velocity: 𝑆𝑉 = if v0 is used to calculate SV, the value obtained is equal to
𝑉 𝜏
𝑉 𝑑
Mean residence time: actual time spent in the reactor For CSTR: 𝑅𝑇 = For PFR: 𝑅𝑇 = 𝐹 ∫0 )
𝑣(−
Solutions:
1
𝜏= = 1 min
𝑆𝑉
Analysis:
Based on the kinetic data obtained in constant volume batch reactor, we can obtain rate law. Since rate law is specific to
nature of reaction, concentration and temperature, but not related to the type of reactor, we can use the same rate law to
solve the design equation to get V for CSTR.
𝐹 𝑥
𝑉=
−𝑟
Solutions:
Kinetic analysis suggests it is a first
order reaction:
𝑟 =− 0.0095𝑐
28
1−𝑥
𝑐 = 𝑐 = (1 − 𝑥 ) × 32.67 × 0.8
1+𝜀 𝑥
20% inert in feed
Data from Excel
100
𝐹 𝑥 𝐹 𝑥 × 0.95
𝑉= = = 3600 = 2.13 m
−𝑟 0.0095𝑐 0.0095 × (1 − 0.95) × 32.67 × 0.8