Fluid Flow and Piping: Section 17
Fluid Flow and Piping: Section 17
FIG. 17-1
Nomenclature
A = pipe cross sectional area, ft2 (A = πD2/4) NE = abscissa of Eaton correlation, Fig. 17-20
c = sum of allowances for corrosion, erosion, etc., in., NLv = liquid velocity number, from Eq 17-53
Fig. 17-23 Ngv = gas velocity number, from Eq 17-54
C = design parameter used in Hazen and Williams Nd = pipe diameter number, from Eq 17-55
formula, Eq. 17-33 NL = liquid viscosity number, from Eq 17-56
C1 = discharge factor from Fig. 17-8 P = pressure, psia
C2 = size factor from Fig. 17-9 P1 = inlet pressure, psia
d = internal diameter of pipe, in. P2 = outlet pressure, psia
do = outside pipe diameter, in. Pavg = average pressure, psia, from Eq 17-16
D = internal diameter of pipe, feet Pb = base absolute pressure, psia (ANSI 2530 specifi-
E = pipeline efficiency factor (fraction) cation: Pb = 14.73 psia)
E′ = longitudinal weld joint factor from ANSI B31.3, Pi = internal design pressure, psig
Fig. 17-23 ∆P100 = pressure drop, psi/100 ft equivalent pipe length
E′′ = longitudinal joint factor from ANSI B31.8, ∆Pe = elevation component of pressure drop, psi
Fig. 17-24 ∆Pf = frictional component of pressure drop, psi
ff = Fanning friction factor ∆Pt = total pressure drop, psi
fm = Moody friction factor (fm = 4.0 ff) q = flow rate, gal./min
fn = single phase friction factor for Dukler calculation, Q = flow rate of gas, cubic feet per day at base condi-
from Eq 17-44 tions
ftpr = friction factor ratio for Dukler calculation, QL = liquid volumetric flow rate at flowing conditions,
Fig. 17-17 ft3/sec
F′′ = construction type design factor used in Qg = gas volumetric flow rate at flowing conditions,
ANSI B31.8, Fig. 17-24 ft3/sec
Fpv = volume correction for a non-ideal fluid due to Re = Reynolds number
compressibility from Eq 17-13 Rey = mixture Reynolds number for Dukler calculation,
√1/ff = transmission factor from Eq 17-45
g = acceleration due to gravity, 32.2 ft/sec2 S = specific gravity of flowing gas (air = 1.0)
gc = gravitational constant, 32.2 (ft • lbm)/(lbf • sec2) S′ = allowable stress, psi, Fig. 17-23
hL = loss of static pressure head due to fluid flow, feet S′′ = specified minimum yield strength, psi, Fig. 17-24
of fluid t = thickness, in., Figs. 17-23, 17-24
H = total energy of a fluid at a point above a datum, tm = minimum required wall thickness, in., Fig. 17-23
from Eq 17-1 T = absolute temperature of flowing gas, °R
HLd = liquid holdup fraction (Dukler), Fig. 17-18 T′′ = temperature derating factor used in ANSI B31.8,
HLe = liquid holdup fraction (Eaton), Fig. 17-20 Fig. 17-24
HLf = liquid holdup fraction (Flanigan), Fig. 17-19 Tavg = average temperature, °R, [Tavg = 1/2 (Tin + Tout)]
IL = liquid inventory in pipe, ft3, from Eq 17-57 Tb = base absolute temperature, °R (ANSI 2530 speci-
L = length of line, feet fication: Tb = 520°R)
Lm = length of line, miles V = single phase fluid velocity, ft/sec
MW = molecular weight Vsg = superficial gas velocity, ft/sec, from Eq 17-36
Nx = Fig. 17-16 horizontal coordinate, ft/sec VsL = superficial liquid velocity, ft/sec, from Eq 17-35
Ny = Fig. 17-16 vertical coordinate, ft/sec Vm = mixture velocity, ft/sec, from Eq 17-46
17-1
FIG. 17-1 (Cont’d)
Nomenclature
17-2
Pressure Loss Due to Friction overall pressure drop is the sum of the ∆ Pf values calculated
for the individual segments. For gas applications the segmen-
Flow is always accompanied by friction. This friction results tal length may be relatively short, as compared to liquid ap-
in a loss of energy available for work. A general equation for plications, since many gas applications involve compressible
pressure drop due to friction is the Darcy-Weisbach2 (often gases where gas densities vary with pressure.
referred to as simply the Darcy) equation. This equation can
be rationally derived by dimensional analysis, with the excep- Friction Factor and Effect of Pipe Roughness
tion of the friction factor, fm, which must be determined ex-
perimentally. Expressed in feet of fluid this equation is: When the fluid flow is laminar (Re<2000), the friction factor
has a direct relationship to the Reynolds number, such that:
2
fm L V
hL = Eq 17-6 fm = 64 /Re or ff = 16 /Re Eq 17-8
2gD
Pipe roughness has no effect on the friction factor in laminar
Converting to pounds per square inch, the equation be- flow.
comes:
Substitution of the formula for Reynolds number, Eq 17-4,
ρ fm L V2
∆Pf = Eq 17-7 into Eq 17-8, yields the following:
(144 ) D (2gc)
64 µe 64 µ 12
It should be noted that the Moody friction factor3, fm, is used fm = = Eq 17-9
in the equations above. Some equations are shown in terms of
DVρ Vρ 1488 d
the Fanning friction factor, ff, which is one fourth of fm (fm = This expression can then be substituted for the friction fac-
4.0 ff). A graph of both Fanning and Moody friction factors as tor in Eq 17-7, resulting in the following formula for pressure
a function of Reynolds number appears in Fig. 17-2. loss in pounds per square inch:
The Darcy-Weisbach equation is valid for both laminar and µLV
∆Pf = 0.000668 Eq 17-10
turbulent flow of any liquid, and may also be used for gases d2
with certain restrictions. When using this equation, changes Eq 17-10 is commonly known as Poiseuille’s law for laminar
in elevation, velocity, or density must be accounted for by ap- flow.
plying Bernoulli’s theorem. The Darcy-Weisbach equation
must be applied to line segments sufficiently short such that When the flow is turbulent, the friction factor depends on
fluid density is essentially constant over that segment. The the Reynolds number and the relative roughness of the pipe,
FIG. 17-2
Friction Factors5
For more accurate values of Z, refer to Section 23. For more The average pressure in the line can be computed by:
approximate calculations, the value of the average compressi- 2 P1 P2
Pavg = P1 + P2 −
3 P1 + P2
bility factor, Zavg, may be calculated from the following equa- Eq 17-16
tions:
In the absence of field data indicating otherwise, an effi-
1
Zavg = Eq 17-12 ciency factor, E, of 1.0 is usually assumed.
(Fpv)2
The AGA Equations — The AGA Equations were devel-
and oped to approximate partially and fully turbulent flow using
0.5 two different transmission factors. The fully turbulent flow
(Pavg ) (3.444 ) (105) (10(1.785) (S)) equation accounts for the relative pipe roughness, ε/D, based
Fpv = 1 + Eq 17-13
T3.825
avg on the rough-pipe law.4 This equation uses the following trans-
mission factor:
Fig. 17-5 contains a plot of the deviation factor, Fpv2.
3.7 D
An estimate for Zavg at pressures below 100 psi is: 1/ff = 4 log10
√ Eq 17-17
ε
1
Zavg = Eq 17-14 When the transmission factor for fully turbulent flow is sub-
1 + 0.0002 Pavg stituted in the general energy equation (Eq 17-15), the AGA
Equation for fully turbulent flow becomes:
0.5
SINGLE PHASE FLOW Tb 3.7 D P21 − P22 2.5
Q = 38.77 E 4 log10 d
Pb ε S Lm Tavg Zavg
Transmission Line Gas Flow
Eq 17-18
Isothermal Flow — The steady-state, isothermal flow be-
The partially turbulent flow equation is based on the
havior of gas in pipelines is defined by a general energy equa-
smooth-pipe law4 and is modified to account for drag-inducing
tion of the form:
elements. The transmission factor for this equation is:
0.5
Tb 1 P21 − P22 2.5
Q = 38.77 E
P f
√
S L T Z
d Eq 17-15 √
1/ff = 4 log10
Re
1/ff
√
− 0.6 Eq 17-19
b f m avg avg
17-4
FIG. 17-3
Relative Roughness of Pipe Materials and Friction Factors for Complete Turbulence5
17-5
FIG. 17-4
Equivalent Length of Valves and Fittings in Feet
Short Long
Nominal Pipe size in.
ell T T
Globe valve or
ell ell
Angle valve
Plug cock
Sudden Std. red. Sudden Std. red.
Threaded
Threaded
Threaded
Threaded
Threaded
2 miter
3 miter
4 miter
Welded
Welded
Welded
Welded
Welded
Equiv. L in terms of small d
d/D = 1/4
d/D = 1/2
d/D = 3/4
d/D = 1/2
d/D = 3/4
d/D = 1/4
d/D = 1/2
d/D = 3/4
d/D = 1/2
d/D = 3/4
11⁄2 55 26 13 7 1 1 2 3 5 2 3 8 9 2 3 5 3 1 4 1 3 2 1 1 –
2 70 33 17 14 2 2 3 4 5 3 4 10 11 3 4 7 4 1 5 1 3 3 1 1 –
21⁄2 80 40 20 11 2 2 – 5 – 3 – 12 – 3 – 8 5 2 6 2 4 3 2 2 –
3 100 50 25 17 2 2 6 4 14 4 10 6 2 8 2 5 4 2 2 –
4 130 65 32 30 3 3 7 5 19 5 12 8 3 10 3 6 5 3 3 –
6 200 100 48 70 4 4 11 8 28 8 18 12 4 14 4 9 7 4 4 1
8 260 125 64 120 6 6 15 9 37 9 25 16 5 19 5 12 9 5 5 2
10 330 160 80 170 7 7 18 12 47 12 31 20 7 24 7 15 12 6 6 2
12 400 190 95 170 9 9 22 14 55 14 28 21 20 37 24 8 28 8 18 14 7 7 2
14 450 210 105 80 10 10 26 16 62 16 32 24 22 42 26 9 – – 20 16 8 – –
16 500 240 120 145 11 11 29 18 72 18 38 27 24 47 30 10 – – 24 18 9 – –
18 550 280 140 160 12 12 33 20 82 20 42 30 28 53 35 11 – – 26 20 10 – –
20 650 300 155 210 14 14 36 23 90 23 46 33 32 60 38 13 – – 30 23 11 – –
22 688 335 170 225 15 15 40 25 100 25 52 36 34 65 42 14 – – 32 25 12 – –
24 750 370 185 254 16 16 44 27 110 27 56 39 36 70 46 15 – – 35 27 13 – –
30 – – – 312 21 21 55 40 140 40 70 51 44
36 – – – 25 25 66 47 170 47 84 60 52
42 – – – 30 30 77 55 200 55 98 69 64
48 – – – 35 35 88 65 220 65 112 81 72
54 – – – 40 40 99 70 250 70 126 90 80
60 – – – 45 45 110 80 260 80 190 99 92
Substituting √
1/ff from Eq 17-19 into Eq 17-15 does not pro- correction factors determined from the system to which it is to
vide an equation which can be solved directly. For partially be applied. Graphs showing gas flow calculations based on the
turbulent flow a frictional drag factor must also be applied to Weymouth equation are shown in Fig. 17-6a and 17-6b.
account for the effects of pipe bends and irregularities. These Panhandle A Equation — In the early 1940s Panhandle
calculations are beyond the scope of this book and the AGA Eastern Pipe Line Company developed a formula for calcula-
“Steady Flow in Gas Pipelines”6 should be consulted for a de- tion of gas flow in transmission lines which has become known
tailed treatment of partially turbulent flow. as the Panhandle A Equation. This equation uses the following
The Weymouth Equation — The Weymouth Equation, expressions of Reynolds number and transmission factor.
published in 19127, evaluated the coefficient of friction as a QS
function of the diameter. Re = 1.934 Eq 17-23
d
0.008
ff = Eq 17-20 0.07305
d1/3 QS
√1/ff = 7.211
= 6.872 (Re)0.07305 Eq 17-24
d
1/ff = 11.18 d1/6
√ Eq 17-21
The transmission factor assumes a Reynolds number value
When the friction factor, ff, is substituted in the general en- from 5 to 11 million based on actual metered experience.
ergy equation, Weymouth’s Equation becomes:
0.5 Substituting Eq 17-24 for √
1/ff in the general energy equa-
Tb P21 − P22 2.667 tion (Eq 17-15), the Panhandle A Equation becomes:
Q = (433.5) E d Eq 17-22
Pb S Lm Tavg Zavg 1.0788
0.5392
Tb P21 − P22
The Weymouth formula for short pipelines and gathering Q = 435.87 E 0.853 d2.6182
systems agrees more closely with metered rates than those Pb S Lm Tavg Zavg
calculated by most other formulae. However, the degree of er-
ror increases with pressure. If the Q calculated from the Wey- Eq 17-25
mouth formula is multiplied by √
, where Z is the
1/Z This equation was intended to reflect the flow of gas through
compressibility factor of the gas, the corrected Q will closely smooth pipes. When “adjusted” with an efficiency factor, E, of
approximate the metered flow. Fig. 17-5 shows a plot of the about 0.90, the equation is a reasonable approximation of the
deviation factor, √
, of a common gas and can be used safely
1/Z partially turbulent flow equation. The equation becomes less
if exact data is not available. accurate as flow rate increases. Many users of the Panhandle
A Equation assume an efficiency factor of 0.92.
The equation cannot be generally applied to any variety of
diameters and roughness, and in the flow region of partially Panhandle B Equation — A new or revised Panhandle
developed turbulence, it is not valid. The Weymouth Equation Equation was published in 1956. This revised equation is
may be used to approximate fully turbulent flow by applying known as the Panhandle B Equation and is only slightly
17-6
compensation for discrepancies through the use of adjustment
FIG. 17-5
factors, usually termed “efficiencies.” These efficiencies are
Deviation Factors8 frequently found in practice by determining the constant re-
quired to cause predicted gas equation behavior to agree with
flow data. As a result, the values of these factors are specific
to particular gas flow equations and field conditions and, un-
der many circumstances, vary with flow rate in a fashion that
obscures the real nature of flow behavior in the pipe.
The Reynolds number dependent equations, such as the
Panhandle equations, utilize a friction factor expression which
yields an approximation to partially turbulent flow behavior
in the case of the Panhandle A equation, and an approximation
to fully turbulent behavior in the case of the Panhandle B.
These equations suffer from the substitution of a fixed gas
viscosity value into the Reynolds number expression, which,
in turn, substituted into the flow equation, results in an ex-
pression with a preconditioned bias.
Regardless of the merits of various gas flow equations, past
practices may dictate the use of a particular equation to main-
tain continuity of comparative capacities through application
of consistent operating policy. A summary of comparisons be-
tween transmission factors used in the above gas equations
are shown in Fig. 17-7.
Reference should be made to “Steady Flow in Gas Pipe-
lines”6, published by American Gas Association, for a complete
analysis of steady flow in gas pipelines.
Low Pressure Gas Flow
Gas gathering often involves operating pressures below 100
psi. Some systems flow under vacuum conditions. For these
low pressure conditions, equations have been developed that
give a better fit than the Weymouth or Panhandle equations.
Two such formulas are:
The Oliphant Formula9 for gas flow between vacuum and
100 psi:
Note: Refer to Section 23 for more accurate compressibility factors. d3 14.4 Tb
Q = 42 (24) d2.5 +
Reynolds number dependent. Therefore, it more nearly ap- 30 Pb 520
proximates fully turbulent flow behavior. The transmission 1⁄ 2
factor used here is: 0.6 520 P21 − P22
Eq 17-28
S T Lm
0.01961
QS
1/ff = 16.70
√ = 16.49 (Re)0.01961 Eq 17-26 The Spitzglass Formula for gas flow below 1 psig at 60°F:
d 1⁄ 2
Substituting Eq 17-23 for √
1/ff in the general energy equa- 27.69 (P1 − P2) d5
tion (Eq 17-15), the Panhandle B Equation becomes: Q = (24) (3550) Eq 17-29
3.6
1.02 0.51 SL 1 + + 0.03 d
Tb P21 − P22 d
Q = 737 E 0.961 d
2.53
Eq 17-27
P
b S L m Tavg Zavg Plant Piping Gas Flow
The equation can be adjusted through the use of an effi- For estimating pressure drop in short runs of gas piping,
ciency term that makes it applicable across a relatively limited such as within plant or battery limits, a simplified formula for
range of Reynolds numbers. Other than this, however, there compressible fluids is accurate for fully turbulent flow, assum-
are no means for adjustment of the equation to correct it for ing the pressure drop through the line is not a significant frac-
variations in pipe surface. Adjusted to an average flowing tion of the total pressure (i.e. no more than 10%).
Reynolds number, the equation will predict low flow rates at
low Reynolds numbers, and high flow rates at high Reynolds The following method is a simplification of the Darcy for-
numbers, as compared to a fully turbulent flow equation. Ef- mula, which eliminates calculation of fm, the Moody friction
ficiencies based on the Panhandle B equation decrease with factor. This simplification was checked over a wide range of
increasing flow rate for fully turbulent flow. The efficiency fac- flows and densities for pressure drops of 0.25 to 1.5 psi/100 ft.
tor, E, used in the Panhandle B equation generally varies be- Density was varied over a range of 100 to 1; flows varied over
tween about 0.88 and 0.94. a range of 75 to 1. Pressure variation was from atmospheric to
1000 psia. The error from using the simplified approach as
Conclusions — The successful application of these trans- compared to the actual friction factor calculated in the Darcy
mission line flow equations in the past has largely involved formula was from zero to 5%, with the simplified approach
17-7
FIG. 17-6a
Gas Flow Based On Weymouth Formula
17-8
FIG. 17-6b
Gas Flow Based On Weymouth Formula
17-9
FIG. 17-7
Comparison of Gas Equation Transmission Factors for Nominal 20 Inch Pipe
giving consistently lower calculated pressure drop for a given perature is 60°F and pressure is 750 psia. The compressibility
flow. factor is 0.905 (from Fig. 23-4).
The Darcy formula can be written in the simple form: Solution Steps
17-10
FIG. 17-8 Water — A graph showing pressure drop for water per 100
5 feet as a function of flow rate in gallons per minute and pipe
Simplified Flow Formula for Compressible Fluids
Values of C1 size is shown in Fig. 17-10. These data are based on the Hazen
and Williams empirical formula10 using a “C” constant of 100
which is commonly used for design purposes in welded and
seamless steel pipe.
Hazen and Williams formula for flow of water:
0.54
P1 − P2
q = 0.442 d 2.63
C Eq 17-33
L
Where:
C = 140 for new steel pipe
C = 130 for new cast iron pipe
C = 100 is often used for design purposes to account for
pipe fouling, etc.
Hydrocarbon — A graph showing pressure drop for hy-
drocarbons per 100 feet as a function of flow rate in gallons
per minute and pipe size is shown in Fig. 17-11. This graph
assumes a specific gravity of 1.0 (water). To correct for differ-
ent liquid densities, the value read from Fig. 17-11 must be
multiplied by the actual specific gravity to obtain the correct
pressure loss.
Steam Flow
Fig. 17-12 contains a graphical representation of Fritzsche’s
formula11 for calculating pressure drop in steam lines.
Fritzsche’s formula and instructions for the chart usage are
given in Fig. 17-12.
The Babcock formula5 for steam flow is:
d + 3.6 W2 L
∆Pf = 3.63 (10−8) 6 Eq 17-34
d ρ
Fire Stream Flow
Fig. 17-13 is a table permitting rapid computation of the
behavior of various sized fire nozzles. The table also includes
the estimated pressure drop in 100 feet of 21⁄2 inch diameter
fire hose.
17-11
FIG. 17-9
Simplified Flow Formula for Compressible Fluids5
Values of C2
17-12
FIG. 17-10
Pressure Drop for Flowing Water
17-13
FIG. 17-11
Pressure Drop for Hydrocarbon Liquids in Smooth Pipe
Vsg = superficial gas velocity cause dispersed flow. A flow regime map generated by Taitel
and Dukler13 contains explicit inclination effects and should
= Qg/A Eq 17-36 be used for inclined pipes. This latter method also contains
Mandhane proposed a fluid property correction to the super- explicit pipe diameter effects but lacks any way of accounting
ficial velocities, but concluded that the fluid property effects for surface tension.
are insignificant compared to the errors in the empirical map. For vertical flow, the stratified flow regime cannot exist as
The map reports the flow regimes: stratified, wavy, annular there is no preferred direction for the liquid to settle. An em-
mist, bubble, slug, and dispersed. Care should be taken in the pirical flow regime map developed by Aziz14 for vertical up-
interpretation of these flow maps as the regime boundaries ward flow is shown in Fig. 17-16. The coordinates for this flow
are strongly affected by pipe inclination. In particular, hori- map are the same as for the Mandhane map in Fig. 17-15 ex-
zontal flow regime maps must not be used for vertical flow, and cept that fluid property corrections are used. The coordinates
vertical flow regime maps must not be used for horizontal flow. used in the Aziz vertical map are:
The Mandhane map given in Fig. 17-15 was developed for hori-
zontal lines flowing air and water at near atmospheric pres- Nx = Vsg XA Eq 17-37
sure. Inclinations in the range of 0.1-1.0 degrees can cause
substantial regime boundary movement. In addition, flow re- Ny = VsL YA Eq 17-38
gime boundary adjustment has been observed due to fluid 0.333
pressure, pipe diameter, and surface tension.13,22 The gas den- ρg
XA = YA Eq 17-39
sity increase caused by high pressure acts to move the slug- ρa
mist boundary to lower superficial gas velocities, while
0.25
increased pipe diameter acts to increase the stratified wavy ρLσwa
flow regime at the expense of the slug flow regime. In addition, YA = Eq 17-40
foamy fluids having a high surface tension have been observed ρwσ
to flow in the dispersed flow regime even though Mandhane For pipe inclinations greater than 10-20 degrees, flow re-
would have predicted superficial liquid velocities too low to gime patterns resemble those of vertical flow more than those
17-14
FIG. 17-12
Pressure Drop in Steam Pipes by Fritzsche’s Formula11
17-15
FIG. 17-13
Table of Effective Fire Streams
Using 100 ft of 21⁄2 in. ordinary best quality rubber lined hose between nozzle and hydrant or pump.
The vertical and horizontal distance in above table are "effective" streams. The "effective stream" is one that has not broken up into a spray
and which will project three fourths of the water through a circle 10 in. in diameter, at the distance indicated.
The bulk of the stream and extreme drops will carry a greater distance, but the stream is spread out too much to penetrate a hot fire and
reach the burning materials before evaporation. The vertical height "h" in feet, which the bulk of the stream will carry, may be computed
V2
by formula h = 0.91
2g
0.321 × gpm
The gallons/minute discharge in the above table check approximately with the formula: V = C√
2 gH =
A
VA
then gpm = V = Velocity in feet per second g = Acceleration due to gravity, = 32.16 ft/sec2
0.321
A = Area of nozzle, in square inches C = Coefficient for smooth nozzles = 0.98 H = Head in feet = pressure x 2.31
Reference: Fire Protection Handbook, 14th Edition, 1976, National Fire Protection Association
Slug
Bubble
Annular Mist
Stratified Wave
17-16
FIG. 17-16 Nx = (1.43) (3.91) = 5.60 ft/sec using Eq 17-37
Aziz et al. Vertical Up-Flow Regime Map14 Ny = (0.49) (1.32) = 0.64 ft/sec using Eq 17-38
Fig. 17-16 shows that this flow is in the slug flow regime.
Pressure Drop Calculation
Calculation of pressure drop in two-phase flow lends itself
better to computer calculation than to hand calculation. Sev-
eral two-phase pressure drop correlations are available for
both horizontal and vertical flows.15,16,17 Due to the complexity
of two-phase flow, uncertainties associated with pressure drop
calculations are much greater than uncertainties in single-
phase pressure drop calculations. As a result, errors in calcu-
lated two-phase pressure drops in the order of plus or minus
twenty percent may normally be anticipated, especially in cir-
cumstances where fluid velocities are unusually high or low,
where terrain is rugged, or where fluid properties are inade-
quately known. In addition, different two-phase flow correla-
tions may give significantly different pressure drops. In order
to evaluate these differences, several correlations should be
used. A method suggested by the American Gas Association18
can serve as a basis for hand calculation generated by Dukler19
and an elevation pressure drop correlation by Flanigan.20
Frictional Component18 — Using the Dukler frictional
pressure drop calculation method, the frictional pressure drop
is given by the equation:
fn ftpr ρk V2m Lm
∆Pf = Eq 17-41
of horizontal flow, and the Aziz vertical map should be used. (0.14623)d
For vertical downward flow, the flow patterns can be more com- where
plicated. A generally accepted vertical down-flow map is not ρL λ2 ρg (1 − λ)2
available. ρk = + Eq 17-42
HLd (1 − HLd)
Although the designer frequently does not have the choice,
and
avoidance of the slug flow regime in horizontal flow and the
slug and froth flow regimes in vertical flow is desirable. Slug QL
λ = Eq 17-43
flow introduces a flow rate and pressure intermittency that QL + Qg
may be troublesome to process control. In some cases, slug flow
The single phase friction factor, fn, can be obtained from the
may be avoided by the choice of smaller pipe sizes. Of course,
correlation:19
frictional pressure drop may be increased by use of this
smaller pipe. In vertical upflow, oversizing the pipe may result fn = 0.0056 + 0.5 (Rey)−0.32 Eq 17-44
in a bubble flow regime, and a large liquid inventory. This liq- The mixture Reynolds number, Rey, is calculated according
uid inventory may cause excessive hydrostatic presure drops. to the equation:
Example 17-3 — A vapor-liquid mixture is flowing verti- (124.0) ρk Vm d
cally upward in a pipe having an inside diameter of 8.0 inches. Rey = Eq 17-45
µn
The fluid is a hydrocarbon liquid-hydrocarbon vapor mixture.
The liquid density is 52 lb/ft3 and the vapor density is 2.0 lb/ Calculation of this Reynolds number requires determina-
ft3. The interfacial surface tension is 20 dynes/cm. The liquid tion of mixture velocity, Vm, and mixture viscosity, µn. These
volumetric flow rate is 0.17 ft3/sec and the vapor flow rate is quantities can be determined according to:
0.5 ft3/sec measured at acutal conditions. What flow regime is Vm = VsL + Vsg Eq 17-46
to be expected?
µn = µL λ + µg (1 − λ) Eq 17-47
Solution Steps
The two-phase friction factor ratio, ftpr, representing a two-
0.25 phase frictional “efficiency” can be determined by reference to
(52.0) (72.4)
YA = = 1.32 using Eq 17-40 Fig. 17-17 or by the equation:
(62.4 (20) y
ftpr = 1 + 4
(2.0)
0.333
1.281 − 0.478y + 0.444y 2
− 0.094y3
+ 0.00843y
XA = (1.32) = 3.91 using Eq 17-39 Eq 17-48
0.0764
where y = −ln(λ).
0.5
Vsg = = 1.43 ft/sec using Eq 17-36 The remaining quantity to be calculated in the Dukler
(π/4) (8/12)2
scheme is an estimate of the liquid holdup, HLd. This holdup
can be estimated using Fig. 17-18. This figure gives liquid hold-
0.17
VsL = = 0.49 ft/sec using Eq 17-35 up as a function of λ and Rey. Since Rey is itself a function of
(π/4) (8/12)2
17-17
liquid holdup, the calculation is, in general, iterative. For most ventory in the pipe, IL, the pipe internal volume is multiplied
calculations, however, the Rey line can be used for a first esti- by this holdup fraction.
mate.
IL = (28.80) HLed2Lm Eq 17-57
Elevation Component — The elevation component of Since holdup fractions may change along the length of the
pressure drop can be found using the Flanigan method. In this pipe, a segmental calculation is more accurate.
method, the elevation component is calculated using the equa-
tion: Example 17-4 — A pipeline segment with a 6-inch inside
ρL HLf diameter, 0.75 miles long, transports a mixture of gas and oil.
∆Pe = Σ Ze Eq 17-49 The pipeline has a gradual upward slope and rises 100 feet
144 over the 0.75 mile length. The inlet pressure of the pipeline is
where HLf is determined from Fig. 17-19 or calculated accord- 400 psia, liquid viscosity is 20 cp, the vapor viscosity is 0.015
ing to the formula: cp, and the interfacial surface tension is 15 dynes/cm. The liq-
uid flow rate is 10 ft3/min and the vapor flow rate is 250 actual
1 ft3/min. The density of the liquid phase is 55 lb/ft3, and the
HLf = Eq 17-50
1 + 0.3264 (Vsg)1.006 density of the gas phase is 1.3 lb/ft3 at operating conditions.
What is the pressure at the downstream end of the line seg-
The term Ze is the vertical elevation rise of a hill. The rises ment, and what is the liquid inventory of the line?
are summed. No elevation drops are considered. This is tanta-
mount to ignoring any possible hydrostatic pressure recover- Solution Steps
ies in downhill sections of pipeline and may lead to a Calculate the flowing liquid volume fraction using Eq 17-43.
considerable error in the pressure drop analysis. 10
λ = = 0.038
Once the frictional component or pressure drop is found us- 250 + 10
ing the Dukler method, and the elevation component is found
Calculate the mixture viscosity, µn using Eq 17-47.
using the Flanigan method, the overall two-phase pressure
drop is found by summing the friction and elevation compo- µn = (20) (0.038) + (0.015) (1.0 − 0.038) = 0.774 cp
nents. For a first guess, assume HLd = λ and estimate ρk using
∆Pt = ∆Pe + ∆Pf Eq 17-51 Eq 17-42.
Since fluid properties and liquid holdups can change rapidly (55) (0.038)2 (1.3) (1.0 − 0.038)2
in a two-phase line, accuracy is improved if this AGA calcula- ρk = + = 3.341 lb/ft3
0.038 (1.0 − 0.038)
tion procedure is performed segmentally. The need for segmen-
tal calculations is one of the reasons why two-phase Calculate the superficial velocities and the mixture velocity.
calculations are best suited for computer calculation. 10.0
VsL = = 0.849 ft/sec using Eq 17-35
Liquid Holdup — The liquid holdup correlation given in (π/4) (6/12)2 (60)
Fig. 17-18 is intended only for use in the Dukler friction pres-
sure drop calculation. A correlation by Eaton et al.21 is better 250
Vsg = = 21.22 ft/sec using Eq 17-36
suited for liquid holdup determination in liquid inventory cal- (π/4) (6/12)2 (60)
culations.
Vm = 0.849 + 21.22 = 22.07 ft/sec using Eq 17-46
The Eaton21 holdup correlation is shown in Fig. 17-20. In
this figure, the holdup fraction, HLe, is plotted directly as a Calculate an estimate of the mixture Reynolds number, Rey,
function of the dimensionless group, NE. This dimensionless using Eq 17-45
group is of the form: (124.0) (3.341) (22.07) (6.0)
Rey = = 70,878
Pavg
0.05
(0.774)
1.84 (NLv)0.575 (NL)
0.1
ρL
0.25 Using this improved HLd, recalculate ρk using Eq 17-42
NLv = 1.938 VsL Eq 17-53
(55) (0.038)2 (1.3) (1.0 − 0.038)2
σ ρk = + = 2.029 lb/ft3
0.25 (0.12) 1.0 − 0.12
ρL
Ngv = 1.938Vsg Eq 17-54 Using this improved ρk recalculate Rey using Eq 17-45
σ (124.0 ) (2.029 ) (22.07) (6.0)
0.50 Rey = = 43,044
ρL (0.774)
Nd = 10.073 d Eq 17-55
σ From Fig. 17-18 with λ = 0.038 and Rey = 43,044, HLd = 0.16.
0.25
Another iteration using HLd = 0.16 indicates Rey = 40,923 and
1 HLd = 0.16.
NL = 0.15726 µL 3 Eq 17-56
ρLσ Calculate the single phase friction factor with Rey = 40,923,
The Eaton correlation has been found reasonably accurate using Eq 17-44
by several investigators, particularly for low holdup flows. fn = 0.0056 + 0.5 (40,923)-0.32 = 0.0223
The liquid holdup fraction, HLe, is the fraction of the flow Determine the two-phase friction factor, ftpr, from Fig. 17-17
area of the pipe occupied by liquid. To calculate the liquid in- using λ= 0.038
17-18
FIG. 17-17
Two-Phase Friction Factor Ratio19
FIG. 17-18
Liquid Holdup Correlation19
17-19
FIG. 17-19
Flanigan Liquid Holdup Correlation20
17-20
FIG. 17-21 FIG. 17-22
Multiple Pipe Slug Catcher Example Line Drip
GRADE
GAS FLOW
DRIP
BELOW
LINE
DRIP LIQUID
DRIP
VESSEL
DRIP
VALVE NO
HLe = 0.14 inventory in the pipeline decreases, and the excess liquid
Note that this estimate is close to the HLf predicted in forms a slug or series of slugs. Finally, pigging can cause very
Fig. 17-19 for elevation pressure drop determination. It also large liquid slugs as the entire liquid inventory of the line is
coincides closely with the value of 0.16 from Fig. 17-18. swept ahead of the pig. Of the four mechanisms described,
wave growth normally produces the shortest slugs, followed in
Calculate the pipeline segment liquid inventory from Eq 17- 57 length by terrain generated slugs. Methods for calculating
IL = (28.80) (0.14) (6.0)2 (0.75) = 108.86 ft3 wave induced slugs were described by Greskovich and
Shrier22, and by Brill et al.23 A preliminary scheme for calcu-
The pipeline segment contains 108.86 cubic feet of liquid at lating terrain generated slugs was reported by Schmidt.24
any instant. Analytical methods for determining inlet flow rate generated
Liquid Slugging slugs were given by Cunliffe,25 and a method of analyzing pig-
ging dynamics was given by McDonald and Baker.26
Purpose of Separators — The slug flow regime is fre-
quently encountered for pipe sizes and flow rates used in proc- Slug Catchers — Slug catchers are devices at the down-
ess and transmission piping. Liquid slugging introduces an stream end or other intermediate points of a pipeline to absorb
additional design and operational difficulty as liquid and va- the fluctuating liquid inlet flow rates through liquid level fluc-
por must generally be separated at the downstream end of the tuation. Slug catchers may be either a vessel or constructed of
two-phase flow line. The downstream separator serves both as pipe. All size specifications discussed in Section 7 to provide
a liquid-vapor disengaging device and as a surge vessel to ab- residence time for vapor-liquid disengagement also apply to
sorb the fluctuating liquid flow rates caused by slugging. In vessels used as slug catchers. In addition, sufficient volume
order to size the separator or slug catcher, the length of the must be provided for liquid level fluctuation. Particularly for
incoming slugs must be determined. Slug length calculation high pressure service, vessel separators may require very
methods are not well developed, and there is large uncertainty thick walls.
in slug length determination.
In order to avoid thick wall vessels, slug catchers are fre-
Mechanisms of Slug Generation — Liquid slug lengths quently made of pipe. Lengths of line pipe tens or hundreds of
are difficult to determine in part because there are at least feet long are used as long, slender horizontal separators. The
four identifiable mechanisms for liquid slug generation. Slugs pipe is generally inclined from one to ten degrees and banks
can form as the result of wave formation at the liquid-gas in- of these slightly inclined pipes are frequently manifolded to-
terface in a stratified flow. When the liquid waves grow large gether. Pipe type slug catchers are frequently less expensive
enough to bridge the entire pipe diameter, the stratified flow than vessel type slug catchers of the same capacity due to thin-
pattern breaks down and a slug flow is formed. ner wall requirements of smaller diameter pipe. The manifold
Slugs can also form due to terrain effects. Liquid collects at nature of multiple pipe slug catchers also makes possible the
a sag in the pipeline and blocks the gas flow. The pressure in later addition of additional capacity by laying more parallel
this blocked gas rises until it blows the accumulated liquid in pipes. A schematic of a multiple pipe (harp) slug catcher ap-
the sag out as a slug. Changes in pipeline inlet flow rate can pears in Fig. 17-21. Different pipe inclinations and different
also cause slugs. When the inlet flow rate increases, the liquid manifolding arrangements are favored by different designers.
17-21
An example of a line drip catcher is shown in Fig. 17-22. A pig may become stuck in the line and require an expensive
drip vessel is connected to the incoming pipeline and often laid shutdown for location and removal.
beneath it. A flow line from the drip vessel is used to blow the
liquids out to a storage or surge vessel as they accumulate.
Pigging — Pipelines are pigged for several reasons. If PIPE AND FLANGE DATA
water is present in the line, it must be removed periodically in
order to minimize corrosion. This water accumulates in sags
in the pipeline, and these low spots are particularly suscepti- The Petroleum Refinery Piping Code (ANSI B31.3) is used
ble to corrosion. Pipelines are also pigged to improve pressure to determine the allowable pressure limits for piping inside
drop-flow rate performance. Water or hydrocarbon liquids that refineries and other processing facilities. Refer to Fig. 17-23
settle in sags in the pipeline constitute partial blockages that for the calculation method.
increase pressure drop. Pigging can remove these liquids and
improve pipeline efficiency. Pigging can also be used as a FIG. 17-24
means of limiting the required slug catcher size. By pigging at
frequent intervals, liquid inventory buildup in a pipeline can Working Pressures Transmission Lines
be reduced, and the maximum slug size can be limited. The
required downstream slug catcher size must take into account
To determine allowable internal working pressures for pip-
pigging frequency.
ing outside of refineries and other processing facilities in
Operational hazards are associated with pigging. The very accordance with ANSI B31.8-1999, “Code of Pressure Pip-
large slugs swept ahead of the pig may overwhelm inade- ing, Gas Transmission and Distribution Piping,” use the
quately sized downstream facilities. Pigs may also occasion- following:
ally be destroyed in the pipeline and the resulting debris may 2 S′′ t
damage downstream fittings or equipment. Even worse, the Pi = (F′′) (E′′) (T′′)
do
FIG. 17-23 Where
Working Pressures Refinery Piping Pi = Design pressure, psig
S′′ = Specified minimum yield strength, psi
To determine allowable internal working pressure for do = Nominal outside diameter, in.
straight sections of pipe in accordance with ANSI B31.3, t = Nominal wall thickness, in.
“Code for Pressure Piping, Petroleum Refinery Piping”, use F′′ = Construction type design factor,
the following: Table 841.114A and Par. 840.2 (see note)
tm = t + c Location Class F′′
Pido 2 tS′E′ 1 Div 1 .80
t = or Pi =
2(S′E′ + PiY′) do − 2 tY′ Div 2 .72
2 .60
do Pi
For t > or > 0.385 calculation of pressure design 3 .50
6 S′ E′
4 .40
thickness requires special consideration.
where: Complete details are covered in Par. 841.
tm = minimum required thickness, satisfying E′′ = Longitudinal joint factor, Table 841.115A
requirements for pressure, and mechanical, Normally a factor of 1.0 is used for seamless
corrosion, and erosion allowances, inches. and welded pipe except for the following:
Selected thickness shall be no less than tm Fusion Welded A 134 and A 139 0.80
plus the manufacturer’s negative tolerance. Spiral Welded A 211 0.80
The typical negative tolerance is 12.5%.
Furnace Butt Welded ASTM-A53, API-5L 0.60
t = pressure design thickness, in.
c = the sum of the mechanical allowances T′′ = Temperature derating factor, Table 841.116A
(thread depth and groove depth), corrosion, Temp, °F Factor T′′
and erosion allowances, in.
250 or less 1.000
Pi = internal design pressure, psig.
300 0.967
do = outside diameter of pipe, in.
350 0.933
S′ = allowable stresses, psi.
E′ = longitudinal weld joint factor: 400 0.900
Seamless = 1.000, ERW = 0.85 450 0.867
Y′ = coefficient having values for ferritic steels
For intermediate temperatures, interpolate for
as follows: 0.4 up to and including 900°F
derating factor.
0.5 for 950°F
Note: Factor reflecting location of line, proximity to roads,
0.7 for 1,000°F and above
public or private land, etc.
17-22
FIG. 17-25
Representative Allowable Stresses in Tension for Materials
(Developed from ASME B31.3, 2002, Tables A-1, A-1B)
17-23
FIG. 17-26
Design Properties and Allowable Working Pressures for Piping
ASTM A106, grade B seamless pipe—Petroleum Refinery Piping Code
for Pressure Piping ANSI B31.3-2002—Corrosion allowance = 0.05 inch
Nom Allowable working pressures for temperatures (in °F) not to exceed.
Weight of Wall ID Flow
pipe Sch. O.D.
pipe thk. (d) area –20 to
size No. in. 200 300 400 500 600 700
lb./ft. in. in. sq ft 100
in.
1/2 S40 .851 .840 .109 .622 .00211 2258 2258 2258 2258 2134 1953 1863
3/4 S40 1.131 1.050 .113 .824 .00371 1933 1933 1933 1933 1827 1672 1595
X80 1.474 .154 .742 .00300 3451 3451 3451 3451 3261 2985 2847
1 S40 1.679 1.315 .133 1.049 .00600 2103 2103 2103 2103 1988 1819 1735
X80 2.172 .179 .957 .00499 3468 3468 3468 3468 3277 3000 2861
160 2.844 .250 .815 .00362 5720 5720 5720 5720 5405 4948 4719
XX 3.659 .358 .599 .00196 9534 9534 9534 9534 9010 8247 7866
1-1/2 S40 2.718 1.900 .145 1.610 .01414 1672 1672 1672 1672 1580 1446 1379
X80 3.632 .200 1.500 .01225 2777 2777 2777 2777 2625 2402 2291
160 4.866 .281 1.338 .00976 4494 4494 4494 4494 4247 3887 3707
XX 6.409 .400 1.100 .00660 7228 7228 7228 7228 6831 6253 5963
2 S40 3.653 2.375 .154 2.067 .02330 1469 1469 1469 1469 1388 1270 1212
X80 5.022 .218 1.939 .02050 2488 2488 2488 2488 2351 2152 2053
160 7.445 .343 1.687 .01556 4600 4600 4600 4600 4347 3979 3795
XX 9.030 .436 1.503 .01232 6284 6284 6284 6284 5939 5436 5185
3 S40 7.58 3.500 .216 3.068 .05130 1640 1640 1640 1640 1550 1419 1353
X80 10.25 .300 2.900 .04587 2552 2552 2552 2552 2412 2207 2105
160 14.33 .438 2.624 .03755 4122 4122 4122 4122 3895 3566 3401
XX 18.58 .600 2.300 .02885 6089 6089 6089 6089 5754 5267 5024
4 S40 10.79 4.500 .237 4.026 .08840 1439 1439 1439 1439 1360 1244 1187
X80 14.99 .337 3.826 .07986 2275 2275 2275 2275 2150 1968 1877
160 22.51 .531 3.438 .06447 3978 3978 3978 3978 3760 3441 3282
XX 27.54 .674 3.152 .05419 5307 5307 5307 5307 5015 4590 4378
6 S40 18.98 6.625 .280 6.065 .2006 1205 1205 1205 1205 1139 1042 994
X80 28.58 .432 5.761 .1810 2062 2062 2062 2062 1948 1783 1701
160 45.30 .718 5.187 .1469 3753 3753 3753 3753 3546 3246 3097
XX 53.17 .864 4.897 .1308 4659 4659 4659 4659 4403 4030 3844
8 S40 28.56 8.625 .322 7.981 .3474 1098 1098 1098 1098 1037 950 906
X80 43.4 .500 7.625 .3171 1864 1864 1864 1864 1761 1612 1537
XX 72.4 .875 6.875 .2578 3554 3554 3554 3554 3359 3074 2932
160 74.7 .906 6.813 .2532 3699 3699 3699 3699 3496 3200 3052
10 S40 40.5 10.750 .365 10.020 .5475 1022 1022 1022 1022 966 884 843
X60 54.7 .500 9.750 .5185 1484 1484 1484 1484 1403 1284 1224
160 115.7 1.125 8.500 .3941 3736 3736 3736 3736 3531 3232 3082
12 S 49.6 12.750 .375 12.000 .7854 888 888 888 888 839 768 732
X 65.4 .500 11.750 .7528 1245 1245 1245 1245 1177 1077 1027
160 160.3 1.312 10.126 .5592 3699 3699 3699 3699 3496 3200 3052
14 10 36.7 14.000 .250 13.500 .9940 486 486 486 486 460 421 401
S30 54.6 .375 13.250 .9575 807 807 807 807 763 698 666
X 72.1 .500 13.000 .9211 1132 1132 1132 1132 1069 979 934
16 10 42.1 16.000 .250 15.500 1.310 425 425 425 425 402 368 351
S30 62.6 .375 15.250 1.268 705 705 705 705 666 609 581
S40 82.8 .500 15.000 1.227 987 987 987 987 933 854 815
18 10 47.4 18.000 .250 17.500 1.670 377 377 377 377 357 326 311
S 70.6 .375 17.250 1.622 625 625 625 625 591 541 516
X 93.5 .500 17.000 1.575 876 876 876 876 828 757 722
20 10 52.7 20.000 .250 19.500 2.074 339 339 339 339 321 293 280
S20 78.6 .375 19.250 2.021 562 562 562 562 531 486 464
X30 104.1 .500 19.000 1.969 787 787 787 787 743 680 649
24 10 63.4 24.000 .250 23.500 3.012 282 282 282 282 267 244 233
S20 94.6 .375 23.250 2.948 468 467 467 467 442 404 386
X 125.5 .500 23.000 2.883 660 654 654 654 618 565 539
Note: The above allowable working pressures are calculated from Fig. 17-23 using a reduction in tm to 87.5% of the wall thickness shown above to recognize mill
wall tolerance of 12.5%.
17-24
FIG. 17-27
Gas Transmission and Distribution Piping
Code for Pressure Piping ANSI B31.8-1999
Carbon Steel and High Yield Strength Pipe
(Values apply to A106, API 5L and API 5LX pipe having the same specified minimum yield strength as shown)
17-25
FIG. 17-27 (Cont’d.)
Gas Transmission and Distribution Piping
17-26
FIG. 17-28
Pressure-Temperature Ratings for Pipe Flanges and Flanged Fittings from ANSI B16.5-1996
Notes: (1) Upon prolonged exposure to temperatures above about 800°F (425°C), the carbide phase of carbon steel may
be converted to graphite: permissible but not recommended for prolonged use above 800°F
(2) Not to be used over 850°F
(3) Not to be used over 700°F
(4) Not to be used over 800°F
(5) At temperatures over 1000°F (540°C), use only when the carbon content is 0.04 percent or higher
(6) For temperatures above 1000°F (540°C), use only if the material is heat treated by heating it to a temp-
erature of at least 1900°F (1040°C) and quenching in water or rapidly cooling by other means
17-27
The table of allowable stresses, Fig. 17-25, is extracted from 16. Beggs, H. Dale, and James P. Brill. “A Study of Two-Phase Flow
ANSI B31.3. The designer is strongly urged to consult the lat- in Inclined Pipes” Trans. AIME, May 1973, pp. 606-617.
est ANSI B31.3 publication for full description of the code. A 17. Orkiszewski, J. “Predicting Two-Phase Pressure Drops in Verti-
tabular compilation of maximum allowable working pressures cal Pipe” J. Pet. Tech, AIME, June 1967, pp. 829-838.
calculated according to ANSI B31.3 appears in Fig. 17-26.
18. Baker, O., et al. “Gas-Liquid Flow in Pipelines, II. Design Manual”
For piping outside of refineries and other processing facili- AGA-API Project NX-28, October 1970.
ties, a separate code applies for determining allowable pres- 19. Dukler, A.E., Moye Wicks, III, and R.G. Cleveland. “Frictional
sure limits, ANSI B31.8. Refer to Fig. 17-24 for the calculation Pressure Drop in Two-Phase Flow: B. An Approach through Simi-
method. larity Analysis” AIChE Journal, Vol. 10, No. 1, January 1964, pp.
The designer is encouraged to refer to the latest ANSI B31.8 44-51.
standard for comprehensive code description. A tabular com- 20. Flanigan, Orin. “Effect of Uphill Flow on Pressure Drop in Design
pilation of maximum allowable working pressures computed of Two-Phase Gathering Systems” Oil and Gas Journal, March
according to ANSI B31.8 appears in Fig. 17-27. 10, 1958, pp. 132-141.
Fig. 17-28 provides pressure ratings for steel flanges and 21. Eaton, Ben A., et al. “The Prediction of Flow Patterns, Liquid
flanged fittings. Holdup and Pressure Losses Occurring During Continuous Two-
Phase Flow in Horizontal Pipelines” J. Pet. Tech. AIME, June
1967, pp. 815-828.
22. Greskovich, Eugene J., and Adam L. Shrier. “Slug Frequency in
REFERENCES Horizontal Gas-Liquid Slug Flow” Ind. Eng. Chem. Process De-
1. “Fluid Mechanics” by N. DeNevers, Addison-Wesley, Reading, velop., Vol. 11, No. 2, 1972, pp. 3;17-318.
Mass., 1970, p. 127. 23. Brill, James P., et al. “Analysis of Two-Phase Tests in Large Di-
2. “Experimental Research on the Flow of Water in Pipes” by H. ameter Prudhoe Bay Flowlines” SPE 8305, AIME, 1979.
Darcy, Mem. Acad. Sci. Inst. Imp. Fr., Vol. 15, 1858, p. 141 (in 24. Schmidt, Zelmir, and James P. Brill. “Experimental Study of Se-
French), “Mechanics of Engineering” by J. Weisbach, Van Nos- vere Slugging in a Two-Phase Flow Pipeline-Riser System” SPE
trand, 1872 (translated from German). 8306 AIME, 1979.
3. “Friction Factors for Pipe Flow” by L. Moody, Trans ASME 66, 25. Cunliffe, R.S. “Condensate Flow in Wet Gas Lines can be pre-
1944. dicted” Oil and Gas Journal, October 30, 1978, pp. 100-108.
4. C.F. Colebrook, “Turbulent Flow in Pipes with Particular Refer- 26. McDonald, Alvis E., and Ovid Baker. “Multiphase Flow in Pipe-
ence to the Transition Region Between the Smooth and Rough lines” Oil and Gas Journal, June 15, June 22, June 29, and July
Pipe Law,” J. Inst. Civil Engineers, London, 1939. 6, 1964
5. “Flow of Fluids through Valves, Fittings and Pipe,” by the Crane 27. Manual of Petroleum Measurement, Chapter 14.3, Orifice Meter-
Co., Chicago, Technical Paper No. 410, 1988. ing of Natural Gas and Other Related Hydrocarbon Fluids, Sec-
6. “Steady Flow in Gas Pipelines,” American Gas Association, IGT ond Edition, September, 1985 (AGA Report #3)(GPA 8185-85)
Technical Report 10, Chicago, 1965. (ANSI/API 2530).
7. T.R. Weymouth, Transactions of the American Society of Me- 28. Chemical Plant and Petroleum Refinery Piping, ANSI/ASME
chanical Engineers, Vol. 34, 1912. B31.3.
8. F.F. Richardson of Foreman-Richardson-Baird, The Petroleum 29. Gas Transmission and Distribution Piping Systems,
Engineer, July, 1945. ANSI/ASME B31.8.
9. F.N. Oliphant, “Production of Natural Gas,” Report of USGS, 30. Chen, N.H., "An Explicit Equation for Friction Factor in Pipe,"
1902. Ind. Eng. Chem. Fund., 18(3), 296 (1979).
10. “Fluid Flow, Pumps, Pipes and Channels” by N. Cheremisinoff,
Ann Arbor Science, Ann Arbor, Michigan, 1981, p. 184. BIBLIOGRAPHY
11. Piping Handbook, Sabin Crocker, McGraw-Hill Publishing Co.,
Inc., 1945. 1. R.W. Leach & W.P. Redmond, “How a Computer is Applied to a
Specific Problem in Pipeline Design,” Oil & Gas Journal, Flow
12. Mandhane, J.M., G.A. Gregory, and K. Aziz. “A Flow Pattern Map
Calculations in Pipelining, 1960.
for Gas-Liquid Flow in Horizontal Pipes” Int. J. Multiphase Flow,
Pergamon Press, Vol. 1, 1974, pp. 537-553. 2. Wilson, McAdams & Seltzer, “Flow of Fluids Through Commer-
cial Pipelines,” Journal of Industrial and Engineering Chemistry,
13. Taitel, Yehuda, and A.E. Dukler. “A Model for Predicting Flow
February, 1922.
Regime Transitions in Horizontal and Near Horizontal Gas-Liq-
uid Flow” AIChE Journal, Vol. 22, No. 1, Jan. 1976, pp. 47-55. 3. T.R. Aude, “Suggested Formula for Calculating Capacity of Prod-
ucts Pipe Lines,” 1943.
14. Aziz, K., G.W. Govier, and M. Fogarasi. “Pressure Drop in Wells
Producing Oil and Gas” J. Cdn. Pet. Tech., July-Sept 1972, pp. 4. Wilson & McAdams, “Contribution No. 19" from the Department
38-48. of Chemical Engineering, Massachusetts Institute of Technology.
15. Lockhart, R.W., and R.C. Martinelli, “Proposed Correlation of 5. W.G. Heltzel, “Fluid Flow and Friction in Pipelines,” Oil & Gas
Data for Isothermal Two-Phase, Two-Component Flow in Pipes” Journal, June 5, 1930.
Chemical Engineering Progress, Vol. 45, No. 1, January 1949, pp. 6. T.R. Young, “Digital Simulation of Crude Oil Pipelines,” API Pipe-
39-48. line report, May 2, 1960.
17-28