smb_mol.m

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

clc
clear;
tic;


colordef white;
colordef white;
global Pe poro K1 K2 v L tc int

K1=0.56;
K2=0.26;
coi=[50 50];
tc=10.3;
L=47.5; % L cm
R=0.46; %R radius of column  cm
poro=0.51;
Pe=[500;500];
v=[14.0442 10.5294 13.0332 10.0314];
cycle=20;
for i=1:cycle
    if i==1
        int=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(4);
        cloumn(4)=cloumn(6);
        cloumn(6)=cloumn(8);
        cloumn(8)=temp2;
        int=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;

    i
    
    
end
figure;
history1=[];history2=[];
for i=1:cycle
    history1=[history1;record{i}(:,end,1)];
    history2=[history2;record{i}(:,end,2)];
end
plot(history1,'r');hold on;
plot(history2,'k');

figure;
aa=record{cycle};
plot([aa(end,:,1),aa(end,:,3),aa(end,:,5),aa(end,:,7)],'mo');hold on;
plot([aa(end,:,2),aa(end,:,4),aa(end,:,6),aa(end,:,8)],'b+')    

toc