source: branches/new_gui/eml/streams/sinks.mso @ 889

#*-------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
*
* This LIBRARY is free software; you can distribute it and/or modify
* it under the therms of the ALSOC FREE LICENSE as available at
* http://www.enq.ufrgs.br/alsoc.
*
* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
* from http://www.rps.eng.br Copyright (C) 2002-2004.
* All rights reserved.
*
* EMSO is distributed under the therms of the ALSOC LICENSE as
* available at http://www.enq.ufrgs.br/alsoc.
*
*--------------------------------------------------------------------
* Model of sink streams
*----------------------------------------------------------------------
* Author: Paula B. Staudt and Rafael de P. Soares
* $Id: sink.mso 757 2009-06-03 20:07:22Z bicca $
*---------------------------------------------------------------------*#

using "streams";

Model sink
        ATTRIBUTES
        Pallete = true;
        Icon = "icon/Sink";
        Brief = "Material stream sink";
        Info = "
        This model should be used for boundary streams when additional
        information about the stream is desired.

        Some of the additional informations calculated by this models are:
         * Mass density
         * Mass flow
         * Mass compostions
         * Specific volume
         * Vapour fraction
         * Volumetric flow
         * Liquid and Vapour compositions
        ";

        PARAMETERS
        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1); 
                  M(NComp)      as molweight    (Brief = "Component Mol Weight");
        
        SET

        M   = PP.MolecularWeight();
        
        VARIABLES
        in Inlet                as stream               (Brief = "Inlet Stream", PosX=0, PosY=0.5308, Protected=true,Symbol="_{in}");
        v                               as fraction             (Brief = "Vapourization fraction",Hidden=true);
        x(NComp)                as fraction             (Brief = "Liquid Molar Fraction",Hidden=true);
        y(NComp)                as fraction             (Brief = "Vapour Molar Fraction",Hidden=true);
        zmass(NComp)    as fraction             (Brief = "Mass Fraction");
        Mw                              as molweight    (Brief = "Average Mol Weight");
        vm                              as volume_mol   (Brief = "Molar Volume");       
        rho                             as dens_mass    (Brief = "Stream Mass Density");
        rhom                    as dens_mol             (Brief = "Stream Molar Density");
        Fw                              as flow_mass    (Brief = "Stream Mass Flow");
        Fvol            as flow_vol     (Brief = "Volumetric Flow");
        T_Cdeg                  as temperature  (Brief = "Temperature in °C", Lower=-200);

        EQUATIONS
        "Flash Calculation"
        [v, x, y] = PP.FlashPH(Inlet.P, Inlet.h, Inlet.z);
        
        "Average Molecular Weight"
        Mw = sum(M*Inlet.z);

        "Molar Density"
                rhom * vm = 1;
                
        "Mass or Molar Density"
        rhom * Mw = rho;

        "Flow Mass"
        Fw      =  Mw*Inlet.F;

        "Molar Volume"
        vm = (1-v)*PP.LiquidVolume(Inlet.T, Inlet.P, x) + v*PP.VapourVolume(Inlet.T,Inlet.P,y);
        
        "Volumetric Flow"
        Fvol = Inlet.F*vm ;
        
        "Mass Fraction"
        zmass = M*Inlet.z / Mw;
        
        "Temperature in °C"
        T_Cdeg = Inlet.T - 273.15 * 'K';

end

Model sink2
        ATTRIBUTES
        Pallete = true;
        Icon = "icon/Sink2";
        Brief = "Material stream sink";
        Info = "
        This model should be used for boundary streams when additional
        information about the stream is desired.

        Some of the additional informations calculated by this models are:
         * Mass density
         * Mass flow
         * Mass compostions
         * Specific volume
         * Vapour fraction
         * Volumetric flow
         * Liquid and Vapour compositions
        ";

        PARAMETERS
        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1); 
                  M(NComp)      as molweight    (Brief = "Component Mol Weight");
        
        SET

        M   = PP.MolecularWeight();
        
        VARIABLES
        in Inlet                as stream               (Brief = "Inlet Stream", PosX=1, PosY=0.5308, Protected=true,Symbol="_{in}");
        v                               as fraction             (Brief = "Vapourization fraction",Hidden=true);
        x(NComp)                as fraction             (Brief = "Liquid Molar Fraction",Hidden=true);
        y(NComp)                as fraction             (Brief = "Vapour Molar Fraction",Hidden=true);
        zmass(NComp)    as fraction             (Brief = "Mass Fraction");
        Mw                              as molweight    (Brief = "Average Mol Weight");
        vm                              as volume_mol   (Brief = "Molar Volume");       
        rho                             as dens_mass    (Brief = "Stream Mass Density");
        rhom                    as dens_mol             (Brief = "Stream Molar Density");
        Fw                              as flow_mass    (Brief = "Stream Mass Flow");
        Fvol            as flow_vol     (Brief = "Volumetric Flow");
        T_Cdeg                  as temperature  (Brief = "Temperature in °C", Lower=-200);

        EQUATIONS
        "Flash Calculation"
        [v, x, y] = PP.FlashPH(Inlet.P, Inlet.h, Inlet.z);
        
        "Average Molecular Weight"
        Mw = sum(M*Inlet.z);

        "Molar Density"
                rhom * vm = 1;
                
        "Mass or Molar Density"
        rhom * Mw = rho;

        "Flow Mass"
        Fw      =  Mw*Inlet.F;

        "Molar Volume"
        vm = (1-v)*PP.LiquidVolume(Inlet.T, Inlet.P, x) + v*PP.VapourVolume(Inlet.T,Inlet.P,y);
        
        "Volumetric Flow"
        Fvol = Inlet.F*vm ;
        
        "Mass Fraction"
        zmass = M*Inlet.z / Mw;
        
        "Temperature in °C"
        T_Cdeg = Inlet.T - 273.15 * 'K';

end

Model simple_sink
        ATTRIBUTES
        Pallete = true;
        Icon = "icon/Sink";
        Brief = "Simple material stream sink";
        Info = "
        This model should be used for boundary streams when no additional
        information about the stream is desired.
        ";
        
        VARIABLES
        in Inlet                as stream       (Brief = "Inlet Stream", PosX=0, PosY=0.5308, Protected=true,Symbol="_{in}");
end

Model simple_sink2
        ATTRIBUTES
        Pallete = true;
        Icon = "icon/Sink2";
        Brief = "Simple material stream sink";
        Info = "
        This model should be used for boundary streams when no additional
        information about the stream is desired.
        ";
        
        VARIABLES
        in Inlet                as stream       (Brief = "Inlet Stream", PosX=1, PosY=0.5308, Protected=true,Symbol="_{in}");
        
end

Model sinkNoFlow
        ATTRIBUTES
        Pallete = true;
        Icon = "icon/SinkNoFlow";
        Brief = "Simple material stream sink";
        Info = "
        This model should be used for seal an outlet material stream port.
        ";
        
        VARIABLES
        in Inlet                as stream       (Brief = "Inlet Stream", PosX=0, PosY=0.5308, Protected=true,Symbol="_{in}");

EQUATIONS
"Stream Molar Flow"
        Inlet.F = 0 * 'kmol/h';
        
end