Changeset 507 for branches/rate/eml/stage_separators

Timestamp:

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

Author:

Argimiro Resende Secchi

Message:

Updating tray rate model.

File:

• branches/rate/eml/stage_separators/tray.mso (modified) (3 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);