source: branches/new_gui/eml/streams/streams.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 basic streams
*----------------------------------------------------------------------
* Author: Paula B. Staudt and Rafael de P. Soares
* $Id: streams.mso 757 2009-06-03 20:07:22Z bicca $
*---------------------------------------------------------------------*#

using "types";

Model stream
        ATTRIBUTES
        Pallete = false;
        Brief = "General Material Stream";
        Info =
        "This is the basic building block for the EML models.
        Every model should have input and output streams derived
        from this model.";
        
        PARAMETERS
        outer NComp as Integer (Brief = "Number of chemical components", Lower = 1); 

        VARIABLES
        F as flow_mol                   (Brief = "Stream Molar Flow Rate");
        T as temperature                (Brief = "Stream Temperature");
        P as pressure                   (Brief = "Stream Pressure");
        h as enth_mol                   (Brief = "Stream Enthalpy");
        v as fraction                   (Brief = "Vapourization fraction");
        z(NComp) as fraction    (Brief = "Stream Molar Fraction");
end

Model liquid_stream as stream
        ATTRIBUTES
        Pallete = false;
        Brief = "Liquid Material Stream";
        Info =
        "Model for liquid material streams.
        This model should be used only when the phase of the stream
        is known ''a priori''.";

        PARAMETERS
        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
        
        EQUATIONS
        "Liquid Enthalpy"
        h = PP.LiquidEnthalpy(T, P, z);
        "Liquid stream"
        v = 0;
end

Model vapour_stream as stream
        ATTRIBUTES
        Pallete = false;
        Brief = "Vapour Material Stream";
        Info =
        "Model for vapour material streams.
        This model should be used only when the phase of the stream
        is known ''a priori''.";

        PARAMETERS
        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
        
        EQUATIONS
        "Vapour Enthalpy"
        h = PP.VapourEnthalpy(T, P, z);
        "Vapour stream"
        v = 1;
end

Model streamPH as stream
        ATTRIBUTES
        Brief = "Stream with built-in flash calculation";
        Info = "
        This model should be used when the vaporization fraction
        is unknown.
        
        The built-in flash calculation will determine the stream
        state as a function of the overall composition '''z''', the
        pressure '''P''' and the enthalpy '''h'''.
        
        Additionally, the liquid composition '''x''' and the vapor
        composition '''y''' are calculated.     
        ";
        Pallete = false;
        
        PARAMETERS
        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
        
        VARIABLES
        x(NComp) as fraction    (Brief = "Liquid Molar Fraction",Hidden=true);
        y(NComp) as fraction    (Brief = "Vapour Molar Fraction",Hidden=true);

        EQUATIONS
        "Flash Calculation"
        [v, x, y] = PP.FlashPH(P, h, z);
        
        "Enthalpy"
        h = (1-v)*PP.LiquidEnthalpy(T, P, x) + v*PP.VapourEnthalpy(T, P, y);
        
end

Model streamPHS as streamPH
        ATTRIBUTES
        Brief = "Stream with built-in flash calculation";
        Info = "
        This model should be used when the vaporization fraction
        is unknown.
        
        The built-in flash calculation will determine the stream
        state as a function of the overall composition '''z''', the
        pressure '''P''' and the enthalpy '''h'''.
        
        Additionally, the liquid composition '''x''', the vapor
        composition '''y''' and the stream entropy are calculated.      
        ";
        Pallete = false;
        
PARAMETERS
        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
        
VARIABLES
        s as entr_mol   (Brief = "Stream Entropy");

EQUATIONS

"Entropy"
        s = (1-v)*PP.LiquidEntropy(T, P, x) +   v*PP.VapourEntropy(T, P, y);
        
end