Changeset 541


Ignore:
Timestamp:
Jun 18, 2008, 7:13:11 PM (14 years ago)
Author:
gerson bicca
Message:

updated pipe model (testing)

File:
1 edited

Legend:

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

    r540 r541  
    107107        dPdLfric(N+1)           as positive             (Brief = "Friction Gradient",Lower = 0, Unit = 'Pa/m', DisplayUnit = 'kPa/m');
    108108        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');
    110110       
    111111#Beggs&Brill Method
     
    120120        Y(N+1)                                  as positive             (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow");
    121121        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
     126XLM(N+1)                                        as Real                         (Brief="Lockhart and Martinelli  Parameter - Square", Lower=1E-6);
     127PhiLM(N+1)                              as Real                         (Brief="Lockhart and Martinelli  Two phase Multiplier - Square", Lower=1E-3);
     128
    123129EQUATIONS
    124130
     
    160166        Liquid.Holdup(i) = 1;
    161167
     168"Lockhart-martinelli Two Phase Multiplier"
     169        XLM(i)  = 1;
     170
     171        PhiLM(i) = 1;
     172       
    162173"Beggs and Brill Line 1"
    163174        BBLine1(i) = 1;
     
    200211        Liquid.Holdup(i) = 0;
    201212
     213"Lockhart-martinelli Two Phase Multiplier"
     214        XLM(i)  = 1;
     215
     216        PhiLM(i) = 1;
     217       
    202218"Beggs and Brill Line 1"
    203219        BBLine1(i) = 1;
     
    240256       
    241257        case "Beggs-Brill":
     258
     259"Lockhart-martinelli Two Phase Multiplier"
     260        XLM(i)  = 1;
     261
     262        PhiLM(i) = 1;
    242263       
    243264"Beggs and Brill Line 1"       
     
    351372        BBLine4 (i)= 1;
    352373
     374"Lockhart-martinelli 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
    353379"Liquid Holdup"
    354         Liquid.Holdup(i)*NFR(i)^0.0868 = 0.98*Liquid.Vsfrac(i)^(0.4846);#usando beggs-brill
    355         Atrans(i) = 1;
     380        Liquid.Holdup(i) = (1/abs(PhiLM(i)))^(1/3);
     381       
     382        Atrans(i) = 1;
     383       
    356384        Map(i)=20;
     385       
    357386        Y(i)=1;
    358387       
     
    417446"Incremental Elevation Pressure Drop"
    418447        dPelv= rho*g*Lincr*sin(Hrise);
     448       
     449        dPelvdL= rho*g*sin(Hrise);
    419450
    420451"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);
    423456       
    424457for i in [1:N+1]
     
    467500
    468501"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)));
    471502        dPacc(i+1) = (Fw/Apipe)^2*(1/rho(i+1)-1/rho(i));
    472         #dPacc(i+1) = dPacc(1);
    473 
     503       
    474504"Incremental Length"
    475505        Lincr(i+1) = i*abs(Lpipe)/N;
     
    551581        Feed.Outlet.F = 10 * 'kmol/h';
    552582        Feed.Outlet.P = 3* 'atm';
    553         Feed.Outlet.T = (30+273.15) * 'K';
     583        Feed.Outlet.T = (60+273.15) * 'K';
    554584        Feed.Outlet.z = [0.5,0.2,0.3];
    555585
     
    558588        Tube.Hrise =    0*'deg'; # working only in the horizontal , must be updated
    559589        Tube.Correlation = "Beggs-Brill";
     590        #Tube.Correlation = "Lockhart-Martinelli";
    560591       
    561592
Note: See TracChangeset for help on using the changeset viewer.