Changeset 613 for branches


Ignore:
Timestamp:
Sep 1, 2008, 7:45:47 PM (14 years ago)
Author:
gerson bicca
Message:

updated column model

Location:
branches/gui
Files:
4 edited

Legend:

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

    r612 r613  
    6666
    6767PARAMETERS
    68         outer PP                        as Plugin       (Brief = "External Physical Properties", Type="PP");
    69         outer NComp     as Integer;
    70         NTrays                          as Integer      (Brief="Number of trays", Default=2);
    71         FeedTray                        as Integer      (Brief="Number of Feed tray", Default=2);
    72         topdown                         as Integer      (Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1,Hidden=true);
    73         top                                     as Integer      (Brief="Number of top tray",Hidden=true);
    74         bot                                     as Integer      (Brief="Number of bottom tray",Hidden=true);
    75         g                                               as acceleration         (Default=9.81,Hidden=true);
     68        outer PP                                                        as Plugin       (Brief = "External Physical Properties", Type="PP");
     69        outer NComp                                     as Integer;
     70        NTrays                                                          as Integer      (Brief="Number of trays", Default=2);
     71        SpecialTrayIndex(NTrays)        as Integer      (Brief="Number of trays", Default=0,Hidden=true);
     72        FeedTray                                                        as Integer      (Brief="Number of Feed tray", Default=2);
     73        topdown                                                         as Integer      (Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1,Hidden=true);
     74        top                                                                     as Integer      (Brief="Number of top tray",Hidden=true);
     75        bot                                                                     as Integer      (Brief="Number of bottom tray",Hidden=true);
     76        g                                                                               as acceleration         (Default=9.81,Hidden=true);
     77        Mw(NComp)                                               as molweight    (Brief="Component Mol Weight");
     78       
     79
     80        Ah                              as area                                 (Brief="Total holes area");
     81        lw                              as length                       (Brief="Weir length");
     82        hw                              as length                       (Brief="Weir height");
     83        beta                            as fraction             (Brief="Aeration fraction");
     84        alfa                            as fraction             (Brief="Dry pressure drop coefficient");
     85        w                                       as Real                         (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1);
     86        btray                   as Real                         (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1);
     87        fw                              as Real                         (Brief="Olsens correlation coefficient", Default=1);
     88        Np                              as Real                         (Brief="Number of liquid passes in the tray", Default=1);
    7689       
    7790        VapourFlow                              as Switcher     (Valid = ["on", "off"], Default = "on",Hidden=true);
     
    8396        top = (NTrays-1)*(1-topdown)/2+1;
    8497        bot = NTrays/top;
     98        SpecialTrayIndex(FeedTray) =1;
     99        Mw = PP.MolecularWeight();
    85100
    86101VARIABLES
    87         zero                    as stream                               (Brief="Dummy stream",Hidden=true);
    88102        in      Inlet           as stream                               (Brief="Feed stream", PosX=0, PosY=0.50);
    89103        trays(NTrays) as trayTeste;
     
    94108        out     Lout            as liquid_stream                (Brief="Feed stream", PosX=0.75, PosY=1);
    95109        out     Vout            as vapour_stream        (Brief="Feed stream", PosX=0.55, PosY=0);
     110       
     111        LiquidConnector         as stream               (Brief="Feed stream", PosX=0.75, PosY=1,Hidden=true);
     112        VapourConnector as stream       (Brief="Feed stream", PosX=0.55, PosY=0,Hidden=true);
    96113
    97114CONNECTIONS
     115
    98116        trays([top+topdown:topdown:bot]).OutletV        to trays([top:topdown:bot-topdown]).InletV;
    99117        trays([top:topdown:bot-topdown]).OutletL        to trays([top+topdown:topdown:bot]).InletL;
    100118       
     119        LiquidConnector to trays(1).InletL;
     120        VapourConnector to trays(NTrays).InletV;
     121       
    101122EQUATIONS
    102123
    103 # Zero stream definition
    104         zero.F = 0 * 'kmol/h';
    105         zero.T = 300 * 'K';
    106         zero.P = 1 * 'atm';
    107         zero.z =  1/NComp;
    108         zero.v = 0;
    109         zero.h = 0 * 'J/mol';
    110 
    111 # Connecting Feed Tray
    112         Inlet.F = trays(FeedTray).Inlet.F;
    113         Inlet.T = trays(FeedTray).Inlet.T;
    114         Inlet.P = trays(FeedTray).Inlet.P;
    115         Inlet.z = trays(FeedTray).Inlet.z;
    116         Inlet.v = trays(FeedTray).Inlet.v;
    117         Inlet.h = trays(FeedTray).Inlet.h;
    118 
    119 if FeedTray equal 1
    120 
    121         then
    122 
    123 for i in [2:NTrays] do
    124        
    125         zero.F = trays(i).Inlet.F;
    126         zero.T = trays(i).Inlet.T;
    127         zero.P = trays(i).Inlet.P;
    128         zero.z = trays(i).Inlet.z;
    129         zero.v = trays(i).Inlet.v;
    130         zero.h = trays(i).Inlet.h;
    131 
    132 end
    133 
    134         else
    135 
    136 if FeedTray equal NTrays
    137        
    138         then
    139        
    140 for i in [1:NTrays-1] do
    141 
    142         zero.F = trays(i).Inlet.F;
    143         zero.T = trays(i).Inlet.T;
    144         zero.P = trays(i).Inlet.P;
    145         zero.z = trays(i).Inlet.z;
    146         zero.v = trays(i).Inlet.v;
    147         zero.h = trays(i).Inlet.h;
    148 
    149 end
    150 
    151         else
    152 
    153 for i in [2:FeedTray-1] do
    154 
    155         zero.F = trays(i).Inlet.F;
    156         zero.T = trays(i).Inlet.T;
    157         zero.P = trays(i).Inlet.P;
    158         zero.z = trays(i).Inlet.z;
    159         zero.v = trays(i).Inlet.v;
    160         zero.h = trays(i).Inlet.h;
    161 
    162 end
    163 
    164 for i in [FeedTray+1: NTrays-1] do
    165 
    166         zero.F = trays(i).Inlet.F;
    167         zero.T = trays(i).Inlet.T;
    168         zero.P = trays(i).Inlet.P;
    169         zero.z = trays(i).Inlet.z;
    170         zero.v = trays(i).Inlet.v;
    171         zero.h = trays(i).Inlet.h;
    172 
    173 end
    174 
    175 end
    176 
    177 end
     124        LiquidConnector.F= Lin.F;
     125        LiquidConnector.T = Lin.T;
     126        LiquidConnector.P = Lin.P;
     127        LiquidConnector.z = Lin.z;
     128        LiquidConnector.v = Lin.v;
     129        LiquidConnector.h = Lin.h;
     130       
     131        VapourConnector.F= Vin.F;
     132        VapourConnector.T = Vin.T;
     133        VapourConnector.P = Vin.P;
     134        VapourConnector.z = Vin.z;
     135        VapourConnector.v = Vin.v;
     136        VapourConnector.h = Vin.h;
     137       
     138        Lout.F= trays(NTrays).OutletL.F;
     139        Lout.T = trays(NTrays).OutletL.T;
     140        Lout.P = trays(NTrays).OutletL.P;
     141        Lout.z = trays(NTrays).OutletL.z;
     142        #Lout.v = trays(NTrays).OutletL.v;
     143#       Lout.h = trays(NTrays).OutletL.h;
     144       
     145        Vout.F= trays(1).OutletV.F;
     146        Vout.T = trays(1).OutletV.T;
     147        Vout.P = trays(1).OutletV.P;
     148        Vout.z = trays(1).OutletV.z;
     149        #Vout.v = trays(1).OutletV.v;
     150        #Vout.h = trays(1).OutletV.h;
     151
     152# Connecting Trays
     153        Inlet.F*SpecialTrayIndex= trays.Inlet.F;
     154        Inlet.T = trays.Inlet.T;
     155        Inlet.P = trays.Inlet.P;
     156        Inlet.z = trays.Inlet.z;
     157        Inlet.v = trays.Inlet.v;
     158        Inlet.h = trays.Inlet.h;
    178159
    179160for i in [1:NTrays] do
     
    184165                                case "default":
    185166                                "Francis Equation"
    186                                 trays(i).OutletL.F*trays(i).vL = 1.84*'1/s'*trays(i).lw*((trays(i).Level-(trays(i).beta*trays(i).hw))/(trays(i).beta))^2;
     167                                trays(i).OutletL.F*trays(i).vL = 1.84*'1/s'*lw*((trays(i).Level-(beta*hw))/(beta))^2;
    187168                       
    188169                                case "Wang_Fl":
    189                                 trays(i).OutletL.F*trays(i).vL = 1.84*'m^0.5/s'*trays(i).lw*((trays(i).Level-(trays(i).beta*trays(i).hw))/(trays(i).beta))^1.5;
     170                                trays(i).OutletL.F*trays(i).vL = 1.84*'m^0.5/s'*lw*((trays(i).Level-(beta*hw))/(beta))^1.5;
    190171                       
    191172                                case "Olsen":
    192                                 trays(i).OutletL.F / 'mol/s'= trays(i).lw*trays(i).Np*trays(i).rhoL/sum(trays(i).Mw*trays(i).OutletV.z)/(0.665*trays(i).fw)^1.5 * ((trays(i).ML*sum(trays(i).Mw*trays(i).OutletL.z)/trays(i).rhoL/trays(i).Ap)-trays(i).hw)^1.5 * 'm^0.5/mol';
     173                                trays(i).OutletL.F / 'mol/s'= lw*Np*trays(i).rhoL/sum(Mw*trays(i).OutletV.z)/(0.665*fw)^1.5 * ((trays(i).ML*sum(Mw*trays(i).OutletL.z)/trays(i).rhoL/trays(i).Ap)-hw)^1.5 * 'm^0.5/mol';
    193174                       
    194175                                case "Feehery_Fl":
    195                                 trays(i).OutletL.F = trays(i).lw*trays(i).rhoL/sum(trays(i).Mw*trays(i).OutletL.z) * ((trays(i).Level-trays(i).hw)/750/'mm')^1.5 * 'm^2/s';
     176                                trays(i).OutletL.F = lw*trays(i).rhoL/sum(Mw*trays(i).OutletL.z) * ((trays(i).Level-hw)/750/'mm')^1.5 * 'm^2/s';
    196177                       
    197178                                case "Roffel_Fl":
    198                                 trays(i).OutletL.F = 2/3*sqrt(2*g)*trays(i).rhoL/sum(trays(i).Mw*trays(i).OutletL.z)*trays(i).lw*(2*trays(i).btemp-1)*(trays(i).ML*sum(trays(i).Mw*trays(i).OutletL.z)/(trays(i).Ap*1.3)/trays(i).rhoL/(2*trays(i).btemp-1))^1.5;
     179                                trays(i).OutletL.F = 2/3*sqrt(2*g)*trays(i).rhoL/sum(Mw*trays(i).OutletL.z)*lw*(2*trays(i).btemp-1)*(trays(i).ML*sum(Mw*trays(i).OutletL.z)/(trays(i).Ap*1.3)/trays(i).rhoL/(2*trays(i).btemp-1))^1.5;
    199180                        end
    200                 when trays(i).Level < (trays(i).beta * trays(i).hw) switchto "off";
     181                when trays(i).Level < (beta *hw) switchto "off";
    201182               
    202183                case "off":
    203184                "Low level"
    204185                trays(i).OutletL.F = 0 * 'mol/h';
    205                 when trays(i).Level > (trays(i).beta * trays(i).hw) + 1e-6*'m' switchto "on";
     186                when trays(i).Level > (beta * hw) + 1e-6*'m' switchto "on";
    206187        end
    207188       
    208         trays(i).btemp = 1 - 0.3593/'Pa^0.0888545'*abs(trays(i).OutletV.F*sum(trays(i).Mw*trays(i).OutletV.z)/(trays(i).Ap*1.3)/sqrt(trays(i).rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
     189        trays(i).btemp = 1 - 0.3593/'Pa^0.0888545'*abs(trays(i).OutletV.F*sum(Mw*trays(i).OutletV.z)/(trays(i).Ap*1.3)/sqrt(trays(i).rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
    209190
    210191switch VapourFlow
     
    212193                        switch VapourFlowModel
    213194                                case "Reepmeyer":
    214                                 trays(i).InletV.F*trays(i).vV = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/(trays(i).rhoV*trays(i).alfa))*trays(i).Ah;
     195                                trays(i).InletV.F*trays(i).vV = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/(trays(i).rhoV*alfa))*Ah;
    215196                       
    216197                                case "Feehery_Fv":
    217                                 trays(i).InletV.F = trays(i).rhoV/trays(i).Ap/trays(i).w/sum(trays(i).Mw*trays(i).OutletV.z) * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-(trays(i).rhoV*g*trays(i).ML*trays(i).vL/trays(i).Ap))/trays(i).rhoV);
     198                                trays(i).InletV.F = trays(i).rhoV/trays(i).Ap/w/sum(Mw*trays(i).OutletV.z) * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-(trays(i).rhoV*g*trays(i).ML*trays(i).vL/trays(i).Ap))/trays(i).rhoV);
    218199                       
    219200                                case "Roffel_Fv":
    220                                 trays(i).InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (trays(i).InletV.P - trays(i).OutletV.P)*1e5 - (trays(i).beta*sum(trays(i).M*trays(i).Mw)/(trays(i).Ap*1.3)*g*1e5) * (trays(i).rhoV*trays(i).Ah/sum(trays(i).Mw*trays(i).OutletV.z))^1.08 * 'm^1.08/mol^1.08';
     201                                trays(i).InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (trays(i).InletV.P - trays(i).OutletV.P)*1e5 - (beta*sum(trays(i).M*Mw)/(trays(i).Ap*1.3)*g*1e5) * (trays(i).rhoV*Ah/sum(Mw*trays(i).OutletV.z))^1.08 * 'm^1.08/mol^1.08';
    221202                       
    222203                                case "Klingberg":
     
    224205                       
    225206                                case "Wang_Fv":
    226                                 trays(i).InletV.F * trays(i).vV = trays(i).Ap * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-trays(i).rhoL*g*trays(i).Level)/trays(i).rhoV*trays(i).alfa);
     207                                trays(i).InletV.F * trays(i).vV = trays(i).Ap * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-trays(i).rhoL*g*trays(i).Level)/trays(i).rhoV*alfa);
    227208                               
    228209                                case "Elgue":
    229                                 trays(i).InletV.F  = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/trays(i).btray);
     210                                trays(i).InletV.F  = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/btray);
    230211                        end
    231212                when trays(i).InletV.F < 1e-6 * 'kmol/h' switchto "off";
     
    240221end
    241222
     223
    242224#*----------------------------------------------------------------------
    243225* Model of a  distillation column containing:
     
    247229*       - a splitter which separate reflux and distillate;
    248230*       - a pump in reflux stream;
    249 *---------------------------------------------------------------------*#
     231*---------------------------------------------------------------------
    250232Model Distillation_kettle_condTeste
    251233        ATTRIBUTES
     
    339321*       - a pump in reflux stream.
    340322*
    341 * ------------------------------------------------------------------*#
     323* ------------------------------------------------------------------
    342324Model Distillation_thermosyphon_subcoolingTeste
    343325        ATTRIBUTES
     
    435417*       - a pump in reflux stream.
    436418*
    437 * ------------------------------------------------------------------*#
     419* ------------------------------------------------------------------*
    438420Model Distillation_thermosyphon_condTeste
    439421        ATTRIBUTES
     
    526508*       - a pump in reflux stream.
    527509*
    528 * ------------------------------------------------------------------*#
     510* ------------------------------------------------------------------*
    529511Model Distillation_kettle_subcoolingTeste
    530512        ATTRIBUTES
     
    613595*       - a splitter which separate reflux and distillate;
    614596*       - a pump in reflux stream;
    615 *---------------------------------------------------------------------*#
     597*---------------------------------------------------------------------*
    616598Model RectifierTeste
    617599        ATTRIBUTES
     
    693675*       - a pump in reflux stream.
    694676*
    695 * ------------------------------------------------------------------*#
     677* ------------------------------------------------------------------*
    696678Model Rectifier_subcoolingTeste
    697679        ATTRIBUTES
     
    773755*       - a splitter which separate reflux and distillate;
    774756*       - a pump in reflux stream;
    775 *---------------------------------------------------------------------*#
     757*---------------------------------------------------------------------*
    776758Model Refluxed_StrippingTeste
    777759        ATTRIBUTES
     
    853835*       - a pump in reflux stream.
    854836*
    855 * ------------------------------------------------------------------*#
     837* ------------------------------------------------------------------*
    856838Model Refluxed_Stripping_subcoolingTeste
    857839        ATTRIBUTES
     
    933915*       - a splitter which separate reflux and distillate;
    934916*       - a pump in reflux stream;
    935 *---------------------------------------------------------------------*#
     917*---------------------------------------------------------------------*
    936918Model Refluxed_AbsorptionTeste
    937919        ATTRIBUTES
     
    1014996*       - a pump in reflux stream.
    1015997*
    1016 * ------------------------------------------------------------------*#
     998* ------------------------------------------------------------------*
    1017999Model Refluxed_Absorption_subcoolingTeste
    10181000        ATTRIBUTES
     
    10941076*       - a kettle reboiler;
    10951077*
    1096 * ------------------------------------------------------------------*#
     1078* ------------------------------------------------------------------*
    10971079Model Reboiled_Stripping_kettleTeste
    10981080        ATTRIBUTES
     
    11541136*       - steady state reboiler (thermosyphon);
    11551137*
    1156 * ------------------------------------------------------------------*#
     1138* ------------------------------------------------------------------*
    11571139Model Reboiled_Stripping_thermosyphonTeste
    11581140        ATTRIBUTES
     
    12171199*       - a kettle reboiler;
    12181200*
    1219 * ------------------------------------------------------------------*#
     1201* ------------------------------------------------------------------*
    12201202Model Reboiled_Absorption_kettleTeste
    12211203        ATTRIBUTES
     
    12771259*       - steady state reboiler (thermosyphon);
    12781260*
    1279 * ------------------------------------------------------------------*#
     1261* ------------------------------------------------------------------*
    12801262Model Reboiled_Absorption_thermosyphonTeste
    12811263        ATTRIBUTES
     
    13361318*  Reactive Distillation Column
    13371319*
    1338 * ------------------------------------------------------------------*#
     1320* ------------------------------------------------------------------*
    13391321Model ReactiveDistillationTeste
    13401322        ATTRIBUTES
     
    14231405*       - NStages = theoretical number of equilibrium stages.
    14241406*
    1425 *---------------------------------------------------------------------*#
     1407*---------------------------------------------------------------------*
    14261408Model Packed_Section_ColumnTeste
    14271409        ATTRIBUTES
     
    14791461*       - a splitter which separate reflux and distillate;
    14801462*       - a pump in reflux stream;
    1481 *---------------------------------------------------------------------*#
     1463*---------------------------------------------------------------------*
    14821464Model PackedDistillation_kettle_cond as Packed_Section_ColumnTeste
    14831465        ATTRIBUTES
     
    15421524
    15431525end
     1526
  • branches/gui/eml/stage_separators/tray.mso

    r522 r613  
    197197        end
    198198       
    199         btemp = 1 - 0.3593/'Pa^0.0888545'*(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
     199        btemp = 1 - 0.3593/'Pa^0.0888545'*abs(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
    200200
    201201        switch VapourFlow
  • branches/gui/eml/stage_separators/trayTeste.mso

    r612 r613  
    137137";     
    138138
    139 PARAMETERS
    140         Ah                              as area                                 (Brief="Total holes area");
    141         lw                              as length                       (Brief="Weir length");
    142         hw                              as length                       (Brief="Weir height");
    143         beta                            as fraction             (Brief="Aeration fraction");
    144         alfa                            as fraction             (Brief="Dry pressure drop coefficient");
    145         w                                       as Real                         (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1);
    146         btray                   as Real                         (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1);
    147         fw                              as Real                         (Brief="Olsens correlation coefficient", Default=1);
    148         Np                              as Real                         (Brief="Number of liquid passes in the tray", Default=1);
    149         Mw(NComp)       as molweight    (Brief="Component Mol Weight");
    150        
    151 #       VapourFlow                              as Switcher     (Valid = ["on", "off"], Default = "on",Hidden=true);
    152 #       LiquidFlow                              as Switcher     (Valid = ["on", "off"], Default = "on",Hidden=true);
    153 #       VapourFlowModel         as Switcher     (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer");
    154 #       LiquidFlowModel                 as Switcher     (Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default");
    155 
    156 SET
    157         Mw = PP.MolecularWeight();
    158        
    159139VARIABLES
    160140        rhoL as dens_mass;
     
    171151        rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
    172152
    173 #*      switch LiquidFlow
    174                 case "on":
    175                         switch LiquidFlowModel
    176                                 case "default":
    177                                 "Francis Equation"
    178                                 OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2;
    179                        
    180                                 case "Wang_Fl":
    181                                 OutletL.F*vL = 1.84*'m^0.5/s'*lw*((Level-(beta*hw))/(beta))^1.5;
    182                        
    183                                 case "Olsen":
    184                                 OutletL.F / 'mol/s'= lw*Np*rhoL/sum(Mw*OutletV.z)/(0.665*fw)^1.5 * ((ML*sum(Mw*OutletL.z)/rhoL/Ap)-hw)^1.5 * 'm^0.5/mol';
    185                        
    186                                 case "Feehery_Fl":
    187                                 OutletL.F = lw*rhoL/sum(Mw*OutletL.z) * ((Level-hw)/750/'mm')^1.5 * 'm^2/s';
    188                        
    189                                 case "Roffel_Fl":
    190                                 OutletL.F = 2/3*sqrt(2*g)*rhoL/sum(Mw*OutletL.z)*lw*(2*btemp-1)*(ML*sum(Mw*OutletL.z)/(Ap*1.3)/rhoL/(2*btemp-1))^1.5;
    191                         end
    192                 when Level < (beta * hw) switchto "off";
    193                
    194                 case "off":
    195                 "Low level"
    196                 OutletL.F = 0 * 'mol/h';
    197                 when Level > (beta * hw) + 1e-6*'m' switchto "on";
    198         end
    199        
    200         btemp = 1 - 0.3593/'Pa^0.0888545'*abs(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
    201 
    202         switch VapourFlow
    203                 case "on":
    204                         switch VapourFlowModel
    205                                 case "Reepmeyer":
    206                                 InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;
    207                        
    208                                 case "Feehery_Fv":
    209                                 InletV.F = rhoV/Ap/w/sum(Mw*OutletV.z) * sqrt(((InletV.P - OutletV.P)-(rhoV*g*ML*vL/Ap))/rhoV);
    210                        
    211                                 case "Roffel_Fv":
    212                                 InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (InletV.P - OutletV.P)*1e5 - (beta*sum(M*Mw)/(Ap*1.3)*g*1e5) * (rhoV*Ah/sum(Mw*OutletV.z))^1.08 * 'm^1.08/mol^1.08';
    213                        
    214                                 case "Klingberg":
    215                                 InletV.F * vV = Ap * sqrt(((InletV.P - OutletV.P)-rhoL*g*Level)/rhoV);
    216                        
    217                                 case "Wang_Fv":
    218                                 InletV.F * vV = Ap * sqrt(((InletV.P - OutletV.P)-rhoL*g*Level)/rhoV*alfa);
    219                                
    220                                 case "Elgue":
    221                                 InletV.F  = sqrt((InletV.P - OutletV.P)/btray);
    222                         end
    223                 when InletV.F < 1e-6 * 'kmol/h' switchto "off";
    224                
    225                 case "off":
    226                 InletV.F = 0 * 'mol/s';
    227                 when InletV.P > OutletV.P + Level*g*rhoL + 1e-1 * 'atm' switchto "on";
    228         end
    229 *#
    230153end
    231154
  • branches/gui/sample/stage_separators/sample_columnTeste.mso

    r612 r613  
    2828using "stage_separators/columnTeste";
    2929
    30 # column section with 2 trays
    3130FlowSheet SectionColumn_Test_with2tray
    32        
     31
    3332PARAMETERS
    3433        PP      as Plugin(Brief="Physical Properties",
     
    4140
    4241DEVICES
    43         sec             as Section_ColumnTeste;
    44         feed            as liquid_stream;
    45         reb             as vapour_stream;
    46         cond    as liquid_stream;
    47 
     42        sec     as Section_ColumnTeste;
     43        feed    as liquid_stream;
     44        reb     as vapour_stream;
     45        cond as liquid_stream;
     46       
    4847SET
    4948        NComp = PP.NumberOfComponents;
    50         sec.FeedTray = 2;
     49        sec.FeedTray=1;
     50        sec.NTrays = 2;
     51        #COLUMN
     52        sec.trays.V = 4 * 'ft^3';
     53        sec.Ah = 0.394 * 'ft^2';
     54        sec.lw = 20.94 * 'in';
     55        sec.hw = 0.125 * 'ft';
     56        sec.trays.Q = 0 * 'kW';
     57        sec.beta = 0.6;
     58        sec.alfa = 4;
     59        sec.trays.Ap = 3.94 * 'ft^2';
    5160       
    5261CONNECTIONS
    53         feed            to sec.Inlet;#Feed at stage 2
    54 
     62        feed to sec.Inlet;
     63        reb to sec.Vin;
     64        cond to sec.Lin;
     65       
    5566SPECIFY
    5667        feed.F = 113.4 * 'kmol/h';
     
    7182        sec.trays(1).OutletV.F = 150 * 'kmol/h';
    7283
    73 SET
    74         sec.NTrays = 2;
    75 #COLUMN
    76         sec.trays.V = 4 * 'ft^3';
    77         sec.trays.Ah = 0.394 * 'ft^2';
    78         sec.trays.lw = 20.94 * 'in';
    79         sec.trays.hw = 0.125 * 'ft';
    80         sec.trays.Q = 0 * 'kW';
    81         sec.trays.beta = 0.6;
    82         sec.trays.alfa = 4;
    83         sec.trays.Ap = 3.94 * 'ft^2';
    84 
    8584INITIAL
    8685        sec.trays.OutletL.T = 270 *'K';
    87         sec.trays.Level = 0.9 * sec.trays.hw;
     86        sec.trays.Level = 0.9 * sec.hw;
    8887        sec.trays.OutletL.z(1) = 0.5;
    8988       
    90 OPTIONS
     89        OPTIONS
    9190        TimeStep = 10;
    9291        TimeEnd = 1000;
    9392       
    9493end
    95 
Note: See TracChangeset for help on using the changeset viewer.