Question

Problem 4: Consider the following system differential equations for a plug flow reactor ar,=-2k,C,'; F,(0) = 2 dV dV dV where k 0.25, k2 0.05 and Cr-1. If you were to solve this system of differential equation using the Euler Method determine the values of F, FB, F for the first two iterations of the procedure. Use a sample period of 0

PLEASE USE MATLAB! Thank you in advance!

Answer 1

%%% Matlab code %%%

clc;
close all;
clear all;
Fa(1)=2;
Fb(1)=0;
Fc(1)=6;
k1=0.25;
k2=0.05;
Ct=1;
h=0.1;

%%%% Eular method
for n=1:2
Ca=Fa(n)/(Fa(n)+Fb(n)+Fc(n))*Ct;
Cb=Fb(n)/(Fa(n)+Fb(n)+Fc(n))*Ct;
Cc=Fc(n)/(Fa(n)+Fb(n)+Fc(n))*Ct;
Fa(n+1)=Fa(n)-h*2*k1*(Ca)^2;
Fb(n+1)=Fb(n)+2*k1*(Ca)^2;
Fc(n+1)=Fc(n)-k2*Cb*Cc;
fprintf('After iteration %d : Fa= %f Fb= %f Fc=%f \n ',n,Fa(n+1),Fb(n+1),Fc(n+1));
  
  
  
end

OUTPUT:

After iteration 1 : Fa= 1.996875 Fb= 0.031250 Fc=6.000000
After iteration 2 : Fa= 1.993782 Fb= 0.062185 Fc=5.999855
>>