Changeset 170 for branches/newlanguage/sample/reactors

Timestamp:

Mar 1, 2007, 3:35:36 PM (17 years ago)

Author:

gerson bicca

Message:

updated pfr reactor model for the new language

File:

• branches/newlanguage/sample/reactors/sample_pfr.mso (modified) (1 diff)

branches/newlanguage/sample/reactors/sample_pfr.mso

@@ -36 +36 @@
 
 FlowSheet PFR_AceticAnhydride
-        PARAMETERS
-        PP   as CalcObject(Brief="Physical Properties",File="vrpp");
-        NComp as Integer;
+        
+PARAMETERS
 
-        DEVICES
-        s1   as streamTP;
+        PP              as Plugin (Brief="Physical Properties",File="vrpp");
+        NComp   as Integer;
+
+DEVICES
+
+        s1      as source;
         Reac as pfr;
 
-        SET
-        PP.Components = ["acetone", "acetic anhydride", "methane" ];
-        PP.LiquidModel = "PR";
-        PP.VapourModel = "PR";
+SET
+
+        PP.Components   = ["acetone", "acetic anhydride", "methane" ];
+        PP.LiquidModel  = "PR";
+        PP.VapourModel  = "PR";
         
         NComp = PP.NumberOfComponents;
         
-        Reac.NDisc = 10;
-        Reac.Across = 0.7 * "in^2";
-        Reac.L = 2.28 * "m";
-        Reac.NReac = 1;
+        Reac.NDisc      = 10;
+        Reac.Across     = 0.7 * 'in^2';
+        Reac.L                  = 2.28 * 'm';
+        Reac.NReac      = 1;
 
         # Reaction 1: A -> B + C
         Reac.stoic(:,1) = [-1, 1, 1];
         
-        EQUATIONS
+EQUATIONS
 
-        for z in [1:Reac.NDisc]
-                "Reaction Rate = k*C(1)"
-                Reac.r(1,z)  = exp(34.34 - (34222 * "K") / Reac.str(z).T)*"1/s" * Reac.C(1,z);
+for z in [1:Reac.NDisc]
+        
+"Reaction Rate = k*C(1)"
+        Reac.r(1,z)  = exp(34.34 - (34222 * 'K') / Reac.str(z).T)*'1/s' * Reac.C(1,z);
 
-                "Heat of reaction"
-                Reac.Hr(1,z) = -80.77 * "kJ/mol";
+"Heat of reaction"
+        Reac.Hr(1,z) = -80.77 * 'kJ/mol';
 
-                "Pressure Drop (no pressure drop)"
-                Reac.str(z+1).P = Reac.str(z).P;
-        end
+"Pressure Drop (no pressure drop)"
+        Reac.str(z+1).P = Reac.str(z).P;
 
-        SPECIFY
-        s1.F = 138/1000 * "kmol/h";
-        s1.T = 1035 * "K";
-        s1.P = 1.6  * "atm";
-        s1.z = [1, 0, 0];
+end
 
-        "Adiabatic"
-        Reac.q     = 0     * "J/s";
+SPECIFY
+        s1.Outlet.F = 138/1000 * 'kmol/h';
+        s1.Outlet.T = 1035 * 'K';
+        s1.Outlet.P = 1.6  * 'atm';
+        s1.Outlet.z = [1, 0, 0];
+
+"Adiabatic"
+        Reac.q     = 0     * 'J/s';
         
-        CONNECTIONS
-        s1  to  Reac.Inlet;
+CONNECTIONS
+
+        s1.Outlet  to  Reac.Inlet;
         
-        INITIAL
-        for z in [2:Reac.NDisc+1]
-                Reac.str(z).T = Reac.Inlet.T;
-                Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
-        end
+INITIAL
+
+for z in [2:Reac.NDisc+1]
         
-        OPTIONS
+        Reac.str(z).T = Reac.Inlet.T;
+        Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
+
+end
+        
+OPTIONS
         # This model can be solved in both steady or dynamic
-        #mode = "steady";
-        time = [0:0.05:3] * "s";
-        NLASolver = "sundials";
+        Dynamic         = true;
+        TimeStep        = 0.05;
+        TimeEnd         = 3;
+        NLASolver       = "sundials";
         #DAESolver = "dassl";
 end