Changeset 813 for branches/gui/sample

Timestamp:

Aug 5, 2009, 4:00:16 PM (14 years ago)

Author:

Rafael de Pelegrini Soares

Message:

Fixed PFR model and sample

File:

• branches/gui/sample/reactors/sample_pfr.mso (modified) (2 diffs)

branches/gui/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;
 
-        s1      as simple_source;
-        Reac as pfr;
+CONNECTIONS
+        s1.Outlet  to  Reac.Inlet;
+        
+SET
+        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;
 
-SET
-        NComp = PP.NumberOfComponents;
-        
-        s1.ValidPhases = "Vapour-Liquid";
-        #s1.CompositionBasis = "Molar";
-        
-        Reac.NDisc      = 10;
-        Reac.Dpipe      = 3 * 'in';
-        Reac.L                  = 15 * 'm';
+        # Only one reaction
         Reac.NReac      = 1;
-        Reac.Roughness  = 4.572E-5*'m';
+        # Reaction 1: acetone -> ketene + methane
+        Reac.stoic(:,1) = [-1, 1, 1];
 
-        # Reaction 1: A -> B + C
-        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
@@ -77 +87 @@
         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.MolarComposition = [1, 0, 0];
-
-        Reac.q     = 0     * 'J/s';
-        
-CONNECTIONS
-
-        s1.Outlet  to  Reac.Inlet;
-        
 INITIAL
 
 for z in [2:Reac.NDisc+1] do
+        Reac.str(z).T = Reac.Inlet.T;
+        Reac.str(z).z(1:NComp) = Reac.Inlet.z(1:NComp);
+end
         
-        Reac.str(z).T = Reac.Inlet.T;
-        Reac.vel(z) = Reac.vel(1);
-        #Reac.str(z).F = Reac.str(1).F;
-        Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
-
-end
-
-VARIABLES
-
-SOMA as Real;
-
-EQUATIONS
-
-SOMA = sum(Reac.Outlet.z);
-
 OPTIONS
-        # This model can be solved in both steady or dynamic
         Dynamic         = true;
         TimeStep        = 0.05;
-        TimeEnd         = 3;
-        DAESolver(File="sundials", RelativeAccuracy=1e-2);
-        NLASolver(RelativeAccuracy=1e-4);
+        TimeEnd         = 10;
         #InitialFile = "PFR_AceticAnhydride";
-        #Dynamic = false;
-        #DAESolver = "dassl";
+        #Dynamic = false; # Requires a GuessFile from a dynamic simulation
 end