# Changeset 541 for branches

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

updated pipe model (testing)

File:
1 edited

Unmodified
Removed
• ## branches/tests/eml/pressure_changers/pipe.mso

 r540 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 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 Liquid.Holdup(i) = 1; "Lockhart-martinelli Two Phase Multiplier" XLM(i)  = 1; PhiLM(i) = 1; "Beggs and Brill Line 1" BBLine1(i) = 1; Liquid.Holdup(i) = 0; "Lockhart-martinelli Two Phase Multiplier" XLM(i)  = 1; PhiLM(i) = 1; "Beggs and Brill Line 1" BBLine1(i) = 1; case "Beggs-Brill": "Lockhart-martinelli Two Phase Multiplier" XLM(i)  = 1; PhiLM(i) = 1; "Beggs and Brill Line 1" 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; "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] "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; 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]; Tube.Hrise =    0*'deg'; # working only in the horizontal , must be updated Tube.Correlation = "Beggs-Brill"; #Tube.Correlation = "Lockhart-Martinelli";
Note: See TracChangeset for help on using the changeset viewer.