Changeset 530


Ignore:
Timestamp:
May 30, 2008, 5:21:43 PM (14 years ago)
Author:
Argimiro Resende Secchi
Message:

Fix some examples.

Location:
trunk
Files:
10 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/reactors/equil.mso

    r471 r530  
    112112
    113113        for i in [1:NComp]
    114           if (Outletm.z(i) > 0) then
     114          if (Outletm.z(i) > 1e-16) then
    115115            "Molar conversion"
    116116            Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
     
    203203
    204204        for i in [1:NComp]
    205           if (Outletm.z(i) > 0) then
     205          if (Outletm.z(i) > 1e-16) then
    206206            "Molar conversion"
    207207            Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
  • trunk/eml/reactors/gibbs.mso

    r471 r530  
    108108                G(i) + sumt(lambda*na(:,i)) = -Rg*Outlet.T*ln(activ(i));
    109109
    110           if (Outletm.z(i) > 0) then
     110          if (Outletm.z(i) > 1e-16) then
    111111            "Molar conversion"
    112112            Fi(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
     
    195195                G(i) + sumt(lambda*na(:,i)) = -Rg*Outlet.T*ln(activ(i));
    196196
    197           if (Outletm.z(i) > 0) then
     197          if (Outletm.z(i) > 1e-16) then
    198198            "Molar conversion"
    199199            Fi(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
  • trunk/eml/reactors/stoic.mso

    r426 r530  
    7373       
    7474        for i in [1:NComp]
    75           if (Outletm.z(i) > 0) then
     75          if (Outletm.z(i) > 1e-16) then
    7676            "Molar conversion"
    7777            Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
     
    124124
    125125        for i in [1:NComp]
    126           if (Outletm.z(i) > 0) then
     126          if (Outletm.z(i) > 1e-16) then
    127127            "Molar conversion"
    128128            Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
  • trunk/eml/reactors/yield.mso

    r426 r530  
    8585       
    8686        for i in [1:NComp]
    87           if (Outletm.z(i) > 0) then
     87          if (Outletm.z(i) > 1e-16) then
    8888            "Molar conversion"
    8989            Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
  • trunk/sample/controllers/sample_PIDIncr_MIMO.mso

    r342 r530  
    264264        col.trays.Emv = 1;
    265265#       col.sptop.Outlet2.F = 75 * 'kmol/h'; # reflux
     266    col.alfaTopo = 2;
    266267
    267268        xtop_sp = 1e-3; # benzene
     
    337338        TimeUnit = 'h';
    338339        InitialFile = "Column_MPID.rlt";
     340        DAESolver(File="dassl");
    339341        #GuessFile = "Column_MPID.rlt";
    340342        #Dynamic = false;
  • trunk/sample/controllers/sample_PID_MIMO.mso

    r340 r530  
    265265        col.trays.Emv = 1;
    266266#       col.sptop.Outlet2.F = 75 * 'kmol/h'; # reflux
     267    col.alfaTopo = 2;
    267268
    268269        xtop_sp = 1e-3; # benzene
     
    338339        TimeUnit = 'h';
    339340        InitialFile = "Column_ctrl.rlt";
     341        DAESolver(File="dassl");
    340342        #GuessFile = "Column_ctrl.rlt";
    341343        #Dynamic = false;
  • trunk/sample/costs/sample_distillation_kettle_cond_cost.mso

    r331 r530  
    8080        col.pump1.dP = 16 * 'kPa';
    8181        col.trays.Emv = 1;
     82        col.alfaTopo = 2;
    8283       
    8384        SET
     
    127128
    128129        # reboiler
    129         col.reb.OutletL.T = 330 *'K';
     130        col.reb.OutletL.T = 300 *'K';
    130131        col.reb.Level = 1 * 'm';
    131132        col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
    132133
    133134        # column trays
    134         col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
     135        col.trays.OutletL.T = [290:(300-290)/(col.NTrays-1):300] * 'K';
    135136        col.trays.Level = 1.2 * col.trays.hw;
    136137        col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
  • trunk/sample/miscellaneous/tenprobs/prob03.mso

    r283 r530  
    170170#               File = "ipopt_emso",
    171171                File = "complex",
    172                 ObjTol = 1e-6,
    173                 AcceptableObjTol = 1e-6);
     172                ObjTol = 1e-6);
    174173        Dynamic = false;
    175174end
     
    197196                File = "ipopt_emso",
    198197#               File = "complex",
    199                 ObjTol = 1e-6,
    200                 AcceptableObjTol = 1e-6);
     198                ObjTol = 1e-6);
    201199        Dynamic = false;
    202200end
     
    225223                File = "ipopt_emso",
    226224#               File = "complex",
    227                 ObjTol = 1e-6,
    228                 AcceptableObjTol = 1e-6);
    229         Dynamic = false;
    230 end
     225                ObjTol = 1e-6);
     226        Dynamic = false;
     227end
  • trunk/sample/stage_separators/sample_column.mso

    r515 r530  
    9494
    9595        INITIAL
    96         sec.trays.OutletL.T = 290 *'K';
     96        sec.trays.OutletL.T = 270 *'K';
    9797        sec.trays.Level = 0.9 * sec.trays.hw;
    9898        sec.trays.OutletL.z(1) = 0.5;
     
    178178
    179179        INITIAL
    180         sec.trays.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *'K';
     180        sec.trays.OutletL.T = [290:(300-290)/(sec.NTrays-1):300] *'K';
    181181        sec.trays.Level = 0.3 * sec.trays.hw;
    182182        sec.trays.OutletL.z(1) = 0.5;
     
    267267
    268268        # reboiler
    269         col.reb.OutletL.T = 330 *'K';
     269        col.reb.OutletL.T = 350 *'K';
    270270        col.reb.Level = 1 * 'm';
    271271        col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
     
    499499        TimeUnit = 'h';
    500500        InitialFile = "Column_ctrl.rlt";
     501        DAESolver(File="dassl");
    501502        #GuessFile = "Column_ctrl.rlt";
    502503        #Dynamic = false;
  • trunk/sample/stage_separators/sample_columnReact.mso

    r346 r530  
    125125        PIDLreb.Parameters.gamma=1;
    126126        PIDLreb.Parameters.beta=1;
    127         PIDLreb.Parameters.gain=1;
    128         PIDLreb.Parameters.intTime=10*'s';
    129         PIDLreb.Parameters.derivTime=1*'s';
     127        PIDLreb.Parameters.gain=0.9;
     128        PIDLreb.Parameters.intTime=20*'s';
     129        PIDLreb.Parameters.derivTime=0*'s';
    130130        PIDLreb.Options.action=-1;
    131131        PIDLreb.Options.clip=1;
     
    174174        PIDTcond.Parameters.gamma=1;
    175175        PIDTcond.Parameters.beta=1;
    176         PIDTcond.Parameters.gain=1;
    177         PIDTcond.Parameters.intTime=10*'s';
     176        PIDTcond.Parameters.gain=0.9;
     177        PIDTcond.Parameters.intTime=100*'s';
    178178        PIDTcond.Parameters.derivTime=1*'s';
    179179        PIDTcond.Options.action=1;
Note: See TracChangeset for help on using the changeset viewer.