Module 7 - Isothermal CSTRs

vandss.m Generate the steady-state input/output curves for the van de Vusse reactor

To use this m-file to generate Figure M7.8, enter the following at the MATLAB command window prompt:

   k1 = 5/6;
   k2 = 5/3;
   k3 = 1/6;
   cafs = 10;
   [fsov,cbs,cas] = vandss(k1,k2,k3,cafs);

to plot Figure M7.8a, enter

   plot(fsov,cas)

to plot Figure M7.8b, enter

   plot(fsov,cbs)

vanmax.m The the maximum of the concentration of component B as a function of F/V for the van de Vusse reactor

Use the command:

  fmin('vanmax',0,2)

ans =

    1.2921

vanvusse.m The function routine to be used by ode45 to integrate the two differential equations which describe the dynamic behavior of the van de Vusse reactor

Use the commands:

  x0 = [3;1.117]

  [t,x] = ode45('vanvusse',[0 9],x0);
% [t,x] = ode45('vanvusse',0,9,x0); if you are using MATLAB 4.x

  plot(t,x)

to find that you are at a steady-state for F/V = 4/7. Modify the function file if you wish to perform step changes, for example.