Docstoc

_ Plotting Contour Arrhenius Eqn

Document Sample
_ Plotting Contour Arrhenius Eqn Powered By Docstoc
					% Plotting Contour Arrhenius Eqn
clear
clc
% Model k=ko*exp(-deltaE/(RT))
% k=First Order Rate Constant 1/second
% ko=PreExponential Constant 1/secong
% deltaE=Activation Energy Calories/Mole
% T=Absolute Temperature Degrees Kelvin
% R=1.987 Calories/Mole/Degree K
R=1.987;
ko=1.0e3
deltaEmin=1000;
deltaEmax=5000;
incdeltaE=1000;
Tmin=400;
Tmax=500;
Tinc=20;
disp('deltaE(Cal/Mole) Temp(Deg K) k(1/seconds)')
i=0;
    for deltaE=deltaEmin:incdeltaE:deltaEmax
        i=i+1;
        j=0;
        for T=Tmin:Tinc:Tmax
             j=j+1;
             k(j,i)=ko*exp((-1/R)*deltaE./T);
             fprintf('    %8.2e      %8.2e
%8.2e\n',deltaE,T,k(j,i))
        end
    end
contour([deltaEmin:incdeltaE:deltaEmax],[Tmin:Tinc:Tmax],k)
,xlabel('deltaE Cal/Mole'),ylabel('Temp Deg K'),zlabel('k
1/sec')

				
DOCUMENT INFO