No BSD License  

Highlights from
Wei-Prater Mechanism

image thumbnail

Wei-Prater Mechanism

by

 

24 Mar 2005 (Updated )

Computes reaction paths in a three-species reaction mixture.

main_wei_prater.m
% Wei-Prater Mechanism 
% Author's Data: Housam BINOUS
% Department of Chemical Engineering
% National Institute of Applied Sciences and Technology
% Tunis, TUNISIA
% Email: binoushousam@yahoo.com 

% Main program calls simulink file wei_prater

tf=50;

k12=0.1;
k23=0.5;
k13=0.1;
k31=2;
k32=1;
k21=k12*k31*k23/(k13*k32);

figure(1);
axis off;
hold on;

set(line([0,1],[0,0]),'Color',[1 0 0]);
set(line([0,0.5],[0,0.866025]),'Color',[1 0 0]);
set(line([1,0.5],[0,0.866025]),'Color',[1 0 0]);

set(text(-0.05,-0.05,'A 1'),'Color',[0 0.75 0]);
set(text(0.5,0.916025,'A 3'),'Color',[0 0.75 0]);
set(text(1.05,-0.05,'A 2'),'Color',[0 0.75 0]);

for i=1:11,

a=0;
b=(i-1)*0.1;
c=1-(i-1)*0.1;

sim('wei_prater');

P=[1 0.5;0 0.866025];
 
x=[c2 c3]*P';

plot(x(:,1),x(:,2))

end

for i=1:11,

b=0;
a=(i-1)*0.1;
c=1-(i-1)*0.1;

sim('wei_prater');

P=[1 0.5;0 0.866025];

x=[c2 c3]*P';

plot(x(:,1),x(:,2))

end

for i=1:11,

c=0;
b=(i-1)*0.1;
a=1-(i-1)*0.1;

sim('wei_prater');

P=[1 0.5;0 0.866025];
 
x=[c2 c3]*P';

plot(x(:,1),x(:,2))

end

Contact us