simulated moving bed process.gpj

gprom语言

<?xml version="1.0" encoding="ISO-8859-1"?>

<gMB:GpromsProject xmlns:gMB="http://www.psenterprise.com/gMB" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.psenterprise.com/gMB http://www.psenterprise.com/gPJ.xsd" version="2" gMB_Version="2.3.4" gMB_Build="2773">
    <Properties>
        <StringProperty name="CreatedBy">aa</StringProperty>
        <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty>
        <StringProperty name="ModifiedBy">Administrator</StringProperty>
        <DateProperty name="ModifiedDate">2007-07-05T01:30:58Z</DateProperty>
        <BooleanProperty name="DisplayEmptyGroups">false</BooleanProperty>
        <BooleanProperty name="LibraryProject">false</BooleanProperty>
    </Properties>
    <ImageCache>
    </ImageCache>
    <ModelReferenceCache>
    </ModelReferenceCache>
    <Group name="Variable Types">
        <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-28T03:48:28Z</DateProperty>
        </Properties>
        <VariableTypeEntity name="Axial_dispersion">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:12:12Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-28T03:48:28Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
            </Properties>
            <DefaultValue>1.0</DefaultValue>
            <MinValue>-1.0E30</MinValue>
            <MaxValue>1.0E30</MaxValue>
            <Units></Units>
        </VariableTypeEntity>
        <VariableTypeEntity name="Concentration">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:09:57Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-22T07:09:57Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
            </Properties>
            <DefaultValue>1.0</DefaultValue>
            <MinValue>-1.0E30</MinValue>
            <MaxValue>1.0E30</MaxValue>
            <Units></Units>
        </VariableTypeEntity>
        <VariableTypeEntity name="Counter">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:11:17Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-22T07:11:17Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
            </Properties>
            <DefaultValue>1.0</DefaultValue>
            <MinValue>-1.0E30</MinValue>
            <MaxValue>1.0E30</MaxValue>
            <Units></Units>
        </VariableTypeEntity>
        <VariableTypeEntity name="Main_l_c1_c2">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:24:58Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-22T07:24:58Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
            </Properties>
            <DefaultValue>1.0</DefaultValue>
            <MinValue>-1.0E30</MinValue>
            <MaxValue>1.0E30</MaxValue>
            <Units></Units>
        </VariableTypeEntity>
        <VariableTypeEntity name="Total_volume_flow">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:12:29Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-22T07:12:29Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
            </Properties>
            <DefaultValue>1.0</DefaultValue>
            <MinValue>-1.0E30</MinValue>
            <MaxValue>1.0E30</MaxValue>
            <Units></Units>
        </VariableTypeEntity>
        <VariableTypeEntity name="Velocity">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:11:54Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-22T07:11:54Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
            </Properties>
            <DefaultValue>1.0</DefaultValue>
            <MinValue>-1.0E30</MinValue>
            <MaxValue>1.0E30</MaxValue>
            <Units></Units>
        </VariableTypeEntity>
    </Group>
    <Group name="Stream Types">
        <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-28T02:39:18Z</DateProperty>
        </Properties>
        <StreamTypeEntity name="Connect_VCC" type="psenterprise.gpromskit.models.StreamTypeEntity" version="1">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:25:22Z</DateProperty>
                <StringProperty name="ModifiedBy">Administrator</StringProperty>
                <DateProperty name="ModifiedDate">2007-07-05T01:30:58Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">true</BooleanProperty>
            </Properties>
            <Body><![CDATA[IS Total_volume_flow, Concentration, Concentration]]></Body>
        </StreamTypeEntity>
    </Group>
    <Group name="Connection Types">
        <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="Models">
        <Properties>
            <StringProperty name="CreatedBy">aa</StringProperty>
            <DateProperty name="CreatedDate">2006-11-22T07:05:55Z</DateProperty>
            <StringProperty name="ModifiedBy">LJQ</StringProperty>
            <DateProperty name="ModifiedDate">2006-11-27T15:06:52Z</DateProperty>
        </Properties>
        <ModelEntity name="Column" type="psenterprise.gpromskit.models.ModelEntity" version="4">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:16:09Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-28T03:49:44Z</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 length 
    ColumnLength  AS     REAL  DEFAULT  1.53    # 3,53  # m                         
    # column cross section area 
    Fl            AS    REAL  DEFAULT  0.002123        # m^2                  
    # porosity (here external porosity)
    ExPorosity    AS    REAL  DEFAULT  0.34
    # axial dispersion coefficient
    D_ax_koeff    AS     REAL  DEFAULT  0                        
    # molecular diffusion as an additional axial dispersion coefficient
    D_add         AS     REAL  DEFAULT  0.0
    # Henry coefficient
    #K_1          AS     REAL  DEFAULT  12.32
    # Langmuir coefficient for component 1 
    a1            AS     REAL  DEFAULT  72.97297    #72.97297
    # langmuir coefficient for component 1
    b1            AS     REAL  DEFAULT  0.4840135  #0.4840135
    # Langmuir coefficient for component 1 
    a2            AS     REAL  DEFAULT  0.32731563
    # Henry coefficient for component 2
    K_2           AS     REAL  DEFAULT  1.3066   
    # mass transfer coefficient for component 1
    ks1           AS     REAL  DEFAULT 0.34191*60      # h^-1
    # mass transfer coefficient for component 2
    ks2           AS     REAL  DEFAULT 0.34555*60      # h^-1
