Changeset 516

Timestamp:

May 15, 2008, 6:37:10 PM (15 years ago)

Author:

Argimiro Resende Secchi

Message:

Updating tray rate.

File:

• branches/rate/eml/stage_separators/tray.mso (modified) (5 diffs)

branches/rate/eml/stage_separators/tray.mso

@@ -494 +494 @@
         
         VARIABLES
-        NL(NComp) as flux_mol   (Brief = "Stream Molar Flux Rate on Liquid Phase");
-        NV(NComp) as flux_mol   (Brief = "Stream Molar Flux Rate on Vapour Phase");
+        NL(NComp) as flow_mol   (Brief = "Stream Molar Rate on Liquid Phase");
+        NV(NComp) as flow_mol   (Brief = "Stream Molar Rate on Vapour Phase");
         T as temperature                (Brief = "Stream Temperature");
         P as pressure                   (Brief = "Stream Pressure");
@@ -503 +503 @@
         htL as heat_trans_coeff (Brief = "Heat Transference Coefficient on Liquid Phase");
         htV as heat_trans_coeff (Brief = "Heat Transference Coefficient on Vapour Phase");      
-        E_liq as heat_flux      (Brief = "Liquid Energy on interface");
-    E_vap as heat_flux      (Brief = "Vapour Energy on interface");     
+        E_liq as heat_rate      (Brief = "Liquid Energy Rate at interface");
+    E_vap as heat_rate      (Brief = "Vapour Energy Rate at interface");        
         hL as enth_mol          (Brief = "Liquid Molar Enthalpy");
         hV as enth_mol          (Brief = "Vapour Molar Enthalpy");
@@ -565 +565 @@
         "Component Molar Balance"
         diff(M_liq)=Inlet.F*Inlet.z + InletL.F*InletL.z 
-        - OutletL.F*OutletL.z + interf.a*interf.NL;
+        - OutletL.F*OutletL.z + interf.NL;
         
         diff(M_vap)=InletFV.F*InletFV.z + InletV.F*InletV.z 
-        - OutletV.F*OutletV.z - interf.a*interf.NV;
+        - OutletV.F*OutletV.z - interf.NV;
         
         "Energy Balance"
         diff(E_liq) = Inlet.F*Inlet.h + InletL.F*InletL.h
-                - OutletL.F*OutletL.h  + Q + interf.a*interf.E_liq;
+                - OutletL.F*OutletL.h  + Q + interf.E_liq;
         
         diff(E_vap) = InletFV.F*InletFV.h + InletV.F*InletV.h
-                - OutletV.F*OutletV.h  - interf.a*interf.E_vap;
+                - OutletV.F*OutletV.h  - interf.E_vap;
         
         "Molar Holdup"
@@ -587 +587 @@
         E_vap = MV*(OutletV.h - OutletV.P*vV);
         
-        "Energy on interface"
-        interf.E_liq = interf.htL*(interf.T-OutletL.T)+sum(interf.NL)*interf.hL;        
-        
-        interf.E_vap = interf.htV*(OutletV.T-interf.T)+sum(interf.NV)*interf.hV;
+        "Energy Rate through the interface"
+        interf.E_liq = interf.htL*interf.a*(interf.T-OutletL.T)+sum(interf.NL)*interf.hL;       
+        
+        interf.E_vap = interf.htV*interf.a*(OutletV.T-interf.T)+sum(interf.NV)*interf.hV;
         
         "Mass Conservation"
-        interf.NL=interf.NV;
+        interf.NL = interf.NV;
         
         "Energy Conservation"
@@ -619 +619 @@
         Level = ML*vL/Ap;
 
-        "Total Mass Transfer Fluxes"
-        interf.NL(1:NC1)=sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
+        "Total Mass Transfer Rates"
+        interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
                 OutletL.z(1:NC1)*sum(interf.NL);
         
-        interf.NV(1:NC1)=sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
+        interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
                 OutletV.z(1:NC1)*sum(interf.NV);