Ignore:
Timestamp:
Aug 25, 2008, 2:11:12 AM (14 years ago)
Author:
gerson bicca
Message:

fixed some equations (compressor model)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/pressure_changers/compressor.mso

    r606 r608  
    3939VARIABLES
    4040
    41         PolyCoeff       as positive             (Brief = "Polytropic Coefficient", Lower=0.2);
    42         IseCoeff        as positive     (Brief = "Isentropic Coefficient", Lower=0.2);
    4341        Pratio          as positive             (Brief = "Pressure Ratio", Symbol ="P_{ratio}");       
    4442        Pdrop           as press_delta  (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
    4543        Pincrease       as press_delta  (Brief = "Pressure Increase", DisplayUnit = 'kPa', Symbol ="P_{incr}");
    4644
    47         Head                    as energy_mass  (Brief = "Actual Head");
     45        Head                    as energy_mass  (Brief = "Actual Head",Hidden=true);
    4846        HeadIsentropic  as energy_mass  (Brief = "Isentropic Head");
    4947        HeadPolytropic  as energy_mass  (Brief = "Polytropic Head");
    50         Tisentropic             as temperature  (Brief = "Isentropic Temperature");
    5148       
    52         PolytropicEff                           as efficiency   (Brief = "Polytropic efficiency");
    53         IsentropicEff                           as efficiency   (Brief = "Isentropic efficiency");
    54         EfficiencyOperation             as efficiency   (Brief = "Compressor efficiency - Polytropic or Isentropic (See Compressor Type)");
    55         MechanicalEff                           as efficiency   (Brief = "Mechanical efficiency");
     49        EfficiencyOperation     as positive     (Brief = "Compressor efficiency - Polytropic or Isentropic (See Compressor Type)",Lower=1E-3,Upper=1);
     50        MechanicalEff                   as positive     (Brief = "Mechanical efficiency",Lower=1E-3,Upper=1);
    5651       
    57         FluidPower      as power                        (Brief = "Fluid Power");
    58         BrakePower      as power                        (Brief = "Brake Power");
    59         PowerLoss       as power                        (Brief = "Power Losses");
    60         Mwm                     as molweight            (Brief = "Mixture Molar Weight");
    61         rho                     as dens_mass            (Brief = "Mass Density");
    62         Zfac_in         as fraction             (Brief = "Compressibility factor at inlet");
    63         Zfac_out        as fraction             (Brief = "Compressibility factor at outlet");
     52        FluidPower      as power        (Brief = "Fluid Power");
     53        BrakePower      as power        (Brief = "Brake Power");
     54        PowerLoss       as power        (Brief = "Power Losses");
     55       
     56        PolyCoeff                       as positive             (Brief = "Polytropic Coefficient", Lower=0.2,Protected=true);
     57        IseCoeff                        as positive     (Brief = "Isentropic Coefficient", Lower=0.2,Protected=true);
     58        PolytropicEff           as positive     (Brief = "Polytropic efficiency",Lower=1E-3,Upper=1,Protected=true);
     59        IsentropicEff           as positive     (Brief = "Isentropic efficiency",Lower=1E-3,Upper=1,Protected=true);
     60
     61        Tisentropic     as temperature          (Brief = "Isentropic Temperature",Protected=true);
     62        Mwm                     as molweight            (Brief = "Mixture Molar Weight",Protected=true);
     63        rho                     as dens_mass            (Brief = "Mass Density",Protected=true);
     64        Zfac_in         as fraction             (Brief = "Compressibility factor at inlet",Protected=true);
     65        Zfac_out        as fraction             (Brief = "Compressibility factor at outlet",Protected=true);
    6466
    6567in      Inlet   as stream       (Brief = "Inlet stream", PosX=0.437, PosY=1, Symbol="_{in}");
     
    112114        PP.VapourEntropy(Tisentropic, Outlet.P, Outlet.z) = PP.VapourEntropy(Inlet.T, Inlet.P, Inlet.z);
    113115
    114 "Fluid Power"
    115         FluidPower = Head*sum(Mw*Inlet.z)*Inlet.F;
    116 
    117116"Brake Power"
    118117        BrakePower = -WorkIn;
    119118
    120119"Brake Power"
    121         BrakePower = FluidPower/MechanicalEff;
     120        BrakePower*MechanicalEff = FluidPower;
    122121
    123122"Power Loss"
     
    139138        HeadPolytropic = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(PolyCoeff/(PolyCoeff-1.001))*Rgas*Inlet.T*((Outlet.P/Inlet.P)^((PolyCoeff-1.001)/PolyCoeff) - 1);
    140139
     140"Fluid Power"
     141        FluidPower = Head*sum(Mw*Inlet.z)*Inlet.F*EfficiencyOperation+PowerLoss;
     142
    141143switch CompressorType
    142144
    143145        case "Isentropic Operation":
    144146
    145 "Discharge Temperature"
     147"Efficiency"
    146148        EfficiencyOperation = IsentropicEff;
    147 
     149       
    148150        case "Polytropic Operation":
    149151
    150 "Discharge Temperature"
     152"Efficiency"
    151153        EfficiencyOperation = PolytropicEff;
    152154
Note: See TracChangeset for help on using the changeset viewer.