Ignore:
Timestamp:
Mar 15, 2007, 9:52:28 AM (16 years ago)
Author:
Argimiro Resende Secchi
Message:

Remove convergence problems of some samples.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/eml/pressure_changers/turbine.mso

    r196 r210  
    3737outer PP                as Plugin                       (Brief = "External Physical Properties", Type="PP");
    3838        Mw(NComp)       as molweight            (Brief = "Molar Weight");
    39         Eff     as positive                     (Default = 0.72, Brief = "Pump efficiency");
    40         Meff    as positive                     (Default = 1.0, Brief = "Brake efficiency");
    41         Beta    as positive                     (Default = 0, Brief = "Volumetric expansivity", Unit = '1/K');
    42        
    4339       
    4440        VARIABLES
     41        Eff     as efficiency           (Brief = "Pump efficiency");
     42        Meff    as efficiency           (Brief = "Brake efficiency");
     43        Beta    as positive             (Brief = "Volumetric expansivity", Unit = '1/K');
    4544        Head    as head                         (Brief = "Head Developed");
    4645        FPower  as power                        (Brief = "Fluid Power");
     
    5958        #Mixtures Properties
    6059        "Calculate Mwm for Inlet Mixture"
    61         Mwm = sum(Mw([1:NComp])*Inlet.z([1:NComp]));
    62        
     60        Mwm = sum(Mw*Inlet.z);
     61
    6362        "Calculate rho using a External Physical Properties Routine"
    6463        rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
     
    7877       
    7978        "Calculate Fluid Power"
    80         FPower = Pdiff * Inlet.F / (rho/Mwm);
     79        FPower * rho = Pdiff * Inlet.F * Mwm;
    8180       
    8281        "Calculate Brake Power"
     
    8786       
    8887        "Calculate Outlet Temperature"
    89         (Outlet.T - Inlet.T)*Cp = (Outlet.h - Inlet.h) -  Pdiff /(rho/Mwm) * (1-Beta*Inlet.T);
     88        (Outlet.T - Inlet.T) * rho * Cp = (Outlet.h - Inlet.h) * rho
     89        -  Pdiff * Mwm * (1-Beta*Inlet.T);
    9090       
    9191        "Calculate Outlet Enthalpy"
    92         Outlet.h - Inlet.h =  Pdiff / (rho/Mwm);
     92        (Outlet.h - Inlet.h) * rho =  Pdiff * Mwm;
    9393       
    9494        "Molar Balance"
     
    9999
    100100        "Calculate Head"
    101         Head = Pdiff/rho;
    102 
     101        Head * rho = Pdiff;
    103102end
Note: See TracChangeset for help on using the changeset viewer.