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/reboiler.mso

    r1 r38  
    141141end
    142142
     143#*-------------------------------------------------------------------
     144* Model of a dynamic reboiler with reaction
     145*-------------------------------------------------------------------*#
     146Model reboilerReact
     147        PARAMETERS
     148ext PP as CalcObject;
     149ext NComp as Integer;
     150        Across as area (Brief="Cross Section Area of reboiler");
     151        V as volume (Brief="Total volume of reboiler");
     152
     153        stoic(NComp) as Real(Brief="Stoichiometric matrix");
     154        Hr as energy_mol;
     155        Pstartup as pressure;
     156
     157        VARIABLES
     158in      Inlet as stream;                        #(Brief="Feed Stream");
     159in      InletL as stream;                       #(Brief="Liquid inlet stream");
     160out     OutletL as stream_therm;        #(Brief="Liquid outlet stream");
     161out     OutletV as stream_therm;        #(Brief="Vapour outlet stream");
     162
     163        Q as heat_rate (Brief="Heat supplied");
     164        M(NComp) as mol (Brief="Molar Holdup in the tray");
     165        ML as mol (Brief="Molar liquid holdup");
     166        MV as mol (Brief="Molar vapour holdup");
     167        E as energy (Brief="Total Energy Holdup on tray");
     168        vL as volume_mol (Brief="Liquid Molar Volume");
     169        vV as volume_mol (Brief="Vapour Molar volume");
     170        Level as length (Brief="Level of liquid phase");
     171        Vol as volume;
     172        startup as Real;
     173        rhoV as dens_mass;
     174        r as reaction_mol (Brief = "Reaction resulting ethyl acetate", Unit = "mol/l/s");
     175        C(NComp) as conc_mol (Brief = "Molar concentration", Lower = -1);
     176
     177        EQUATIONS
     178        "Molar Concentration"
     179        OutletL.z = vL * C;
     180       
     181        "Component Molar Balance"
     182        diff(M)= Inlet.F*Inlet.z + InletL.F*InletL.z
     183                - OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r*ML*vL;
     184       
     185        "Energy Balance"
     186        diff(E) = Inlet.F*Inlet.h + InletL.F*InletL.h
     187                - OutletL.F*OutletL.h - OutletV.F*OutletV.h + Q + Hr * r * vL*ML;
     188       
     189        "Molar Holdup"
     190        M = ML*OutletL.z + MV*OutletV.z;
     191       
     192        "Energy Holdup"
     193        E = ML*OutletL.h + MV*OutletV.h - OutletL.P*V;
     194       
     195        "Mol fraction normalisation"
     196        sum(OutletL.z)=1.0;
     197       
     198        "Liquid Volume"
     199        vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     200        "Vapour Volume"
     201        vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z); 
     202        "Vapour Density"
     203        rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z);
     204       
     205        "Level of liquid phase"
     206        Level = ML*vL/Across;
     207
     208        Vol = ML*vL;
     209       
     210        "vaporization fraction "
     211        OutletV.v = 1.0;
     212        OutletL.v = 0.0;
     213       
     214        "Mechanical Equilibrium"
     215        OutletL.P = OutletV.P;
     216       
     217        "Thermal Equilibrium"
     218        OutletL.T = OutletV.T; 
     219       
     220        "Geometry Constraint"
     221        V = ML*vL + MV*vV;             
     222
     223        "Chemical Equilibrium"
     224        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =
     225        PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
     226
     227        sum(OutletL.z)=sum(OutletV.z);
     228       
     229end
Note: See TracChangeset for help on using the changeset viewer.