Changeset 664


Ignore:
Timestamp:
Oct 10, 2008, 1:45:07 AM (14 years ago)
Author:
gerson bicca
Message:

added info_stream model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/streams.mso

    r651 r664  
    588588
    589589end
     590
     591Model info_stream
     592       
     593ATTRIBUTES
     594        Pallete = true;
     595        Icon = "icon/Info_Stream";
     596        Brief = "Material stream information";
     597        Info = "
     598        This model should be used for middle streams when additional
     599        information about the stream is desired.
     600
     601        Some of the additional informations calculated by this models are:
     602         * Mass density
     603         * Mass flow
     604         * Mass compostions
     605         * Specific volume
     606         * Vapour fraction
     607         * Volumetric flow
     608         * Liquid and Vapour compositions
     609         * Viscosity
     610         * Heat Capacity
     611         * Thermal Conductivity
     612         * Temperature in Celsius Degrees
     613        ";
     614
     615PARAMETERS
     616        outer PP                        as Plugin                       (Brief = "External Physical Properties", Type="PP");
     617        outer NComp     as Integer                      (Brief = "Number of chemical components", Lower = 1);
     618                  M(NComp)      as molweight    (Brief = "Component Mol Weight");
     619       
     620SET
     621
     622        M   = PP.MolecularWeight();
     623       
     624VARIABLES
     625
     626        in      Inlet           as stream               (Brief = "Inlet Stream", PosX=0, PosY=0.5308, Protected=true , Symbol="_{in}");
     627        out     Outlet          as stream               (Brief = "Outlet Stream", PosX=1, PosY=0.5308, Protected=true , Symbol="_{out}");
     628       
     629        v                                                       as fraction                     (Brief = "Vapourization fraction",Hidden=true);
     630        x(NComp)                        as fraction                     (Brief = "Liquid Molar Fraction",Hidden=true);
     631        y(NComp)                        as fraction                     (Brief = "Vapour Molar Fraction",Hidden=true);
     632       
     633        F(NComp)        as flow_mol             (Brief = "Component Molar Flow",Protected=true);
     634        FwTotal         as flow_mass            (Brief = "Total Mass Flow",Protected=true);
     635        Fw(NComp)               as flow_mass            (Brief = "Component Mass Flow",Protected=true);
     636        FvolTotal           as flow_vol         (Brief = "Total Volumetric Flow",Protected=true);
     637        T_Cdeg                          as temperature          (Brief = "Temperature in °C", Lower=-200,Protected=true);
     638
     639        Mu      as viscosity            (Brief="Stream Viscosity",Lower=0.0001, Symbol = "\mu",Protected=true);
     640        Cp              as cp_mol                       (Brief="Stream Molar Heat Capacity", Upper=1e10,Protected=true);       
     641        K               as conductivity         (Brief="Stream Thermal Conductivity", Default=1.0, Lower=1e-5, Upper=500,Protected=true);
     642        Mw                                              as molweight            (Brief = "Average Mol Weight",Protected=true);
     643        vm              as volume_mol   (Brief = "Molar Volume",Protected=true);       
     644        rho             as dens_mass            (Brief = "Stream Mass Density",Protected=true);
     645        rhom    as dens_mol             (Brief = "Stream Molar Density",Protected=true);
     646        s               as entr_mol             (Brief = "Stream Entropy",Protected=true);
     647        zmass(NComp) as fraction                        (Brief = "Mass Fraction",Protected=true);
     648       
     649EQUATIONS
     650
     651"Flash Calculation"
     652        [v, x, y] = PP.FlashPH(Inlet.P, Inlet.h, Inlet.z);
     653       
     654"Average Molecular Weight"
     655        Mw = sum(M*Inlet.z);
     656
     657"Mass Density"
     658        rho * ((1-v)/PP.LiquidDensity(Inlet.T,Inlet.P,x) + v/PP.VapourDensity(Inlet.T,Inlet.P,y)) = 1;
     659       
     660"Mass or Molar Density"
     661        rhom * Mw = rho;
     662
     663"Total Flow Mass"
     664        FwTotal =  Mw*Inlet.F;
     665
     666"Component Flow Mass"
     667        Fw      =  FwTotal*zmass;
     668
     669"Molar Volume"
     670        vm = (1-v)*PP.LiquidVolume(Inlet.T, Inlet.P, x) + v*PP.VapourVolume(Inlet.T,Inlet.P,y);
     671       
     672"Total Volumetric Flow"
     673        FvolTotal = Inlet.F*vm ;
     674       
     675"Mass Fraction"
     676        zmass = M*Inlet.z / Mw;
     677
     678"Stream Heat Capacity"
     679        Cp      =       (1-v)*PP.LiquidCp(Inlet.T, Inlet.P, x) + v*PP.VapourCp(Inlet.T,Inlet.P,y);
     680
     681"Stream Viscosity"
     682        Mu      =       (1-v)*PP.LiquidViscosity(Inlet.T, Inlet.P, x) + v*PP.VapourViscosity(Inlet.T,Inlet.P,y);
     683
     684"Stream ThermalConductivity"
     685        K       =       (1-v)*PP.LiquidThermalConductivity(Inlet.T, Inlet.P, x) + v*PP.VapourThermalConductivity(Inlet.T,Inlet.P,y);
     686
     687"Stream Overall Entropy"
     688        s = (1-v)*PP.LiquidEntropy(Inlet.T, Inlet.P, x) + v*PP.VapourEntropy(Inlet.T, Inlet.P, y);
     689       
     690"Temperature in °C"
     691        T_Cdeg = Inlet.T - 273.15 * 'K';
     692
     693"Outlet Flow"
     694        Outlet.F = Inlet.F;
     695
     696"Component Molar Flow"
     697        F = Inlet.F*Inlet.z;
     698
     699"Outlet Temperature"
     700        Outlet.T = Inlet.T;
     701
     702"Outlet Pressure"
     703        Outlet.P = Inlet.P;
     704
     705"Outlet Vapour Fraction"
     706        Outlet.v = Inlet.v;
     707
     708"Outlet Enthalpy"
     709        Outlet.h = Inlet.h;
     710
     711"Outlet Composition"
     712        Outlet.z= Inlet.z;
     713end
Note: See TracChangeset for help on using the changeset viewer.