Process Engineering Calculations (Part 1)
Process Engineering Calculations (Part 1)
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