Changeset 530

Ignore:
Timestamp:
May 30, 2008, 5:21:43 PM (15 years ago)
Message:

Fix some examples.

Location:
trunk
Files:
10 edited

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

 r471 for i in [1:NComp] if (Outletm.z(i) > 0) then if (Outletm.z(i) > 1e-16) then "Molar conversion" Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i)); for i in [1:NComp] if (Outletm.z(i) > 0) then if (Outletm.z(i) > 1e-16) then "Molar conversion" Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
• trunk/eml/reactors/gibbs.mso

 r471 G(i) + sumt(lambda*na(:,i)) = -Rg*Outlet.T*ln(activ(i)); if (Outletm.z(i) > 0) then if (Outletm.z(i) > 1e-16) then "Molar conversion" Fi(i) = Outletm.F*Outletm.z(i)*(1 - conv(i)); G(i) + sumt(lambda*na(:,i)) = -Rg*Outlet.T*ln(activ(i)); if (Outletm.z(i) > 0) then if (Outletm.z(i) > 1e-16) then "Molar conversion" Fi(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
• trunk/eml/reactors/stoic.mso

 r426 for i in [1:NComp] if (Outletm.z(i) > 0) then if (Outletm.z(i) > 1e-16) then "Molar conversion" Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i)); for i in [1:NComp] if (Outletm.z(i) > 0) then if (Outletm.z(i) > 1e-16) then "Molar conversion" Outlet.F*Outlet.z(i) = Outletm.F*Outletm.z(i)*(1 - conv(i));
• trunk/eml/reactors/yield.mso

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

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

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

 r331 col.pump1.dP = 16 * 'kPa'; col.trays.Emv = 1; col.alfaTopo = 2; SET # reboiler col.reb.OutletL.T = 330 *'K'; col.reb.OutletL.T = 300 *'K'; col.reb.Level = 1 * 'm'; col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01]; # column trays col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K'; col.trays.OutletL.T = [290:(300-290)/(col.NTrays-1):300] * 'K'; col.trays.Level = 1.2 * col.trays.hw; col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
• trunk/sample/miscellaneous/tenprobs/prob03.mso

 r283 #               File = "ipopt_emso", File = "complex", ObjTol = 1e-6, AcceptableObjTol = 1e-6); ObjTol = 1e-6); Dynamic = false; end File = "ipopt_emso", #               File = "complex", ObjTol = 1e-6, AcceptableObjTol = 1e-6); ObjTol = 1e-6); Dynamic = false; end File = "ipopt_emso", #               File = "complex", ObjTol = 1e-6, AcceptableObjTol = 1e-6); Dynamic = false; end ObjTol = 1e-6); Dynamic = false; end
• trunk/sample/stage_separators/sample_column.mso

 r515 INITIAL sec.trays.OutletL.T = 290 *'K'; sec.trays.OutletL.T = 270 *'K'; sec.trays.Level = 0.9 * sec.trays.hw; sec.trays.OutletL.z(1) = 0.5; INITIAL sec.trays.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *'K'; sec.trays.OutletL.T = [290:(300-290)/(sec.NTrays-1):300] *'K'; sec.trays.Level = 0.3 * sec.trays.hw; sec.trays.OutletL.z(1) = 0.5; # reboiler col.reb.OutletL.T = 330 *'K'; col.reb.OutletL.T = 350 *'K'; col.reb.Level = 1 * 'm'; col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01]; TimeUnit = 'h'; InitialFile = "Column_ctrl.rlt"; DAESolver(File="dassl"); #GuessFile = "Column_ctrl.rlt"; #Dynamic = false;
• trunk/sample/stage_separators/sample_columnReact.mso

 r346 PIDLreb.Parameters.gamma=1; PIDLreb.Parameters.beta=1; PIDLreb.Parameters.gain=1; PIDLreb.Parameters.intTime=10*'s'; PIDLreb.Parameters.derivTime=1*'s'; PIDLreb.Parameters.gain=0.9; PIDLreb.Parameters.intTime=20*'s'; PIDLreb.Parameters.derivTime=0*'s'; PIDLreb.Options.action=-1; PIDLreb.Options.clip=1; PIDTcond.Parameters.gamma=1; PIDTcond.Parameters.beta=1; PIDTcond.Parameters.gain=1; PIDTcond.Parameters.intTime=10*'s'; PIDTcond.Parameters.gain=0.9; PIDTcond.Parameters.intTime=100*'s'; PIDTcond.Parameters.derivTime=1*'s'; PIDTcond.Options.action=1;
Note: See TracChangeset for help on using the changeset viewer.