Changeset 208

Timestamp:

Mar 15, 2007, 2:18:39 AM (17 years ago)

Author:

Argimiro Resende Secchi

Message:

Fix more sample problems to new language.

Location:

branches/newlanguage

Files:

• eml/reactors/cstr.mso (modified) (3 diffs)
• eml/stage_separators/condenser.mso (modified) (1 diff)
• eml/stage_separators/flashPH.mso (modified) (1 diff)
• sample/heat_exchangers/Sample_Heater.mso (modified) (1 diff)
• sample/mixers_splitters/sample_mixer.mso (modified) (1 diff)
• sample/reactors/sample_cstr.mso (modified) (3 diffs)
• sample/stage_separators/sample_column.mso (modified) (4 diffs)
• sample/stage_separators/sample_columnReact.mso (modified) (1 diff)
• sample/stage_separators/sample_condenser.mso (modified) (2 diffs)
• sample/stage_separators/sample_flashPH.mso (modified) (1 diff)
• sample/stage_separators/sample_reboiler.mso (modified) (2 diffs)
• sample/stage_separators/sample_tray.mso (modified) (1 diff)

branches/newlanguage/eml/reactors/cstr.mso

@@ -48 +48 @@
 EQUATIONS
 
+"Molar Concentration"
+        Outlet.z * M = C * Vr;
+
 "Molar fraction"
-        Outlet.z * M = C * Vr;
-        
-"Vapourisation Fraction"
-        Outlet.v = Inlet.v;
-
+    sum(Outlet.z) = 1;
 end
 
@@ -85 +84 @@
         Info            =
         "write some information";
-        
+
 VARIABLES
         
@@ -97 +96 @@
 
 "Reactor Energy Balance"
-        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +sum(Hr*sum(stoic*r))*Vr - q;
+        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +sum(Hr*sum(stoic*r))*Vr - q; #FIXME sum(sum())
 
 end

branches/newlanguage/eml/stage_separators/condenser.mso

@@ -128 +128 @@
         "Pressure"
         DP = InletV.P - OutletL.P;
-
-        "Vapourisation Fraction"
-        OutletL.v = 0.0;
 end
 

branches/newlanguage/eml/stage_separators/flashPH.mso

@@ -49 +49 @@
         vfrac as fraction(Brief="Real vaporization fraction");
         vsat as Real(Lower=-5, Upper=5, Brief="Vaporization fraction if saturated");
-        Tsat as temperature(Lower=173, Upper=473, Brief="Temperature if saturated");
+        Tsat as temperature(Lower=173, Upper=1473, Brief="Temperature if saturated");
         xsat(NComp) as Real(Lower=-5, Upper=5, Brief="Liquid composition if saturated");
         ysat(NComp) as Real(Lower=-5, Upper=5, Brief="Vapour composition if saturated");

branches/newlanguage/sample/heat_exchangers/Sample_Heater.mso

@@ -58 +58 @@
         streamcold1.Outlet.F            = 100   * 'kmol/h';
         streamcold1.Outlet.P            = 1     * 'atm';
-        streamcold1.Outlet.T            = 100   * 'K';
+        streamcold1.Outlet.T            = 300   * 'K';
     streamcold1.Outlet.z                = [1,0];
 
         streamcold2.Outlet.F            = 20*'kmol/h';
         streamcold2.Outlet.P            = 1     * 'atm';
-        streamcold2.Outlet.T            = 100   * 'K';
+        streamcold2.Outlet.T            = 300   * 'K';
     streamcold2.Outlet.z                = [0,1];
         
-        Heater.Outlet.T         = 70*'K';
+        Heater.Outlet.T         = 330*'K';
         Heater.Outlet.P   = 1 * 'atm';
 

branches/newlanguage/sample/mixers_splitters/sample_mixer.mso

@@ -46 +46 @@
         
         CONNECTIONS
-        stream1 to mixer1.Inlet_mixer(1);
-        stream2 to mixer1.Inlet_mixer(2);
+        stream1.Outlet to mixer1.Inlet_mixer(1);
+        stream2.Outlet to mixer1.Inlet_mixer(2);
         
         SPECIFY

branches/newlanguage/sample/reactors/sample_cstr.mso

@@ -39 +39 @@
 SET
 
-        PP.Components   = ["1,2-propylene oxide", "methanol","water","sulfuric acid", "1,2-propylene glycol" ];
+        PP.Components   = ["1,2-propylene oxide", "methanol","water","1,2-propylene glycol" ];
         PP.LiquidModel  = "PR";
         PP.VapourModel  = "PR";
@@ -45 +45 @@
         NComp                   = PP.NumberOfComponents;
         Reac.NReac      = 1;
-        Reac.stoic              = [-1, 0, -1, 0, 1];
+        Reac.stoic              = [-1, 0, -1, 1];
         
 EQUATIONS
 
