📄 simulated moving bed process.gpj
字号:
Vo_ex AS REAL DEFAULT 1800e-6 # m^3/h # volumic velocity for raffinate Vo_raff AS REAL DEFAULT 2700e-6 # m^3/h # volumic velocity for eluent Vo_el AS REAL DEFAULT 3600e-6 #2700e-6 # m^3/h NC AS INTEGER DEFAULT 10 VARIABLE C1_in AS Concentration C2_in AS Concentration C1_out AS Concentration C2_out AS Concentration C1_raff AS Concentration C2_raff AS Concentration C1_ex AS Concentration C2_ex AS Concentration F_str AS Counter El_str AS Counter EX_str AS Counter R_str AS Counter zwstrom AS ARRAY(NC) OF Counter Vo_in AS Total_volume_flow Vo_out AS Total_volume_flow STREAM Inlet : Vo_in, C1_in, C2_in AS Connect_VCC Outlet : Vo_out, C1_out, C2_out AS Connect_VCC EQUATION Vo_out = SIGMA(zwstrom)*Vo_in+El_str*Vo_el+F_str*(Vo_feed+Vo_in)+EX_str*(-Vo_ex+Vo_in)+R_str*(-Vo_raff+Vo_in); C1_out = SIGMA(zwstrom)*C1_in+F_str*((Vo_in*C1_in+Vo_feed*C1_f)/(Vo_in+Vo_feed))+El_str*(Vo_in*C1_in/(Vo_in+Vo_el))+R_str*C1_in+Ex_str*C1_in; C2_out = SIGMA(zwstrom)*C2_in+F_str*((Vo_in*C2_in+Vo_feed*C2_f)/(Vo_in+Vo_feed))+El_str*(Vo_in*C2_in/(Vo_in+Vo_el))+R_str*C2_in+Ex_str*C2_in; C1_ex = Ex_str*C1_out; C2_ex = Ex_str*C2_out; C1_raff= R_str*C1_out; C2_raff= R_str*C2_out; ]]></Body> </ModelEntity> <ModelEntity name="ExterColumn" type="psenterprise.gpromskit.models.ModelEntity" version="4"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:17:48Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-28T03:50:06Z</DateProperty> <BooleanProperty name="Win32LineSeparation">false</BooleanProperty> <StringProperty name="DefaultTab">gPROMS language</StringProperty> <BooleanProperty name="DisplayedInPalette">true</BooleanProperty> </Properties> <Icon>-199773695</Icon> <Body><![CDATA[PARAMETER # column cross section area ExColumnArea AS REAL DEFAULT 0.00002826 # m^2 # column length ExLength AS REAL DEFAULT 0.12 # m # axail dispersion coefficient D_ax_koeff AS REAL DEFAULT 0DISTRIBUTION_DOMAIN Axial_2 AS [ 0 : ExLength ] VARIABLE v_vel AS Velocity L_in AS Total_volume_flow C1_in AS Concentration C2_in AS Concentration C1_out AS Concentration C2_out AS Concentration D_ax AS Axial_dispersion C1 AS DISTRIBUTION(Axial_2) OF Concentration C2 AS DISTRIBUTION(Axial_2) OF ConcentrationSTREAM Inlet : L_in, C1_in, C2_in AS Connect_VCC Outlet : L_in, C1_out, C2_out AS Connect_VCCSET Axial_2 := [BFDM,2,10]; BOUNDARY C1(0) = C1_in ; C2(0) = C2_in ; C1(ExLength) = C1_out; C2(ExLength) = C2_out;EQUATION v_vel=L_in/(ExColumnArea); D_ax=v_vel*D_ax_koeff; FOR z := 0|+ TO ExLength DO $C1(z)=D_ax_koeff*PARTIAL(C1(z),Axial_2,Axial_2)-(v_vel*PARTIAL(C1(z),Axial_2)); $C2(z)=D_ax_koeff*PARTIAL(C2(z),Axial_2,Axial_2)-(v_vel*PARTIAL(C2(z),Axial_2)); END]]></Body> </ModelEntity> </Group> <Group name="Tasks"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> <Group name="Processes"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:19:10Z</DateProperty> </Properties> <ProcessEntity name="SMBSimulation" type="psenterprise.gpromskit.models.ProcessEntity" version="1"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:19:10Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-28T03:45:06Z</DateProperty> <BooleanProperty name="Win32LineSeparation">true</BooleanProperty> </Properties> <Body><![CDATA[PARAMETER NoColumn AS INTEGER DEFAULT 10 # the number of column in your SMB system INTERVALL AS REAL DEFAULT 21/60 # the switching time # hUNIT V AS ARRAY(NoColumn) of ExterColumn # the dead volumne model M AS ARRAY(NoColumn) of Connect # connect columns with flowrate and concentration T AS ARRAY(NoColumn) of Column # the main models:the single column model EQUATION FOR i := 1 to NoColumn DO IF INT(time/INTERVALL)-10*INT(INT(time/INTERVALL)/10)+1=i then M(i).El_str=1; M(i+2-10*INT((i+1)/10)).Ex_str=1; M(i+5-10*INT((i+4)/10)).F_str=1; M(i+7-10*INT((i+6)/10)).R_str=1; M(i).zwstrom(i)=0; M(i+2-10*INT((i+1)/10)).zwstrom(i)=0; M(i+5-10*INT((i+4)/10)).zwstrom(i)=0; M(i+7-10*INT((i+6)/10)).zwstrom(i)=0; M(i+1-10*INT((i)/10)).zwstrom(i)=1; M(i+3-10*INT((i+2)/10)).zwstrom(i)=1; M(i+4-10*INT((i+3)/10)).zwstrom(i)=1; M(i+6-10*INT((i+5)/10)).zwstrom(i)=1; M(i+8-10*INT((i+7)/10)).zwstrom(i)=1; M(i+9-10*INT((i+8)/10)).zwstrom(i)=1; ELSE M(i).El_str=0; M(i+2-10*INT((i+1)/10)).Ex_str=0; M(i+5-10*INT((i+4)/10)).F_str=0; M(i+7-10*INT((i+6)/10)).R_str=0; M(i).zwstrom(i)=0; M(i+2-10*INT((i+1)/10)).zwstrom(i)=0; M(i+5-10*INT((i+4)/10)).zwstrom(i)=0; M(i+7-10*INT((i+6)/10)).zwstrom(i)=0; M(i+1-10*INT((i)/10)).zwstrom(i)=0; M(i+3-10*INT((i+2)/10)).zwstrom(i)=0; M(i+4-10*INT((i+3)/10)).zwstrom(i)=0; M(i+6-10*INT((i+5)/10)).zwstrom(i)=0; M(i+8-10*INT((i+7)/10)).zwstrom(i)=0; M(i+9-10*INT((i+8)/10)).zwstrom(i)=0; END # END IF END # END FOR FOR j := 1 to NoColumn DO M(j).outlet = T(j).inlet; # Co. Connect and Column T(j).outlet = V(j+1-10*INT((j)/10)).inlet; # Co. Column and ExterColumn V(j+1-10*INT((j)/10)).outlet =M(j+1-10*INT((j)/10)).inlet; # Co. ExterColumn and Connect END # END FORINITIAL# set the initial values FOR k := 1 to NoColumn DO FOR z := 0|+ to T(1).ColumnLength DO WITHIN T(k) DO C1(z) = 0.0; C2(z) = 0.0; END # END WITHIN END FOR z := 0 to T(1).ColumnLength DO WITHIN T(k) DO Q1(z) = 0.0; Q2(z) = 0.0; END # END WITHIN END FOR z := 0|+ to V(1).ExLength DO WITHIN V(k) DO C1(z) = 0.0; C2(z) = 0.0; END # END WITHIN END END # END FOR SOLUTIONPARAMETERS ReportingInterval :=0.035 # hSCHEDULE SEQUENCE Continue for 105 # h END # END SEQUENCE ]]></Body> </ProcessEntity> </Group> <Group name="Optimisations"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> <Group name="Estimations"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> <Group name="Experiment Designs"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> <Group name="Experiments"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> <Group name="... for parameter estimation"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> <Group name="... for experiment design"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> </Group> <Group name="Saved Variable Sets"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group> <Group name="Miscellaneous Files"> <Properties> <StringProperty name="CreatedBy">aa</StringProperty> <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty> <StringProperty name="ModifiedBy">aa</StringProperty> <DateProperty name="ModifiedDate">2006-11-22T07:05:55Z</DateProperty> </Properties> </Group></gMB:GpromsProject>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -