Changeset 235


Ignore:
Timestamp:
Apr 11, 2007, 10:16:05 PM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Added orientation type (vertical and horizontal) into the flash model. "Across" is now a variable and "diameter" is a new parameter.

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/stage_separators/flash.mso

    r118 r235  
    4141Model flash
    4242        PARAMETERS
    43         outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
     43        outer PP as Plugin (Brief = "External Physical Properties", Type="PP");
    4444        outer NComp as Integer (Brief = "Number of chemical components", Lower = 1);
    45         V as volume(Brief="Total Volume of the flash");
     45        V as volume (Brief="Total Volume of the flash");
    4646        Mw(NComp) as molweight;
    47         Across as area  (Brief="Flash Cross section area");
    48        
     47        orientation as Switcher (Valid=["vertical","horizontal"],Default="vertical");
     48        diameter as length (Brief="Vessel diameter");
     49
    4950        SET
    5051        Mw=PP.MolecularWeight();
    51        
     52
    5253        VARIABLES
    5354        in      Inlet as stream(Brief="Feed Stream");
     
    6364        vV as volume_mol (Brief="Vapour Molar volume");
    6465        Level as length (Brief="liquid height");
    65        
     66        Across as area (Brief="Flash Cross section area");
     67
    6668        EQUATIONS
    6769        "Component Molar Balance"
     
    98100       
    99101        "Geometry Constraint"
    100         V = ML* vL + MV*vV;
     102        V = ML * vL + MV * vV;
     103
     104        switch orientation
     105        case "vertical":
     106        "Cross Section Area"
     107                Across = 0.5 * asin(1) * diameter^2;
     108       
    101109        "Liquid Level"
    102         ML* vL = Across * Level;
     110                ML * vL = Across * Level;
     111
     112        case "horizontal":
     113        "Cylindrical Side Area"
     114                Across = 0.25*diameter^2 * (asin(1) - asin((diameter - 2*Level)/diameter)) +
     115                                (Level - 0.5*diameter)*sqrt(Level*(diameter - Level));
     116
     117        "Liquid Level"
     118                0.5 * asin(1) * diameter^2 * ML* vL = Across * V;
     119        end
    103120end
    104121
  • trunk/sample/stage_separators/sample_flash.mso

    r213 r235  
    6565       
    6666        SET
    67         fl.V = 2000 * 'm^3';
    68         fl.Across = 1 * 'm^2';
    69        
     67        fl.V = 50 * 'm^3';
     68        fl.diameter = 2 * 'm';
     69        fl.orientation = "vertical";
     70#       fl.orientation = "horizontal";
     71
    7072        INITIAL
    7173        fl.OutletL.T = 338 * 'K';
    72         fl.Level = 1 * 'm';
     74        fl.Level = 0.4 * 'm';
    7375
    7476        fl.OutletL.z(1) = 0.1;
     
    8385        Dynamic = true;
    8486        TimeStep = 0.1;
    85         TimeEnd = 20;
     87        TimeEnd = 5;
    8688        TimeUnit = 'h';
    8789end
Note: See TracChangeset for help on using the changeset viewer.