CSTR with Cooled Jacket (ODEs System solved by MATLAB)
#CSTR #reactor #chemical #engineering
clc;
clear;
% System Parameters
V = 1.0; % Reactor volume (m^3)
F = 0.1; % Inlet flow rate (m^3/s)
CA_in = 1.0; % Inlet concentration of A (mol/m^3)
T_in = 300; % Inlet temperature (K)
U = 500; % Heat transfer coefficient (W/m^2-K)
A = 10; % Heat transfer area (m^2)
Cp = 4.18; % Heat capacity (kJ/kg-K)
rho = 1000; % Liquid density (kg/m^3)
k0 = 1.0e3; % Pre-exponential reaction rate constant (1/s)
Ea = 20000; % Activation energy (J/mol)
R = 8.314; % Universal gas constant (J/mol-K)
% Reaction Parameters
Hr = -50000; % Reaction enthalpy (J/mol)
% Cooling Jacket
T_c = 290; % Cooling jacket temperature (K)
