smb_edmodel.m

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

function []=smb_EDmodel
%   the programe for simulated moving bed with MOL method 
%   Author: Jianqi Liu
%   Southern yangtze university, Wuxi, PRC
%   $Date: 2006/05/01 $
%   illumination: this programe 
clear;
Q1_s=0.0079146;  % the flow rates through the Zone 1 [m/s];
Q2_s=0.0054106;  % the flow rates through the Zone 2 [m/s];
Q3_s=0.0059175;  % the flow rates through the Zone 3 [m/s];
Q4_s=0.0049273;  % the flow rates through the Zone 4 [m/s];

Q1=Q1_s*60;  % Q1 the flow rates through the Zone 1 [m/min];
Q2=Q2_s*60;  % Q2 the flow rates through the Zone 2 [m/min];
Q3=Q3_s*60;  % Q3 the flow rates through the Zone 3 [m/min];
Q4=Q4_s*60;  % Q4 the flow rates through the Zone 4 [m/min];

Qd=Q1-Q4; %  Qd the desorbent glow rate;
Qe=Q1-Q2; %  Qe the extract glow rate;
Qf=Q3-Q2; %  Qf the feed glow rate;
Qr=Q3-Q4; %  Qr the raffinate glow rate;

Cf=2.9; %  Cf the feed concentration [kg/m3]; 
Ce=0; %  Cf the extract concentration [kg/m3]; 
Cd=0; %  Cf the desorbent concentration [kg/m3]; 
Cr=0; %  Cf the raffinate concentration [kg/m3]; 

C1_in=0;  % C1_in the concentration rates through the Zone 1 at the inlet [kg/m3];
C1_ou=0   % C1_ou the concentration rates through the Zone 1 at the outlet [kg/m3]; 
C2_in=0;  % C2_in the concentration rates through the Zone 2 at the inlet  [kg/m3];
C2_ou=0;  % C2_ou the concentration rates through the Zone 2 at the outlet [kg/m3];
C3_in=0;  % C3_in the concentration rates through the Zone 3 at the inlet [kg/m3];
C3_ou=0;  % C3_ou the concentration rates through the Zone 3 at the outlet [kg/m3];
C4_in=0;  % C4_in the concentration rates through the Zone 4 at the inlet [kg/m3];
C4_ou=0;  % C4_ou the concentration rates through the Zone 4 at the outlet [kg/m3];

%  the node model 
C3_in=(C2_ou*Q2+Cf*Qf)/Q3;

C4_in=C3_ou;

C1_in=(C2_ou*Q4+Cd*Qd)/Q1;

C2_in=C1_ou;