Changeset 697 for branches/gui


Ignore:
Timestamp:
Dec 11, 2008, 10:23:51 AM (14 years ago)
Author:
Rafael de Pelegrini Soares
Message:

Finished high index Flash problem

File:
1 edited

Legend:

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

    r696 r697  
    9595       
    9696        EQUATIONS
    97         if time < 1 * 'h' then
    98                 fl.z = [0.45, 0.35, 0.2];
    99         else
    100                 fl.z = [0.55, 0.25, 0.2];
    101         end
    102 
     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
    103103
    104104        SPECIFY
    105         # Feed conditions
     105        # Feed condition
    106106        fl.F = 1 * 'kmol/h';
     107        # Steady-state 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';
    107121       
    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)
    115128        #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
    118132        INITIAL
    119133        #fl.T = (80+273.15) * 'K';
    120134        #fl.nt = 1 * 'kmol';
    121         fl.x(1:2) = [0.55, 0.25];
     135
     136        # Steady-state 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';
    122142       
    123143        OPTIONS
    124         TimeEnd = 2;
    125         TimeStep = 0.01;
     144        TimeEnd = 4;
     145        TimeStep = 0.05;
    126146        TimeUnit = 'h';
    127147
    128         DAESolver(File="dasslc");
     148        #DAESolver(File="dasslc"); # slow integration
     149        DAESolver(File="mebdf");  # much faster
    129150       
    130151        Dynamic = true;
Note: See TracChangeset for help on using the changeset viewer.