Question 9.5.3: Energy Balance on an Adiabatic Reactor The dehydrogenation o......

Energy Balance on an Adiabatic Reactor
The dehydrogenation of ethanol to form acetaldehyde

C_{2}H_{5}OH(v) → CH_{3}CHO(v) + H_{2}(g)

is carried out in a continuous adiabatic reactor. Ethanol vapor is fed to the reactor at 400°C, and a conversion of 30% is obtained. Calculate the product temperature (a) using heat capacity formulas from Table B.2, and (b) using APEx.

Step-by-Step
The 'Blue Check Mark' means that this solution was answered by an expert.
Learn more on how do we answer questions.

Basis: 100 mol Feed
Material balances lead to the flowchart shown here.

Since only one reaction occurs, we could equally well choose the reactants and products [C_{2}H_{5}OH(v) , CH_{3}CHO(v) , H_{2}(g)] or their elemental constituents [C(s), H_{2}(g), O_{2}(g)] as references for enthalpy calculations. Let us choose the molecular species.

The open system energy balance neglecting kinetic and potential energy changes and shaft work and setting Q = 0 for this adiabatic reactor is

\Delta H= \xi \Delta H^{°}_{r} + \sum{n_{out}\hat{H}_{out}} – \sum {n_{in}\hat{H}_{in}} =0

\begin{matrix}\Delta \dot{H}=\xi \Delta \dot{H}^{°}_{r} + \sum{\dot{n}_{out}\hat{H}_{out}} – \sum {\dot{n}_{in}\hat{H}_{in}}    (single reaction)    \pmb{(9.5-1a)}\end{matrix}

The expression for ΔH is that of Equation 9.5-1a, with the dots above H and n omitted since a molar amount and not a flow rate was chosen as a basis of calculation.
(a) Manual Calculation
Extent of Reaction:
We could use any reactant or product as a basis for calculating ξ. Let us use acetaldehyde. From Equation 9.1-3,

