Changeset 109


Ignore:
Timestamp:
Jan 11, 2007, 8:03:27 PM (16 years ago)
Author:
Argimiro Resende Secchi
Message:

corrected the water specific density at 60F and uncommented some equations.

Location:
mso
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • mso/eml/pressure_changers/valve.mso

    r108 r109  
    2727*               * Isentalpic
    2828*               * Liquid
    29 *       
     29*
    3030*---------------------------------------------------------------------
    3131* Author: Estefane Horn, Núbia do Carmo Ferreira
     
    5454       
    5555        SET
    56         rho60F = 99.022 * "kg/m^3";
     56        rho60F = 999.2  * "kg/m^3";
    5757       
    5858        EQUATIONS
     
    6969        Outlet.z = Inlet.z;
    7070
    71         "Valve Equation"
    72         #if Pdiff >= 0 then
    73         #       "Flow"
     71        if Pdiff > 0 then
     72                "Valve Equation - Flow"
    7473                Qv = fc*cv*sqrt(Pdiff/Gf);     
    75         #else
    76         #       "Closed"
    77         #       Pdiff < 0 * "kmol/h";
    78         #end
     74        else
     75                "Valve Equation - Closed"
     76                Qv = 0 * "m^3/h";
     77        end
    7978       
    8079        "Calculate Gf"
     
    9291end
    9392
    94 #*----------------------------------------------------------------------
    95 * Model of a  valve_linear
    96 *---------------------------------------------------------------------*#
    97 #*      - Specify:
    98 *               * Opening fraction (x)
    99 *
    100 *---------------------------------------------------------------------*#
    101 
    102        
    10393Model valve_linear as valve_basic
    10494
     
    113103end
    114104
    115 #*----------------------------------------------------------------------
    116 * Model of a valve_parabolic
    117 *---------------------------------------------------------------------*#
    118 #*      - Specify:
    119 *               * Opening fraction (x)
    120 *
    121 *---------------------------------------------------------------------*#
    122 
    123105Model valve_parabolic as valve_basic
    124106
     
    136118end
    137119
    138 #*----------------------------------------------------------------------
    139 * Model of a  valve_equal
    140 *---------------------------------------------------------------------*#
    141 #*      - Specify:
    142 *               * Opening fraction (x)
    143 *
    144 *---------------------------------------------------------------------*#
    145 
    146         Model valve_equal as valve_basic
     120Model valve_equal as valve_basic
    147121
    148122        PARAMETERS
     
    158132
    159133end
    160 
    161 #*----------------------------------------------------------------------
    162 * Model of a  valve_quick
    163 *---------------------------------------------------------------------*#
    164 #*      - Specify:
    165 *               * Opening fraction (x)
    166 *
    167 *---------------------------------------------------------------------*#
    168134
    169135Model valve_quick as valve_basic
     
    245211                (1-Inlet.v)*PP.LiquidVolume((Inlet.T+Outlet.T)/2, (Inlet.P+Outlet.P)/2, Outlet.z);
    246212
    247         #if Inlet.P > Outlet.P then
    248         #       "Flow"
    249                 Outlet.F * v = cv*x*sqrt((Inlet.P - Outlet.P)*rho_ref / rho ) ;
    250         #else
    251         #       "Closed"
    252         #       Outlet.F = 0 * "kmol/h";
    253         #end
    254 end
     213        if Inlet.P > Outlet.P then
     214                "Flow"
     215                Outlet.F * v = k*x*sqrt((Inlet.P - Outlet.P)*rho_ref / rho ) ;
     216        else
     217                "Closed"
     218                Outlet.F = 0 * "kmol/h";
     219        end
     220end
  • mso/sample/pressure_changers/sample_valve.mso

    r81 r109  
    5959        OPTIONS
    6060        relativeAccuracy = 1e-6;
    61         time = [0:20:1000];
    6261        mode = "steady";
    63        
    6462end
    6563
     
    9795        OPTIONS
    9896        relativeAccuracy = 1e-6;
    99         time = [0:20:1000];
    10097        mode = "steady";
    101        
    10298end
    10399
     
    128124        V1.Qv = 113.56                  * "m^3/h";
    129125        V1.x = 0.5;
    130         V1.cv = 3.997           * "m^3/h/kPa^0.5";
     126        V1.cv = 10                              * "m^3/h/kPa^0.5";
     127       
     128        OPTIONS
     129        relativeAccuracy = 1e-6;
     130        mode = "steady";
    131131end
Note: See TracChangeset for help on using the changeset viewer.