Changeset 697 for branches

Timestamp:

Dec 11, 2008, 10:23:51 AM (15 years ago)

Author:

Rafael de Pelegrini Soares

Message:

Finished high index Flash problem

File:

• branches/gui/sample/miscellaneous/sample_high_index_optimal.mso (modified) (1 diff)

branches/gui/sample/miscellaneous/sample_high_index_optimal.mso

@@ -95 +95 @@
         
         EQUATIONS
-        if time < 1 * 'h' then
-                fl.z = [0.45, 0.35, 0.2];
-        else
-                fl.z = [0.55, 0.25, 0.2];
-        end
-
+        # Disturb
+        #if time < 0.5 * 'h' then
+        #       fl.z = [0.45, 0.35, 0.2];
+        #else
+        #       fl.z = [0.55, 0.25, 0.2];
+        #end
 
         SPECIFY
-        # Feed conditions
+        # Feed condition
         fl.F = 1 * 'kmol/h';
+        # Steady-state feed
+        #fl.z = [0.45, 0.35, 0.2];
+        # Disturb on feed composition
+        fl.z = [0.55, 0.25, 0.2];
+
+        # Desired production of y1 (index 2 - will determine the "perfect" heat profile)
+        fl.V = 0.1 * 'kmol/h';
+        fl.y(1) = 0.7;
+        fl.nt = 1 * 'kmol';
+
+        # Default specification (index 1 - heat is given)
+        #fl.V = 0.1 * 'kmol/h';
+        #fl.Q = 0.37 * 'kW';
+        #fl.nt = 1 * 'kmol';
         
-        fl.V = 0.1 * 'kmol/h';
-        #fl.L = 0.2 * 'kmol/h';
-        fl.nt = 1 * 'kmol';
-        
-        # Operating conditions
-        #fl.T = (80+273.15) * 'K';
-        #fl.Q = 2.9 * 'kW';
+        # Fixed Temperature (index 2 - will determine the "perfect" heat profile)
+        #fl.T = (90+273.15) * 'K';
+        #fl.V = 0.1 * 'kmol/h';
+        #fl.nt = 1 * 'kmol';
+
+        # Fixed Pressure (index 2 - will determine the "perfect" heat profile)
         #fl.P = 1.2 * 'atm';
-        fl.y(1) = 0.8;
-        
+        #fl.V = 0.1 * 'kmol/h';
+        #fl.nt = 1 * 'kmol';
+
         INITIAL
         #fl.T = (80+273.15) * 'K';
         #fl.nt = 1 * 'kmol';
-        fl.x(1:2) = [0.55, 0.25];
+
+        # Steady-state composition with feed = [0.45, 0.35, 0.2];
+        fl.x(1) = 0.4222;
+        fl.x(2) = 0.3622;
+
+        # steady state compositions
+        # diff(fl.n(1:2)) = 0 * 'mol/s';
         
         OPTIONS
-        TimeEnd = 2;
-        TimeStep = 0.01;
+        TimeEnd = 4;
+        TimeStep = 0.05;
         TimeUnit = 'h';
 
-        DAESolver(File="dasslc");
+        #DAESolver(File="dasslc"); # slow integration
+        DAESolver(File="mebdf");  # much faster
         
         Dynamic = true;