Question

Component A reacts irreversibly and at a specific reaction rate k to form product, component B via the first order reaction A as shown in figure 2 B in a three CSTRs connected in series Let the concentration of component A in the inflowing feed stream be Co. (kmoles of A per unit volume) and in the reactor CA. The resulting govering equations are given below decat (C0CH) -RCHI den (C1-C12) -kCz2 des C12-602) – kad Given that at t 0 C 1 = 0.4 kmole/mº, CA2 = 0.2 kmole/mº. C13 = 0.1 kmole/ m, and Co = 1.8 kmole/m". The value of reaction rate constant can be taken as k=0.5, and t = 2 min. Let C(1) Cai, C(2)** CA2, C(3) Caz. Write a matlab code to find the following: a. The outlet concentration of components Car, CA, and Cx from the reactor and plot the concentration after 3 minutes. b. Generate algebraic function of the ode using deval c. At what time the concentration of component Cxo is consumed by half. d. Study the effect of varying the ratio of inlet on the reactor output, from equimolar to 2:1 of component C to A. plot the outlet concentration ys to A ratio,

by using MATLAB answer is Code

Answer 1

clc

clear all

close all

format long

CA0=1.8;

CA1=0.4;

CA2=0.2;

CA3=0.1;

k=0.5;

tau=2;

f=@(t,y) [1/tau*(CA0-y(1))-k*y(1);1/tau*(y(1)-y(2))-k*y(2);1/tau*(y(2)-y(3))-k*y(3)];

[T,C]=ode45(f,[0,3],[CA1,CA2,CA3]);

plot(T,C);

legend('CA1(T)','CA2(T)','CA3(T)')