Changeset 540 for branches


Ignore:
Timestamp:
Jun 18, 2008, 3:26:58 PM (14 years ago)
Author:
gerson bicca
Message:

updated pipe model

File:
1 edited

Legend:

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

    r531 r540  
    11
    22using "streams";
     3# Not Finished yet !!!!
     4# Do Not Use it
    35
    46Model props
     
    1012PARAMETERS
    1113
    12 outer   N                               as Integer      (Brief = "Number of Profile Intervals", Default = 1, Lower = 1, Upper = 100);
     14outer   N               as Integer      (Brief = "Number of Profile Intervals", Default = 1, Lower = 1, Upper = 100);
    1315outer   NComp   as Integer      (Brief = "Number of chemical components", Lower = 1);
    1416       
    1517VARIABLES
    1618
    17         Fw(N+1)                                 as flow_mass    (Brief = "Mass Flow Profile" , Symbol = "F_w");
    18         rho(N+1)                                as dens_mass    (Brief = "Mass Density Profile" , Symbol = "\rho");
    19         mu(N+1)                                 as viscosity            (Brief = "Viscosity Profile" , Symbol = "\mu");
    20 
    21         Vel(N+1)                                        as velocity     (Brief = "Velocity Profile");
    22         vm(N+1)                                 as vol_mol      (Brief = "Molar Volume Profile");
    23         frac(N+1)                                       as fraction             (Brief = "Molar Fraction Profile");
    24         Vsfrac(N+1)                             as fraction             (Brief = "No Slip Volume Fraction Profile",Lower=0);
    25         Holdup(N+1)                             as fraction             (Brief = "Holdup Profile",Lower=0);
    26         Mfrac(N+1,NComp)        as fraction             (Brief = "Phase Molar Fraction Profile");
    27         h(N+1)                                          as enth_mol     (Brief = "Molar Enthalpy profile");
    28        
     19        Fw(N+1)                         as flow_mass    (Brief = "Mass Flow Profile" , Symbol = "F_w");
     20        rho(N+1)                        as dens_mass    (Brief = "Mass Density Profile" , Symbol = "\rho");
     21        mu(N+1)                         as viscosity            (Brief = "Viscosity Profile" , Symbol = "\mu");
     22
     23        Vel(N+1)                        as velocity     (Brief = "Velocity Profile");
     24        vm(N+1)                         as vol_mol      (Brief = "Molar Volume Profile");
     25        frac(N+1)                       as fraction     (Brief = "Molar Fraction Profile");
     26        Vsfrac(N+1)                     as fraction     (Brief = "No Slip Volume Fraction Profile",Lower=0);
     27        Holdup(N+1)                     as fraction     (Brief = "Holdup Profile",Lower=0);
     28        Mfrac(N+1,NComp)        as fraction     (Brief = "Phase Molar Fraction Profile");
     29        h(N+1)                          as enth_mol (Brief = "Molar Enthalpy profile");
    2930       
    3031end
     
    9192       
    9293        Lincr(N+1)              as length                       (Brief = "Length Points", Symbol = "L_{incr}");
    93         f(N+1)                          as fricfactor           (Brief = "Friction Factor");
     94        fns(N+1)                        as fricfactor           (Brief = "No Slip Friction Factor");
     95        ftp(N+1)                        as fricfactor           (Brief = "Two phase Friction Factor");
    9496        Mw                                      as molweight    (Brief = "Average Mol Weight");
    9597        zmass(NComp)    as fraction                     (Brief = "Mass Fraction");
     
    103105        h(N+1)                          as enth_mol             (Brief = "Molar Enthalpy profile");
    104106       
    105 #       dPdLfric(N+1)           as positive             (Brief = "Friction Gradient",Lower = 0, Unit = 'Pa/m', DisplayUnit = 'kPa/m');
    106 #       dPelvdL(N+1)            as positive             (Brief = "Elevation Gradient", Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m');
    107 #       dPaccdL(N+1)            as positive             (Brief = "Acceleration Gradient",Lower = 0 , Unit = 'Pa/m', DisplayUnit = 'kPa/m');
     107        dPdLfric(N+1)           as positive             (Brief = "Friction Gradient",Lower = 0, Unit = 'Pa/m', DisplayUnit = 'kPa/m');
     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');
    108110       
    109111#Beggs&Brill Method
    110112#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    111         NFR(N+1)                                        as positive             (Brief="Froude Number", Lower=1E-10,Symbol = "N_{FR}");
     113        NFR(N+1)                                as positive             (Brief="Froude Number", Lower=1E-10,Symbol = "N_{FR}");
    112114        BBLine1(N+1)                    as positive             (Brief="Beggs and Brill Correlation Parameter 1");
    113115        BBLine2(N+1)                    as positive             (Brief="Beggs and Brill Correlation Parameter 2");
    114116        BBLine3(N+1)                    as positive             (Brief="Beggs and Brill Correlation Parameter 3");
    115117        BBLine4(N+1)                    as positive             (Brief="Beggs and Brill Correlation Parameter 4");
    116         Map(N+1)                                        as positive             (Brief="Beggs and Brill Flag ");#apenas verificação do algoritmo
     118        Map(N+1)                                as positive             (Brief="Beggs and Brill Flag ");#apenas verificação do algoritmo
    117119        Atrans(N+1)                             as positive             (Brief="Beggs and Brill Correction When Flow is in Transition Pattern");
    118         Y(N+1)                                          as positive             (Brief="Beggs and Brill Correction for Friction Factorin Two Phase Flow");
     120        Y(N+1)                                  as positive             (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow");
     121        S(N+1)                                  as Real                         (Brief="Beggs and Brill Correction for Friction Factor in Two Phase Flow");
    119122       
    120123EQUATIONS
     
    173176Y(i) = 1;
    174177
     178ftp(i)= 1;
     179       
     180        S(i)= 1;
     181       
    175182else if Liquid.frac(i) equal 0 then # Only Vapour phase
    176183
     
    209216Y(i) = 1;
    210217
     218ftp(i)= 1;
     219       
     220        S(i)= 1;
     221       
    211222else # Two phase
    212223
     
    225236"Vapour Mass Flow"
    226237        Vapour.Fw(i) =  Inlet.F*Vapour.frac(i)*sumt(M(1:NComp)*Vapour.Mfrac(i,1:NComp));
    227 
     238       
    228239switch Correlation
    229240       
     
    244255Y(i)*Liquid.Holdup(i)*Liquid.Holdup(i) = Liquid.Vsfrac(i);
    245256
     257        ftp(i)= fns(i)*exp(S(i));
     258       
     259        S(i)= ln(abs(2.2*Y(i)-1.2));
     260       
    246261# Find the Flow Pattern for Beggs-Brill Method
    247262#++++++++++++++++++++++++++++++++++++++++++
     
    324339        case "Lockhart-Martinelli":
    325340       
     341"Do Not Use Variables related to Beggs and Brill Method"
    326342        BBLine1(i) = 1;
    327343
     344"Do Not Use Variables related to Beggs and Brill Method"
    328345        BBLine2 (i)= 1;
    329346
     347"Do Not Use Variables related to Beggs and Brill Method"
    330348        BBLine3(i) = 1;
    331349
     350"Do Not Use Variables related to Beggs and Brill Method"
    332351        BBLine4 (i)= 1;
    333352
     
    337356        Map(i)=20;
    338357        Y(i)=1;
     358       
     359        ftp(i)= 1;
     360       
     361        S(i)= 1;
    339362
    340363end
     
    389412        Vel= Vapour.Vel+Liquid.Vel;
    390413
    391 "Incremental Friction Pressure Drop"
    392         #dPfric(1:N+1)  = (0.5*f(1:N+1)*Lincr(1:N+1)*rho(1:N+1)*Vel(1:N+1)^2/Dpipe);
    393         dPfric  = (0.5*f*Lincr*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel*Vel/Dpipe);
    394        
    395         #DPDLfric       = (0.5*f*(Liquid.rho*Liquid.Vsfrac+Vapour.rho*Vapour.Vsfrac)*Vel^2/Dpipe);
    396 
    397414"Incremental Acceleration Pressure Drop at Pipe Entry"
    398415        dPacc(1) = 0*'Pa';
     
    400417"Incremental Elevation Pressure Drop"
    401418        dPelv= rho*g*Lincr*sin(Hrise);
     419
     420"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);
    402423       
    403424for i in [1:N+1]
     
    450471        dPacc(i+1) = (Fw/Apipe)^2*(1/rho(i+1)-1/rho(i));
    451472        #dPacc(i+1) = dPacc(1);
    452        
     473
    453474"Incremental Length"
    454475        Lincr(i+1) = i*abs(Lpipe)/N;
     
    479500
    480501"Friction Factor for Pressure Drop - laminar Flow"
    481         f(i)*Re(i) = 16;
     502        fns(i)*Re(i) = 16;
    482503
    483504        when Re(i) > 2300 switchto "turbulent";
     
    486507
    487508"Friction Factor  for Pressure Drop - Turbulent Flow"
    488         1/sqrt(f(i))= -2*log(Roughness/Dpipe/3.7+2.51/Re(i)/sqrt(f(i)));
     509        1/sqrt(fns(i))= -2*log(Roughness/Dpipe/3.7+2.51/Re(i)/sqrt(fns(i)));
    489510
    490511        when Re(i) <= 2300 switchto "laminar";
    491512
    492513end
    493 
    494 end
    495 
    496 end
     514 
     515end
     516
     517end
     518
    497519
    498520FlowSheet Pipe_simples
     
    515537SET
    516538        NComp  = PP.NumberOfComponents;
    517         Tube.N  = 2;
     539        Tube.N  = 5;
    518540       
    519541        #Tube.Thermal ="Constant Temperature";
     
    539561
    540562        Tube.Dpipe                      = 3*'in';
    541         Tube.Roughness  = 4.572E-05*'m';
     563        Tube.Roughness  = 4.572E-5*'m';
    542564       
    543565        OPTIONS
Note: See TracChangeset for help on using the changeset viewer.