# Changeset 697 for branches/gui/sample

Ignore:
Timestamp:
Dec 11, 2008, 10:23:51 AM (14 years ago)
Message:

Finished high index Flash problem

File:
1 edited

Unmodified
Added
Removed
• ## branches/gui/sample/miscellaneous/sample_high_index_optimal.mso

 r696 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;
Note: See TracChangeset for help on using the changeset viewer.