Changeset 38 for mso/eml/stage_separators/reboiler.mso

Timestamp:

Oct 23, 2006, 5:26:39 PM (17 years ago)

Author:

Paula Bettio Staudt

Message:

Included reactive distillation column

File:

• mso/eml/stage_separators/reboiler.mso (modified) (1 diff)

mso/eml/stage_separators/reboiler.mso

@@ -141 +141 @@
 end
 
+#*-------------------------------------------------------------------
+* Model of a dynamic reboiler with reaction
+*-------------------------------------------------------------------*#
+Model reboilerReact
+        PARAMETERS
+ext PP as CalcObject;
+ext NComp as Integer;
+        Across as area (Brief="Cross Section Area of reboiler");
+        V as volume (Brief="Total volume of reboiler");
+
+        stoic(NComp) as Real(Brief="Stoichiometric matrix");
+        Hr as energy_mol;
+        Pstartup as pressure;
+
+        VARIABLES
+in      Inlet as stream;                        #(Brief="Feed Stream");
+in      InletL as stream;                       #(Brief="Liquid inlet stream");
+out     OutletL as stream_therm;        #(Brief="Liquid outlet stream");
+out     OutletV as stream_therm;        #(Brief="Vapour outlet stream");
+
+        Q as heat_rate (Brief="Heat supplied");
+        M(NComp) as mol (Brief="Molar Holdup in the tray");
+        ML as mol (Brief="Molar liquid holdup");
+        MV as mol (Brief="Molar vapour holdup");
+        E as energy (Brief="Total Energy Holdup on tray");
+        vL as volume_mol (Brief="Liquid Molar Volume");
+        vV as volume_mol (Brief="Vapour Molar volume");
+        Level as length (Brief="Level of liquid phase");
+        Vol as volume;
+        startup as Real;
+        rhoV as dens_mass;
+        r as reaction_mol (Brief = "Reaction resulting ethyl acetate", Unit = "mol/l/s");
+        C(NComp) as conc_mol (Brief = "Molar concentration", Lower = -1);
+
+        EQUATIONS
+        "Molar Concentration"
+        OutletL.z = vL * C;
+        
+        "Component Molar Balance"
+        diff(M)= Inlet.F*Inlet.z + InletL.F*InletL.z
+                - OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r*ML*vL;
+        
+        "Energy Balance"
+        diff(E) = Inlet.F*Inlet.h + InletL.F*InletL.h
+                - OutletL.F*OutletL.h - OutletV.F*OutletV.h + Q + Hr * r * vL*ML;
+        
+        "Molar Holdup"
+        M = ML*OutletL.z + MV*OutletV.z; 
+        
+        "Energy Holdup"
+        E = ML*OutletL.h + MV*OutletV.h - OutletL.P*V;
+        
+        "Mol fraction normalisation"
+        sum(OutletL.z)=1.0;
+        
+        "Liquid Volume"
+        vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
+        "Vapour Volume"
+        vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);  
+        "Vapour Density"
+        rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z);
+        
+        "Level of liquid phase"
+        Level = ML*vL/Across;
+
+        Vol = ML*vL;
+        
+        "vaporization fraction "
+        OutletV.v = 1.0;
+        OutletL.v = 0.0;
+        
+        "Mechanical Equilibrium"
+        OutletL.P = OutletV.P;
+        
+        "Thermal Equilibrium"
+        OutletL.T = OutletV.T;  
+        
+        "Geometry Constraint"
+        V = ML*vL + MV*vV;              
+
+        "Chemical Equilibrium"
+        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = 
+        PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
+
+        sum(OutletL.z)=sum(OutletV.z);
+        
+end