Changeset 38 for mso/eml/stage_separators/column.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/column.mso (modified) (1 diff)

mso/eml/stage_separators/column.mso

@@ -909 +909 @@
         spbottom.Outlet2 to reb.InletL;
 end
+
+# Reactive Distillation!!
+Model ReactiveDistillation
+        PARAMETERS
+ext PP as CalcObject;
+ext NComp as Integer;
+        NTrays as Integer(Brief="Number of trays", Default=2);
+        
+        VARIABLES
+        trays(NTrays) as trayReact;
+        cond as condenserReact;
+        reb as reboilerReact;
+        sp as splitter;
+        p as pump;
+        
+        EQUATIONS
+        if ( reb.OutletV.P > 1 * "atm" ) then
+                "Pressure Drop through the tray"
+                reb.OutletV.F = trays(1).Ah/reb.vV * sqrt((reb.OutletV.P - 1*"atm") / (0.15*reb.rhoV) );
+        else
+                "Prato selado"
+                reb.OutletV.F = 0.0 * "mol/s";
+        end
+        
+        CONNECTIONS
+        #vapor
+        reb.OutletV to trays([NTrays]).InletV;
+        trays([2:NTrays]).OutletV to trays([1:NTrays-1]).InletV;
+        trays(1).OutletV to cond.InletV;
+
+        #liquid
+        cond.OutletL to sp.Inlet;       
+        sp.Outlet2 to p.Inlet;
+        p.Outlet to trays(1).InletL;
+        trays([1:NTrays-1]).OutletL to trays([2:NTrays]).InletL;
+        trays(NTrays).OutletL to reb.InletL;
+end