Changeset 697
 Timestamp:
 Dec 11, 2008, 10:23:51 AM (14 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/gui/sample/miscellaneous/sample_high_index_optimal.mso
r696 r697 95 95 96 96 EQUATIONS 97 if time < 1 * 'h' then98 fl.z = [0.45, 0.35, 0.2];99 else100 fl.z = [0.55, 0.25, 0.2];101 end102 97 # Disturb 98 #if time < 0.5 * 'h' then 99 # fl.z = [0.45, 0.35, 0.2]; 100 #else 101 # fl.z = [0.55, 0.25, 0.2]; 102 #end 103 103 104 104 SPECIFY 105 # Feed condition s105 # Feed condition 106 106 fl.F = 1 * 'kmol/h'; 107 # Steadystate feed 108 #fl.z = [0.45, 0.35, 0.2]; 109 # Disturb on feed composition 110 fl.z = [0.55, 0.25, 0.2]; 111 112 # Desired production of y1 (index 2  will determine the "perfect" heat profile) 113 fl.V = 0.1 * 'kmol/h'; 114 fl.y(1) = 0.7; 115 fl.nt = 1 * 'kmol'; 116 117 # Default specification (index 1  heat is given) 118 #fl.V = 0.1 * 'kmol/h'; 119 #fl.Q = 0.37 * 'kW'; 120 #fl.nt = 1 * 'kmol'; 107 121 108 fl.V = 0.1 * 'kmol/h'; 109 #fl.L = 0.2 * 'kmol/h'; 110 fl.nt = 1 * 'kmol'; 111 112 # Operating conditions 113 #fl.T = (80+273.15) * 'K'; 114 #fl.Q = 2.9 * 'kW'; 122 # Fixed Temperature (index 2  will determine the "perfect" heat profile) 123 #fl.T = (90+273.15) * 'K'; 124 #fl.V = 0.1 * 'kmol/h'; 125 #fl.nt = 1 * 'kmol'; 126 127 # Fixed Pressure (index 2  will determine the "perfect" heat profile) 115 128 #fl.P = 1.2 * 'atm'; 116 fl.y(1) = 0.8; 117 129 #fl.V = 0.1 * 'kmol/h'; 130 #fl.nt = 1 * 'kmol'; 131 118 132 INITIAL 119 133 #fl.T = (80+273.15) * 'K'; 120 134 #fl.nt = 1 * 'kmol'; 121 fl.x(1:2) = [0.55, 0.25]; 135 136 # Steadystate composition with feed = [0.45, 0.35, 0.2]; 137 fl.x(1) = 0.4222; 138 fl.x(2) = 0.3622; 139 140 # steady state compositions 141 # diff(fl.n(1:2)) = 0 * 'mol/s'; 122 142 123 143 OPTIONS 124 TimeEnd = 2;125 TimeStep = 0.0 1;144 TimeEnd = 4; 145 TimeStep = 0.05; 126 146 TimeUnit = 'h'; 127 147 128 DAESolver(File="dasslc"); 148 #DAESolver(File="dasslc"); # slow integration 149 DAESolver(File="mebdf"); # much faster 129 150 130 151 Dynamic = true;
Note: See TracChangeset
for help on using the changeset viewer.