smb_mol.asv

模拟移动床色谱的计算机模拟

clc
clear;
tic;
colordef white;
colordef white;
global alpha Pe poro K1 K2 v L tc initial
K1=0.56;
K2=0.23;
coi=[50 50];
tc=10.3;
L=47.5; % L cm
R=0.46; %R radius of column  cm
poro=0.85;
Pe=[500;500];
v=[14.0442 10.5294 13.0332 10.0314];
cycle=300;
for i=1:cycle
    if i==1
        initial=zeros(8,1);
    else
       temp1=cloumn(1);
       temp2=cloumn(2);
       cloumn(1)=cloumn(3);
       cloumn(3)=cloumn(5);
       cloumn(5)=cloumn(7);
       cloumn(7)=temp1;
       cloumn(2)=cloumn(3);
       cloumn(4)=cloumn(6);
       cloumn(6)=cloumn(8);
       cloumn(8)=temp2;
       initial=cloumn;
    end
    sol=edmodel;
    cloumn=[sol(end,end,1);sol(end,end,2);sol(end,end,3);sol(end,end,4);sol(end,end,5);sol(end,end,6);sol(end,end,7);sol(end,end,8)]
    record{i}=sol;
end

history=[];
for i=1:cycle
history=[history;record{i}(:,end,1)];
end
plot(history);


aa=record{cycle};
plot([aa(end,:,1),aa(end,:,3),aa(end,:,5),aa(end,:,7)],'mo')    


toc