Changeset 859

Timestamp:

Oct 9, 2009, 5:15:53 PM (13 years ago)

Author:

gerson bicca

Message:

updated samples

Location:

trunk/sample

Files:

• processes/Sample_Process.mso (modified) (1 diff)
• reactors/sample_pfr.mso (modified) (1 diff)

trunk/sample/processes/Sample_Process.mso

@@ -126 +126 @@
         
         Reac.str(z).T = Reac.Inlet.T;
-        Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
+        Reac.str(z).z(1:NComp) = Reac.Inlet.z(1:NComp);
 
 end

trunk/sample/reactors/sample_pfr.mso

@@ -39 +39 @@
 PARAMETERS
 
-        PP              as Plugin (Brief="Physical Properties",
+
+
+        PP as Plugin (Brief="Physical Properties",
+
                 Type = "PP",
-                Components = ["acetone", "acetic anhydride", "methane" ],
+
+                Components = ["acetone", "ketene", "methane" ],
+
                 LiquidModel = "PR",
+
                 VapourModel = "PR"
+
         );
+
         NComp   as Integer;
+
+SET
+
+        NComp = PP.NumberOfComponents;
+
+
 
 DEVICES
 
         s1      as source;
-        Reac as pfr;
+
+        Reac    as pfr;
+
+
+
+CONNECTIONS
+
+        s1.Outlet  to  Reac.Inlet;
+
+        
 
 SET
-        NComp = PP.NumberOfComponents;
-        
-        Reac.NDisc      = 10;
-        Reac.Across     = 0.7 * 'in^2';
-        Reac.L                  = 2.28 * 'm';
+
+        Reac.NDisc      = 15;
+
+        Reac.L          = 2.28 * 'm';
+
+        # Area was chosen to match the original example data
+
+        Reac.Across = 1.27*'m^3'/Reac.L;
+
+
+
+        # Only one reaction
+
         Reac.NReac      = 1;
 
-        # Reaction 1: A -> B + C
+        # Reaction 1: acetone -> ketene + methane
+
         Reac.stoic(:,1) = [-1, 1, 1];
+
+
+
+SPECIFY
+
+        # Feed specification
+
+        s1.Fw = 8000 * 'kg/h';
+
+        s1.T = 1035 * 'K';
+
+        s1.P = 1.6  * 'atm';
+
+        s1.Composition = [0.98, 0.01, 0.01];
+
+
+
+        # Adiabatic
+
+        Reac.q     = 0     * 'J/s';
+
         
+
 EQUATIONS
 
-for z in [1:Reac.NDisc]
+
+
+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';
 
+
+
 "Pressure Drop (no pressure drop)"
+
         Reac.str(z+1).P = Reac.str(z).P;
+
+
 
 end
 
-SPECIFY
-        s1.F = 138/1000 * 'kmol/h';
-        s1.T = 1035 * 'K';
-        s1.P = 1.6  * 'atm';
-        s1.Composition = [1, 0, 0];
 
-"Adiabatic"
-        Reac.q     = 0     * 'J/s';
-        
-CONNECTIONS
 
-        s1.Outlet  to  Reac.Inlet;
-        
 INITIAL
 
-for z in [2:Reac.NDisc+1]
-        
+
+
+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);
+
+        Reac.str(z).z(1:NComp) = Reac.Inlet.z(1:NComp);
 
 end
+
         
+
 OPTIONS
-        # This model can be solved in both steady or dynamic
+
         Dynamic         = true;
+
         TimeStep        = 0.05;
-        TimeEnd         = 3;
-        #Dynamic = false;
-        #DAESolver = "dassl";
+
+        TimeEnd         = 10;
+
+        GuessFile = "PFR_AceticAnhydride";
+
+        Dynamic = false; # Requires a GuessFile from a dynamic simulation
 end