Ignore:
Timestamp:
Oct 23, 2006, 5:26:39 PM (16 years ago)
Author:
Paula Bettio Staudt
Message:

Included reactive distillation column

File:
1 edited

Legend:

Unmodified
Added
Removed
  • mso/eml/stage_separators/column.mso

    r37 r38  
    909909        spbottom.Outlet2 to reb.InletL;
    910910end
     911
     912# Reactive Distillation!!
     913Model ReactiveDistillation
     914        PARAMETERS
     915ext PP as CalcObject;
     916ext NComp as Integer;
     917        NTrays as Integer(Brief="Number of trays", Default=2);
     918       
     919        VARIABLES
     920        trays(NTrays) as trayReact;
     921        cond as condenserReact;
     922        reb as reboilerReact;
     923        sp as splitter;
     924        p as pump;
     925       
     926        EQUATIONS
     927        if ( reb.OutletV.P > 1 * "atm" ) then
     928                "Pressure Drop through the tray"
     929                reb.OutletV.F = trays(1).Ah/reb.vV * sqrt((reb.OutletV.P - 1*"atm") / (0.15*reb.rhoV) );
     930        else
     931                "Prato selado"
     932                reb.OutletV.F = 0.0 * "mol/s";
     933        end
     934       
     935        CONNECTIONS
     936        #vapor
     937        reb.OutletV to trays([NTrays]).InletV;
     938        trays([2:NTrays]).OutletV to trays([1:NTrays-1]).InletV;
     939        trays(1).OutletV to cond.InletV;
     940
     941        #liquid
     942        cond.OutletL to sp.Inlet;       
     943        sp.Outlet2 to p.Inlet;
     944        p.Outlet to trays(1).InletL;
     945        trays([1:NTrays-1]).OutletL to trays([2:NTrays]).InletL;
     946        trays(NTrays).OutletL to reb.InletL;
     947end
Note: See TracChangeset for help on using the changeset viewer.