Ignore:
Timestamp:
Jan 14, 2007, 11:01:46 PM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Add perturbation in the Column_ctrl example.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • mso/sample/stage_separators/sample_column.mso

    r111 r112  
    312312        Pad as Real (Brief="Dimensionless pressure");
    313313        Tad as Real (Brief="Dimensionless temperature");
    314        
     314        RR      as positive (Brief="Reflux ratio");
     315
    315316        SET
    316317        PP.Components = [ "isobutane", "n-pentane", "propylene",
     
    367368        PC.Options.clip = 1;
    368369        PC.Options.autoMan = 0;
    369         PC.Parameters.intTime = 6*"s";
     370        PC.Parameters.intTime = 50*"s";
    370371        PC.Parameters.gain = 0.5;
    371372        PC.Parameters.derivTime = 1*"s";
     
    411412        col.reb.OutletL.F = Frmin + (Frmax-Frmin) * LCbot.Ports.output;
    412413
     414        RR * (col.cond.OutletV.F + col.sptop.Outlet1.F) = col.sptop.Outlet2.F;
     415       
     416        if time < 1 * "h" then
     417                col.sptop.Outlet2.F = 70 * "kmol/h"; # reflux
     418        else
     419                col.sptop.Outlet2.F = 85 * "kmol/h"; # reflux
     420        end
     421
    413422        SPECIFY
    414423        feed.F = 113.4 * "kmol/h";
     
    424433        zero.h = 0 * "J/mol";
    425434       
    426         col.sptop.Outlet2.F = 85 * "kmol/h"; # reflux
    427435        Qr = 3e6 * "kJ/h";
    428436        col.pump1.dP = 16 * "kPa";
Note: See TracChangeset for help on using the changeset viewer.