Changeset 297 for trunk


Ignore:
Timestamp:
Jun 28, 2007, 9:27:02 AM (16 years ago)
Author:
Argimiro Resende Secchi
Message:

Added switcher to select a model to evaluate the mixture density.

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/eml/streams.mso

    r290 r297  
    126126        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
    127127        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1);
    128                         M(NComp)  as molweight  (Brief="Component Mol Weight");
     128                  M(NComp)      as molweight    (Brief = "Component Mol Weight");
     129                  rhoModel              as Switcher             (Brief = "Density model", Valid = ["volume", "correlation"], Default="volume");
    129130       
    130131        SET
    … …  
    133134
    134135        VARIABLES
    135         out Outlet as stream;
    136         x(NComp) as fraction(Brief = "Liquid Molar Fraction");
    137         y(NComp) as fraction(Brief = "Vapour Molar Fraction");
    138         hl as enth_mol;
    139         hv as enth_mol;
    140         zmass(NComp)    as fraction                     (Brief = "Mass Fraction");
    141         Mw                                      as molweight            (Brief="Average Mol Weight");
    142         vm                                      as volume_mol   (Brief="Molar Volume");
    143         rho                                     as dens_mass            (Brief="Stream Density");
    144         Fw                                      as flow_mass            (Brief="Stream Mass Flow");
    145         Fvol                            as flow_vol         (Brief = "Volumetric Flow");
     136        out Outlet                      as stream;
     137        x(NComp)                        as fraction                     (Brief = "Liquid Molar Fraction");
     138        y(NComp)                        as fraction                     (Brief = "Vapour Molar Fraction");
     139        hl                                      as enth_mol;
     140        hv                                      as enth_mol;
     141        zmass(NComp)            as fraction                     (Brief = "Mass Fraction");
     142        Mw                                      as molweight            (Brief = "Average Mol Weight");
     143        vm                                      as volume_mol           (Brief = "Molar Volume");       
     144        rho                                     as dens_mass            (Brief = "Stream Mass Density");
     145        rhom                            as dens_mol                     (Brief = "Stream Molar Density");
     146        Fw                                      as flow_mass            (Brief = "Stream Mass Flow");
     147        Fvol                    as flow_vol         (Brief = "Volumetric Flow");
    146148       
    147149        EQUATIONS
    … …  
    159161        Mw = sum(M*Outlet.z);
    160162
     163        switch rhoModel
     164                case "volume":
     165        "Molar Density"
     166                rhom * vm = 1;
     167               
     168                case "correlation":
    161169        "Mass Density"
    162         rho =   (1-Outlet.v)*PP.LiquidDensity(Outlet.T,Outlet.P,x) + Outlet.v*PP.VapourDensity(Outlet.T,Outlet.P,y);
     170                rho*((1-Outlet.v)/PP.LiquidDensity(Outlet.T,Outlet.P,x) + Outlet.v/PP.VapourDensity(Outlet.T,Outlet.P,y)) = 1;
     171        end
     172       
     173        "Mass or Molar Density"
     174        rhom * Mw = rho;
    163175
    164176        "Flow Mass"
    … …  
    171183        Fvol = Outlet.F*vm ;
    172184       
    173 "Mass Fraction"
     185        "Mass Fraction"
    174186        zmass = M*Outlet.z / Mw;       
    175187       
    … …  
    197209        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
    198210        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1);
    199                         M(NComp)  as molweight  (Brief="Component Mol Weight");
     211                  M(NComp)      as molweight    (Brief = "Component Mol Weight");
     212                  rhoModel              as Switcher             (Brief = "Density model", Valid = ["volume", "correlation"], Default="volume");
    200213       
    201214        SET
    … …  
    204217       
    205218        VARIABLES
    206         in Inlet as stream;
    207         v as fraction;
    208         x(NComp)        as fraction                     (Brief = "Liquid Molar Fraction");
    209         y(NComp)        as fraction                     (Brief = "Vapour Molar Fraction");
    210         zmass(NComp)    as fraction     (Brief = "Mass Fraction");
    211         Mw                      as molweight            (Brief="Average Mol Weight");
    212         vm                      as volume_mol   (Brief="Molar Volume");
    213         rho                     as dens_mass            (Brief="Stream Density");
    214         Fw                      as flow_mass            (Brief="Stream Mass Flow");
    215         Fvol          as flow_vol          (Brief = "Volumetric Flow");
     219        in Inlet                as stream;
     220        v                               as fraction;
     221        x(NComp)                as fraction             (Brief = "Liquid Molar Fraction");
     222        y(NComp)                as fraction             (Brief = "Vapour Molar Fraction");
     223        zmass(NComp)    as fraction             (Brief = "Mass Fraction");
     224        Mw                              as molweight    (Brief = "Average Mol Weight");
     225        vm                              as volume_mol   (Brief = "Molar Volume");       
     226        rho                             as dens_mass    (Brief = "Stream Mass Density");
     227        rhom                    as dens_mol             (Brief = "Stream Molar Density");
     228        Fw                              as flow_mass    (Brief = "Stream Mass Flow");
     229        Fvol            as flow_vol     (Brief = "Volumetric Flow");
    216230       
    217231        EQUATIONS
    … …  
    222236        Mw = sum(M*Inlet.z);
    223237
     238        switch rhoModel
     239                case "volume":
     240        "Molar Density"
     241                rhom * vm = 1;
     242               
     243                case "correlation":
    224244        "Mass Density"
    225         rho =   (1-v)*PP.LiquidDensity(Inlet.T,Inlet.P,x) + v*PP.VapourDensity(Inlet.T,Inlet.P,y);
     245                rho * ((1-v)/PP.LiquidDensity(Inlet.T,Inlet.P,x) + v/PP.VapourDensity(Inlet.T,Inlet.P,y)) = 1;
     246        end
     247       
     248        "Mass or Molar Density"
     249        rhom * Mw = rho;
    226250
    227251        "Flow Mass"

  • trunk/sample/stage_separators/sample_batch_dist_PID.mso

    r247 r297  
    4747        reflux  as stream;
    4848        feed    as stream;
    49         pidP    as PID_Ideal_AWBT;
     49        pidP    as PID;
    5050       
    5151        CONNECTIONS
    … …  
    9696
    9797        SET
     98        pidP.PID_Select = "Ideal_AWBT";
    9899        batch.V = 3 * 'm^3';
    99100        batch.Across = 1 * 'm^2';

  • trunk/sample/stage_separators/sample_column.mso

    r270 r297  
    422422        feed as source;
    423423        zero as stream;
    424         TCcond as PIDIncr_Ideal_AW;
    425         LCtop as PIDIncr_Ideal_AW;
    426         LCbot as PIDIncr_Ideal_AW;
    427         PC as PIDIncr_Ideal_AW;
     424        TCcond as PIDIncr;
     425        LCtop as PIDIncr;
     426        LCbot as PIDIncr;
     427        PC as PIDIncr;
    428428
    429429        CONNECTIONS
    … …  
    548548        col.reb.Across = 1 * 'm^2';
    549549
     550        # Controllers type
     551        TCcond.PID_Select = "Ideal_AW";
     552        PC.PID_Select = "Ideal_AW";
     553        LCtop.PID_Select = "Ideal_AW";
     554        LCbot.PID_Select = "Ideal_AW";
     555
    550556        Qrmax = 5e6 * 'kJ/h';
    551557        Qrmin = 1e6 * 'kJ/h';
Note: See TracChangeset for help on using the changeset viewer.