Changeset 218


Ignore:
Timestamp:
Mar 15, 2007, 11:55:04 PM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Remove more convergence problems in samples.

Location:
branches/newlanguage
Files:
1 added
5 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/eml/stage_separators/flashPH.mso

    r210 r218  
    4848in      Q as heat_rate (Brief="Rate of heat supply");
    4949        vfrac as fraction(Brief="Real vaporization fraction");
    50         vsat as Real(Lower=-5, Upper=5, Brief="Vaporization fraction if saturated");
     50        vsat as Real(Lower=-0.1, Upper=1.1, Brief="Vaporization fraction if saturated");
    5151        Tsat as temperature(Lower=173, Upper=1473, Brief="Temperature if saturated");
    52         xsat(NComp) as Real(Lower=-5, Upper=5, Brief="Liquid composition if saturated");
    53         ysat(NComp) as Real(Lower=-5, Upper=5, Brief="Vapour composition if saturated");
     52        xsat(NComp) as Real(Lower=0, Upper=1, Brief="Liquid composition if saturated");
     53        ysat(NComp) as Real(Lower=0, Upper=1, Brief="Vapour composition if saturated");
    5454       
    5555        zero_one as fraction(Brief="Regularization Variable");
  • branches/newlanguage/sample/miscellaneous/sample_galvanic.mso

    r182 r218  
    4444        VARIABLES
    4545
    46         y1  as  current(Default=0.5, Brief="Mole fraction of NiOOH");
    47         y2  as  current(Default=0.5, Brief="Potential difference at the solid-liquid interface");
     46        y1  as  Real(Default=0.5, Brief="Mole fraction of NiOOH");
     47        y2  as  Real(Default=0.5, Brief="Potential difference at the solid-liquid interface");
    4848
    4949
     
    5353       
    5454        "Equation 1(a)"
    55         pho*V/W*diff(y1) = i01*(2*(1-y1)*exp(0.5*F/(R*T)*(y2-phi1))-2*y1*exp(-0.5*F/(R*T)*(y2-phi1)))/F;
     55        pho*V/W*diff(y1)*'s' = i01*(2*(1-y1)*exp(0.5*F/(R*T)*(y2-phi1))-2*y1*exp(-0.5*F/(R*T)*(y2-phi1)))/F;
    5656       
    5757        "Equation 1(b)"
  • branches/newlanguage/sample/miscellaneous/sample_streams.mso

    r182 r218  
    4747FlowSheet SampleStreamsTherm
    4848        PARAMETERS
    49         PP   as Plugin(Brief="Physical Properties", File="vrpp");
     49        PP   as Plugin(Brief="Physical Properties", Type="PP",
     50                Components = [ "isobutane", "benzene", "ethanol" ],
     51                LiquidModel = "IdealLiquid",
     52                VapourModel = "Ideal");
    5053        NComp as Integer;
    5154
    5255        SET
    53         PP.LiquidModel = "IdealLiquid";
    54         PP.VapourModel = "Ideal";
    55         PP.Components = [ "isobutane", "benzene", "ethanol" ];
    5656        NComp = PP.NumberOfComponents;
    5757
  • branches/newlanguage/sample/optimization/ammonia.mso

    r216 r218  
    108108FlowSheet Ammonia
    109109        PARAMETERS
    110         PP      as Plugin(Brief="Physical Properties", Type="PP",Components = ["hydrogen", "nitrogen", "argon", "methane", "ammonia"],
    111         LiquidModel = "APR",    VapourModel = "APR");
     110        PP      as Plugin(Brief="Physical Properties", Type="PP",
     111                Components = ["hydrogen", "nitrogen", "argon", "methane", "ammonia"],
     112                LiquidModel = "APR",
     113                VapourModel = "APR");
    112114        NComp   as Integer;
    113115        SET
     
    182184        RelativeAccuracy = 1e-5;
    183185end
    184 
  • branches/newlanguage/sample/stage_separators/sample_column.mso

    r213 r218  
    196196        PP      as Plugin(Brief="Physical Properties",
    197197                Type="PP",
    198                 Components = [ "isobutane", "benzene"],
     198                Components = [ "isobutane", "n-pentane", "propylene",
     199                "benzene", "isobutene" ],
    199200                LiquidModel = "PR",
    200201                VapourModel = "PR"
     
    390391        LCtop.Parameters.intTime = 10*'s';
    391392        LCtop.Parameters.gain = 1;
    392         LCtop.Parameters.derivTime = 1*'s';
     393        LCtop.Parameters.derivTime = 0*'s';
    393394        LCtop.Ports.setPoint = (1.0 * 'm' - Hmint)/(Hmaxt-Hmint);
    394395        LCtop.Ports.input = Had_top;
     
    408409        LCbot.Parameters.intTime = 100*'s';
    409410        LCbot.Parameters.gain = 1;
    410         LCbot.Parameters.derivTime = 1*'s';
     411        LCbot.Parameters.derivTime = 0*'s';
    411412        LCbot.Ports.setPoint = (1.0 * 'm' - Hminb)/(Hmaxb-Hminb);
    412413        LCbot.Ports.input = Had_bot;
     
    481482        col.reb.Level = 1 * 'm';
    482483        col.reb.OutletL.z([1:4]) = [0.1, 0.4, 0.1, 0.3];
    483        
     484
    484485        # column trays
    485486        col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
     
    488489
    489490        OPTIONS
    490         RelativeAccuracy = 1e-3;
    491         TimeStep = 0.01;
    492         TimeEnd = 2;
     491        TimeStep = 0.1;
     492        TimeEnd = 5;
    493493        TimeUnit = 'h';
    494         #InitialFile = "Column_ctrl.rlt";
     494        InitialFile = "Column_ctrl.rlt";
    495495        #GuessFile = "Column_ctrl.rlt";
    496496        #Dynamic = false;
Note: See TracChangeset for help on using the changeset viewer.