Ignore:
Timestamp:
Jul 25, 2008, 5:03:50 PM (14 years ago)
Author:
gerson bicca
Message:

updated source_testing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/streams.mso

    r575 r576  
    534534Model source_testing
    535535
    536 #Testing a new way to specify stream composition !!!!
    537 # Must be added a set of equations to solve the  "Switch CompostionBasis" !!!!!!
    538 
    539         ATTRIBUTES
     536ATTRIBUTES
    540537        Pallete = true;
    541538        Icon = "icon/Source";
     
    550547         * Temperature
    551548         * Pressure
    552          * Molar (mass or volumetric) composition
     549         * Molar or mass composition
    553550       
    554551        No matter the specification set, the model will calculate some
     
    563560        ";
    564561
    565         PARAMETERS
    566         outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
    567         outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1);
    568                   M(NComp)      as molweight    (Brief = "Component Mol Weight");
    569                   rhoModel              as Switcher             (Brief = "Density model", Valid = ["volume", "correlation"], Default="volume");
     562PARAMETERS
     563        outer PP                                                as Plugin                       (Brief = "External Physical Properties", Type="PP");
     564        outer NComp                             as Integer                      (Brief = "Number of chemical components", Lower = 1);
     565                  M(NComp)                              as molweight    (Brief = "Component Mol Weight");
     566                  rhoModel                              as Switcher             (Brief = "Density model", Valid = ["volume", "correlation"], Default="volume");
    570567                  CompostionBasis               as Switcher             (Brief = "Molar or Mass Compostion", Valid = ["Molar", "Mass"], Default="Molar");
    571568       
    572         Composition(NComp) as fraction  (Brief = "Stream Composition");
    573        
    574         SET
     569
     570SET
    575571
    576572        M   = PP.MolecularWeight();
    577573
    578         VARIABLES
    579         out Outlet                      as stream                       (Brief = "Outlet stream", PosX=1, PosY=0.5256, Symbol="_{out}");
     574VARIABLES
     575
     576        out Outlet                      as stream                       (Brief = "Outlet stream", PosX=1, PosY=0.5256, Symbol="_{out}",Protected=true);
     577       
     578        Composition(NComp) as fraction                  (Brief = "Stream Composition");
     579        F                                                               as flow_mol             (Brief = "Stream Molar Flow Rate");
     580        Fw                                                      as flow_mass            (Brief = "Stream Mass Flow");
     581        Fvol                                    as flow_vol        (Brief = "Volumetric Flow");
     582        T                                                               as temperature  (Brief = "Stream Temperature");
     583        T_Cdeg                                          as temperature  (Brief = "Temperature in °C", Lower=-200);
     584        P                                                               as pressure             (Brief = "Stream Pressure");
     585       
    580586        x(NComp)                        as fraction                     (Brief = "Liquid Molar Fraction",Hidden=true);
    581587        y(NComp)                        as fraction                     (Brief = "Vapour Molar Fraction",Hidden=true);
    582         hl                                      as enth_mol                     (Brief = "Liquid Enthalpy");
    583         hv                                      as enth_mol                     (Brief = "Vapour Enthalpy");
    584         s                                       as entr_mol                     (Brief = "Stream Entropy");
    585         sl                                      as entr_mol                     (Brief = "Liquid Entropy");
    586         sv                                      as entr_mol                     (Brief = "Vapour Entropy");     
    587         zmass(NComp)            as fraction                     (Brief = "Mass Fraction",Protected=true);
    588         Mw                                      as molweight            (Brief = "Average Mol Weight");
    589         vm                                      as volume_mol           (Brief = "Molar Volume");       
    590         rho                                     as dens_mass            (Brief = "Stream Mass Density");
    591         rhom                            as dens_mol                     (Brief = "Stream Molar Density");
    592         Fw                                      as flow_mass            (Brief = "Stream Mass Flow");
    593         Fvol                    as flow_vol         (Brief = "Volumetric Flow");
    594         T_Cdeg                          as temperature          (Brief = "Temperature in °C", Lower=-200);
     588       
     589        Mw                                              as molweight                    (Brief = "Average Mol Weight",Protected=true);
     590        vm                                              as volume_mol           (Brief = "Molar Volume",Protected=true);       
     591        rho                                             as dens_mass                    (Brief = "Stream Mass Density",Protected=true);
     592        rhom                                            as dens_mol                     (Brief = "Stream Molar Density",Protected=true);
     593       
     594        zmass(NComp)            as fraction                             (Brief = "Mass Fraction",Protected=true);
    595595       
    596596        EQUATIONS
     
    617617        Mw = sum(M*Outlet.z);
    618618       
    619         "Flash Calculation"
     619"Flash Calculation"
    620620        [Outlet.v, x, y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
    621        
    622         "Overall Enthalpy"
    623         Outlet.h = (1-Outlet.v)*hl + Outlet.v*hv;
    624 
    625         "Liquid Enthalpy"
    626         hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, x);
    627 
    628         "Vapour Enthalpy"
    629         hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
    630 
    631         "Overall Entropy"
    632         s = (1-Outlet.v)*sl + Outlet.v*sv;
    633 
    634         "Liquid Entropy"
    635         sl = PP.LiquidEntropy(Outlet.T, Outlet.P, x);
    636        
    637         "Vapour Entropy"
    638         sv = PP.VapourEntropy(Outlet.T, Outlet.P, y);
    639 
    640         switch rhoModel
    641                 case "volume":
    642         "Molar Density"
     621
     622"Overall Enthalpy"
     623        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, x) + Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
     624
     625switch rhoModel
     626
     627        case "volume":
     628"Molar Density"
    643629                rhom * vm = 1;
    644630               
    645                 case "correlation":
     631case "correlation":
    646632        "Mass Density"
    647633                rho = (1-Outlet.v)*PP.LiquidDensity(Outlet.T,Outlet.P,x) + Outlet.v*PP.VapourDensity(Outlet.T,Outlet.P,y);
    648         end
    649        
    650         "Mass or Molar Density"
     634end
     635       
     636"Mass or Molar Density"
    651637        rhom * Mw = rho;
    652638
    653         "Flow Mass"
     639"Flow Mass"
    654640        Fw      =  Mw*Outlet.F;
    655641
    656         "Molar Volume"
     642"Molar Volume"
    657643        vm = (1-Outlet.v)*PP.LiquidVolume(Outlet.T, Outlet.P, x) + Outlet.v*PP.VapourVolume(Outlet.T,Outlet.P,y);
    658644       
    659         "Volumetric Flow"
     645"Volumetric Flow"
    660646        Fvol = Outlet.F*vm ;
    661647       
    662         "Temperature in °C"
     648"Temperature in °C"
    663649        T_Cdeg = Outlet.T - 273.15 * 'K';
    664650
    665 end
     651"Equate Flow"
     652        Outlet.F = F;
     653
     654"Equate Pressures"
     655        Outlet.P = P;
     656
     657"Equate Temperatures"
     658        Outlet.T = T;
     659
     660end
Note: See TracChangeset for help on using the changeset viewer.