Changeset 805 for branches/gui/eml/stage_separators/reboiler.mso

Timestamp:

Jul 24, 2009, 2:03:40 PM (14 years ago)

Author:

gerson bicca

Message:

updates - added thermosyphon reboiler / fixed sump tank

File:

• branches/gui/eml/stage_separators/reboiler.mso (modified) (5 diffs)

branches/gui/eml/stage_separators/reboiler.mso

@@ -30 +30 @@
 * perfect mixing of both phases;
 * no thermodynamics equilibrium;
-* no liquid entrainment in the vapour stream.
 
 == SET ==
@@ -50 +49 @@
         outer NComp     as Integer              (Brief="Number of Components");
         Pdrop                   as press_delta  (Brief="Pressure Drop in the reboiler", Symbol = "\Delta P");
-
+        FlowConstant                    as Real (Brief = "Flow Constant");
+        k                       as Real (Brief = "Flow Constant", Hidden = true, Unit='mol^3/(kg*m^2)');
+
+SET
+
+        k = 1*'mol^3/(kg*m^2)';
+        
 VARIABLES
         in      InletLiquid     as stream                               (Brief="Liquid inlet stream", PosX=0.345, PosY=1, Symbol="_{inL}");
-        out     OutletVapour    as vapour_stream                (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
+        out     OutletVapour    as streamPH                             (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
         in      InletQ                  as power                                (Brief="Heat supplied", PosX=1, PosY=0.08, Symbol="Q_{in}", Protected = true);
         vV                              as volume_mol                   (Brief="Vapour Molar volume", Protected = true);
@@ -86 +91 @@
 "Pressure indicator"
         PI * 'atm' = OutletVapour.P;
+        
+"Flow through the reboiler"
+        OutletVapour.F^3 = FlowConstant*k*InletQ;
 
 end
 
+Model thermosyphon
+
+ATTRIBUTES
+        Pallete         = true;
+        Icon            = "icon/Thermosyphon"; 
+        Brief           = "Model of a  Steady State reboiler thermosyphon.";
+        Info            =
+"== ASSUMPTIONS ==
+* perfect mixing of both phases;
+* no thermodynamics equilibrium;
+
+== SET ==
+* the pressure drop in the reboiler;
+
+== SPECIFY ==
+* the InletLiquid stream;
+* the InletQ (the model requires an energy stream, also you can use a controller for setting the heat duty using the heat_flow model)
+OR the outlet temperature (OutletVapour.T);
+
+== OPTIONAL ==
+* the reboiler model has two control ports
+** TI OutletVapour Temperature Indicator;
+** PI OutletVapour Pressure Indicator;
+";      
+
+PARAMETERS
+        outer PP                as Plugin               (Brief = "External Physical Properties", Type="PP");
+        outer NComp     as Integer              (Brief="Number of Components");
+        Pdrop                   as press_delta  (Brief="Pressure Drop in the reboiler", Symbol = "\Delta P");
+        FlowConstant                    as Real (Brief = "Flow Constant");
+        k                       as Real (Brief = "Flow Constant", Hidden = true, Unit='mol^3/(kg*m^2)');
+
+SET
+
+        k = 1*'mol^3/(kg*m^2)';
+        
+VARIABLES
+        in      InletLiquid     as stream                               (Brief="Liquid inlet stream", PosX=0.44, PosY=1, Symbol="_{inL}");
+        out     OutletVapour    as streamPH                             (Brief="Vapour outlet stream", PosX=0, PosY=0.09, Symbol="_{outV}");
+        in      InletQ                  as power                                (Brief="Heat supplied", PosX=1, PosY=0.77, Symbol="Q_{in}", Protected = true);
+        
+        out     TI as control_signal    (Brief="Temperature  Indicator of Reboiler", Protected = true, PosX=1, PosY=0.57);
+        out     PI as control_signal    (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=1, PosY=0.35);
+
+EQUATIONS
+
+"Molar Flow Balance"
+        InletLiquid.F = OutletVapour.F;
+
+"Molar Composition Balance"
+        InletLiquid.z = OutletVapour.z;
+        
+"Energy Balance"
+        InletLiquid.F*InletLiquid.h + InletQ = OutletVapour.F*OutletVapour.h;
+        
+"Pressure Drop"
+        OutletVapour.P = InletLiquid.P - Pdrop;
+
+"Temperature indicator"
+        TI * 'K' = OutletVapour.T;
+
+"Pressure indicator"
+        PI * 'atm' = OutletVapour.P;
+        
+"Flow through the thermosyphon reboiler"
+        OutletVapour.F^3 = FlowConstant*k*InletQ;
+
+end
+
 Model reboilerSteady_fakeH
+
         ATTRIBUTES
         Pallete         = true;
@@ -113 +191 @@
 VARIABLES
         in      InletLiquid     as stream                               (Brief="Liquid inlet stream", PosX=0.345, PosY=1, Symbol="_{inL}");
-        out     OutletVapour    as vapour_stream                (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
+        out     OutletVapour    as stream                               (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
         in      InletQ                  as power                                (Brief="Heat Duty", PosX=1, PosY=0.08, Symbol="Q_{in}", Protected = true);
 
@@ -131 +209 @@
 
 "Fake Vapourisation Fraction"
-        OutletVapour.v = 1.0;
-        
+        OutletVapour.v = 0.6;
 "Fake output temperature"
         OutletVapour.T = 300*'K';
         
-"Pressure Drop through the reboiler"
+        #OutletVapour.v = PP.VapourFraction(OutletVapour.T, OutletVapour.P, OutletVapour.z);
+        #OutletVapour.h = OutletVapor.v * PP.VapourEnthalpy(OutletVapour.T, OutletVapour.P, OutletVapour.z)
+        #       +  (1-OutletVapor.v) * PP.LiquidEnthalpy(OutletVapour.T, OutletVapour.P, OutletVapour.z)
+
+"Flow through the reboiler"
         OutletVapour.F = FlowConstant*k*InletQ;