Changeset 507 for branches/rate

Timestamp:

Apr 28, 2008, 4:17:48 PM (15 years ago)

Author:

Argimiro Resende Secchi

Message:

Updating tray rate model.

Location:

branches/rate

Files:

• eml/stage_separators/tray.mso (modified) (3 diffs)
• sample/stage_separators/sample_tray.mso (modified) (2 diffs)

branches/rate/eml/stage_separators/tray.mso

@@ -507 +507 @@
         hL as enth_mol          (Brief = "Liquid Molar Enthalpy");
         hV as enth_mol          (Brief = "Vapour Molar Enthalpy");
-        ConcL as conc_mol       (Brief = "Liquid Molar Concentration");
-        ConcV as conc_mol       (Brief = "Vapour Molar Concentration");
         kL(NC1,NC1) as velocity (Brief = "Mass Transfer Coefficients");
         kV(NC1,NC1) as velocity (Brief = "Mass Transfer Coefficients");
@@ -600 +598 @@
         interf.E_liq = interf.E_vap;
         
-        
         "Mol fraction normalisation"
         sum(OutletL.z)= 1.0;
@@ -621 +618 @@
         "Level of clear liquid over the weir"
         Level = ML*vL/Ap;
-        
+
         "Total Mass Transfer Fluxes"
-        interf.NL(1:NC1)=interf.ConcL*sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))+
+        interf.NL(1:NC1)=sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
                 OutletL.z(1:NC1)*sum(interf.NL);
         
-        interf.NV(1:NC1)=interf.ConcV*sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))+
+        interf.NV(1:NC1)=sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
                 OutletV.z(1:NC1)*sum(interf.NV);
 

branches/rate/sample/stage_separators/sample_tray.mso

@@ -271 +271 @@
         t1.interf.kL=0.001*'m/s';
         t1.interf.kV=0.001*'m/s';
-        t1.ML=3*'mol';
-        t1.MV=3*'mol';
 
         SET
@@ -285 +283 @@
         t1.Q = 0 * 'kW';
         
-        
         INITIAL
         t1.OutletL.T = 290 *'K';
         t1.OutletV.T = 290 *'K';
-        #t1.Level = 0.9*'m' ;
+        t1.Level = 0.9*'m' ;
         t1.OutletL.z(1) = 0.5;
         t1.OutletV.z(1) = 0.5;
-        #t1.MV = 0.01*'mol' ;
-        
-        OPTIONS
+        t1.MV = 0.01*'mol' ;
+        
+        OPTIONS
+        TimeStep = 0.1;
         TimeEnd = 100;
-end
+        TimeUnit = 's';
+#       GuessFile = "tray_Test.rlt";
+end