\xi = \left. \Large{\{} \right. \begin{matrix} Batch: \frac{(n_{i,out} – n_{i,in})(mol  i)}{\nu _{i}(mol  i)} = \frac{n_{ir}}{\left|\nu _{i}\right|} \\Continuous: \frac{(\dot{n}_{i,out} – \dot{n}_{i,in})(mol  i/ time)}{\dot{\nu} _{i}(mol   i/ time)} = \frac{\dot{n}_{ir}}{\left|\dot{\nu} _{i}\right|} \end{matrix}                 (9.1-3)

\xi = \frac{\left|(n_{CH_{3}CHO})_{out}  –  (n_{CH_{3}CHO})_{in} \right| }{\left|\nu_{CH_{3}CHO} \right| } = \frac{\left|30.0  mol –  0  mol\right| }{1   mol} =30.0

Standard Heat of Reaction:
From Equation 9.3-1 and Table B.1 (heats of formation),

\boxed{ΔH^{°}_{r} = \sum\limits_{i}{\nu _{i}\Delta \hat{H}^{°}_{fi}} = \sum\limits_{products}{\left|\nu _{i}\right| \Delta \hat{H}^{°}_{fi}}- \sum\limits_{reactants}{\left|\nu _{i}\right| \Delta \hat{H}^{°}_{fi}}}              (9.3-1)
\begin{matrix}ΔH^{°}_{r} = \sum{\nu _{i}\Delta \hat{H}^{°}_{f}} = (-1  mol  C_{2}H_{5}OH)(\Delta \hat{H}^{°}_{f})_{C_{2}H_{5}OH(v)} + (1  mol CH_{3}CHO)(\Delta \hat{H}^{°}_{f})_{CH_{3}CHO(v)} \\ + ( 1  mol  H_{2})(\Delta \hat{H}^{°}_{f})_{H_{2}(g)} \\= [(-1)(-235.31)+(1)(- 166.2)+(1)(0)]  kJ = 69.11  kJ\end{matrix}

Inlet Enthalpy:

\hat{H}_{1} = \int_{25°C}^{400°C}{(C_{p})_{C_{2}H_{5}OH}} \overset{C_{p}  from  Table  B.2 }{\left. \Large{\Longrightarrow } \right.} \hat{H}_{1} = 33.79  kJ/mol

Outlet Enthalpies:
The heat capacities of ethanol vapor and hydrogen are given in Table B.2. For acetaldehyde vapor, the heat capacity is given by Poling, Prausnitz, and O’Connell^{6} :

\begin{matrix}(C_{p})_{CH_{3}CHO} \left(\frac{kJ}{mol\cdot °C} \right) = 0.05048 + 1.326 \times 10^{-4}T  –  8.050 \times 10^{-8}T^{2} + 2.380 \times 10^{-11}T^{3}\end{matrix}

where T is in °C. For the three species in the product stream,

\hat{H}_{i} = \int_{25°C}^{T_{ad}}C_{pi}(T)  dT,    i=1,2,3

If the heat capacity formulas for the three species are substituted in this expression and the integrals are evaluated, the results are three fourth-order polynomial expressions for \hat{H}_{2}(T_{ad}) , \hat{H}_{3}(T_{ad}) , and \hat{H}_{4}(T_{ad}) :

\begin{matrix} \pmb{C_{2}H_{5}OH}:    \hat{H}_{2}(kJ/mol) = 4.958 \times 10^{-12}T^{4}_{ad}  –  2.916 \times 10^{-8}T^{3}_{ad}\\ +7.860 \times 10^{-5}T^{2}_{ad}+ 0.06134T_{ad}  –  1.582\\ \pmb{CH_{3}CHO:}   \hat{H}_{3}(kJ/mol) = 5.950 \times 10^{-12}T^{4}_{ad}  –  2.683 \times 10^{-8}T^{3}_{ad} \\ +6.630 \times 10^{-5}T^{2}_{ad} + 0.05048T_{ad}  –  1.303 \\  \pmb{H_{2}:}    \hat{H}_{4}(kJ/mol) = – 0.2175 \times 10^{-12}T^{4}_{ad}  + 0.1096 \times 10^{-8}T^{3}_{ad} \\ +0.003825 \times 10^{-5}T^{2}_{ad} + 0.02884T_{ad}  – 0.7210\end{matrix}

Solve the Energy Balance for T_{ad}

\begin{matrix}\Delta H= \xi \Delta H^{°}_{r} + (70.0  mol)\hat{H}_{2} + (30.0  mol)\hat{H}_{3} + (30.0  mol)\hat{H}_{4}  –  (100.0  mol)\hat{H}_{1} =0 \\ \left. \Large{\Downarrow} \right. Substitute  for  \xi (=30.0), \Delta H^{°}_{r} (=69.11  kJ),\hat{H}_{1}(=33.79  kJ/mol), and  \hat{H}_{2}  through  \hat{H}_{4} \\ \Delta H= 5.190 \times 10^{-10} T^{4}_{ad}+ 2.813\times 10^{-6} T^{3}_{ad} +7.492 \times 10^{-3} T^{2}_{ad} +6.673 T_{ad}  –  1477 =0\end{matrix}

This equation can be solved using an equation-solving program or a spreadsheet.^{7} The solution is

\boxed{T_{ad} = 185°C}

(b) Spreadsheet Calculation. Replicate the following solution for yourself in Excel with APEx to make sure you understand it. Each step matches a calculation in Part (a). The numbers and formulas entered in Column B—including an initially guessed temperature of 100°C—are shown below the spreadsheet.

[B3] 30
[B4] =DeltaHfg(“acetaldehyde”) – DeltaHfg(“ethanol”)
[B5] =Enthalpy(“ethanol”,25,400,,“g”)
[B6] 100
[B7] =Enthalpy(“ethanol”,25,B$6,,“g”)
[B8] =0.05048*(B$6 – 25) + 1.326e – 4*(B$6^2 – 25^2)/2

– 8.05e-8*(B$6^3 – 25^3)/3 + 2.38.e-11*(B$6^4 – 25^4)/4

[B9] =Enthalpy(“hydrogen”,25,B$6)
[B10] = B3*B4 + 70*B7 + 30*B8 + 30*B9 –  100*B5
Once the cell entries have been made, use Excel’s Solver or Goal Seek to set the value in Cell B10 (which equals ΔH) equal to 0 by varying the value in Cell B6 (the temperature). When you click “Solve” or “OK,” the value in Cell B6 changes to 185, so that the solution is \boxed{185°C} and the value in Cell B10 changes to – 2.4 × 10^{-7} , more than close enough to zero for our purposes.
Contrast the ease of this calculation with the cumbersome set of calculations in Part (a). This example would be even easier if the heat capacity formula for acetaldehyde were listed in Table B.2 so that the Enthalpy function of APEx could be used to determine \hat{H}_{3}. Since it isn’t, the formula for the integral of the given formula must be entered manually in [B8].

^{6}B. E. Poling, J. M. Prausnitz, and J. P. O’Connell, The Properties of Gases and Liquids, 5th Edition, McGraw-Hill, New York, 2001. The formula given has been derived from the one shown in this reference, which is for the heat capacity in J/(mol·K) with the temperature expressed in kelvin

^{7}To obtain the solution using a spreadsheet, put a guessed value of T_{ad} in one cell and the expression for ΔH in an adjacent cell, and use the Goal Seek tool to determine the value of T_{ad} for which the expression for ΔH equals zero. A first guess might be the value of T_{ad} obtained by dropping all higher-order terms in the expression, leaving 6.673T_{ad}  –  1134 = 0 \Longrightarrow T_{ad}\approx 170°C.

References: C_{2}H_{5}OH (v), CH_{3}CHO (v), H_{2} (g) at 25°C and 1 atm

Substance n_{in}

(mol)

\hat{H}_{in}

(kJ/mol)

n_{out}

(mol)

\hat{H}_{out}

(kJ/mol)

C_{4}H_{5}OH 100.0 \hat{H}_{1} 70.0 \hat{H}_{2}
CH_{3}CHO 30.0 \hat{H}_{3}
H_{2} 30.0 \hat{H}_{4}
A B
1 Example 9.5-3
2
3 Xi 30
4 DHr 69.11
5 H1 33.79029
6 T 100
7 H2(T) 5.309161
8 H3(T) 4.381741
9 H4(T) 2.164416
10 DH -737.703
2

Related Answered Questions