Dynamic Simulation of a Non-isochoric & Non-isothermal CSTR

A simple dynamic simulation of a CSTR with varying volume and temperature conditions, solved both in the t & s domain (with simple controls)

A CSTR with two inlet streams and one outlet stream. A reaction of Acetic Acid and Ethanol to produce Ethyl acetate is occuring in the CSTR. It provides heat due to its reaction. 5 state equations were developed, with 4 input variables and 5 output variables. The t-domain solution is in the m files: 'cstr1.m' and 'plot_cstr1.m'. The 5 state equations were linearized using Jacobian function in the cstr1_linearized.m file, with an ode15s evaluation. The s-domain solution, where the parameters of the transfer functions come from the z matrix of coefficients from the cstr1_linearized.m file, is given in the simulink file 'cstr1_s_block.slx'. A simple feedback and a simple cascade control scheme is implemented in the system using the Ziegler-Nichols Manual Tuning Rules, found in 'cstr1_s_block_control.slx file'. Note that the transfer functions in both simulink files are just the same. The one in 'cstr1_s_block_control.slx' is only sub-grouped together. A documentation of the whole coding and illustration of the problem compliments the m and slx files.

This is now updated which includes 2 CSTRs in series.