Question 19.8: A batch reactor converts component A into B, which in turn d...

A batch reactor converts component A into B, which in turn decomposes into C :

\stackrel{k_{1}}{ A \rightarrow B } \stackrel{k_{2}}{\rightarrow C }

where k_{1}=k_{10} e^{-E_{1} / R T} and k_{2}=k_{20} e^{-E_{2} / R T}.

The concentrations of A and B are denoted by x_{1} and x_{2}, respectively. The reactor model is

\begin{aligned}&\frac{d x_{1}}{d t}=-k_{10} x_{1} e^{-E_{1} / R T} \\&\frac{d x_{2}}{d t}=k_{10} x_{1} e^{-E_{1} / R T}-k_{20} x_{2} e^{-E_{2} / R T}\end{aligned}

Thus, the ultimate values of x_{1} and x_{2} depend on the reactor temperature as a function of time. For

\begin{array}{ll}k_{10}=1.335 \times 10^{10} min ^{-1}, & k_{20}=1.149 \times 10^{17} min ^{-1} \\E_{1}=75,000 J / g mol , & E_{2}=125,000 J / g mol \\R=8.31 J /( g mol K ) & x_{10}=0.7 mol / L , \quad x_{20}=0\end{array}

Find the constant temperature that maximizes the amount of B, for 0 \leq t \leq 6  min.

The Blue Check Mark means that this solution has been answered and checked by an expert. This guarantees that the final answer is accurate.
Learn more on how we answer questions.

The reactor equations are:

\frac{d x_{1}}{d t}=-k_{1} x_{1}       (1)

\frac{d x_{2}}{d t}=k_{1} x_{1}-k_{2} x_{2}      (2)

Where k_{1}=1.335 * 10^{10} e^{-75000 /(8.31 * T)} ; k_{2}=1.149 * 10^{17} e^{-125000 /(8.31 * T)}

By using MATLAB, this differential equation system can be solved using the command “ode45”. Furthermore, we need to apply the command “fminsearch” ino order to optimize the temperature. In doing so, the results are:

T_{o p}=360.92 K ; x_{2, \max }=0.343

The 'Blue Check Mark' means that either the MATLAB code/script/answer provided in the answer section has been tested by our team of experts; or the answer in general has be fact checked.

Learn more on how do we answer questions.

%% Exercise 19.8
function y = Exercise_19_8(T)
k10 = 1.335*10^10; % min^(-1)
k20 = 1.149*10^17; % min^(-1)
E1 = 75000; % J/(g.mol)
E2 = 125000; % J/(g.mol)
R = 8.31; % J/(g.mol.K)
x10 = 0.7; % mol/L
x20 = 0; % mol/L
k1 =k10*exp(-E1/(R*T));
k2 = k20*exp(-E2/(R*T));
time = [0,6]; % Time period;
initial_val = [x10, x20];
options = odeset(‘RelTol’,1e-4,’AbsTol’,[1e-4 1e-4]);
[~,X] = ode45(@reactor, time,initial_val, options);
y = -X(end,2); % Because of fminsearch, has to be opposite
function dx = reactor(t,x)
dx = zeros(2,1); % A column vector
dx(1) = -k1*x(1);
dx(2) = k1*x(1)-k2*x(2);
end
end
%% Exercise 22.10 main
clear all;clc; close all;
T_range = [200, 500];
T = fminsearch(@Exercise_19_8, 200);
x2_max =-Exercise_19_8(T);

Related Answered Questions