clc
clear all
close all
V= 2; %cm3
Cawe= 5.0; %mol/cm3;
Fwe= 0.01; %cm3/s;
k0= 8.94e12; %1/s;
E= 22500.0; %cal/mol;
Twe= 298; %K;
ro=0.9; %g/cm3;
cp= 1; %cal/gK;
deltaH= -10.0; %cal/mol;
Tc= 305.0; %K;
R= 1.997; %cal/molK;
UA= 0.001; %cal/Ks;
T=300:1:350;
F=Fwe;
Ca=Cawe./(k0.*exp(-E./(R.*T)).*V/F+1);
r=k0.*exp(-E./(R.*T)).*Ca;
Qw=-r.*V.*deltaH;
plot(T,Qw)