Ignore:
Timestamp:
Jul 24, 2009, 2:03:40 PM (13 years ago)
Author:
gerson bicca
Message:

updates - added thermosyphon reboiler / fixed sump tank

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/stage_separators/reboiler.mso

    r794 r805  
    3030* perfect mixing of both phases;
    3131* no thermodynamics equilibrium;
    32 * no liquid entrainment in the vapour stream.
    3332
    3433== SET ==
     
    5049        outer NComp     as Integer              (Brief="Number of Components");
    5150        Pdrop                   as press_delta  (Brief="Pressure Drop in the reboiler", Symbol = "\Delta P");
    52 
     51        FlowConstant                    as Real (Brief = "Flow Constant");
     52        k                       as Real (Brief = "Flow Constant", Hidden = true, Unit='mol^3/(kg*m^2)');
     53
     54SET
     55
     56        k = 1*'mol^3/(kg*m^2)';
     57       
    5358VARIABLES
    5459        in      InletLiquid     as stream                               (Brief="Liquid inlet stream", PosX=0.345, PosY=1, Symbol="_{inL}");
    55         out     OutletVapour    as vapour_stream                (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
     60        out     OutletVapour    as streamPH                             (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
    5661        in      InletQ                  as power                                (Brief="Heat supplied", PosX=1, PosY=0.08, Symbol="Q_{in}", Protected = true);
    5762        vV                              as volume_mol                   (Brief="Vapour Molar volume", Protected = true);
     
    8691"Pressure indicator"
    8792        PI * 'atm' = OutletVapour.P;
     93       
     94"Flow through the reboiler"
     95        OutletVapour.F^3 = FlowConstant*k*InletQ;
    8896
    8997end
    9098
     99Model thermosyphon
     100
     101ATTRIBUTES
     102        Pallete         = true;
     103        Icon            = "icon/Thermosyphon";
     104        Brief           = "Model of a  Steady State reboiler thermosyphon.";
     105        Info            =
     106"== ASSUMPTIONS ==
     107* perfect mixing of both phases;
     108* no thermodynamics equilibrium;
     109
     110== SET ==
     111* the pressure drop in the reboiler;
     112
     113== SPECIFY ==
     114* the InletLiquid stream;
     115* the InletQ (the model requires an energy stream, also you can use a controller for setting the heat duty using the heat_flow model)
     116OR the outlet temperature (OutletVapour.T);
     117
     118== OPTIONAL ==
     119* the reboiler model has two control ports
     120** TI OutletVapour Temperature Indicator;
     121** PI OutletVapour Pressure Indicator;
     122";     
     123
     124PARAMETERS
     125        outer PP                as Plugin               (Brief = "External Physical Properties", Type="PP");
     126        outer NComp     as Integer              (Brief="Number of Components");
     127        Pdrop                   as press_delta  (Brief="Pressure Drop in the reboiler", Symbol = "\Delta P");
     128        FlowConstant                    as Real (Brief = "Flow Constant");
     129        k                       as Real (Brief = "Flow Constant", Hidden = true, Unit='mol^3/(kg*m^2)');
     130
     131SET
     132
     133        k = 1*'mol^3/(kg*m^2)';
     134       
     135VARIABLES
     136        in      InletLiquid     as stream                               (Brief="Liquid inlet stream", PosX=0.44, PosY=1, Symbol="_{inL}");
     137        out     OutletVapour    as streamPH                             (Brief="Vapour outlet stream", PosX=0, PosY=0.09, Symbol="_{outV}");
     138        in      InletQ                  as power                                (Brief="Heat supplied", PosX=1, PosY=0.77, Symbol="Q_{in}", Protected = true);
     139       
     140        out     TI as control_signal    (Brief="Temperature  Indicator of Reboiler", Protected = true, PosX=1, PosY=0.57);
     141        out     PI as control_signal    (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=1, PosY=0.35);
     142
     143EQUATIONS
     144
     145"Molar Flow Balance"
     146        InletLiquid.F = OutletVapour.F;
     147
     148"Molar Composition Balance"
     149        InletLiquid.z = OutletVapour.z;
     150       
     151"Energy Balance"
     152        InletLiquid.F*InletLiquid.h + InletQ = OutletVapour.F*OutletVapour.h;
     153       
     154"Pressure Drop"
     155        OutletVapour.P = InletLiquid.P - Pdrop;
     156
     157"Temperature indicator"
     158        TI * 'K' = OutletVapour.T;
     159
     160"Pressure indicator"
     161        PI * 'atm' = OutletVapour.P;
     162       
     163"Flow through the thermosyphon reboiler"
     164        OutletVapour.F^3 = FlowConstant*k*InletQ;
     165
     166end
     167
    91168Model reboilerSteady_fakeH
     169
    92170        ATTRIBUTES
    93171        Pallete         = true;
     
    113191VARIABLES
    114192        in      InletLiquid     as stream                               (Brief="Liquid inlet stream", PosX=0.345, PosY=1, Symbol="_{inL}");
    115         out     OutletVapour    as vapour_stream                (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
     193        out     OutletVapour    as stream                               (Brief="Vapour outlet stream", PosX=0.17, PosY=0, Symbol="_{outV}");
    116194        in      InletQ                  as power                                (Brief="Heat Duty", PosX=1, PosY=0.08, Symbol="Q_{in}", Protected = true);
    117195
     
    131209
    132210"Fake Vapourisation Fraction"
    133         OutletVapour.v = 1.0;
    134        
     211        OutletVapour.v = 0.6;
    135212"Fake output temperature"
    136213        OutletVapour.T = 300*'K';
    137214       
    138 "Pressure Drop through the reboiler"
     215        #OutletVapour.v = PP.VapourFraction(OutletVapour.T, OutletVapour.P, OutletVapour.z);
     216        #OutletVapour.h = OutletVapor.v * PP.VapourEnthalpy(OutletVapour.T, OutletVapour.P, OutletVapour.z)
     217        #       +  (1-OutletVapor.v) * PP.LiquidEnthalpy(OutletVapour.T, OutletVapour.P, OutletVapour.z)
     218
     219"Flow through the reboiler"
    139220        OutletVapour.F = FlowConstant*k*InletQ;
    140221
Note: See TracChangeset for help on using the changeset viewer.