Changeset 690

Timestamp:

Nov 21, 2008, 3:02:51 PM (15 years ago)

Author:

jjunges

Message:

valve gas update

File:

• branches/gui/eml/pressure_changers/valve.mso (modified) (7 diffs)

branches/gui/eml/pressure_changers/valve.mso

@@ -1 @@
-#*-------------------------------------------------------------------
+ #*-------------------------------------------------------------------
 * EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
 *
@@ -176 +176 @@
         ValveType       as Switcher             (Valid = ["linear", "parabolic", "equal", "quick", "hyperbolic"], Default = "linear");
         ValidPhases     as Switcher             (Brief = "Valid Phases for Flash Calculation", Valid = ["Vapour-Only", "Liquid-Only"], Default="Liquid-Only");
-        Tau                     as time_sec             (Brief="valve time constant");
+#       Tau                     as time_sec             (Brief="valve time constant");
         rho60F                  as dens_mass            (Brief = "Water Mass Density at 60 F",Hidden=true);     
 
 VARIABLES
-
+        W as flow_mass(DisplayUnit='kg/s');
         Pdrop                   as press_delta          (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
         Fvol                            as flow_vol                     (Brief = "Volumetric Flow");
         fc                              as positive                     (Brief = "Opening Function",Hidden=true);
         Cv                              as Real                                 (Brief="Valve Flow Coefficient", Unit='gal/min/psi^0.5');
+        Cv1                             as Real                                 (Brief="Valve Flow Coefficient", Unit='m^2');
+        Cg                              as Real                                 (Brief="Valve Gas Flow Coefficient", Unit='ft^3/h/psi');
+        C                               as Real                                 (Brief="Liquid-gas Coefficient Ratio", Unit='(ft^3/gal)*(min/h)/(psi^.5)');
         StemPosition    as fraction                     (Brief = "Actual valve stem position");
-        
-
+        a as Real;
+        #b as Real (Brief="d", Unit='1/(psi^.5)');
         vm                              as vol_mol                      (Brief = "Mixture Molar Volume");
         rho                             as dens_mass            (Brief = "Mixture Mass Density");       
-        vsp                             as fraction                     (Brief = "Valve stem position",Hidden=true);
+#       vsp                             as fraction                     (Brief = "Valve stem position",Hidden=true);
         
 in              Inlet   as stream                       (Brief = "Inlet stream", PosX=0, PosY=0.7365, Symbol="_{in}");
 out     Outlet  as streamPH             (Brief = "Outlet stream", PosX=1, PosY=0.7365, Symbol="_{out}");
-in              vsignal as fraction             (Brief = "Flow Signal", PosX=0.5, PosY=0);
+#in             vsignal as fraction             (Brief = "Flow Signal", PosX=0.5, PosY=0);
 
 SET
@@ -202 +205 @@
 EQUATIONS
 
-"First order valve dynamics"
-        Tau*diff(StemPosition) = vsp-StemPosition;
-
-"Flow Signal"
-        vsp = vsignal;
+#"First order valve dynamics"
+#       Tau*diff(StemPosition) = vsp-StemPosition;
+
+#"Flow Signal"
+#       vsp = vsignal;
 
 "Pressure Drop"
@@ -219 +222 @@
 "Outlet Composition"
         Outlet.z = Inlet.z;
-
-switch ValidPhases      
-        
+        
+        Cv1=Cv*'1/(gal/min/psi^0.5)'*2.3837e-5*'m^2';
+
+switch ValidPhases
+                
+#############################################################
+
         case "Liquid-Only":
         
@@ -228 +235 @@
 "Valve Equation - Liquid Flow"
         Fvol = fc*(Cv/sqrt(1/rho60F))*sqrt(Pdrop/rho);
-        
+        "Liquid-gas Coefficient Ratio"
+        C*Cv=Cg;
+a=1/(1.6764e-2*C*'1/((ft^3/gal)*(min/h)/(psi^.5))')*sqrt(Pdrop/Inlet.P);        
 else
 
 "Valve Equation - Liquid Flow"
         Fvol = fc*(Cv/sqrt(1/rho60F))*sqrt(Pdrop/rho);
-
+        "Liquid-gas Coefficient Ratio"
+        C*Cv=Cg;
+a=1/(1.6764e-2*C*'1/((ft^3/gal)*(min/h)/(psi^.5))')*sqrt(Pdrop/Inlet.P); 
 end
         
@@ -241 +252 @@
 "Liquid Molar Volume"
         vm = PP.LiquidVolume(Inlet.T,Inlet.P,Inlet.z);
+        
+
+        
+############################################################
 
         case "Vapour-Only":
         
 if Pdrop > 0 then #Update for gas flow !!!!
+        
+"Liquid-gas Coefficient Ratio"
+        C*Cv=Cg;
+        a=1/(1.6764e-2*C*'1/((ft^3/gal)*(min/h)/(psi^.5))')*sqrt(Pdrop/Inlet.P); 
+
+        if 1.5708 > a then
+        "Valve Equation - Vapour Flow"
+                #Fvol = fc*Cg*sqrt(Inlet.P/1000*rho60F/rho);####rho60f/rho ok!!!!
+                #Fvol = fc*Cv*sqrt(Pdrop/1000*rho60F/rho);
+                #W = fc*Cv1*sqrt(Pdrop/1000*rho);
+                Fvol = fc*0.13446*'psi^.5'*Cg*sqrt(Inlet.P/1000*rho60F/rho)*sin(a*'rad');
+                
+        else
+        "Valve Equation - Vapour Flow"
+                Fvol = fc*0.13446*Cv*sqrt(Inlet.P*rho60F/rho);
+        end
+else
 
 "Valve Equation - Vapour Flow"
         Fvol = fc*(Cv/sqrt(1/rho60F))*sqrt(Pdrop/rho);
-        
-else
-
-"Valve Equation - Vapour Flow"
-        Fvol = fc*(Cv/sqrt(1/rho60F))*sqrt(Pdrop/rho);
-
+
+"Liquid-gas Coefficient Ratio"
+        C*Cv=Cg;
+        a=1/(1.6764e-2*C*'1/((ft^3/gal)*(min/h)/(psi^.5))')*sqrt(Pdrop/Inlet.P); 
 end
         
 "Vapour Mass Density"
         rho = PP.VapourDensity(Inlet.T,Inlet.P,Inlet.z);
-
+        #rho=3.708741*'kg/m^3';
 "Vapour Molar Volume"
         vm = PP.VapourVolume(Inlet.T,Inlet.P,Inlet.z);
@@ -264 +294 @@
 end
 
+######################################################
+
 "Calculate Mass Flow"
-        Fvol = Inlet.F*vm;      
+        Fvol = Inlet.F*vm;
+        
+        W=Fvol*rho;
+        
         
 switch ValveType #Update the valve Type