-#       Reac.Hr = -36400 * "Btu/lbmol" - 7 * "Btu/lbmol/degR" * (Reac.Fout.T - 528 * "degR");
-#       Reac.r = 16.96e12 / "h" * exp(- 32400 * "Btu/lbmol"/R/Reac.Fout.T)
-#               * Reac.C(1);
-        Reac.Hr = 0 * 'Btu/lbmol';
-        Reac.r  = 0 * 'mol/l/s';
+        Reac.Hr = -36400 * 'Btu/lbmol' - 7 * 'Btu/lbmol/degR' * (Reac.Outlet.T - 528 * 'degR');
+        Reac.r = 16.96e12 / 'h' * exp(- 32400 * 'Btu/lbmol'/R/Reac.Outlet.T)
+                * Reac.C(1);
 
-        SPECIFY
+SPECIFY
         s1.Outlet.F = 52126 * 'mol/h';
         s1.Outlet.P = 1     * 'atm';
         s1.Outlet.T = 534.67* 'degR';
-        s1.Outlet.z = [0.0468, 0.0783, 0.874, 0.0009, 0];
-        s1.Outlet.v = 0.0;
+        s1.Outlet.z = [0.0468, 0.0783, 0.8749, 0];
 
         Reac.Vr                 = 300   * 'gal';
         Reac.Outlet.F   = Reac.Inlet.F;
 
-"Adiabatic"
+        "Adiabatic"
         Reac.q                  = 0     * 'J/s';
         Reac.Outlet.P   = 1     * 'atm';
@@ -75 +72 @@
 INITIAL
 
-        Reac.Outlet.T = Reac.Inlet.T;
-        Reac.Outlet.z = Reac.Inlet.z;
+        Reac.Outlet.T = 530 * 'degR';
+        Reac.M = 100 * 'kmol';
+        Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874];
 
 OPTIONS

branches/newlanguage/sample/stage_separators/sample_column.mso

@@ -105 +105 @@
 # can be obtained by using the results from that transient:
 
-        #guessFile="SectionColumn_Test_with2tray.rlt";
-        #mode = "steady";
+        #GuessFile="SectionColumn_Test_with2tray.rlt";
+        #Dynamic = false;
 end
 
@@ -184 +184 @@
         TimeStep = 1;
         TimeEnd = 100;
-        #guessFile="SectionColumn_Test_with8tray.rlt";
-        #mode = "steady";
+        #GuessFile="SectionColumn_Test_with8tray.rlt";
+        #Dynamic = false;
 end
 
@@ -277 +277 @@
         TimeEnd = 50;
         #time = [0:0.01:1, 2:50];
-        #guessFile="Distillation_kettle_cond_Test.rlt";
-        #mode = "steady";       
+        #GuessFile="Distillation_kettle_cond_Test.rlt";
+        #Dynamic = false;       
 end
 
@@ -488 +488 @@
         TimeEnd = 2;
         TimeUnit = 'h';
-        #time = [0:0.01:0.1, 0.11:0.01:2]*"h";
-        #initialFile = "Column_ctrl.rlt";
-        #guessFile = "Column_ctrl.rlt";
-        #mode = "steady";
+        #InitialFile = "Column_ctrl.rlt";
+        #GuessFile = "Column_ctrl.rlt";
+        #Dynamic = false;
 end
 

branches/newlanguage/sample/stage_separators/sample_columnReact.mso

@@ -62 +62 @@
 
         CONNECTIONS
-        feed to col.trays(5).Inlet;
+        feed.Outlet to col.trays(5).Inlet;
         zero to col.reb.Inlet;
         zero to col.trays([1:4]).Inlet;

branches/newlanguage/sample/stage_separators/sample_condenser.mso

@@ -47 +47 @@
         
         CONNECTIONS
-        s1 to c1.InletV;
+        s1.Outlet to c1.InletV;
         c1.OutletL to sp.Inlet;
         Q to c1.Q;
@@ -96 +96 @@
         
         CONNECTIONS
-        s1 to c1.InletV;
+        s1.Outlet to c1.InletV;
         Q to c1.Q;
         

branches/newlanguage/sample/stage_separators/sample_flashPH.mso

@@ -44 +44 @@
         
         CONNECTIONS
-        s1 to fl.Inlet;
+        s1.Outlet to fl.Inlet;
         Q to fl.Q;
         

branches/newlanguage/sample/stage_separators/sample_reboiler.mso

@@ -46 +46 @@
         CONNECTIONS
         feed to r1.Inlet;
-        s1 to r1.InletL;
+        s1.Outlet to r1.InletL;
         Q to r1.Q;
         
@@ -102 +102 @@
 
         CONNECTIONS
-        s1 to r1.InletL;
+        s1.Outlet to r1.InletL;
         Q to r1.Q;
         

branches/newlanguage/sample/stage_separators/sample_tray.mso

@@ -43 +43 @@
         
         CONNECTIONS
-        feed to t1.Inlet;
+        feed.Outlet to t1.Inlet;
         inL to t1.InletL;
         inV to t1.InletV;