Ur Process Pocket
Ur Process Pocket
This tool was developed for process engineers to quickly estimate line pressure drops (single-phase gas
or liquid) estimate equivalent length for pipelines, and access a list of standards checklists on site. As a
part of a larger project to develop what is similar to Carl Branan’s book “process engineers Pocket
Handbook” these tools would allow a process engineer to quickly calculate/estimate equipment
efficiencies or sizing using standardized calculations
The aim here is to take little-known data from the field (flow, pressures, temperatures,
compositions..etc.) and use it as input for a rough estimation without having to return to the office to
use commercial software or calculations Excel sheets to validate or to calculate. Additionally, These tools
may also serve as a gathered data validation tool.
Nomenclature
Input Required
Table 1: Note 1: Input required depending on the case
Output Obtained
s: specific gravity
u: average velocity (ft/s)
L: Equivalent length (ft)
d: pipe diameter (in)
f: fannings friction factor (equals ¼ moody’s friction factor)
𝑅𝑒 : Reynold’s number
Z: viscosity (cP)
Note: fannings equations isn’t valid for gases when pressure drop is larger than 10% of
downstream pressure
Total equivalent length = actual length + fittings estimated length
Nelson’s fittings equivalent length[6]
Equivalent length table P.396
Reynold’s number
𝝆𝒗𝒅
𝑹𝒆 =
𝝁
μ: Viscosity (kg/m.s) or 0.001*cP
v: fluid velocity (m/s)
Darcy friction factor
You either opt to use the famous moody’s chart or use the Colebrook-white
equation
Laminar flow
𝟔𝟒
𝒇=
𝑹
Turbulent flow
Colebrook-White equation
𝟏 𝒆 𝟐. 𝟓𝟏
= −𝟐 𝐥𝐨𝐠 𝟏𝟎 ( + )
√𝒇 𝟑. 𝟕𝑫 𝑹√𝒇
When the piping layout is not available, the equivalent length (Le) of the piping will be
estimated based on the straight length (Ls) as follows:
Note that Ls is the sum of XYZ coordinate length. For large size or high pressure piping, it is
recommended to estimate the number of elbows tees and valves and evaluate the equivalent
length, assuming piping layout. [5]
2. GPSA table[1]:
This table is the table used to estimate your equivalent length in the tool developed
16 𝑄 1.852 𝑆. 𝑔
∆𝑃𝑓 = 0.98066 ∗ 1.1101 ∗ 10 ∗( )
𝐶 𝐷4.87
Where: ∆𝑃𝑓 : frictional pressure drop (kg/cm2)
C: Hazen-williams C factor
Hazen-williams C factor
Pipe material Hazen-williams C factor
Smooth pipes (all metals) 130-140
Cast iron (old) 100
Iron (worn/pitted) 60-80
Polyvinyl chloride (PVC) 150
Brick 100
Smooth wood 120
Smooth masonry 120
Vitrified clay 110
Note: This equation is yet to be utilized in our tool
1.2 Gases pressure drop calculations
Several equations are available that relate the gas flow rate with gas properties, pipe diameter and
length, and upstream and downstream pressures. These equations are listed as follows:
1. General Flow equation
2. Weymouth equation
3. Panhandle A equation
4. Panhandle B equation
Other equations which aren’t discussed here: Oliphant, Fritzsche, Muller, IGT, Spitzglass high, Spitzglass
low [4]
Where:
Q: gas flow rate (stdm3/hr)
Tb: base temperature K (15.5556 oC +273.15 in tool)
Pb: Base pressure kPa (101.325 Kpa in tool)
Tf: average flow temperature K
P1: Upstream pressure kPa
P2: Downstream pressure kPa
L: Pipe length (km)
Le: Equivalent length of pipe (km)
H1: Upstream elevation (m)
H2: Downstream elevation (m)
G: gas gravity (air =1)
f: Darcy’s friction factor
D: pipe inside diameter, mm.
Z: Gas compressibility
The equivalent length, Le, and the term 𝑒 𝑠 take into account the elevation difference
between the upstream and downstream ends of the pipe segment. The parameter s depends upon the
gas gravity, gas compressibility factor, the flowing temperature, and the elevation difference.
𝑷𝒃 𝑮𝑸
𝑹𝒆 = 𝟎. 𝟓𝟏𝟑𝟒 ∗ 𝟐𝟒 ∗ 𝟏𝟎𝟎 ∗ ∗
𝑻𝒃 𝝁𝑫
𝝁: viscosity in cP
Z is determined based on average P where:
𝟐 𝑷𝟏𝑷𝟐
𝑷𝒂𝒗𝒈. = (𝑷𝟏 + 𝑷𝟐 − )
𝟑 𝑷𝟏 + 𝑷𝟐
0.5
−3
𝑇𝑏 𝑃12 − 𝑃22
𝑄 = 3.7435 ∗ 10 ∗ 24 ∗ 𝐸 ∗ ( ) ∗ 𝐷2.667
𝑃𝑏 𝐺 𝑇𝑓 𝐿𝑒 𝑍
remaining units are the same for the general flow equation
Pipleline efficiency E
Note: API RP 14E doesn't use E for the Weymouth equation.
0.5394
𝑇𝑏 1.0788
𝑃12 − 𝑃22
𝑄 = 4.5965 ∗ 10−3 ∗ 24 ∗ 𝐸 ∗ ( ) ( ) ∗ 𝐷2.6182
𝑃𝑏 𝐺0.8539 𝑇𝑓 𝐿𝑒 𝑍
remaining units are the same for the general flow equation
The successful application of these transmission line flow equations in the past has largely involved
compensation for discrepancies through the use of adjustment factors, usually termed “efficiencies.”
These efficiencies are frequently found in practice by determining the constant required
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, under 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 expression
with a preconditioned bias. Regardless of the merits of various gas flow equations, past
practices may dictate the use of a particular equation to maintain continuity of comparative capacities
through application of consistent operating policy. A summary of comparisons between transmission
factors used in the above gas equations are shown in Fig. 17-7. [4]
In Figure 2.5, we consider a pipeline 100 mi long, NPS 16 with 0.250 in. wall thickness, operating at a
flow rate of 100 MMSCFD. The gas flowing temperature is 80°F. With the upstream pressure fixed at
1400 psig, the downstream pressure was calculated using the different flow equations. By examining
Figure 2.5, it is clear that the highest pressure drop is predicted by the Weymouth equation and the
lowest pressure drop is predicted by the Panhandle B equation. It must be noted that we used a pipe
roughness of 700 μin. for both the AGA and Colebrook equations, whereas a pipeline efficiency of 0.95
was used in the Panhandle and Weymouth equations. Figure 2.6 shows a comparison of the flow
equations from a different perspective. In this case, we calculated the upstream pressure required for
an NPS 30 pipeline, 100 miles long, holding the delivery pressure constant at 800 psig. The upstream
pressure required for various flow rates, ranging from 200 to 600 MMSCFD, was calculated using the five
flow equations. Again it can be seen that the Weymouth equation predicts the highest upstream
pressure at any flow rate, whereas the Panhandle A equation calculates the least pressure. We
therefore conclude that the most conservative flow equation that predicts the highest pressure drop is
the Weymouth equation and the least conservative flow equation is Panhandle A.
Remaining calculations
Liquid Calculations
NPSHa = Psuction – Pv = P1 - ∆𝑃𝑓 + H – Pv
(P1 and Pv variables in kg/cm2.a and P1, Pv & - ∆𝑃𝑓 converted to m.liq by *10/S.Gr)
Gas Calculations
𝑷𝒃 (𝒌𝒈/𝒄𝒎𝟐 ) (𝟏. 𝟎𝟑𝟑) 𝑻 (°𝑲)
𝑸𝒂𝒄𝒕𝒖𝒂𝒍 = 𝑸𝒔𝒕𝒅𝒎𝟑 /𝒉𝒓 × ×
𝑷 (𝒌𝒈/𝒄𝒎𝟐 ) 𝑻𝒃 (°𝑲)(𝟏𝟓. 𝟓𝟔 °𝑪)
Similarly, 𝑸𝒏𝒐𝒓𝒎𝒂𝒍 (𝒏𝒎𝟑/𝒉𝒓) can be calculated at T = 0 +273.15 K and P=Pb
𝑸𝒂𝒄𝒕𝒖𝒂𝒍
𝑽𝒆𝒍𝒐𝒄𝒊𝒕𝒚 =
𝑨𝒄𝒓𝒐𝒔𝒔 𝒔𝒆𝒄𝒕𝒊𝒐𝒏𝒂𝒍
𝟗. 𝟖𝟏 ∗ 𝒌 ∗ 𝑹 ∗ 𝑻
𝑺𝒐𝒏𝒊𝒄 𝒗𝒆𝒍𝒐𝒄𝒊𝒕𝒚 𝒗𝒔 = √
𝑴𝒘𝒕
𝒗
𝒎𝒂𝒄𝒉 𝒏𝒖𝒎𝒃𝒆𝒓 =
𝒗𝒔
𝑪𝒆
𝑬𝒓𝒓𝒐𝒔𝒊𝒐𝒏𝒂𝒍 𝒗𝒆𝒍𝒐𝒄𝒊𝒕𝒚 𝑽𝒆 =
√𝝆
Table 2: Source: JGC
Ce Services
180-240 Solid-free fluids where corrosion is not anticipated or when corrosion is controlled by
inhibition or by employing corrosion resistant alloys, continuous.
300 Same as above but intermittent
120 Solids and/or corrosives are expected
Table 3: C = 100 for continous and C=125 for itermittent services are conservative [8]
Where:
R: gas constant 847.9 (kgf/m2)(m3)/(kg-mole)(°K)
K: Cp/Cv
If [5]:
1 The flow is subsonic
Mach number <
√𝑘
1 The flow is sonic
Mach number =
√𝑘
1 The flow is supersonice
Mach number >
√𝑘
𝑇𝑐 = ∑ 𝑦𝑖 ∗ 𝑇𝑐𝑖
𝑖=1
𝑛
𝑃𝑐 = ∑ 𝑦𝑖 ∗ 𝑃𝑐𝑖
𝑖=1
𝑇 𝑃
𝑇𝑟 = , Pr =
𝑇𝑐 𝑃𝑐
Where Tci and Pci are each component’s critical temperature and pressure
Naturally, The Standing and Katz method is applied where Z is obtained as a function of Tr and Pr
through the famous chart (figure no. ). Alternatively, Dranchuk, Purvis, and Robinson method is used. In
this method the Benedict-Webb-Rubin equation of state is used to correlate the Standing-Katz Z
factor chart. Eight coefficients A1, A2, etc., are used in this equation as shown [1]:
𝑨𝟐 𝑨𝟑 𝑨𝟓 𝑨𝟓 𝑨𝟔 𝝆𝒓 𝟓 𝑨𝟕 𝝆𝒓 𝟑
𝒁 = 𝟏 + (𝑨𝟏 + + ) 𝝆𝒓 + (𝑨𝟒 + ) 𝝆𝒓 𝟐 + +
𝑻𝒑𝒓 𝑻𝟑 𝑻𝒑𝒓 𝑻𝒑𝒓 𝑻𝒑𝒓 𝟑 (𝟏 + 𝑨𝟖 𝝆𝒓 𝟐 )𝒆𝒙𝒑(−𝑨𝟖 𝝆𝒓 𝟐 )
Where:
𝟎. 𝟐𝟕 ∗ 𝑷𝝆𝒓
𝝆𝒓 =
𝒁𝑻𝝆𝒓
K (Cp/Cv) [1]
𝐶𝑝 𝑀𝐶𝑝
𝑘= = , 𝑅 = 1.986 𝐵𝑡𝑢/(𝑙𝑏𝑚𝑜𝑙 ˚𝐹)
𝐶𝑣 𝑀𝐶𝑝 − 𝑅
∑𝑛𝑖=1 𝑦𝑖 ∗ 𝑀𝐶𝑝
𝑘=
(∑𝑛𝑖=1 𝑦𝑖 ∗ 𝑀𝐶𝑝) − 1.986
MCp values at different temperatures for a variety of compounds are found in figure 1 [1]
Figure 4: Mcp, Tc and Pc values [1]
Example No.1: Liquid pressure drop calculations (Units converted to match table units)
Computer Methods for chemical engineers (2nd edition Page no. 37) (input highlighted)
Table 4: Note: Software pgk utilized is PR, using ASME stem is recommended though for this particular case as fluid is Water!
Manual calculations:
𝑄 8∗4
𝑣= = = 4 𝑚/𝑠
𝐴 3600 ∗ 3.14 ∗ 0.02662
𝜌𝑣𝑑 1007 ∗ 4 ∗ 0.0266
𝑅= = = 120333.33
𝜇 0.8904 ∗ 0.001
𝜀 0.00005
The relative roughness of the smooth pipe is =
𝑑 0.0266
Turbulent flow, Iterate through the following function to obtain f using the goal
seek function in excel or use moody’s chart
Colebrook-White equation
1 𝑒 2.51
= −2 log10 ( + )
√𝑓 3.7𝐷 𝑅√𝑓
f = 0.0245
𝑓 ∗ 𝜌 ∗ 𝑢2 ∗ 𝐿 1.01972 ∗ 10−5 ∗ 0.0245 ∗ 1007 ∗ 42 ∗ 10
−5
∆𝑃𝑓 = 1.01972 ∗ 10 ∗ =
𝑑2 0.0266 ∗ 2
= 0.751 𝑘𝑔/𝑐𝑚2
P2 = P1 - ∆𝑃𝑓 = 0.282 kg/cm2
Calculations Examples using commercial software, hand written calculations and the tool developed:
Example No.2: Gas pressure drop calculations (Units converted to match table units)
Computer Methods for chemical engineers (2nd edition Page no. 43) (input highlighted)
Hand Tool
Commercial
Property Unit Calculations Developed
software[Note]
in book (NiqSim)
P1 Kg/cm2 3.518 3.518 3.518
P2 Kg/cm2 3.384 3.334 3.376
o
T C 25 25 25
L m 20000 20000 20000
D m 0.154 0.154 0.154
flow_rate Std_m3/hr 424.1 363 kg/hr 424.1
Composition Vol% 85 CH4 15 CO2 85 CH4 15 CO2 85 CH4 15 CO2
flow_rate Nm3/hr - - 402.02
Actual_flow_rate m3/hr - - 133.99
Dp% to % 3.8
Upstream 4.022
pressure
dp/100 m Kg/cm2/100m 0.00067 - 0.0007
friction factor/E - 0.024 -
Reynolds 70480 71400
71670
number
Out velocity m/s 1.98 1.98 1.99
Sonic Velocity m/s 398.8 - 399.24
Out Mach 0.005 -
0.005
number
ꝬV^2 - - 11.527
mol.weight 20.24 - 20.238
z 0.9931 - 0.993
k 1.308 - 1.30
1/root(k) - - 0.876
viscosity cP 0.0118 0.012 0.0116
Table 5: Note: Gas pipe for compressible gases were used, adiabatic pipe could also be used with similar results
Manual Calculations
General gas equation
K, Pr and Tr calculations
Carbon
Gas Methane dioxide total
Mol wt 16.043 44.01 Tc 207.59 K
Pc 667 1070 Pc 51.144952 kg/cm2.a
Tc -116.66 87.76 M.wt 20.23805 g/gmol
0F 8.23 8.38
50 F 8.42 8.7 T1 298.15 K
60 F 8.46 8.76 76.37 F
100 F 8.65 9 P1 3.518 kg/cm2.a
150 F 8.95 9.29
200 F 9.28 9.56
250 F 9.64 9.81
300 F 10.01 10.05
mol fraction 0.85 0.15 1
Tc -85.997 F P1
Pc 727.45 Psia Tr 1.436
M.wt 20.23805 g/gmol Pr 0.068
P1 MCp 8.505
K1 1.304648
Viscosity calculations
Viscosity at 1 atm µa (cP) Fig. 32.-22 0.0104
Viscosity addition due to CO2 15 mol% Fig. 32.-22 0.0008
Total µa (cP) 0.0112
µ/µa at Pr = 0.068 and Tr = 1.54 Fig. 32.-24 1.05
µ (cP) 1.05*0.0112 0.01176
Z calculations
𝑨𝟐 𝑨𝟑 𝑨𝟓 𝟐
𝑨𝟓 𝑨𝟔 𝝆𝒓 𝟓 𝑨𝟕 𝝆𝒓 𝟑
𝒁 = 𝟏 + (𝑨𝟏 + + ) 𝝆 + (𝑨𝟒 + )𝝆 + +
𝑻𝒑𝒓 𝑻𝟑 𝒓 𝑻𝒑𝒓 𝒓 𝑻𝒑𝒓 𝑻𝒑𝒓 𝟑 (𝟏 + 𝑨𝟖 𝝆𝒓 𝟐 )𝒆𝒙𝒑(−𝑨𝟖 𝝆𝒓 𝟐 )
Where:
𝟎. 𝟐𝟕 ∗ 𝑷𝝆𝒓
𝝆𝒓 =
𝒁𝑻𝝆𝒓
(Use the Goal seek function in an Excel Sheet, Change your assumed Z till ABS(Z_assumed –
Z_calculated)*1000 = 0 for instance)
Z = 0.992
𝑤𝑡 20.24
𝐺 = 𝑀. = = 0.698
29 29
𝑷𝒃 𝑮𝑸 𝟏𝟎𝟏.𝟑𝟐𝟓 𝟎.𝟔𝟗𝟖∗𝟒𝟐𝟒.𝟏
𝑹𝒆 = 𝟎. 𝟓𝟏𝟑𝟒 ∗ 𝟐𝟒 ∗ 𝟏𝟎𝟎 ∗ ∗ = 𝟎. 𝟏𝟓𝟑𝟒 ∗ 𝟐𝟒 ∗ 𝟏𝟎𝟎 ∗ ∗ =70714
𝑻𝒃 𝝁𝑫 𝟐𝟕𝟑+𝟏𝟓.𝟓𝟓𝟔 𝟎.𝟎𝟏𝟏𝟕𝟔∗𝟏𝟓𝟒
Turbulent flow, Iterate through the following function to obtain f using the goal seek function in excel or
use moody’s chart
Colebrook-White equation
1 𝑒 2.51
= −2 log10 ( + )
√𝑓 3.7𝐷 𝑅√𝑓
f = 0.0206
0.5
11.4946∗ 10−4 𝑇𝑏 𝑃12 −𝑒 𝑠 𝑃22
𝑄= ∗ ( ) ∗ 𝐷2.5 , where s = 0 as no elevation change assumed
24√𝑓 𝑃𝑏 𝐺 𝑇 𝑓 𝐿𝑒 𝑍
0.5
11.4946 ∗ 10−4 273 + 15.56 3.5182 − 𝑃22
424.1 = ∗ ( ) ∗ 1542.5
24√0.0206 101.325 0.698 ∗ (273.15) ∗ 20000 ∗ 0.992
Line sizing is an essential aspect of designing chemical plants. It involves determining the appropriate
diameter of pipes to transport fluids (gases, liquids, or multiphase mixtures) within a chemical plant.
Proper line sizing ensures that the plant operates efficiently, safely, and with minimal maintenance
while meeting the desired process parameters.
1. Pressure drop: A key consideration is minimizing pressure drop across the pipeline. High
pressure drops can lead to insufficient flow rates, increased energy consumption, and reduced
plant efficiency. Engineers use various correlations and models to calculate pressure drops and
select appropriate pipe diameters.
2. Flow regime: The flow characteristics (laminar, turbulent, or transitional) depend on the fluid's
properties, flow rate, and pipe diameter. Engineers should ensure that the desired flow regime
is maintained to avoid issues such as poor heat and mass transfer, increased pressure drop, or
flow instabilities.
3. Velocity limits: Fluid velocities in pipes should be maintained within specific limits to avoid
issues such as erosion, noise, and vibration. Higher velocities can cause erosion and damage to
pipes, while low velocities can lead to sedimentation or other issues.
4. Economic factors: The cost of piping and installation should also be considered when sizing
pipelines. Larger diameter pipes are more expensive and may require more support structures,
while smaller pipes may be cheaper but can result in higher pressure drops and energy
consumption.
1. Accuracy of models and correlations: The models and correlations used for pressure drop
calculations and flow regime predictions are often based on empirical data and may not
accurately represent the specific conditions of a plant. Engineers must exercise judgment when
applying these tools and may need to consult experimental data or conduct simulations to
validate their assumptions.
2. Complexity of multiphase flows: In cases where multiple phases (gas, liquid, solid) are
present, predicting flow behavior and pressure drops can be more challenging. Engineers may
need to rely on experimental data, simulations, or specialized models for multiphase flow
systems.
3. Changes in operating conditions: Plant operating conditions can change over time due to
factors such as equipment wear, fouling, or changes in feedstock properties. These changes can
affect the flow behavior and pressure drops in pipelines, potentially requiring adjustments to
pipe diameters or other design parameters.
In summary, line sizing in chemical plants involves balancing various factors such as pressure drop, flow
regime, velocity limits, and economic considerations. Engineers must also consider the limitations of the
models and correlations used for line sizing, and be prepared to adapt to changing operating conditions.
Proper line sizing ensures that chemical plants operate safely, efficiently, and with minimal
maintenance. The basic principle to determine the line sizes shall be based on economical aspects, i.e.,
minimizing the sum of operational costs and investment costs. However, line sizes should not
exceed the limitations. In some instances, the process requirements will take precedence over
the economical aspects; for example, in the case of pump suction lines where the NPSH is the
main concern. [5]
References
1. GPSA Engineering Data Book Section 23 (Physical properties) and section 17 (fluid flow and
piping).
2. Computer methods for chemical engineering, Second Edition, Chapter 2: Fluid flow in pipes.
3. Piping Calculations manual, E. Shashi Menon. Chapter 7 Gas piping system (equations Q units
were converted to stdm3/hr instead of stdm3/day).
4. Gas Hydraulics, E. Shashi Menon. Chapter 2 pressure drop due to friction.
5. JGC Hydraulics standard JGS 210-120-1-37E
6. Petroleum refinery Engineering, W. L. Nelson.
7. American national standard for rotodynamic pumps for pump piping ANSI/HI 9.6.6 -2009
8. API RP 14E fifth edition, recommended practice for design and installation of offshore
production platform piping systems
1. JGC Hydraulics standard JGS 210-120-1-37E (focuses on calculations and the economic aspect of
line sizing, presents shortcut calculations and charts for quick calculations and line size
estimation)
2. Elements of chemical process engineering Chapter 1 P. 49: 1.10 Hydraulic analysis of a system (If
you are interested in a complete pumping system hydraulic analysis and you don’t mind US
Customary units)
3. Line sizing guidelines for worley parsons PTD-DGS-130 (thorough guidelines for various cases in
line sizing)
4. Chemical process engineering by A. Kayode Coker (present the manual calculations of pressure
drop in different cases and compares results with commercial software. Overall the an excellent
book, worth checking out and should be added to your library)
5. Chemical process design and simulation by Juma Hayday. Chapter 4: Pressure changing
equipment. (Basic calculations and correlations used in commercial software along with aspen
plus and aspen hysys case studies)
Pump hydraulic system analysis and specifications
Link: https://processpocket.streamlit.app/
This tool was developed for process engineers to quickly estimate pump system analysis on-site
or check your pump specs. As a part of a larger project to develop what is similar to Carl
Branan’s book “process engineers Pocket Handbook” these tools would allow a process
engineer to quickly calculate/estimate equipment efficiencies or sizing using standardized
calculations
The aim here is to take little-known data from the field (flow, pressures, temperatures,
compositions..etc.) and use it as input for a rough estimation without having to return to the
office to use commercial software or calculations Excel sheets to validate or to calculate.
Additionally, these tools may also serve as a gathered data validation tool.
Nomenclature
Input Required
Output Obtained
Corrected BHP kW
Selected motor/motor η kW / %
• Use Darcy equation (check previous calculations documentation for liquid pressure
drop)
𝑓 ∗ 𝜌 ∗ 𝑢2 ∗ 𝐿
∆𝑃𝑓 = 1.01972 ∗ 10−5 ∗
2∗𝑑
Where: ∆Pf : pressure drop (kg/cm2)
Reynold’s number
ρvd
Re =
μ
μ: Viscosity (kg/m.s) or 0.001*cP
Laminar flow
64
f=
R
Turbulent flow
Colebrook-White equation
1 e 2.51
= −2 log10 ( + )
√f 3.7D Re√f
Where: f: Darcy friction factor
D: pipe inside diameter, in (or mm)
e: absolute pipe roughness, in (or mm) (0.045-0.05 for commercial steel)
• Convert ∆P to head loss
𝑛 10
ℎ𝑡𝑜𝑡𝑎𝑙 (𝑚) = ∑ ∆P𝑖 ∗
𝑖=0 𝑆. 𝐺
Fitting K1 Kꚙ
Elbow, 90°, Standard (R/D = 1), Screwed 800.0 0.40
Elbow, 90°, Standard (R/D = 1), Flanged/welded 800.0 0.25
Elbow, 90°, Long-radius (R/D = 1.5), All types 800.0 0.20
Elbow, 90°, Mitered (R/D = 1.5), 1 weld (90° angle) 1000.0 1.15
Elbow, 90°, Mitered (R/D = 1.5), 2 weld (45° angle) 800.0 0.35
Elbow, 90°, Mitered (R/D = 1.5), 3 weld (30° angle) 800.0 0.30
Elbow, 90°, Mitered (R/D = 1.5), 4 weld (22.5° angle) 800.0 0.27
Elbow, 90°, Mitered (R/D = 1.5), 5 weld (18° angle) 800.0 0.25
Elbow, 45°, Standard (R/D = 1), All types 500.0 0.20
Elbow, 45°, Long-radius (R/D 1.5), All types 500.0 0.15
Elbow, 45°, Mitered (R/D=1.5), 1 weld (45° angle) 500.0 0.25
Elbow, 45°, Mitered (R/D=1.5), 2 weld (22.5° angle) 500.0 0.15
Elbow, 45°, Standard (R/D = 1), Screwed 1000.0 0.70
Elbow, 180°, Standard (R/D = 1), Flanged/welded 1000.0 0.35
Elbow, 180°, Long-radius (R/D = 1.5), All types 1000.0 0.30
Elbow, Used as, Standard, Screwed 500.0 0.70
Elbow, Elbow, Long-radius, Screwed 800.0 0.40
Elbow, Elbow, Standard, Flanged/welded 800.0 0.80
Elbow, Elbow, Stub-in type branch 1000.0 1.00
Tee, Run, Screwed 200.0 0.10
Tee, Through, Flanged or welded 150.0 0.05
Tee, Tee, Stub-in type branch 100.0 0.00
Valve, Gate, Full line size, Beta = 1 300.0 0.10
Valve, Ball, Reduced trim, Beta = 0.9 500.0 0.15
Valve, Plug, Reduced trim, Beta = 0.8 1000.0 0.25
Valve, Globe, Standard 1500.0 4.00
Valve, Globe, Angle or Y-type 1000.0 2.00
Valve, Diaphragm, Dam type 1000.0 2.00
Valve, Butterfly, 800.0 0.25
Valve, Check, Lift 2000.0 10.00
Valve, Check, Swing 1500.0 1.50
Valve, Check, Tilting-disc 1000.0 0.50
• Convert K to head loss
𝑛 𝑣𝑖2
ℎ𝑡𝑜𝑡𝑎𝑙 (𝑚) = ∑ 𝐾𝑖 ∗
𝑖=0 2𝑔
For control valves
• Given the control valve diameter and if the control valve Cv isn’t available, use Fisher’s
table for approximate Cv estimation[3]:
−1/2 𝑞
∑ 𝐾 𝐶𝑣 2 𝐶𝑣 =
𝐹𝑃 = [1 + ( ) ]
𝑁2 𝑑2 ∆𝑃𝑠𝑖𝑧𝑖𝑛𝑔
𝑁1 𝐹𝑃 √
𝜌1 /𝜌𝑜
“the sum of the following pressure drops at maximum flow may be used for this purpose:
1. 20% of the friction drop in the circuit (excluding the valve). (A circuit generally includes
all equipment between the discharge of the pump, compressor or vessel and the next
point downstream of which pressure is controlled. In most cases this latter point is a
vessel)
2. 10% of the static pressure of the vessel into which the circuit discharges up to the
pressure of 200 psig, 20 psig from 200 psig to 400 psig, and 5% above 400 psig.”
𝑆. 𝐺
𝐶𝑣 = 𝑄𝐿 √
△𝑃
Where:
QL : liquid flow in gpm
△ 𝑃 : pressure drop across the valve in psi
“In order to ensure that the valve is in a controlling position at maximum flow rate, the valve Cv
is the maximum process value determined above divided by 0.8. The reasons for using this
factor are that:
1. It is not desirable to have the valve fully open at maximum flow since it’s not then in a
controlling position.
2. The valves supplied by a single manufacturer often vary as much as 10-20% in Cv.
3. Allowance must be made for pressure drop, flow rate, etc, values which differ from
design.”
After estimating valve Cv, valve sizing can be estimated using the following equations:
𝑣𝑎𝑙𝑣𝑒 𝐶𝑣 1 Single-seated valves
𝑆 (𝑖𝑛𝑐ℎ𝑒𝑠) = ( )2
9
𝑣𝑎𝑙𝑣𝑒 𝐶𝑣 1 Double-seated valves
𝑆 (𝑖𝑛𝑐ℎ𝑒𝑠) = ( )2
12
𝑣𝑎𝑙𝑣𝑒 𝐶𝑣 1 Butterfly valve sizes
𝑆 (𝑖𝑛𝑐ℎ𝑒𝑠) = ( )2
20
Opening percentage = Cv desired / Cv max. for the valve
η= 80 – 0.2855F+3.78*10-4-FG-2.38*10-7FG2 + 5.39*10-4*F2-6.39*10-7*F2G+4*10-10*F2G2
Where:
η: pump percentage efficiency
F: developed head (ft)
G: flow (GPM)
Range of applicability:
F: 50 – 300 ft
G: 100 – 1000 GPM
For G: 25 – 99:
η= η (at G=100 GPM) – 0.35% * G * (100-G)
BHP = hydraulic power / ηestimated
Drawing the Hydraulic Gradient [6]
Effect of elevation
There are situations where the ground elevation may have drastic peaks and valleys that require
careful consideration of the pipeline topography. In some instances, the total pressure required
to transport a given volume of liquid through a long pipeline may depend more on the ground
elevation profile than on the actual frictional pressure drop. If we calculated the total pressure
required for a flat pipeline as 825 psi and an uphill pipeline to be 1002 psi. In the uphill case the
static elevation difference contributed to 17 percent of the total pressure required. Thus the
frictional component was much higher than the elevation component.
When there are drastic elevation differences in a long pipeline, sometimes the last section of
the pipeline toward the delivery terminus may operate in an open-channel flow. This means that
the pipeline section will not be full of liquid and there will be a vapor space above the liquid.
Such situations are acceptable in ordinary petroleum liquid (gasoline, diesel, and crude oil)
pipelines compared to high vapor pressure liquids such as liquefied petroleum gas (LPG) and
propane. To prevent such open-channel flow or slack line conditions, we must pack the line by
providing adequate back pressure at the delivery terminus as illustrated in Fig. 6.15.
Hydraulic Gradient
The graphical representation of the pressures along the pipeline, as shown in Fig. 6.16, is the
hydraulic gradient. Since elevation is measured in feet, the pipeline pressures are converted to
feet of head of liquid and plotted against the distance along the pipeline superimposed on the
elevation profile. If we assume a beginning elevation of 100 ft, a delivery terminus elevation of
500 ft, a total pressure of 1000 psi required at the beginning, and a delivery pressure of 25 psi at
the terminus, we can plot the hydraulic pressure gradient graphically by the following method.
The first example of a branched pipe system is a circulation system. In such a case Hstat = 0. The
volume flow Qp, which
passes through the pump, is divided at the branch point. Continuity dictates that
Qp=Qa+Qb
Both branches A and B have their own system curves, the sum of which determines the
resultant system curve. The intersection point of the resultant system curve and the pump curve
determines the pump’s optimum operating point.
Figure 9.2a also shows how large a portion of the pump flow
near to the pump i.e. that flow losses between the pump and
the losses in the main pipeline up to the branch point. This is Figure 1: h is the difference between Pb and Pa in meters
done by reducing the pump curve by the magnitude of the losses up to the branch point. The
“reduced” curve is then matched to the remaining system curve as before.
In the third example the pump curve is first reduced to the branch point. Thereafter the system
curves A+B are added to the resultant system curve Hsyst A+B. The intersection point between
the reduced pump curve and the resultant system curve determines the pump operational
point, as shown in figure 9.2c.
In the fourth example there is a level difference (in this example a supply head) between the
supply vessel on the suction side and the branch point. The system curve for the supply pipeline
exhibits therefore, a static head (Hstat < 0).
The method is the same as in the previous examples. First reduce the pump curve with the
supply line curve up to the branch point. Then determine the system curves for the branches
A+B from the branch point. The resultant branch curve is then matched to the pump's reduced
curve as shown in figure 9.2d.
This method can be applied in principle to calculate any pipework system, however complicated
and containing any number of branch points. The first step is to determine the resultant system
curve for the branch furthest from the pump in relation to the furthest branch point. The next
step is to deal with the next branch point and so on until the branch point nearest to the pump
is reached.
A common problem is to dimension a piping network with many branches and many tapping
points in such a way as to guarantee a predetermined capacity at each tapping point.
Dimensioning is carried out so that the pipeline requiring the most pressure, usually the pipeline
to the furthest tapping point, is dimensioned first. The branch lines are then calculated on the
basis of the excess pressure available at the branch points. The excess pressure at the branch
points can be such that it is necessary to equip a branch line with extra throttling to prevent it
from "stealing" capacity from the other tapping points.
The first step toward the goal of selecting the proper pump for a service involves four areas of
analysis:
Pump boundary conditions must be defined clearly and with as much certainty as is possible. A
comprehensive knowledge of all boundary parameters is essential to correct user specifications
and vendor bid preparation.
• Does the calculated value for NPSHa allow for increases in system resistance in the
suction piping due to fouling?
• Where suction is from a vessel, was the minimum possible operating level used in
calculating NPSHa? The minimum level at which maximum flow is expected should
provide the limiting guidelines here.
• What margin between NPSHa and NPSHr is considered as acceptable?
• Has the maximum expected flow been used to set NPSHa?
Min./Normal/rated flow should be determined. The rated flow should reflect the maximum flow
the system can envisage under current consideration, but must also consider the long range
outlook. Minimum flow requirements can conflict with rated requirements and recirculation
facilities may be required.
When determining these design flows, care must be taken to avoid an extremely conservative
approach. Higher than expected flow will result in a larger than required pump as in head
consideration. A more simple single stage, overhung pump application may require a double-
suction between bearings design under increased flow requirements.
While it’s of prime importance to define maximum and minimum flow requirements properly,
it’s also important to clarify the percentage of time at which the pump will operate at minimum,
normal and rated (or maximum) flows. Where a pump is used for two different services, the
lower flow may require excellent turndown while the higher flow will impose more stringent
NPSHa restrictions. Long term operation at lower flow can mean higher maintenance costs due
to higher bearing loads and shaft deflections.
the points of low flow cavitation and high temperature rise are
only valid considerations where extremely low flows are
considered probable for short periods which may cause severe
cavitation and eventual vapor locking of the pump. Such events
will lead to to mechanical seal failure and require protection
against even short duration of one or two minutes, where
volatile
liquids close to their vapor points being pumped.
1.3 Fluid specifications
System analysis must include a clear definition of the fluid to be
pumped and show all variations expected in fluid quality.
Fluid temperature variations must be defined. This will assist in
Specifying the NPSHa and must include heating and cooling
which may occur between the vessel or drum being drawn from
the
pump suction, in the case of vapor-liquid interface.
Higher temperatures will necessitate bearing housing cooling
and
require that mechanical seal arrangements be suitably designed.
Corrosion due to chemical due to chemical attack or oxidation
Must be considered. Compatibility of materials to resist
electrolytic reaction is important.
Erosion due to high percentage of particulate matter may
cause a premature performance decline. Large particles may
necessitate an open faced impeller. It may prove necessary to
specify wear plates or a hard coating to prolong life. Wear ring
flushing from an external source may be required.
Fluid toxicity may dictate the use of dual seals as may very high
temperatures, high flammability, and/or high vapor pressures.
Carcinogenic, strongly acid or strongly alkaline fluid impose
similar needs for more elaborate sealing and often required an
external clean fluid supply.
Entrained gases may cause cavitation and may have a very
negative effect on pump’s ability to produce the required
differential pressure.
1.4 Criticality of the service
Check if pump is spared or un-spared?
Loss of flow process debits (where a specific time can be tolerated for pump outage for
maintainability needs and reparability)
Safety Consideration (where fluid is toxic, highly flammable or carcinogenic special care may be
required in materials, mechanical seals and structural areas)
Continuous operation is normally viewed as the prime criteria for sparing, but intermittent
operation may also demand a spare pump. The definition of intermittent is important. A pump
may be required to operate one week in four, but 24-hout service during that may be crucial,
requiring a spare for high reliability of service.
Liquid specification In addition to the parameters outlined in API 610 datasheet,
comment must be made on solids content, toxicity & setup
temperature
Operating and site conditions Define min., normal and maximum flows. Remember to add a
table of time of which these flows are expected to operate in
Suction pressure min., normal and maximum
The maximum discharge pressure The max. pressure that will be encountered under conditions of
max. flow and min. suction pressure
Site conditions These would influence items such as electrical or steam tracing
requirements, lubricant quality, type of lubrication, motor
protection, etc..
Performance Important performance variables such as specific speed, suction
specific speed, NPSHr, MCF, POR, AOR.
API 610
“Rated flow shall be within the region of 80% to 110% of the best efficiency flowrate of the
pump as furnished.”
“the BEP flowrate for the pump as furnished should preferably be between the rated point and
the normal point”
For slurry pumps handling high concentrations of abrasive solids, the generally acceptable range
of rates of flow is primarily determined by wear considerations of the wear components. Refer
to ANSI/HI 12.1-12.6 Rotodynamic (Centrifugal) slurry pumps for nomenclature, definitions,
Applications, and Operation for information.
API 610
“pumps shall have a preferred operating region of 70% to 120% of the best efficiency flowrate
of the pump as furnished.” Rated flow shall be within the region of 80% to 110% of the best
efficiency flowrate of the pump as furnished.
AOR depends on a large number of factors, some of which are application-specific. This
discussion is limited to thos factors related to operating rates of flow and pump design. For the
purpose of this guideline it’s assumed that the pumped liquid is a Newtonian pure liquid with no
vapor, gas, suspended solids or abrasives.
Certain fluid mixtures may have other characteristics that affect the AOR. for example, the
minimum rate of flow when pumping a fluid that contains entrained gasses may be limited bt
gas accumulation at the eye of the impeller. These considerations are outside the scope of this
guideline;
Factors such as proper equipment selection, installation, maintenance and operation which also
greatly affect pump reliability are not factors in defining the AOR.
Bearing life will be reduced and noise, vibration and components stresses will be increased
when a pump is operated outside its POR. As a result, service life of a pump operated within the
AOR may be lower than within the POR. While the predicted bearing life will vary significantly
over the AOR. It’s recommended that the calculated bearing life be a minimum of two years of
operation in this range (basic rating life L10 equal to or greater than 17500 hours)
AOR: MCF in ANSI table below () to 120% of the flow at the best efficiency point (BEP)
AOR: MCF in ANSI table below () to 120% of the flow at the best efficiency point (BEP)
Determine the normally attainable efficiency of an API-type, single-stage, end suction process
pump driven by a 50 Hz motor when pumping clear water at 30 oC given Q equals 700 m3/hr or
0.19 m3/s. H equals 130 m, NPSHa equals 7.5 m, and the customer has specified a maximum
pump rotational speed of 1500 rpm.
1470 ∗ 0.190.5
𝑛𝑠 = = 17
1300.75
Because this calculated value is already below the specific speed for optimum efficiency, no
consideration will be given to selections at speeds below 1500 rpm.
From chart, at 0.19 m3/s and curve F, efficiency corresponding to the optimum specific speed for
API end suction process type pumps = 85.1%
From chart (fig. 20.e) at 0.19 m3/s, the normal deviation is ±3%
After calculating BHP estimated for your pump at max Q multiply it by one of the following
factors mentioned in API 610
Kw HP
Estimate motor efficiency from the following table and Calculate rated motor output
5.1 Minimum flow rate at max. impeller to avoid suction circulation (ANSI/HI 9.6.3-2012)
After calculating suction specific speed & estimating BEP obtain result from the following curve
1. Check Approximate Hydraulic coverage table from the ASME standard and select the
appropriate pump
2. Check MCSF percentage for your selected pump in MCF table using affinity laws
𝑄𝐵𝐸𝑃_𝑛𝑒𝑤
𝑀𝐶𝐹% 𝑎𝑡 𝑟𝑎𝑡𝑒𝑑 ℎ𝑒𝑎𝑑 = 𝑀𝐶𝐹% ∗
𝑄𝐵𝐸𝑃 𝑏𝑒𝑓𝑜𝑟𝑒 𝑡𝑟𝑖𝑚𝑚𝑖𝑛𝑔
Note: This tool trims the impeller diameter to achieve the rated head, QBEP here is QBEP
for the pump after adjusting pump head for rated head not Q rated in specs input
Vertical in-line centrifugal pumps (ASME B73.2-2016)
6. Predicted/required NPSHr
NPSHr Desired is either NPSHa/NPSH margin ratio or NPSHa – NSPH margin (or whichever
greater)
As per API 610, if the parallel operation is specified and the pumps are not individually flow
controlled, Head rise from rated point to shutoff shall be at least 10%.
Generally if specified, pump head curve shall be continuously rising from rated point to shutoff.
HI (ANSI/HI 9.6.1-2012)
Application Margins
Chemical process pumps SSS < 210 metric or 11000 US units SSS >= 210 metric
Min. margin ratio of 1.1 or a margin of 0.6 m Pumps operating within POR
Electric power plant pumps Power (kW/Stage) POR NPSH Ratio AOR NPSH ratio/Margin
Boiler feed > 225 and < 500 1.2 (ratio) 1.5 (ratio)
Nuclear power plant pumps Power (kW/Stage) POR NPSH Ratio AOR NPSH ratio/Margin
(1.0 m minimum ) (1.5 m minimum )
Waste water (cast iron impeller) <45 1.1 (ratio) 1.2 (ratio)
Waste water (stainless steel impeller) <45 1.05 (ratio) 1.1 (ratio)
Waste water (cast iron impeller) > 45 1.2 (ratio) 1.3 (ratio)
Waste water (stainless steel impeller) > 45 1.1 (ratio) 1.2 (ratio)
Pulp and paper stock pumps SSS < 145 metric or 7500 US units SSS > 145 metric
Min. margin ratio of 1.1 throughout AOR Min. margin ratio of 1.2 or a
margin of 1.0 m
Building Services SSS < 145 metric or 7500 US units SSS > 145 metric
Min. margin ratio of 1.0 throughout AOR Min. margin ratio of 1.1 or a
margin of 0.6 m
General industrial pumping applications Pumps operating within POR Pumps operating within AOR
Min. margin ratio of 1.05 or a margin of 0.6 Min. margin ratio of 1.1 or a
m margin of 1.0 m
An operating NPSH margin is necessary to ensure satisfactory operation. A minimum margin of 3 ft (0.9m) or a margin ratio of 1.2 (whichever
yields a higher NPSH requirement) should be made available. This margin should be increased if variables exist that will increase the NPSHR of
the pump
JGC NPSHr Charts[9]
𝑓 ∗ 𝜌 ∗ 𝑢2 ∗ 𝐿
−5
∆𝑃𝑓 = 1.01972 ∗ 10 ∗
2∗𝑑
0.01607 ∗ 900 ∗ 1.432 ∗ 25.908
= 1.01972 ∗ 10−5 ∗ = 0.0153
2 ∗ 0.254
Psuction = Pvessel - ∆𝑃𝑓 +Helevation = 0.0532 – 0.0153 + 1.2375 = 1.2754 kg/cm2.a = 0.24 kg/cm2.g
NPSHa = Psuction – Pvapor = 14.14 – 0 = 14.14 m
𝐷2 (25.4∗0.001∗6)2
𝐴=𝜋 = 3.14 ∗ = 0.018232
4 4
ρvd 900∗5.68∗0.0254∗4
Re = μ
= 1.08∗0.001
=534294.2
Using moody’s Chart, At Re = 534294.2, f = 0.0168
𝑓 ∗ 𝜌 ∗ 𝑢2 ∗ 𝐿
−5
∆𝑃𝑓 = 1.01972 ∗ 10 ∗ =
2∗𝑑
0.0168 ∗ 900 ∗ 2.5242 ∗ 36.82
= 1.01972 ∗ 10−5 ∗ = 0.125
2 ∗ 0.1524
Pressure at tee point
△P total= Pfittings + ∆𝑃𝑓 +∆𝑃𝑐𝑜𝑛𝑡𝑟𝑜𝑙 𝑣𝑎𝑙𝑣𝑒 + Helevation + ∆𝑃𝑚𝑖𝑠𝑐 =
0.033+0.125+2.412+(13.72*0.9/10) + (8.22*0.9/10)= 4.55 Kg/cm2
substituting for
F = 206 ft
G = 1153.7 gpm
Gives, η = 59.95 %
Power Calculations
𝑄∗𝐻∗𝑆𝐺 262∗62.88∗0.9
Hydraulic power = 367.46 = 367.46 = 40.35 Kw
Rated power = Hydraulic power/pump efficiency = 40.35/0.5995 = 67.306 Kw
How
You can remove any entry simply by selecting it and clicking “remove selected entry”
5. Check your calculations output or download the Calculations sheet
E
6. Go to the “Pump specs Checking” Tab and fill out the table below or click on “Use
pump system data” to auto fill
7. Check your output table
Note 1: you can skip the first tab and directly use the pump specs checking tab if you already have your pump
specifications
References
This tool was developed for the APRCO process team to quickly estimate compressor efficiency and its
anticipated power consumption. As a part of a larger project to develop what is similar to Carl Branan’s
book “process engineers Pocket handbook” these tools would allow a process engineer to quickly
calculate/estimate equipment efficiencies or sizing using standardized calculations
The aim here is to take little-known data from the field (flow, pressures, temperatures,
compositions..etc.) and use it as input for a rough estimation without having to return to the office to
use commercial software or calculations excel sheets to validate or to calculate. Additionally, These
tools may also serve as a gathered data validation tool.
Nomenclature
Input Required
Output Obtained
Polytropic exponent
1
𝑛=
ln 𝑇𝑑/𝑇𝑠
1−
ln 𝑃𝑑/𝑃𝑠
Adiabatic Temperature [1]
𝑘
𝑃𝑑 𝑘−1
𝑇𝑎𝑑𝑖𝑎𝑏𝑎𝑡𝑖𝑐 = 𝑇𝑠𝑢𝑐𝑡𝑖𝑜𝑛 𝑥 ( )
𝑃𝑠
Adiabatic Efficiency [1]
𝑘
𝑃𝑑 𝑘−1
𝑇𝑠𝑢𝑐𝑡𝑖𝑜𝑛 𝑥 (( 𝑃𝑠 ) − 1)
ƞ𝑎𝑑 =
𝑇𝐷𝑖𝑠𝑐ℎ𝑎𝑟𝑔𝑒 − 𝑇𝑠𝑢𝑐𝑡𝑖𝑜𝑛
𝑀. 𝑤𝑡 = ∑ 𝑦𝑖 ∗ 𝑀. 𝑤𝑡𝑖
𝑖=1
K (Cp/Cv) [1]
𝐶𝑝 𝑀𝐶𝑝
𝑘= = , 𝑅 = 1.986 𝐵𝑡𝑢/(𝑙𝑏𝑚𝑜𝑙 ˚𝐹)
𝐶𝑣 𝑀𝐶𝑝 − 𝑅
∑𝑛𝑖=1 𝑦𝑖 ∗ 𝑀𝐶𝑝
𝑘=
(∑𝑛𝑖=1 𝑦𝑖 ∗ 𝑀𝐶𝑝) − 1.986
MCp values at different temperatures for a variety of compounds are found in figure 1 [1]
Compressibility factor Z [3]
𝑇𝑐 = ∑ 𝑦𝑖 ∗ 𝑇𝑐𝑖
𝑖=1
𝑛
𝑃𝑐 = ∑ 𝑦𝑖 ∗ 𝑃𝑐𝑖
𝑖=1
𝑇 𝑃
𝑇𝑟 = , Pr =
𝑇𝑐 𝑃𝑐
Where Tci and Pci are each component’s critical temperature and pressure
Naturally, The Standing and Katz method is applied where Z is obtained as a function of Tr and Pr through the
famous chart (figure no. ). Alternatively, Dranchuk, Purvis, and Robinson method is used. In this method the
Benedict-Webb-Rubin equation of state is used to correlate the Standing-Katz Z
factor chart. Eight coefficients A1, A2, etc., are used in this equation as shown [1]:
𝑨𝟐 𝑨 𝟑 𝑨𝟓 𝑨𝟓 𝑨𝟔 𝝆𝒓 𝟓 𝑨𝟕 𝝆𝒓 𝟑
𝒁 = 𝟏 + (𝑨𝟏 + + ) 𝝆𝒓 + (𝑨𝟒 + ) 𝝆𝒓 𝟐 + +
𝑻𝒑𝒓 𝑻𝟑 𝑻𝒑𝒓 𝑻𝒑𝒓 𝑻𝒑𝒓 𝟑 (𝟏 + 𝑨𝟖 𝝆𝒓 𝟐 )𝒆𝒙𝒑(−𝑨𝟖 𝝆𝒓 𝟐 )
Where:
𝟎. 𝟐𝟕 ∗ 𝑷𝝆𝒓
𝝆𝒓 =
𝒁𝑻𝝆𝒓
Calculations Examples using commercial software, hand written calculations and the tool developed:
Example No.1: Unknown composition, Known physical properties and operating conditions
Hand Calculations
Computer Methods for chemical engineers (2nd edition Page no. 78) (input highlighted)
𝑨𝟐 𝑨 𝟑 𝑨𝟓 𝑨𝟓 𝑨𝟔 𝝆𝒓 𝟓 𝑨𝟕 𝝆𝒓 𝟑
𝒁 = 𝟏 + (𝑨𝟏 + + ) 𝝆𝒓 + (𝑨𝟒 + ) 𝝆𝒓 𝟐 + +
𝑻𝒑𝒓 𝑻𝟑 𝑻𝒑𝒓 𝑻𝒑𝒓 𝑻𝒑𝒓 𝟑 (𝟏 + 𝑨𝟖 𝝆𝒓 𝟐 )𝒆𝒙𝒑(−𝑨𝟖 𝝆𝒓 𝟐 )
Where:
𝟎. 𝟐𝟕 ∗ 𝑷𝝆𝒓
𝝆𝒓 =
𝒁𝑻𝝆𝒓
For both suction and discharge Tr and Pr (obtained above) then average your result
(Use the Goal seek function in an Excel Sheet, Change your assumed Z till ABS(Z_assumed –
Z_calculated)*1000 = 0 for instance)
Z = 0.993
𝑘
1 𝑘 𝑍𝑅𝑇𝑠𝑢𝑐𝑡𝑖𝑜𝑛 𝑃𝑑 𝑘−1
𝐵𝐻𝑃 = ∗𝐺∗ ∗( ) ∗ (( ) − 1)
(ƞ𝑎𝑑 ∗ 36) 𝑘−1 𝑀. 𝑤𝑡 𝑃𝑠
1 1.26 0.993𝑥8.3143𝑥303.15 1.26
= 𝑥 100 𝑥 𝑥 𝑥 ((3.33)0.26 − 1) = 6.62 𝐾𝑤
74.8 𝑥 36 0.26 19.44
Tool input and output
References:
1. JGC Compressors standard JGS 210-120-1-28E (Calculations using Schultz method and PROII)
2. Elements of chemical process engineering Chapter 5: Compressors (If you are interested in
specifying compressors and you don’t mind US Customary units)
3. PTC-10: Performance test code on compressors and exhausters (A more detailed calculations,
examples and test parameters using the Schultz method that the commercial software are
using)
Shell and Tube Heat Exchanger Calculations
This tool was developed for process engineers to quickly estimate heat exchangers (single-phase gas or
liquid) efficiencies and create preliminary designs using kern, bell Delaware and NTU methods, and
access a list of useful tables on site. As a part of a larger project to develop what is similar to Carl
Branan’s book “process engineers Pocket Handbook” these tools would allow a process engineer to
quickly calculate/estimate equipment efficiencies or sizing using standardized calculations
The aim here is to take little-known data from the field (flow, pressures, temperatures,
compositions..etc.) and use it as input for a rough estimation without having to return to the office to
use commercial software or calculations Excel sheets to validate or to calculate. Additionally, these tools
may also serve as a gathered data validation tool.
Nomenclature
Kern & Bell Delaware Rating Calculations Kern & Bell Delaware Preliminary Calculations
Flow rates kg/hr Flow rates kg/hr
o o
Inlet/outlet temperatures C Inlet/outlet temperatures C
Inlet/outlet pressures Kg/cm2.a Inlet/outlet pressures Kg/cm2.a
Average Densities [1] Average Densities [1]
Average Heat Capacities [1] Kcal/kg. oC Average Heat Capacities [1] Kcal/kg. oC
Average Heat Thermal Conductivities[1] W/m. °C Average Heat Thermal W/m. °C
Conductivities[1]
Average Viscosities [1] cP Average Viscosities [1] cP
Alternative option: Alternative option:
Liquid/gas: Composition Liquid/gas: Composition Vol% or
Vol% or mol%
Liquid oil fraction: SG. & temperatures Liquid oil fraction: SG. & mol%
- /°C
temperatures - /°C
Fouling factors W/m. K Fouling factors W/m. K
Exchanger Geometry Desired features in the design
Tube Diameter & thickness mm Bell Delaware
Pitch type and pitch -/mm Pitch type and pitch ratio -/-
Shell diameter mm Kern
Number of tubes / number of passes Assumed Uservice W/m2. °C
Tube Length mm Length mm
Baffle spacing mm Pitch type and pitch ratio -/-
Baffle cut %
Output Obtained
Duty KCal/hr Q W
∆TLMTD UA W/K
Ft Cr W/K
Corrected LMTD Cmin W/K
Surface Area m² Cmax W/K
Tube Heat transfer Coef. W/m²-°C effectiveness
Shell Heat transfer Coef. W/m²-°C NTU
Uclean W/m²-°C Thi °C
Udirty W/m²-°C Tho °C
Uservice W/m²-°C Tci °C
Over Design % Tco °C
Over Surface %
Shell Reynolds Number
Tube Reynolds Number
Tube Velocity m/s
Shell Velocity m/s
Shell D mm
Baffle Spacing mm
Number of baffles
Do mm
Di mm
Length mm
Number of tubes
Number of passes
Tube pitch mm
pitch type
baffle cut %
Heat Exchangers calculations
Shell and tube heat exchanger rating is an essential process in evaluating the performance and
efficiency of a shell and tube heat exchanger. Heat exchanger rating involves determining the actual
heat transfer rate, pressure drop, and overall thermal performance of the heat exchanger based on
its operating conditions and design specifications.
During the rating process, several key factors are taken into account. These include the inlet and
outlet temperatures of both the hot and cold fluids, the flow rates of the fluids, the physical
properties of the fluids (such as density, specific heat capacity, thermal conductivity, and viscosity),
and the heat exchanger's geometric and design parameters (such as tube diameter, tube length,
shell diameter, and number of tube passes).
1. Heat Transfer Calculation: This calculation estimates the actual heat transfer rate occurring in
the heat exchanger.
2. Pressure Drop Calculation: This calculation determines the pressure drop experienced by the
fluids as they flow through the heat exchanger.
By performing heat transfer and pressure drop calculations, the rating process allows engineers to
evaluate the actual performance of the heat exchanger and compare it against the design
specifications.
Additionally, heat exchanger rating also provides valuable insights into the heat exchanger's thermal
performance over time, considering factors like fouling, scaling, or degradation. This information
aids in determining maintenance schedules, cleaning intervals, and potential efficiency
improvements.
Where
S: Or P, Temperature efficiency
Using Fig. 18 in reference [1]. We get the correction factor Ft or from [2]
√𝑅 + 1 ∗ ln( )
𝐹𝑡 =
( √ )
(𝑅 − 1) ∗ ln
( √ )
n: no. of units
of times the shell fluid passes the tube bundle. n number of units in series.
Tube side
𝜋𝑑 𝑁
𝐴 = ×
4 𝑁
Nt: no. of tubes
Np: no. of passes
Fluid mean velocity
𝑚̇
𝑈 =
𝜌 ×𝐴
ρ: fluid density kg/m3
𝑅 = , µ in N.s/m2
µ
𝑓 = (1.58 ln 𝑅 − 3.28)
Pr: Prandtl number Cpµ/k
𝑁𝑢 . 𝑘
ℎ =
𝑑
Tube side pressure drop [3]
4𝑓𝐿𝑁 ⍴𝑢
∆𝑝 = 𝑛. + 4𝑁
𝑑 2
n: number of units in series
1.2 Heat exchangers rating (Bell-Delaware Method) [note]
In Bell Delaware method, the fluid flow in the shell is divided into a number of individual
streams. Each of these streams introduces a correction factor which is used to correct heat
transfer coefficient and pressure drop across the shell. This article gives step-by-step guidance on
doing heat exchanger rating analysis based on Bell-Delware method.
Note:
This section is mostly the Che Guide article with minor clarifications Link
Cross flow area, Sm in m2 is the minimum flow area in one baffle space at the center of the tube
bundle. It is calculated by following equation:
DotL = Ds - (12.5+( Ds /200))
Sm = B[(Ds - DOTL) + (DOTL - Do)(PT - Do)/PT,eff ]*10-6
where, PT is tube pitch, B is central baffle spacing, DOTL is outer tube limit diameter, Ds is shell
diameter and Do is tube outside diameter. (all units in mm)
GS = 𝑚 /Sm ∗ 3600
where, mS is shell side mass flow rate. Shell side Reynolds number Re S is then calculated from
ReS = Do.GS / μS *10
where, μS is the shell side fluid dynamic viscosity at average bulk temperature in cP.
where, CP,S is the shell side fluid specific heat (kcal/kg.oC) and kS is the shell side fluid thermal
conductivity (W.m/k).
1.33
𝑗=𝑎 . (𝑅𝑒 )
𝑎
𝑎=
1 + 0.14(𝑅𝑒 )
where a1, a2, a3 and a4 are the correlation constants listed below.
The ideal tube bank based coefficient (W/m2.k) is calculated from –
µ .
ℎ = 1.163 ∗ 𝑗. 𝐶 , . 𝐺 . (𝑃𝑟 ) . ( )
µ,
The factor JC accounts for heat transfer in the baffle windows. It has a value of 1.0 for exchanger
with no tubes in the windows.
JC = 0.55 + 0.72FC
FC = 1 - 2FW
FW = (θCTL - Sin(θCTL))/2π
θCTL = 2cos-1(Ds(1 - 2*Bc/100)/DCTL)
DCTL = DOTL - Do
where, Bc is segemental baffle cut in %. Fc is Fraction of tubes in cross flow. Fw is Fraction of tubes in
windows flow. Dctl is central tube limit diameter in mm (the diameter of the circle that passes through the
centers of the outermost tubes in the bundle)
Figure 1: Baffle geometry used for calculating the fraction, Fw, of tubes
in one baffle window
a.2 Correction factor for Baffle Leakage, JL
The correction factor JL considers the effects of the tube-to-baffle and shell-to-baffle leakage
streams on heat transfer.
DSB = 3.1+0.004* Ds
JL = 0.44(1-rS) + (1-0.44(1-rS))exp(-2.2rL) LTB = 0.4 mm
𝑟 = Ssb /(Ssb + Stb) If the longest unsupported
𝑟 = (Ssb + Stb)/ Sm
Ssb = Ds*DSB(π - 0.5θDS)*10 tube length is < 3 ft,
Stb = (π/4)((Do+LTB)2 - Do2)Nt(1-FW)*10 otherwise: 0.8 mm
θDS = 2cos-1(1 - 2Bc/100)
Figure 2: Ssb Shell-to-baffle leakage
area.
Where, Nt is number of tubes, DSB is diametral clearance between shell & baffle and LTB is
diametral clearance between tube and baffle (mm). rS is ratio of shell-baffle to total area and rL is
ratio of leakage to cross flow. θDS is Baffle window angle, Ssb is Shell to baffle leakage area and
Stb is Tube to baffle leakage area. (m2)
Bundle bypass correction factor JB accounts for the bundle bypass stream flowing in the gap
between the outermost tubes and the shell. The number of effective rows crossed in one cross
flow section, Ntcc between the baffle tips is provided by following equation.
The bundle bypass flow area, Sb (m2) is defined as the area between the outermost tubes and the shell
at the shell centerline in one central baffle space. It is part of the cross-flow area , Sm.
Sb = B(Ds - DOTL - Do/2)*10
The factor JR accounts for the decrease in the heat transfer coefficient with downstream distance
in laminar flow.
Ntcw = (0.8/Pp)(Ds(Bc/100) - (Ds-(DOTL-Do))/2 )
NB = 1 + int((L - 2Ls - LBIn - LBOut)/B)
NC = (Ntcw + Ntcc)(1 + NB)
JRL = (10/NC)0.18
JR = 1, ReS > 100
JR = JRL + (20-ReS)(JRL - 1)/80, ReS <= 100, ReS > 20
JR = JRL, ReS <= 20
where, L is tube length, Ls is tubesheet thickness, LBIn is inlet baffle spacing and LBOut is outlet baffle
spacing. All units in mm. Ntcw is the number of tube rows crosses, NB is number of baffles and Nc is Tube
rows crossed in entire exchanger.
Ls ~ DS √P/C (mm)
a.5 Correction factor for unequal baffle spacing, J S P : Design pressure (kPa or
Kg/cm2G) of shell or tube
n1 = 0.6, ReS >= 100 side, whichever is greater.
n1 = 1/3, ReS < 100 C : Constant (578 for kPaG,
JS = ((NB-1)+(LBIn/B)1-n1 + (LBOut/B)1-n1)/((NB-1)+(LBIn/B) + (LBOut/B)) 58.3 for Kg/cm2G)
Or whichever greater than:
Shell side heat transfer coefficient is calculated as Ls = 50 mm Ds ≥ 500mm
Ls = 0.1* Ds for Ds < 500mm
hs = hIdeal(JC.JL.JB.JS.JR)
LBIn , LBOut is assumed = B
b. Shell Side Pressure Drop, ΔPs
1.33
𝑓=𝑏 . (𝑅𝑒 )
𝑏
𝑏=
1 + 0.14(𝑅𝑒 )
RL = exp(-1.33(1+rS)rLp)
p = 0.8 - 0.15(1+rS)
RS = 0.5((B/LBIn)2-n + (B/LBOut)2-n)
n = 0.2, ReS >= 100
n = 1.0, ReS < 100
where, Di is tube inside diameter mm, v is velocity m/s, ρ t is density, μt is viscosity in cP, kt is
thermal conductivity W/m.K and Cp,t is specific heat for fluid on tube side kcal/kg.oC.
For laminar flow, ReT < 2300, Sieder and Tate correlation is used for Nusselt's nubmer.
Nu = 1.86(ReT.PrT.Di/Leff)1/3
Leff = L - 2*Ls , where Leff is effective tube length in mm for fixed tube sheet, see [2] for others
For turbulent flow, ReT > 10,000, following equation developed by Petukhov-Kirillov can be
used.
Rtube = Do*0.001/(2.ln(Do/Di).ktube)
where, ktube is thermal conductivity of tube material W/m.k. Overall clean heat transfer
coefficient, UClean is calculated as per below equation
Overall dirty heat transfer coefficient, UDirty is calculated as per below expression
where, fshell & ftube are fouling factors for shell and tube side in m2.k/W.
where, Q is heat duty in W, A is heat transfer area and LMTD corrected is corrected logarithmic
mean temperature difference.
Accuracy of Correlations Between Kern’s Method and the Bell–Delaware’s Method [4]
Heat exchangers designs have been arrived at using the Kern’s, Tinker’s, or the Bell’s method. The
suitability and accuracy of these have been reviewed in the literature. Kern’s method cannot be applied
to a TEMA type T floating head heat exchanger without sealing strips or with unsealed pass partition
lanes. Whitley presented a study of the errors found in heat transfer coefficient and pressure drop
predictions obtained with the Kern and Bell methods. Palen and Taborek show that the Bell–Delaware
method allows the prediction of shell-side film coefficients in the range from 50% lower to 100% higher
than the real values. Table below shows a comparison of the Kern, Bell–Delaware, and the Tinker’s
methods.
Kern Bell-Delaware Tinker
1. Ease of use Simple More involved More involved
2. Resultant design Very conservative Not conservative Not conservative
3. Cost of equipment Very high Relatively low Relatively low
4. Does it account for
a. By pass and leakage streams No Yes Yes
b. Inlet and outlet baffle
spacing being different than
the central one No Yes Yes
c. Number of tube rows being
different in inlet and outlet zones
than in the center. No Yes No
d. Seal strips No Yes Yes
e. Different tube layouts and
baffle cuts No Yes Yes
f. Effectiveness of tube rows in
window by a separate
calculation. No Yes Yes
Original plots
extended.
Method remains Δp calculated lam Method remains
g. Laminar flow unchanged. differently. unchanged.
Accounts for it by
Assumes shell full number of tubes in Accounts for it by mean
h. Size of tube bundle of tubes cross flow bundle width.
i. Pressure drop in the nozzle. Yes No No
Yes (takes an
j. Δp due to gradual fouling of average
heat exchanger working unit) No Yes
5. Basis of Reynolds number Equivalent
calculation diameter Tube O.D. Tube O.D.
Cross-flow mass velocity
Geometric mean of multiplied by a factor
Cross flow over the cross-flow and for tubes in baffle
6. Which mass velocity used? tube bundle window flow. window.
Source: Petroleum Refining Design and applications, vol. 4, A. Kayode coker
The NTU method is based on the concept of heat transfer effectiveness, which represents the ratio of
the actual heat transfer in the heat exchanger to the maximum possible heat transfer. It takes into
account the temperature difference between the hot and cold fluids, the heat exchanger configuration,
and the fluid flow rates.
In the NTU method, the heat exchanger is typically represented as a series of heat transfer units.
Each heat transfer unit consists of a hot fluid stream and a cold fluid stream. The NTU value
represents the number of these heat transfer units in the heat exchanger
1. Determining the NTU Value: The NTU value is calculated based on the heat exchanger geometry,
flow rates, and heat transfer characteristics. It can be determined using equations or graphical
methods.
2. Calculating the Heat Transfer Effectiveness: The effectiveness ( ε) of the heat exchanger is
determined using the NTU value and the heat exchanger configuration. The effectiveness
represents the ratio of the actual heat transfer to the maximum possible heat transfer. It is a
measure of how efficiently heat is transferred between the hot and cold fluids.
3. Estimating the Heat Transfer Rate: Once the effectiveness is known, the heat transfer rate in the
heat exchanger can be calculated. This involves multiplying the maximum possible heat transfer
rate (which is determined by the temperature difference between the hot and cold fluids) by the
effectiveness.
4. Estimating outlet temperatures by simple heat balance in shell and tube
It should be noted that the NTU method assumes certain simplifications and may not account for all
complex factors, such as fouling or non-uniform fluid flow distribution. In such cases, more detailed
analysis and modeling may be required.
Capacity ratio
ṁ
Capacity ratio=
ṁ
With the right tools it may be possible to manually adjust exchanger outlet temperatures to come up
with acceptable heat balances across the preheat train system; however, this is not a trivial undertaking
due to the complexity of interactions between exchangers. The final solution is also subjective as it
involves a set of sequential decisions as to which exchanger to start with and whether to adjust the
tube- or shell-side temperature in order to achieve a heat balance. There are many possible heat
balance solutions. The question is: Which is the most representative solution?
Mathematical data reconciliation can effectively answer this question. It is a vital
component of any successful fouling monitoring application. Reconciliation is a least squares
optimization process where the objective function minimizes the deviation between reconciled, heat
balanced data and the raw plant data. Temperature, fluid flow, and flow split data can be included in the
data reconciliation model. The relative weight the reconciliation model places on different data types is
governed by assigning trust factors or typical errors to each plant measurement. Occasionally, though
not often, it may be necessary to include local exchanger bypasses in the reconciliation problem.
Given that local bypass flows are never measured, bypasses must be estimated and thus very low trust
factors (or high typical errors) should be assigned to these variables.
A measure of “acceptability” should be established for any optimized solution, to determine whether to
accept or reject the current data set. If reconciled plant data deviations on a particular plant
measurement are consistently higher than expected instrumentation errors, the flow meter or
temperature indicator should be checked.
To track the amount of fouling in an exchanger one can calculate Uo frequently (weekly) and plot Uo versus
time. The pattern of the fouling curve will usually repeat itself after each cleaning. 1/Uo = R= RCLEAN +RDIRT
If conditions have changed significantly from design condition, h i and ho can be adjusted to current
conditions using the equations from the Natural Gasoline Processing Handbook summarized in
Key Formulas. RCLEAN can then be recalculated. The current RDIRT can be compared to the design RDIRT.
Many heat exchange services are regulated by instrumentation to provide a stream to a downstream
process at a consistent temperature. This temperature is accomplished by bypassing part of the
exchanger feed around the exchanger and mixing this bypassed material with the exchanger effluent.
The degree of bypassing is regulated by a temperature controller which senses the downstream mixture
temperature. When the exchanger is clean, the maximum amount of feed is bypassed; and when the
exchanger is fouled to the maximum, no feed is bypassing the exchanger. By recording the clean
bypassing condition (percent of feed bypassed), you can easily estimate the amount of exchanger
fouling by observing the amount of feed bypassed with the partially fouled exchanger.
The following equations which can be readily derived from the
above and are often quoted in the literature:
Fraction dirty = rn/rd = (Ud / UN)[UC – Un ] / (UC – Ud )
where:
rd = RDIRT design
rn = RDIRT current (now)
UC = Clean U with no fouling factors
UD = Design U which includes fouling factors
Un = Current U (now)
3. Shell &Tube Preliminary sizing
The preliminary sizing stage focuses on establishing initial estimates and configurations before
proceeding to detailed design calculations. During preliminary sizing, several key parameters need to be
considered. These include the heat duty or the amount of heat to be transferred, the temperature
difference between the two fluids, the flow rates of the fluids, and the physical properties of the fluids,
such as density and specific heat capacity.
The goal of preliminary sizing is to determine the approximate size and layout of the shell and tube heat
exchanger, including the number of tubes, tube diameter, tube length, and shell diameter. This is
achieved by employing empirical correlations, design guidelines, and engineering experience to estimate
the required surface area for heat transfer.
Additionally, factors like pressure drop, fouling, material selection, and construction constraints are also
taken into account during the preliminary sizing phase.
The results obtained from preliminary sizing serve as a foundation for detailed design and analysis,
where more precise calculations and considerations are made to optimize the heat exchanger's
performance, efficiency, and cost.
In preliminary sizing of shell and tube heat exchangers, there are certain aspects and considerations that
are generally not accounted for. These may include:
1. Vibration Analysis
2. Detailed Pressure Drop Calculations.
3. Thermal Stress Analysis
While preliminary sizing may not incorporate these detailed analyses and considerations, it is important
to note that they are crucial in the subsequent stages of design and engineering to ensure the heat
exchanger's reliability, performance, and longevity.
A = area (m2)
d = tube diameter (m)
Ds = shell diameter (m)
Ft = ΔTLMTD correction factor (non-counter-current flow)
h = film heat transfer coefficient (W/m2 oK)
K = dimensional constant, solely dependent on physical properties, volumetric flow rate, and a
single
characteristic dimension.
R = fouling resistance (m2 K/W)
Δp = pressure drop (kPa)
Subscripts
i = inside surface
o = outside surface
s = shell side
t = tube side
Source 7: Rapid Design Algorithms for Shell-and-tube and Compact Heat Exchangers, G. T. Polley
Walkthrough
a. Initialization
a. Assume an exchanger with shell diameter Ds = 387 mm baffle cut = 25%, 6 m length, 2 passes, same
tube diameter and same pitch ratio and type as specified.
b. Estimate the number of tubes that could fit in your shell diameter
c. Estimate shell Reynold’s number
µ .
ℎ = 1.163 ∗ 𝑗. 𝐶 , . 𝐺 . (𝑃𝑟 ) . ( )
µ,
1.33
𝑓=𝑏 . (𝑅𝑒 )
𝑃
𝐷𝑜
𝑏
𝑏=
1 + 0.14(𝑅𝑒 )
RL = exp(-1.33(1+rS)rLp)
p = 0.8 - 0.15(1+rS)
b.2 Correction factor for Bundle Bypass effect, R B
RB = exp(-Cr(Sb / Sm)(1 - (2rss)1/3)) for rss < 0.5
RB = 1 for rss >= 0.5
Cr = 4.5 for ReS < 100
Cr = 3.7 for ReS >= 100
ReT = Di.v.ρt/μt*10
PrT = 1.163*Cp,t.μt/kt
For laminar flow, ReT < 2300, Sieder and Tate correlation is used for Nusselt's nubmer.
For transient flow, Nusselt number can be interpolated from Nu Laminar & Nu Turbulent.
ht = Nu.(kt/Di)(μt/μt, w)0.14
.
∆𝑃 = 𝑘 ∗ 𝐴 ∗ ℎ
.
∆𝑃 = 𝑘 ∗ 𝐴 ∗ ℎ
.
∆𝑃 = 𝑘 ∗ 𝐴 ∗ ℎ
.
∆𝑃 = 𝑘 ∗ 𝐴 ∗ ℎ
Where ∆𝑃 and ∆𝑃 are shell and tube allowable pressure drops
Solve the following equations to obtain tube Reynold’s number and friction factor
𝑣 = Di.Re.ρt/μt
at = mt/(⍴t*vt*3600)
*
Nt = at/ π. ∗
For laminar flow, ReT < 2300, Sieder and Tate correlation is used for Nusselt's nubmer.
Leff = L - 2*Ls , where Leff is effective tube length
Nu = 1.86(ReT.PrT.Di/Leff)1/3
For turbulent flow, ReT > 10,000, following equation developed by Petukhov-Kirillov can be
used.
Nu = (f/2)ReT.PrT/(1.07+12.7(f/2)0.5(PrT2/3-1))
f = (1.58 ln(ReT) - 3.28)-2
For transient flow, Nusselt number can be interpolated from Nu Laminar & Nu Turbulent.
ht = Nu.(kt/Di)(μt/μt, w)0.14
if ht = ht.solved stop iteration, else change Re accordingly and start over.
Now you’ve calculated Tube length, Number of baffles, Ret and Vt and tube count.
c. Iterate!
1. Estimate Shell diameter from the tube count. If it didn’t change proceed, if changed repeat from
initialization with a new shell diameter accordingly.
2. Calculate ∆𝑃 If it didn’t change from your allowable shell pressure drop proceed, if it changed repeat
from initialization with a new baffle cut accordingly.
3. If no baffle cut % achieved ∆𝑃 start from initialization with a new number of passes.
c. Calculate ho from
.
𝐷𝑒 𝐷𝑒. 𝐺𝑠 𝐶𝑝. µ µ .
ℎ ∗ = 0.36 . .( )
𝑘 µ 𝑘 µ
d. If ho appears too low, assume closer baffle spacing, up to 1/5 of the shell diameter and
recalculate Gs and ho. If this second trial is obviously too low, then a larger shell size
may be indicated; therefore, return to step 3, re-evaluate the assumed U to be certain
that is attainable.
11. If the ho appears to have possibilities of satisfying the design, continue to a conclusion by
assuming the tube-side and shell-side fouling and calculate U overall.
12. Compare values calculated in steps 10 and 5. If the calculated Uoverall is too small, re-assume a
new lower U for step 5 or try closer baffle spacing in step 9 but do not get baffles closer than 1/5
the shell I.D.
13. Calculate the percent of over Design. A reasonable figure is 10–20%.
14. Calculate the shell-side pressure drop. If p is too high, reassume U (step 3).
15. Calculate the tube-side pressure drop.
If the tube-side pressure drop exceeds a critical allowable value for the process system, then recheck by
either lowering the flow rate and changing the temperature levels or reassuming a unit with fewer
passes on tube-side or more tubes per pass. The unit must then be rechecked for the effect of changes
on heat transfer performance. The following figure illustrates the influence of various geometrical
parameters on heat exchanger heat transfer and pressure drop.
Need to Increase U
Increase F
[1] Use consistent units for any one variable in both cases
[2] F is the ratio of the new value to the old value for a given variable. The overall f is the
product of the individual fs.
[3] Number of rows of tubes exposed to cross flow (as opposed to parallel flow). This
Number is determined by baffle and bundle geometry.
[4] (hi/h2) = r2/ r1
.
ℎ 𝑟 𝜇 . 𝑘 . 𝐶 𝐺 . 𝐷 .
= = . . . .
ℎ 𝑟 𝜇 𝑘 𝐶 𝐺 𝐷
. . .
∆𝑃 𝜇 𝜌 𝑁 𝐺 𝐷
== . . . .
∆𝑃 𝜇 𝜌 𝑁 𝐺 𝐷
Note: Check for the lowest h in your calculations it should be the dominating factor
and the most factor that would influence your overall U.
Methods to improve the flow fraction
To increase “B” fraction, “A”, “C”, “E” and “F” fractions should be decreased by the following methods:
Method \ Fraction A C E F
Enlargement of baffle spacing Effective Effective
Increasing baffle cut Effective Effective
Decreasing tube-to-baffle clearance(*1) Effective
Installation of sealing strips Effective
Installation of seal rods Effective
Change of tube passes arrangement Effective
Change of baffle cut orientation Effective
NOTE(*1): Refer to TEMA RCB-4.2 with regard to standard its clearance
The recommended limitation about flow fractions are shown hereunder:
B stream : For turbulent flow; Min. 50%, For laminar flow; Min. 40%
C stream : Max. 10%
E stream : Max. 15%
F stream : Max 20%, preferably 15% or below
Guidelines to be kept in mind
The dynamic head, ρV2(*1), is a limitation to prevent or minimize erosion of the tube bundle,
and is used as a guideline for the necessity for impingement protection devices. Impingement
protection devices should be provided when ρV2 through the inlet nozzle exceeds the
following:
Non corrosive, non-abrasive, single-phase fluid: 2230 (kg/m-sec 2)
All other liquids, including a liquid at its bubble point: 740 (kg/m-sec 2)
Note(*1): ρ is the fluid density, and V is the linear velocity of fluid. For two-phase fluid, the
mean density should be calculated assuming a homogeneous vapor-liquid mixture.
The allowable baffle cut range for each baffle style is as follows:
Segmental: 10~49% (single phase: recommended 25%, 10-15 accepted)
Double segmental: 10~30%
Triple segmental: 10~40%
Segmental NTIW: 15~30%
Minimum spacing[2]
Segmental baffles normally should not be spaced closer than 1/5 of the shell inside diameter, or
50 mm, whichever is greater. This is the design practice specified by TEMA.
Max. spacing[2]
Baffle or support plates should be spaced such that the unsupported tube span does not exceed
the values
Table: TEMA Maximum Unsupported Tube Spans
NOTE(*1): If sea water is not treated with an effective iron ion, the maximum water
velocity should be less than the following conditions: 1.2m/sec for admiralty brass,
1.8m/sec for Aluminum brass. Engineers should confirm the sea water supply
conditions.
(c) To reduce the tendency to foul, the maximum allowable outlet temperature for
cooling water is limited as follows:
SW: 50°C
CTW: 60°C
(d) Another important criterion is the film temperature on the tube surface including a
fouling layer. If the film temperature exceeds the following conditions, severe fouling
could result.
SW: 65°C
CTW: 70°C
Slurry handling[3]
25.0 mm O.D. or larger tubes should be used to prevent plugging of the tubes with
slurry “chunks”, and not less than 12 BWG thickness because of the possibility of
erosion that may occur at the higher velocities.
4. Tool assumptions
1. Flow is steady and isothermal, and fluid properties are independent of time.
2. Fluid density is dependent on the local temperature only or is treated as constant
3. The pressure at a point in the fluid is independent of direction.
4. There are no energy sinks or sources along the streamline; flow stream mechanical
energy dissipation is idealized as zero.
5. The friction factor is considered constant with passage flow length.
6. Tube wall temperature viscosity effect is negligible.
7. Single phases for both shell and tube (two phases to be added soon)
1. Leff = L - 2*Ls , where Leff is effective tube length in mm for fixed tube sheet
2. Ls = 0.1* Ds
𝑁𝑡𝑐𝑐
3. Nss = 𝑖𝑛𝑡𝑒𝑔𝑒𝑟
6
4. DSB = 3.1+0.004* Ds
5. LTB = 0.4 mm
6. DotL = Ds - (12.5+( Ds /200))
References
1. Process Heat transfer, Kern.
2. JGC standard practice “Shell & tube Exchangers” JGS 210-120-1-24E
3. TOTAL process engineering design manual
4. Petroleum Refining Design and applications, vol. 4, A. Kayode Coker
5. GPSA Engineering Data Book Section 9 (Heat exchangers)
6. Rapid Design Algorithms for Shell-and-tube and Compact Heat Exchangers, Graham
Polley.
7. Effectively design shell-and-tube heat exchanges, Rajiv Mukherjee.
8. Energy and Process Optimization for the Process Industries, Frank Zhu.
9. Energy Management and Efficiency for the Process Industries, Alan P. Rossiter, Beth P.
Jones.
10. Applied Heat Transfer and Heat Exchangers, Aramco Engineering Encyclopedia
11. Process Heat Transfer Principles and Applications, Robert W. serth.
Other Recommended Readings and resources
1. Chemical process engineering by A. Kayode Coker vol 2 (present the manual calculations
of heat exchangers in different cases and compares results with commercial software.
Overall, an excellent book, worth checking out and should be added to your library)
2. Monitoring the Thermal Efficiency of Fouled Heat Exchangers: A Simplified Method, M.
A. S. Jerónimo et al
3. Data Reconciliation and Fouling Analysis in Heat Exchanger Network, Ahmad Nuruddin
bin Abdul Aziz.
4. Optimization of operating conditions for mitigating fouling in heat exchanger networks,
C. Rodriguez and R. Smith.
Resources
1. Bell-Delaware Calculations excel sheet
2. Kern Calculations excel sheet
3. NTU calculations excel sheet
Preliminary Design Example
Kern Method [4] page: 291
Use Kern’s Method to obtain a preliminary design for the following conditions
Kerosene API 42, 25,000 kg/hr (Shell side)
Crude oil API 34, 85,000 kg/hr (Tube side)
Kerosene Inlet Average Outlet
o
Temperature C 200 145 90
o
Specific heat kJ/kg. C 2.72 2.47 2.26
Thermal conductivity W/m. oC 0.13 0.132 0.135
Density Kg/m3 690 730 770
-2
Viscosity mNsm 0.22 0.43 0.8
Crude Oil
o
Temperature C 79 59.5 40
Specific heat kJ/kg. oC 2.09 2.05 2.01
o
Thermal conductivity W/m. C 0.133 0.134 0.135
Density Kg/m3 800 820 840
-2
Viscosity mNsm 2.4 3.2 4.3
Assumptions:
Uservice = 350 W/m2. oC
Length = 5000 mm, Do = 19.05 mm, Di = 14.83 mm, triangular pitch 30o, pitch ratio = 1.25
Allowable ∆Ps = 0.9 bar, ∆Pt = 0.9 bar,
Output
Table: input for tool was only S.G and temperatures.
P= = 0.2437, R = = 2.82
F = 0.8673 , ∆TCMTD= F * ∆TLMTD = 69.43
Q=UA∆TCMTD , A = Q/U∆TCMTD = 1886.8 x 103/(350x69.43) = 77.64 m2
Nt = A/πDoL = 77.64/(3.14x0.01905x5) = 259.46 (Say 260)
Tube cross sectional area = Ac = πDi2/4 = π(0.01483)2/4 = 0.0001727 m2
Area per pass = Nt/pnAc =(260/2x0.0001727)= 0.02246 m2
Gs = Wh/ρ = 85,000/3600*820= 0.02879 m3/s
Tube velocity = ut = 0.02879/0.02246 = 1.28 m/s
Estimate Shell diameter
From table 8.46 for two passes K1 = 0.249 , n1 = 2.207
.
𝐷 =𝐷 ∗ = 19.05 ∗ = 444 mm
.
Clearance is between 50 to 80 mm for split-ring floating head
Choose C= 56 mm, Ds = 444 +56 = 500 mm
∗ . ∗ .
𝑅 = = = 4864
. ∗
. ∗ ∗ . ∗
Pr= Cpμ/k = =48.96
.
Re= 4864, L/di =5000/14.84 = 337, from figure 8.55 the friction factor Jh
Jh = 3.5 x 10
Nu = Jh Re Pr 0.33(μ/μw)0.14
.
= 𝐽 𝑅𝑒𝑃𝑟 .
, assume negligible μ/μw
ℎ = 3.5 x 10 (0.134/0.01483)(4864)(48.96)0.33=555.48 W/m2. oC
In order to increase hi , tube passes is doubled, this will halve cross sectional area and double
the tube velocity.
ut = 2 x 1.28 = 2.56 m/s
Re = 2 x 4864 = 9728, Jh = 4.9 x 10
hi= 4.9 x 10 (0.134/0.01483)(9728)(48.96)0.33= 1555 W/m2. oC √√
From table 8.46 for two passes K1 = 0.175 , n1 = 2.285
.
𝐷 =𝐷 ∗ = 19.05 ∗ = 465.6 mm
.
Ds = 465.6 +56 = 521.6 mm
Baffle spacing trials
B = Ds/5 = 521.6/5 = 104.32
( ) . .
Shell flow area = = 𝑥522𝑥144 = 10853.04 mm2 =0.01853 m2
.
.
Equivalent shell diameter = De = (𝑝 − 0.913𝑑 ) =13.53 mm
Gs = Wh/ρ = 25,000/3600*730= 0.00951 m3/s
Shell velocity = us = 0.00951/0.01853 = 0.51 m/s
∗ . ∗ .
𝑅 = = = 11,714
. ∗
. ∗ ∗ . ∗
Pr= Cpμ/k = = 8.17
.
-3
jh= 5.8 x 10
.
Nu = = 𝐽 𝑅𝑒𝑃𝑟 .
𝑑
1 1 𝑑 1 𝑑 𝑑 ln 1
𝑑
= + + + +𝑅
𝑈 ℎ 𝑑 ℎ 𝑑 2𝑘 ℎ
19.05
1 1 19.05 1 19.05 0.001905 ln 1
= + 0.00035 + + 14.83 + + 0.0002
𝑈 1555 14.83 ℎ 14.83 2 ∗ 55 1325.69
= 440 W/m2. oC
Over Design = ∗ 100% = 25.7%
Pressure drops
Tube pressure drop
𝐿 μ ρv
∆Pt = 𝑁 8𝐽 + 2.5
𝐷 μ 2
=4[8x4.8x10 x337.15+2.5]820x2.56 /2 =96443.9 N/m2 = 0.96 bar
-3 2
a. Estimate the number of tubes that could fit in your shell diameter
b. Estimate shell Reynold’s number
θDS = 2cos-1(1 - 2Bc/100)
θCTL = 2cos-1(Ds(1 - 2*Bc/100)/DCTL)
FW = (θCTL - Sin(θCTL))/2π
SWG = (Ds²/8)(θDS - Sin(θDS))
SWT = Nt.FW(πDo²/4)
SW = SWG - SWT
Based on the assumption that baffle spacing is same as baffle cut area
Sm = SW
GW = mS/(Sm.SW)0.5
GS = mS/Sm
ReS = Do.GS / μS
c. Calculate shell heat transfer coefficient and pressure drop
Sm = B[(Ds - DOTL) + (DOTL - Do)(PT - Do)/PT,eff ] Di = Do – 2*thickness
PT,eff = PT for 30° and 90° layouts Ltb = 0.4
PT,eff = 0.707*PT for 45° layout
Ls= 0.1* Ds
GS = mS/Sm
ReS = Do.GS / μS Dsb = 3.1+0.004* Ds
PrS = CP,S.μS / kS Leff = L- 2* Ls
Dotl= Ds - (12.5+( Ds/200))
1.33
𝑗=𝑎 . (𝑅𝑒 )
𝑃
𝐷𝑜
𝑎
𝑎=
1 + 0.14(𝑅𝑒 )
µ .
ℎ = 𝑗. 𝐶 , . 𝐺 . (𝑃𝑟 ) . ( )
µ,
1.33
𝑓=𝑏 . (𝑅𝑒 )
𝑃
𝐷𝑜
𝑏
𝑏=
1 + 0.14(𝑅𝑒 )
ΔPIdeal = 2f(GS²/ρS)(μS/μS,W)0.14 Ntcc
For turbulent flow, ReT > 10,000, following equation developed by Petukhov-Kirillov
can be used.
Nu = (f/2)ReT.PrT/(1.07+12.7(f/2)0.5(PrT2/3-1))
f = (1.58 ln(ReT) - 3.28)-2
For transient flow, Nusselt number can be interpolated from Nu Laminar & Nu Turbulent.
ht = Nu.(kt/Di)(μt/μt, w)0.14
∆𝑃 = 𝑘 ∗ ℎ .
→ 11.66 = kt*(ht )3.5→ ht = 908 W/m².°C
∆𝑃 = 𝑘 ∗ ℎ .
→ 13.7 = ks*(hs )4.412→ hs = 6335 W/m².°C
A = 64.74 m2
For laminar flow, ReT < 2300, Sieder and Tate correlation is used for Nusselt's nubmer.
Nu = 1.86(ReT.PrT.Di/Leff)1/3
Leff = L - 2*Ls , where Leff is effective tube length
For turbulent flow, ReT > 10,000, following equation developed by Petukhov-Kirillov
can be used.
Nu = (f/2)ReT.PrT/(1.07+12.7(f/2)0.5(PrT2/3-1))
f = (1.58 ln(ReT) - 3.28)-2
For transient flow, Nusselt number can be interpolated from Nu Laminar & Nu Turbulent.
ht = Nu.(kt/Di)(μt/μt, w)0.14
if ht = ht.solved stop iteration, else change Re accordingly and start over.
∆𝑃 = 0.196 Kg/cm2
NTU Method
For the previous example:
Fluid Data Hot Cold
Flowrate kg/s 22.4 78.0
Inlet Temperature °C 100.0 7.1
Specific Heat kJ/kg.°K 2.407 4.19
Exchanger Data
Heat Exchanger Type Shell & Tube - 1 Shell - 2,4, .. Tube Passes
Heat Transfer Coefficient kW/m².°K 1.09
Area m² 52.30
Solution
UA = 52.3*1.088 = 56.9 kW/°C
Cmin = 22.4*3600*2.407* 0.239006*1.163 = 53952.94 W/oC
Cmax = 78*3600*4.19*0.239006 *1.163 = 327039.1 W/ oC
Cr = Cmin/Cmax = 0.1652
NTU = UA/Cmin = 56.9/53.95 = 1.06
Qmax = Cmin(T1-t1) = (22.4*3600*2.407* 0.239006) * (100-7.1) = 4309741 kCal
2
𝜀= = 0.6177
1 + exp (−𝑁𝑇𝑈 1 + 𝐶 )
1+𝐶 + 1+𝐶
1 − exp (−𝑁𝑇𝑈 1 + 𝐶 )
𝜔 = 𝑦𝜔 GPSA[2]
Using Charts
𝑍 = 0.291 − 0.08𝜔
1. Use fig 23-31 at operating
𝑀 = 𝑦𝑀 temperature to find KA
𝑇 =𝑇 2. Calculate Tr and Pr
3. Use fig. 23-32 to find thermal
𝑇 = (𝑇 𝑇 )
conductivity ratio
𝑉 =𝑉
4. K = ratio * KA
(𝑉+𝑉 ) Stiel and thodos Method is recommended
𝑉 =
8 For mixtures of simple gases, it doesn’t
∑ ∑ 𝑦𝑦𝑉 𝑇 apply to mixtures containing high CO2[2]
𝑇 =
𝑉 Miscellaneous gases can be obtained from
𝑉 = 𝑦𝑦𝑉 fig 23-33 and hydrocarbon gases from fig
23-34
𝑃 =𝑍 𝑅𝑇 /𝑉
/ Liquid Thermal conductivity
𝑇 𝑀 Ludwig’s [3] [Notes]
𝛤 = 210
𝑃 KL = A+ BT + CT2
𝑃 = 𝑉 /𝑉 KL: µcal/s.cm.oC
For organic compounds, the correlation for thermal
if ρrm<0.5: conductivity of
liquid as a function of temperature is [2]:
(λ−λ∘)ΓZcm5=1.22×10−2[exp(0.535ρr)−1] log10 kliq
if 0.5< ρrm <2.0: Thermal conductivity for liquid mixtures
(λ−λ∘)ΓZcm5=1.14×10−2[exp(0.67ρrm)−1.069] ∑ . .
if 2.0< ρrm <2.8: = ∑
or = ∑𝑤 , where
.
(λ−λ∘)ΓZcm5=2.6×10−3[exp(1.155ρrm)−2.016] Xi: Mole fraction of component i
Wi: weight fraction of component i
Tcm, Pcm, Zcm, Vcm are critical temperature, Mi: molecular weight of component i
pressure, compressibility factor and specific From graph 23-35 liquid paraffin
volume of mixture in oK, bar, -, cm3/mol hydrocarbons
Vm, Mm, λ, λ∘ are specific volume at 1 atm, The Sato-Riedel and Missenard methods
molecular weight of mixture g/mol, are-recommended.[6] [Notes]
thermal conductivity and thermal Missenard method
conductivity at 1 atm in W/m.k k/k∗=1+QPr0.7
yi, yj are mol. Fractions of mixture k*: k at low pressure (W/m.K)
components
𝛤 is the reduced, inverse thermal
conductivity in [W/m.K]-1
Density
Volume
Component Mole fraction Mol wt Weight, lb (60F),
, cu ft
lb/cf
Methane 0.20896 16.043 3.352 - -
Carbon dioxide 0.39730 44.010 17.485 51.016 0.3427
Ethane 0.01886 30.070 0.567 - -
Propane 0.02387 44.097 1.053 31.619 0.0333
n-Butane 0.03586 58.123 2.084 36.423 0.0572
n-Pentane 0.02447 72.150 1.766 39.36 0.0449
n-Hexane 0.01844 86.177 1.589 41.4 0.0384
n-Heptane 0.02983 100.204 2.989 42.92 0.0696
n-Octane 0.02995 114.231 3.421 44.09 0.0776
n-Decane 0.18208 142.285 25.907 45.79 0.5658
n-Tetradecane 0.03038 198.394 6.027 47.85 0.1260
Total 1.00000 66.241
Given:
Pressure = 300 psia
Temperature = 50 deg F
Methane = 80 mol %
Nitrogen = 15 mol %
Carbon dioxide = 5 mol %
Solution
Steps
Tcm = 335.9 deg R (from Table above)
3
Vcm = 1.558 ft /lb (from Table above)
Pcm = 0.2877*10.73*335.9/1.558=665.6 psia
Zcm = 0.2877 dimensionless
5.4402*(335.9)0.1666/(665.60.6666)/(19.23)0.5
ξ = dimensionless (from Eq. 23-20)
=0.043
Tr = 1.518 dimensionless
ξμA = 0.000489 dimensionless (from Eq. 23-21)
μA = 0.000489/0.043=0.014 cP (from Eq. 23-22)
Thermal Conductivity
Example 23-12 Find the thermal conductivity of a natural gas.
Given:
Pressure = 700 psia
Temperature = 300 deg F
Tc = 440 deg R
Pc = 660 psia
Solution Steps
kA = 0.0248 Btu/[(hr-sq.ft.-deg (from Fig.23-35)
F)/ft]
Tr = (300+459.67)/440=1.73 dimensionless
Pr = 700/660=1.06 dimensionless
k/kA = 1.15 dimensionless (from Fig.23-36)
k = 1.15*0.0248=0.0285 Btu/[(hr-sq.ft.-deg
F)/ft]
Example 23-13 -- Find the thermal conductivity of the gaseous mixture shown in Fig. 23-41 at
200°F and one atmosphere.
Given:
Component Mole Thermal conductivity, MW MW-3 (yi)MW-3 (yi)(ki)MW-3
Fraction Btu/(hr-ft-degF)
CO2 0.10 0.0127 44.01 3.530 0.3530 0.00448
H2S 0.20 0.0136 34.076 3.242 0.6483 0.00882
N2 0.05 0.0175 28.013 3.037 0.1518 0.00266
CH4 0.60 0.0258 16.043 2.522 1.5131 0.03904
C2H6 0.05 0.0176 30.07 3.109 0.1555 0.00274
Total 1.00 2.8218 0.05773
Solution
km = 0.05773/2.8218=0.0205 Btu/[(hr-sq.ft.-degF)/ft]
Calculate Naphtha properties (S.G =0.727) at 25oC and vapor naphtha at 2 kg/cm2.a & 150o C
density 0.727*1000*(1-0.0005*((25*1.8+32)-60)) 720.8 Kg/m3
Cp = (1/ 𝑠𝑔))*(0.388+0.00045*(t*1.8+32)) 0.49569 Kcal/kg.C
Cv Cp-(0.09/sg) 0.3719 Kcal/kg.C
thermal conductivity = (0.813/sg)*(1-0.0003*((t*1.8+32)-32)) *0.1442279 0.15911 W/m.k
viscosity 3.518-0.01591*t+1.734*10-5t2 2.39574 cP
Mol. Weight 1/(0.0001644 ∗ 𝐴𝑃𝐼 − 0.000972) 106.29846 -
16.018463 ∗ 𝑚_𝑤𝑡 ∗ 𝑝 ∗ 14.7 5.92723 Kg/m3
Gas Density =
10.732 ∗ 𝑧 ∗ (𝑡 + 459.67)
Gas thermal ( -0.000500777+1.0906*10-5*(t*1.8+32)+0.061137256/√M.wt + 0.02316 W/m.k
conductivity 0.000158966*(t*1.8+32)/√M.wt)*1.730735
=-0.0092696 +T0.5(0.001383-5.9712x10-5M0.5) 0.00771 cP
Gas Viscosity +1.1249x10-5M
Source 1: Bureau of standards report
Report Conclusions
Specific volume of vapor: The data given in Table 9 for products of gravity 50°
to 150° A. P. I. appear to be sufficiently reliable for most industrial purposes. The
data given for products of gravity 20° to 50° A. P. I. are admittedly approximate and
are included mainly for the benefit of those users of petroleum products who have no
information at hand as to the source or the volatility of these products. Additional
experimental data on molecular weights of products of gravity 10° to 40° A. P. I.
appear desirable.
Thermal conductivity: The data given in Table 10 appear to be sufficient for all
practical uses of such data.
Specific heat: The data given in Tables 12 and 13 appear to be sufficiently
reliable for all practical applications wherein only moderate pressures (less than 50
lbs. /in.2) are involved. Additional experimental data on petroleum vapors and also
on petroleum liquids at high temperatures and pressures are desirable.
Heat content: The data given in Tables 16 to 19 should provide simplicity and
convenience in calculations of the quantities of heat involved in the heating and
cooling of petroleum products.
Figure 23-15: Density
correction for
compressibility of
hydrocarbon liquids
Figure 23-14: Pseudo Density of mixtures
Figure 23-17: Density Correction for Thermal Expansion of Hydrocarbon Liquids
FIG. 23-33
Thermal Conductivity
of Miscellaneous
Gases
at One Atmosphere
FIG. 23-34
Thermal Conductivity
of H.Cs Gases
at One Atmosphere
Figure 23-4: Standing-Katz Chart
Figure 23-4: Standing-Katz Chart
Supplementary Material
Condensing Heat
Transfer
Steam, ammonia No noncondensable 8,000-12,000
Pure component, 0.1 bar abs, no
Light organics 2,000-5,000
noncondensable
Light
multicomponent
Medium condensing range, 1 bar abs 1,000-2,500
mixture, all
condensable
Medium
multicomponent
Medium condensing range, 1 bar abs 600-1,500
mixture, all
condensable
Heavy
multicomponent
Medium condensing range, 1 bar abs 300-600
mixture, all
condensable
Vaporizing Heat
Transfer
Fixed Tubes One-Pass ¾ in. on 15/ 16 in ∆ 33 69 105 135 193 247 307 391 481 553 663 763 881 1019 1143 1,269
Fixed Tubes One-Pass ¾ in. on 1in. ∆ 33 57 91 117 157 217 277 343 423 493 577 667 765 889 1007 1,127
Fixed Tubes One-Pass ¾ in. on 1in ◽ 33 53 85 101 139 183 235 287 355 419 495 587 665 765 865 965
Fixed Tubes One-Pass 1 in. on ¼ in. ∆ 15 33 57 73 103 133 163 205 247 307 361 427 481 551 633 699
Fixed Tubes One-Pass 1 in. on 1 ¼ in ◽ 17 33 45 65 83 111 139 179 215 255 303 359 413 477 545 595
Fixed Tubes Two-Pass ¾ in. on 15/ 16 in ∆ 32 58 94 124 166 228 300 370 452 528 626 734 846 964 1088 1,242
Fixed Tubes Two-Pass ¾ in. on 1 in. ∆ 28 56 90 110 154 208 264 326 398 468 556 646 746 858 972 1,088
Fixed Tubes Two-Pass ¾ in. on 1 in ◽ 26 48 78 94 126 172 222 280 346 408 486 560 644 746 840 946
Fixed Tubes Two-Pass 1 in. on 1 ¼ in ∆ 16 32 52 62 92 126 162 204 244 292 346 410 462 530 608 688
Fixed Tubes Two-Pass 1 in. on 1 ¼ in ◽ 12 26 40 56 76 106 136 172 218 248 298 348 402 460 522 584
U Tubes Two-Pass ¾ in. on 15/ 16 in ∆ 8 34 64 94 134 180 234 304 398 460 558 648 768 882 1008 1,126
U Tubes Two-Pass ¾ in. on 1in. ∆ 8 26 60 72 108 158 212 270 336 406 484 566 674 772 882 1,000
U Tubes Two-Pass ¾ in. on 1in ◽ 12 30 52 72 100 142 188 242 304 362 436 506 586 688 778 884
U Tubes Two-Pass 1 in. on 1 ¼ in ∆ XX 8 26 42 58 84 120 154 192 234 284 340 396 466 532 610
U Tubes Two-Pass 1 in. on 1 ¼ in ◽ XX 12 22 38 58 76 100 134 180 214 256 304 356 406 464 526
Fixed Tubes Four-Pass ¾ in. on 15/ 16 in ∆ XX 48 84 108 154 196 266 332 412 484 576 680 788 904 1024 1,072
Fixed Tubes Four-Pass ¾ in. on 1in. ∆ XX 44 72 96 134 180 232 292 360 424 508 596 692 802 912 1,024
Fixed Tubes Four-Pass ¾ in. on 1in ◽ XX 48 72 88 126 142 192 242 308 366 440 510 590 688 778 880
Tubesheet type and no. of
passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
Fixed Tubes Four-Pass 1 in. on 1 ¼ in ∆ XX 24 44 60 78 104 138 176 212 258 308 368 422 486 560 638
Fixed Tubes Four-Pass 1 in. on 1 ¼ in ◽ XX 24 40 48 74 84 110 142 188 214 260 310 360 414 476 534
U Tubes Four-Pass ¾ in. on 15/ 16 in ∆ XX 28 56 84 122 166 218 286 378 438 534 622 740 852 976 1,092
U Tubes Four-Pass ¾ in. on 1in. ∆ XX 20 52 64 98 146 198 254 318 386 462 542 648 744 852 968
U Tubes Four-Pass ¾ in. on 1in ◽ XX 24 44 64 90 130 174 226 286 342 414 482 560 660 748 852
U Tubes Four-Pass 1 in. on 1 ¼ in ∆ XX XX 20 36 50 74 110 142 178 218 266 322 376 444 508 584
U Tubes Four-Pass 1 in. on 1 ¼ in ◽ XX XX 16 32 50 66 90 122 166 198 238 286 336 384 440 500
Fixed Tubes Six-Pass ¾ in. on 15/ 16 in ∆ XX XX XX 80 116 174 230 294 372 440 532 632 732 844 964 1,106
Fixed Tubes Six-Pass ¾ in. on 1in. ∆ XX XX XX 66 104 156 202 258 322 388 464 548 640 744 852 964
Fixed Tubes Six-Pass ¾ in. on 1in ◽ XX XX XX 54 78 116 158 212 266 324 394 460 536 634 224 818
Fixed Tubes Six-Pass 1 in. on 1 ¼ in ∆ XX XX XX 34 56 82 112 150 182 226 274 338 382 442 514 586
Fixed Tubes Six-Pass 1 in. on 1 ¼ in ◽ XX XX XX XX 44 66 88 116 154 184 226 268 318 368 430 484
U Tubes Six-Pass ¾ in. on 15/ 16 in ∆ XX XX XX 74 110 156 206 272 358 416 510 596 716 826 944 1,058
U Tubes Six-Pass ¾ in. on 1in. ∆ XX XX XX 56 88 134 184 238 300 366 440 518 626 720 826 940
U Tubes Six-Pass ¾ in. on 1in ◽ XX XX XX 56 80 118 160 210 268 322 392 458 534 632 718 820
U Tubes Six-Pass 1 in. on 1 ¼ in ∆ XX XX XX 30 42 68 100 130 168 206 252 304 356 426 488 562
U Tubes Six-Pass 1 in. on 1 ¼ in ◽ XX XX XX XX 42 60 80 110 152 182 224 268 316 362 420 478
Tubesheet type and no. of
passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
Fixed Tubes eight-Pass ¾ in. on 15/ 16 in ∆ XX XX XX XX 94 140 198 258 332 398 484 576 682 790 902 1,040
Fixed Tubes eight-Pass ¾ in. on 1in. ∆ XX XX XX XX 82 124 170 224 286 344 422 496 588 694 798 902
Fixed Tubes eight-Pass ¾ in. on 1in ◽ XX XX XX XX XX 94 132 174 228 286 352 414 490 576 662 760
Fixed Tubes eight-Pass 1 in. on 1 ¼ in ∆ XX XX XX XX XX 66 90 120 154 190 240 298 342 400 466 542
Fixed Tubes eight-Pass 1 in. on 1 ¼ in ◽ XX XX XX XX XX XX 74 94 128 150 192 230 280 334 388 438
U Tubes eight-Pass ¾ in. on 15/ 16 in ∆ XX XX XX 68 102 142 190 254 342 398 490 578 688 796 916 1,032
U Tubes eight-Pass ¾ in. on 1in. ∆ XX XX XX 52 82 122 170 226 286 350 422 498 600 692 796 908
U Tubes eight-Pass ¾ in. on 1in ◽ XX XX XX 48 70 106 146 194 254 306 374 438 512 608 692 792
U Tubes eight-Pass 1 in. on 1 ¼ in ∆ XX XX XX 24 38 58 90 118 154 190 238 290 340 404 464 540
U Tubes eight-Pass 1 in. on 1 ¼ in ◽ XX XX XX XX 34 50 70 98 142 170 206 254 300 344 396 456
TEMA P or S One-Pass 1.25 in. on 1 9/16 in ◽ 12 21 29 38 52 70 85 108 136 154 184 217 252 289 329 372
TEMA P or S Two-Pass 1.25 in. on 1 9/16 in ◽ 12 12 28 34 48 66 84 108 128 154 180 212 248 276 316 368
TEMA P or S Four-Pass 1.25 in. on 1 9/16 in ◽ 4 8 16 34 44 56 70 100 128 142 158 204 234 270 310 354
Tubesheet type and no.
of passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
TEMA P or S Six-Pass 1.25 in. on 1 9/16 in ◽ 0 12 18 24 48 50 80 96 114 136 172 198 236 264 304 340
U Tubes Two-Pass 1.25 in. on 1 9/16 in ◽ 4 12 26 36 44 60 82 100 128 154 176 212 242 280 324 358
U Tubes Four-Pass 1.25 in. on 1 9/16 in ◽ 4 8 20 28 44 60 76 100 120 148 172 204 240 280 312 352
U Tubes Six-Pass 1.25 in. on 1 9/16 in ◽ 6 12 12 15 32 56 79 100 120 130 160 198 234 274 308 350
TEMA L or M One-
Pass 1.25 in. on 1 9/16 in ∆ 15 27 38 55 66 88 117 136 170 198 237 268 312 357 417 446
TEMA L or M Two-
Pass 1.25 in. on 1 9/16 in ∆ 10 22 36 44 64 82 106 134 164 188 228 266 304 346 396 446
TEMA L or M Four-
Pass 1.25 in. on 1 9/16 in ∆ 8 16 26 42 52 78 98 124 146 166 208 242 284 322 372 422
TEMA L or M Six-Pass 1.25 in. on 1 9/16 in ∆ 12 12 24 40 50 68 96 108 148 168 192 236 276 324 364 408
TEMA P or S One-
Pass 1.25 in. on 1 9/16 in ∆ 13 18 33 38 57 81 100 126 159 183 208 249 291 333 372 425
TEMA P or S Two-
Pass 1.25 in. on 1 9/16 in ∆ 10 20 26 44 58 72 94 120 146 172 206 238 282 326 368 412
TEMA P or S Four-
Pass 1.25 in. on 1 9/16 in ∆ 4 8 26 34 48 62 86 116 132 150 190 224 262 298 344 394
TEMA P or S Six-Pass 1.25 in. on 1 9/16 in ∆ 0 12 18 24 44 68 80 102 132 148 180 220 256 296 336 384
U Tubes Two-Pass 1.25 in. on 1 9/16 in ∆ 6 14 28 34 52 72 90 118 148 172 200 242 282 326 362 416
Tubesheet type and no.
of passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
U Tubes Four-Pass 1.25 in. on 1 9/16 in ∆ 4 12 20 28 48 68 84 112 132 160 188 228 264 308 344 396
U Tubes Six-Pass 1.25 in. on 1 9/16 in ∆ 6 12 18 30 40 64 78 102 120 152 180 216 250 292 336 384
TEMA P or S One-
Pass 5/8 in. on 13/16 in ◽ 55 88 140 178 245 320 405 502 610 700 843 970 1127 1288 1479 1647
TEMA P or S Two-
Pass 5/8 in. on 13/16 in ◽ 48 78 138 172 232 308 392 484 584 676 812 942 1096 1250 1438 1604
TEMA P or S Four-
Pass 5/8 in. on 13/16 in ◽ 34 62 112 146 208 274 352 442 536 618 742 868 1014 1172 1330 1520
TEMA P or S Six-Pass 5/8 in. on 13/16 in ◽ 24 56 100 136 192 260 336 424 508 600 716 840 984 1148 1308 1480
U Tubes Two-Pass 5/8 in. on 13/16 in ◽ 52 90 140 180 246 330 420 510 626 728 856 998 1148 1318 1492 1684
U Tubes Four-Pass 5/8 in. on 13/16 in ◽ 40 80 128 164 232 312 388 488 596 692 816 956 1108 1268 1436 1620
U Tubes Six-Pass 5/8 in. on 13/16 in ◽ 32 74 108 148 216 292 368 460 562 644 780 920 1060 1222 1388 1568
TEMA L or M One-
Pass 3/4 in. on 15/16 in ∆ 64 85 122 151 204 264 332 417 495 579 676 785 909 1035 1164 1304
TEMA L or M Two-
Pass 3/4 in. on 15/16 in ∆ 48 72 114 142 192 254 326 396 478 554 648 762 878 1002 1132 1270
TEMA L or M Four-
Pass 3/4 in. on 15/16 in ∆ 34 52 94 124 166 228 290 364 430 512 602 704 814 944 1062 1200
TEMA L or M Six-Pass 3/4 in. on 15/16 in ∆ 24 50 96 112 168 220 280 348 420 488 584 688 792 920 1036 1168
TEMA P or S One-
Pass 3/4 in. on 15/16 in ∆ 34 60 109 126 183 237 297 372 450 518 618 729 843 962 1090 1233
Tubesheet type and no. of
passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
TEMA P or S Two-Pass 3/4 in. on 15/16 in ∆ 32 62 98 120 168 228 286 356 430 498 602 708 812 934 1064 1196
TEMA P or S Four-Pass 3/4 in. on 15/16 in ∆ 16 52 78 106 146 202 258 324 392 456 548 650 744 868 990 1132
TEMA P or S Six-Pass 3/4 in. on 15/16 in ∆ 18 44 68 100 136 192 248 316 376 444 532 624 732 840 972 1100
U Tubes Two-Pass 3/4 in. on 15/16 in ∆ 32 64 98 126 180 238 298 370 456 534 628 736 846 978 1100 1238
U Tubes Four-Pass 3/4 in. on 15/16 in ∆ 24 52 88 116 160 224 280 352 428 500 600 696 812 928 1060 1200
U Tubes Six-Pass 3/4 in. on 15/16 in ∆ 24 52 78 108 148 204 262 334 408 474 570 668 780 904 1008 1152
TEMA P or S One-Pass 3/4 in. on 1 in ◽ 28 52 80 104 136 181 222 289 345 398 477 554 637 730 828 937
TEMA P or S Two-Pass 3/4 in. on 1 in ◽ 26 48 76 90 128 174 220 272 332 386 456 532 624 712 812 918
TEMA P or S Four-Pass 3/4 in. on 1 in ◽ 16 44 66 70 128 154 204 262 310 366 432 510 588 682 780 882
TEMA P or S Six-Pass 3/4 in. on 1 in ◽ 12 24 56 80 114 160 198 260 308 344 424 496 576 668 760 872
U Tubes Two-Pass 3/4 in. on 1 in ◽ 28 52 78 96 136 176 224 284 348 408 480 562 648 748 848 952
U Tubes Four-Pass 3/4 in. on 1 in ◽ 24 44 72 92 132 176 224 280 336 392 468 548 636 728 820 932
U Tubes Six-Pass 3/4 in. on 1 in ◽ 12 32 70 90 120 160 224 274 328 378 460 530 620 718 816 918
TEMA L or M One-Pass 3/4 in. on 1 in ∆ 42 73 109 136 183 237 295 361 438 507 592 692 796 909 1023 1155
Tubesheet type and no.
of passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
TEMA L or M Two-
Pass 3/4 in. on 1 in ∆ 40 66 102 128 172 228 282 346 416 486 574 668 774 886 1002 1124
TEMA L or M Four-
Pass 3/4 in. on 1 in ∆ 26 52 88 112 146 208 258 318 382 448 536 632 732 836 942 1058
TEMA L or M Six-Pass 3/4 in. on 1 in ∆ 24 44 80 102 148 192 248 320 372 440 516 604 708 812 920 1032
TEMA P or S One-
Pass 3/4 in. on 1 in ∆ 31 56 88 121 159 208 258 320 400 450 543 645 741 843 950 1070
TEMA P or S Two-
Pass 3/4 in. on 1 in ∆ 26 48 78 106 148 198 250 314 384 442 530 618 716 826 930 1052
TEMA P or S Four-
Pass 3/4 in. on 1 in ∆ 16 42 62 94 132 182 228 290 352 400 488 574 666 760 878 992
TEMA P or S Six-Pass 3/4 in. on 1 in ∆ 12 40 68 88 132 180 220 276 336 392 468 556 648 740 856 968
U Tubes Two-Pass 3/4 in. on 1 in ∆ 32 52 84 110 152 206 266 330 400 472 554 648 744 852 974 1092
U Tubes Four-Pass 3/4 in. on 1 in ∆ 24 48 76 100 140 188 248 316 384 440 528 616 716 816 932 1056
U Tubes Six-Pass 3/4 in. on 1 in ∆ 24 40 74 98 136 182 234 296 356 424 502 588 688 788 908 1008
TEMA P or S One-
Pass 1 in. on 1.25 in ◽ 17 30 52 61 85 108 144 173 217 252 296 345 402 461 520 588
TEMA P or S Two-
Pass 1 in. on 1.25 in ◽ 12 30 48 56 78 108 136 166 208 240 280 336 390 452 514 572
Tubesheet type and no. of
passes tube size and pitch 8 10 12 13¼ 15¼ 17¼ 19¼ 21¼ 23¼ 25 27 29 31 33 35 37
TEMA P or S Four-Pass 1 in. on 1.25 in ◽ 8 16 42 52 62 104 130 154 194 230 270 310 366 432 494 562
TEMA P or S Six-Pass 1 in. on 1.25 in ◽ 12 18 24 50 64 96 114 156 192 212 260 314 368 420 484 548
U Tubes Two-Pass 1 in. on 1.25 in ◽ 14 30 44 60 80 104 132 172 212 244 290 340 400 456 518 584
U Tubes Four-Pass 1 in. on 1.25 in ◽ 8 24 40 48 72 100 132 168 204 240 284 336 384 444 504 576
U Tubes Six-Pass 1 in. on 1.25 in ◽ 6 12 32 44 74 100 120 148 198 230 274 328 372 440 502 566
TEMA L or M One-Pass 1 in. on 1.25 in ∆ 27 42 64 81 106 147 183 226 268 316 375 430 495 579 645 729
TEMA L or M Two-Pass 1 in. on 1.25 in ∆ 26 40 66 74 106 134 176 220 262 302 360 416 482 554 622 712
TEMA L or M Four-Pass 1 in. on 1.25 in ∆ 8 34 52 62 88 124 150 204 236 274 336 390 452 520 586 662
TEMA L or M Six-Pass 1 in. on 1.25 in ∆ 12 24 44 56 92 114 152 186 228 272 324 380 448 504 576 648
TEMA P or S One-Pass 1 in. on 1.25 in ∆ 18 33 51 73 93 126 159 202 249 291 345 400 459 526 596 672
TEMA P or S Two-Pass 1 in. on 1.25 in ∆ 14 28 48 68 90 122 152 192 238 278 330 388 450 514 584 668
TEMA P or S Four-Pass 1 in. on 1.25 in ∆ 8 16 42 52 78 112 132 182 216 250 298 356 414 484 548 626
TEMA P or S Six-Pass 1 in. on 1.25 in ∆ 12 18 44 44 76 102 136 172 212 240 288 348 400 464 536 608
U Tubes Two-Pass 1 in. on 1.25 in ∆ 14 28 52 64 90 122 152 196 242 286 340 400 456 526 596 668
U Tubes Four-Pass 1 in. on 1.25 in ∆ 12 24 40 56 80 112 140 180 224 264 320 380 436 504 572 636
U Tubes Six-Pass 1 in. on 1.25 in ∆ 6 24 40 52 78 102 136 176 216 246 300 352 414 486 548 614
Source : Source: Process Heat Transfer Principles and Applica ons & Petroleum Refining Design and applica ons, vol. 4, A. Kayode coker
Birmingham Wire Gage (inches) (mm)
BWG
7 0.180 4.572
8 0.165 4.191
9 0.148 3.759
10 0.134 3.404
11 0.120 3.048
12 0.109 2.769
13 0.095 2.413
14 0.083 2.108
15 0.072 1.829
16 0.065 1.651
17 0.058 1.473
18 0.049 1.245
19 0.042 1.067
20 0.035 0.889
21 0.032 0.813
22 0.028 0.711
23 0.025 0.635
24 0.022 0.559
25 0.020 0.508
26 0.018
27 0.016
28 0.014
29 0.013
30 0.012
31 0.010
32 0.009
33 0.008
34 0.007
35 0.005
36 0.004
Source: Heat Exchangers: Selec on, Ra ng, and Thermal Design
Source: Heat Exchangers: Selec on, Ra ng, and Thermal Design
Fired Heater Efficiency Calculations
This tool was developed for process engineers to quickly estimate heaters efficiency using heat loss
simplified method and the infamous API 560 also indirect heat loss method and access a list of useful
tables on site. As a part of a larger project to develop what is similar to Carl Branan’s book “process
engineers Pocket Handbook” these tools would allow a process engineer to quickly calculate/estimate
equipment efficiencies or sizing using standardized calculations
The aim here is to take little-known data from the field (flow, pressures, temperatures, compositions..etc.)
and use it as input for a rough estimation without having to return to the office to use commercial
software or calculations Excel sheets to validate or to calculate. Additionally, these tools may also serve as
a gathered data validation tool.
Nomenclature
1. Simplified Equation[2]
1. Collect data
Fluids properties
• Fuel composition (S.G., CH ratio & LHV if fuel, Complete composition & LHV if gas)
• Air & Fuel inlet conditions (humidity% and T)
• Calorific values (preferred for fuel oil but can be calculated for gas)
𝑚 = ∑𝑛𝑖=0 𝑥𝑖 . 𝑚𝑖 where:
m: amount formed or required of a flue gas (Air, N2, H2O, CO2 or SO2) per lb fuel.
xi : mass fraction of fuel gas component.
mi: amount formed by the fuel gas component per lb fuel.
n: number of fuel gas components.
3. Calculate pounds of moisture per pound of fuel = Pounds of wet air per pound of fuel
- Air required:
a. Moisture in Air
𝐴𝑖𝑟 𝑟𝑒𝑞𝑢𝑖𝑟𝑒𝑑
Pounds of wet air per pound of fuel required=
1−𝑀𝑜𝑖𝑠𝑡𝑢𝑟𝑒 𝑖𝑛 𝑎𝑖𝑟
= Pounds of wet air per pound of fuel - Air required per pound of fuel
5. Calculate total pounds of H2O per pound of fuel (corrected for excess air)
𝑝𝑒𝑟𝑐𝑒𝑛𝑡 𝑒𝑥𝑐𝑒𝑠𝑠 𝑎𝑖𝑟
Total lbs of H2O /lb of fuel (corrected for excess air) = 100
× 𝑝𝑜𝑢𝑛𝑑𝑠 𝑜𝑓 𝑚𝑜𝑖𝑠𝑡𝑢𝑟𝑒 𝑝𝑒𝑟 𝑝𝑜𝑢𝑛𝑑 𝑜𝑓 𝑓𝑢𝑒𝑙
7. Calculate radiant section Heat Loss (Assume 2% heat loss of total LHV)
hr = 2*0.01*hl
8. Calculate hf Fuel sensible Heat correction BTU/lb of fuel
hf = Cpf.(Tf – Td)
where Cpf is fuel heat capacity, Tf: fuel temperature and Td design datum temperature
9. Calculate ha air sensible heat correction BTU/lb of fuel
ha = Cpair (Tambient – Td).ma/mf
ma/mf is total air per lb fuel
10. Calculate atomizing medium heat correction hm BTU/lb of fuel
hm = Cpm(Tm – Td).mm/mf or mm/mf . ∆Hm
where Cpm is steam heat capacity, Tf: steam temperature and Td design datum temperature
11. Calculate Efficiencies
a. Net Thermal Efficiency
𝐻𝑒𝑎𝑡 𝑈𝑠𝑒𝑑
𝜂=
𝐻𝑒𝑎𝑡 𝑐𝑜𝑛𝑡𝑒𝑛𝑡 𝑓𝑜𝑟 𝑎𝑙𝑙 𝑖𝑛𝑝𝑢𝑡
((𝐻𝑒𝑎𝑡 𝐹𝑟𝑜𝑚 𝐹𝑢𝑒𝑙 + 𝐴𝑖𝑟 𝐻𝑒𝑎𝑡 + 𝐹𝑢𝑒𝑙 𝐻𝑒𝑎𝑡 + 𝐴𝑡𝑜𝑚𝑖𝑧𝑖𝑛𝑔 𝑠𝑡𝑒𝑎𝑚 ℎ𝑒𝑎𝑡)‐ (𝑟𝑎𝑑𝑖𝑎𝑛𝑡 𝐿𝑜𝑠𝑠 + 𝑆𝑡𝑎𝑐𝑘 𝐿𝑜𝑠𝑠))
=
(𝐻𝑒𝑎𝑡 𝐹𝑟𝑜𝑚 𝐹𝑢𝑒𝑙 + 𝐴𝑖𝑟 𝐻𝑒𝑎𝑡 + 𝐹𝑢𝑒𝑙 𝐻𝑒𝑎𝑡 + 𝐴𝑡𝑜𝑚𝑖𝑧𝑖𝑛𝑔 𝑠𝑡𝑒𝑎𝑚 ℎ𝑒𝑎𝑡)
(ℎ𝐿 + ℎ𝑎 + ℎ𝑓 + ℎ𝑚 ) − (ℎ𝑟 − ℎ𝑠 )
=
ℎ𝐿 + ℎ𝑎 + ℎ𝑓 + ℎ𝑚
b. Fuel Efficiency
2. Use Carbon ratio and other Oil content (H2O, ash, sulfur, Sodium and others) to estimate C% and
H% and other contents
References
1. API 560 ANNEX G.
2. Aramco stds: calculating process heater thermal efficiency
Recommended Readings
1. Aramco stds: calculating process heater thermal efficiency (extended document on the
various methods of estimating thermal efficiencies (direct & indirect methods)
2. Correctly Modeling and Calculating Combustion Efficiencies In Fired Equipment,
David Schmitt
3. Get the most from your fired heater, ashutosh Garg and H. Ghosh chemical engineering.
4. Optimize fired heater operations to save money , A. Garg, Hydrocarbon processing.
Column Column Column Column
Column 1 Column 2 3 4 Column 5 Column 6 Column 7 Column 8 Column 9 Column 10 Column 11 12 Column 13 14
Net SO2
Heating Dry Air H2O N2
Heating CO2 formed formed
Fuel Component total Value required Dry Air CO2 formed H2O formed
Volume value (pounds of N2 formed (pounds
M.Wt weight (british (pounds of required formed (pounds of formed (pounds
fraction (british CO2 per (lbs) of SO2
(lbs) thermal air per (lbs) (lbs) H2O per (lbs) of N2 per
thermal pound) per
units) pound) pound) pound)
units) pound)
Hydrogen 73.98 2.016 1.49 51600 76958.14 34.29 51.14 0 0 8.94 13.33344 26.36 39.314274 0
Methane 5.42 16 0.87 21500 18644.8 17.24 14.95 2.74 2.376128 2.25 1.9512 13.25 11.4904 0
Ethane 3.72 30.1 1.12 20420 22864.68 16.09 18.02 2.93 3.28078 1.8 2.015496 12.37 13.8509364 0
Ethylene 28.1 0.00 20290 0.00 14.79 0.00 3.14 0 1.28 0 11.36 0 0
Propane 4.36 44.1 1.92 19930 38320.61 15.68 30.15 2.99 5.749052 1.63 3.134099 12.05 23.169258 0
Propylene 42.1 0.00 19690 0.00 14.79 0.00 3.14 0 1.28 0 11.36 0 0
Butane 4.89 58.1 2.84 19670 55884.24 15.46 43.92 3.03 8.608503 1.55 4.40369 11.88 33.7521492 0
Butylene 56.1 0.00 19420 0.00 14.79 0.00 3.14 0 1.28 0 11.36 0 0
Pentane 3.93 72.1 2.83 19500 55253.84 15.33 43.44 3.05 8.642267 1.5 4.250295 11.78 33.3789834 0
Hexane 0.23 86.2 0.20 19390 3844.26 15.24 3.02 3.06 0.606676 1.46 0.28946 11.71 2.3216246 0
Hydrogen Sulfide 0.68 34.1 0.23 6550 1518.81 6.08 1.41 0 0 0.53 0.122896 4.68 1.0851984 1.88
Water 18 0 0 0 1 0 0
Total per pound of fuel 100.00 1.00 22242.01763 16.7696231 2.381641 2.40094 12.9521469
EXCESS AIR AND RELATIVE HUMIDITY WORKSHEET
Atomizing steam:
0 Pounds per pound of fuel (assumed or measured)
= 0.076328 50% 18
14.696 100 28.85
where:
Pvapor = vapor pressure of water at ambient temperature in psia (from steam tables)
16.76962
1 -0.000016
16.770
Pounds of moisture per pound of fuel = Pounds of wet air per pound of fuel - Air required
= 16.770 -16.7696
= 0.00027
Pounds of H2O per pound of fuel = H2O formed + Pounds of moisture per pound of fuel +
atomizing steam
= 2.400944 +0.00027 +0
= 2.401
= 3.761878
16.76962
= 22.43%
EXCESS AIR AND RELATIVE HUMIDITY WORKSHEET
total pounds of H2O per pound of fuel (corrected for excess air)
= 2.401
NOTE: All values used in the calculations above shall be on per-pound-of-fuel basis. Numbers in parentheses indicate values to be
taken from the "total per pound fuel" line of the combustion work sheet, and letters in parentheses indicate values to be taken
from the corresponding line of this work sheet.
FLUE GAS WORK SHEET
150 C
302 F
Enthalpy at
Component pounds of
Tc (Btu per Heat Content (BTU per pound of
Wet Mol. % Wet Wt% component T/year
pound fuel)
formed/hr
pounds of component formed per formed)
pound of fuel
Carbon
dioxide 2.38 6.93% 11.1% 17482.30 62936.29 52.81542 125.7874
water vapor 2.40 17.07% 11.2% 17626.44 63455.17 101.6308 244.0436
Nitrogen 15.92 72.77% 74.0% 116889.47 420802.09 57.49809 915.6013
Oxygen 0.79 3.16% 3.7% 53.63188 42.36888
Air 3.76 27613.86 99409.89
SO2 0.0355 0.07% 0.2% 39.15918 1.389333
total 21.53 100.00% 100.00% 179612.07 - 1327.801
Design Basis
Design datum temperature, Td = 60 oF
Ambient air temperature, Ta,a = 28 oF
Air temperature, Ta = 28 oF
Relative Humidity = 50% %
Fuel:
Av specific heat capacity of the air, over dataum cpa = 0.24 Btu/lb.oF
Air sensible massic heat correction, Δha = -158 Btu/lb
Av specific heat capacity of the fuel, over datum cpf = 0.48 Btu/lb.oF
Fuel sensible massic heat correction, Δhf = 19.2 Btu/lb
Steam temperature 0.0
Atomizing medium sensible massic heat correction, Δhm = 0.00 Btu/lb
Efficiency calculations
Net thermal efficiency, expressed as a %, e
% = 91.5%
Gross thermal efficiency, expressed as a %, eg
Link: https://processpocket.streamlit.app
This tool was developed for process engineers to quickly estimate Heaters’ maximum skin temperature
using the method described in API 530 and ISO 13704 on site. As a part of a larger project to develop
what is similar to Carl Branan’s book “process engineers Pocket Handbook” these tools would allow a
process engineer to quickly calculate/estimate equipment efficiencies or sizing using standardized
calculations.
The aim here is to take little-known data from the field (flow, pressures, temperatures,
compositions..etc.) and use it as input for a rough estimation without having to return to the office to
use commercial software or calculations Excel sheets to validate or to calculate. Additionally, these tools
may also serve as a gathered data validation tool.
Nomenclature
Summary Table
o
Max. Skin temperature C
o
Mean bulk temperature C
Total mass flow rate Kg/hr
Heat Duty Kcal/hr
Heat Transfer coefficient W/m.k
Max heat flux W/m2
Vapor fraction -
Average Heat Capacities Kcal/kg. oC
Average Heat Thermal Conductivities W/m. °C
Average Viscosities cP
Reynold’s number -
Prandtl Number -
The maximum heat flux density at any point in a coil can be estimated as follows:
Where
FL is the factor accounting for the effect of tube metal temperature on the radiant heat flux density
FT is the factor accounting for the effect of tube metal temperature on the radiant heat flux density
The circumferential variation factor, Fcir, is given as a function of tube spacing and coil geometry in
Figure 2. The factor given by this figure is the ratio of the maximum local heat flux density at the fully
exposed face of a tube to the average heat flux density around the tube. This figure was developed from
considerations of radiant heat transfer only.
As mentioned above, influences such as conduction around the tube and flue gas convection act to
reduce this factor. Since these influences are not included in this calculation, the calculated value will be
somewhat higher than the actual maximum heat flux density.
The longitudinal variation factor, FL, is not easy to quantify. Values between 1 .0 and 1.5 are most often
used. In a firebox that has a very uniform distribution of heat flux density, a value of 1.0 can be
appropriate. Values greater than 1.5 can be appropriate in a firebox that has an extremely uneven
distribution of heat flux density (for example, a long or a tall, narrow firebox with burners in one end
only).
FL = c10 x (1 + H / c11)
Where H = Inside firebox height, feet (mm) c10 = 1.00 and 1.10 for fuel gas and fuel oil/combination
firing, respectively c11 = 100 and 30500 for H in feet and mm, respectively
The tube metal temperature factor, FT, will be less than 1.0 near the coil outlet or in areas of maximum
tube metal temperature. It will be greater than 1.0 in areas of lower tube metal temperatures. For most
applications, the factor can be approximated as follows:
4 4
𝑇𝑔,𝑎𝑣𝑔 − 𝑇𝑡𝑚
𝐹𝑇 = ( 4 4 )
𝑇𝑔,𝑎𝑣𝑔 − 𝑇𝑡𝑚,𝑎𝑣𝑔
∗
𝑇𝑔,𝑎𝑣𝑔 : is the average flue-gas temperature, expressed in kelvins (degrees Rankine), in the radiant
section
∗
𝑇𝑡𝑚 : is the tube metal temperature, expressed in kelvins (degrees Rankine), at the point under
consideration
∗
𝑇𝑡𝑚,𝑎𝑣𝑔 : is the average tube metal temperature, expressed in kelvins (degrees Rankine), in the radiant
section.
The convective heat flux density in most parts of a radiant section is usually small compared with the
radiant heat flux density. In the shock section, however, the convective heat flux density can be
significant; it should therefore be added to the radiant heat flux density when the maximum heat flux
density in the shock section is estimated.
In addition to the heat-transfer coefficient and the maximum heat flux density, the temperature profile
of the fluid in the coil is necessary for calculating the maximum tube metal temperature in the radiant
section of the heater. This profile, which is often calculated by the heater supplier, defines the variation
of the bulk fluid temperature through the heater coil. For operation at or near design, the design profile
can be used. For operation significantly different from design, a bulk temperature profile shall be
developed. Once the bulk fluid temperature is known at any point in the coil, the maximum tube metal
temperature can be calculated as follows:
∆Tff is the temperature difference across the fluid film, expressed in degrees Celsius
∆Tcoke is the temperature difference across coke or scale, expressed in degrees Celsius
∆Ttw is the temperature difference across the tube wall, expressed in degrees Celsius
QR,max is the maximum radiant heat flux density, in W/(m2-K) for the outside surface
𝜆𝑡𝑚 is the thermal conductivity, in W/(m-K) of the tube metal. (from figure 1)
A value necessary for calculating the maximum tube metal temperature is the fluid heat-transfer
coefficient at the inside wall of the tube. Although the following correlations are extensively used and
accepted in heater design, they have inherent inaccuracies associated with all simplified correlations
that are used to describe complex relationships.
For single-phase fluids, the heat-transfer coefficient is calculated by one of the two equations below, in
which Re is the Reynolds number and Pr is the Prandtl number. No correlation is included for the heat-
transfer coefficient in laminar flow, since this flow regime is rare in process heaters. There is inadequate
information for reliably determining the inside coefficient in laminar flow for oil in tube sizes that are
normally used in process heaters.
𝜆𝑓,𝑇𝑏 𝑇𝑏 0.5
ℎ𝑣 = 0.021 ∗ ( ) ∗ 𝑅𝑒0.8 ∗ 𝑃𝑟0.4 ∗ ( )
𝐷𝑖 𝑇𝑤
𝐷𝑖 ∗ 𝑞𝑚𝐴
𝑅𝑒 =
𝜇𝑓,𝑇𝑏
𝐶𝑝 𝜇𝑓,𝑇𝑏
𝑃𝑟 =
𝜆𝑓,𝑇𝑏
For two-phase flows, the heat-transfer coefficient can be approximated using the following equation:
hf = hv * wv + hL * wL
The liquid and vapor heat-transfer coefficients, hL and hv, should be calculated using the mixed-phase
area mass flow rate but using the liquid and vapor material properties, respectively.
Iteration part
𝑇𝑏 0.5
( ) = 0.91
𝑇𝑤
1. Solve for hL and hv
2. Calculate ∆𝑇𝑓𝑓 , ∆𝑇𝑐𝑜𝑘𝑒 and ∆𝑇𝑡𝑤
3. Calculate mean tube wall temperature & oil film temperature & maximum tube temperature
Tfilm =Tb + ∆𝑇𝑓𝑓 (estimate μfilm at Tfilm)
∆𝑇𝑡𝑤
Tm = Tb + ∆𝑇𝑓𝑓 +
2
Tmax = Tb + ∆𝑇𝑓𝑓 + ∆𝑇𝑡𝑤 + ∆𝑇𝑐𝑜𝑘𝑒
0.14
𝜇𝑓,𝑇𝑏 𝑇 0.5
4. Calculate ( ) , ( 𝑏) and calculate your error percentage in assumptions
𝜇𝑓,𝑇𝑤 𝑇𝑤
0.14 0.5
𝜇 𝑇
5. Re-assume (𝜇 𝑓,𝑇𝑏 ) , (𝑇𝑏 ) accordingly until your assumption matches your calculated
𝑓,𝑇𝑤 𝑤
values!
6. Check mean tube temperature calculated against temperature assumed for tube thermal
conductivity
Calculating total mass flow rate of feed and light gases
Light gases
𝑚
𝑀𝑤𝑡 = ∑ 𝑀𝑤𝑡𝑖 ∗ 𝑦𝑖
1
𝑀𝑤𝑡𝑜𝑖𝑙 = ∑ 𝑀𝑤𝑡𝑖 ∗ 𝑦𝑖
1
𝑐𝑜𝑚𝑝𝑜𝑛𝑒𝑛𝑡 𝑚𝑜𝑙𝑠
𝑡𝑜𝑡𝑎𝑙 𝑓𝑒𝑒𝑑 𝑚𝑜𝑙 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛𝑠 =
𝑛𝑡𝑜𝑡𝑎𝑙
NOTE: from composition calculated physical properties could be calculated for both light ends and oil
fractions. Oil characterization is beyond the sope of this document thought (sorry)
Figure 1: GPSA metal thermal conductivity
Figure 2: Fcir as function of tube spacing/D
In the heater under consideration, the medium-carbon-steel tubes are in a single row against the wall. Other
aspects of the heater configuration are as follows:
The flow in the tubes is two-phase with 10 % mass vapour. Other operating conditions are as follows:
Flow rate (total liquid plus vapour) is 6.3 kg/s
Tb = 271 °C
qR.ave = 31 546 W/m2
From the inside diameter, the flow area is equal to 8.107 10 3 m2 (0.087 3 ft2). Using the total flow rate:
qmA = 6.3/(8.107 x 10-3),
qmA 777.1 kg/(m2 s).
0.0346
ℎ𝑣 = 0.021 ∗ (0.1016) ∗ (1.12 ∗ 107 )0.8 ∗ (0.486)0.33 ∗ 0.91 = 2126 W/m2.K
The temperature difference through each part of the system can now be calculated for the
fluid film:
𝑄𝑚𝑎𝑥 𝐷𝑜 66278 114.3
∆𝑇𝑓𝑓 = ( )= ( ) = 113 𝐾
ℎ𝑓 𝐷𝑖 − 2𝛿𝑐𝑜𝑘𝑒 659.9 101.6
Checking the assumed viscosity ratio, at the oil-film temperature calculated above, 271+ 113 = 384 °C, the viscosity is 1.1
mPa.s. So, for the liquid:
0.14
𝜇𝑓,𝑇𝑏 0.002 0.14
( ) =( ) = 1.820.14 = 1.09
𝜇𝑓,𝑇𝑤 0.0011
For the vapour:
𝑇𝑏 0.5 270 + 273 0.5
( ) = ( ) = 0.91
𝑇𝑤 384 + 273
Both values are close to the values assumed for the calculation of Kl and Kv. so no additional work is needed.
This is close to the temperature assumed for the tube conductivity, so no additional work is required.
Fabricated Example to explore differences between a commercial software and tool
Commercial
Parameter Tool Units
software
o
Max Skin Temperature 408.557 380 C
o
Mean Bulk temperature 399.507 363 C
Total Mass flow rate 19964.63 20,000 kg/hr
Heat Duty 819933 824150 Kcal/hr
Heat transfer coefficient 1339.157 1404.1 W/m oC
Max Heat Flux 126555.7 122200 W/m oC
o
Inlet Temperature 250 250 C
o
Outlet Temperature 290 290 C
Inlet Pressure 4 4 Kg/cm2.g
Vapor Fraction 1 1 -
Reynold's Number V/L [3124642.0, 'NA'] -
Prandtl Number V/L [0.273, 'NA'] 0.46 -
Viscosity V/L [0.012, 'NA'] 0.0221 cP
Cp 1.04 1.0479 Kcal/kg. oC
Thermal Conductivity V/L [0.188, 'NA'] 0.2239 W/m.K
Metal Thermal Conductivity 43.375 W/m.C
Fc 1.686 1.678 -
Fl 0.955 0.9551 -
Ft 1 1 -
References
This tool was developed for process engineers to quickly estimate the outcome of the safety relief valve sizing
equations on site. As a part of a larger project to develop what is similar to Carl Branan’s book “process
engineers Pocket Handbook” these tools would allow a process engineer to quickly calculate/estimate
equipment efficiencies or sizing using standardized calculations
The aim here is to take little-known data from the field (flow, pressures, temperatures, compositions..etc.) and
use it as input for a rough estimation without having to return to the office to use commercial software or
calculations Excel sheets to validate or to calculate. Additionally, these tools may also serve as a gathered data
validation tool.
Nomenclature
1. Introduction
Sizing pressure relief valves (PSV) requires a vivid understanding of the process where the safety valves are
installed, the sizing equations, and the numerous relevant standards. In this document, we’ll be exploring the
hand calculations required to preliminary size different safety valves. First, we’ll explore how to estimate the
relieving and set pressures according to standards. Secondly, we’ll explore the basic safety orifices sizing
equations (for gases, steam, or liquids) and learn how to estimate the required coefficients. Finally, we’ll take a
glance at how to estimate the required relief loads for different processes and contingencies as an introduction
to this complex matter. However, the back pressure calculations were not discussed in this document. The
backpressure used in safety calculations is the sum of superimposed backpressure and built-up back pressure.
The superimposed backpressure is the static pressure and the relief header before the valve opening, while the
built-up back pressure is the pressure increase due to the valve opening building up extra pressure to overcome
the outlet friction losses. The estimation of the built-up back pressure could be discussed later in the context of
flare systems calculations.
This document did not discuss the various safety valve types' working mechanisms, the rupture disc types and
sizing procedures, or guidelines for sizing the inlet/outlet piping of PSV. However, you may find the mentioned
topics in references and recommended reads.
2. Set/relief Pressure Calculations
Table 2.1: Relief Pressure estimation table (percentages are applied to MAWP’s gauge pressures) [1]
Contingency Installation Single-Valve Installations Multiple-Valve Installations
Set Pressure Max. Accumulated Set Pressure Max. Accumulated
(%) pressure (%) (%) pressure (%)
Non-fire only First valve 100 110 [1] 100 116 [2]
𝑃1 190.5 ∗ 𝑊 11.78𝑄√𝐺
Acalc 𝑇= ∗𝑇 𝐴= 𝐴=
𝑃𝑛 𝑛 𝑃1 𝐾𝑑 𝐾𝑠ℎ 𝐾𝑁 𝐾𝑏 𝐾𝑐 𝐾𝑑 𝐾𝑤 𝐾𝑣 𝐾𝑐 √𝑃1 − 𝑃2
𝑃𝐶𝑓 2 𝑘⁄(𝑘−1)
=[ ]
𝑃1 𝑘+1
2 𝑘+1/𝑘−1
𝐶 = 0.03948√𝑘 [ ]
𝑘+1
𝑘 1 − 𝑟 (𝑘−1)/𝑘
𝐹2 = √[ ] 𝑟 2/𝑘 [ ]
𝑘−1 1−𝑟
A: mm2
A: mm2 A: mm2
W: kg/h
W: kg/h Q: L/min
Units Q: Nm3/min
P1,P2: kPa P1,P2: kPag
P1,P2: kPa
T: Kelvin T: Kelvin
T: Kelvin
W (kg/h ) = 0.044 * M.wt * Q (Nm3/hr)
Q (Nm3/min) = Q (Nm3/hr)/60
Conversion Q (L/min) = Q (m3/hr) *16.67
Notes P (kPa) = P (kg/cm2) * 98.0665
A (cm2) = A (mm2)*0.01
A (in2) = A (cm2) * 0.155
Table 3.1: Table of single-phase relief valve sizing equations [1]
NOTE 1: The Napier coefficient KN must be considered when P1 ≥ 106.5 kg/cm2a.
NOTE 2: for gas equations
Kd = effective coefficient of discharge. For preliminary sizing, use the following values:
= 0.975 when a pressure relief valve is installed with or without a rupture disk in combination.
= 0.62 when a pressure relief valve is not installed and sizing is for rapture disk
Wrated = maximum rated flow rate through the valve
3.1 Gas PSV Sizing[1]
1. Estimate if flow is critical or sub-critical from the following equation
2 𝑘/𝑘+1
𝑃𝑐 = 𝑃1 ∗ [ ] , where
𝑘+1
Pc is the critical pressure, P1 is the upstream pressure (relieving), both in absolute units
k is the specific heats ratio for any ideal gas
2. Based on the flow condition, you can proceed with one of the following equations:
Kd: effective coefficient of discharge that should be obtained from the valve manufacturer. For preliminary sizing
estimation, a discharge coefficient of 0.65 can be used.
Kw: correction factor for back pressure. If the backpressure is atmospheric, Kw =1.0.
Balanced-bellow valves in backpressure service will require the correction factor obtained from
Figure 2 Conventional valves require no special correction.
Kv: correction factor for viscosity as determined from the following equation:
2.878 342.75 −1 170 −0.5
𝑘𝑣 = (0.9935 + + ) or 𝑘𝑣 = ( + 1)
𝑅𝑒0.5 𝑅𝑒1.5 𝑅𝑒
When a pressure relief valve is sized for viscous liquid service, it should first be sized as if it were for non-viscous
type application (i.e., Kv = 1.0) so that a preliminary required discharge area, A, can be obtained from the previous
equation. From API standard orifice sizes (Table 3.3 may be used for preliminary estimation), the next larger
orifice size should be used in determining the Reynolds number, Re, the following equation:
𝑄∗(18800∗𝐺)
𝑅𝑒 = ,
𝜇√𝐴
Inlet piping
inlet piping pressure losses due to friction should not be higher than 3% of set pressure and calculated using the
maximum rated capacity of the pressure relief device [1][3]. Exceptions were mentioned in recommended reads[II]
When two or more pressure relief devices are placed on one connection, the inlet piping internal cross-sectional
area shall be at least equal to the combined inlet areas of the pressure relief devices connected to it
Selection Criteria of PSV
1. Conventional PSVs
2. The sum of the maximum variable superimposed back pressure plus the built-up back pressure is less than
10% of the set pressure.
3. Fouling or corrosive conditions are not expected.
2. Bellows Type PSVs
1. The sum of the variable superimposed back pressure plus built-up back pressure exceeds 10% of set
pressure.
2. Fouling or corrosive conditions are expected and protection cannot be afforded by using alternative
materials or devices.
More on the advantages and disadvantages of different PSV types can be found at recommended reads [II] [III]
4. Examining possible cases (Contingencies)
Estimating the required relief loads could be a grueling task that requires extensive heat and mass balance
calculations at relief conditions. It could be best to use simulation tools for this one. However, this section
aims to familiarize readers with the basis and variables needed for different scenarios. Moreover, manual
calculations could be an alternative option where minimum data are available for simulation in some cases.
All causes of overpressure, or contingencies, must be evaluated for each PSV installation in terms of the
pressures generated and the rates at which fluids must be relieved. Causes of overpressure in process
equipment can range from a single event to a complex combination of events. The basis for calculating a
valve size follows calculations of valid contingencies. The contingency that requires the largest effective area
dictates the size of the PSV.
More on “guidelines on estimating relief loads for different contingencies” can be found on recommended
reads [V][I][III]
Table 4.1: Guidelines on estimating relief loads for different contingencies [2]
Item Pressure Relief Device Pressure Relief Device
No. Condition (Liquid Relief) (Vapor Relief)*
1 Closed outlets on vessels Maximum liquid pump- Total incoming steam and vapor plus that
in rate generated therein at relieving conditions
2 Cooling water failure to condenser — Total vapor to condenser at relieving conditions
3 Top-tower reflux failure — Total incoming steam and vapor plus that
generated therein at relieving conditions less
vapor condensed by side-stream reflux
4 Side-stream reflux failure — Difference between vapor entering and leaving
section at relieving conditions
5 Lean oil failure to absorber — None, normally
6 Accumulation of noncondensables — Same effect in towers as found for Item 2; in
other vessels, same effect as found for Item 1
7 Entrance of highly volatile — For towers, usually not predictable
material: Water in hot oil or Light
hydrocarbons in hot oil
8 Overfilling storage or surge vessel Maximum liquid —
Pump-in rate
9 Failure of automatic controls — Must be analyzed on a case-by-case basis
10 Abnormal heat or vapor input — Estimated maximum vapor generation including
non-condensables from overheating
11 Split exchanger tube — Steam or vapor entering from twice the cross-
sectional area of one tube; also same effects
found in Item 7 for exchangers
12 Internal Explosions — Not controlled by conventional relief devices
but by avoidance of circumstances
13 Chemical reaction — Estimated vapor generation from both normal
and uncontrolled conditions
14 Hydraulic expansion Use Liquid expansion —
1. Cold fluid shut in equation (table 4.7)
2. Line outside process area
3. Shut in
15 Exterior fire — Estimate by external fire equations (table 4.2)
16 Power failure (steam, electric, or — Study the installation to determine the effect of
other) power failure; size relief valve for the worst
condition that can occur
1. Fractionators — All pumps could be down, with the result that
reflux and cooling water would fail
2. Reactors — Consider failure of agitation or stirring, quench
or retarding steam; size valves for vapor
generation from a runaway reaction
3. Air-cooled exchangers — Fans would fail; size valves for the difference
between normal and emergency duty
4. Surge vessels Maximum liquid inlet —
rate
* Considerations may be given to the suppression of vapor production as the result of the device’s relieving
pressure being above operating pressure, assuming constant heat input. (Procedures for sizing pressure relief devices
are presented in Section 4 of API-RP-520.)
After the specifications are determined, the next activity is to calculate a preliminary effective discharge area.
Figure 4.1 shows the inputs used to determine the basis for calculating the size of a PSV.
An overpressure problem does not occur on If the design pressure for the compressor discharge system is
the discharge line of a centrifugal pump higher than the pressure of surge point at maximum speed,
because the centrifugal pump discharge overpressure does not occur. If the design pressure is lower
system normally has a design pressure equal than that, overpressure protection should be considered. In
to or higher than the pump shut-off this case, the relieving load should be the flow rate (FD) at the
pressure. head equivalent to the design pressure (PD) at maximum
However, when the design pressure is set speed or should be the anti-surge flow (FS) at maximum
lower than the pump shut-off pressure, an speed, whichever is greater. That value is usually obtained
overpressure problem will occur, and the from such compressor performance curve as shown below:
flow rate at the head, equivalent to the set
pressure of the PSV minus the maximum
suction pressure, should be read from the
pump performance curve for the relieving
load.
The relieving load should be determined
based on the pump’s maximum speed when
the steam turbine driver is provided.
If the vapor pressure of the cold medium at the When the outlet of the cold side is blocked off, and the hot
inlet temperature of hot side is more than 1.3 medium continues to flow, overpressure will
times the design pressure of the cold side, occur, and the relieving load should be calculated based on the
overpressure protection due to vaporization is vaporizer heat duty at the relieving condition obtained by the
needed. If it is assumed that the outlet of the cold procedure shown in Table 4.3
side is blocked off, and the hot medium continues
to flow, the relieving load on the cold side should
be determined as follows :
𝑄 (𝑇1 − 𝑡𝑏𝑝 )
𝑊 = .
𝜆 (𝑇1 − 𝑡𝑎𝑣 ) Steam Vaporizer
𝑡1 + 𝑡2 (1) Inlet and outlet temperatures of process fluid to/from
𝑡𝑎𝑣 =
2 reboiler at relieving pressure will increase, because
W : relieving load (kg/hr) the boiling temperature of process fluid increases, when the
Q : normal heat exchanger duty (kcal/hr) operating pressure reaches the relieving pressure.
λ : latent heat of vaporization at tbp (kcal/kg) (2) When heating steam is supplied under steam flow control,
T1: hot side inlet temperature (oC) the saturated temperature of steam at the
Tbp: cold side inlet temperature (oC) pressure of supply header should be applied to ΔT calculation,
tav: average cold side temperature (normal because the flow control valve tend to open to maintain the
operation) (oC) steam flow at a constant. If the other control system is applied
t1/t2: cold side inlet/outlet temperature (oC) to the steam supply, the normal operating temperature of
steam may be used.
(3) Ucalc (calculated overall heat transfer coefficient ) should be
used for recalculation.
Table 4.7: Relief loads guidelines for Liquid Thermal Expansion and Solar radiation
Inadvertent Control Valve Opening
Flashing Liquid [4]
NOTE 2: Shortcut mentioned in ref [8] were used in the tool (Cv =10*d2), d in inches
Wetted Area Calculations [5]
To determine vaporization rate, the surface
area wetted by a vessel's internal liquid and is Semi-Ellipsoidal/Elliptical Heads
up to 7.6 m (25 feet) above grade, denoted as Or Torispherical/Dished Heads
"fire elevation", needs to be considered. The
term "grade" usually refers to ground level,
but may be at any level at which a sizable fire
could be sustained.
z
D
Vessels
Cylindrical part
h
Vertical
𝐴𝑤
𝐷2 ℎ
= 𝜋 (( − 0.5) 𝐵 + 1
8 𝐷
ℎ
L 1 4𝜀 ( − 0.5) + 𝐵
+ ln( 𝐷 ))
D 4𝜀 2 − √3
2
ℎ
√
𝐵 = 1 + 12 ( − 0.5)
𝐷
4𝑧 2
𝜀 = √1 −
𝐷2
D
z: inside dish depth
h
𝜀 : Eccentricity of elliptical heads
L
𝐴𝑤 = 𝐿𝐷 cos −1 (1 − (2ℎ/𝐷))
Heads
Hemispherical Heads
𝐷
𝐴𝑤 = 𝜋ℎ
2
Shortcut Calculations for various wetted Areas [4]
Output
Relieving pressure kPa / kg/cm2 670
Z - 0.9 0.98
Critical Pressure kPa / kg/cm2 392 -
Cp/Cv - 1.11 1.05
Kb - 1 1
Kc - 1 1
kd - 0.975 0.975
Acalculated cm2/ in2 36.98 / 5.73 39.5 / 6.12
Aselected cm2/ in2 41.16 / 6.38 41.16 / 6.38
Orifice Designation P P
1. Estimating Relief Pressure and constants
Relief pressure for one relief valve, no rupture disc, and non-fire case will equal 1.1 the equipment
design pressure
Pr = 1.1 * 517 = 670 kPa (6.8321 kg/cm2)
Kb = 1 as backpressure is atmospheric
Kc = 1, as there’s no rupture disc installed
Kd = 0.975
1. Checking whether the flow is critical
2 𝑘⁄(𝑘−1) 2 1.11⁄(1.11−1)
𝑃𝐶𝑓 = 𝑃1 ∗ [ ] = 670 ∗ [ ] = 670 ∗ 0.5825 = 390.33
𝑘+1 1.11 + 1
The PRV sizing is based on the critical flow equation since the backpressure (0 kPag) is less than the
critical flow pressure (291 kPag).
2 (𝑘+1)/(𝑘−1) 2 (1.11+1)/(1.11−1)
𝐶 = 0.03948√𝑘 [ ] =0.03948*√1.11 [ ] = 0.03948 ∗ 0.63045 =
𝑘+1 1.11+1
0.02489
Using The C graph provided earlier (Figure 2.7), C = 0.0249 for Cp/Cv of 1.11
𝑊 ∗ √𝑇𝑍 24270 348 ∗ 0.9
𝐴= = √ = 3698 𝑚𝑚2
𝐶𝐾𝑑 𝑃1 𝐾𝑏 𝐾𝑐 √𝑀 0.0249 ∗ 0.975 ∗ 1 ∗ 1 ∗ 670 51
See API 526 for the selection of the proper orifice size. API 526 provides standard effective
orifice areas in terms of letter designations. For this example, a “P” size orifice should be selected since it
has an effective orifice area of 6.38 in.2 (4116 mm2)
Step 1: Calculate Relief Pressure Step 2: Input mass flow rate and composition
Step 3: Input relief/back conditions (T&P) Step 4: Download your calculations table!
API RP 520 Example No. 4
Input Units API RP 520 Part I Tool
Mass flow Kg/hr 69615 69615
Relieving Temperature C 427 427
Design pressure of equipment kPag / kg/cm2g 11032 112.5
Backpressure kPag / kg/cm2g 101.325 1.033
Overpressure % 10 10
Output
Relieving pressure kPa / kg/cm2 12236 (124.77) 124.77
Kb - 1 1
Kc - 1 1
kd - 0.975 0.975
KN - 1.01 1.01
KSH - 0.855 0.862
Acalculated cm2/ in2 12.87 / 1.995 12.75 / 1.97
Aselected cm2/ in2 18.41 / 2.853 18.41 / 2.853
Orifice Designation L L
1. Estimating Relief Pressure and constants
Relief pressure for one relief valve, no rupture disc, and non-fire case will equal 1.1 the equipment
design pressure
Pr = 1.1 * 11032 + 101.325 = 12236.5 kPa (124.77 kg/cm2)
Kb = 1 for conventional valve discharging to atmosphere
Kc = 1, as there’s no rupture disc installed
Kd = 0.975
KSH : 0.855 from tables
KN : P1 (124.77 kg/cm2.a) is > 106.5 kg/cm2.a
KN = (0.1906 *1774.7 - 1000)/(0.2292 *1774.7- 1061) = 1.01
orifice areas in terms of letter designations. For this example, a “L” size orifice should be selected since it
Step 3: Input relief/back conditions (T&P) Step 4: Download your calculations table!
API RP 520 Example No. 5
Input Units API RP 520 Part I Tool
Mass flow L/min 6814 368010 (kg/hr)
Specific Gravity 0.9 0.9
Viscosity Saybolt / cP 2000 / 388.5 388.5
Design pressure of equipment kPag / kg/cm2g 1724 17.58
Backpressure kPag / kg/cm2g 345 3.52 (4.55 abs)
Overpressure % 10 10
Output
Relieving pressure kPag / kg/cm2 1896 (19.33 kg/cm2.g) 20.37 (19.337 gauge)
Kc - 1.0 1
Kd - 0.65 0.65
Kw - 1 0.97
Kv - 1 (initially) / 0.982 0.99
Acalculated cm2/ in2 31.22 / 4.84 30.71 / 4.83
Aselected cm2/ in2 41.16 / 6.38 41.16 / 6.38
Orifice Designation P P
1. Estimating Relief Pressure and constants
Relief pressure for one relief valve, no rupture disc, and non-fire case will equal 1.1 the equipment
design pressure
Pr = 1.1 * 1724 = 1896.4 kPa (kg/cm2)
Kw = 1 as backpressure is atmospheric
Kc = 1 as there’s no rupture disc installed
Kd= 0.975
Kv: initially assumed 1
2. Calculating Orifice Area
See API 526 for the selection of the proper orifice size. API 526 provides standard effective
orifice areas in terms of letter designations. For this example, a “P” size orifice should be selected since it
has an effective orifice area of 6.38 in.2 (4116 mm2)
3. Check Kv and recalculate Acalculated
𝑄 ∗ (18800 ∗ 𝐺) 6814*18800*0.9
𝑅𝑒 = = = 4625
𝜇√𝐴 388.5√4116
−0.5
170
𝑘𝑣 = ( + 1) = 0.982
𝑅𝑒
𝐴𝑅 3066
𝐴= = = 3122 𝑚𝑚2
𝐾𝑣 0.982
Step 1: Calculate Relief Pressure Step 2: Input mass flow rate and composition
Step 3: Input relief/back conditions (T&P) Step 4: Download your calculations table!
Fire Case Example: Benzene in a Horizontal Vessel
Input Units Reference [1] Tool
Units Converted
Fluid Benzene Benzene
Molecular weight - 78.11 78.11
Relieving Temperature K 37.78 + 273.15 [Note] 250 + 273.15 [Note 2]
Design pressure of equipment kPag / kg/cm2.g 1379 14.06 kg/cm2.g (1379 kPa)
Backpressure kPa / kg/cm2 Atmospheric Atmospheric
Overpressure % 21 21
Vessel Geometry data Bare Vessel Bare Vessel
Horizontal + Horizontal +
spherical heads Cylindrical heads
Diameter m 4.572 4.572
T-T Length m 9.144 9.144
Level (Normal liquid level) m 3.734 3.734
Height above grade m 4.572 4.572
Output
Relieving pressure kPa / kg/cm2 1769.88 18.04 kg/cm2.g (1769.7 kPa)
ΔH kCal/kg 95.62(input) 146 [Note 1]
Wetted Area (Aw) m2 83.7 121
Mass flow Kg/hr 14665 13020.15
Z - 1 (input) 0.81
Critical Pressure kPa / kg/cm2
Flow condition - Critical Critical
Cp/Cv - 1.12 (input) 0.97
Kb - 1 1
Kc - 1 1
kd - 0.975 0.975
C
Acalculated cm2/ in2 6.787 / 1.052 in2 8.20495 / 1.27177 in2
Aselected cm2/ in2 8.303 / 1.287 in2 8.303 / 1.287 in2
Orifice Designation J (nozzles 2”/3”) J (nozzles 2”/3”)
Notes:
1. Heat of vaporization is higher in the tool as a result of adding the heat required to bring the fluid
from the subcooled phase to saturation temperature
2. A temperature of 250 was assumed to obtain benzene at the vapor phase; assuming relived and
vaporized benzene at 37.78 C & 18 kg/cm2 is illogical. This, however, will increase Acalculaed.
3. Estimate Vessel wetted area
2ℎ
𝐴𝑤 (𝐶𝑦𝑙𝑖𝑛𝑑𝑟𝑖𝑐𝑎𝑙 𝑝𝑎𝑟𝑡) = 𝐿𝐷 cos −1 (1 − ( ))
𝐷
3.03
= 9.1 ∗ 4.57 𝑐𝑜𝑠 −1 (1 − (2 ∗ )) = 587 ∗ cos −1 (−0.326) = 79.14 𝑚2
4.57
𝐴𝑤 (2 ℎ𝑒𝑎𝑑𝑠) = 𝜋ℎ𝐷 = 3.14 ∗ 3.03 ∗ 4.057 = 38.6 𝑚2
𝐴𝑤 𝑡𝑜𝑡𝑎𝑙 = 79.147 + 38.6 = 117.74 𝑚2
4. Estimate Relief Load
Q = qoFA0.82
Relief pressure for one relief valve, no rupture disc, and fire case will equal 1.21 the equipment design
pressure
Pr = 1.21 * 1379 + 101.325 = 1770 kPa (18.05 kg/cm2)
Kb = 1 as backpressure is atmospheric
Kc = 1, as there’s no rupture disc installed
Kd = 0.975
4. Checking whether the flow is critical
2 𝑘⁄(𝑘−1) 2 1.12⁄(1.12−1)
𝑃𝐶𝑓 = 𝑃1 ∗ [ ] = 1770 ∗ [ ] = 1770 ∗ 0.5805 = 1027.5 𝑘𝑃𝑎
𝑘+1 1.12 + 1
The PRV sizing is based on the critical flow equation since the backpressure (0 kPag) is less than the
critical flow pressure (1027 kPag).
5. Calculating Orifice Area
2 (𝑘+1)/(𝑘−1) 2 (1.12+1)/(1.12−1)
𝐶 = 0.03948√𝑘 [ ] =0.03948*√1.12 [ ] = 0.03948 ∗ 0.4 =
𝑘+1 1.12+1
0.0249
Using The C graph provided earlier (Figure 2.7), C = 0.025 for Cp/Cv of 1.12
See API 526 for the selection of the proper orifice size. API 526 provides standard effective orifice areas
in terms of letter designations. For this example, a “J” size orifice should be selected since it has an
effective orifice area of 1.287 in.2 (830 mm2)
Step 1: Calculate Relief Pressure Step 2: Select estimate relief loads and choose “External Fire
due to vaporization”
Step 3: input relieved stream composition and conditions Step 4: Input wetted Area required inputs
Step 5: check your results and calculations’ matrix Step 6: go to “3. Gas or vapor” and use Calculated Relief Load
Step 7: Input mass flow rate and composition Step 8: Download your calculations table!
Maximum allowable working pressure (MAWP) is the maximum gauge pressure permissible at the top of
a completed vessel in its operating position for a designated temperature. The pressure is based on
calculations for each element in a vessel using nominal thickness, exclusive of additional metal thickness
allowed for corrosion and loading other than pressure. The maximum allowable working pressure is the
basis for the pressure setting of the pressure relief devices that protect the vessel.
Design gauge pressure refers to at least the most severe conditions of coincident temperature and
pressure expected during operation. This pressure may be used in place of the maximum allowable
working pressure in all cases where the MAWP has not been established. The design pressure is equal to
or less than the MAWP.
Accumulation is the pressure increase over the maximum allowable working pressure of the vessel
during discharge through the pressure relief device, expressed in pressure units or as a percent.
Maximum allowable accumulations are established by applicable codes for operating and fire
contingencies.
Overpressure is the pressure increase over the set pressure of the pressure relief device, expressed in
pressure units or as a percent. It is the same as accumulation when the relieving device is set at the
maximum allowable working pressure of the vessel.
Rated relieving capacity is the portion of the measured relieving capacity permitted by the applicable
code regulation to be used as a basis for the application of a pressure relief device. Stamped capacity is
the rated relieving capacity that appears on the device nameplate.
The stamped capacity is based on the set pressure plus the allowable overpressure for compressible
fluids and the differential pressure for incompressible fluids. The stamped capacity shall not exceed 90%
of the average capacity of the valves tested.[ASME VIII Div. 1 UG-13 1 (d)(l)]
Set pressure is the inlet gauge pressure at which the pressure relief device is set to open under service
conditions.
Pressure differential is the difference between the set pressure of the pressure relief device and the
operating pressure of the protected vessel.
Cold differential test pressure is the pressure at which the pressure relief valve is adjusted to open on the
test stand. The cold differential test pressure includes corrections for the service conditions of back
pressure or temperature or both.
Back pressure is the pressure that exists at the outlet of pressure relief devices as a result of the pressure
in the discharge system. It is the sum of the superimposed and built-up back pressure.
Superimposed back pressure: The static back pressure that exists at the outlet of a pressure relief device
at the time the device is required to operate. It is result of pressure in the discharge system corning from
other sources and may be constant or variable. Built-up back pressure: The built-up back pressure is the
increase in pressure in the discharge header that develops as a result of flow after the pressure relief
device or devices open.
The built-up back pressure is caused by flow from the particular device and others, if any, which
simultaneously discharge into the disposal system. This type of back pressure is variable. The built-up
back pressure shall be less than the allowable back pressure, 10% of the set pressure for the
conventional type or 50% of the set pressure for the balanced-bellows type.
Constant Back Pressure: The static back pressure that exists under normal operation where no relief
device is operated; i.e. constant back pressure of the superimposed back pressure. The constant back
pressure is used to determine the spring set pressure of conventional type pressure relief valve as a
difference between the set pressure and the constant back pressure. Therefore, the opening pressure
will vary depend on the built-up back pressure. The opening pressure of balanced-bellows type is the set
pressure of the valve and is independent of any back pressure.
Blowdown is the difference between the set pressure and the closing pressure of a pressure relief valve,
expressed as a percent of the set pressure or in pressure units.
Opening pressures the value of increasing inlet static pressure at which there is a measurable lift of the
disc or at which discharge of the fluid becomes continuous.
Closing pressure is the value of decreasing inlet static pressure at which the valve disc reestablishes
contact with the seat or at which lift becomes zero.
Simmer is the audible or visible escape of compressible fluid between the seat and disc at an inlet static
pressure below the set pressure and at no measurable capacity
Leak-test pressure is the specified inlet static pressure at which a seat leak test is performed. The term
relieving conditions is used to indicate the inlet pressure and temperature on a pressure relief device at a
specific overpressure. The relieving pressure is equal to the valve set pressure plus the overpressure. The
temperature of the flowing fluid at relieving conditions may be higher or lower than the operating
temperature.
Chatter refers to the motion that causes the disc to contact the seat and damage the valve and
associated piping. Chattering may result in lowered capacity and damage to the seating surfaces.
Flutter refers to the abnormally rapid reciprocating motion of the movable parts of a pressure relief valve
in which the disc does not contact the seat
Atmospheric discharge is the release of vapors and gases from pressure relief or depressing devices to
the atmosphere.
Flare system is a means for the safe disposal of waste gasses by closed pipeline and combustion system.
With an elevated mre the combustion is carried out at the top of a pipe or stack where the burner and
igniter are located. A ground flare is similarly equipped except that combustion is carried out at or near
ground level. A burn pit differs from a flare in that it is normally designed to handle both liquids and
vapors.
Vent stack is the elevated vertical termination of a disposal system that discharges vapors into the
atmosphere without combustion or chemical conversion of the relieved fluid.
References
[1] API RP 520 Part I, 10th edition.
[2] Engineering Encyclopedia: Aramco desktop standards: sizing and selecting Pressure relief valves
[3] Pressure relief devices rev. 5, JGC standard practice 2002 JGS 210-120-1-61E
[4] Pressure Relieving design rev. 5, JGC standard practice 2008 JGS 210-120-1-40E
[5] Volume and Wetted Area of Partially Filled Horizontal Vessels, Link:
https://neutrium.net/equipment/volume-and-wetted-area-of-partially-filled-horizontal-vessels/
[6] Sizing Pressure-Relief Devices, Daniel A. Crowl and Scott A. Tipler. (AIChE)
[7] Elements of chemical process engineering chapter 1, Basic Process Engineering principles
[8] Chemical process engineering Vol. 1, Chapter 5 P. 246 By Kayode Coker
This tool was developed for process engineers to es mate a control valve Cv/Size on site
quickly. As a part of a larger project to develop what is similar to Carl Branan’s book “process
engineers Pocket Handbook” these tools would allow a process engineer to quickly
calculate/es mate equipment efficiencies or sizing using standardized calcula ons
The aim here is to take li le-known data from the field (flow, pressures, temperatures,
composi ons..etc.) and use it as input for a rough es ma on without having to return to the
office to use commercial so ware or calcula ons Excel sheets to validate or to calculate.
Addi onally, these tools may also serve as a gathered data valida on tool.
Nomenclature
Summary Table
Calculation - D inlet assumed in.
Inlet pressure kg/cm².a D outlet assumed in.
Flow Rate Nm³/hr FL -
Outlet pressure kg/cm².a Fd -
Inlet temperature ºC xT -
Molecular weight - FP -
Compressibility factor - xTP -
Avergage k (Cp/Cv) - FR -
Viscosity cP choked -
Shortcut Cv gpm/psi Y -
Kv m³.hr-1/bar Rev -
Cv gpm/psi laminar -
Control Valve size in.
Control Valve Specifications sheet Table of matching valves
Control valve sizing
1. Rules of thumb for Estimating △Psizing and QSizing & Other Guidelines
2. Different options for Calculating Valve Cv
3. Selecting Valve Size and calculating the opening percentage (or lift %).
4. General guidelines when selecting a control valve
1. Rules of thumb for Estimating △Psizing and QSizing[3][4]
“The sum of the following pressure drops at maximum flow may be used for this purpose:
1. 20% of the friction drop in the circuit (excluding the valve). (A circuit generally
includes all equipment between the discharge of the pump, compressor, or vessel
and the next point downstream of which pressure is controlled. In most cases this
latter point is a vessel)
2. 10% of the static pressure of the vessel into which the circuit discharges up to the
pressure of 200 psig (14.06 kg/cm2.g),
20 psig from 200 psig to 400 psig (1.406 from 14.06 to 28.12 kg/cm 2.g)
and 5% above 400 psig.” [Note]
“In order to ensure that the valve is in a controlling posi on at the maximum flow rate, the
valve Cv is the maximum process value determined above divided by 0.8. The reasons for
using this factor are that:
1. It is not desirable to have the valve fully open at maximum flow since it’s not then in
a controlling position.
2. The valves supplied by a single manufacturer often vary as much as 10-20% in Cv.
3. Allowance must be made for pressure drop, flow rate, etc. values which differ from
design.” [1]
Another Rule of thumb is to use the greater of the following [2]
1. QSizing = 1.3* QNormal or 2*QNormal [4]
2. QSizing = 1.1* QMaximum
For saturated steam keep the control valve outlet velocity below 0.25 mach. Branan Also suggests
keeping the inlet velocity below 300 /s for 2” pipes and below 200 /s for larger sizes. [2]
The above figures are guidelines for typical applica ons.
In general, smaller-sized valves handle slightly higher veloci es and large valves handle lower
veloci es. Special applica ons have par cular velocity requirements;
a few of which are provided below.
Liquid applica ons – where the fluid temperature is close to the satura on point – should be
limited to 30 feet per second (9.1 m/s) to avoid reducing the fluid pressure
below the vapor pressure. This is also an appropriate limit for applica ons designed to pass the
full flow rate with a minimum pressure drop across the valve.
Valves in cavita ng service should also be limited to 30 feet per second (9.1 m/s) to minimize
damage to the downstream piping. This will also localize the pressure recovery
which causes cavita on immediately downstream from the vena contracta.
In flashing services, veloci es become much higher due to the increase in volume resul ng from
vapor forma on.
For most applica ons, it is important to keep veloci es below 500 feet per second. (152.4 m/s)
Expanded outlet-style valves – such as the Mark One-X – help to control outlet veloci es on such
applica ons. Erosion damage can be limited by using chrome-moly body material and hardened
trim. On smaller valve applica ons that remain closed for most of the me – such as heater drain
valves – higher veloci es of 800 to 1500 feet per second (243.84 - 457.2 m/s) may be acceptable
with appropriate materials.
Gas applica ons where special noise a enua on trim is used should be limited to approximately
0.33 Mach. In addi on, pipe veloci es downstream from the valve are cri cal to the overall noise
level. Experimenta on has shown that veloci es around 0.5 Mach can create substan al noise
even in a straight pipe. The addi on of a control valve to the line will increase the turbulence
downstream, resul ng in even higher noise levels.
Noise
Typical Noise limita on: 85-90 dB(A) [4][5]
Cau on: Noise levels over 110 dBA may cause vibra on in valves/piping resul ng in equipment
damage. [8]
P in bar
Non-choked flow (with fi ngs) (ΔP < (FLP/Fp) 2(P1-FFPv) )
ρ in kg/m3
𝑸 𝝆𝟏 ⁄𝝆𝒐
𝐶𝑣 = Eq.2 ρ0 : 1000 kg/m3
𝑵𝟏 𝑭 𝑷 ∆𝑷
𝑸 𝝆𝟏 ⁄𝝆𝒐
𝐶𝑣 = Eq.3
𝑵𝟏 𝑭𝑳 𝑷𝟏 − 𝑭𝑭 𝑷𝒗
𝑸 𝝆 𝟏 ⁄𝝆 𝒐
𝐶𝑣 = Eq.4
𝑵𝟏 𝑭𝑳𝑷 𝑷𝟏 − 𝑭𝑭 𝑷𝒗
𝑸 𝝆𝟏 ⁄𝝆𝒐
𝐶𝑣 = Eq.5
𝑵𝟏 𝑭 𝑹 ∆𝑷
Specify the Variables Required to Size the Valve[5] The values of FL for the globe, v-notch ball, and bu erfly valves can
• Process fluid (e.g. oil and water); be extracted from Table 2. Ball valves with standard design are not
• Service condi ons and fluid physical proper es such as volume rate suitable for flow control; however, a v-notch ball valve with a unique
of flow (q), mass rate of flow (w), Upstream and downstream v-shaped hole in the ball allows precise flow control. Furthermore,
pressure values (P1 and P2), pressure drop across the valve (ΔP = P1 − bu erfly valves can be used for fluid control or thro ling in addi on
P2), upstream temperature (T1), Specific gravity (SG), Viscosity (µ), to fluid isola on.
Vapor pressure (Pv), Cri cal pressure (Pc), and solid%. If the valve is installed with fi ngs such as reducers a ached to it,
• Constants FL & Fd: A er selec ng required Valve type, trim type and the combined liquid pressure recovery factor and piping geometry
flow direc on, these valves can be es mated using Table 2. factor of a valve with a ached fi ngs known as parameter FLP is
applied instead of FL in Eq. (35) for maximum flow rate (qmax)
Determine Piping Geometry Factor (FP)
calcula on. FLP is calculated as per Eq. (30).
FP is a correc on factor that accounts for the pressure losses due to
piping fi ngs that are used to change the direc on or size of the Maximum Pressure Drop Calcula on for Valves Without Any
piping or take branches from the headers such as tees, elbows, and A ached Fi ngs
reducers. If fi ngs are used before or a er the valve, the effect of ΔPmax = FL2(P1-FFPv)
piping geometry must be considered for the valve sizing. However, if Maximum Pressure Drop Calcula on for Valves With A ached
no fi ng is a ached to the valve, FP is equal to one and does not Fi ngs
impact the valve sizing process. Equa on (36) shows how to calculate ΔPmax = (FLP/Fp) 2(P1-FFPv)
the piping geometry.
a) Solve for Flow Coefficient (valve diameter known)
Determine the Maximum Flow Rate (qmax) and Maximum Pressure Drop
Equa on (1 & 2) is used to calculate the modified flow efficient (Cv)
(ΔPmax)
value if ΔPmax > ΔP meaning that the flow will not be choked.
As its name implies, the maximum flow rate, also called limi ng flow
Equa on (3 & 4) applies to adjusted flow coefficient calcula on if
rate (qmax) or choked flow calculated from Eq. (35), cannot be
ΔPmax < ΔP when the flow will be choked.
increased even by increasing the pressure drop across the valve. For
1. Select FL from Table 2
example, flashing the vapor from liquids when the sta c valve
2. Calculate FF from Eq. 31
pressure drops below the liquid vapor pressure can cause a
3. Calculate ∆Pmax and determine if flow is choked or not.
maximum flow rate. It is important to calculate the maximum
4. Calculate Ci by Eq. 1 (non-choked) or Eq. 3 (Choked)
allowable pressure drop (ΔPmax) for selec ng the correct size of the
5. Calculate Rev using Eq. 28, Obtain Fd from Table 2.
valve if it is possible for the choked flow to develop. The calculated
6. If Rev > 10000, flow is turbulent and If valve size = pipe size
ΔPmax value is compared with the actual pressure drop across the
Use Ci previously Calculated. If not, Calculate Fp and
valve specified by the process engineer, and the lesser of these two is
subsequently its Cv using the iteration method described in
applied to the valve sizing process. If ΔPmax > ΔP, the flow will not be
flowchart.
choked and there is no need to determine the maximum flow rate as
7. If Rev < 10000, flow is non-turbulent. Calculate FR and its
per step 4. However,
subsequently its Cv using the iteration method described in
if a choked flow condi on exists (ΔPmax > ΔP), then step 5 for sizing
flowchart.
valves for liquids shall be changed by replacing the actual service
pressure differen al (ΔP = P1 − P2) with ΔPmax.
Select the Correct Valve Size
Maximum Flow Rate (qmax)
As a general rule, the smallest valve with a flow coefficient larger
𝑃 −𝐹 𝑃 Eq.35 than that calculated based on either Eq. (1 & 2) or Eq.(3 & 4),
𝑄 =𝑁𝐹𝐶 according to the informa on provided in Fisher representa ve Table
𝑆𝐺
shall be selected. Select the desired flow Characteris c and calculate
the es mated opening for the chosen valve.
The liquid cri cal pressure ra o factor (FF) can be calculated by using
Eq. (31), or Figure 5 for water. When a liquid pressure drops below its
b) Alterna ve op on for valve sizing (valve diameter isn’t known):
vapor pressure, the fluid state alters from a liquid to gas, that is called
Use the simplified equa ons described earlier then use the sizing
choked flow or vena contracta pressure condi on. FF is the ra o of
equa on in the following sec on (Another method for es ma ng
the vena contracta pressure at choked flow condi on to the vapor
control valve size (if tables not available))
pressure of the liquid at inlet temperature.
Check Noise levels, velocity and Cavita on index
𝑃 −𝑃
𝐶𝑎𝑣𝑖𝑡𝑎𝑡𝑖𝑜𝑛 𝑖𝑛𝑑𝑒𝑥 (𝛿) =
𝑃 −𝑃
Cavita on index: Ra o of pressure drop across the valve divided by
the difference between the inlet pressure and the vapor pressure of
the liquid. Valve selec on to ensure opera on above the cavita on
point is the primary use.
Turbulent flow
Non-choked flow (without fi ngs) (ΔP/P1 < Fϒ.xT )
𝑾 N6 : 27.3
𝐶𝑣 = Eq.6 N7 : 417
𝑵𝟔 𝒀 𝒙𝑷𝟏 𝝆𝟏
N8 : 94.8
N9 : 2120
𝑾 𝑻𝟏 𝒁 W in kg/hr
𝐶𝑣 = Eq.7
𝑵𝟖 𝑷𝟏 𝒀 𝒙𝑴 Q in Nm3/hr
P in bar
𝑸 𝑴𝑻𝟏 𝒁 Density in kg/m3
𝐶𝑣 = Eq.8a T in kelvin
𝑵𝟗 𝑷𝟏 𝒀 𝒙
𝑸 𝑮𝒈𝑻𝟏 𝒁
𝐶𝑣 = Eq.8b
𝑵𝟕 𝑷𝟏 𝒀 𝒙
𝑾
𝐶𝑣 = Eq.9
𝑵𝟔 𝑭𝑷 𝒀 𝒙𝑷𝟏 𝝆𝟏
𝑾 𝑻𝟏 𝒁
𝐶𝑣 = Eq.10
𝑵𝟖 𝑭𝑷 𝑷𝟏 𝒀 𝒙𝑴
𝑸 𝑴𝑻𝟏 𝒁
𝐶𝑣 = Eq.11a
𝑵𝟗 𝑭𝑷 𝑷𝟏 𝒀 𝒙
𝑸 𝑮𝒈𝑻𝟏 𝒁
𝐶𝑣 = Eq.11b
𝑵𝟕 𝑭𝑷 𝑷𝟏 𝒀 𝒙
Choked flow (without fi ngs) (ΔP/P1> Fϒ.xT )
𝑾
𝐶𝑣 = Eq.12
𝟎. 𝟔𝟔𝟕 ∗ 𝑵𝟔 𝐅ϒ 𝒙𝑻 𝑷𝟏 𝝆𝟏
𝑾 𝑻𝟏 𝒁
𝐶𝑣 = Eq.13
𝟎. 𝟔𝟔𝟕 ∗ 𝑵𝟖 𝑷𝟏 𝐅ϒ 𝒙𝑻 𝑴
𝑸 𝑴𝑻𝟏 𝒁
𝐶𝑣 = Eq.14a
𝟎. 𝟔𝟔𝟕 ∗ 𝑵𝟗 𝑷𝟏 𝐅ϒ 𝒙𝑻
𝑸 𝑮𝒈𝑻𝟏 𝒁
𝐶𝑣 = Eq.14b
𝟎. 𝟔𝟔𝟕 ∗ 𝑵𝟕 𝑷𝟏 𝐅ϒ 𝒙𝑻
𝑾
𝐶𝑣 = Eq.15
𝟎. 𝟔𝟔𝟕 ∗ 𝑭𝑷 𝑵𝟔 𝐅ϒ 𝒙𝑻 𝑷𝟏 𝝆𝟏
𝑾 𝑻𝟏 𝒁
𝐶𝑣 = Eq.16
𝟎. 𝟔𝟔𝟕 ∗ 𝑭𝑷 𝑵𝟖 𝑷𝟏 𝐅ϒ 𝒙𝑻 𝑴
𝑸 𝑴𝑻𝟏 𝒁
𝐶𝑣 = Eq.17a
𝟎. 𝟔𝟔𝟕 ∗ 𝑭𝑷 𝑵𝟗 𝑷𝟏 𝐅ϒ 𝒙𝑻
𝑸 𝑮𝒈𝑻𝟏 𝒁
𝐶𝑣 = Eq.17b
𝟎. 𝟔𝟔𝟕 ∗ 𝑭𝑷 𝑵𝟕 𝑷𝟏 𝐅ϒ 𝒙𝑻
𝑾 𝑻𝟏 Eq.18
𝐶𝑣 = N27: 0.67
𝑵𝟐𝟕 𝑭𝑹 𝚫𝐏(𝑷𝟏 + 𝑷𝟐 )𝑴
N22: 1500
𝑸 𝑴𝑻𝟏
𝐶𝑣 = Eq.19
𝑵𝟐𝟐 𝑭𝑹 𝚫𝐏(𝑷𝟏 + 𝑷𝟐 )
Specify the Variables Required to Size the Valve [5] a) Solve for the Required Flow Coefficient (Cv) (valve diameter known)
• Process fluid (e.g. oil and water); Flow coefficient value is different from volumetric flow rate than
• Service condi ons and fluid physical proper es such as volume rate mass flow rate units. When the gas-specific gravity, Gg, is known for
of flow (q), mass rate of flow (w), Upstream and downstream volumetric flow rate units,
pressure values (P1 and P2), pressure drop across 1. the flow coefficient is calculated as per Eq. (15 & 16).
the valve (ΔP = P1 − P2), upstream temperature (T1), gas-specific 2. If the molecular weight, M, is known for volumetric flow rate
gravity that is defined as a dimensionless number equal to the ra o units, the flow coefficient is calculated as per Eq. (17a).
of the flowing gas density to the density of air with both at standard 3. If the specific weight at the valve inlet, γ, is known for mass flow
condi ons (Gg), molecular weight (M), the ra o of specific heats (k), rate units, the flow coefficient is calculated as per Eq. (17b).
gas compressibility factor (Z), viscosity (µ) and solid% . 4. If the molecular weight, M, is known for mass flow rate units, the
• Constants FL, Fd and xT: A er selec ng required Valve type, trim flow coefficient is calculated as per Eq. (16).
type and flow direc on, these valves can be es mated using Table 2. 1. Select xT from Table 2
2. Calculate Fγ from Eq. 34
3. Calculate FγxT and determine if flow is choked or not. (x <
Determine Piping Geometry Factor (FP) FγxT)
FP is a correc on factor that accounts for the pressure losses due to 4. If flow is choked Y = 0.667 if not Calculate Y.
piping fi ngs that are used to change the direc on or size of the 5. Calculate Ci by Eq. 6,7,8 (non-choked) or Eq. 12,13,14
piping or take branches from the headers such as tees, elbows, and (Choked)
reducers. If fi ngs are used before or a er the valve, the effect of 6. Calculate Rev using Eq. 28, Obtain Fd from Table 2.
piping geometry must be considered for the valve sizing. However, if 7. If Rev > 10000, flow is turbulent and If valve size = pipe size
no fi ng is a ached to the valve, FP is equal to one and does not Use Ci previously Calculated. If not, Calculate Fp and
impact the valve sizing process. Equa on (36) shows how to calculate subsequently its Cv using the iteration method described in
the piping geometry. flowchart.
8. If Rev < 10000, flow is non-turbulent. Calculate FR and its
Determine the Expansion Factor (Y) subsequently its Cv using the iteration method described in
When sizing components such as valves in the piping system handling flowchart.
compressible gas, it is essen al to account for all factors, including
the expansion factor that would affect the flow rate and pressure Select the appropriate flow coefficient table and the calculated flow
drop in the piping system. As compressible gas flows through each coefficient valve to select the suitable valve size.
component in the piping system, the head loss causes expansion as Use fisher tables (or any representa ve tables) to select a valve with
well as change in density, temperature, and velocity. The gas Cv larger than your calculated Cv. Select the desired flow
expansion process adds resistance to the flow rate, resul ng in a flow Characteris c and calculate the es mated opening for the chosen
rate reduc on for a given pressure drop. Thus, the gas expansion valve.
factor, Y, shall be taken into account.
𝑥 b) An alterna ve op on for valve sizing (valve diameter isn’t known):
𝑌 = 1−
3Fϒ 𝑥 Use the simplified equa ons described earlier then use the sizing
x: The pressure drop ra o= ΔP/P1; If x < 0.2, the change in gas density equa on in the following sec on (Another method for es ma ng
and expansion is small, the fluid can be assumed incompressible, and control valve size (if tables not available))
the parameter Y can be neglected. For x > 0.4, the expansion factor Check Noise levels, velocity and Mach Number
shall be considered as the fluid can be assumed compressible.
xT: The pressure drop ra o factor for valves installed without
a ached fi ngs can be extracted from Table 2. More precisely, xT is
the pressure drop ra o required to produce cri cal or maximum flow
through the valve when Fγ = 1. It should be noted that the value of Y
cannot be less than 0.667. When the calculated Y value is less than
0.667, the Y value shall be assumed to be 0.667.
Suppose the control valve is installed with fi ngs such as reducers or
elbows a ached to it. In that case, the effect of fi ngs shall be taken
into account by replacing XT by a factor of XTP that can be calculated
from Eq. (33) as follows.
𝑥
𝐹
𝑥 = Eq.33
𝑥 ∑𝑘 𝐶
1+ 𝑁 ( )
𝑑
N5: 0.00241, d in mm
Determina on of correc on factors[6]
𝐾 =1− Eq. 22
N1 : 0.865
N2: 0.00214
Q in m3/hr
P in bar
∑ / Density in kg/m3
𝐹 = 1+ ( ) Eq. 36
ρ in kg/m3
ρ0 : 1000 kg/m3
D and D in mm
If Rev > 10,000 and Valve Size doesn’t equal pipe size:
Assume Fp =1 and Cv = Ci (your calculated Cv previously) and re-calculate Fp, check whether
flow is choked or not and re-calculate Cv accordingly, if Ci/C > 0.99 you can use your
calculated Cv. If not, re-assume Ci to be your calculated Cv and Fp to be your re-calculated Fp.
Reynolds number factor FR 1. For Transi onal flow
The Reynolds number factor FR is required when .
0.33𝐹 𝑅𝑒
non-turbulent flow condi ons are established 𝐹 =1+( ) log ( )
through a control valve because of a low 𝑛 . 10000
pressure differen al, a high viscosity, a very Where:
small flow coefficient, or a combina on thereof.
The FR factor is determined by dividing the flow 𝑁
𝑛 =
rate when non-turbulent flow condi ons exist by 𝐶
( )
the flow rate measured in the same installa on 𝑑
under turbulent condi ons. 2. For laminar flow
Tests show that FR can be determined from the
.
curves given in Figure 3a or 3b using a valve 𝐹 = 𝑛 𝑅𝑒 (not to exceed FR =1)
Reynolds number calculated from the following
equa on: For reduced trim valves where Ci /d2 at rated
travel is less than 0.016 N18 and Rev ≥ 10,
𝑁𝐹 𝑄 𝐹 𝐶 calculate F from the following equa ons:
𝑅𝑒 = ( + 1) .
Eq.28
𝑣 𝐶𝐹 𝑁 𝐷 1. For Transi onal flow
.
0.33𝐹 𝑅𝑒
This calcula on will require itera on. Proceed by 𝐹 =1+( ) log ( )
𝑛 . 10000
calcula ng the flow coefficient C for turbulent
flow. The valve style modifier Fd converts the n2= 1 + N32 ( ) /
Next, establish Ci as per Equa on 26. Liquid pressure recovery factor without
Ci = 1.3C (Equa on 26) a ached fi ngs FL
Apply Ci as per Equa on 26 and determine FR. FR FL is the liquid pressure recovery factor of the
is determined from Figure 3a for full-size trim valve without a ached fi ngs. This factor
valves. accounts for the influence of the valve internal
FR is determined from Figure 3b for reduced trim geometry on the valve capacity at choked flow. It
valves where Ci /d 2 at rated travel is less than is defined as the ra o of the actual maximum
0.016 N18. flow rate under choked flow condi ons to a
theore cal, non-choked flow rate which would
C/FR ≤ Ci Eq. 29 be calculated if the pressure differen al used
If the condi on of Equa on 29 is sa sfied, then was the difference between the valve inlet
use the Ci established from Equa on 26. If the pressure and the apparent vena contracta
condi on of Equa on 29 is not met, then repeat pressure at choked flow condi ons. The factor FL
the above equa on by again increasing Ci by 30 may be determined from tests in accordance
percent. This may require several itera ons un l with ANSI/ISA-75.02-1996. Typical values of FL
the condi ons required in Equa on 29 are met. versus percent of rated flow coefficient are
The equa ons defining non-turbulent flow for shown in Figure 4a/b and Table 2.
full size and reduced trim valves are stated in
Combined liquid pressure recovery factor
Annex G in ANSI/ISA 75.01.01. (Summarized
below, Also check flowchart below) and piping geometry factor with a ached
fi ngs FLP
For full size trim where Ci /d2 ≥ 0.016 N18
FLP is the combined liquid pressure recovery
and Rev ≥ 10, calculate FR from the following
factor and piping geometry factor for a control
equa ons:
valve with a ached fi ngs. It is obtained in the d) Reynolds number; and
same manner as FL. e) specific heat ra o .
To meet a devia on of ±5% for FLP, FLP shall be The influence of items a), b), c), and e) is
determined by tes ng. When es mated values accounted for by the pressure differen al ra o
are permissible, the following equa on shall be factor xT.
used: The Reynolds number is the ra o of iner al to
𝐹 viscous forces at the control valve orifice. In the
𝐹 = case of compressible flow, its value is generally
𝐹 𝐶 Eq.30 beyond the range of influence, except where the
1+ 𝐾( ) flow rate or the Cv is very low or a combina on
𝑁 𝑑
of both exist .
The pressure differen al ra o xT is influenced by
the specific heat ra o of the fluid.
Liquid cri cal pressure ra o factor FF
FF is the liquid cri cal pressure ra o factor. This 𝑥
factor is the ra o of the apparent vena contracta
𝑌 =1−
3F 𝑥 Eq.32
pressure at choked flow condi ons to the vapor
pressure of the liquid at inlet temperature. At
vapor pressures near zero, this factor is 0.96. The value of x for calcula on purposes shall not
Values of FF may be determined from the curve exceed Fϒ.xT. If x > Fϒ.xT, then the flow becomes
in Figure 5 or approximated from the following choked and Y = 0.667.
equa on:
Pressure differen al ra o factor without fi ngs
xT
𝑃
𝐹 = 0.96 − 0.28 Eq.31 xT is the pressure differen al ra o factor of a
𝑃
control valve installed without reducers or other
fi ngs. If the inlet pressure P1 is held constant
and the outlet pressure P2 is progressively
Specific heat ra o factor Fϒ lowered, the mass flow rate through a valve will
The factor xT is based on air near atmospheric increase to a maximum limit, a condi on
pressure as the flowing fluid with a specific heat referred to as choked flow. Further reduc ons in
ra o of 1.40. If the specific heat ra o for the P2 will produce no further increase in flow rate.
flowing fluid is not 1.40, the factor Fis used to This limit is reached when the pressure
adjust xT. Use the following equa on to calculate differen al x reaches a value of Fϒ.xT. The
the specific heat ra o factor: limi ng value of x is defined as the cri cal
𝛾 differen al pressure ra o. The value of x used in
Fϒ = Eq.34 any of the sizing equa ons and in the
1.4
rela onship for Y (Equa on 32) shall be held to
this limit even though the actual pressure
Expansion factor Y differen al ra o is greater. Thus, the numerical
value of Y may range from 0.667, when x = Fϒ.xT,
The expansion factor Y accounts for the change
to 1.0 for very low differen al pressures.
in density as the fluid passes from the valve inlet
to the vena contracta (the loca on just
Representa ve values of xT for several types of
downstream of the orifice where the jet stream
control valves with full size trim and at full rated
area is a minimum). It also accounts for the openings are given in Table 2. Cau on should be
change in the vena contracta area as the exercised in the use of this informa on. When
pressure differen al is varied. precise values are required, they should be obtained
Theore cally, Y is affected by all of the following: by test.
a) ra o of port area to body inlet area;
b) shape of the flow path;
c) pressure differen al ra o x ;
Pressure differen al ra o factor with
a ached fi ngs xTP
N4 : 0.076
N32 : 127
N5 : 0.00241
D or d in mm
Q in m3/hr
Kinema c viscosity
(𝑣) in m2/s
Control Valve Sizing flow charts
Li ra o of the control valve is preferred to be between 50% - 70% for normal flow condi ons.
(50% ≤ Li % ≤ 70%)
(1) The flow characteristic of the control valve will be “Linear”, except for the following
items (2) or (3).
Also ‘Linear’ could be selected when the change in valve opening causes a maximum of
10% change in ΔP
(2) Equal percentage will be applied for the following services:
- For temperature control
- For the services where the pressure drop across control valve or flow rate varies
significantly
- For service with relatively small ϒo
(3) Quick opening (denoted on-off) will be applied for the following services:
- Self-actuated pressure regulator.
- On-off control service.
- For service with relatively large friction, such as product rundown line.
(1) General features and typical application for control valves are summarized as follows:
Globe body valve Flow characteristic: Optional Most common, provided in the middle of
Rangeability (effective): 10 straight line.
Angle body valve Flow characteristic: Optional Application to flashing, erosive, high AP, or
Rangeability (effective): 10 slurry services
Available single seat only
Solid contained less than 0.05 kg/cm3
Butterfly valve Flow characteristic: Optional Applicable to high flow rate and low pressure
Rangeability (effective): 15 drop services
Sampled and economical, use for 4” and larger
Gas service with low pressure drop
High viscous, slurry service
Max. opening is usually limited to 60 degrees
for throttling
Don’t use with opening of 10% or less as a
control valve.
Ball valve Flow characteristic: Essentially EQ% Applicable to high flow rate and high shut-off
Rangeability (effective): 50 pressure services
Low resistance at full-open
Suitable as a shut-off valve
Solid contained service
High rangeability
Eccentric plug- Flow characteristic: Globe/Ball Advantages over butterfly valves
valve (Camflex) Rangeability (effective): 100
Calcula ons
𝑃 70.1
𝐹 = 0.96 − 0.28 = 0.96 − 0.28 = 0.944
𝑃 22120
∆Pmax = 𝐹𝐿 (𝑃 − 𝐹 𝑃 ) = 497.2 𝑘𝑃𝑎
Which is more than the differen al pressure (∆P = 460 kPa); therefore, the flow is non-choked,
and the flow coefficient C is calculated using Equa on 1.
𝑄 𝜌 ⁄𝜌
𝐶𝑣 = = 165
𝑁 ∆𝑃
Where Q = 360 m3 / h; N1 = 1 × 10–1 from Table 1; ρ1/ρo = 0.965; and ∆P = 460 kPa. Cv is in Kv’s
units.
Calculate Rev
𝑁𝐹 𝑄 𝐹 𝐶
𝑅𝑒 = ( + 1) . = 2.967 ∗ 10
𝑣 𝐶𝐹 𝑁 𝐷
N2 = 1.60 × 10–3 from Table 1; N4 = 7.07 × 10–2 from Table 1;
Fd = 0.46; Q = 360 m3 /h; ν = 3.26 × 10–7 m 2 /s; Ci = C = Kv = 165 m3/ h;
FL = 0.90; and D = 150 mm.
Since the valve Reynolds number is greater than 10,000, the flow is turbulent, and the flow
coefficient C as calculated above is correct.
Example no. 2 : Compressible flow – non-choked flow with a ached fi ngs
Process data:
Fluid: carbon dioxide
Inlet temperature: T1 = 433 K
Molecular mass: M = 44.01 kg/kmol
Kinema c viscosity: ν = 1.743 × 10–5 m 2 /s
Specific heat ra o: γ = 1.30
Compressibility factor: Z = 0.988
Inlet absolute pressure: P1 = 680 kPa
Outlet absolute pressure: P2 = 310 kPa
Flow rate: Q = 3,800 standard m3 /h at 101.325 kPa and 0°C
Inlet pipe size: D1 = 80 mm
Outlet pipe size: D2 = 100 mm
Reducers: short length, concentric
Valve data:
Valve style: rotary
Trim: eccentric rotary plug
Flow direc on: flow-to-open Valve size: d = 50 mm
Pressure differen al ra o factor: xT = 0.60 (from Table 2)
Liquid pressure recovery factor: FL = 0.85 (from Table 2)
Valve style modifier: Fd = 0.42 (from Table 2)
Calcula ons:
ϒ 1.3
Fϒ = = = 0.929
1.4 1.4
𝛥𝑃 680 − 310
𝑥= = = 0.544
𝑃 680
Which is less than Fϒ.xT = 0.557; therefore, the flow is non-choked and the flow coefficient is
calculated from Equa on 11. Next, Y is calculated from Equa on 32;
.
𝑌 =1− =1− = 0.674
ϒ ∗ . ∗ .
. ∗ ∗ .
𝐶𝑣 = = = 62.7 𝑓𝑜𝑟 𝑘
. ∗ ∗ ∗ . .
Calculate Rev
𝑁𝐹 𝑄 𝐹 𝐶 0.0707 ∗ 0.42 ∗ 3800 0.85 62.7
𝑅𝑒 = ( + 1) . = ( + 1) . = 8.96 ∗ 10
𝑣 𝐶𝐹 𝑁 𝐷 1.743 ∗ 10 √62.7 ∗ 0.85 80
Since the valve Reynolds number is greater than 10,000, the flow is turbulent. Now, calculate the
effect of the inlet and outlet reducers on C.
Since both reducers are concentric, short length, the velocity head loss coefficients can be
calculated as follows:
𝑑
𝐾 = 0.5(1 − ) = 0.186
𝐷
𝐾 =1− = 0.847
𝐾 =1− = 0.938
∑ 𝐾 = 𝐾 + 𝐾 + 𝐾 − 𝐾 =0.658
Now, the effect of the reducers is calculated by itera on, star ng with Ci = C and FP(1) = 1.
/ /
∑𝐾 𝐶 0.658 62.7
𝐹 = 1+ ( ) = 1+ ( ) = 0.891
𝑁 𝑑 1.6 ∗ 10 50
Since FP(2) /FP(1) = 0.891/1 < 0.99, one more itera on step shall be done
𝐶 62.7 m
𝐶 = = = 70.4
𝐹 ( ) 0.891 h
/ /
∑𝐾 𝐶 0.658 70.4
𝐹 ( ) = 1+ ( ) = 1+ ( ) = 0.868
𝑁 𝑑 1.6 ∗ 10 50
Now, FP(3)/FP(2) = 0.868/0.891 > 0.99 so FP(3) will be used as FP for the final calcula on.
𝑥 0.6
𝐹 0.868
𝑥 = = = 0.626
𝑥 𝐶 0.6 70.4
1+ ( ) 1+ ( )
𝑁 𝑑 1.8 ∗ 10 50
and with this Fϒ xTP = 0.582, which is greater than x = 0.544.
Input steps
Recommended Reads
This tool was developed for process engineers to estimate a vessel Size on site quickly and compare
various methods from the literature. As a part of a larger project to develop what is similar to Carl
Branan’s book “process engineers’ Pocket Handbook” these tools would allow a process engineer to
quickly calculate/estimate equipment efficiencies or sizing using standardized calculations.
The aim here is to take little-known data from the field (flow, pressures, temperatures,
compositions..etc.) and use it as input for a rough estimation without having to return to the office to
use commercial software or calculations Excel sheets to validate or to calculate. Additionally, these tools
may also serve as a gathered data validation tool.
Nomenclature*
Property Table
Gas Mass Fraction -
Gas Mass Flow rate kg/hr
Liquid Mass Flow rate kg/hr
Gas Flow rate Nm³/hr
Liquid Flow rate m³/hr
Stream Pressure kg/cm²
Stream Temperature C
Gas Density kg/m³
Liquid Density kg/m³
Gas Viscosity cP
Liquid Viscosity cP
Gas Molecular weight -
Gas Compressibility -
2.7𝛾𝑃(𝑃𝑠𝑖𝑎) 𝑙𝑏
𝜌𝑔 = ( 3)
𝑇(ᵒ𝑅)𝑍 𝑓𝑡
𝑙𝑏
𝜌𝐿 = 62.4 ∗ 𝑆. 𝐺( )
𝑓𝑡 3
6. Calculate Actual Volume flow rate of gas
TOTAL • Try to avoid vessels with wall thickness
Horizontal or vertical design (notes) [2] greater than 100 mm as these require
• Provided a sufficient L/D ratio is selected, a special fabrication and can prove expensive.
horizontal separator is more efficient than a Calculation theory and equations
vertical for the same flow area.
• Vapor velocity in a horizontal drum can • Valid only for pure gravity settlers with no
exceed the liquid settling velocity provided internals to enhance separation)
L/D > 1. For vertical drums, the velocity
Liquid-vapour settling velocity
cannot.
• Horizontal drums are more effective and 1. Vs = K* ((ρL- ρv)/ (ρv))1/2 m/s
0.5
geometrically more practical for heavy 2. K = 0.003615 (D/C) m/s
liquid phase removal than vertical drums. 3. UT = 0.003615 (D/C) * ((ρL- ρv)/ (ρv))1/2
0.5
• A rising liquid level in a vertical drum does K in m/s, Densities in kg/m3, D is droplet size in
not alter the vapor flow area. Consequently, microns and C is drag coefficient.
vertical drums are preferred for compressor • For medium and low pressure’ with gases of
and fuel gas KO drums. Vertical drums viscosity less than 0.01 cP Figure 2 can be
utilize a smaller plot area and are easier to used to estimate Vs.
instrument with alarms and shutdown • For higher pressures (> 50 bar) or viscosities
controls. For floating installations are over 0.01 cP it is necessary to calculate Vs.
preferred as less “sloshing” occurs. The drag coefficient is calculated using
• For high volume flow rates a split flow Figure 3 (curve 2) where:
horizontal drum is preferable as smaller CRe = (0.3072*10-11*ρv*D3*(ρL - ρv) / µ2
drum diameters can be used. The preferred
split flow arrangement is a single centre Equation in step 3 is then used to calculate Vs.
entry nozzle with two end exists. Head exits Liquid-liquid settling velocity (based on Stoke's
can be used where plot space is limited. law of terminal settling)
• Each design case must be evaluated
The following equation can be used for
separately but in general the following can
calculating the settling velocity of water in oil or
be used as a guideline:
the upward “settling” of oil in water. The
Vertical drums important fact is to use the viscosity of the
Compressor KO drums continuous phase i.e : for oil settling upwards
De-gassing boots through water use the water viscosity, for water
Fuel gas KO drums settling in oil use the oil viscosity.
Absorber feed KO drums Ut = (g*D2*(ρH – ρL)) / 18*µc m/s
Floating installations
Horizontal drums Ut, = terminal velocity m/s
Production separators HP g: gravitation acceleration m/s2
Reflux drums ρH = density heavy fluid kg/m3
3-phase separation ρL = density light fluid kg/m3
Flare KO drums µc = viscosity (continuous) kg/m.s
Setting the particle size to 125 microns and eliminator. It is recommended to install
using more useful units gives: a mist eliminator for most
applications. if not calculate Vs using
Ut = 0.513*(ρH – ρL) / µc
Ut mm/min µc in cP the equation in step 3 earlier.
3. De-rate the calculated settling velocity
• The above equation is valid for Reynold’s no by 85 % design margin to give a
of 0.1 – 0.3 If the calculated settling velocity maximum allowable vapor velocity.
is > 250 mm/min use 250 max 4. Calculate the drum internal diameter
and round to the nearest 50 mm. (note
Vessel volumes
further adjustment of ID: OD can be
Partial volumes of a horizontal cylinder can be made to suit standard head
calculated using the partial volume charts in dimensions).
Figure 2 or estimated using the following 5. Check if wall thickness is less than 100
equations : (for vessels with a diameter < 1.2 m mm.
ignore head volumes)
𝐷2 𝐷−2ℎ 𝐷
𝐴𝐿 = cos−1 [ ] − ( − ℎ)(𝐷ℎ − ℎ2 )0.5 m2
4 𝐷 2
Horizontal cylinder H1
VC = AL* L m in radians H2
2 dished head H3
Vdh = 0.21543 h2 (1.5 D – h) m3
ɸ
2 elliptical heads H4
Vel = 0.52 194 h² (1.5 D – h) m³ (most common)
2 hemispherical heads
HLA
Vhh = 1.047 h2 (1.5 D -h) m³ (gives extra vol) H5 HLL
Volume up to baffle
H6 NLL
For depth h
(0.52194 h2 (1.5D-h))/2+ AL*B (Elliptical heads) H7
LLL
Calculation procedure horizontal vessel (Vapor- • For high volumetric flows of gas with small
liquid) liquid volumes consider using a split flow
L’ arrangement. Design is as above but with
half vapor volume flow.
L • Normal design is with top entry and exit
nozzles. However, if space is limiting
(primarily offshore) head mounted nozzles
can be used to increase flow path.
ɸ HLL • L is designated as the flow path length i.e
NLL distance between inlet and outlet nozzle. L’
LLL or TTL is the tangent-tangent length. For 1st
estimates L’ = L + 1.5 ɸi + 1.5 ɸ2
2 conditions
Sizing Procedure
Find V* qd and qf
Fd = F g
Vf.Af = qf qs. qs =
m = qs.pf
In summary, the size (diameter and seam-to-
seam length or height) of a vertical separator is
determined as follows:
Find V*
Vf.Af = qf Vmax = Vf
qs = qs.pf/ps 1. Determine the minimum allowable vessel
c(po – pg/pi)**0.5 diameter.
0.5
𝑇𝑍 (ρo -ρg ) dm
D2 = 5.508 ∗ 𝑄𝑔 ∗ ∗ [ . ] in2
𝑃 ρg Cd
Find V* qs = qf By ideal gas law change
(ps/ps) 2. For diameters larger than the minimum, the
Vf.Af = qf (pf/ps)(Ts/Tf)z As f(T,P)
following is used to determine combinations
of D and H.
D2H=8.565*Qo*t in3
Cd could only be determined by an iterative 3. The seam-to-seam length, Ls, for each
procedure as follows: combination of D and H is determined using
1. Assume a value for Cd (a value of 0.34 one of the following expressions as
could be used as a first assumption). appropriate:
2. Calculate the velocity, u, from: If D < 36 in, then Ls =(H+76)/12 ft
0.5 If D> 36 in, then Ls =(H+D+40)/12 ft
(ρo -ρg ) dm
u=0.01186. [ . ] 4. For each combination of D and Ls, the
ρg Cd
slenderness ratio, SR, defined as the ratio of
3. Calculate Re from
length to diameter, is determined.
Re = 0.0049 ρg dm u/µg
Separators with SR between 3 and 4 are
4. Calculate Cd and compare to the assumed
commonly selected.
value from:
Cd = 0.34 + 3/Re0.5 +24/Re Notes: dm in micron, densities in lb/ft3, Qg in
5. If no match is obtained, use the calculated MMSCFD and Qo in bbl/d, t in mins, H and D in
value of Cd and repeat steps 2–4 until inches, T in oR and P in Psia.
convergence is obtained.
Sizing Horizontal Gas–Oil Separators volumetric flow rate, Qg, by one-half of the
separator cross-sectional area, A, that is,
ug = Qg /(0.5[0.25*π*D2]) ft/s
Qg is usually reported in units of MMSCFD
and should, therefore, be converted into
actual cubic feet per second (ft3 /s); also D,
which is usually given in inches, should be
converted into feet in order to obtain the
velocity in units of feet per second (ft/s). The
preceding equation, therefore, becomes
ug = 120*Qg*TZ/(D2*P) ft/s
The gas travels horizontally along the
effective length of the separator, L (ft), in a
The horizontal separator is determined by its time, tg, that is given by:
required capacity for gas and oil. It has been tg = L/ug s
shown that the gas capacity constraint for This time must, at least, be equal to the time
vertical separators determines the minimum it takes the smallest oil droplet to be
allowable vessel diameter. For horizontal removed from the gas to travel a distance of
separators, however, the gas capacity constraint D/2 to reach the gas–oil interface. This
yields, as shown in the following section, a
settling time, ts, is obtained by dividing the
relationship between the diameter and effective
distance (D/2) by the settling velocity from u
length of the separator. This along with a similar
equation; therefore,
relationship derived from the liquid capacity
0.5
constraint are used in determining the size of 𝐷 (ρo -ρg ) dm
𝑡𝑠 = . 0.01186. [ . ]
the separator. In reality, either the gas capacity 2𝑥12 ρg Cd
constraint or the liquid capacity constraint Equating both above equation, substituting for
governs the design and only one of the two ug from ug = 120*Qg*TZ/(D2*P), and solving for
constraints equations is used in determining the
the product LD, we obtain
size. 0.5
𝑄𝑔 𝑇𝑍 (ρo -ρg ) dm
𝐿𝐷 = 422. 𝑃
.[ ρg
.C ] ft /in
d
In the following discussion, it is assumed that
This equation provides a relationship between
each of the gas and oil phases occupies 50% of
the effective separator volume. Similar the vessel diameter and effective length that
equations as those derived in the following satisfies the gas capacity constraint. Any
could be obtained for other situations where combination of D and L satisfying the above
either of the two phases occupies more or less equation ensures that all oil droplets having
than 50% of the separator effective volume. diameter dm and larger will settle out of the gas
flowing at a rate of Qg MMSCFD into the
When Sizing a separator, two constraints should
separator that is operating at P psia and T°R.
be accounted for:
3. Liquid Capacity Constraint: A gas–oil separator
2. Gas Capacity Constraint: Since the gas has to have a sufficient volume to retain the
occupies the top half of the separator, its liquid for the specified retention time before it
average flowing velocity within the leaves the separator. For a horizontal separator
separator, ug, is obtained by dividing the
that is half full of liquid, the volume occupied by The preceding equations and sizing procedures
the liquid is given by are very sufficient for the determination of
Vo = 0.5*0.25* π*(D/12)2.L ft3 separator diameter and length as well as for the
performance evaluation of existing separators.
Substituting in Qo (bbl/day) =257* Vo / t, the Students and practicing engineers should be
following equation is obtained: familiar with such design equations and
D2L = 1.428Qot ft3 procedures before attempting the use of
commercially available software. The retention
The equation provides another relationship time is an important parameter in designing
between D and L that satisfies the liquid gas–oil separators. It is best obtained from
capacity (retention) time constraint. laboratory tests that simulate the field
Sizing Procedure operating conditions. This, however, may not
always be available. In such cases, experience
For a given set of operating conditions and data from offset fields, if available, will be
(pressure, temperature, gas and oil flow rates, very valuable.
gas and oil properties, and oil retention time),
the size (diameter and seam-to-seam length) of
a horizontal separator is determined as follows:
Ls = Lg + D/12 ft
2 6686.𝑄𝑜 .𝜇𝑜
𝐷𝑚𝑖𝑛 = (𝛾 2 in2
0 −𝛾𝑊 ).𝑑𝑚
0.5
2 𝑇𝑍 (ρo -ρg ) dm
𝐷𝑚𝑖𝑛 = 5508 ∗ 𝑄𝑔 ∗ ∗ [ . ]
𝑃 ρg Cd
Where Δγ = γw – γo. γo and γw are the specific 2. Determine the minimum diameter that
gravity of oil and water, respectively. satisfies the gas capacity constraint from
0.5
2 𝑇𝑍 (ρo -ρg ) dm
2. Water Droplets Settling Constraint: The 𝐷𝑚𝑖𝑛 = 5508 ∗ 𝑄𝑔 ∗ ∗ [ . ]
𝑃 ρg Cd
condition for the settling and separation of
3. The larger of the two minimum diameters
water droplets from the oil is established by
determined in steps 1 and 2 is then
equating the average upward velocity of the
considered as the minimum allowable vessel
oil phase, uo, to the downward settling
diameter.
velocity of the water droplets of a given size,
uw. The average velocity of the oil is obtained
4. For various values of diameter larger than can, therefore, be summarized in the following
the minimum allowable vessel diameter, steps:
Determine combinations of diameters and
1. Determine the value of Aw/A from:
liquid heights. 𝐴𝑤 𝑄𝑤 𝑡𝑤
= 0.5
(H𝑜 + H𝑤 ). 𝐷 2 = 8.576 ∗ (𝑄𝑜 𝑡0 + 𝑄𝑤 𝑡𝑤 ) in3 𝐴 𝑄𝑜 𝑡𝑜 +𝑄𝑤 𝑡𝑤
5. For each combination, determine the seam- 2. Use Figure 1 to determine the value of Ho/D
to-seam length from the following: for the calculated value of Aw/A.
If D < 36 in,
Figure 1: Ho/D ratio as a function Aw/Atotal
then Ls =(1/12)*(Ho+ Hw+76) ft
If D > 36 in,
then Ls =(1/12)*(Ho+ Hw+D+40) ft
Notes:
T in absolute R
P in Psia
Densities in lb/ft3
Specific gravities are dimensionless
Qo and Qw in bbl/day
Qg in MMSCFD
Droplet sizes dm,w and dm,o in microns
Diameter (D) in inches
LD in in.ft and D2L in in2.ft
Ho in inches
Areas in ft2
L in ft
VH = TH * QL ft3
24 in x 5 ft
24 in x 7½ ft
24 in x 10 ft
30 in x 5 ft
30 in x 7½ ft
30 in x 10 ft
36 in x 5 ft
36 in x 7½ ft
36 in x 10 ft
36 in x 15 ft
42 in x 7½ ft
42 in x 10 ft
42 in x 15 ft
Figure 7: G values as a function of density difference and the downcomer
allowable flow
48 in x 7½ ft
48 in x 10 ft
48 in x 15 ft
54 in x 7½ ft
54 in x 10 ft
54 in x 15 ft
60 in x 7½ ft
60 in x 10 ft
60 in x 15 ft
Example No. 1: Two-Phase Vertical Separator[1]
Step 2: Determine Cd
assumed drag coefficient of Cd = 1.21
0.5
𝜌𝐿 − 𝜌𝑔 𝑑𝑚
𝑢 = 0.01186 [( ) ] = 0.41 𝑓𝑡/𝑠
𝜌𝑔 𝐶𝑑
𝜌𝑔 𝑑𝑚 𝑢
𝑅𝑒 = 0.0049 = 52
𝜇𝑔
3 24
𝐶𝑑 = 0.34 + + = 1.22, Acceptable Cd assumption
𝑅𝑒 0.5 𝑅𝑒
Step 4: Check for Gas Capacity Constraint
𝑇𝑍 𝜌𝑔 𝐶𝑑 0.5
𝐷𝑚𝑖𝑛 2 = 5058𝑄𝑔 ( ) ( ) = 2059
𝑃 𝜌𝐿 − 𝜌𝑔 𝑑𝑚
∴ 𝐷𝑚𝑖𝑛 = 45.4 𝑖𝑛
Step 5: Check for liquid capacity constraint
D2H = 8.565*Qot = 25695 in2
The solution selected has the highest SR ratio in the range of 2-4 as per GPSA
TOTAL method
Liquid-vapour settling velocity and Dmin
Step 3: Check Operating conditions and physical props output Step 4: Check TOTAL output image and table
Step 5: Check Hussien’s output table Step 6: Check Svrcek output image and table
1. Calculate The vapor volumetric flow rate Assumptions for the Surge and Holdup
145600 𝑙𝑏/ℎ𝑟 times (Tables 6 and 7):
𝑄𝑣 = 𝑠 = 10.09 𝑓𝑡 3 /𝑠
3600 ∗ 4.01 𝑙𝑏/𝑓𝑡 3
ℎ Service: Unit Feed Drum
2. Calculate The liquid volumetric flow rate
Surge 10 mins, Holdup 5 mins
46100 𝑙𝑏/ℎ𝑟
𝑄𝐿 = = 19.79 𝑓𝑡 3 /𝑚𝑖𝑛
𝑚𝑖𝑛 3 Personal, Factor = 1
60 ∗ 38.83 𝑙𝑏/𝑓𝑡
ℎ
Well-Instrumented, Factor = 1
3. Calculate the vertical terminal velocity
For other examples the factors are
38.83 − 4.01 assumed both 1 and 1 too.
𝑈𝑡 = 0.13√ = 0.38 𝑓𝑡/𝑠
4.01
K = 0.13 (GPSA Value divided by two since “no mist eliminator is used”)
Uv = 0.75 Ut = 0.29 ft/s
𝑉𝐻 + 𝑉𝑠 197.9 + 98.95
𝐿= = = 29.3, 𝑠𝑎𝑦 29.5 𝑓𝑡
𝐴 𝑇 − 𝐴𝑣 − 𝐴𝐿𝐿𝐿 19.63 − 7.34 − 2.16
2 𝑓𝑡
∅= = 6.9 𝑠
0.29 𝑓𝑡/𝑠
10.09 𝑓𝑡 3 /𝑠
𝑈𝑉𝐴 = = 1.37 𝑓𝑡/𝑠
7.34 𝑓𝑡 2
1072 ∗ 60 1
𝑡𝑠 = + = 2.322 𝑖𝑛
2 ∗ 17500 ∗ 0.85 − 1.2 ∗ 1072 16
• tH = 2.25 in
1072 ∗ 60 1
𝑡𝐻 = + = 2.240 𝑖𝑛
2 ∗ 17500 ∗ 0.85 − 0.2 ∗ 1072 16
16.Calculate the surface area of the shell and heads according to table 9:
As = π*5 ft*29.5 ft = 463.38 ft2 and
AH = 1.09 * (5 ft)2 = 27.25 ft2
17.Calculate the approximate vessel weight = 50,224 lb
𝑙𝑏 2.375 𝑖𝑛
𝑊 = 490 (463.38 𝑓𝑡 2 + 2 ∗ 27.25)
𝑓𝑡 3 12 𝑖𝑛
𝑓𝑡
18.Try D = 5.5 ft and repeat calculations until minimum weight of shell and head is obtained
API Method
Drum Sizing
For first trial tres = 15 mins volume required =15 * QL =15*0.01*60= 9 m3
Trial 1 2
Selected h/D 0.7 assumed 0.6
modified
2 0.71
Vapor Area Av m Step 1 1.57
Step 6
% Total Area 70.5 Step 2 61.83
Step 5
2
Total Area Ar m 1.017 Step 3 2.54
Step 4
Liquid Area Al m2 0.305 Step 4 0.97
Step 7
Calculated drum φ mm 1130 Step 5
Selected Drum D mm 1200 Step 6 1800 ↗ Step 1
L/D (3 – 4 ) 4 assumed 4
Flowpath length L mm 4800 Step 7 7200 Step 2
T-T length L’ mm 5325 Step 8 7725 Step 3
HLL height mm 840 Step 9 1080 Step 8
3
Vol @ HLL m 4.8 Step 10 12.86 Step 9
LLL height mm 300 Step 11 300 Step 10
3
Vol @ LLL m 1.24 Step 12 2.22 Step 11
∆Vol m3 3.56 Step 13 10.64 Step 12
Calculated tres mins 5.93 Step 14 17.73 Step 13
Notes Small vessel Good enough!
tan-tan length L’ = L + 1.5φ1 + 1.5φ2 = 7200 + 1.5 * (150+200) = 7725 mm
Selected Diameter: 1650 mm x 7725 mm tan-tan length
Step 1: Select Desired Calculations and modify assumptions Step 2: Input operating conditions and physical props
Step 3: Check Operating conditions and physical props output Step 4: Check TOTAL output image and table
Step 5: Check Hussien’s output table Step 6: Check Svrcek output image and table
1. Calculate the vertical terminal velocity. Using table 5, calculate K using the york demister
equations
53.95 − 0.6973
𝑈𝑡 = 0.313√ = 2.74 𝑓𝑡/𝑠
0.6973
4∗165.32 𝑓𝑡 3 /𝑠
𝐷𝑉𝐷 = √ = 10.13 𝑓𝑡 , Use DVD = 10.5 ft
𝜋∗2.05 𝑓𝑡/𝑠
4. Calculate the settling velocity of the heavy liquid out of the light liquid. Using Table 12 ks = 0.163
Then:
UHL = 0.163*(62.11-53.95)/0.63 = 2.11 in/min
5. Calculate the settling velocity of the light liquid out of the heavy liquid.
ULH = 0.163*(62.11-53.95)/0.764 = 1.74 in/min
6. Calculate the light and heavy liquid volumetric flow rates
165000 𝑙𝑏/ℎ𝑟
𝑄𝐿𝐿 = = 5.10 𝑓𝑡 3 /𝑚𝑖𝑛
𝑚𝑖𝑛
60 ∗ 53.95 𝑙𝑏/𝑓𝑡 3
ℎ
1300 𝑙𝑏/ℎ𝑟
𝑄𝐻𝐿 = = 0.35 𝑓𝑡 3 /𝑚𝑖𝑛
𝑚𝑖𝑛 3
60 ∗ 62.11 𝑙𝑏/𝑓𝑡
ℎ
7. Assume HL =1 ft and calculate the time for heavy liquid to settle out of the light liquid phase:
tHL = 12*1.0/2.11 = 5.7 min
8. Assume HH = 1 ft and Calculate the time for the light liquid to rise out of the heavy liquid phase:
tLH = 12 * 1.0 /1.74 = 6.9 min
9. Calculate the baffle plate area
ρL- ρv = 53.95-0.6973 = 53.25 lb/ft3
Assume HR = 12 in, HL + HR = 24 in
Using Figure 3. G = 9800 gph/ft3
AD=7.48 gal/ft3*60 min/h*(5.10+0.35) (ft3)/9800 gph/ft3
Assume WD = 4 in
WD/D = 4/(12*10.5)=0.0317
Using table 3, AD/A = 0.0095
A = 0.25π*10.5 ft2=86.59 ft2
AD=0.0095*86.59=0.82 ft2
10. Calculate the residence time of each phase
ϴLL=1 ft * 85.77 ft2/5.10 ft3/min = 16.8 min
ϴHL=1 ft * 85.77 ft2/0.35 ft3/min = 247.4 min
11. Calculate the height of the light liquid above the outlet, based on holdup:
HR=5.10 ft3/min*25 min / 85.77 ft2= 1.5 ft
HS = (5.10+0.35) (ft3/min) (5 min)/86.59 ft2 = 0.31
Use HS=0.5 ft
12. Calculate dN according to table 5:
λ = QL/ (QL+ Qv) = (5.10 + 0.35)/(5.10+0.35+165.32*60)=0.0006
16500 1300
𝜌𝐿 = 𝑥53.95 + ∗ 62.11 = 54.55 𝑙𝑏/𝑓𝑡 3
17800 17800
ρM= ρL λ+ρv(1-λ)=54.55*0.0006+0.6973*(1-0.0006)=0.73
QM=165.32+(5.10+0.35)/60 = 165.41 ft3/s
4 𝑥 165.41 0.5
dN > ( ) = 1.73𝑓𝑡 𝑜𝑟 21 𝑖𝑛, 𝑈𝑠𝑒 24 𝑖𝑛
60𝜋/√0.73
Calculate HD = 24 + 24/2 = 36 in = 3.0 ft (min.)
Use HD=5.5 ft from figure 1, HT = 1 ft and s = 0.5 ft. Calculate HBN
HBN=0.5* 2 ft+ 2 ft= 3 ft
Set HA=0.5 ft
Final Dimensions
D = 10.5 ft, HH = 1.0 ft, HL = 1.0 ft, HR = 1.5 ft, HA=0.5 ft, HBN=3.0 ft, and HD=5.5 ft. Add 1.5 ft for
the mist eliminator.
HT = 14 ft
HT/D = 14/10.5 = 1.3
Add 2 ft to HT to be 16 ft(HR=2 ft, HD=7.0 ft) so that HT/D = 1.52 (HT/D should be in the range of 1.5 to 6.0)
Step 3: Check Operating conditions and physical props output Step 4: Check Hussien’s output tables
Step 5: Check Svrcek output image and table Step 6: Check Svrcek output image and table (Cont.)
40.5 − 0.19
𝑈𝑡 = 0.175√ = 2.55 𝑓𝑡/𝑠
0.19
6. Since the mass rate of vapor is about 82% of the loading, set Hv to be much greater than the
minimum. Assume Hv = 0.70D = (0.70)*(11.0 ft) = 7.70 ft. Using Table 3, Av/AT = 0.748, Av= 71.08
ft2
7. HLLL = 0.5 * 11.0 + 7 = 12.5 in, use 13 in.
HLLL/D = 13/(11*12) = 0.098, Using Table 3, ALLL/AT = 0.051
ALLL=0.051*95.03 ft2 = 4.85 ft2
8. Hw = 11.0 - 7.70 = 3.30 ft
9. Calculate L2
277.8
𝐿2 = = 14.54 𝑓𝑡, Use L2 = 15 ft
95.03−71.08−4.85
10. HHL = HLL= 3.30/2 = 1.65 ft
11. HHL/D= 1.65/11= 0.15, From Table 11, AHL/AT = 0.094
AHL=0.094*95.03 ft2 = 8.93 ft2
ALL=95.03 – 71.08 – 8.93 = 15.02 ft2
12. From Table 12, ks=0.333 as the light phase is the oil with SG < 0.85 and the heavy phase is water
UHL= (0.333)(62 - 40.5)/0.24= 29.83 in./min
Use 10 in./min (maximum)
ULH= (0.333)(62 - 40.5)/0.682 = 10.50 in./min, Use 10 in./min (maximum)
13. tHL = 12 in./ft)(1.65 ft)/10 in./min = 1.98 min, use 2.0 mins
tLH = tHL = 2.0 min
14. Calculate L1
𝑡𝐿𝐻 𝑄𝐻𝐿 2 𝑚𝑖𝑛𝑠 ∗ 2.02 𝑓𝑡 3 /𝑚𝑖𝑛
= = 0.45 𝑓𝑡
𝐴𝐻𝐿 8.93 𝑓𝑡 2
𝑡𝐻𝐿 𝑄𝐿𝐿 2 𝑚𝑖𝑛𝑠 ∗ 18.52 𝑓𝑡 3 /𝑚𝑖𝑛
= = 2.47 𝑓𝑡
𝐴𝐿𝐿 15.02 𝑓𝑡 2
Use L1 = 3 ft
15. L = 3.0 + 15.0 = 18.0 ft
16. Ф = 7.7 ft/1.91 ft/s = 4.03 s
17. UVA = Qv /Av = (343.57 ft3/s )/ (71.08 ft2) = 4.83 ft/s
18. Lmin = (4.83 ft/s) * (4.03 s) = 19.5 ft
19. Since L < Lmin set L = 19.5 ft (set L2 =16.0 ft, L1=3.5 ft)
20. L/D = 19.5/11.0 = 1.78
21. Assume dished heads per Table 9.
Calculate the thickness of the shell and heads according to table 9:
• Assume E = 0.85
• Assume SA 516 70 Carbon steel
Design Temp. = 650 oF
• S = 17500 psi
• CA = 1/16 in (Assumed)
• P = 25 + 30 = 55psig
• tS = 3/8 in
55 ∗ 132 1
𝑡𝑠 = + = 0.307 𝑖𝑛
2 ∗ 17500 ∗ 0.85 − 1.2 ∗ 55 16
• tH = 0.5 in
0.885 ∗ 55 ∗ 132 1
𝑡𝐻 = + = 0.495 𝑖𝑛
2 ∗ 17500 ∗ 0.85 − 0.1 ∗ 55 16
22. Calculate the surface area of the shell and heads according to table 7:
As = π*11 ft*19.5 ft = 673.87 ft2
and
AH = 0.842 * (11 ft)2 = 101.88 ft2
23.Calculate the approximate vessel weight = 17920 lb
𝑙𝑏 0.5 𝑖𝑛
𝑊 = 490 (673.87 𝑓𝑡 2 + 2 ∗ 101.88 𝑓𝑡 2 ) = 17920 𝑙𝑏𝑠
𝑓𝑡 3 12 𝑖𝑛
𝑓𝑡
24. In this example, calculations were performed for only one diameter. However, nearly the minimum
L/D corresponded to a diameter of 11.0 ft, therefore. the next diameter should be smaller, resulting
in a larger L/D. Also, calculations should be performed using a diameter of 11.5 ft.
25. For the light liquid compartment
HHLL = Hw = 3.3 ft ~ 3 ft, 4 in.
ANLL = 4.85 + 185.200/16 = 16.43 ft2
ANLL/AT = 16.43/95.03 = 0.173
Using Table 11, HNLL/D = 0.229
HNLL = 0.229*11 = 2.52 ft ~ 2 ft, 6 in.
HLLL = 13 in.
Comment: Due to the small amount of heavy liquid and a large amount of vapor, a better design
would have used a boot. A vertical vessel should be compared, As well.
Step 1: Select Desired Calculations and modify assumptions Step 2: Input operating conditions and physical props
Step 3: Check Operating conditions and physical props output Step 4: Check Svrceks’s output image and table
References
Resources
$/time
Fixed
Cost
The most economical velocity in piping costs
Pumping
continues to shift downwards power
cost
over the last 40 years Optimum
P
ipe sizing calculations include Carbon Stainless Carbon Stainless
Alumi- Brass
certain values that are time Terms Steel Steel Steel Steel
num 2008 2008
dependent, such as costs. While 1998 1998 2008 2008
some costs rise, others decrease n 1.35 0.7793 1.472 0.924 0.769 0.907
over time. However, there will always x 29.52 130 6.607 30.7 22.26 32.3
be a minimum total. For flow of flu- Le’ 2.74 2.74 2.74 2.74 2.74 2.74
ids calculations, heuristic criteria are M 0.102 0.102 0.064 0.064 0.064 0.064
often used, which is the main reason E 0.5 0.5 0.5 0.5 0.5 0.5
why these values must be updated P 150 150 150 150 150 150
from time to time. In this paper rec- K 0.04 0.04 0.07 0.07 0.07 0.07
ommended velocities are updated for Y 365 365 365 365 365 365
several materials of construction, and
& 0.55 0.55 0.55 0.55 0.55 0.55
other velocities for new materials are
Z 0.1 0.1 0.1 0.1 0.1 0.1
calculated, using the Generaux Equa-
F 6.7 7.5 6.5 7.4 7.1 7.2
tion [1].
a+b 0.2 0.2 0.2 0.2 0.2 0.2
Background a’ + b’ 0.4 0.4 0.4 0.4 0.4 0.4
Process design is one of the most im-
portant tasks for engineering. In the An example of this phenomena is costs of 1998 derived by the Marshall
development, optimal use of resources the determination of the recommended and Swift Index [2]. Similarly recom-
must be accomplished in order to velocity for a fluid in pipe sizing. The mended velocities are updated in this
reach equilibrium between minimum velocity is intrinsically related to fluid paper with costs for the current year.
cost and maximum efficiency. properties and a host of economic pa- We compare 1998 values with the ac-
Engineers tend to use heuristic for- rameters, including energy and pipe tual calculated ones (2008) and calcu-
mulas to make quick estimates. Such material costs (It’s important to re- late new values for other pipe materi-
calculated estimations are accurate member that pipe diameter is directly als such as brass and aluminum.
enough for grasping a quick overview proportional to its cost).
of the desired process. Since 1968 prices have changed Defining optimum
Some of these calculations involve due to inflation. In 1998 it was dem- For pipe design, several criteria have
certain factors that changing over the onstrated that the most economical to be taken into account. Common cri-
years. Nevertheless, the values fre- velocity in piping had shifted down- teria used are recommended velocity,
quently used as inputs in these equa- wards over the preceding 30 years [6]. economical criteria and other opera-
tions are over thirty years old, driving In that paper, recommended velocities tive criterion.
estimations even further away from were recalculated using the Gener- However, an economical criterion
reality [6]. aux equation [1] with parameters and is the decisive factor for determin-
2 CHEMICAL ENGINEERING WWW.CHE.COM JANUARY 2010
Recommended Fluid Velocities, 1998 to 2008 Costs
Stainless
velocity, ft/s
steel
43 Aluminium
Optimum equation
23 Brass Although there are several rules for
Carbon calculating the optimum pipe diame-
3 steel
0 0.2 0.4 0.6 0.8 1 1.2 ter, this paper only considers the Gen-
Density, lb/ft3 @ µ = 0.2 cp; D = 1 ft eraux equation [1], shown in the box
on p. TK. This equation determines
FIGURE 3. Recommended velocities for gases (a) and liquids (b) according to 1998 the fluid velovcities required to obtain
and 2008 prices. Recommended velocities tend to fall down as years pass by
the most economic pipe diameter.
Recommended Fluid Velocities at different Diameters for Carbon steel, 2008 Costs Ten years ago this equation was
8 used to obtain updated recommended
7 D = 6 in velocities because the 1968 values
were still used in that time [6]. For
Recommended
6 D = 8 in
velocity, ft/s
Circle XX on p. 62 or go to adlinks.che.com/23022-XX
CHEMICAL ENGINEERING WWW.CHE.COM JANUARY 2010 5
View publication stats