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

Timestamp:

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

Author:

Paula Bettio Staudt

Message:

Included reactive distillation column

File:

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

mso/eml/stage_separators/condenser.mso

@@ -122 +122 @@
         OutletL.v = 0.0;
 end
+
+#*-------------------------------------------------------------------
+* Condenser with reaction in liquid phase
+*--------------------------------------------------------------------*#
+Model condenserReact
+        PARAMETERS
+ext PP as CalcObject;
+ext NComp as Integer;
+        V as volume (Brief="Condenser total volume");
+        Across as area (Brief="Cross Section Area of reboiler");
+
+        stoic(NComp) as Real(Brief="Stoichiometric matrix");
+        Hr as energy_mol;
+        Pstartup as pressure;
+
+        VARIABLES
+in      InletV as stream;                       #(Brief="Vapour inlet stream");
+out     OutletL as stream_therm;        #(Brief="Liquid outlet stream");
+out     OutletV as stream_therm;        #(Brief="Vapour outlet stream");
+
+        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");
+        Q as heat_rate (Brief="Heat supplied");
+        Vol as volume;
+        r as reaction_mol (Brief = "Reaction rate", 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) = InletV.F*InletV.z - OutletL.F*OutletL.z
+                                - OutletV.F*OutletV.z + stoic*r*ML*vL;
+
+        "Energy Balance"
+        diff(E) = InletV.F*InletV.h - OutletL.F*OutletL.h
+                                - OutletV.F*OutletV.h + Q + Hr * r * ML*vL;
+
+        "Molar Holdup"
+        M = ML*OutletL.z + MV*OutletV.z; 
+        
+        "Energy Holdup"
+        E = ML*OutletL.h + MV*OutletV.h - OutletV.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);
+
+        "Thermal Equilibrium"
+        OutletL.T = OutletV.T;
+
+        "Mechanical Equilibrium"
+        OutletV.P = OutletL.P;
+
+        "Geometry Constraint"
+        V = ML*vL + MV*vV;
+
+        Vol = ML*vL;
+        
+        "Level of liquid phase"
+        Level = ML*vL/Across;
+        
+        "Vapourisation Fraction"
+        OutletL.v = 0.0;
+        OutletV.v = 1.0;
+        
+        "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