Changeset 188 for branches/newlanguage/sample/reactors/fogler/chap9/cstr_startup.mso

Timestamp:

Mar 7, 2007, 1:53:12 PM (17 years ago)

Author:

Rodolfo Rodrigues

Message:

Updated reactor samples for the new language

File:

• branches/newlanguage/sample/reactors/fogler/chap9/cstr_startup.mso (modified) (7 diffs)

branches/newlanguage/sample/reactors/fogler/chap9/cstr_startup.mso

@@ -57 +57 @@
         NComp           as Integer              (Brief="Number of components");
         stoic(NComp)as Real             (Brief="Stoichiometric number");
-        UA              as Real             (Brief="Exchange heat", Unit="Btu/h/degR");
+        UA              as Real             (Brief="Exchange heat", Unit='Btu/h/degR');
         V               as volume           (Brief="Volume of the reactor");
         Ta1             as temperature  (Brief="Cooling temperature");
@@ -66 +66 @@
         ko                      as frequency    (Brief="Frequency factor");
         E                       as energy_mol   (Brief="Activation energy");
-        R                       as Real                 (Brief="Universal gas constant", Unit="Btu/lbmol/degR", Default=1.987);
+        R                       as Real                 (Brief="Universal gas constant", Unit='Btu/lbmol/degR', Default=1.987);
         
         VARIABLES
-        C(NComp)        as conc_mol     (Brief="Molar concentration", Unit="lbmol/ft^3", Lower=0);
-        N(NComp)        as mol                  (Brief="Molar holdup", Unit="lbmol");
-        Co(NComp)       as conc_mol     (Brief="Initial molar concentration", Unit="lbmol/ft^3", Lower=0);
-        Fo(NComp)       as flow_mol     (Brief="Initial molar flow", Unit="lbmol/h");
-        T               as temperature  (Brief="Reactor temperature", Unit="degR");
-        To              as temperature  (Brief="Initial reactor temperature", Unit="degR");
-        Ta2             as temperature  (Brief="Temperature of heat exchange", Unit="degR");
-        r(NComp)        as reaction_mol (Brief="Rate of reaction", Unit="lbmol/ft^3/h");
-        k               as frequency    (Brief="Specific rate of reaction", Unit="1/h", Upper=1e5);
-        mc              as flow_mol     (Brief="Molar flow of cooling water", Unit="lbmol/h");
-        Q               as heat_rate    (Brief="Heat exchange", Unit="Btu/h");
+        C(NComp)        as conc_mol     (Brief="Molar concentration", DisplayUnit='lbmol/ft^3', Lower=0);
+        N(NComp)        as mol                  (Brief="Molar holdup", DisplayUnit='lbmol');
+        Co(NComp)       as conc_mol     (Brief="Initial molar concentration", DisplayUnit='lbmol/ft^3', Lower=0);
+        Fo(NComp)       as flow_mol     (Brief="Initial molar flow", DisplayUnit='lbmol/h');
+        T               as temperature  (Brief="Reactor temperature", DisplayUnit='degR');
+        To              as temperature  (Brief="Initial reactor temperature", DisplayUnit='degR');
+        Ta2             as temperature  (Brief="Temperature of heat exchange", DisplayUnit='degR');
+        r(NComp)        as reaction_mol (Brief="Rate of reaction", DisplayUnit='lbmol/ft^3/h');
+        k               as frequency    (Brief="Specific rate of reaction", DisplayUnit='1/h', Upper=1e5);
+        mc              as flow_mol     (Brief="Molar flow of cooling water", DisplayUnit='lbmol/h');
+        Q               as heat_rate    (Brief="Heat exchange", DisplayUnit='Btu/h');
         Theta(NComp)as Real                     (Brief="Parameter Theta");
-        vo              as flow_vol     (Brief="Volumetric flow", Unit="ft^3/h");
-        tau             as time_h       (Brief="Residence time", Unit="h");
+        vo              as flow_vol     (Brief="Volumetric flow", DisplayUnit='ft^3/h');
+        tau             as time_h       (Brief="Residence time", DisplayUnit='h');
 
         EQUATIONS
@@ -122 +122 @@
         stoic = [-1, -1, 1, 0]; # A + 2B -> C + D
         
-        UA      = 16000*"Btu/h/degR";
-        V       = 500*"gal";
-        Ta1     = (60 + 460)*"degR";
-        DH      =-3.6e4*"Btu/lbmol";
-        rho     = [0.923, 3.45, 1, 1.54]*"lbmol/ft^3";
-        cp      = [35, 18, 46, 19.5]*"Btu/lbmol/degR";
-        
-        ko = 16.96e12*"1/h";
-        E  = -32400*"Btu/lbmol";
+        UA      = 16000*'Btu/h/degR';
+        V       = 500*'gal';
+        Ta1     = (60 + 460)*'degR';
+        DH      =-3.6e4*'Btu/lbmol';
+        rho     = [0.923, 3.45, 1, 1.54]*'lbmol/ft^3';
+        cp      = [35, 18, 46, 19.5]*'Btu/lbmol/degR';
+        
+        ko = 16.96e12*'1/h';
+        E  = -32400*'Btu/lbmol';
 end
 
