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/compressor.mso

    r196 r210  
    3333outer PP                as Plugin               (Brief = "External Physical Properties", Type="PP");
    3434outer NComp     as Integer              (Brief = "Number of chemical components", Lower = 1);
    35         Effs            as positive     (Default = 0.72, Brief = "Isentropic efficiency", Lower = 0, Upper = 1);
    3635        R                       as positive     (Default = 8.31451, Brief = "Constant of Gases", Unit= 'kJ/kmol/K');
    3736        Mw(NComp)       as molweight    (Brief = "Molar Weight");
    38        
     37
    3938        VARIABLES
    40         n                       as positive             (Brief = "Politropic Coefficient");
    41         k                       as positive     (Brief = "Isentropic Coefficient");
     39        n                       as positive             (Brief = "Politropic Coefficient", Lower=0);
     40        k                       as positive     (Brief = "Isentropic Coefficient", Lower=0);
    4241        Cp              as cp_mol               (Brief = "Heat Capacity");
    4342        Cv                      as cv_mol               (Brief = "Heat Capacity");
    … …  
    4746        Ws                      as energy_mol   (Brief = "Isentropic Head");
    4847        Tiso            as temperature  (Brief = "Isentropic Temperature");
    49         Effp            as positive     (Brief = "Politropic efficiency", Lower = 0, Upper = 1);
     48        Effp            as efficiency   (Brief = "Politropic efficiency");
     49        Effs            as efficiency   (Brief = "Isentropic efficiency");
    5050        FPower          as power                (Brief = "Fluid Power");
    5151        Mwm                     as molweight    (Brief = "Mixture Molar Weight");
    52        
     52
    5353        SET
    5454        Mw = PP.MolecularWeight();
    … …  
    5757       
    5858        "Calculate Mwm for Inlet Mixture"
    59         Mwm = sum(Mw([1:NComp])*Inlet.z([1:NComp]));
     59        Mwm = sum(Mw*Inlet.z);
    6060
    6161        "Calculate Outlet Stream Pressure"
    … …  
    7272       
    7373        "Calculate Isentropic Coeficient"
    74         k = Cp/Cv;
     74        k * Cv = Cp;
    7575       
    7676        "Calculate Isentropic Head"
    … …  
    8787       
    8888        "Calculate Politropic Efficiency"
    89         Effp = (n/(n-1))/(k/(k-1));
     89        Effp * (n-1) * k = n * (k-1);
    9090       
    9191        "Calculate Politropic Head"
    … …  
    9999       
    100100        Outlet.z = Inlet.z;
    101 
    102 
    103101end
Note: See TracChangeset for help on using the changeset viewer.