Changeset 207 for trunk

Timestamp:

Mar 15, 2007, 2:17:35 AM (16 years ago)

Author:

Argimiro Resende Secchi

Message:

Fix initial conditions of sample_cstr.

Location:

trunk

Files:

• eml/reactors/cstr.mso (modified) (2 diffs)
• sample/reactors/sample_cstr.mso (modified) (4 diffs)

trunk/eml/reactors/cstr.mso

@@ -41 +41 @@
 
         EQUATIONS
+        "Molar concentration"
+        Outlet.z * M = C * Vr;
+
         "Molar fraction"
-        Outlet.z * M = C * Vr;
+        sum(Outlet.z) = 1;
+
         "Vapourisation Fraction"
         Outlet.v = Inlet.v;
@@ -65 +69 @@
         EQUATIONS
         "Component Molar Balance"
-        diff(Outlet.z*M) = (Inlet.F*Inlet.z - Outlet.F*Outlet.z) + sumt(stoic*r)*Vr;
+        diff(Outlet.z*M) = (Inlet.F*Inlet.z - Outlet.F*Outlet.z) + sumt(stoic*r)*Vr; #FIXME sumt()
 
         "Reactor Energy Balance"
         diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +
-                sum(Hr*sum(stoic*r))*Vr - q;
+                sum(Hr*sum(stoic*r))*Vr - q; # FIXME sum()
 end

trunk/sample/reactors/sample_cstr.mso

@@ -38 +38 @@
         SET
         PP.Components = ["1,2-propylene oxide", "methanol",
-                "water", "sulfuric acid", "1,2-propylene glycol" ];
+                "water", "1,2-propylene glycol" ];
         PP.LiquidModel = "PR";
         PP.VapourModel = "PR";
@@ -44 +44 @@
         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
@@ -57 +55 @@
         s1.P = 1     * "atm";
         s1.T = 534.67* "degR";
-        s1.z = [0.0468, 0.0783, 0.874, 0.0009, 0];
+        s1.z = [0.0468, 0.0783, 0.874, 0];
         s1.v = 0.0;
 
@@ -70 +68 @@
         
         INITIAL
-        Reac.Outlet.T = Reac.Inlet.T;
-        Reac.Outlet.z = Reac.Inlet.z;
+        Reac.Outlet.T = 500*"degR";
+        Reac.M = 100 * "kmol";
+        Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874];
 
         OPTIONS