Project -

For this project, you are required to simulate the CSTR described in example #6 of chapter 10 in the course textbook (page 177), but with the temperature controlled by a cooling jacket. Our goal is to change the cooling fluid  temperature  so  that  the  reactor  temperature  is held  at  38^oC,  i.e., SP =  38^oC. The CSTR has a valve controlling the temperature of the cooling water being fed to a cooling jacket with the following properties: h = 5kJ/(min·m²·^oC), and A = 4.1m². You should use the same values for all other parameters as described in example #6 except the rate constant, k = 200·e^-2000/T min^-1, and T is in K.

Begin by downloading project17.m and project17driver.m "attached". For part 1 of the project, you need to modify project17driver.m so that it uses a PI control algorithm to set the requested cooling fluid temperature, T_r, based  on  the  measured  temperature  from  the  temperature  sensor  (i.e.,  the  third  unknown  in  the simulation).  Hints and suggestions for part 1:

• The set point (SP) for this process is an outflow temperature of 38^oC. The initial conditions should be unchanged from example #6 (and unchanged from those included in project17driver.m).
• I highly recommend using the velocity form for the PI control equation
• You should calculate the error and new CO value each second (i.e., Δt = 1.0 sec.). This means that you need to call the ODE solver (ode45) many times and it should simulate the next second every time it is called.
• You will need to store 2 error values if you use the velocity form of the PI equation.
• Generate at least 2 plots. The first plot should display the temperature and the measured (or sensed) temperature each second. The second plot should display the requested coolant temperature each second.
• I found it easiest to just have the controller output equal to the requested temperature. You may also use a non-dimensional control equation, but it is not required and may be slightly more difficulty.

It is critical that part 1 be successfully completed in order to move on to parts 2 and 3.

For part 2 you need to determine the PI parameters, K_c and τ_I, using any one of the tuning methods covered in this course.  Please show plots of the process response using the parameters you obtained.

For part 3, you are concerned about the accuracy of the two constants in the rate constant, k = 200·e^-2000/T min^-1. Pick one of the two constants and analyze how the tuned PI controller constants change when the rate constant parameter is increased and decreased.  Use the same tuning method as part 2.

Assignment Files -

https://www.dropbox.com/s/6z9vhonc2j6biqr/Assignment%20Files.rar?dl=0