Changeset 57 for mso/eml/pressure_changers/valve.mso

Timestamp:

Nov 16, 2006, 2:46:48 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

Added new pressure_changers models

File:

• mso/eml/pressure_changers/valve.mso (modified) (2 diffs)

mso/eml/pressure_changers/valve.mso

@@ +1 @@
+#*-----------------------------------------------------------------------------
+* Author:       Estefane Horn                                   Date:06/10/2006
+*                       
+*$Id:$                                                                  
+*-----------------------------------------------------------------------------
+* Assumptions
+* -> Steady State
+* -> Isentalpic
+* -> Liquid
+* 
+*#
+
+using "streams";
+using "pressure_changers/flux_machine_basic";
+        
+
+Model valve_basic as flux_machine_basic_TP
+        
+        PARAMETERS
+ext PP                  as CalcObject   (Brief = "External Physical Properties", File = "vrpp");
+ext NComp       as Integer              (Brief = "Number of chemical components", Lower = 1);
+        rho60F  as dens_mass;
+
+        VARIABLES
+        Pdiff   as press_delta          (Brief = "Pressure Increase", Unit = "kPa");
+        Qv              as flow_vol                     (Brief = "Volumetric Flow");
+        fc              as positive                     (Brief = "Opening Function");
+        cv              as positive                     (Brief = "Valve Coefficient", Unit = "m^3/h/kPa^0.5");
+        Gf              as positive                     (Brief = "Specific Gravity");
+        rho     as dens_mass;   
+        vm              as vol_mol                      (Brief = "Mixture Molar Volume", Unit = "m^3/kmol");    
+        
+        SET
+        rho60F = 99.022 * "kg/m^3";
+        
+        EQUATIONS
+        "Calculate Outlet Stream Pressure"
+        Inlet.P - Outlet.P = Pdiff;
+        
+        "Enthalpy Balance"
+        Outlet.h = Inlet.h;
+        
+        "Molar Balance"
+        Outlet.F = Inlet.F;
+        
+        "Calculate Outlet Composition"
+        Outlet.z = Inlet.z;
+
+        "Valve Equation"
+        Qv = fc*cv*sqrt(Pdiff/Gf);
+        
+        "Calculate Gf"
+        Gf = rho/rho60F;
+        
+        "Calculate Specific Mass"
+        rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
+        
+        "Calculate Mass Flow"
+        Qv = Inlet.F*vm;        
+        
+        "Calculate Liquid Molar Volume"
+        vm = PP.LiquidVolume(Inlet.T,Inlet.P,Inlet.z);
+        
+end
+
+Model valve_linear as valve_basic
+
+        VARIABLES
+        x               as fraction (Brief = "Opening");
+
+        EQUATIONS
+        
+        "Opening Equation"
+        fc = 100*x;
+
+end
+
+Model valve_parabolic as valve_basic
+
+        PARAMETERS
+        n               as positive (Brief = "Constant", Lower = 1.4, Upper = 2.6);
+
+        VARIABLES
+        x               as fraction (Brief = "Opening");
+        
+        EQUATIONS
+        
+        "Opening Equation"
+        fc = 100*x^n;
+
+end
+
+Model valve_equal as valve_basic
+
+        PARAMETERS
+        a               as Real (Brief = "Constant", Default = 100);
+
+        VARIABLES
+        x               as fraction (Brief = "Opening");
+
+        EQUATIONS
+        
+        "Opening Equation"
+        fc = 100*a^(x-1);
+
+end
+
+Model valve_quick as valve_basic
+
+        PARAMETERS
+        a               as positive     (Brief = "Constant", Default = 0.05);
+        n               as positive (Brief = "Constant", Default = 5);
+
+        VARIABLES
+        x               as fraction (Brief = "Opening");
+
+        EQUATIONS
+        
+        "Opening Equation"
+        fc = 100*(1-(a*(1-x)-(a-1)*(1-x)^n));
+
+end
+
 #*-------------------------------------------------------------------
 * Model of a valve
@@ -27 +150 @@
 * $Id$
 *--------------------------------------------------------------------*#
-
-using "streams";
 
 Model valve