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

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/tests/eml/pressure_changers/pipe.mso
r540 r541 107 107 dPdLfric(N+1) as positive (Brief = "Friction Gradient",Lower = 0, Unit = 'Pa/m', DisplayUnit = 'kPa/m'); 108 108 dPelvdL(N+1) as positive (Brief = "Elevation Gradient", Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m'); 109 dPaccdL(N+1) as positive (Brief = "Acceleration Gradient",Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m');109 #dPaccdL(N+1) as positive (Brief = "Acceleration Gradient",Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m'); 110 110 111 111 #Beggs&Brill Method … … 120 120 Y(N+1) as positive (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow"); 121 121 S(N+1) as Real (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow"); 122 122 123 #Lockhart&Martinelli Method 124 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 125 126 XLM(N+1) as Real (Brief="Lockhart and Martinelli Parameter  Square", Lower=1E6); 127 PhiLM(N+1) as Real (Brief="Lockhart and Martinelli Two phase Multiplier  Square", Lower=1E3); 128 123 129 EQUATIONS 124 130 … … 160 166 Liquid.Holdup(i) = 1; 161 167 168 "Lockhartmartinelli Two Phase Multiplier" 169 XLM(i) = 1; 170 171 PhiLM(i) = 1; 172 162 173 "Beggs and Brill Line 1" 163 174 BBLine1(i) = 1; … … 200 211 Liquid.Holdup(i) = 0; 201 212 213 "Lockhartmartinelli Two Phase Multiplier" 214 XLM(i) = 1; 215 216 PhiLM(i) = 1; 217 202 218 "Beggs and Brill Line 1" 203 219 BBLine1(i) = 1; … … 240 256 241 257 case "BeggsBrill": 258 259 "Lockhartmartinelli Two Phase Multiplier" 260 XLM(i) = 1; 261 262 PhiLM(i) = 1; 242 263 243 264 "Beggs and Brill Line 1" … … 351 372 BBLine4 (i)= 1; 352 373 374 "Lockhartmartinelli Two Phase Multiplier" 375 XLM(i)*Vapour.rho(i)*Vapour.Vel(i)^2 = Liquid.rho(i)*Liquid.Vel(i)^2; 376 377 PhiLM(i) = 1+20/sqrt(XLM(i))+1/XLM(i); 378 353 379 "Liquid Holdup" 354 Liquid.Holdup(i)*NFR(i)^0.0868 = 0.98*Liquid.Vsfrac(i)^(0.4846);#usando beggsbrill 355 Atrans(i) = 1; 380 Liquid.Holdup(i) = (1/abs(PhiLM(i)))^(1/3); 381 382 Atrans(i) = 1; 383 356 384 Map(i)=20; 385 357 386 Y(i)=1; 358 387 … … 417 446 "Incremental Elevation Pressure Drop" 418 447 dPelv= rho*g*Lincr*sin(Hrise); 448 449 dPelvdL= rho*g*sin(Hrise); 419 450 420 451 "Incremental Friction Pressure Drop" 421 dPfric = (0.5*fns*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe); 422 #dPfric = (0.5*ftp*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe); 452 #dPfric = (0.5*fns*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe); 453 dPfric = PhiLM*(0.5*fns*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe); 454 455 dPdLfric = (0.5*fns*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe); 423 456 424 457 for i in [1:N+1] … … 467 500 468 501 "Incremental Acceleration Pressure Drop" 469 #dPacc(i+1) = (Fw/Apipe)^2*(1/(Liquid.rho(i+1)*Liquid.Holdup(i+1)+Vapour.rho(i+1)*Vapour.Holdup(i+1))1/470 #(Liquid.rho(i)*Liquid.Holdup(i)+Vapour.rho(i)*Vapour.Holdup(i)));471 502 dPacc(i+1) = (Fw/Apipe)^2*(1/rho(i+1)1/rho(i)); 472 #dPacc(i+1) = dPacc(1); 473 503 474 504 "Incremental Length" 475 505 Lincr(i+1) = i*abs(Lpipe)/N; … … 551 581 Feed.Outlet.F = 10 * 'kmol/h'; 552 582 Feed.Outlet.P = 3* 'atm'; 553 Feed.Outlet.T = ( 30+273.15) * 'K';583 Feed.Outlet.T = (60+273.15) * 'K'; 554 584 Feed.Outlet.z = [0.5,0.2,0.3]; 555 585 … … 558 588 Tube.Hrise = 0*'deg'; # working only in the horizontal , must be updated 559 589 Tube.Correlation = "BeggsBrill"; 590 #Tube.Correlation = "LockhartMartinelli"; 560 591 561 592
Note: See TracChangeset
for help on using the changeset viewer.