Changeset 541 for branches

Timestamp:

Jun 18, 2008, 7:13:11 PM (15 years ago)

Author:

gerson bicca

Message:

updated pipe model (testing)

File:

• branches/tests/eml/pressure_changers/pipe.mso (modified) (10 diffs)

branches/tests/eml/pressure_changers/pipe.mso

@@ -107 +107 @@
         dPdLfric(N+1)           as positive             (Brief = "Friction Gradient",Lower = 0, Unit = 'Pa/m', DisplayUnit = 'kPa/m');
         dPelvdL(N+1)            as positive             (Brief = "Elevation Gradient", Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m');
-        dPaccdL(N+1)            as positive             (Brief = "Acceleration Gradient",Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m');
+        #dPaccdL(N+1)           as positive             (Brief = "Acceleration Gradient",Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m');
         
 #Beggs&Brill Method
@@ -120 +120 @@
         Y(N+1)                                  as positive             (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow");
         S(N+1)                                  as Real                         (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow");
-        
+
+#Lockhart&Martinelli Method
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+
+XLM(N+1)                                        as Real                         (Brief="Lockhart and Martinelli  Parameter - Square", Lower=1E-6);
+PhiLM(N+1)                              as Real                         (Brief="Lockhart and Martinelli  Two phase Multiplier - Square", Lower=1E-3);
+
 EQUATIONS
 
@@ -160 +166 @@
         Liquid.Holdup(i) = 1;
 
+"Lockhart-martinelli Two Phase Multiplier"
+        XLM(i)  = 1;
+
+        PhiLM(i) = 1;
+        
 "Beggs and Brill Line 1"
         BBLine1(i) = 1;
@@ -200 +211 @@
         Liquid.Holdup(i) = 0;
 
+"Lockhart-martinelli Two Phase Multiplier"
+        XLM(i)  = 1;
+
+        PhiLM(i) = 1;
+        
 "Beggs and Brill Line 1"
         BBLine1(i) = 1;
@@ -240 +256 @@
         
         case "Beggs-Brill":
+
+"Lockhart-martinelli Two Phase Multiplier"
+        XLM(i)  = 1;
+
+        PhiLM(i) = 1;
         
 "Beggs and Brill Line 1"        
@@ -351 +372 @@
         BBLine4 (i)= 1;
 
+"Lockhart-martinelli Two Phase Multiplier"
+        XLM(i)*Vapour.rho(i)*Vapour.Vel(i)^2    = Liquid.rho(i)*Liquid.Vel(i)^2;
+        
+        PhiLM(i) = 1+20/sqrt(XLM(i))+1/XLM(i);
+
 "Liquid Holdup"
-        Liquid.Holdup(i)*NFR(i)^0.0868 = 0.98*Liquid.Vsfrac(i)^(0.4846);#usando beggs-brill
-        Atrans(i) = 1;
+        Liquid.Holdup(i) = (1/abs(PhiLM(i)))^(1/3);
+        
+        Atrans(i) = 1;
+        
         Map(i)=20;
+        
         Y(i)=1;
         
@@ -417 +446 @@
 "Incremental Elevation Pressure Drop"
         dPelv= rho*g*Lincr*sin(Hrise);
+        
+        dPelvdL= rho*g*sin(Hrise);
 
 "Incremental Friction Pressure Drop"
-        dPfric = (0.5*fns*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe);
-        #dPfric = (0.5*ftp*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe);
+        #dPfric = (0.5*fns*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe);
+        dPfric = PhiLM*(0.5*fns*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe);
+        
+        dPdLfric = (0.5*fns*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe);
         
 for i in [1:N+1]
@@ -467 +500 @@
 
 "Incremental Acceleration Pressure Drop"
-        #dPacc(i+1) = (Fw/Apipe)^2*(1/(Liquid.rho(i+1)*Liquid.Holdup(i+1)+Vapour.rho(i+1)*Vapour.Holdup(i+1))-1/
-        #(Liquid.rho(i)*Liquid.Holdup(i)+Vapour.rho(i)*Vapour.Holdup(i)));
         dPacc(i+1) = (Fw/Apipe)^2*(1/rho(i+1)-1/rho(i));
-        #dPacc(i+1) = dPacc(1);
-
+        
 "Incremental Length"
         Lincr(i+1) = i*abs(Lpipe)/N;
@@ -551 +581 @@
         Feed.Outlet.F = 10 * 'kmol/h';
         Feed.Outlet.P = 3* 'atm';
-        Feed.Outlet.T = (30+273.15) * 'K';
+        Feed.Outlet.T = (60+273.15) * 'K';
         Feed.Outlet.z = [0.5,0.2,0.3];
 
@@ -558 +588 @@
         Tube.Hrise =    0*'deg'; # working only in the horizontal , must be updated
         Tube.Correlation = "Beggs-Brill";
+        #Tube.Correlation = "Lockhart-Martinelli";