@@ -140 +140 @@
         SPECIFY
         "Inlet molar flow"
-        CSTR.Fo = [80, 1000, 0, 100]*"lbmol/h";
+        CSTR.Fo = [80, 1000, 0, 100]*'lbmol/h';
         "Inlet reactor temperature"
-        CSTR.To = (75 + 460)*"degR";
+        CSTR.To = (75 + 460)*'degR';
         "Molar flow of cooling water"
-        CSTR.mc = 1e3*"lbmol/h";
+        CSTR.mc = 1e3*'lbmol/h';
 
         INITIAL
         "Molar concentration"
-        CSTR.C  = [0, 3.45, 0, 0]*"lbmol/ft^3";
+        CSTR.C  = [0, 3.45, 0, 0]*'lbmol/ft^3';
         "Reactor temperature"
         CSTR.T  = CSTR.To;
 
         OPTIONS
-        time = [0:0.05:4]*"h";
+        TimeStep = 0.05;
+        TimeEnd = 4;
+        TimeUnit = 'h';
 end
 
@@ -185 +187 @@
         SPECIFY
         "Inlet molar flow"
-        CSTR.Fo = [80, 1000, 0, 100]*"lbmol/h";
+        CSTR.Fo = [80, 1000, 0, 100]*'lbmol/h';
         "Inlet reactor temperature"
-        CSTR.To = (70 + 460)*"degR"; # Reduction of temperature: 75°F to 70°F
+        CSTR.To = (70 + 460)*'degR'; # Reduction of temperature: 75°F to 70°F
         "Molar flow of cooling water"
-        CSTR.mc = 1e3*"lbmol/h";
+        CSTR.mc = 1e3*'lbmol/h';
 
         INITIAL
         "Molar concentration"
-        CSTR.C = [0.039, 2.12, 0.143, 0.226]*"lbmol/ft^3"; # Final values of SS in 9-4
+        CSTR.C = [0.039, 2.12, 0.143, 0.226]*'lbmol/ft^3'; # Final values of SS in 9-4
         "Reactor temperature"
-        CSTR.T = (138.5 + 460)*"degR"; # Final values of SS in 9-4
+        CSTR.T = (138.5 + 460)*'degR'; # Final values of SS in 9-4
 
         OPTIONS
-        time = [0:0.05:4]*"h";
+        TimeStep = 0.05;
+        TimeEnd = 4;
+        TimeUnit = 'h';
 end
 
@@ -210 +214 @@
         Tsp                     as temperature  (Brief="Reference temperature");
         mco             as flow_mol     (Brief="Molar flow of cooling water");
-        kc                      as Real                 (Brief="Gain", Unit="lbmol/degR/h");
+        kc                      as Real                 (Brief="Gain", Unit='lbmol/degR/h');
         
         VARIABLES
-        I                       as Real                 (Brief="Integral action", Unit="degR*h");
+        I                       as Real                 (Brief="Integral action", Unit='degR*h');
         X                       as fraction             (Brief="Fraction conversion", Lower=0);
         
@@ -233 +237 @@
         
         SET
-        CSTR.mco = 1000*"lbmol/h";
-        CSTR.Tsp = (138 + 460)*"degR";
-        CSTR.kc = 8.5;
+        CSTR.mco = 1000*'lbmol/h';
+        CSTR.Tsp = (138 + 460)*'degR';
+        CSTR.kc = 8.5*'lbmol/degR/h';
 
         SPECIFY
         "Inlet molar flow"
-        CSTR.Fo = [80, 1000, 0, 100]*"lbmol/h";
+        CSTR.Fo = [80, 1000, 0, 100]*'lbmol/h';
         "Inlet reactor temperature"
-        CSTR.To = (70 + 460)*"degR";
+        CSTR.To = (70 + 460)*'degR';
 
         INITIAL
         "Molar concentration"
-        CSTR.C = [0.03789, 2.12, 0.143, 0.2265]*"lbmol/ft^3";
+        CSTR.C = [0.03789, 2.12, 0.143, 0.2265]*'lbmol/ft^3';
         "Reactor temperature"
-        CSTR.T = (138.53 + 460)*"degR";
+        CSTR.T = (138.53 + 460)*'degR';
         "Integral action"
-        CSTR.I = 0;
+        CSTR.I = 0*'degR*h';
 
         OPTIONS
-        time = [0:0.01:4]*"h";
-end
+        TimeStep = 0.01;
+        TimeEnd = 4;
+        TimeUnit = 'h';
+end