Changeset 46

Timestamp:

Nov 7, 2006, 2:47:55 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

Fixed some sample files

Location:

mso

Files:

• eml/stage_separators/column.mso (modified) (2 diffs)
• eml/stage_separators/reboiler.mso (modified) (1 diff)
• sample/controllers/sample_tanks_pid.mso (modified) (1 diff)
• sample/miscellaneous/sample_pend.mso (modified) (1 diff)
• sample/reactors/sample_pfr.mso (modified) (1 diff)
• sample/stage_separators/sample_column.mso (modified) (1 diff)
• sample/stage_separators/sample_tank.mso (modified) (6 diffs)

mso/eml/stage_separators/column.mso

@@ -167 +167 @@
         "Pressure Drop through the tray"
         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 * 
-                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
+                trays(top).Ah/trays(top).vV /2* sqrt(2*(trays(top).OutletV.P -
                 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
                 
@@ -173 +173 @@
             (trays([top+topdown:topdown:bot]).OutletV.P -
             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
-            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
+            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV /2*
             sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /

mso/eml/stage_separators/reboiler.mso

@@ -139 +139 @@
         "Vapourisation Fraction"
         OutletV.v = 1.0;
+end
+
+Model reboilerSteady_fakeH
+        PARAMETERS
+ext PP as CalcObject;
+ext NComp as Integer;
+        DP as press_delta (Brief="Pressure Drop in the reboiler");
+        k as Real (Brief = "Flow Constant", Unit="mol/J");
+        
+        VARIABLES
+in      InletL as stream; #(Brief="Liquid inlet stream");
+out     OutletV as stream; #(Brief="Vapour outlet stream");
+in      Q as heat_rate (Brief="Heat supplied");
+
+        EQUATIONS
+        "Molar Balance"
+        InletL.F = OutletV.F;
+        InletL.z = OutletV.z;
+        
+        "Energy Balance"
+        InletL.F*InletL.h + Q = OutletV.F*OutletV.h;
+        
+        "Pressure"
+        DP = InletL.P - OutletV.P;
+
+        "Fake Vapourisation Fraction"
+        OutletV.v = 1.0;
+        
+        "Fake output temperature"
+        OutletV.T = 300*"K";
+        
+        "Pressure Drop through the reboiler"
+        OutletV.F = k*Q;
 end
 

mso/sample/controllers/sample_tanks_pid.mso

@@ -127 +127 @@
         
         OPTIONS
+        NLASolver = "sundials";
         time = [0:0.1:30] * "h" ;
 end

mso/sample/miscellaneous/sample_pend.mso

@@ -64 +64 @@
         time = [0:0.1:30];
         integration = "original"; # original, index0 or index1
-        accuracyType = "scalar"; # scalar or vector     
         relativeAccuracy = 1e-8;
         absoluteAccuracy = 1e-9;

mso/sample/reactors/sample_pfr.mso

@@ -78 +78 @@
         #mode = "steady";
         time = [0:0.05:3] * "s";
+        NLASolver = "sundials";
+        #DAESolver = "dassl";
 end

mso/sample/stage_separators/sample_column.mso

@@ -73 +73 @@
         OPTIONS
         #relativeAccuracy = 1e-3;
+        NLASolver = "sundials";
         time = [0:10:1000];
 end

mso/sample/stage_separators/sample_tank.mso

@@ -11 +11 @@
         NComp = PP.NumberOfComponents;
         
+        VARIABLES
+        Qtank as heat_rate (Brief="Heat rate supplied to tank");
+        
         DEVICES
         t as tank;
@@ -17 +20 @@
         CONNECTIONS
         s to t.Inlet;
+        Qtank to t.Q;
         
         SPECIFY
@@ -25 +29 @@
         s.v = 0.698; 
         
-        t.Q = 0 * "J/s";
+        Qtank = 0 * "J/s";
         t.Outlet.F = 179 * "kmol/h";
         
@@ -50 +54 @@
         PP.VapourModel = "PR";
         NComp = PP.NumberOfComponents;
+
+        VARIABLES
+        Qtank as heat_rate (Brief="Heat rate supplied to tank");
         
         DEVICES
@@ -57 +64 @@
         CONNECTIONS
         s to t.Inlet;
+        Qtank to t.Q;
         
         SPECIFY
@@ -65 +73 @@
         s.v = 0.368;
 
-        t.Q = 0 * "J/s";
+        Qtank = 0 * "J/s";
         t.Outlet.F = 490 * "kmol/h";