Ignore:
Timestamp:
Aug 14, 2009, 5:20:34 PM (13 years ago)
Author:
gerson bicca
Message:

updates

File:
1 edited

Legend:

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

    r829 r831  
    110110        VapourSideStreamLocation                        as Integer                      (Brief="Vapour Side Stream Location", Default=2);
    111111        Pi                                                                      as constant             (Brief="Pi Number",Default=3.14159265, Symbol = "\pi",Hidden=true);
    112         g                                                                       as acceleration         (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
     112        Gconst                                                          as acceleration         (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
    113113        Mw(NComp)                                                       as molweight            (Brief="Component Mol Weight",Hidden=true);
     114        zero_flow                                                       as flow_mol             (Brief = "Stream Flow closed",Default = 0, Hidden=true);
     115        low_flow                                                        as flow_mol             (Brief = "Low stream Flow",Default = 1E-6, Hidden=true);
    114116
    115117        VapourFlowModel         as Switcher     (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer");
    116         LiquidFlowModel         as Switcher     (Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default");
     118        LiquidFlowModel         as Switcher     (Valid = ["Francis", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "Francis");
    117119
    118120        TrayDiameter                    as length               (Brief="Tray Diameter",Default=1.600);
     
    120122        Fraction_HoleArea               as fraction     (Brief="Fraction of the active area that is occupied by the holes with respect to the total tray area",Default=0.10);
    121123        Fraction_DowncomerArea  as fraction     (Brief="Fraction of the downcomer area with respect to the total tray area",Default=0.20);
    122         WeirLength                              as length               (Brief="Weir length");
    123         WeirHeight                              as length               (Brief="Weir height");
    124         TrayLiquidPasses                as Real                 (Brief="Number of liquid passes in the tray", Default=1);
    125         HeatSupply                              as heat_rate    (Brief="Rate of heat supply");
     124        WeirLength                              as length               (Brief="Weir length", Default = 1);
     125        WeirHeight                              as length               (Brief="Weir height", Default= 0.05);
     126        TrayLiquidPasses                as positive     (Brief="Number of liquid passes in the tray", Lower = 1,Default=1);
     127        HeatSupply                              as heat_rate    (Brief="Rate of heat supply",Default = 0);
    126128        AerationFraction                as Real                 (Brief="Aeration fraction", Default = 1);
    127129        DryPdropCoeff                   as Real                 (Brief="Dry pressure drop coefficient", Default= 0.60);
     
    135137        OlsenCoeff              as Real         (Brief="Olsens correlation coefficient", Default=1,Hidden=true);
    136138       
    137         VapourFlow      as Switcher     (Valid = ["on", "off"], Default = "off",Hidden=true);
    138         LiquidFlow      as Switcher     (Valid = ["on", "off"], Default = "off",Hidden=true);
     139        VapourFlow      as Switcher     (Brief="Flag for Vapour Flow condition",Valid = ["on", "off"], Default = "off",Hidden=true);
     140        LiquidFlow      as Switcher     (Brief="Flag for Liquid Flow condition",Valid = ["on", "off"], Default = "off",Hidden=true);
    139141
    140142SET
     
    142144        LiqSideTrayIndex(LiquidSideStreamLocation) =1;
    143145        Mw = PP.MolecularWeight();
     146        zero_flow = 0 * 'kmol/h';
     147        low_flow = 1E-6 * 'kmol/h';
    144148
    145149        PlateArea = 0.25*Pi*(TrayDiameter^2)*(1-Fraction_DowncomerArea);
     
    231235                case "on":
    232236                        switch LiquidFlowModel
    233                                 case "default":
     237                                case "Francis":
    234238                                "Francis Equation"
    235239                                TRAYS(i).OutletLiquid.F*TRAYS(i).vL = 1.84*'1/s'*WeirLength*((TRAYS(i).Level-(AerationFraction*WeirHeight))/(AerationFraction))^2;
     
    245249                       
    246250                                case "Roffel_Fl":
    247                                 TRAYS(i).OutletLiquid.F = 2/3*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletLiquid.z)*WeirLength*(TRAYS(i).ML*sum(Mw*TRAYS(i).OutletLiquid.z)/(PlateArea*1.3)/TRAYS(i).rhoL)^1.5*sqrt(2*g/
     251                                TRAYS(i).OutletLiquid.F = 2/3*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletLiquid.z)*WeirLength*(TRAYS(i).ML*sum(Mw*TRAYS(i).OutletLiquid.z)/(PlateArea*1.3)/TRAYS(i).rhoL)^1.5*sqrt(2*Gconst/
    248252                                                        (2*(1 - 0.3593/'Pa^0.0888545'*abs(TRAYS(i).OutletVapour.F*sum(Mw*TRAYS(i).OutletVapour.z)/(PlateArea*1.3)/sqrt(TRAYS(i).rhoV))^0.177709)-1)); #/'(kg/m)^0.0888545/s^0.177709';
    249253                        end
     
    251255               
    252256                case "off":
     257               
    253258                "Low level"
    254                 TRAYS(i).OutletLiquid.F = 0 * 'mol/h';
     259                TRAYS(i).OutletLiquid.F = zero_flow;
     260               
    255261                when TRAYS(i).Level > (AerationFraction * WeirHeight) + 1e-6*'m' switchto "on";
     262               
    256263        end
    257264       
     
    263270                       
    264271                                case "Feehery_Fv":
    265                                 TRAYS(i).InletVapour.F = TRAYS(i).rhoV/PlateArea/FeeheryCoeff/sum(Mw*TRAYS(i).OutletVapour.z) * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-(TRAYS(i).rhoV*g*TRAYS(i).ML*TRAYS(i).vL/PlateArea))/TRAYS(i).rhoV);
     272                                TRAYS(i).InletVapour.F = TRAYS(i).rhoV/PlateArea/FeeheryCoeff/sum(Mw*TRAYS(i).OutletVapour.z) * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-(TRAYS(i).rhoV*Gconst*TRAYS(i).ML*TRAYS(i).vL/PlateArea))/TRAYS(i).rhoV);
    266273                       
    267274                                case "Roffel_Fv":
    268                                 TRAYS(i).InletVapour.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)*1e5 - (AerationFraction*sum(TRAYS(i).M*Mw)/(PlateArea*1.3)*g*1e5) * (TRAYS(i).rhoV*HolesArea/sum(Mw*TRAYS(i).OutletVapour.z))^1.08 * 'm^1.08/mol^1.08';
     275                                TRAYS(i).InletVapour.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)*1e5 - (AerationFraction*sum(TRAYS(i).M*Mw)/(PlateArea*1.3)*Gconst*1e5) * (TRAYS(i).rhoV*HolesArea/sum(Mw*TRAYS(i).OutletVapour.z))^1.08 * 'm^1.08/mol^1.08';
    269276                       
    270277                                case "Klingberg":
    271                                 TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*g*TRAYS(i).Level)/TRAYS(i).rhoV);
     278                                TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*Gconst*TRAYS(i).Level)/TRAYS(i).rhoV);
    272279                       
    273280                                case "Wang_Fv":
    274                                 TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*g*TRAYS(i).Level)/TRAYS(i).rhoV*DryPdropCoeff);
     281                                TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*Gconst*TRAYS(i).Level)/TRAYS(i).rhoV*DryPdropCoeff);
    275282                               
    276283                                case "Elgue":
    277284                                TRAYS(i).InletVapour.F  = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/ElgueCoeff);
    278285                        end
    279                 when TRAYS(i).InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
     286                when TRAYS(i).InletVapour.F < low_flow switchto "off";
    280287               
    281288                case "off":
    282                 TRAYS(i).InletVapour.F = 0 * 'mol/s';
    283                 when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P + TRAYS(i).Level*g*TRAYS(i).rhoL + 1e-1 * 'atm' switchto "on";
     289                TRAYS(i).InletVapour.F = zero_flow;
     290               
     291                when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P + TRAYS(i).Level*Gconst*TRAYS(i).rhoL + 1e-1 * 'atm' switchto "on";
     292               
    284293        end
    285294
Note: See TracChangeset for help on using the changeset viewer.