DISTRIBUTION_DOMAIN
    Axial         AS  [0 : ColumnLength]  
VARIABLE
    # intersticial velocity in column   
    v_int          AS     Velocity
    # total axial dispersion coefficient
    D_ax          AS      Axial_dispersion
    # inlet concentration of desorbent
    C1_in         AS      Concentration
    C2_in         AS      Concentration
    # outlet concentration of desorbent
    C1_out        AS      Concentration
    C2_out        AS      Concentration
    # flow rate in the column
    L_in          AS      Total_volume_flow
    C1            AS      DISTRIBUTION(Axial)   OF   Concentration
    Q1_equ        AS      DISTRIBUTION(Axial)   OF   Concentration
    Q1            AS      DISTRIBUTION(Axial)   OF   Concentration
    C2            AS      DISTRIBUTION(Axial)   OF   Concentration
    Q2_equ        AS      DISTRIBUTION(Axial)   OF   Concentration
    Q2            AS      DISTRIBUTION(Axial)   OF   Concentration
    
STREAM 
    Inlet:   L_in, C1_in, C2_in     AS Connect_VCC
    Outlet : L_in, C1_out, C2_out   AS Connect_VCC
SET 
     Axial := [BFDM,2,40];
BOUNDARY
     C1(0) = C1_in ;
     C2(0) = C2_in ;
     C1(ColumnLength) = C1_out;
     C2(ColumnLength) = C2_out;
EQUATION
    v_int=L_in/(ExPorosity*Fl);
    D_ax=v_int*D_ax_koeff+D_add;
    # mass balance using Transport-Dispersion (TD) model in the flow phase
    FOR z := 0|+ TO ColumnLength DO
        $C1(z)=(D_ax_koeff*PARTIAL(C1(z),Axial,Axial)-(v_int*PARTIAL(C1(z),Axial))-((1-ExPorosity)/ExPorosity*ks1*(Q1_equ(z)-Q1(z)))); # the first component
        $C2(z)=(D_ax_koeff*PARTIAL(C2(z),Axial,Axial)-(v_int*PARTIAL(C2(z),Axial))-((1-ExPorosity)/ExPorosity*ks2*(Q2_equ(z)-Q2(z)))); # the second component
    END
    # linear driving force model to descrip the adsorbent phase mass balance for tow components    
    FOR z := 0 TO ColumnLength DO
        $Q1(z)=ks1*(Q1_equ(z)-Q1(z));
        $Q2(z)=ks2*(Q2_equ(z)-Q2(z));
        Q1_equ(z)=a1*C1(z)/(1+b1*C1(z))+a2*C1(z); # the modified langmuir isotherms
        Q2_equ(z)=K_2*C2(z);  # the linear isotherms
    END
]]></Body>
        </ModelEntity>
        <ModelEntity name="Connect" type="psenterprise.gpromskit.models.ModelEntity" version="4">
            <Properties>
                <StringProperty name="CreatedBy">aa</StringProperty>
                <DateProperty name="CreatedDate">2006-11-22T07:08:59Z</DateProperty>
                <StringProperty name="ModifiedBy">aa</StringProperty>
                <DateProperty name="ModifiedDate">2006-11-28T03:52:08Z</DateProperty>
                <BooleanProperty name="Win32LineSeparation">false</BooleanProperty>
                <StringProperty name="DefaultTab">gPROMS language</StringProperty>
                <BooleanProperty name="DisplayedInPalette">true</BooleanProperty>
            </Properties>
            <Icon>-199773695</Icon>
            <Body><![CDATA[PARAMETER
  
    # feed concentration for component 1
    C1_f       AS REAL    DEFAULT  585.8         # g/L        
    # feed concentration for component 2                                            
    C2_f       AS REAL    DEFAULT  50            # g/L  
    # volumic velocity for feed                                            
    Vo_feed    AS REAL    DEFAULT 900e-6        # m^3/h                         
    # volumic velocity for extract