Pressure Drop2016
Pressure Drop2016
Pressure Drop2016
|2016|
a
y
It expands to fill the vessel that contains
r
the gas.
A
h
s
r
a
d
A
.
r
D
Gases are Compressible.
July 2014
Aug
2014
PES
2016
Properties of Gases
1. Mass
2. Volume
3. Specific Gravity
4. Viscosity
d
6. Compressibility factor . A
r
7. Average Pressure D
Calculations
5. Average Molecular Weight
h
s
ar
a
y
r
A
8. Heating Value
July 2014
Aug
2014
PES
2016
MASS
Mass is the quantity of matter.
For representing units, two systems will be used :
a
y
r
A
d
A
.
r kg)
(1 slug =D
14.5
h
s
ar
July 2014
Aug
2014
PES
2016
Standard Conditions
Standard conditions (also called base conditions) of
a
y
r
A
h
s
r are 15
Standard conditions in SI system
a
d
A
.
KPa.
r
D
July 2014
C & 101.325
Aug
2014
PES
2016
Volume
Volume of gas is measured in ft3 in USCS units and m3 in SI
units.
a
y
Other units for volume In USCS include thousand
r
ft (Mft
A
h
and million ft (MMft )
s
r
a
d
and thousand m (km ) and
million
m (Mm ) in SI units
A
.
r
DMillion standard cubic feet per day.
MMSCFD stands for
3
July 2014
Aug
2014
PES
2016
a
y
r
In SI :
k (kilo) represents a thousand)
A
h
s
r
M refers to million.
a
d
A
.
Therefore, 500 MSCFD
in USCS units refers to 500 thousand
r
D
MM refers to million
July 2014
Aug
2014
PES
2016
SPECIFIC GRAVITY
g
G
air
h
s
ar
a
y
r
A
d
M
G r. A
DM 28.9625
Mg
air
July 2014
Aug
2014
PES
2016
VISCOSITY
a
y
calculate the viscosity from the viscosities ofrcomponent gases:
A
h
s
r
a
d
A
.
r
D i yi M i
July 2014
Mi
Aug
2014
PES
2016
1 y1 M 1 2 y2 M 2 3 y3
y1
July 2014
a
y
r
M
...
y
A
h
s
M y M ry M
a
d
A
.
r
D
1
n n
Mn
... yn M n
Aug
2014
PES
2016
a
y
r
A
h
s
ar
d
A
.
r
D
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
M a yi M i
h
s
ar
a
y
r
A
d
A
.
r
M = Molecular Weight
of
individual
gas
component.
D
yi=Mole fraction of gas component.
i
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
COMPRESSIBILITY FACTOR
It is a measure of how much the properties of a real gas deviate
from ideal gas.
a
y
The compressibility factor is defined as the rratio of the actual
A
h
gas volume to the volume the gas s
would
occupy if it were an
r
a
ideal gas at the same temperature
and pressure.
d
A
.
r
D
July 2014
Aug
2014
PES
2016
a
y
r
A
h
PV ZnRT
s
r
a
d
A
.
r
D
July 2014
Aug
2014
PES
2016
h
s
ar
a
y
r
A
d
A
.
d. CNGA method (Most
commonly used )
r
D
c. AGA method
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
A
h
s
ar
d
A
.
r
T = average
D gas temperature, R (1R=0.556 K)
Pavg = average gas pressure, psig
f
July 2014
Aug
2014
PES
2016
a
y
r
A
Pav
July 2014
P
2
h
i
j
s
r
Aid Pj
3 r.
Pi Pj
Aug
2014
PES
2016
a
y
r
A
and a temperature of 70F, using the CNGA
method.
h
s
r
a
d
A
.
r
D
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
HEATING VALUE
a
y
r
A
h
s
ar
d
A
.
r
D
gases and their mole fractions using the following equation:
H m Hi M i
July 2014
Aug
2014
PES
2016
Where:
a
y
r
y = mole fraction or percent of gas component
i
A
h
s
Hi = heating value of gas component,
Btu/ft
r
a
d
A
.
r
D
Hm = gross heating value of mixture, Btu/ft3
i
July 2014
Aug
2014
PES
2016
h
s
ar
a
y
r
A
d
A
DN is nominal diameter
in mm.
.
r
D
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
FLOW EQUATIONS
1. General Flow equation.
Colebrook-White equation
Modified Colebrook-White equation
AGA equation
2. Weymouth equation.
d
A
.
r
D
3. Panhandle A equation
4. Panhandle B equation
5. IGT equation
6. Spitzglass equation
h
s
ar
a
y
r
A
7. Mueller equation
8. Fritzsche equation
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Tb
Q 77.54
Pb
P P
GT
LZf
f
2
1
d
A
.
r
D
July 2014
2
2
0.5
h
s
ar
a
y
r
A
D 2.5
Aug
2014
PES
2016
Tb
Q 1.1494*10
Pb
3
P P
GT f LZf
d
A
.
r
D
July 2014
2
1
2
2
h
s
ar
0.5
a
y
r
A
D 2.5
Aug
2014
PES
2016
Transmission Factor
h
2
s
r
a
F
d
A
.
r
f
D
July 2014
a
y
r
A
Aug
2014
PES
2016
In USCS
a
y
r
A
Tb P12 P22
Q 38.77 F rsh
a
P
GT
LZ
f
b d
0.5
2.5
A
.
r
D
July 2014
Aug
2014
PES
2016
In SI System
a
0.5
y
2 r 2
Tb P1 A P2
4
2.5
Q 5.747*10 F s h
*D
P
GT
LZ
a
f
b
d
A
.
r
D
July 2014
Aug
2014
PES
2016
In USCS
a
y
0.5
r
s A2
Tb P12 h
e P2
2.5
s
Q 38.77 F
D
a
P
GT
L
Z
b A d f e
.
r
D
July 2014
Aug
2014
PES
2016
In SI.
a
0.5
y
2
s r 2
P1 h A
e P2
4 Tb
2.5
Q 5.747*10 F rs
*D
a
P
GT
L
Z
f e
Adb
.
r
D
July 2014
Aug
2014
PES
2016
Where
L(e 1)rya
A
Le
h
s
s
ar
s
d
A
.
r
D
July 2014
Aug
2014
PES
2016
In USCS
H 2 H1
a
s 0.0375G
y
r
hT fA
s
r
a
d
A
.
r
D
July 2014
Aug
2014
PES
2016
H 2 H1
a
s 0.0684G
y
Tf Z A
r
s
r
a
d
A
.
r
D
July 2014
Aug
2014
PES
2016
a
y
r
A
h
s
r of length L has a series
If, however, the pipe segment
a
d
A
.
of slopes, then werintroduce a parameter j as follows
D
for each individual pipe sub segment that constitutes
the pipe length from point 1 to point 2.
July 2014
Aug
2014
PES
2016
a
y
r
A
e h 1
ji ars
d
A
s
.
i
r
D
si
July 2014
Aug
2014
PES
2016
a
y
r
A
h
s
ar
d
A
.
r
The equivalent length
D term L is calculated as:
e
Le j1 L1 j2 L2 e j3 L3e ...
s1
July 2014
s2
Aug
2014
PES
2016
a
y
r
A
h
s
ar
d
A
.
segment, from ther pipeline inlet to the end of each
D
segment.
July 2014
Aug
2014
PES
2016
Qb
u 0.002122 2
D
d
A
.
r
D
July 2014
Pb ZT
Tb P
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Qb Pb ZT a
u 14.7359 2 ry
D Tbh
AP
s
r
a
d
A
.
r
D
July 2014
Aug
2014
PES
2016
EROSIONAL VELOCITY
It is always tried to keep the flow rate of gas and hence the
velocity as high.
a
y
But a high velocity in pipeline leads to vibration
r
and noise.
A
h
s
In addition, higher velocities will
cause erosion of the pipe
r
a
interior over a long periodA
ofd
time.
.
r
D
July 2014
Aug
2014
PES
2016
umax
100
a
y
r
A
h
s
r
a
d
A
.
r
D
Where
July 2014
Aug
2014
PES
2016
h
s
ar
a
y
r
A
d
A
.
r
D
ZRT
umax 100
July 2014
29GP
Aug
2014
PES
2016
Where
Z = compressibility factor of gas, dimensionless
R = gas constant = 10.73 ft3 psia/lb-moleR
T = gas temperature, R
G = gas gravity (air = 1.00)
d
A
.
r
D
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
A
h
s
r
a
pipe if the inlet pressure is 1000 psig and the outlet pressure is
d
A
.
850 psig. The base pressure
and base temperature are 14.7 psia
r
D
and 60F, respectively. Assume compressibility factor Z = 1.00.
flow, calculate the velocity of gas at the inlet and outlet of the
What is the erosional velocity for this pipeline based on the above
data and a compressibility factor Z = 0.90?
July 2014
Aug
2014
PES
2016
Soln.
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Reynolds No.
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Ques.
A natural gas pipeline, NPS 20 with 0.500 in. wall thickness,
a
y
r
A
0.6 and viscosity is 0.000008 lb/ft-s. Calculate
the value of the
h
s
r
a
Reynolds number of flow. Assume the base temperature and
d
A
.
r
base pressure are 60F
D and 14.7 psia, respectively.
transports 100 MMSCFD of gas. The specific gravity of gas is
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
FRICTION FACTOR
The term friction factor is a dimensionless parameter that
depends upon the Reynolds number of flow.
a
y
Two types of friction factor are commonly used
r
:
A
h
s
i) Darcy friction factor.
ii) Fanning friction factor.
r
a
d
A
Both friction factors are .co-related
by the following equation :
r
D
fD
ff
4
July 2014
Aug
2014
PES
2016
h
s
ar
July 2014
d
A
.
r
D
a
y
r
A
Aug
2014
PES
2016
Laminar Flow
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Turbulent Flow
For turbulent flow, many empirical relationships for calculating
f have been put forth by researchers.
h
s
ar
d
Modified The Colebrook-White
A
.
r
D
AGA equations
i) The Colebrook-White
iii)
a
y
r
A
July 2014
ii)
Aug
2014
PES
2016
The Colebrook-White
Turbulent flow in pipes (Re > 4000) is subdivided into three
separate regions as follows:
h
s
ar
a
y
r
A
d
A
.
r
D
3. Turbulent flow in fully rough pipes.
2. Turbulent flow in smooth pipes.
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Where
f = friction factor, dimensionless
h
s
r
D = pipe inside diameter, in.
a
d
A
.
r
Droughness, in.
e = absolute pipe
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
A
h
s
ar
July 2014
d
A
.
r
D
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
A
use for many years in both liquid flow and gas flow.
h
s
rform of the Colebrookthat introduced a modified
a
d
A
.
White equation. r
D
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Ques.
Q. A gas pipeline, NPS 20 with 0.500 in. wall thickness, flows
a
y
r
A
0.000008 lb/ft-s. Using the modified Colebrook-White
equation,
h
s
r
a
calculate the pressure drop in a 50 mi segment of pipe, based on
d
A
.
an upstream pressureD
ofr1000 psig.
200 MMSCFD gas of specific gravity 0.6 and viscosity of
July 2014
Aug
2014
PES
2016
a
y
r
psia, respectively. Neglect elevation
effects and use
A
h
s
60F for gas flowing temperature
and compressibility
r
a
d
A
factor Z = 0.88.
.
r
D
July 2014
Aug
2014
PES
2016
Sol.
Inside diameter of pipe = 20 2 . 0.5 = 19.0 in.
The base temperature = 60 + 460 = 520 R
Gas flow temperature = 60 + 460 = 520 R
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
By Successive iterations:
F=19.81
Now using Generalized flow equation,
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
P2=853.23psia
= 838.5psig
Hence pressure drop in the pipeline
= 1014.73-853.23
=161.5psia.
July 2014
d
A
.
r
D
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
transmission factor for gas pipelines to beA
used in the General
h
s
r
a
Flow equation.
d
A
.
r
This is sometimes D
referred to as the AGA NB-13 method.
Aug
2014
PES
2016
h
s
ar
a
y
r
A
d
A
.
r
D
July 2014
Aug
2014
PES
2016
a
y
r
A
d
A
.
r
D
July 2014
h
s
ar
Aug
2014
PES
2016
For the partially turbulent zone, F is calculated from the following equations
a
y
r
A
using the Reynolds number, a parameter Df known as the pipe drag factor,
and the Von Karman smooth pipe transmission factor Ft:
d
A
.
r
D
July 2014
h
s
ar
Aug
2014
PES
2016
The pipe drag factor Df is a parameter that takes into account the
a
y
r
A
number of bends, fittings, etc. Its value ranges from 0.90 to 0.99.
h
s
ar
The Bend index is the sum of all the angles and bends in the pipe
segment, divided by the total length of the pipe section under
consideration
July 2014
d
A
.
r
D
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Ques.
Using the AGA method, calculate the transmission factor and
a
y
r
0.500 in. wall thickness. The flow rate isA
200 MMSCFD, gas
h
s
r
a
gravity = 0.6, and viscosity = 0.000008 lb/ft-sec. The absolute
d
A
.
pipe roughness is D
700r in. Assume a bend index of 60, base
friction factor for gas flow in a bare steel NPS 20 pipeline with
July 2014
Aug
2014
PES
2016
Sol.
Inside diameter of pipe = 20 2 . 0.5 = 19.0 in.
The base temperature = 60 + 460 = 520R
d
A
.
r
D
h
s
ar
a
y
r
A
July 2014
Aug
2014
PES
2016
h
s
ar
a
y
r
A
d
A
.
r
D
July 2014
Aug
2014
PES
2016
a
y
r
A
h
s
Therefore, the corresponding friction
factor f is found
r
a
d
A
.
r
D
July 2014
Aug
2014
PES
2016
Ques.
Using the AGA method, calculate the transmission factor and friction factor
for gas flow in a DN 500 pipeline with 12 mm wall thickness. The flow rate
a
y
r
A
is 6 Mm3/day, gas gravity = 0.6, and viscosity = 0.00012 Poise. The absolute
h
s
ar
pipe roughness is 0.02 mm. Assume a bend index of 60, base pressure of
d
A
.
r
D
101 kPa, and base temperature of 15C. For a 60 km pipe length, calculate
the upstream pressure needed to hold a downstream pressure of 5 MPa
(absolute). Assume flow temperature = 20C and compressibility factor Z =
0.85. Neglect elevation effects. (DN 500= 500mm)
July 2014
Aug
2014
PES
2016
Sol.
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
WEYMOUTH EQUATION
Gas flows in pipelines with diameters of 15 inches or below are
a
y
r
A
h
s
ar
d
A
.
r
D T P
e
P
b
Q 433.5 E
P
GT
L
Z
f e
b
2
1
July 2014
2
2
0.5
*D
2.667
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
In SI units:
a
y
r
T P eAP
Q 3.7435*10 E sh
rGT L Z
P
d
a
A
.
r
D
3
July 2014
2
1
2
2
0.5
*D
2.667
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
PANHANDLE A EQUATION
Incorporates an efficiency factor instead of pipe friction.
a
y
r
A
h
s
ar
d
A
.
r
D
USCS units as follows:
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
In SI system
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Ques.2
Using the Panhandle A equation, calculate the outlet pressure in
a
y
r
A
miles long. The gas flow rate is 100 MMSCFD
at an inlet
h
s
r
a
pressure of 1000 psia. The gas gravity = 0.6 and viscosity =
d
A
.
0.000008 lb/ft-sec.Dr
July 2014
Aug
2014
PES
2016
a
y
r
A
h
s
ar
0.92.
d
A
.
r
(Accuracy to be maintained
D for downstream pressure= 0.5 psia
As first approximation use a downstream pressure of 800psia
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
PANHANDLE B EQUATION
The Panhandle B equation, also known as the revised Panhandle equation, is
used in large diameter, high pressure transmission lines.
a
y
r
A
h
s
ar
d
A
.
r
D
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Ques.3
Using the Panhandle B equation, calculate the outlet pressure in
a
y
r
A
miles long. The gas flow rate is 100 MMSCFD
at 1000 psia
h
s
r
a
inlet pressure. The gas gravity = 0.6 and viscosity = 0.000008
d
A
.
r
lb/ft-sec. The average
D gas temperature is 80F. Assume base
a natural gas pipeline, NPS 16 with 0.250 in. wall thickness, 15
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
A
h
s
ar
d
A
.
r
D
pressure was calculated using the different flow equations
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Q1. Using the IGT equation, calculate the flow rate in a natural
a
y
r
The inlet and outlet pressure are 1000A
psig and 800 psig,
h
s
r
a
respectively. The gas gravity = 0.6 and viscosity = 0.000008
d
A
.
r
lb/ft-s. The average
gas
D temperature is 80F, base pressure =
gas pipeline, NPS 16 with 0.250 in. wall thickness, 15 mi long.
Aug
2014
PES
2016
SOL.
Inside diameter of pipe = 16 2 . 0.25 = 15.5 in.
The pressures given are in psig, and they must be converted to
absolute pressures.
P1 = 1000 + 14.7 = 1014.7 psia
P2 = 800 + 14.7 = 814.7 psia
d
A
.
r
D
T = 80 + 460 = 540R
Tb = 60 + 460 = 520R
h
s
ar
a
y
r
A
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
SPITZGLASS EQUATION
Used in fuel gas piping calculations
a
y
r
One equation is for low pressure (less than
or
equal to 1 psig)
A
h
s
and another is for high pressure (more
than 1 psig).
r
a
d
A
.
r
D
There are two versions of the Spitzglass equation.
July 2014
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
a
y
r
A
ft. The flowing temperature of fuel gas is
60F, and the inlet
h
s
r
a
pressure is 1.0 psig. In the pipeline the pressure decreases to
d
A
.
0.0252 psig and D
the rspecific gravity of gas = 0.6. Assume
inside diameter of 6.065 in. and a total equivalent length of 180
Aug
2014
PES
2016
Base Temperature=60+460=5200R
Gas Temperature= 520 R
Pressure Drop=1-0.025=0.9747psi.
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
MUELLER EQUATION
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
FRITZSCHE EQUATION
d
A
.
r
D
July 2014
h
s
ar
a
y
r
A
Aug
2014
PES
2016
Conclusion
Highest pressure drop is predicted by the Weymouth equation
a
y
r
A
h
s
ar
d
A
.
r
D drop is the Weymouth equation and the
the highest pressure
July 2014
Aug
2014
PES
2016