Changeset 372 for trunk/eml/stage_separators/flash.mso

Timestamp:

Sep 21, 2007, 7:17:36 PM (15 years ago)

Author:

Argimiro Resende Secchi

Message:

Adding pressure drop variable (Pdrop) and pressure ratio (Pratio) in some models. Need to do the same for other models.

File:

• trunk/eml/stage_separators/flash.mso (modified) (8 diffs)

trunk/eml/stage_separators/flash.mso

@@ -64 +64 @@
         Across as area (Brief="Flash Cross section area");
         vfrac as positive (Brief="Vapourization fraction", Symbol="\phi");
+        Pratio as positive      (Brief = "Pressure Ratio", Symbol ="P_{ratio}");        
+        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
 
         EQUATIONS
@@ -102 +104 @@
         "Mechanical Equilibrium"
         OutletV.P = OutletL.P;
-        
+
+        "Pressure Drop"
+        OutletL.P  = Inlet.P - Pdrop;
+
+        "Pressure Ratio"
+        OutletL.P = Inlet.P * Pratio;
+
         "Geometry Constraint"
         V = ML * vL + MV * vV;
@@ -151 +159 @@
 in      InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7559, Symbol="_{in}");
         vfrac as fraction (Brief="Vapourization fraction", Symbol="\phi");
+        Pratio as positive      (Brief = "Pressure Ratio", Symbol ="P_{ratio}");        
+        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
 
         EQUATIONS
@@ -170 +180 @@
         "Mechanical Equilibrium"
         OutletV.P = OutletL.P;
+
+        "Pressure Drop"
+        OutletL.P  = Inlet.P - Pdrop;
+
+        "Pressure Ratio"
+        OutletL.P = Inlet.P * Pratio;
 end
 
@@ -203 +219 @@
         vfrac as fraction(Brief="Vaporization fraction", Symbol="\phi");
         h as enth_mol(Brief="Mixture enthalpy");
+        Pratio as positive (Brief = "Pressure Ratio", Symbol ="P_{ratio}");     
+        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
 
         EQUATIONS
@@ -222 +240 @@
         "Mechanical Equilibrium"
         OutletV.P = OutletL.P;
+
+        "Pressure Drop"
+        OutletL.P  = Inlet.P - Pdrop;
+
+        "Pressure Ratio"
+        OutletL.P = Inlet.P * Pratio;
 end
 
@@ -264 +288 @@
         xsat(NComp) as Real(Lower=0, Upper=1, Brief="Liquid composition if saturated");
         ysat(NComp) as Real(Lower=0, Upper=1, Brief="Vapour composition if saturated");
+        Pratio as positive (Brief = "Pressure Ratio", Symbol ="P_{ratio}");     
+        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
         
         zero_one as fraction(Brief="Regularization Variable");
@@ -295 +321 @@
         "Mechanical Equilibrium"
         OutletV.P = OutletL.P;
-        
+
+        "Pressure Drop"
+        OutletL.P  = Inlet.P - Pdrop;
+
+        "Pressure Ratio"
+        OutletL.P = Inlet.P * Pratio;
+
         # regularization functions
         zero_one = (1 + tanh(B * vsat))/2;