Processengineeringfundamentals Researchgate
Processengineeringfundamentals Researchgate
Processengineeringfundamentals Researchgate
net/publication/371378445
CITATIONS READS
0 2,076
1 author:
SEE PROFILE
All content following this page was uploaded by Kevin Graham Harding on 08 June 2023.
Kevin Harding
Acknowledgements
The lecturers from the University of the Witwatersrand, Johannesburg, whose notes, time
and/or other materials were considered in compiling this document, as well as the many,
many students that supplied, reviewed and checked solutions:
Disclaimer: This document is a draft of lecture notes compiled from various sources for use
in several chemical engineering courses at the University. The contents are continually
updated and are in no way guaranteed to be correct. Not all sections will be examinable in
the courses this document is made available to, and material not in these notes, but still
lectured, may be examined.
ACKNOWLEDGEMENTS I
TABLE OF CONTENTS II
LIST OF EXAMPLES IX
LIST OF VIDEOS XI
1 NUMERICAL METHODS 1
1.1 INTERPOLATION 1
1.1.1 INTRODUCTION 1
1.1.2 PIECEWISE CONSTANT INTERPOLATION 1
1.1.3 LINEAR INTERPOLATION 2
1.1.4 POLYNOMIAL INTERPOLATION 6
1.1.5 SPLINE INTERPOLATION 7
1.1.6 OTHER FORMS OF INTERPOLATION 8
1.2 EXTRAPOLATION 8
1.2.1 INTRODUCTION 8
1.2.2 LINEAR EXTRAPOLATION 9
1.3 ROOT FINDING 9
1.3.1 INTRODUCTION 9
1.3.2 DIRECT SUBSTITUTION 10
1.3.3 BISECTION METHOD 11
1.3.4 REGULA FALSI 12
1.3.5 ITP (INTERPOLATE, TRUNCATE, AND PROJECT) 16
1.3.6 NEWTON-RAPHSON (DERIVATIVE BASED METHODS) 16
1.3.7 SECANT METHOD 20
1.3.8 MS-EXCEL BASED METHODS 20
1.3.9 OTHER FORMS OF ROOT FINDING. 21
1.4 NUMERICAL INTEGRATION 22
1.4.1 INTRODUCTION 22
1.4.2 PLOT AND COUNT 22
Kevin Harding (kevin.harding@wits.ac.za) ii
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
1.4.3 RECTANGULAR RULE 24
1.4.4 TRAPEZOIDAL RULE 25
1.4.5 SIMPSON’S 1/3 RULE 26
1.4.6 SIMPSON’S 3/8 RULE 27
PROBLEMS 29
2 MASS BALANCES 31
2.1 INTRODUCTION 31
2.2 BATCH, CONTINUOUS AND SEMI-BATCH PROCESSES 31
2.3 STEADY STATE VS UNSTEADY STATE 32
2.4 BASIC MASS BALANCE EQUATION 33
2.5 MATERIAL BALANCES WITHOUT REACTION 35
2.6 MATERIAL BALANCES WITH REACTION 35
2.7 MATERIAL BALANCE FOR MULTI-COMPONENT SYSTEMS 36
2.7.1 RECYCLE SYSTEMS 36
2.7.2 PURGE STREAMS 36
2.7.3 CASCADES 37
2.8 SYSTEMS NEEDING ITERATIVE SOLUTIONS 37
2.9 UNSTEADY STATE MASS BALANCE 38
2.10 DEGREE OF FREEDOM ANALYSIS 39
2.11 GENERAL PROCEDURE TO SOLVE A MASS BALANCE 39
ADDITIONAL READING 40
PROBLEMS 42
3 HEAT TRANSFER 48
3.1 INTRODUCTION 48
3.1.1 WHAT IS HEAT TRANSFER? 48
3.1.2 CONDUCTION 48
3.1.3 CONVECTION 48
3.1.4 RADIATION 49
3.2 CONDUCTIVE HEAT TRANSFER 50
3.2.1 FOURIER’S LAW 50
3.2.2 THERMAL CONDUCTIVITY 52
3.2.3 CONDUCTION THROUGH A WALL 53
3.2.4 CONDUCTION THROUGH A CYLINDER 60
3.2.5 CONDUCTION THROUGH AN IRREGULAR SHAPE 63
3.2.6 CONDUCTION WHEN THERMAL CONDUCTIVITY IS NOT CONSTANT 65
3.2.7 THERMAL RESISTANCE 68
3.3 CONVECTIVE HEAT TRANSFER 70
3.3.1 NEWTON’S LAW OF COOLING 70
3.3.2 CONVECTION FROM A WALL 71
3.3.3 CONVECTION (AND CONDUCTION) THROUGH A CYLINDER 77
Kevin Harding (kevin.harding@wits.ac.za) iii
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
3.3.4 LOG MEAN RADIUS AND LOG MEAN AREA 82
3.4 OVERALL HEAT TRANSFER COEFFICIENT 83
3.5 HEAT TRANSFER THROUGH COMPOSITE WALLS 84
3.5.1 HEAT TRANSFER THROUGH SERIES CONFIGURATION 84
3.5.2 HEAT TRANSFER THROUGH PARALLEL CONFIGURATION 85
3.5.3 CONTACT RESISTANCE 88
3.5.4 POROUS MATERIALS 89
3.6 HEAT TRANSFER FROM EXTENDED SURFACES 90
3.6.1 INTRODUCTION 90
3.6.2 FINS OF UNIFORM CROSS-SECTIONAL AREA 91
3.7 RADIATIVE HEAT TRANSFER 94
BLACK BODY: A PERFECT IDEALIZED PHYSICAL BODY WHICH ABSORBS ALL INCIDENT ELECTROMAGNETIC RADIATION AND
IS ALSO THE BEST POSSIBLE EMITTER OF THERMAL RADIATION. 94
3.8 HEAT EXCHANGERS 95
3.8.1 INTRODUCTION 95
3.8.2 TYPES OF SHELL AND TUBE HEAT EXCHANGERS 97
3.8.3 CO-CURRENT HEAT EXCHANGERS 98
3.8.4 COUNTER CURRENT HEAT EXCHANGERS 98
3.8.5 DERIVATION OF HEAT EXCHANGER EQUATION 99
3.9 HEAT GENERATION 104
3.9.1 NUCLEAR REACTIONS 104
3.9.2 CATALYSTS 104
3.9.3 ELECTRIC CURRENT 105
3.10 UNSTEADY STATE HEAT TRANSFER 107
ADDITIONAL READING 109
PROBLEMS 110
6 PUMPING 175
List of Figures
Figure 1.1: Graphical representation of Piecewise Constant Interpolation (Part 1) .............................. 1
Figure 1.2: Graphical representation of Piecewise Constant Interpolation (Part 2).............................. 2
Figure 1.3: Graphical representation of Piecewise Constant Interpolation (Part 3).............................. 2
Figure 1.4: Graphical representation of linear interpolation................................................................. 3
Figure 1.5: Graphical representation of double linear interpolation .................................................... 4
Figure 1.6: Graphical representation of Polynomial Interpolation (Part 1) ........................................... 6
Figure 1.7: Graphical representation of Polynomial Interpolation (Part 2) ........................................... 7
Figure 1.8: Graphical representation of Spline Interpolation (Part 1) ................................................... 7
Figure 1.9: Graphical representation of Spline Interpolation (Part 2) ................................................... 7
Figure 1.10: Graphical representation of Spline Interpolation (Part 3) ................................................. 8
Figure 1.11: Graphical representation of linear extrapolation .............................................................. 9
Figure 1.12: Graphical representation for bisection method (Part 1) ................................................. 11
Figure 1.13: Graphical representation for bisection method (Part 2) ................................................. 12
Figure 1.14: Graphical representation for Regula-falsi root finding method (Part 1) ......................... 13
Figure 1.15: Graphical representation for Regula-falsi root finding method (Part 2) ......................... 14
Figure 1.16: Graphical representation for Newton Raphson root finding method ............................. 17
Figure 1.17: Situation when Newton-Raphson is not practical ........................................................... 20
Figure 1.18: Screenshot of the Goal Seek tool in MS-Excel ................................................................. 20
Figure 1.19: Screenshot of the Solver tool in MS-Excel ....................................................................... 21
Figure 2.1: Batch process. No material entering or exiting. ............................................................... 31
Figure 2.2: Continuous reactor. Both an inlet and outlet. .................................................................. 32
Figure 2.3: Semi-batch process. In this example, only an inlet. .......................................................... 32
Figure 2.4: Simple representation of a mass balance without reaction ............................................. 35
Figure 2.5: Simple representation of a mass balance with reaction................................................... 36
Figure 2.6: Simple representation of a mass balance with reaction and recycle ............................... 36
Figure 2.7: Simple representation of a mass balance with reaction, recycle and purge .................... 37
Figure 2.8: Simple representation of co-current cascade system with three units ............................ 37
Figure 2.9: Simple representation of counter current cascade system with three units ................... 37
Figure 3.1: Conduction – heat moves from the hot end to the cold end of a metal pole ................... 48
Figure 3.2: Convection – movement of water as it heats in a pot ....................................................... 49
Figure 3.3: Thermal conductivities of various materials ...................................................................... 53
Figure 3.4: Some examples of fins: Rectangular, Pins, Triangular, Parabolic, Annular. ....................... 90
Figure 3.5: Graphical representation of a rectangular fin protruding from a wall .............................. 91
Figure 3.6: Graphical representation of a circular pin protruding from a wall .................................... 91
Figure 3.7: Plate and frame heat exchanger (metal plates transfer heat between two fluids)
(https://commons.wikimedia.org/wiki/File:Plate_frame_1.svg) ........................................................ 96
List of Videos
Video 1.1: Interpolation (YouTube) ....................................................................................................... 1
Video 1.2: Linear interpolation (YouTube) ............................................................................................ 3
Video 1.3: Double interpolation (YouTube) ........................................................................................... 5
Video 1.4: Linear extrapolation (YouTube) ............................................................................................ 9
Video 1.5: Root finding – Direct substitution (YouTube) ..................................................................... 10
Video 1.6: Root finding – Regula Falsi (YouTube) ................................................................................ 15
Video 1. 7: Root finding – Newton Raphson (YouTube) ...................................................................... 18
Video 2.1: Mass balance – General equation (YouTube) ..................................................................... 35
Video 2.2: Mass balance – Ink in a tank example (YouTube) ............................................................... 39
Video 3.3: Heat transfer – Introduction (YouTube) ............................................................................. 49
Video 3.4: Heat transfer – Conduction through slabs (YouTube) ........................................................ 58
Video 3.5: Heat transfer – Conduction through cylinders (YouTube).................................................. 62
Video 3.6: Heat transfer – Parallel configurations (YouTube) ............................................................. 88
Video 3.7: Heat transfer – Conduction through slabs (YouTube) ........................................................ 94
Video 3.8: Heat transfer – Heat exchangers introduction (YouTube) ................................................. 99
Video 3.9: Heat transfer – Heat transfer around a pipe (YouTube) ................................................... 111
Video 3.10: Heat transfer – Counter current heat exchanger (YouTube)........................................... 112
Video 4.1: Intensive vs extensive properties (YouTube).................................................................... 115
Video 4.2: Energy balance – Introduction (YouTube) ........................................................................ 119
Video 4.3: Energy balance – Steam tables (Part 1) (YouTube).......................................................... 121
Video 4.4: Energy balance – Steam tables (Part 2) (YouTube).......................................................... 122
Video 4.5: Energy balance – Mollier diagrams (YouTube) ................................................................ 125
Video 4.6: Energy balance – General equation (YouTube) ............................................................... 128
Video 4.7: Energy balance – Heat of mixing (YouTube) .................................................................... 129
Video 4.8: Energy balance – Heat of reaction (YouTube) ................................................................. 130
Video 4.9: Energy balance – Joule Thomson expansion (YouTube) .................................................. 131
Video 5.1: Fluid flow – U-tube manometer (YouTube) ..................................................................... 138
Video 5.2: Fluid flow – Laminar vs turbulent flow (YouTube)........................................................... 143
Video 5.3: Fluid flow – Reynolds number (YouTube)........................................................................ 144
Video 5.4: Fluid flow – Fanning friction factor (YouTube) ................................................................ 160
Video 5.5: Fluid flow – Pipe roughness (YouTube) ........................................................................... 161
1 Numerical methods
1.1 Interpolation
1.1.1 Introduction
Interpolation is a way of estimating values between a given range of data.
Between each of these data points, individual portions can be distinguished for the x-axis
values closest to these points.
It is then assumed that the y-value in each portion is equal to the data point value in each
of these sub-sections (red lines).
What is the value of H3, at a given temperature (T3) between the given enthalpy values at
T1 and T2?
Temperature Enthalpy
T1 *
H1*
T3 # H3^
T2* H2*
* Given, e.g., from steam tables
# Value to be used between known values, e.g., from the question
^ Unknown to be calculated
For simplicity, and with no other relationship given, we assume linearity across the system
(Linear Interpolation).
" ! #!
! "" " # #!
! ""
= 1.1
$! #$" $# #$"
" ! #!
! ""
$% = %
H $&
& . (T% − T& ) + H 1.2
$ ! #$"
(T2, Ĥ2)
(T3, Ĥ3)
(T1, Ĥ1)
Figure 1.4: Graphical representation of linear interpolation.
NOTE: Points 1 & 2: Given data. Curve: Unknown function connecting known points. Straight line: Linear
approximation. Point 3: Calculated approximation
What is the value of H at P3, T3 (i.e. H7)? Only H at points 1-4 are known.
P H@T1 H@T3 H@T2
P1 H1 H3 = unknown H2
P3 H8 H7 = unknown H9
P2 H4 H6 = unknown H5
(T2, Ĥ2)P1
(T3, Ĥ3)P1
(T3, Ĥ6)P2
(T1, Ĥ4)P2
Figure 1.5: Graphical representation of double linear interpolation
NOTE: Outer points: Given data. Curves: Unknown function connecting known points. Straight lines: Linear
approximations. Middle points (top and bottom): Calculated approximations from curves. Centre point:
Linear approximation from previous points
NOTE: It would also have been acceptable to calculate H8 (from H1 and H4) and H9 (from
H2 and H5), as in the table above, and used this data to calculate H7.
where:
f(_,_): values given at the extremes of the interpolation
x: variables along x-axis
y: variables along y-axis
Example 1.1: Interpolation for internal energy of water at different temperatures and constant pressure
Find the internal energy and specific volume for water at:
i) 6 bar and 330°C
ii) 6 bar and 297°C
Solution:
i) Internal energy: Linear interpolation between 300 and 330°C
(U300 − U350)/(T300 − T350) = (U300 − U330)/(T300 − T330)
(2801 − 2881)/(300 − 350) = (2801 – U330)/(300 − 330)
U300 = 2849kJ/kg
Calculate the specific volume of superheated steam at 8.16 bar and 512°C.
Solution:
The specific volume can be determined by the double interpolation formula:
$! %$ $ %$" ) %) $ %$ $ %$" )%)"
𝑓(𝑥, 𝑦) = ()
$! %$"
* 𝑓(𝑥& , 𝑦& ) + )$ * 𝑓(𝑥& , 𝑦( ), ) !%) + ()$ !%$ * 𝑓(𝑥( , 𝑦& ) + )$ * 𝑓(𝑥( , 𝑦( ), )
! %$" ! " ! " ! %$" ! %)"
where:
x = Temperature
y = Pressure
f(~,~) = Specific volume
For these points, a single polynomial can be fitted that connects as many points as possible.
From this curve, f(x), data points between the original data (triangles) can then be found.
For each data point, 5 individual intervals can be obtained between each point.
Using a cubic spline, for interval 1 to interval 5, 5 independent cubic functions can be
obtained, using the datapoints closest to each interval.
Kevin Harding (kevin.harding@wits.ac.za) 7
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Numerical methods – Extrapolation
1.2 Extrapolation
1.2.1 Introduction
Extrapolation is a way of estimating values beyond a given range of data.
" "
$ % = %!! #!" & . (T% − T& ) + H
H $& 1.5
$ #$ ! "
(T3, Ĥ3)
(T2, Ĥ2)
(T1, Ĥ1)
Figure 1.11: Graphical representation of linear extrapolation
NOTE: Points 1 & 2: Given data. Curve: Unknown function connecting known points. Straight line: Linear
approximation. Point 3: Calculated approximation
methods to solve. However, for more complex equations, a root finding technique may be
needed to. Most of the techniques are iterative in nature, and involve will some error,
depending on the number of iterations performed.
1) Rework the equation such that the left-hand side (LHS) contains the unknown
variable ONLY; to the power one. Move everything else to the RHS, which should
also contain at least one occurrence of the unknown variable.
2) Guess a value of the unknown: x0
3) Substitute this value in the right-hand side (RHS) of the equation and calculate the
unknown on the LHS – x1
4) Use this new value to substitute in the right-hand side and again calculate the LHS
value.
5) Repeat the substitution of the x (xn) into the RHS to calculate a new x (xn+1) until
the xn-1 – xn = 0
The bisection method is a form of a bracketing method of root finding. Two points are
chosen, with successive iterations continually finding a smaller and smaller interval
between each point.
f(x)
+ve
x1
x
x2
-ve
Method:
1) Choose arbitrary value x1 à find f(x1)
2) Choose arbitrary value x2 so that f(x2) is the opposite sign of f(x1)
3) Find x3, equal to (x1 + x2)/2
4) Discard either x1 or x2 based on the sign of f(x3). [Keep x3; discard either of x1 or x2 that has
a corresponding f(x1) or f(x2) with the same sign as f(x3)]
5) Repeat the procedure using the two remaining x values until f(xn) = 0 (or as close as possible)
If f(x3) has the same sign as f(x1), discard x1; else discard x2
If x1 is discarded, work with x2 and x3 (replace all x1 with x3)
If x2 is discarded, work with x1 and x3 (replace all x2 with x3)
f(x)
+ve
x1 x3
x
x2
-ve
-ve
f(x) = [Function]
Iteration Difference
number xn xn + 1 New xn+2 yn yn + 1 yn + 2 New xn New xn + 1 in x-values
[initial
guess to
give
opposite
[initial signs of yn
1 guess] and yn + 1 ] =(B4+C4)/2 = f(xn ) = f(xn + 1 ) = f(xn + 2 ) =D4 =IF(G4*F4 > 0; B4; C4) =ABS(H4-C4)
2 =H4 =I4 =(B5+C5)/2 = f(xn ) = f(xn + 1 ) = f(xn + 2 ) =D5 =IF(G5*F5 > 0; B5; C5) =ABS(H5-C5)
3 [copy row above and paste]
4 [copy row above and paste]
5 [copy row above and paste]
6 [Repeat for as many rows as needed]
Key:
From the question
User input
Final output
f(x)
+ve
x1
x
x2
-ve
Figure 1.14: Graphical representation for Regula-falsi root finding method (Part 1)
Method:
1) Choose arbitrary value x1 à find f(x1)
2) Choose arbitrary value x2 so that f(x2) is the opposite sign of f(x1)
3) Do linear interpolation between these two points [f(x1) and f(x2)]
4) From the linear interpolation, find x3, the point of intersection with the x-axis à find f(x3)
5) Discard either x1 or x2 based on the sign of f(x3). [Keep x3; discard either of x1 or x2 that has
a corresponding f(x1) or f(x2) with the same sign as f(x3)]
6) Repeat the procedure using the two remaining x values until f(xn) = 0 (or as close as possible)
If f(x3) has the same sign as f(x1), discard x1; else discard x2
If x1 is discarded, work with x2 and x3 (replace all x1 with x3)
If x2 is discarded, work with x1 and x3 (replace all x2 with x3)
f(x)
+ve
x1 x3
x
-ve -ve x2
Figure 1.15: Graphical representation for Regula-falsi root finding method (Part 2)
Linear interpolation:
'()! )#'()" ) '()# )#'()" )
= 1.6
)! #)" )# #)"
But f(x3) = y3 = 0
OR
) #)
∴ x% = x& − y& . 2+! #+" 3 1.11
! "
Solution:
x1 = 1 f(x1) = -3
x2 = 2 f(x2) = 5
Since f(x3) is negative, discard x1 (also negative) and use x2 and x3 to find x4
Solution:
a) Starting from x1 = 0.2 and x2 = 0.8:
x1 = 0.2 f(x1) = 5.25
x2 = 0.8 f(x2) = -2.25
𝑥$ − 𝑥# 0.8 − 0.2
𝑥% = 𝑥# − 𝑦# . , - = 0.2 − 5.25. , - = 0.62
𝑦$ − 𝑦# −2.25 − 5.25
\ discard x1
x4 = 0.6517 f(x4) = 0.1633
\ discard x3 (use x2 and x4)
x5 = 0.6617 f(x5) = 0.0553
\ discard x4 (use x2 and x5)
x6 = 0.6650 f(x6) = 0.0183
The Newton-Raphson method starts from an arbitrary guess (x0). It then assumes that the
gradient of the curve at this point can be modeled as a linear line tangential to the curve to
calculate a new x value (x1). This new x is then used to determine a new gradient and linear
approximation line tangential to the curve. This is repeated until a root is found.
f(x)
x2 x0
x1 x
Figure 1.16: Graphical representation for Newton Raphson root finding method
Method:
1) Choose arbitrary value x0
2) Calculate f(x0)
3) Calculate f’(x0)
4) Calculate x1 (assuming a linear tangent line intersecting at y = 0; as show below)
5) Use x1 as new starting point for step (2) above
6) Repeat until f(xn) = 0
+" #'(). )
= f′(x/ ) 1.13
)" #).
But y1 = 0
#'(). )
∴ = f′(x/ ) 1.14
)" #).
'(). )
∴) = f′(x/ ) 1.15
. #)"
'() )
∴ x& = x/ − '0(). ) 1.16
.
OR
'() )
x,-& = x, − '0()* ) 1.17
*
It is a good idea to check that values are smaller in absolute value is successive iterations,
i.e., |f(xn+1)| < |f(xn)|
Solution:
f(x) = x3 + x – 5
f’(x) = 3x2 + 1
𝑓(𝑥/ )
𝑥/0& = 𝑥/ −
𝑓′(𝑥/ )
x0 = 1 f(x0) = -3 f’(x0) = 4
x1 = 1.75 f(x1) = 2.109375 f’(x1) = 10.1875
x2 = 1.5429448 f(x2) = 0.2162004 f’(x2) = 8.1420358
x3 = 1.5163912 f(x3) = 0.003245 f’(x3) = 7.893266
x4 = 1.5159803 f(x4) = 7.68 × 10-7 f’(x4) = 7.894589
x5 = 1.5159802 f(x5) = 4.26 × 10-14 f’(x5) = 7.8945882
x6 = 1.5159802 f(x6) = 0 f’(x6) = 7.89458815
⸫ x = 1.1516
Solution:
$ $%(
𝑓(𝑥) = $%&
+1− ($
(&)($%&)%(&)($) (&)(($)%(()($%()
𝑓 1 (𝑥) = ($%&)!
− (($)!
1 (𝑥) $%&%$ ($%($0-
𝑓 = ($%&)!
− -$ !
%& &
𝑓 1 (𝑥) = ($%&)! − $ !
4($ )
𝑥/0& = 𝑥/ − 41($# )
#
a)
x0 = 0.2 f(x0) = 5.25 f’(x0) = -26.5625
x1 = 0.397647 f(x1) = 2.354637 f’(x1) = -9.0803
x2 = 0.65696 f(x2) = 0.107055 f’(x2) = -10.8148
x3 = 0.666859 f(x3) = -0.00216 f’(x3) = -11.2591
x4 = 0.666667 f(x4) = -8.7 × 10-7 f’(x4) = -11.25
x5 = 0.666667 f(x5) = -1.4 ×10-14 f’(x5) = -11.25
x6 = 0.666667 f(x6) = 0 f’(x6) = -11.25
⸫ x = 0.666667
b)
x0 = -1.5 f(x0) =0.433333 f’(x0) = -0.60444
x1 = -0.78309 f(x1) = -0.33782 f’(x1) = -1.94524
x2 = -0.95675 f(x2) = -0.05625 f’(x2) = -1.35362
x3 = -0.99831 f(x3) = -0.00212 f’(x3) = -1.25381
x4 = -1 f(x4) = -3.2 × 10-6 f’(x4) = -1.25001
x5 = -1 f(x5) = -7.4 ×10-12 f’(x5) = -1.25
x6 = -1 f(x6) = 0 f’(x6) = -1.25
⸫ x = -1
It is possible that there is more than one solution.
c)
x0 = 2 f(x0) = 3 f’(x0) = -1.25
x1 = 4.4 f(x1) = 2.02139 f’(x1) = -0.13816
x2 = 19.031 f(x2) = 1.608006 f’(x2) = -0.00854
x3 = 294.5214 f(x3) = 1.506802 f’(x3) = -2.3 × 10-5
x4 = 6524.43 f(x4) = 1.500031 f’(x4) = -4.7 × 10-10
x5 = 3.21 × 109 f(x5) = 1.5 f’(x5) = -1.9 × 10-19
x6 = 7.73 × 1018 f(x6) = 1.5 f’(x6) = -3.3 × 10-38
x7 = 4.48 × 1037 f(x7) = 1.5 f’(x7) = 1 × 10-75
x8 = 1.51 × 1075 f(x8) = 1.5 f’(x8) = -9 × 10-151
No solution is found.
The Newton-Raphson method is not useful for finding roots that have very flat gradients
or for graphs with discontinuities.
Ther also exist other derivative based methods with higher convergence. These include
the Householder and Halley methods.
Goal Seek
In the Office 365 version of MS-Excel: “Data” ribbon – “What-If-Analysis” (about three-
quarters of the way to the right of the screen) – “Goal Seek”.
Solver
Solver is a more robust than Goal Seek and allows for more manipulation of the tool. It is
also not always installed in MS-Excel by default. To ensure it is available, search for “Add-
ins”, and ensure the “Solver” box is ticked.
Once installed, in the Office 365 version of MS-Excel: “Data” ribbon – “Solver” (often
hiding at the very right of the screen).
Solver can not only obtain a specific value (“Value of”), but also determine maximum
(“Max”) and minimum values (“Min”). The user is also able to manipulate several variables
at ones (“By changing Variable Cells”), set constraints, and change the method of solving.
- Broyden’s method
- Steffenson’s method
- Inverse interpolation
- Brent’s method
- Ridder’s method
Various techniques exist to determine the numerical integrals of data, including the
Newton-Cotes formulae, Gaussian Quadrants and Clenshaw-Curtis Quadrants.
The Newton-Cotes formulae make use of equally spaced intervals and are then divided into
both closed and open types. Closed Newton-Cotes formulae use the function points at all
points while open types do not use the function points at the ends.
Newton-Cotes Closed formulae include the Trapezoidal Rule, Simpson’s Rules and Boole’s
Rule. The Open formulae include the Rectangular Rule, Trapezoidal Rule and Milne’s
Rule.
Given the function below, what is the integral over the given range:
Step one:
Plot on graph paper.
Step two:
Count the number of squares below the graph.
Step three:
Calculate the area of a single block. Convert the number of blocks to an equivalent area.
The smaller the blocks, the more accurate the method.
The concern with this method is that near the graph, it is unclear whether to include or
exclude a block in the counting; particularly if the graph dissects a block.
Original graph:
This can then be solved by either determining the area of each rectangle or by the
mathematical formula below. The formula below calculates the height of the rectangle as
the height at the right-hand side of the rectangle. This results in overcalculation on the
increasing portions of the graph, but underrepresentation on the decreasing sections.
Other formulae may be seen where the left corner, or the middle of the rectangle (midpoint
rule) are taken as the height.
2
∫3 f(x) ≈ h[f(x& ∗ ) + f(x. ∗ ) + ⋯ + f(x, ∗ )] 1.18
where:
h = width of intervals = (b-a)/n [NOTE: Intervals are of equal size]
n = number of intervals
f(xi*) = function value at the midpoint of interval i
Original graph:
This can then be solved by either determining the area of each trapezium or by the
mathematical formula below. How is this mathematically different to the rectangular rule
given above? And for the mid-point rule?
Simplifying:
2 & &
∫3 f(x) ≈ h 2. f(a) + f(x& ) + f(x. ) + ⋯ + f(x,#& ) + . f(x, )3 1.21
where:
a = x value at left most point of concern
h = width of intervals [NOTE: Intervals are of equal size]
n = number of intervals
The simplest of these is Simpson’s 1/3 Rule which uses a quadratic formula, which gives
the following integral estimation:
where:
h = width of intervals [NOTE: Intervals are of equal size]
n = number of intervals
NOTE:
Outer function values have co-efficient 1
All even function values have a co-efficient of 4
All odd function values have a co-efficient of 2
For the Simpson’s Rule (also called the 1/3 Simpson’s Rule) we need an even number of
integrals, i.e., an odd number of points.
where:
h = width of intervals [NOTE: Intervals are of equal size]
n = number of intervals
NOTE:
Outer function values have co-efficient 1
Kevin Harding (kevin.harding@wits.ac.za) 27
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Numerical methods – Numerical Integration
For the Simpson’s 3/8 Rule we need the number of integrals to be equal to a multiple of 3.
Problems
Example 1.7: Vapour pressure of 1-chlorotetradecane 1.0bar 1.2bar 1.5bar 𝑓 1(0.971762 … ) = 2𝑒 "($.'2#23"… ) − 1 = 12.9966654 …
300°C 2811 kJ/kg A 2809kJ/kg 0.011565 …
𝑥" = 0.971762 … − = 0.970871 …
12.966654 …
Vapour pressures of 1-chlorotetradecane (C14H29Cl) at several temperatures are 330°C C
tabulated below: 400°C 2968 kJ/kg B 2967 kJ/kg Iteration 3:
𝑓(0.970871 … ) = 𝑒 "($.'2$)2#… ) − 0.970871 … − 6 = 1.11 × 10&*
T (°C) 98.5 131.8 148.2 166.2 199.8 215.5 A: At 300°C: 𝑓 1(0.970871 … ) = 2𝑒 "($.'2$)2#… ) − 1 = 12.94176 …
P* (mm Hg) 1 5 10 20 60 100 (UA − 2811)/(2811) = (1.2 − 1.0)/(1.0 − 1.5) 1.11 × 10&*
𝑥5 = 0.970871 − = 0.970871 …
UA = 2810.2 kJ/kg 12.94176 …
i) Calculate P* for a temperature of 483.15K
ii) 1-chlorotetradecane boils at approximately 140°C and behaves similarly to an B: At 400°C: Since x2 is very close to x3, we can stop iterating: Final x value calculated: 0.970871
ideal gas. Determine the mass of the gas of 1-chlorotetradecane at 142.3°C in a (UB − 2968)/(2968 − 2967) = (1.2 − 1.0)/(1.0 − 1.5)
50 dm3 container. UB = 2967.6 kJ/kg Example 1.10: Enthalpy of steam at 17 bar
Solution: C: At 330°C and 1.2bar: Find the enthalpy (H) of steam at 17 bar and 340°C.
"#$ & #''.)
i) 𝑃∗ = 60 + (100 − 60) -
"#*.* & #''.)
. = 85.99 𝑚𝑔 𝐻𝑔 (UC − 2810.2)/(2810.2 − 2967.6) = (330 − 300)/(300 − 400)
UC = 2857.42 kJ/kg Solution:
ii) P* at T = 142.3°C From steam tables (superheated steam) we can find the enthalpy values as below
Updated table with calculated internal energy values for 1.2 bar (kJ/kg). With x, y & z the values to obtain to solve the problem:
148.2 − 131.8
𝑃∗ = 5 + (10 − 5) 7 ; = 8.20 𝑚𝑚 𝐻𝑔 → 0.0108 𝑎𝑡𝑚 1.0bar 1.2bar 1.5bar
148.2 − 132.8
300°C 2811 kJ/kg UA = 2810.2 kJ/kg 2809kJ/kg 300°C 340°C 350°C
From the ideal gas law: 330°C UC = 2857.42 kJ/kg 15 bar 3039 x 3148
𝑃𝑉 = 𝑛𝑅𝑇 400°C 2968 kJ/kg UB = 2967.6 kJ/kg 2967 kJ/kg 17 bar z
𝑃𝑉 (0.0108). (50)
𝑛= = = 0.0158 𝑚𝑜𝑙 20 bar 3025 y 3138
𝑅𝑇 (8.206 × 10&" )(415.45)
Example 1.9: Newton-Raphson for a function involving e
Converting to mass: At constant 15 bar, interpolate across temperatures:
𝑀𝑜𝑙𝑎𝑟 𝑚𝑎𝑠𝑠 = 14(12.01) + 29(1.008) + 35.45 = 232.822 𝑔/𝑚𝑜𝑙 Using Newton-Raphson method, find all solutions of the following equation: e2x = x + 𝐻5$$℃ − 𝐻5*$℃ 𝐻/ − 𝐻5$$℃
=
𝑇5$$℃ − 𝑇5*$℃ 𝑇57$℃ − 𝑇5$$℃
6. Use an initial value of 1 for x.
Therefore: 3039 − 3148 𝐻/ − 3039
=
𝑀𝑎𝑠𝑠 = (𝑛𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑚𝑜𝑙𝑒𝑠). (𝑚𝑜𝑙𝑎𝑟 𝑚𝑎𝑠𝑠) = (0.0158)(232.822) = 3.688 𝑔 300 − 350 340 − 300
Solution:
-(/ ) 𝐻/ = 3126.2 𝑘𝑔/𝑘𝑔
Newton’s method: 𝑥+,# = 𝑥+ − -" (/! )
Example 1.8: Internal energy of water !
2 Mass balances
2.1 Introduction
The typical difference between the interests of the chemist and the chemical engineer is
that the latter is concerned with continuous processes. These are systems into which there
is a continuous flow of feed material, some (at least) of which is converted to a desired
product, which then flows continuously out of the system, together with any by-products
and unused feed material. The basic calculations performed on such a process (as distinct
from the design of the individual units) are the mass balance, which is considered below,
and the energy balance.
A continuous process is one that allows for inputs and outputs to flow in and out of the
unit all the time, with reaction or transformation happening while material is added (raw
materials) and removed (product). This is a steady state process (Figure 2.2).
Examples: Many large scale industrial production facilities that run 24-hours a day, e.g.
petroleum refining.
A semi-batch process, also called a fed-batch process, is a combination of these two systems,
where consideration is made for either an input or an outputs stream (not both). This
process is also considered unsteady state (Figure 2.3).
:
=0 2.1
:;
where:
t = time
A system is said to be under unsteady or transient state (sometimes also non-steady state),
if any variable is changing at a specific location, with time.
:
≠0 2.2
:<
Solution:
a) Steady state. Assuming he is not changing the flow at the tap.
b) Unsteady state. The volume increases with time, and the temperature in the bath may also change due
to hot water taking time to heat. There is no flow out to balance temperature or volume.
c) Unsteady state: A geyser is designed to keep the temperature constant. However, the typical operation
of the thermostat relies on a constant on/off cycling where the heat is supplied only when the
temperature drops below a certain point and will again switch off when the water reaches the higher
set point. Since there is constant change in temperature over time, unsteady state.
d) Steady state. As long as the inlets do not vary in temperature or flowrate, even though the outlet is a
different temperature, at any specific location, the flow and temperature are constant.
e) If we consider the engine alone.
At startup: Unsteady state. The engine will start to warm up, i.e., temperature changes over time.
After a while: Steady state. Burning constant amount of fuel will result in a constant rotation, with
any temperature generated, dissipated constantly given sufficient time.
If we consider the system as the engine AND the fuel tank.
Unsteady state: The volume of the fuel tank will decrease over time, resulting in the engine cutting
out once the fuel is completely consumed.
In
Out
In = Out 2.3
∴ Ṁ=, = Ṁ>?@ 2.4
where:
Ṁ = mass flow rate (kg/s)
This can also be rewritten for reaction/consumption of a single species in the system.
In Generated Out
Consumed
Rearranging:
ṀE>,F?GBD − ṀAB,BC3@BD = Ṁ=, − Ṁ>?@ 2.6
Where:
ṀE>,F?GBD − ṀAB,BC3@BD = mass change inside the system per time.
D
NOTE: This is no longer a steady state system, i.e., D@ ≠ 0
where:
M = Mass inside the system boundary
t = time
Kevin Harding (kevin.harding@wits.ac.za) 34
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Mass balances – Material balances without reaction
where:
i = stream number. Where i is positive for inlet streams and negative for outlet stream
n = total number of streams
A A, Aʹ
B B, Bʹ
*
A, A
*
B, B
Figure 2.4: Simple representation of a mass balance without reaction
Where: A and B are different non-reacting components, and Aʹ, A*, Bʹ and B* are different
forms (phase, size class or similar) of each respective component
often make it convenient to express the amounts of the various materials in moles rather
than in mass; at the same time, we need to note that there is no conservation of (total)
moles as there is of mass due to the stoichiometry of the chemical reaction that is taking
place. For example, if A + B à C, we have two moles total entering, but only one leaving.
A A
B Reactor B
C
Example 2.2:
Solution:
A A
B Reactor Separator B
C
Figure 2.6: Simple representation of a mass balance with reaction and recycle
Where: A and B are different reacting components to form a single product C (A + B à C).
A, B, C, I
A
A
B
B Reactor Separator
C
I
I
Figure 2.7: Simple representation of a mass balance with reaction, recycle and purge
Where: A and B are different reacting components to form a single product C (A + B à C), and
I is the inert.
2.7.3 Cascades
A cascade is a multi-step process that can act as a reactor or separator. Flows from one unit
proceed to successive units in either a co-current (Figure 2.8) or a counter current fashion
(Figure 2.9).
Figure 2.8: Simple representation of co-current cascade system with three units
Figure 2.9: Simple representation of counter current cascade system with three units
In very complex cases, especially when there are interactions between the mass and energy
balances, we resort to simply calculating "round and round" the various loops until
convergence is obtained. Computer packages, including MS-Excel, are available to perform
these calculations, including the modelling (or design) of the process units.
Water enters and leaves a 5 m3 tank at a rate of 0.12 m3/min. At time t = 0, 1 kg ink is injected into the
tank. Assume that the ink is well mixed, and evenly distributed. Complete the mass balance.
Solution:
Mass balance of the water:
/
𝑑𝑀
= E 𝑀̇5
𝑑𝑡
56&
But at steady state, and only two streams:
𝑀̇5/ − 𝑀̇789 = 0
𝑀̇5/ = 𝑀̇789
𝑤ℎ𝑒𝑟𝑒, 𝑀̇5/ = 0.12 𝑚, /𝑚𝑖𝑛
Mass balance of the ink:
/
𝑑𝑀
= E 𝑀̇5
𝑑𝑡
56&
𝑑𝑀
= 𝑀̇5/ − 𝑀̇789
𝑑𝑡
𝑑𝑀
= 0 − 𝑀̇789
𝑑𝑡
𝑑𝐶. 𝑉
= −𝑀̇789
𝑑𝑡
𝑑𝐶
𝑉. = −(𝑣. 𝐶)789
𝑑𝑡
𝑑𝐶
𝑉. = −𝑣789 𝑑𝑡
𝐶789
𝑑𝐶
𝑉. Q = −𝑣789 Q 𝑑𝑡
𝐶789
𝐶789
𝑉. 𝑙𝑛 = −𝑣789 (𝑡 − 𝑡5 )
𝐶789,5
At i: ti = 0, Ci = 0.2 kg/m3
V = 5 m3
vout = 0.12 m3/min
Substituting values:
𝐶789
5. 𝑙𝑛 = − 0.12(𝑡 − 0)
0.2
𝐶789 0.12
𝑙𝑛 = −0 7.𝑡
0.2 5
𝑙𝑛(𝐶789 ) − 𝑙𝑛(0.2) = −0.024. 𝑡
If DOF > 0, the problem is underspecified. This means that there are more unknown
variables than the number of independent relations and hence, cannot be solved. One way
of fixing this is by merging multiple data into a single solution.
DOF < 0, the problem is over-specified. This means that there are possibly inconsistent
independent relations, therefore, it is unsolvable. Gather additional information until the
problem is well defined.
Decide on a basis
If no basis is given, and the DOF = 1, decide on an appropriate basis of calculation to bring
the DOF to 0. It may be useful to set the basis as the input to the system or the input to the
reactor (after mixing from recyle), but how to do this will be come easier with practice.
Additional Reading
Felder RM, Rousseau RW, 2000. Elementary Principles of Chemical Processes, 3rd ed.,
Wiley.
Himmelblau, DM, 2012. Basic Principles and Calculations in Chemical Engineering, 8th
ed., Prentice-Hall International.
Murphy, RM, 2007. Introduction to Chemical Processes: Principles, Analysis, Synthesis,
McGraw-Hill International Edition.
Reklaitis, GV, 1983. Introduction to Material and Energy Balances, Wiley.
Kevin Harding (kevin.harding@wits.ac.za) 40
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Mass balances – Additional Reading
Skogestad, S, 2009. Chemical and Energy Process Engineering, CRC Press, Taylor & Francis
Group.
Problems
Example 2.4: Combustion in a furnace ii) Sulphur, hydrogen and nitrogen do not react. ∴ Feed = Output
Mass flow rate of sulphur, in stream 1 and 4, using moles of carbon in the feed: Substitute back into Equation 5: E1 = 10(0.71) = 7.1 kmol/hr
Coal (91% carbon, 4% sulphur, 4% hydrogen and 1% oxygen) is mixed with air (78% m1,S = [(m1,C)/(X1,C)] × (X1,S)
nitrogen and 22% oxygen by volume) (as per the diagram), at STP, and is pumped into = [(120 kg/hr)/(0.91)] × (0.04) = 5.27kg/hr = m4,S Example 2.5: Flow in a river
a furnace to produce CO2 gas. Oxygen is 28% in excess and has an overall conversion
of 92%. There is also CO formed in a side reaction with selectivity of 5. Given that the Mass flow rate of hydrogen, in stream 1 and 4, using moles of carbon in the feed: A colleague has been working at a local river collecting data on water quality and other
ratio of oxygen in coal and in air is 1/10, and that 55 kmol/hr of nitrogen is coming out m1,H = [(m1,C)/(X1,C)] × (X1,H) = [(120 kg/hr)/(0.91)] × (0.04) = 5.27kg/hr = m4,H related things. They would like to find out the total flowrate in the river and have
of the furnace, calculate: asked if you can help given the following messy data:
i) Mass of carbon in the coal feed per hour. Mass flow rate of nitrogen, in stream 2 and 4, using n2,N:
ii) Mass fraction of sulphur in stream 4. m2,N = (n2,N) × (MN) = (55 kmol/hr) × (14 kg/kmol) = 770 kg/hr = m4,N [1: 0: 0.00], [2: 9: 0.06], [3: 39: 0.22], [4: 38: 0.16], [5: 50: 0.18], [6: 42: 0.15], [7: 38: 0.12],
iii) Extent of each reaction. [8: 29: 0.10], [9: 23: 0.07], [10: 22: 0.02], [11: 0: 0.00].
Fraction of sulphur in stream 4:
Reactions: C (s) + O2 (g) → CO2 (g) 2C (s) + O2 (g) → 2CO (g) X4,S = (m4,S)/(m4) = (m4,S)/(m4,S + m4,H + m4,N) Each entry in the data set contains a reference point starting at the riverbank, a depth
= (5.27)/(5.27 + 5.27 + 770) = 0.68/100 = 0.68% (cm) and velocity (m/s). Each successive entry is 500 mm away from the previous point
1: Air 3: CO, CO2, O2 across the river. What is the flowrate across the river assuming both the depth and
Furnace
iii) Elemental mole balances to get extents. velocity is constant in each of the divisions above? (Give your answer to the nearest
CO2 mole balance: 0 + E1 = n3,CO2 + 0 à E1 = n3,CO2 (1) litre.)
2: Coal 4: S, H2, N2
CO mole balance: 0 + 2E2 = n3/CO + 0 à 2E2 = n3,CO (2)
Oxygen mole balance: Solution:
Solution: n1,O + n2/O + 0 = n3/O + 2(E1 + E2) à n1/O + n2/O = n3/O + 2(E1 + E2) It was given than each point was 500 mm (0.5 m) from the previous point. Therefore
i) Nitrogen does not react, so n2,N = n4,N = 55 kmol/hr (see table):
For a gas at STP: volume% is the same as mole%: Feed ratio of oxygen to get n1,O: n1,O = (n2/O)/10 = (15.51 kmol/hr)/10 = 1.55 kmol/hr Values in column 2 (Distance from bank) can be calculated by adding 500mm to the
n2,O = (n2,N / X2,N) × (X2,O ) = (55/0.78) × (0.22) = 15.51 kmol/hr previous point, starting at zero (note: it makes no difference which end of the table you
Overall conversion of oxygen to get n3/O: start working this out)
The ratio of carbon to hydrogen reacting can be obtained using the sum of the reactions: Overall Conversion % = (moles reacted)/(moles in feed) Values in column 3 (Gap) were given and added again for simple MS-Excel calculation
C (s) + O2 (g) → CO2 (g) Overall Conversion % = [(moles in feed)-(moles in output)]/(moles in feed) Values in column 4 (Depth – cm) were given
+ 2C (s) + O2 (g) → 2CO (g) (92/100) = [(n1,O + n2,O) – (n3,O)] / (n1,O + n2,O) Values in column 5 (Depth – m) are a unit conversion
3C (s) + 2O2 (g) → CO2 (g) + 2CO (g) = [(1.55 + 15.51) – (n3,O)] / (1.55 + 15.51) Values in column 6 (Velocity) were given
∴ 3 Moles of Carbon reacts with 4 moles of Oxygen. n3,O = 1.36 kmol/hr Values in column 7 (volumetric flow) are the values from column 5 multiplied by
column 6.
Using % excess to get moles of carbon in the feed: Substitute n1,O, n2,O and n3,O into oxygen mole balance:
% Excess = [(1/4)(nfeed) – (1/3)(nstoch)] / (1/3)(nstoch) n1,O + n2,O = n3,O + 2(E1 + E2) Distance Gap Depth Depth Velocity Volumetric
(20/100) = [(1/4)(n1,O + n2,O) – (1/3)(n1,C)] / (1/3)(n1,C) 1.55 + 15.51 = 1.36 + 2(E1 + E2) from (m) (cm) (m) (m/s) flowrate(m3/sec)
7.85 = E1 + E2 (3) riverbank [GIVEN] [GIVEN]
(m)
Using ratio of oxygen in coal and air to get n1,O:
1 5 0.5 0 0 0 0
(n1,O/n2,O) = (1/10) Selectivity = (moles of desired product)/(moles of undesired product)
2 4.5 0.5 9 0.09 0.06 0.0027
n1,O = (n2,O/10) -Substitute this and the value for n2,O into excess equation 5 = (n3,CO2)/(n3,CO) (4) 3 4 0.5 39 0.39 0.22 0.0429
4 3.5 0.5 38 0.38 0.16 0.0304
(20/100) = [(1/4)((n2,O/10) + n2,O) – (1/3)(n1,C)] / (1/3)(n1,C) Substitute Equations 1 and 2 into Equation 4: 5 3 0.5 50 0.5 0.18 0.0450
6 2.5 0.5 42 0.42 0.15 0.0315
(20/100) = [(1/4)((15.51/10) + 15.51) – (1/3)(n1,C)] / (1/3)(n1,C) 5 = (E1)/(2E2)
7 2 0.5 38 0.38 0.12 0.0228
n1,C = 10 kmol/hr E1 = 10E2 (5)
8 1.5 0.5 29 0.29 0.1 0.0145
9 1 0.5 23 0.23 0.07 0.0081
Converting to mass units Substitute Equation 5 into Equation 3: 10 0.5 0.5 22 0.22 0.02 0.0022
m1,C = (n1,C) × (MC) = (10 kmol/hr) × (12 kg/kmol) = 120 kg/hr 7.85 = (10E2) + E2 11 0 0.5 0 0 0 0
E2 = 0.71 kmol/hr
The total flow rate is the sum of all values in column 7. This equals 0.20005m3/s n1,B = X1,B × n1 = 0.27 × 100 = 27 mol/min ṁ3, E = 98 kg/hr
Therefore, converting to litres: 200 l/s
Stoichiometric ratio of A:B is 4:1, and the feed ratio of A:B is 73:27 = 2.70:1 Since weight percent of ethanol is known in stream 3:
Example 2.6: Formation of milk of lime \Fe3O4 is limiting ṁ3, E = xe × ṁ3
98 kg/h = 0.05 × ṁ3
Limestone Ca(OH)2 has a solubility of 2 g/L in water. If 63 kg/hr of hydrated lime is ii) Stoichiometric requirement i.e. nStoich of B to react with all A present ṁ3 = 98 / 0.05 = 1960 kg/hr
added to 520 kg/hr of water to form milk of lime, find the percentage of lime that nStoich = n1,A × 1/4 = 73 × 1/4 = 18.25 mol/min
dissolves and the mass of solid in the milk of lime? %excess of B = (n1B - nStoich B)/ nStoich B = (27 – 18.25) × 100/18.25 = 47.95% And the amount of water in stream 3: ṁ3, W = ṁ3 – ṁ3, E = 1960 – 98 = 1862 kg/hr
If the SG of lime is 2.5, calculate the density of milk of lime.
iii) %conversion of A = 78% Mass flow rate of water in stream 2(water balance):
63 kg/hr lime n2,A = n1,A – 0.78 n1,A = 73 – 0.78(73) = 73 – 56.94 = 16.06 mol/min ṁ1, W = ṁ2, W + ṁ3, W
Milk of lime
ṁ2, W = 2275 – 1862 = 413 kg/hr
Component A balance: nOut A = n2,A = nIn A + nFormed A – nReacted A
520 L/hr water 16.06 = 73 + 0 – 4ε Example 2.9: Production of ethylene
-4ε = 16.06 – 73 = -56.94
Solution: \ε = 14.235 mol/min Ethylene (E) (C2H4) is produced by the following reactions, including hydrogen (H)
Lime dissolved in 520 L/hr: 2 g/L × 520 L/hr = 1.04 kg/hr and methane (M) side reactions:
The percentage of lime that dissolves is: (1.04/63) × 100= 2% iv) Component balance on C: nOut C = nIn C + nFormed C – nReacted C C2H6 à C2H4 +H2 (1)
Therefore, 98% of lime does not dissolve = (98/100) × 63 = 61.7 kg nOut C = 0 + 6ε – 0 = 6 × 14.235 = 85.41 mol/min C2H4 + H2 à 2CH4 (2)
\ Production rate = nOut C = 85.41 mol/min
1/density of mixture (AB) = 1/A + 1/B = 1/2500 + 1/1000 =7/5000 = 714.3 kg/m3 A feed with 89% percent ethane (P) with inerts (I), is fed to the reactor, where the
v) n2 = Σ n2i = n2A + n2B + n2C overall conversion of ethane is 65%. Given a selectivity for ethylene is 0.9; calculate
n2,A = 16.06 mol/min the extent of reactions (ɛ1 and ɛ2) and yields (Y) for both reactions, and the reactor
Example 2.7: Solid phase reaction with no recycle
n2,B = n1B + nFormedB – nReactedB = 27 + 0 – 1 × ε = 27 + 14.235 = 12.765 mol/min product composition.
Iron (III) oxide (Ferric oxide or hematite) has many industrial uses, including jewellery
n2 = n2,A + n2,B + n2,C = 16.06 + 12.765 + 85.41 = 114.235 mol/min Solution:
polishing, pigmentation, magnetic media storage, catalysis, and others. The Gamma
X2,A = n2,A / n2 = 16.06/114.235 = 14.06% Basis: 100 kmol/hr feed
phase of iron (III) oxide (γ-Fe2O3) is the specific type used in magnetic storage media
X2,B = n2,B / n2 = 12.765/14.235 = 11.17%
such as tapes, CD’s, and hard disk drives. This phase of iron (III) oxide can be prepared
X2,C = n2,C / n2 = 85.41/114.235 = 74.77% I balance: n2,I = n1,I = x1,I n1 = (0.11)(100 kmol/hr) = 11 kmol/hr
by the careful oxidation of Fe3O4.
+#,%&+&,%
Overall Reaction: 4Fe3O4(s) + O2(g) à 6Fe2O3(s). The feed stream is 73% Fe3O4 by Example 2.8: Distillation of water and ethanol Overall %Conversion E = = 0.65
+#,%
Stoichiometric equivalent of oxygen = 500/1 = 500 kmol/hr In this process, a feed enters the reactor composed of 35% H2S and CH4. A second feed
Ethylene is the limiting reactant since it has the smaller stoichiometric equivalent. stream enters the reactor composed of pure SO2 gas. The process yields 2500 kg pure 5 n3 H2S = n3 H2O = 51.978 kmol/hr
solid sulphur. n3 H2S = 10.3956 kmol/hr
%excess = [(nO2)feed × (nO2)stoic] / [(nO2)stoic ] × 100
Stoichiometric ratio: [(nO2) stoic] / [(nC2H4)] = ½ 1: 35% H2S n3 SO2 = 2 n3 H2S = 2(10.3956) = 20.7912 kmol/hr
(nO2) stoic = 200/2 = 100 kmol/hr. 65% CH4
3 5: SO2, CH4, H2S, H2O
Therefore, %excess = [(nO2)feed × (nO2)stoic] / [(nO2)stoic ] × 100 = [(500 – 100)/100] × 100 = Sulfur Dioxide balance:
Reactor
400% Out = In + Formed – Reacted
2: 100% SO2
20.7912 = m2 –25.989
2) (i) 4: 2 500 kg S(s) m2 = 46.7802 kmol/hr
O2 balance about the reactor:
Out = In + Formed – Reacted Hydrogen Sulfide balance:
(nO2)out = 500 + 0 – Ɛ i) Perform a Degree of Freedom analysis. Out = In + Formed – Reacted
ii) Find the flow rates of the feed streams and the single pass conversion of the 10.3956 = 0.35m1 – 2(25.989)
C2H4 balance about the reactor: limiting reactant m1 = 178.2111 kmol/hr
Out = In + Formed – Reacted
0 = 200 + 0 – (nC2H4)reacted Solution: ∴ The limiting reactant is SO2
Ɛ = 200/2 =100 kmol/hr (i) Degrees of Freedom:
8 Equations: Single pass conversion = Ɛ/m2 = 25.989/46.7802 = 0.555555 = 55.56%
From O2 balance: (nO2)out = 500 – 200 = 300 kmol/hr. • m1 = m1(x1 H2S) + m1(x1 CH4)
• m2 = (x2 SO2) m2 Example 2.14: Vibrating screen separator
C2H4O balance: • m3SO2 = (x3 SO2)m3
Out = In + Formed – Reacted • m3H2S = (x3 H2S)m3 A vibrating screen separator consists of 2 decks. Granulated ore, consisting of large,
(nC2H4O)out = 0 + (2 × 100) – 0 = 200 kmol/hr. • m3H2O = (x3 H2O)m3 medium, and small particles, is separated into 3 products in a continuous steady state.
• m3CH4 = (x3 CH4)m3 The wire mesh of each deck produces an overflow (particles too large to pass through
(ii) • m4 = (x4 S)m4 the deck’s openings) and an underflow (particles that pass through the deck’s openings).
Fractional conversion = [(nO2) reacted]/ [(nO2) feed] = 0.7 • m1 + m2 = m3 + m4 The feed to the screen is in a ratio of 2:1:1 of the large, medium, and small particles.
Hence, (nO2)reacted = 0.7 × 200 = 140 The vibrating separator is not 100% efficient, as particles that should be part of the
Ɛ = 140/2 =70 kmol/hr 8 Variables: underflow product are often included in the in the overflow products. 89% of the small
• x1 H2S particles report to the second deck’s underflow. 85% of the medium particles report to
C2H4 balance: • x1 CH4 the second deck’s overflow. The composition of the second deck’s overflow consists of
Out = In + Formed – Reacted • x2 SO2 10% smalls.
(nC2H4)OUT =200 + 0 – 140 = 60 kmol/hr • x3 SO2
• x3 H2S Feed
O2 balance: • x3 H2O
50% large (L), 25% medium (M),
Out = In + Formed – Reacted • x3 CH4 25% small (S) 1
(nO2)out = 500 + 0 – 70 = 430 kmol/hr • m3 2 Top Overflow
∴ 0 degrees of freedom
First Deck
Example 2.13: Mass balance with a reaction for removing impurities
Raw gases which contain hydrogen sulphide or carbon dioxide are considered as sour (ii) Basis: 2 500 kg of solid sulfur in an hour -à 77.967 kmol/hr
3 Bottom Overflow
gases. These gases need to go through a process for their impurities to be removed.
Sulfur balance: Second Deck 85% of all medium particles
Mercaptans, hydrogen sulphide, and other sulphur compounds are removed from
natural gas by a “sweetening process” producing less impure sour gases. The desired Out = In + Formed – Reacted X3S = 10%
condition in this specific process is the removal of H2S gas, by using the following 77.967 = 0 + 3Ɛ 4
equation: Ɛ = 25.989 kmol/hr Bottom Underflow 89% of all Small particles
ii) Calculate the percentage composition of the particles in the first deck’s %x2S = (m2S/m2) × 100 = (0.39/54.14) × 100 = 0.72% Reactor H balances:
overflow. Out = In + Formed - Reacted
Example 2.15: Formation of milk of lime 𝑛"N = 𝑛#N + ᶓ" + ᶓ#
Solution: 0.8827*32.5 = 0 + 2.223 + ᶓ#
Basis: 100 tph feed flowrate. A simplified process used to produce monochloroethane and hydrochloric acid is ∴ ᶓ# = 26.465 kmol/min
i) To find: m2L, m2M, m2S shown below. Ethane is chlorinated in a continuous reactor by reaction 1:
Using the large, medium, small ratio of 2:1:1, the flowrates of stream 1 can be C₂H6 + Cl2 → C2H5Cl + HCl. From the given information ethane in the feed stream can be calculated and with the
calculated. reactive balance around the reactor we can find M and E in the reactor product.
X1L = (100/4) × 2 = 50% The reactor is deliberately used for low conversion of ethane to minimise the further Conversion information: 68%
m1L = 100 × 50% = 50 tph chlorination of the product C₂H₅Cl in the reactor by the undesirable side reaction 2:
++,-./,0/'2304.,0
C₂H₅Cl + Cl₂ → C₂H₄Cl₂ + HCl Percentage conversion = +5,,0
X1M = (100/4) × 1 = 25% 0.68 = +
ᶓ#
m1M = 100 × 25% = 25tph The reactor product reports to the first separator where only C₂H₅Cl is removed by #%
"3.73*
𝑛#P = = 38.919 kmol/min
stream 3. The second separator product (stream 6) has a flowrate of 32.5 kmol/min and $.3)
X1S = (100/4) × 1 = 25% its composition is 4.89 mole% chlorine, 6.84 mole% C₂H₄Cl₂ and the rest is HCl.
m1S = 100 × 25% = 25 tph Separator 2 removes all the HCl and C₂H₄Cl₂. The single pass conversion of ethane is Reactor E balance:
68%. Out = In +Formed-Reacted
Using recovery rate of small particles to the second deck underflow: 𝑛"P = 𝑛#P + 0 - ᶓ# = 38.919 – 26.465 = 12.454 kmol/min
m4S = 25 × 89% = 22.25tph Calculate the composition of reactor product.
Reactor M balance:
Total stream 4 flowrate is 22.25tph since: 3 Out = In + Formed –Reacted
m4L = 0 tph 𝑛"Q = 𝑛#Q + ᶓ# - ᶓ" = 0 + 26.465 – 2.223 = 24.242 kmol/min
SEPARATOR
m4M = 0 tph (screen only allows small particles through, medium particles are separated SEPARATOR
1 2
to stream 3)
1 2
4 6 To write a Cl balance around the reactor we need the flowrate of Cl in the feed.
Overall Cl balance:
Using recovery rate of medium particles to the second deck overflow: Out = In +Formed-Reacted
m3M = 25 × 85% =21.25tph REACTOR 𝑛3R? = 𝑛#R? - ᶓ# - ᶓ"
5
0.0489*32.5 = 𝑛#R? – 26.465 – 2.223
∴ 𝑛#R? = 30.277 kmol/min
KEY:
Because composition of m3 consists of 10% smalls and m3L=0tph: • E- Ethane
• Cl- Chlorine
Liquid:
0.04 Benzene
3 Heat transfer
3.1 Introduction
3.1.1 What is heat transfer?
Heat transfer is the process of the movement of energy due to a temperature difference.
The calculations we are interested in typically include determining the final temperatures
of materials or how long it takes for these materials to reach these temperatures. This can
help inform the level of insulation required to ensure heat is not lost from a system.
Typically, heat loss is proportional to a temperature gradient (driving force or potential).
3.1.2 Conduction
Conduction is the form of heat that exists due to direct contact without movement. A
temperature gradient within a substance causes a flow of energy from a hotter to a colder
region. These gradients can exist in solids, liquids, and gases; provided there is no
movement in the fluid phases, i.e., fluids which are not well mixed. Over time the
temperature difference will reduce and approach thermal equilibrium (same temperature).
Conduction occurs in a solid, liquid or gas; provided there is no bulk movement.
Examples of conduction include the end of a metal rod placed in a fire heating up from one
end to the other, hot coffee heating through the mug or the cooling of ice cream in the
bowl it is placed in.
Figure 3.1: Conduction – heat moves from the hot end to the cold end of a metal pole
3.1.3 Convection
Convection is the transfer of heat due to the bulk movement of fluids. As such convection
only applies to heat transfer within a fluid or between a solid and fluid but not the heat
transfer within a solid. This heat transfer is achieved by the movement of molecules within
the fluid. The term convection can refer to either mass transfer and/or heat transfer.
Typically, when referred to as ‘convection’, heat transfer is meant.
Examples of convection include the effect of hot air rising and falling (convection currents)
or the large-scale convection currents of the atmosphere and oceans.
3.1.4 Radiation
Radiation is the transfer of energy due to electromagnetic waves when thermal energy is
converted by the movement of the charges of electrons and protons in the material. When
a body radiates, the energy comes from the entire depth of the body, not just the surface.
Radiation does not require a temperature gradient. A person standing some distance from
the source will still feel the effects of the heat, e.g., a person near a fire is heated by the fire,
not by the air surrounding them.
Examples include infra-red radiation such as, an incandescent light bulb emitting visible
light, the infrared radiation emitted by a common household radiator or electric heater, as
well as the sun heating the earth.
The power is going to be down for the next 5 hours. You have just enough time to boil a cup of water for
a cup of coffee. However, you don’t really want a cup of coffee right now – but would rather have it in
30 min time. Your colleagues come up with several suggestions on how to ensure the coffee is as warm as
possible when you want to drink it. Unfortunately, they are not sure if these will make the coffee hotter
or colder. For each suggestion indicate which would keep the coffee hotter for when you want to drink
it. Explain why.
a) Pouring the (cold) milk in when the coffee has been made OR pouring the milk in just before drinking.
Assume the milk stays at the same temperature even with the power out due to good insulation of the
fridge.
b) Covering the top of the coffee cup with something OR leaving the top of the cup open to the
atmosphere.
c) Using a steel mug OR a ceramic mug.
d) Leaving a spoon in the coffee cup OR taking the spoon out.
e) Leaving the coffee near an open window where there is a slight draft OR where there is no wind.
Solution:
a) Pouring the cold milk in as soon as the coffee has been made will keep it warmer.
Pouring cold milk into the coffee after the coffee has been left to cool for some time results in a larger
temperature gradient and the coffee will cool down faster at the start.
b) Covering the top of the coffee cup will keep it warmer.
This will reduce convective heat transfer between surface of the coffee and the surrounding air.
c) Using a ceramic mug will keep it warmer.
Ceramic is a good insulator, while metal is a very good conductor. Convective heat loss from the coffee
through the walls of the mug to the surrounding air will be faster through metal walls than ceramic
walls.
d) Taking the spoon out of the coffee mug will keep it warmer.
The spoon (assuming it is metal) is a good conductor and acts as a type of fin and increases the surface
area for heat transfer. This would increase the convective heat loss.
e) Leaving the coffee cup where there is no wind will keep it warmer.
A draft around the outside of the coffee mug will increase the rate of convective heat loss from the
walls of the mug and from the liquid surface. The draft moves heat away from the outer edges of the
mug and the liquid surface faster; thereby creating a larger temperature gradient. Heat will then be
transferred from the system more quickly.
∆$
q0 ∝ 3.1
∆)
where:
q’ = Heat flux (W/m2)
DT = Temperature difference (K)
x = direction in which there is a temperature gradient (direction of heat transfer) (m)
where:
k = thermal conductivity (W/(m.K))
Or more conveniently:
D$
q̇ = −k. A D) 3.3
where:
q̇ = Heat (W)
A = Cross sectional area (m2)
Thot
Tcold
q̇
Assumptions:
- Thot is at the same temperature across entire area;
- Use area perpendicular to the direction of heat transfer; and
- x is the direction from hot to cold (direction of temperature gradient).
Th
Tc
k
q̇
L
x=0 x=L
where:
Assumptions:
- Area through slab is constant
- k is constant
- Steady state conditions. “Energy in one face of the slab” = “Energy out other face of
the slab”
Therefore:
)IO K̇ $
∫)I/ dx = ∫$ D dT 3.5
#MN E
Since none of q̇ , k or A are functions of temperature or x, we can remove these from the
integral:
K̇ )IO $
∫ dx = ∫$ D dT 3.6
#MN )I/ E
And integrate:
K̇
(L − 0) = (TE − T4 ) 3.7
#MN
Therefore:
MN($E #$D )
q̇ = 3.8
O
An aluminium plate (k = 215 W/m°C) is heated to 300°C. If the heat flux is 8.6 MW/m2, how hot is the
other face, if the metal is 5 mm thick?
Solution:
dT
𝑞1 = −k
dx
Assuming constant k and A:
𝑘(𝑇F − 𝑇G )
𝑞1 =
𝐿
𝑞1 𝐿
∴ = 𝑇F − 𝑇G
𝑘
𝑞1 𝐿 (8600000)(0.005)
∴ 𝑇G = 𝑇F − = 300 −
𝑘 215
∴ 𝑇G = 100
T1
Th
Tc
k1 k2
q̇
L1 L2
x=0 x = L1 x = L1 + L2
where:
Th: Hot temperature on one side of the slab (K)
T1: Temperature at the end of slab 1 (K)
Tc: Cold temperature on the other side of slabs (K)
q̇ : Heat Flow (W)
A: Area perpendicular to heat flow (m2)
L1: Length of slab 1 (m)
L2: Length of slab 2 (m)
x: Direction in which heat flows (m)
k1: Thermal conductivity of slab 1 (W/(m.K))
k2: Thermal conductivity of the slab 2 (W/(m.K))
Assumptions:
- Area through slab is constant
- k1 ≠ k2
- Steady state conditions; therefore, “Energy in one face of the slab” = “Energy out other face
of the slab”
Therefore:
M" N($E #$" )
Slab 1: q̇ = 3.10
O"
M! N($" #$D )
Slab 2: q̇ = 3.11
O!
Adding these equations, we are left with an expression without T1 and can calculate q̇ :
K̇ O" K̇ O!
T4 − T& + T& − TE = + 3.14
M" N M! N
K̇ O" K̇ O!
T4 − T& + T& − TE = + 3.15
M" N M! N
K̇ O" O
T4 − TE = % + M! & 3.16
N M" !
K̇ $E #$D
= H H 3.17
N P "- !Q
I" I!
To determine T1, either of the original heat transfer equations can then be used.
The wall of a furnace is constructed from an inner steel layer of 0.5 cm (k = 40 W/m.K) and a brick outer layer of
10 cm (k = 2.5 W/mk). The inner surface temperature is 900K and the outside surface temperature is 460K. What is
the temperature between the steel and the brick?
Solution:
Assumptions:
- Steady State (energy in = energy out)
- Linear heat transfer (ignore the effects of edges)
&'
From Fourier’s Law: 𝑞 = −𝑘𝐴 &!
( '!" "'#$%
For 2 slabs this simplifies to: )
= &' &)
*
(' ()
Therefore:
𝑞 (900 − 460) 440
= = = 10965.732 … 𝑊/𝑚$
𝐴 0.005 0.1 (0.000125 + 0.04)
40 +
2.5
If we let 𝑇# be the temperature between the steel and the brick and:
-./'
𝑇# = 𝑇+, − = 900 − (10965)(0.000125) = 898.629 … 𝐾
).0'
T1 T2 T3 Tn-1
Th
q̇
L1 L2 L3 Ln
x=0 x = L1 x = L1 + L2 + L3
x = L1 + L2 x = L1 + L2 + L3 + … + Ln
where:
Th: Hot temperature on one side of the slab (K)
Ti: Temperature at the end of slab i (K)
Tc: Cold temperature on the other side of slabs (K)
q̇ : Heat Flow (W)
A: Area perpendicular to heat flow (m2)
Li: Length of slab i (m)
x: Direction in which heat flows (m)
ki: Thermal conductivity of slab i (W/(m.K))
n: Number of slabs
Assumptions:
- Area through slab is constant
Kevin Harding (kevin.harding@wits.ac.za) 57
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Conductive heat transfer
- ki ≠ ki +1
- Steady state conditions; therefore “Energy in one face of the slab” = “Energy out other face
of the slab”
As before we can solve for the heat flow through each slab:
M" N($E #$" )
Slab 1: q̇ = 3.18
O"
M! N($" #$! )
Slab 2: q̇ = 3.19
O!
M# N($! #$# )
Slab 3: q̇ = 3.20
O#
M* N($*J" #$D )
Slab n: q̇ = 3.21
O*
Rearranging as before:
K̇ O"
T4 − T& = 3.22
M" N
K̇ O!
T& − T. = 3.23
M! N
K̇ O#
T. − T% = 3.24
M# N
K̇ O*
T,#& − TE = 3.25
M* N
Adding these equations, we are left with an expression without Ti and can calculate q̇ :
K̇ O O O O*
T4 − TE = % " + M! + M# + ⋯ + & 3.26
N M " ! # M*
On a chilly winter day, you decide to add on a few layers of clothes to stay warm. What heat flux do you
experience for each successive layer added as below? Assume your skin is at 36°C and the temperature
outside is at 4°C:
1) T-Shirt (k = 0.05 W/mK; Thickness = 0.75 mm)
2) Jersey 1 (k = 0.06 W/mK; Thickness = 2 mm)
3) Jersey 2 (k = 0.05 W/mK; Thickness = 2 mm)
4) Jacket (k = 0.005 W/mK; Thickness = 3 mm)
5) Overcoat (k = 0.02 W/mK; Thickness = 4.5 mm)
It is estimated that the body released between 90-140 W of heat. Further, the body has an average surface
area of between 1.6 and 1.9 m2. Given this information, what would be a reasonable number of layers to
wear?
Solution:
∆O ∆O
q̇ = \] OR 𝑞1 = \]
∑`
]a" ∑`
]a"
^] _ ^]
∆O ,(
1) 𝑞 = 1
\" = b.cd×"bef
= 2133.3 𝑊/𝑚(
^" b.bd
∆O ,(
2) 𝑞1 = \" \ = b.cd×"bef !×"bef
= 662.1 𝑊/𝑚(
0 ! 0
^" ^! b.bd b.bg
∆O ,(
3) 𝑞 = \" 1
\ \ = b.cd×"bef !×"bef !×"bef
= 362.2 𝑊/𝑚(
0 !0 f 0 0
^" ^! ^f b.bd b.bg b.bd
∆O ,(
4) 𝑞1 = \" \ \ \ = b.cd×"bef !×"bef !×"bef f×"bef
= 46.5 𝑊/𝑚(
0 !0 f0 h 0 0 0
^" ^! ^f ^h b.bd b.bg b.bd b.bbd
∆O ,(
5) 𝑞1 = \" \ \ \ \ = b.cd×"bef !×"bef !×"bef f×"bef h.d×"bef
= 35 𝑊/𝑚(
0 !0 f0 h0 d
^" ^! ^f ^h ^d 0 0 0 0
b.bd b.bg b.bd b.bbd b.b!
From the average body data, it can be shown that the body releases between 47 and 87 W/m2. Therefore,
is would be a reasonable assumption that the average person would be most comfortable with four layers
as calculated.
Rearranging:
q̇ L = kA(T4 − TE ) 3.29
K̇ O
= T4 − TE 3.30
MN
#K̇ O
TE = + T4 3.31
MN
LARGE, and the slope is great, i.e., not much temperature change. NOTE: Only applicable
for a constant heat transfer constant.
k
Th
Tc
x
0 L
In the same way, this can be extended for n-slabs in series. The temperature profile in each
slab is linear (provided k is constant), with a different linear slope depending on the value
of k.
k1 k2 k3 ... kn
Th
Tc
x
x1 x2 x3 xn-1 xn
0
Tin
ri
Tout
ro
q̇
where:
Tin: Temperature on the inside of the cylinder (K)
Tout: Temperature on the outside of the cylinder (K)
q̇ : Heat Flow (W)
ri: Inner radius of the cylinder (m)
ro: Outer radius of the cylinder (m)
k: Thermal conductivity of the slab (W/(m.K))
Assumptions:
- Cylinder thickness is constant throughout
- k is constant
- Steady state conditions; therefore, “Energy in one face of the slab” = “Energy out other face
of the slab”
- No axial heat flow
- Area through which heat flows is NOT constant!
But we do not have an x-direction in a cylinder. Therefore, convert Fourier’s Law into
cylindrical equivalent:
D$
q̇ = −k. A DC 3.32
The area is that which is perpendicular to the direction of heat transfer. Therefore:
A = 2π.r.L (circumference of a cylinder × length) 3.33
Kevin Harding (kevin.harding@wits.ac.za) 61
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Conductive heat transfer
Integrating:
K̇ C
ln % CQ & = T>?@ − T=, 3.36
#M.TO K
Or:
K̇ C
T=, − T>?@ = M.TO ln % CQ & 3.37
K
Water enters a pipe at a temperature of 60°C. Assuming the pipe is made of stainless steel (k = 30 W/m2)
and that the pipe has an inner and outer diameter of 20 and 25 cm respectively, what is the outer
temperature of the pipe if heat is lost at a rate of 15 kJ/m.s? Assume the temperature is constant down the
length of the pipe.
Solution:
𝑞̇ 𝑟7
𝑇5/ − 𝑇789 = 𝑙𝑛 0 7
𝑘2𝜋𝐿 𝑟5
𝑞̇ 𝑟7 15000 25
𝑇789 = 𝑇5/ − 𝑙𝑛 0 7 = 60 − 𝑙𝑛 0 7
𝑘2𝜋𝐿 𝑟5 (30)2𝜋 20
𝑇789 = 42.2°𝐶
Liquid ammonia at 0.5 atm (-45°C) is stored in a spherical container of ID 1 foot and OD 18 inches. The
walls of the sphere are filled with insulation (k = 0.0692 W.m-1.K-1), and the outer surface is at room
temperature of 20°C. What is the heat transfer through the sphere? Assuming the temperature in the
sphere remains constant, how much ammonia is lost per hour by evaporation? Hvap NH3: 1390 kJ/kg
Solution:
Fourier’s law:
𝑑𝑇
𝑞̇ = −𝑘𝐴
𝑑𝑟
𝐴 = 𝑠𝑢𝑟𝑓𝑎𝑐𝑒 𝑎𝑟𝑒𝑎 𝑜𝑓 𝑠𝑝ℎ𝑒𝑟𝑒 = 4𝜋𝑟 (
𝑞̇ 𝑑𝑟
∴ = 𝑑𝑇
−4𝜋𝑘 𝑟 (
Ti Ui
𝑞̇ 𝑑𝑟
Q ( = Q 𝑑𝑇
−4𝜋𝑘 𝑟
Tj Uj
𝑞̇ −1 Ti
∴ ? @ = 𝑇5 − 𝑇7
−4𝜋𝑘 𝑟 Tj
𝑞̇ 1 1
0 − 7 = 𝑇5 − 𝑇7
4𝜋𝑘 𝑟5 𝑟7
4𝜋𝑘(𝑇5 − 𝑇7 )
𝑞̇ =
1 1
)𝑟 − 𝑟 *
5 7
4𝜋(0.0692𝑊𝑚%& 𝐾 %& )(228𝐾 − 293𝐾)
𝑞̇ = = −25.84 𝑊 = −25.84 𝐽𝑠 %& = −93.024 𝑘𝐽ℎ%&
1 1
−
)0.1524𝑚 0.2286𝑚*
𝑞̇ 93.024 𝑘𝐽ℎ%&
𝑚̇ = = = 0.06692 𝑘𝑔ℎ%& = 66.92 𝑔ℎ%&
∆𝐻mVWX,YZ, 1390 𝑘𝐽𝑘𝑔%&
Now:
K̇ )IO & $
∫ dx = ∫$ D dT 3.38
#M )I/ N E
In a building in a hot climate (45°C), the inside is cooled to 15°C. The walls of the building include an
irregular shaped steel girder which forms a structural member of the wall. If no heat is lost or gained to
or from the insulation through the sides of the girder, calculate the heat flow from the room from a single
girder.
Solution:
𝑑𝑇
𝑞̇ = −𝑘. 𝐴
𝑑𝑥
𝑑𝑇
𝑞̇ = −𝑘. (−50𝑥 ( + 3)
𝑑𝑥
𝑞̇
𝑑𝑥 = 𝑑𝑇
−𝑘(−50𝑥 ( + 3)
$6[ Uk
𝑞̇
Q 𝑑𝑥 = Q 𝑑𝑇
−𝑘(−50𝑥 ( + 3)
$6\ Ul
$6[ Uk
𝑞̇ 1
Q 𝑑𝑥 = Q 𝑑𝑇
𝑘 50𝑥 ( − 3
$6\ Ul
[
𝑞̇ 2 2 × 3𝑥
? 𝑡𝑎𝑛%& @ = 𝑇F − 𝑇G
𝑘 √4 × 50 × 3 √4 × 50 × 3 \
𝑞̇ 2 6𝑥 \.(
? 𝑡𝑎𝑛%& @ = 45 − 15 = 30
45 √600 √600 \
(30 × 45)
𝑞̇ =
2 6(0.2) 2 6(0)
0 𝑡𝑎𝑛%& 7− 0 𝑡𝑎𝑛%& 7
√600 √600 √600 √600
1350
𝑞̇ =
1.2
)0.081 𝑡𝑎𝑛%& 24.5* − 0
1350
𝑞̇ = = 810.6 𝑊
1.665
The figure shows the plan view of a steel girder forming a structural member in the wall of a cold room.
The thickness of the wall is 30 cm. The width of the girder is 10 cm on the inside and 20 cm on the outside.
Assuming no heat transfer across the insulation on the sides, calculate the heat flow from the cold room
To = 25°C
Ti = 5°C
per metre of the girder. k = 45W.m-1.K-1
Solution:
Drawing the system in 3-D, with heat flowing perpendicular to the diagram, the surface area perpendicular
to the flow of heat will be rectangular, with the area increasing as the heat flows from the cold room to
the outside because of the tapering sides of the girder (the girder boundaries become increasingly further
apart).
20 cm To = 25°C
30 cm
10 cm x
Ti = 5°C
Area
We need an expression to describe how the rectangular surface area perpendicular to the direction of heat
flow changes with distance through the wall. The sides of the girder can be described using the equation
for a straight line. Defining an x-y coordinate system such that at x = 0, y = 10 cm , and at x = 30 cm, y =
20 cm.
Let u = 1/3 x + 10
du = 1/3 dx
3du = dx
And at x = 0, u = 0.10
at x = 0.30, u = 0.20
\.(\ Uj
𝑞̇ 3
∴ Q 𝑑𝑢 = Q 𝑑𝑇
−𝑘 𝑢
\.&\ Ui
3𝑞̇ 0.20
. 𝑙𝑛 0 7 = 𝑇] − 𝑇5
−𝑘 0.10
−𝑘(𝑇] − 𝑇5 ) −45𝑊𝑚%& 𝐾 %& (298𝐾 − 278𝐾)
𝑞̇ = = = −433 𝑊. 𝑚%&
3 × 𝑙𝑛 (2) 3 × 𝑙𝑛 (2)
The sign is negative because heat flow is conventionally down the temperature gradient (i.e. from hot to
cold)
)IO K̇ $
∫)I/ dx = ∫$ D dT 3.5
#MN E
Now:
K̇ )IO $
− N ∫)I/ dx = ∫$ D kdT 3.39
E
Example 3.9: Conduction with variable conductive heat transfer coefficient (1)
A furnace wall has a thickness of 0.5 ft, with inner and outer temperatures of 200°F and 100°F respectively.
The thermal conductivity is given by: K(BTU/h.ft.°F) = 0.1 + 0.001T(°F).
Calculate the heat transfer per square foot through the furnace wall.
Solution:
𝑑𝑇
𝑞̇ = −𝑘𝐴
𝑑𝑥
$ U!
𝑞̇
Q 𝑑𝑥 = − Q 𝑘𝑑𝑇
\ 𝐴 U"
𝑞̇ $ U!
Q 𝑑𝑥 = − Q (0.1 + 0.001 𝑇)𝑑𝑇
𝐴 \ U"
𝑞̇ 0.001𝑇( ( 0.001𝑇& (
(𝑥 − 0) = − r(0.1𝑇( + , − (0.1𝑇& + ,s
𝐴 2 2
𝑞̇ 0.001(100)( 0.001(200)(
(0.5) = − r(0.1(100) + , − (0.1(200) + ,s
𝐴 2 2
𝑞̇ t(10) − (40)u
=−
𝐴 0.5
^̇
`
= 60 BTU/hr.ft2
Example 3.10: Conduction with variable conductive heat transfer coefficient (2)
A house in a cold climate with a flat insulated rook is heated and maintained at 23°C. The roof consists of
three layers that are in contact with one another: ceiling board (1.50 cm thick), insulation (40.0 cm thick)
and a waterproof layer on the outside (0.5 cm thick). The thermal conductivities of the ceiling board,
insulation and waterproof material are k1 = 0.690 W.m-1.K-1, k2 = 85.0 W.m-1.K-1 and k3 = 377 W.m-1.K-1.
Lying on the waterproof material is a layer of snow, which gets thicker as the snow falls. Due to
compression of the snow in contact with the roof, the thermal conductivity is a function of the thickness
of the snow according to:
1
𝑘a/7b =
𝑎𝑙 + 𝑏
Where:
l is the thickness of the snow in m, with a = 15.6 K.W-1 and b = 1.24 m.K.W-1
The film heat transfer coefficient on the inside is 18 W.m-2.K-1 and on the outside it is 50 W.m-2.K-1. If
heat is lost at a rate of 50.0 W.m-2 and the outside temperature is -5°C, what is the thickness of the snow
on the roof and the overall heat transfer coefficient?
Ta
h0 T4
L4
L (m)
L3 = 42cm
3
L2 = 41.5cm
2
1 L1 = 1.5cm
0cm
hi
Solution:
𝑞̇v %(
𝐴 = 50.0 𝑊. 𝑚
𝑇W = −5℃ = 268𝐾
For the inner film:
𝑞̇ 𝑞̇
= ℎ5 (𝑇5 − 𝑇\ ) ⇒ 𝑇\ = 𝑇5 −
𝐴 𝐴ℎ5
∴ 𝑇\ = 293.22𝐾
For the ceiling board:
𝑞̇ 𝑑𝑇 𝑞̇
= −𝑘& ⇒ (𝐿 ) = 𝑇& − 𝑇\
𝐴 𝑑𝐿 −𝑘& 𝐴 &
∴ 𝑇& = 292.13𝐾
For the insulation:
𝑞̇ 𝑑𝑇 𝑞̇
= −𝑘( ⇒ (𝐿 − 𝐿& ) = 𝑇( − 𝑇&
𝐴 𝑑𝐿 −𝑘( 𝐴 (
∴ 𝑇( = 291.89𝐾
For the waterproofing:
𝑞̇ 𝑑𝑇 𝑞̇
= −𝑘, ⇒ (𝐿 − 𝐿( ) = 𝑇, − 𝑇(
𝐴 𝑑𝐿 −𝑘, 𝐴 ,
∴ 𝑇, = 291.89𝐾
For the snow:
𝑞̇ 1 𝑑𝑇
=−
𝐴 (𝑎𝐿 + 𝑏) 𝑑𝐿
𝑞̇
⇒ − (𝑎𝐿 + 𝑏)𝑑𝐿 = 𝑑𝑇
𝐴
Integrating, with L = thickness of the snow:
^̇ W[!
−`) (
+ 𝑏𝐿* = 𝑇- − 𝑇, (1)
For the outer film:
𝑞̇ 𝑞̇
= ℎ\ (𝑇- − 𝑇W ) ⇒ 𝑇- = 𝑇W +
𝐴 𝐴ℎ\
∴ 𝑇- = 269𝐾
From (1):
15.6
−50𝑊. 𝑚%( /0 𝐾𝑊 %& 7 𝐿( + (1.24 𝑚𝐾𝑊 %& )𝐿= = 269𝐾 − 291.89𝐾
2
(7.8 𝐾𝑊 %& )𝐿( + (1.24 𝑚𝐾𝑊 %& )𝐿 − 0.4578 𝑚𝐾. 𝑊 %& = 0
Or:
$E #$D
q̇ = H" H! H H 3.41
P - - # -⋯- * Q
I" L I! L I# L I* L
∆$
q̇ = HK 3.42
∑*
KM" I L
K
where:
DT: Temperature difference (K)
R: Thermal resistance (K/W)
OK
For slab calculations, R = ∑,=I& , however, R can represent the resistance for any range
MK N
of shapes.
R1 R2 R3
A circular pipe of 20 cm is enclosed centrally in a square section insulator of side 36 cm. The thermal
conductivity of the material is given as 8.5W/mK. The inside surface is at 200°C, while the outside is at
30°C. Determine the heat flow for a length of 5 m.
where:
a = side length of length of square section
l = length
r = radius of inner portion
k = 8.5 W/m.K
30°C
0.2 m
b=5m
200°C
a = 0.36m
Solution:
∆𝑇
𝑞̇ =
𝑅
1 1 1.08 × 0.36
𝑅= } 𝑙𝑛 ~ = 0.002489362 °C/W
2𝜋 × 5 8.5 0.2
2. 2
∆𝑇 200 − 30
∴ 𝑞̇ = = = 68290.59 = 68 000 W
𝑅 0.002489 …
A pipe of 30 cm diameter is insulated by a material of thermal conductivity 0.45 W/mK. Due to space
restrictions, the pipe is placed slightly off-centre (5cm), resulting in a portion of insulation thicker than
the rest of the pipe (see diagram). The inner surface is at 250°C and the outer surface is at 60°C. Determine
the heat loss for a 5 m length of pipe.
where:
e = eccentricity, m
60°C
0.15 m
0.05 m
0.3 m
250°C
k = 0.45 W/m.K
Solution:
∆𝑇
𝑞̇ =
𝑅
1 €[(0.15 + 0.3)( − 0.05( ] + €[(0.15 − 0.3)( − 0.05( ]
𝑅= 𝑙𝑛
2𝜋 × 0.45 × 5 €[(0.15 + 0.3)( − 0.05( ] − €[(0.15 − 0.3)( − 0.05( ]
1 0.4472 + 0.14142
𝑅= 𝑙𝑛 = 0.046325 ℃/W
2𝜋 × 0.45 × 5 0.4472 − 0.14142
250 − 60
∴ 𝑞̇ = = 4101 W
0.046325
Ts
Tb
Solid surface
where:
δ: thickness of stationary boundary layer (m)
Ts: Surface temperature (K)
Tb: Bulk temperature (K)
If we assume that there is no movement in a boundary layer between the slab and bulk
fluid (δ – Greek delta), we can approximate the heat transfer to Fourier’s Law:
MN($k #$l )
q̇ = 3.44
W
where:
h: Heat Transfer Coefficient (W/m2.K)
Note:
It is harder to measure wall temperatures compared to the bulk fluid temperatures.
Therefore, the equations are developed to use the bulk fluid temperatures and not
intermediate wall temperatures.
T1
Tc
q̇
hc
where:
T1: Temperature at the end of slab 1 (K)
T2: Temperature at the end of slab 2 (K)
Tc: Cold temperature (K)
q̇ : Heat Flow (W)
A: Area perpendicular to heat flow (m2)
hc: Heat Transfer coefficient on the cold side (W/m2.K)
Assumptions:
- The fluid has a heat transfer coefficient (h)
- Area through convective area is constant
- Steady state conditions, therefore, “Energy in” = “Energy out”
A stainless-steel plate is maintained at 150°C, while the surrounding temperature is 20°C. The heat transfer
coefficient is given as 11.2 W/m2.K. What is the convective heat loss per square meter?
Solution:
𝑞̇ = ℎ𝐴(𝑇( − 𝑇& )
𝑞̇
= ℎ(𝑇( − 𝑇& )
𝐴
𝑞̇
= 11.2(150 − 20) = 1456 𝑊/𝑚( = 1.456 𝑘𝑊/𝑚(
𝐴
T0 T1 T2
Th
Tc
q̇
hh k1 k2 hc
L1 L2
x=0 x = L1 x = L1 + L2
where:
Th: Hot temperature on one side of the slab (K)
T0: Temperature at the start of slab 1 (K)
T1: Temperature at the end of slab 1 (K)
T2: Temperature at the end of slab 2 (K)
Tc: Cold temperature on the other side of slabs (K)
q̇ : Heat Flow (W)
A: Area perpendicular to heat flow (m2)
L1: Length of slab 1 (m)
L2: Length of slab 2 (m)
x: Direction in which heat flows (m)
k1: Thermal conductivity of slab 1 (W/(m.K))
k2: Thermal conductivity of the slab 2 (W/(m.K))
hh: Heat Transfer coefficient on the hot side (W/m2.K)
hc: Heat Transfer coefficient on the cold side (W/m2.K)
Assumptions:
- The slabs have thermal conductivities (k); the fluids have heat transfer coefficients
(h)
Then adding the last four equations, we are left with an expression without Tc and are able
to calculate q̇ :
K̇ K̇ O" K̇ O! K̇
T4 − T/ + T/ − T& + T& − T. + T. −TE = +M +M +4 3.57
4E N "N !N DN
K̇ K̇ O" K̇ O! K̇
T4 − T/ + T/ − T& + T& − T. + T. −TE = +M +M +4 3.58
4E N " N ! N DN
K̇ & O O &
T4 −TE = % + M" + M! + 4 & 3.59
N 4E " ! D
K̇ $E #$D
= " H" H! "
3.60
N Y - - - Z
EE I" I! ED
Given the area, heat transfer coefficients ad thermal conductivities, only need to know
inner and outer temperatures (NOT any of the inside temperatures) to solve for q.
The inner surface of a 23 cm brick furnace wall is kept at 820°C, while the temperature of the outer surface
is 180°C.
a) Calculate the heat loss from the furnace, given a thermal conductivity of the brick of 0.865 W.m-1.K-1.
b) A further 23 cm is added to the furnace wall, consisting of insulating brick with a thermal conductivity
of 0.26W.m-1.K-1. Calculate the reduction in the heat loss, and the maximum temperature to which
the insulating brick is exposed. Assume that the inner temperature remains at 820°C. The temperature
of the surrounding air is 20°C.
c) Re-calculate the heat loss and the maximum temperature to which the insulating is exposed, but allow
for a convective heat transfer coefficient on the outer surface of 11.9 W.m-2.K-1. Assume that the inner
temperature remains at 820°C, and the temperature of the surrounding air is 20°C.
d) How do the answers compare when convection is included?
Solution:
mU
a) 𝑞̇ = −𝑘𝐴 m$
[ Uj
𝑞̇
∴ Q 𝑑𝑥 = Q 𝑑𝑇
−𝑘𝐴
\ Ui
𝑞̇ 𝑘 0.865 𝑊𝑚%& 𝐾 %&
= (𝑇5 − 𝑇\ ) = (1093𝐾 − 453𝐾) = 2407 𝑊. 𝑚%(
𝐴 𝐿 0.23𝑚
b) For wall 1:
𝑑𝑇
𝑞̇ = −𝑘& 𝐴
𝑑𝑥
[" Um
𝑞̇
∴ Q 𝑑𝑥 = Q 𝑑𝑇
−𝑘& 𝐴
\ Ui
𝑞̇ 𝐿&
⇒ = 𝑇5 − 𝑇b
𝑘& 𝐴
^̇ e
`
= [" (𝑇5 − 𝑇b ) (1)
"
For wall 2:
𝑑𝑇
𝑞̇ = −𝑘( 𝐴
𝑑𝑥
^̇ e!
`
= [
(𝑇b − 𝑇\ ) (2)
!
From (2):
𝑞̇ 0.26 𝑊. 𝑚%& . 𝐾 %&
= (908 𝐾 − 293𝐾) = 695.2 𝑊. 𝑚%(
𝐴 0.23 𝑚
∆𝑞 = 2407 𝑊. 𝑚%( − 695 𝑊. 𝑚%( = 1712 𝑊. 𝑚%( reduction in heat loss
c) For wall 1:
𝑑𝑇
𝑞̇ = −𝑘& 𝐴
𝑑𝑥
^̇ e"
`
= [ (𝑇5 − 𝑇b ) (1)
"
For wall 2:
^̇ e!
`
=
[!
t𝑇b − 𝑇4 u (2)
For the film:
𝑞̇ = ℎ𝐴(𝑇4 − 𝑇\ )
^̇
𝑇4 − 𝑇\ =
F`
(3)
Add (1) + (2):
^̇ [ [
𝑇5 − 𝑇4 = ` )e" + e! * (4)
" !
d) The heat loss was 695.2 W/m2 neglecting the convective film heat transfer on the outer surface of the
furnace wall. Taking the convective heat transfer into account (as calculated above), there is a 47.2
W/m2 reduction in heat loss. The film acts as an additional insulating layer.
Reminder:
∆$
q̇ = 3.61
V
where:
Ri : individual thermal resistances for different layers (K/W)
Example 3.15:
Solution
T0 T1 T2 T3 Tn-1 Tn
Th
Tc
…
q̇
hh k1 k2 k3 kn hc
L1 L2 L3 L3
x=0 x = L1 x = L1 + L2 + L3
x = L1 + L2 x = L1 + L2 + L3 + … + Ln
Example 3.16:
Solution
Insulation
Pipe wall
Tout T2 T1 T0 Tin
r2
hin r0
k1
k2 r1
hout
q̇
where:
Tin: Temperature on the inside of the cylinder (K)
Tout: Temperature on the outside of the cylinder (K)
q̇ : Heat Flow (W)
r0: Inner radius of the cylinder (m)
r1: Outer radius of the 1st cylinder (pipe) (m)
r2: Outer radius of the 2nd cylinder (insulation) (m)
k1: Thermal conductivity of cylinder 1 (W/(m.K))
k2: Thermal conductivity of cylinder 2 (W/(m.K))
hin: Heat Transfer coefficient on the inside (inside) (W/m2.K)
hout: Heat Transfer coefficient on the outside (outside) (W/m2.K)
Assumptions:
Kevin Harding (kevin.harding@wits.ac.za) 78
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Convective heat transfer
Area is the area in contact with the film (we don’t know how thick the film is so can’t use
any other area)
∴ q̇ = h=, 2πr/ L(T=, − T/ ) 3.67
K̇
T=, − T/ = 3.68
4K .TC. O
Integrating:
K̇ C
T/ − T& = M ln %C" & 3.71
" .TO .
K̇
T=, − T/ = 3.75
4K* .TC. O
K̇ C"
T/ − T& = M ln %C & 3.76
" .TO .
K̇ C
T& − T. = M ln %C! & 3.77
! .TO "
K̇
T. − T>?@ = 3.78
4QRS .TC! O
Saturated steam at 14 bar (195°C) condenses inside a nickel tube 10 m long with a 2 cm ID and 4 cm OD.
The temperature of the surrounding air is 25°C. The thermal conductivity of nickel is 140 W.m-1.K-1. The
pipe is covered with a 0.5 cm layer of insulation (thermal conductivity 0.21 W.m-1.K-1).
a) Calculate the conductive energy across the pipe and insulation.
b) Given a convective heat transfer coefficient on the outside of the insulation (10 W.m-2.K-1), determine
the overall heat transfer coefficient, Uo
c) Determine the mass of steam which condenses per hour.
d) How do answers compare when convection is included vs when it is not?
Solution:
Ta T3 T2 T1 r3
r1 r2
h
mU
a) 𝑞̇ = −𝑘𝐴 mT
𝐴 = 2𝜋𝑟𝐿
𝑞̇ 𝑑𝑟
∴ = 𝑑𝑇
−2𝜋𝑘𝐿 𝑟
The conductive energy across the pipe and insulation is approx. 10 kW.
For this system, the overall heat transfer co-efficient (U0) is given by:
1
𝑈\ =
𝑙𝑛 (𝑟( ⁄𝑟& ) 𝑙𝑛 (𝑟, ⁄𝑟( ) 1
0 2𝜋𝑘 𝐿 + 2𝜋𝑘 𝐿 + 2𝜋ℎ𝑟 𝐿7
Y q ,
𝑈\
1
=
𝑙𝑛 (2𝑐𝑚⁄1𝑐𝑚) 𝑙𝑛 (2.5𝑐𝑚⁄2𝑐𝑚) 1
0 + + 7
2𝜋 × 140 𝑊𝑚%& 𝐾 %& × 10𝑚 2𝜋 × 0.21 𝑊𝑚%& 𝐾 %& × 10𝑚 2𝜋 × 10𝑊𝑚%( 𝐾 %& × 0.025𝑚 × 10𝑚
& %&
𝑈\ = e" = 12.4 𝑊𝐾
\.\f\rs t u
For saturated steam at 14bar, Ts = 195°C, Hf = 830 kJ/kg and Hg = 2790 kJ/kg
∴ ∆𝐻VWX = 1960 𝑘𝐽/𝑘𝑔
The enthalpy change is the heat of vaporization since the steam is condensing.
2.108 𝑘𝐽𝑠 %&
𝑚̇ = = 1.076 × 10%, 𝑘𝑔𝑠 %& = 3.872 𝑘𝑔ℎ%&
1960 𝑘𝐽𝑘𝑔%&
Kevin Harding (kevin.harding@wits.ac.za) 81
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Convective heat transfer
d) The heat transfer is smaller when convection is included (2.108 kW compared to 10 kW)
where:
rLM = Log mean radius
r2 = outer radius
r1 = inner radius
Note:
C! # C" C! # C"
p ≠ 3.84
[, ! .
p"
where:
ALM = Log mean area
From the heat flow through two cylinders, (conduction and convection, it was shown that:
$K* #$QRS
q̇ = " " p" " p "
3.86
Y - [,P Q- [,P ! Q- Z
EK* !op. H I" !oH p. I! !oH p" EQRS !op! H
From the definition of RLM, and converting wherever possible to surface area instead of
radius, this becomes:
$K* #$QRS
q̇ = " " (p" Jp. ) " (p Jp ) "
3.88
Y - . - . ! "- Z
EK* L. I" !oH pHv,". I! !oH pHv,!" EQRS L!
$K* #$QRS
q̇ = " (p" Jp. ) (p Jp ) "
3.89
Y - - ! " - Z
EK* L. w" LHv,". w! LHv,!" EQRS L!
$K* # $QRS
q̇ = " * ∆xK "
3.90
Y -∑KM" - Z
EK* L. IK LHv,K EQRS L*
where:
DT: Temperature difference (K)
R: Thermal resistance (K/W)
For multiple slabs in series the thermal resistance was shown to be:
OK
R = ∑,=I& 3.92
MK N
A new term is now introduced for the Overall Heat Transfer Coefficient (U):
&
U = N.V 3.95
∴ q̇ = UA(T=, − T>?@ ) 3.96
Note: Thermal resistance (R) includes the term for area (A), while the overall heat transfer
coefficient (A) does not. Further, R and U are inverse relationships.
&
UA = " ∆pK "
3.98
Y -∑ *
KM" I - Z
EK* K EQRS
While we typically use the ALM for these calculations, or assume the thickness is thin
enough not to worry about which radius, and thus aera, to use, we can then define different
U values for different areas:
]N
U/ = For outer cylinder 3.100
N.
]N
U= = For inner cylinder 3.101
NK
Example 3.18: Convective heat transfer, with conduction, across multiple layers (overall heat transfer
coefficient)
A furnace is constructed of a 130 mm steel layer and an outer surface of bricks which is 170 mm thick.
The thermal conductivities of the steel and brick layers are 144 W/m.K and 4.38 W/m.K respectively,
whilst the heat transfer of the inner and outer surface are 18W/𝑚( . 𝐾 and 73W/𝑚( . 𝐾. Given that the
temperature inside the furnace is kept at 650°C and the temperature of the surroundings is 45°C, calculate
the overall heat flow and overall heat transfer coefficient.
Solution:
Because the system is at steady state, with areas and thermal conductivities, it was shown earlier that:
𝑞 𝑇5 − 𝑇W
=
𝐴 1 + 𝐿& + 𝐿( + 1
ℎ5 𝑘& 𝑘( ℎ\
Therefore:
𝑞
= 5552𝑊/𝑚(
𝐴
Now: 𝑞 = 𝑈𝐴(𝑇5 − 𝑇W )
^
𝑈 = `(U %U ) = 9.18 W/𝑚( . 𝐾
i p
Th
Tc
k1 k2 k3
q̇ hh hc
L1 L2 L3
Rh R1 R2 R3 Rhc
h
where:
Ri = Thermal resistance through each layer (including convection and conduction)
In the same way as in an electrical circuit, to get the total resistance in series, add the
resistance.
where:
R = Total thermal resistance
Therefore, from:
∆$
q̇ = 3.91
V
∆$
q̇ = V 3.103
EE -V" -V! -V# -VED
∆$
q̇ = " H H H " 3.104
- "- !- #-
EEE. L I" L I! L I# L EED. L
The sum of resistance in a parallel configuration is equal to the sum of the inverse
resistances, all to the power of negative one.
R2
where:
R = Total thermal resistance
Given the configuration below, heat will travel from the left to the right through the
parallel configuration. In the calculations we have performed so far, we have assumed 1-
dimendionsal heat transfer. However, with the parallel configuration now in place, when
k2 and k3 are not equal, heat will move between these sections. To ensure the 1-
dimensional assumptions in the calculations, two scenarios are presented with simplifying
assumptions.
L2
Th
Tc
k1 k2 k4
q̇ hh hc
k3
L1 L3 L4
Scenario 1: Assuming the surfaces normal to the direction of heat flow are isothermal, i.e.,
the entire length of slabs 1 and 4 have the same temperature, this gives:
Rh R1 R2 R4 Rhc
h
R3
where:
Ri = Thermal resistance through each layer (including convection and conduction)
∆$
q̇ = 3.91
V
In a similar manner to calculating electrical resistance, the resistance through the parallel
section can be written as:
& & #&
R ./% = %V + V & 3.106
! #
M N M# N #&
R ./% = % O! + & 3.107
! O#
Therefore:
∆$
q̇ = " H I L I L J" H "
3.108
- " -P ! - # Q - n -
EEE. L I" L H! H# In L EED. L
Scenario 2: However, if we assume that the surfaces parallel to heat flow are adiabatic, i.e.,
no heat flows from slab 2 to 3, a better approximation is given by:
Rh R1 R2 R4 Rhc
h
Rh R1 R3 R4 Rhc
h
Now:
∆$
q̇ = 3.91
V
Adding all the sections which are in series, inverting to account for addition of multiple
parallel sections and again inverting the sum of this:
∆$
q̇ = J" 3.109
" "
_ " H" H! Hn " - " H" H# H "
`
+ + + + + + + n +
EEE. L I" L I! L In L EED. L EEE. L I" L I# L In L EED. L
Th
Tc
k1 k2
q̇
L1 L2
R1 RC R2
where:
RC = Contact resistance (m2.K/W)
Th
Tc
q̇
Various estimations of the heat transfer through the material are available, including ways
such as dividing the material into parallel and series systems (not presented here), but a
more comprehensive method is given as:
where:
keff: Effective heat transfer coefficient across slid and porous section of slab
(W/(m.K))
where:
kf: Heat transfer coefficient of the fluid in the void (W/(m.K))
ks: Heat transfer coefficient of the solid (W/(m.K))
e: Void fraction (-)
Fins may come in various forms. These could be as straight, annular or pins.
Figure 3.4: Some examples of fins: Rectangular, Pins, Triangular, Parabolic, Annular.
TW TA
Convection
y
Conduction
TS
w
x x + Dx
Convection
TW TA
Convection
Conduction
r
TS
x x + Dx
Simplifying assumption (for rectangular fins of uniform cross-sectional area): Assume thin
slice of fin that has a temperature profile in x-direction only, NOT in y-direction.
By an energy balance:
Energy in (left – conduction) = energy out (right - conduction) + energy out (top + bottom
- convection) 3.112
q̇ |) = q̇ |)-∆) + h(2w∆x)(TF − TN ) 3.113
where:
q̇ : Conductive heat transfer
h: convective heat transfer coefficient
Dx: element through which heat flows
TS surface temperature
TA: Ambient/air temperature
Taking limits as Dx à 0
}~ }~
#M(d.+) e -M(d.+) e
}• •+∆• }• •
lim = −h(2w)(TF − TN ) 3.117
∆)→/ ∆)
D D$
2−k(w. y) D) 3 = −h(2w)(TF − TN ) 3.118
D)
Assuming k is constant:
D! $
−k(w. y) D)! = −h(2w)(TF − TN ) 3.119
Simplifying:
D! $
ky D)! = 2h(TF − TN ) 3.120
And:
D! $
ky D)! = 2h(T − TN ) 3.121
D! $ .4
= M+ (T − TN ) 3.122
D)!
Let:
q = T – TA 3.123
And:
m = i2h/ky 3.124
D! g
= m. θ 3.125
D)!
D. − m. θ = 0 3.126
(D − m)(D + m)θ = 0 3.127
θ = C& eG) + C. e#G) 3.128
From BC1:
q1 = C1 + C2 3.129
From BC2:
Dg
= C& meG) − C. me#G) 3.130
D)
0 = C1memL – C2me-mL 3.131
C2 = C1e2mL 3.132
g"
C& = 3.133
&- B!ۥ
g" .B!ۥ
C. = 3.134
&- B!ۥ
g E>F4 G.(O#))
= 3.138
g" E>F4 GO
where:
TW: Wall temperature (K)
cosh i2h/ky. (L − x)
∴ T = (Ti − TN ) p r + TN
cosh i2h/ky. L
where:
q̇ : Heat Flux (W)
A: Cross Sectional Area (m2)
σ: Stefan Boltzmann Constant = 5.67 × 10-8 W.m-2K-4
2π6 k 5j
=
15h% c .
kB = Boltzmann constant – 1.38 × 10-23 J/K;
Kevin Harding (kevin.harding@wits.ac.za) 94
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Heat Exchangers
What is the radiative heat transfer per square meter for a body at the following temperatures?
a) 100°C
b) 400°C
c) 800°C
Solution:
𝑞̇ /𝐴 = 𝐴𝜎𝑇 -
𝑞̇ /𝐴 = 𝜎𝑇 -
𝑞̇ /𝐴 = 5.76 × 10%f . 𝑇 -
a) 𝑞̇ /𝐴 = 1.12 kW/m2
b) 𝑞̇ /𝐴 = 11.8 kW/m2
c) 𝑞̇ /𝐴 = 76.4 kW/m2
Given that the Stefan Boltzmann constant is very small, radiative heat transfer only plays
a part at high temperatures (approx. greater than 400°C), while conduction and convection
describe heat transfer at lower temperatures.
A heat exchanger is a piece of equipment built for efficient heat transfer from one medium
to another.
Figure 3.7: Plate and frame heat exchanger (metal plates transfer heat between two fluids)
(https://commons.wikimedia.org/wiki/File:Plate_frame_1.svg)
Figure 3.8: Passive heat sink on a motherboard (metal pins act as fins to allow for cooling)
(https://commons.wikimedia.org/wiki/File:Harumphy.dg965.heatsink.jpg)
Figure 3.9: Typical radiator found in a car (liquid inside the mechanism allows for conductive cooling
despite the contradictory naming) (https://commons.wikimedia.org/wiki/File:Automobile_radiator.jpg)
A shell and tube heat exchanger is a class of heat exchanger most commonly used in oil
refineries and other large chemical processes. It consists of a shell with a bundle of tubes
inside. One fluid runs through the tubes, and another fluid flows over the tubes (through
the shell) to transfer heat between the two fluids.
Figure 3.11: Double pipe heat exchanger (one Figure 3.12: Shell and tube heat exchanger (one
pass tube side; one pass shell side; no baffles) pass tube side; one pass shell side; with baffles)
(https://commons.wikimedia.org/wiki/File:Double- (https://commons.wikimedia.org/wiki/File:Straight-
Pipe_Heat_Exchanger.png) tube_heat_exchanger_1-pass.png)
Figure 3.13: Shell and tube heat exchanger (two Figure 3.14: U-tube heat exchanger (two pass
pass tube side; one pass shell side; with baffles) tube side; one pass shell side; with baffles)
(https://commons.wikimedia.org/wiki/File:Straight- (https://commons.wikimedia.org/wiki/File:U-
tube_heat_exchanger_2-pass.png) tube_heat_exchanger.png)
It is common to put the dirtier fluid (fluid which would foul) inside the tubes as the tubes
are easier to clean or replace. The easiest way to increase the area of a heat exchanger could
simply be to in increase the number of tubes.
Temperature
Length
Figure 3.15: Temperature profile through a co-current heat exchanger
Typically, the maximum temperature difference at any point can be no lower than 10°C.
Note in the counter current setup that the outlet temperature difference can be less than
this.
Temperature
Length
Tc Cold
Tc + dTc Tc
dA
Assumptions:
- No phase change
- No reaction
- Constant pressure
- Constant specific heat (Cp)
For the portion over dA (assuming the hot and cold fluids don’t change over the length of
the heat exchanger – because the length is so small):
dq̇ = UdA(T4 − TE ) 3.142
(4) – (5)
dq̇ /ṁE CvpE 1 − dq̇ /ṁ4 Cvp4 1 = −ṁ4 Cvp4 . dT4 /ṁE CvpE 1 + ṁE CvpE . dTE /ṁ4 Cvp4 1 3.149
dq̇ /ṁE CvpE − ṁ4 Cvp4 1 = −ṁ4 Cvp4 . ṁE CvpE (dT4 − dTE ) 3.150
Substituting (1)
−UdA(T4 − TE )/ṁE CvpE − ṁ4 Cvp4 1 = ṁ4 Cvp4 . ṁE CvpE (dT4 − dTE ) 3.151
Rearranging:
" lD #ĠE k
ĠD k " lE D$E #D$D
−UdA " lE .ĠD k
" lD
= 3.152
ĠE k $E #$D
Now:
" lD #ĠE k
ĠD k " lE
Let " lE .ĠD k
" lD
=R (6) 3.153
ĠE k
Or
" lD #ĠE k
ĠD k " lE & &
R= " "
= %Ġ "
− Ġ " lD
& 3.154
ĠE klE .ĠD klD E klE Dk
D$E #D$D
∴ −RUdA = 3.155
$E #$D
D($E #$D )
∴ = −RUdA 3.156
$E #$D
Integrate:
∆$
ln(T4 − TE )∆$!" = −RUA 3.157
∆$
ln ∆$! = −RUA 3.158
"
Or:
∆$
ln ∆$" = RUA 3.159
!
∆~"
[,
∆~!
R= (7) 3.160
]N
(8) – (9)
& &
Q %Ġ "
− Ġ " lD
& = nT4,=, − T4,>?@ o − nTE,>?@ − TE,=, o 3.163
E klE Dk
& &
Q %Ġ "
− Ġ " lD
& = nT4,=, − TE,>?@ o − nT4,>?@ − TE,=, o 3.164
E klE Dk
Substituting (6):
QR = ∆T& − ∆T. 3.165
∆$" #∆$!
Q= 3.166
V
Substituting (7)
∆$" #∆$!
Q = UA p ∆~" r 3.167
[,
∆~!
where:
Kevin Harding (kevin.harding@wits.ac.za) 101
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Heat Exchangers
∆$" #∆$!
p ∆~" r = ∆TOH 3.168
[,
∆~!
∴ Q = UA∆TOH 3.169
Assumptions:
i) Overall heat transfer co-efficient (U) is constant
ii) Heat exchanger is well insulated
iii) Physical properties of each fluid are constant
iv) Steady state
v) No phase change, no reaction, no mixing, constant pressure
Phase change; no
temperature change
DTLM,a
DTLM,b
DTLM,c
CAN NOT use DTLM calculation BUT: If we break the problem into 3,
can use DTLM calculation
Example 3.20: Heat exchanger – log mean temperature vs. temperature difference
Calculate the log mean temperature difference given the following information for a counter current
system:
Temperature hot in: 90°C
Temperature hot out: 65°C
Temperature cold in: 40°C
Temperature cold out: 70°C
Solution:
90°C
Temperature 65°C
DT2
70°C DT1
40°C
Length
∆𝑇& − ∆𝑇( 25 − 20
∆𝑇[ƒ = Š ‹= ƒ „ = 22.4℃
∆𝑇 25
𝑙𝑛 & 𝑙𝑛 ) *
∆𝑇( 20
(s0(\
This is not the same as average temperature: ∆𝑇`„… = ) (
* = 22.5℃
Calculate the log mean temperature difference given the following information for a counter current
system Temperature hot in: 90°C
Kevin Harding (kevin.harding@wits.ac.za) 103
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Heat transfer – Heat Generation
Solution:
90°C
Temperature
DT2 60°C
70°C
DT1
40°C
Length
∆𝑇# − ∆𝑇$ 20 − 20
∆𝑇/1 = C F= G J = 𝑢𝑛𝑑𝑒𝑓𝑖𝑛𝑒𝑑
∆𝑇# 20
𝑙𝑛 ∆𝑇 𝑙𝑛 H20I
$
Since DTLM is undefined, can use DTAVE. Also: since DT is constant across the length of heat exchanger (and also equal
to DTAVE), DTAVE can be used.
20 + 20
∆𝑇)23 = O P = 20℃
2
3.9.2 Catalysts
Metals are expensive. To maximize the surface area to volume ratio of a catalyst, we
typically use metal oxide supports which have a large surface area, e.g., activated carbon.
Reactants adsorb and react in the pores; products desorb and are released.
Consider an electric wire of circular cross section with radius R and electrical conductivity ke ohm-1cm-1.
Through this wire an electric current is passed with a current density of I amp/cm2. This process is an
irreversible process, converting some electrical energy into heat (thermal energy). The rate of heat
production per unit volume is given by the expression:
𝐼(
𝑆i =
𝑘i
The quantity Se is the heat source resulting from electrical dissipation. It is assumed that the temperature
rise in the wire is not so large that the temperature dependence of either the electrical or thermal
conductivity need be considered. The surface of the wire is maintained at temperature T0.
R
L
Dr
Solution:
Taking limits as Dr à 0
(𝑟𝑞T ) |T0∆T − (𝑟𝑞T ) |T
𝑙𝑖𝑚 = 𝑆i 𝑟
∆T→\ ∆𝑟
𝑑
(𝑟𝑞 ) = 𝑆i 𝑟
𝑑𝑟 T
Equating:
𝑑𝑇 𝑆i 𝑟
−𝑘. 0 7 =
𝑑𝑟 2
At r = R; T = T0
𝑆i 𝑅(
𝐶( = + 𝑇\
4𝑘
𝑆i 𝑅( 𝑟 (
𝑇 − 𝑇\ = ?1 − ) * @
4𝑘 𝑅
Maximum temperature is at r = 0
ˆq ‰ ! ˆq ‰ !
𝑇‡W$ − 𝑇\ = -e
OR 𝑇‡W$ = -e
+ 𝑇\
Average temperature:
(c ‰
∫\ ∫\ (𝑇(𝑟) − 𝑇\ )𝑟𝑑𝑟. 𝑑𝜃 𝑆i 𝑅(
〈𝑇〉 − 𝑇\ = (c ‰ =
∫\ ∫\ 𝑟𝑑𝑟. 𝑑𝜃 8𝑘
q̇ &,=, q̇ .,>?@
q̇ &,>?@ q̇ .,=,
1 2
A well-stirred spherical vessel of diameter 6.5 cm containing water at 50°C is placed in a constant
temperature bath at 20°C. The overall heat transfer coefficient is 625 W/m2.K.
a) Assuming a constant heat loss over time, how much energy is lost from the sphere in 2 minutes?
b) What is the temperature of the water in the vessel after 2 minutes?
Solution:
a) 𝑞 = 𝑈\ 𝐴(𝑇5 − 𝑇\ )∆𝑡
𝑞 = 625 𝐽𝑠 %& 𝑚%( 𝐾 %& × 4𝜋(0.0325)( 𝑚( × (323𝐾 − 293𝐾) × 120𝑠
𝑞 = 625 𝐽𝑠 %& 𝑚%( 𝐾 %& × 4𝜋(0.0325)( 𝑚( × (323𝐾 − 293𝐾) × 120𝑠
𝑞 = 29865 𝐽 = 29.865 𝑘𝐽
b) Set up a differential balance describing how the temperature of the water within the vessel (the system)
varies with time. This is a constant pressure system => ∆H = q
𝐻|9 = 𝐻|90∆9 + 𝑞̇ ∆𝑡
𝐻 = 𝑚𝐶“X ∆𝑇 = 𝑚𝐶“X (𝑇 − 𝑇Š )
where Tb is the basis temperature relative to which the Cp is calculated
$ &
Note: ∫$ ! W%$ 𝑑𝑥 = [− 𝑙𝑛(𝑎 − 𝑥)]$$"!
"
U
∴ 𝑚𝐶“X [−𝑙𝑛 (𝑇\ − 𝑇5 )]s\r = 𝑈\ 𝐴𝑡
𝑙𝑛 (𝑇\ − 323)
𝑚𝐶“X = 𝑈\ 𝐴𝑡
𝑙𝑛 (𝑇\ − 𝑇4 )
𝑇\ − 323 𝑈\ 𝐴𝑡
𝑙𝑛 ( ,=
𝑇\ − 𝑇4 𝑚𝐶“X
‹b `9
𝑇\ − 323 Ž‡Œ•
( ,=𝑒 s
𝑇\ − 𝑇4
‹b `9
Ž‡Œ•
𝑇\ − 323 = 𝑒 s × (𝑇\ − 𝑇4 )
‹b `9 ‹b `9
Ž‡Œ• Ž‡Œ•
𝑇\ − 323 = 𝑇\ . 𝑒 s − 𝑇4 . 𝑒 s
‹b `9
Ž‡Œ•
𝑇\ (𝑒 s − 1, + 323
∴ 𝑇4 = ‹b `9
Ž‡Œ•
𝑒 s
r(s×-c(\.\,(s)! ×&(\‘
• \.&-×-&f(’
293 × (𝑒 − 1, + 323
𝑇4 = r(s×-c(\.\,(s)! ×&(\‘
• \.&-×-&f(’
𝑒
𝑇4 = 298.5𝐾 = 25.5℃
Additional Reading
Çengel, YA, Ghajar, AJ, 2011. Heat and Mass Transfer – Fundamentals and Applications,
McGraw Hill.
Coulson, JM, Richardson, JF, Backhurst, JR, Harker, JH, 1999. Coulson & Richardson’s
Chemical Engineering, Volume 1, Elsevier.
Incropera, FP, DeWitt, DP, Bergman and Lavine, 2013. Fundamentals of Mass and Heat
Transfer, 7th Edition, John Wiley & Sons.
Perry, RH, and Green, DW, 1984. Perry's Chemical Engineers Handbook.
Rathakrishnan, E, 2012. Elements of Heat Transfer, CRC Press, Taylor & Francis Group.
Skogestad, S, 2009. Chemical and Energy Process Engineering, CRC Press.
Welty, JR, Wicks, CE, Wilson, RE, Rorrer, GL. 2008. Fundamentals of Momentum, Heat,
and Mass Transfer, Wiley.
Problems
e.) None of the above How is a phase change dealt with in a heat exchanger? In order to calculate the heat
Example 3.24: [MCQ] Conduction: Definition of A exchanger area needed:
Solution: a.) Calculate the log mean temperature for three sections (liquid,
Conduction is given by the equation: q = -k.A.dt/dx, where x is the direction of heat B phase change and gas) and use the average of these temperatures
flow. A is the: in a single calculation.
a.) Average surface area of the object, in m2 Example 3.28: [MCQ] Units of k b.) Use the log mean temperature across the whole exchanger.
b.) Area parallel to heat flow, in m2 c.) Calculate the log mean temperature for three sections (liquid,
c.) Area perpendicular to heat flow, in m2 What are the units of k (thermal conductivity coefficient): phase change and gas) and calculate 3 individual areas.
d.) Amount of material through which heat flows, in kg a.) J/(m.K.s) d.) Divide the calculation into one part for the liquid and one for the
e.) Total amount of material through which heat flows PLUS b.) W/(kg.K) gas.
insulating material which does not allow heat flow, in kg c.) J/(m.K) e.) Use the average (NOT log mean) temperature across the whole
d.) W/K exchanger.
Solution: e.) J/(m2.K)
C Solution:
Solution: C
Example 3.25: [MCQ] Definition of thermal insulator C
Example 3.32: [MCQ] Overall heat transfer coefficient
A material that is considered a thermal insulator will have a: Example 3.29: [MCQ] Conduction through the walls of a cylinder
a.) High thermal conductivity UA is given as 100 W/K. Ao = 20 m2 and Ai = 10m2. The correct Uo and Ui respectively
b.) Low thermal conductivity Heat flows through a cylinder (from inside to outside). Which of the following are (in W/m2.K):
c.) Large cross-sectional area equations would you use to solve heat flow: a.) 2&5
d.) Large surface area a.) 𝑞 = −𝑘2𝜋𝑟𝐿(𝑇L − 𝑇B ) b.) 2 & 10
e.) None of the above b.) 𝑞=
&J"UV(W8&W.)
2
c.) 5 & 10
?+ 6!960,
234/960,
d.) 5 & 20
&J"UV(W8&W.)
Solution: c.) 𝑞 = (; & & e.) 10 & 20
34/960,) &(;6!960,)
B d.) 𝑞=
&J"UV(W8&W.)
;34/960,&;6!960,
Solution:
e.) None of the above
Example 3.26: [MCQ] Definition of 𝒒̇ C
Solution:
𝑞̇ has the units of: Example 3.33: Heat transfer around a pipe
B
a.) kJ
b.) kJ/kg Water flows through a 2” steel pipe (2.067” ID & 2.375” OD). The pipe is heated on
Example 3.30: [MCQ] Heat transfer equipment
c.) kJ/m3 the outside by steam at 230°F. The overall heat transfer coefficient is 250 Btu/hr ft2
d.) kJ/s °F based on the outside area of the pipe. If the inlet and outlet temperatures for a 1 ft
Which of the following is least likely to be classified as heat transfer equipment:
e.) kJ/m2 length of pipe are 80°F and 180°F, what is the water flow rate?
a.) Fin fan cooler
b.) Distillation column
Solution: Solution:
c.) Radiator
D 𝑄 = 𝑈𝐴∆𝑇VQ
d.) Heat Exchanger
𝐴 = 2𝜋𝑟$ 𝐿
e.) Furnace
Example 3.27: [MCQ] Conduction: Proportionality
∆𝑇# − ∆𝑇" }𝑇EA@:=,F+ − 𝑇N"$,>HA ~ − }𝑇EA@:=,>HA − 𝑇N"Y,F+~
Solution: ∆𝑇VQ = =
∆𝑇 }𝑇EA@:=,F+ − 𝑇N"Y,>HA ~
Conductive heat transfer is proportional to: 𝑙𝑛 # 𝑙𝑛
B ∆𝑇"
}𝑇EA@:=,>HA − 𝑇N"Y,F+ ~
a.) Dx/DT
b.) DT/Dx
Example 3.31: [MCQ] Heat exchanger: Phase change Assume that Tsteam,in = Tsteam,out = 230°F
c.) DA/DT
d.) DT/DA
(230℉ − 180℉) − (230℉ − 80℉) Cp, oil = 1.6 kJ/kg.K (1m – 3(5mm)/2 = 0.4925 mm thick. Further assume that the section is 1m deep.
∴ ∆𝑇VQ = = 91℉
(230℉ − 180℉) Cp, steam = 2 kJ/kg.K NOTE: Units DO NOT balance with this assumption and as such, the heat transfer
𝑙𝑛
(230℉ − 80℉)
given needs to be written as 4kW as we are taking a basis of 1 m.
Solution:
𝑄 = 𝑈𝐴∆𝑇VQ
1𝑓𝑡 a) At a viscosity of 6508 cP; x = 119.8653 m [2nd equation] Numbering the sections from 1 – 8 from left to right (top to bottom):
&# &" &#
𝑄 = 250 𝐵𝑡𝑢 ℎ 𝑓𝑡 ℉ × 2𝜋 × 1.1875𝑖𝑛 × × 91℉ This also corresponds to a temperature of 30°C [1st equation]
12𝑖𝑛
&#
&#
1 𝑘𝐽𝑠 Calculating the heat needed to raise the oil from 30°C to 80°C: Li Ki Ai Resistance per
𝑄 = 14 145.35 𝐵𝑡𝑢 ℎ × = 4.146 𝑘𝐽𝑠 &# section
3412 𝐵𝑡𝑢 ℎ&# 𝑄 = 𝑚𝐶𝑝∆𝑇 = 100 𝑘𝑔 ∗ 1.6 𝑘𝐽/𝑘𝑔. 𝐾 ∗ (80 − 30)𝐾 = 7999.529 𝑘𝐽/ℎ𝑟
(mm) (m) (w/mK) (m2) Ri = Li/KiAi
Since the heat gained by the oil must be equal to the heat lost by the hot
1 Ceramic 40 0.04 0.52 1 0.076923077
At T = 27°C, hf = 113.1 kJ/kg water: 2 Al left 20 0.02 200 1 0.0001
𝑄>F? = 𝑄Z:A@; = 7999.529 𝑘𝐽/ℎ𝑟 = 𝑚 ∗ 4.1868 ∗ (10) 3 Al top 500 0.5 200 0.005 0.5
At T = 80°C, hf = 334.9 kJ/kg ∴ 𝑚Z:A@; = 399.976 ≈ 400 𝑘𝑔/ℎ𝑟 4 Air top 500 0.5 0.0234 0.4925 43.38583019
At T = 85°C, hf = 355.9 kJ/kg b) Length at which viscosity is at the limit = 119.86m (from above). Therefore, 5 Al middle 500 0.5 200 0.005 0.5
6 Air bottom 500 0.5 0.0234 0.4925 43.38583019
By interpolation, hf = 343.3 kJ/kg at T = 82°C need a heat exchanger every 119 m 7 Al bottom 500 0.5 200 0.005 0.5
Therefore (3716/119) = 31 8 Al right 24 0.024 200 1 0.00012
Energy balance: Plus adding the initial heat exchanger to the start of the line, a total of 32 heat
∆𝐻 „N"Y = 𝑞… exchangers are needed. Now:
∴ 𝑞̇ = 𝑚̇N"Y ∆𝐻 „N"Y ∆𝑇
𝑞̇ =
𝑞̇ = 𝑚̇N"Y }ℎ†-,)"℃ − ℎ†-,"2℃ ~ Example 3.35: Heat transfer through composite slab 𝑅
𝑞̇ 4.146 𝑘𝐽𝑠 &#
⇒ 𝑚̇N"Y = = = 0.01801 𝑘𝑔𝑠 &#
}ℎ†-,)"℃ − ℎ†-,"2℃ ~ (343.3 𝑘𝐽𝑘𝑔&# − 113.1 𝑘𝐽𝑘𝑔&#) Assuming a section of composite material below, what is the inside wall temperature Or for the system above:
𝑚̇N"Y = 0.01801 𝑘𝑔𝑠 &# × 2.205 𝑙𝑏𝑘𝑔&# × 3600 𝑠 ℎ&# = 143 𝑙𝑏 ℎ&# ∆𝑇
if the outside is maintained at a temperature of 1000°C and the heat transfer is 𝑞̇ =
𝐿# 𝐿 𝑘 𝐴 𝑘 𝐴 𝑘 𝐴 𝑘 𝐴 𝑘 𝐴 &# 𝐿
4 kW/m. Give your answer to the nearest degree Fahrenheit. + " + 7 5 5 + 7 7 + * * + 3 3 + 2 2; + )
𝑘# 𝐴# 𝑘" 𝐴" 𝐿5 𝐿7 𝐿* 𝐿3 𝐿2 𝑘) 𝐴
k-values (W/mK):
Aluminium: 200 Ceramic: 0.52 Air: 0.0234 Which can also be simplified to (because portions of the parallel section are the same
material, i.e., same k, AND same area):
Assumptions to be used: ∆𝑇
𝑞̇ =
- Ignore convection 𝐿# 𝐿" 𝑘5 𝐴5 𝑘7 𝐴7
&#
𝐿
𝑘# 𝐴# + 𝑘" 𝐴" + Ž3 - 𝐿5 . + 2 - 𝐿7 .• + 𝑘 )𝐴
- Assume the surfaces normal to the direction of heat flow are isothermal )
Solution: Solution:
Given: Heat exchangers A process stream is heated or cooled via exchange of heat
It can be assumed that all three “bars” are 5 mm thick. Therefore, the length of the
𝑇 = 25 + 55𝑒 &$.$"/
Units: °C with another process stream
air sections are:
𝑉𝑖𝑠𝑐𝑜𝑠𝑖𝑡𝑦 = 5000 + 5000𝑒 &$.$#/ Units: cP
Where x is the length down the pipe (in metres)
Heaters A process stream is heated via exchange of heat with a Hot always A stream (oil) of 9.2 kg/s is to be heated from 65°C to 95°C in a heat exchanger, using
hotter than
hot service/utility stream (e.g. steam) cold 16.7 kg/s of water entering at 105°C. If the overall heat-transfer coefficient is
Hot out
Coolers A process stream is cooled via exchange of heat with a lower than
1300 J/m2s°C, calculate the area of heat exchanger required if the flows are counter
cold service/utility stream (e.g. cooling water) cold out
current.
Reboilers Typically associated with distillation columns. The
liquid bottoms product is heated and partially or fully Solution:
vaporized via exchange of heat with a hot service/utility mH,in = 16.7kg/s mH,out = 16.7kg/s
TH,in = 105°C TH,out = ? °C
stream (e.g. steam)
Condensers Typically associated with distillation columns. The mC,out = 9.2kg/s mC,in = 9.2kg/s
vapour top product is heated and partially or fully TC,out = 95°C TC,in = 65°C
Example 3.38: Counter current heat exchanger
condensed via exchange of heat with a cold
For the oil stream:
service/utility stream (e.g. cooling water/refrigerant)
A small counter current heat exchanger operates with the following stream 𝑄 = 𝑚𝐶‘G ∆𝑇
Agitated vessels (jacketed heaters) Process fluid is heated or cooled in an
temperatures: Assume vegetable oil or mineral oil (both have the same heat capacity)
agitated tank via exchange of heat with a hot or cold
Cp,oil = 1.67 kJ/kg.K
jacket (steam or cooling water circulated through the
Cold stream in: 20°C 𝑄 = 9.2 𝑘𝑔𝑠 &# × 1.67 𝑘𝐽𝑘𝑔&# 𝐾 &# × (368𝐾 − 338𝐾)
outer jacket of the vessel) 𝑄 = 460.92 𝑘𝐽𝑠 &#
Cold stream out: 100°C
Hot stream in: 120°C
Example 3.37: Shell and tube heat exchanger For the water stream:
Hot stream out 70°C
𝑄 = 𝑚𝐶‘G ∆𝑇
In a shell and tube heat exchanger: 𝑄
The unit has a total area for heat transfer of 1 m² and an overall heat transfer ∴ ∆𝑇 =
𝑚𝐶‘G
coefficient of 500W/m².K. What is the rate of energy transfer?
a) What are the role of baffles?
b) How does fouling influence heat transfer? Solution: Assume Cp,water = 4.19 kJ/kg.K
c) Can the outlet temperature of a hot stream be less than the outlet temperature 460.92 𝑘𝐽𝑠 &#
∆𝑇 =
of the cold stream? Explain and include a sketch to show this. TH,in = 120°C 16.7 𝑘𝑔𝑠 &# × 4.19 𝑘𝐽𝑘𝑔&# 𝐾 &#
TH,out = 70°C ∆𝑇 = 6.6𝐾
Solution: ∆𝑇 = 𝑇N,>HA − 𝑇N,F+
a) Baffles define the flow path of fluid through the heat exchanger, increase the TC,out = 100°C ∴ 𝑇N,>HA = 371.4𝐾 = 98.4℃
TC,in = 20°C
surface area for heat exchange and support the tubes.
𝑄 = 𝑈$ 𝐴∆𝑇VQ 𝑄 = 𝑈$ 𝐴∆𝑇VQ
}𝑇N,F+ − 𝑇R,>HA ~ − }𝑇N,>HA − 𝑇R,F+ ~ (393𝐾 − 373𝐾) − (343𝐾 − 293𝐾) 𝑄
b) Fouling refers to the deposition of solid material on the heat exchange surfaces. ∆𝑇VQ = = ∴𝐴=
}𝑇N,F+ − 𝑇R,>HA ~ (393𝐾 − 373𝐾) 𝑈$ ∆𝑇VQ
All process and utility streams will cause fouling and the extent depends on the 𝑙𝑛 Ž • 𝑙𝑛 7 ;
}𝑇N,>HA − 𝑇R,F+ ~ (343𝐾 − 293𝐾) }𝑇N,F+ − 𝑇R,>HA ~ − }𝑇N,>HA − 𝑇R,F+ ~ (105℃ − 95℃) − (98.4℃ − 65℃)
composition of the fluid. The fouled material acts as a resistance to heat transfer ∆𝑇VQ = =
= 32.74𝐾 }𝑇N,F+ − 𝑇R,>HA ~ (105℃ − 95℃)
and, therefore, reduce the rate of heat transfer. ln Ž • ln 7 ;
}𝑇N,>HA − 𝑇R,F+ ~ (98.4℃ − 65℃)
b) What is the energy input required at the ends of the rod to maintain =0
the specified temperatures T1 and T2? b) For conduction down the length of the rod: 𝑞̇ = ℎ(2𝜋𝑅𝐿)(𝑇" − 𝑇#)
Using the quadratic formula:
𝑞̇ −𝑏 ± √𝑏 " − 4𝑎𝑐
Solution: =
𝐿 2𝑎
a) Energy in|x = Energy out|x+∆x – energy loss due to convection over ∆x Example 3.41:
Note:
𝑞̇ |/ = 𝑞|̇/,∆/ + ℎ(2𝜋𝑅∆𝑥)(𝑇 − 𝑇$ ) A thick-walled pipe of inner radius r1 and outer radius r2 is constructed of material k (W/mK) = a (W/mK) + b (W/mK2) * T (°C)
whose thermal conductivity is a function of the temperature, i.e. k = a + bT, where a 𝑏¢
𝑊𝑖𝑡ℎ 𝑎 = 2
Applying Fourier’s Law: and b are constants. Inside the pipe, there is a liquid at temperature To. Outside, the (2𝜋ℎ" 𝑟")"
𝑑𝑇 𝑑𝑇 pipe is surrounded by air at a temperature TA. The film coefficients on the inside and 0.004 𝑊𝑚&# 𝐾 &"¢
−𝑘(𝜋𝑅" ) — = −𝑘(𝜋𝑅" ) — + 2𝜋ℎ𝑅∆𝑥(𝑇) 2
𝑑𝑥 / 𝑑𝑥 /,∆/ outside of the pipe are h1 and h2 respectively. 𝑎= = 2.0264 × 10&3 𝑊 &# 𝑚
(2𝜋 × 25 𝑊𝑚&" 𝐾 &# × 0.20 𝑚)"
Re-arranging and dividing throughout by ∆x:
𝑑𝑇 𝑑𝑇
𝑘𝜋𝑅" - ˜ − a.) If the coefficient h1 is very large, so that the temperature of the inner surface 𝑟
𝑑𝑥 /,∆/ 𝑑𝑥 ˜/ . ln} "¢𝑟# ~ 𝑎 + 𝑏𝑇b
= 2𝜋ℎ𝑅𝑇 may be assumed to be at To, calculate the steady-state heat flow through the 𝑊𝑖𝑡ℎ 𝑏 = +
∆𝑥 2𝜋 2𝜋ℎ" 𝑟"
Taking lim pipe wall per unit length. 1 }0.1 + 0.004(10℃)~𝑊𝑚&# 𝐾 &#
Žln-0.20 𝑚¢0.05 𝑚. +
∆/→$
𝑑 𝑑𝑇 2ℎ𝑇 b.) Recalculate 𝑞̇ /𝐿 allowing for h1 given in the data below. 𝑏= • = 0.22509
2𝜋 (25 𝑊𝑚&" 𝐾 &# × 0.20 𝑚)
7 ;=
𝑑𝑥 𝑑𝑥 𝑘𝑅
𝑑 " 𝑇 2ℎ𝑇 𝑑 " 𝑇 2ℎ 𝑏
∴ "= 𝑜𝑟 − 𝑇=0 Data:
𝑑𝑥 𝑘𝑅 𝑑𝑥 " 𝑘𝑅 𝑊𝑖𝑡ℎ 𝑐 = 𝑎(𝑇$ − 𝑇b ) + }𝑇$ " − 𝑇b " ~
r1 = 5 cm 2
If D = d/dx and we let θ = (2h/kR)½ , then: 0.004
(𝐷" − 𝜃 " )𝑇 = 0 r2 = 20 cm 𝑐 = 0.1 𝑊𝑚&# 𝐾 &# (100℃ − 10℃) + 𝑊𝑚&# 𝐾 &"(100℃" − 10℃" )
2
(𝐷 + 𝜃)(𝐷 − 𝜃)𝑇 = 0 h1 = 150 W/m2.K = 28.8 𝑊𝑚 &#
⇒ 𝑇 = 𝐴𝑒 &]/ + 𝐵𝑒 ]/ h2 = 25 W/m2.K
To = 100°C `̇ &$.""*$'±d$.""*$'&&7(".$"37×#$;< f;#=)(&").) f=;#)
= = 128 𝑊𝑚&#
TA = 10°C V "(".$"37×#$;< f;#=)
Boundary condition 1: At x = 0, T = T1 = 150°C
150 = 𝐴𝑒 &]($) + 𝐵𝑒 ]($) = 𝐴 + 𝐵 k = 0.1 + 0.004T in W/mK for T in °C
∴ 𝐵 = 150 − 𝐴 b) Film heat transfer from fluid to inner pipe wall:
Boundary condition 2: At x = L, T = T2 = 60°C Solution:
60 = 𝐴𝑒 &]V + (150 − 𝐴)𝑒 ]V a) If h1 is very large, it offers minimal resistance to heat transfer 𝑞̇ = ℎ# (2𝜋𝑟# 𝐿)(𝑇$ − 𝑇#)
`̇
𝐴(𝑒 ]V − 𝑒 &]V ) = 150𝑒 ]V − 60 𝑇# = 𝑇$ − (1)
"UL#;#V
150𝑒 ]V − 60 Conduction through the pipe material:
𝐴 = ]V
𝑒 − 𝑒 &]V 𝑑𝑇
𝑞̇ = −𝑘𝐴 Conduction through the pipe material:
𝑑𝑟 𝑑𝑇
150𝑒 ]V − 60 ;& W& 𝑞̇ = −𝑘𝐴
∴ 𝐵 = 150 − ]V 1 𝑑𝑟
(𝑒 − 𝑒 &]V ) ∴ 𝑞̇ ž 𝑑𝑟 = − ž (𝑎 + 𝑏𝑇)𝑑𝑇 ;& W&
2𝜋𝑟𝐿 1
150(𝑒 ]V − 𝑒 &]V ) − (150𝑒 ]V − 60) ;# W: ∴ 𝑞̇ ž 𝑑𝑟 = − ž (𝑎 + 𝑏𝑇)𝑑𝑇
𝐵= `̇ ;& 9 9 2𝜋𝑟𝐿
(𝑒 ]V − 𝑒 &]V ) ln = 𝑎𝑇$ + 𝑇"− 𝑎𝑇" − 𝑇" " (1)
"UV ;# " $ " ;# W#
60 − 150𝑒 &]V `̇
ln
;& 9
= 𝑎𝑇# + 𝑇# " − 𝑎𝑇" − 𝑇" "
9
(2)
𝐵= "UV ;# " "
𝑒 ]V − 𝑒 &]V
Film heat transfer from outer surface of pipe:
𝑞̇ = ℎ"(2𝜋𝑟" 𝐿)(𝑇" − 𝑇b) Film heat transfer from outer surface of pipe:
150𝑒 ]V − 60 &]/ 60 − 150𝑒 &]V ]/ `̇
𝑇= 𝑒 + ]V 𝑒 𝑇" = 𝑇b + (2) 𝑞̇ = ℎ"(2𝜋𝑟" 𝐿)(𝑇" − 𝑇b) (3)
𝑒 ]V − 𝑒 &]V 𝑒 − 𝑒 &]V "UL&;&V
]/
(150𝑒 ]V − 60)𝑒 &]/ + (60 − 150𝑒 &]V )𝑒
𝑇=
𝑒 ]V − 𝑒 &]V Substitute (2) into (1): Substitute (1) and (3) into (2):
"
𝑞̇ 𝑟" 𝑏 𝑞̇ 𝑏 𝑞̇ " 𝑞̇ 𝑟" 𝑞̇ 𝑏 𝑞̇ 𝑞̇
At x = L/2 ln = 𝑎𝑇$ + 𝑇$ " − 𝑎 Ÿ𝑇b + − Ÿ𝑇b + ln = 𝑎 Ÿ𝑇$ − + Ÿ𝑇$ − − 𝑎 Ÿ𝑇b +
2𝜋𝐿 𝑟# 2 2𝜋ℎ" 𝑟" 𝐿 2 2𝜋ℎ" 𝑟" 𝐿 2𝜋𝐿 𝑟# 2𝜋ℎ# 𝑟# 𝐿 2 2𝜋ℎ# 𝑟# 𝐿 2𝜋ℎ" 𝑟" 𝐿
"
]V
(150𝑒 ]V − 60)𝑒 & ^" + (60 − 150𝑒 &]V )𝑒
]V^
" Re-arranging to solve for q/L: 𝑏 𝑞̇
𝑇= − Ÿ𝑇b +
𝑒 ]V − 𝑒 &]V 2 2𝜋ℎ" 𝑟" 𝐿
Re-arranging to solve for q/L:
𝑏¢ 𝑏¢ 𝑟" a = 0.25
𝑞̇ " 2 2 £ + 𝑞̇ ¤ln} ¢𝑟# ~ + 𝑎 + 𝑏𝑇$ + 𝑎 + 𝑏𝑇b ¥
7 ; ¡ − x1 = 50 mm
𝐿 (2𝜋ℎ" 𝑟" )" (2𝜋ℎ# 𝑟# )" 𝐿 2𝜋 2𝜋ℎ# 𝑟# 2𝜋ℎ" 𝑟"
x2 = 250 mm
𝑏
− Ÿ𝑎(𝑇$ − 𝑇b ) + }𝑇$ " − 𝑇b " ~ = 0 T1 = 600K
2
T2 = 400K
Using the quadratic formula:
k = 3.46 W/m.K
𝑞̇ −𝑏 ± √𝑏 " − 4𝑎𝑐
=
𝐿 2𝑎
Note: The sides are well insulated, while the top surface of the cone (x1) is maintained at T1,
k (W/mK) = a (W/mK) + b (W/mK2) * T (°C) T2 T1
𝑏¢ q and the bottom surface (x2) at T2. What is the rate of the heat transfer across the cone?
𝑊𝑖𝑡ℎ 𝑎 = 2 Ÿ 1 − 1
(2𝜋)" (ℎ" 𝑟" )" (ℎ# 𝑟#)" Ignore convection.
0.004 𝑊𝑚&# 𝐾 &"¢ 1 1
𝑎= 2Ÿ −
(2𝜋)" (25 𝑊𝑚&" 𝐾 &# × 0.20 𝑚)" (150 𝑊𝑚&" 𝐾 &# × 0.05 𝑚)" Given:
𝑎 = 1.1258 × 10&3 𝑊 &# 𝑚 X2 X1
kaluminium: 238 W/m.K
;
With 𝑏 =
ghi &^;# j
+
:,9W:
+
:,9W= a: 1m-1/2
"U "UL#;# "UL&;& Solution:
1 }0.1 + 0.004(100℃)~𝑊𝑚&# 𝐾 &# a) Assuming, steady-state, heat transfer only in x direction and well insulated lateral x1: 0.075 m
𝑏= Žln-0.20 𝑚¢0.05 𝑚. +
2𝜋 (150 𝑊𝑚&" 𝐾 &# × 0.05 𝑚) surface – no heat transferred to wall. T1: 100°C
}0.1 + 0.004(10℃)~𝑊𝑚&# 𝐾 &#
+ • q = -kAdT/dx x2: 0.225 m
(25 𝑊𝑚&" 𝐾 &# × 0.20 𝑚)
𝑏 = 0.23570 T2: 20°C
With: A = π D2/(4)
But D = ax
𝑏
𝑊𝑖𝑡ℎ 𝑐 = 𝑎(𝑇$ − 𝑇b ) + }𝑇$ " − 𝑇b " ~ ⸫ A = (π/4) × a2.x2 Solution:
2
0.004 DW UM&
𝑐 = 0.1 𝑊𝑚&# 𝐾 &# (100℃ − 10℃) + 𝑊𝑚&# 𝐾 &"(100℃" − 10℃" ) From Fourier’s law: 𝑞 = −𝑘𝐴 D/, with 𝐴 = = 𝜋𝑎" /4)𝑥 5
Substituting: q = -k [(π/4).a2x2]dT/dx 7
2
= 28.8 𝑊𝑚&#
Rearranging:
4q/[-π.k.a2x2]dx = dT Therefore:
𝑞̇ −0.23570 ± ¨0.23570" − 4(1.1258 × 10&3 𝑊 &# 𝑚)(−28.8 𝑊𝑚&# )
= 4𝑞𝑑𝑥
𝐿 2(1.1258 × 10&3 𝑊 &# 𝑚) Integrating: = −𝑘𝑑𝑇
𝜋𝑎" 𝑥 5
𝑞̇
= 122 𝑊𝑚&# [4q/- π.k.a2] × (-1/x + 1/x1) = T(x) – T1
𝐿 T(x) = T1 – (4q/π.k.a2) × (1/x1 – 1/x)
Rearranging and integrating:
Example 3.42: Temperature distribution through a solid conical section At time t2: q = (T1 – T2)(π.k.a2)/[4 × (1/x1 – 1/x2)] 4𝑞 −1 𝑥
Ÿ = −𝑘(𝑇# − 𝑇" )
𝜋𝑎" 2𝑥 " 𝑥#
A conical section fabricated in a wall (as below) has a diameter of D = ax. Assuming Substituting q in into T(x): T(x) = T1 + (T1 – T2) × [(1/x – 1/x1)/(1/x1 – 1/x2)]
2𝑞 1 1
the lateral surface is well insulated, derive an expression for the temperature 𝑇 = 𝑇# + " Ÿ " − "
𝜋𝑎 𝑘 𝑥 𝑥#
distribution assuming one dimensional heat flow. b) Substituting all the constants to solve for q: 2.12 W
𝜋𝑎" 𝑘 𝑇" − 𝑇# 𝜋(1)238 (20 − 100)
𝑞= = . = 189.2613 … 𝑊
a) Show an expression for T(x) in terms of T, q, k, a, and x values AND An Example 3.43: Heat transfer through truncated cone 2 Ÿ 1 − 1 2 [(0.225)&" − (0.075)&"]
𝑥" " 𝑥# "
expression for T(x) in terms of T and x values ONLY
b) Given the following information, calculate the heat transfer rate (in Watts A truncated solid aluminium cone of circular cross section area, has a diameter given
to 3 significant figures) through the cone.
by D = ax3/2
4 Energy balances
4.1 Introduction
As for mass balances, the conservation laws extend in an analogous way for energy. This
chapter will cover basic introductory concepts and then energy balances over non-reacting
systems as well as over reacting systems.
4.1.1 Heat
Heat is the disorganized transfer of energy that results due to a temperature gradient. Heat
is not included in the energy of a system since it is the transfer of energy between systems
while energy that part within a system.
For a closed system heat, energy and work are related by:
DE = Q + W 4.1
Where DE = Ef – E0, or the change in energy from the initial to final states.
4.1.2 Work
Work is typically the transfer of energy through movement due to a force. It is typically
broken down into the following subcategories:
- Expansion work (W∆" ) [mechanical work], also called PV Work (WPV), is the work needed
to change a systems volume.
- Flow work (WFlow) is the work performed when a stream exits or enters a system.
- Shaft work (WS) is the work due to mechanical work changing the pressure, e.g. from
pumps compressors, turbines etc.
- Electrochemical work (Wel) is the work performed due to such things as batteries or fuel
cells.
- Other forms of work (Wother) might include electromagnetic work, sound work, etc.
Where:
Wn is all the non-flow work
4.1.3 Energy
Energy is often referred to as the “ability to perform work”. However, the definition needs
to further extend to include the concept of being able to release heat. Therefore, the
definition of energy can be said to be either the ability to perform work or to release heat.
- Kinetic energy (EK) is the energy due to movement of the centre of mass of an object.
- Rotational energy (ER) is the energy when the centre of mass is at rest but there is rotation.
- Potential energy (EP) is the energy due to an elevated height relative to a reference point.
- Electric energy (Eel) is energy stored in such things as a capacitor.
- Surface energy (ES) is the intermolecular energy creating a surface.
- Internal energy (U) is the energy of molecules. It can be seen as the sum of thermal energy,
chemical bonding energy and latent energy on a molecular level.
For a system with no potential, kinetic, rotational etc. energy, this can be simplified to:
E=U 4.7
Therefore:
DU = Q + W 4.8
where the internal energy term is an ordinary differential, while Q and W are partial
differentials.
Since, from the work term below, all W terms are zero
W = W'[>d + W∆o + WF + WB[ + W>@4BC 4.12
4.1.5 Enthalpy
Enthalpy is defined as:
H = U + PV 4.13
where:
H: Enthalpy (J)
U: Internal Energy (J)
P: System pressure (N/m2)
V: System volume (m3)
Enthalpy is measured according to some known value. This is known as the basis.
Typically, the basis is chosen as elements in their standard state (at 1 atm and 298.15K) à
H = 0.
0: standard state
f: formation
Given a stream of CO(g), flowing at 100 mol/min at a temperature of 298.15 K and pressure of 1 bar,
calculate the enthalpy.
Solution:
Basis: Elements in standard state at 298.15K and 1 bar
Enthalpy flow:
𝑚𝑜𝑙 𝑘𝐽
𝐻̇ = 100 × (−110.53)
𝑚𝑖𝑛 𝑚𝑜𝑙
𝑘𝐽
𝐻̇ = −11 053
𝑚𝑖𝑛
Converting:
𝑘𝐽 1000 𝐽 1𝑚𝑖𝑛
𝐻̇ = −11. 053 × 10, × × = −1.84 × 10s 𝐽/𝑠
𝑚𝑖𝑛 𝑘𝐽 60𝑠
𝐻̇ = −1.84 𝑥 10s 𝑊 = −184 𝑘𝑊
The negative value indicates that energy is being released, i.e. exothermic
In addition to 25°Cand 1atm basis entails elements in standard state having H0f = 0 kJ/mol
(this accounts for formation and breaking of bonds in a reactor).
It is possible to measure these properties under other conditions for use as a basis. However,
25°C and 1 atm is the most convenient for measurement.
Definition:
s!
Cvl = %s$ & 4.14
)t,q
Where:
𝑥} = composition
P = Pressure
i.e., constant composition and pressure
Where:
T1 = reference temperature (K)
T2 = new temperature (K)
$! $! s!
∫$" Cvl dt = ∫$" % s$ & dt = H|$! − H|$" 4.16
)t,q
i.e., if we don’t have Cp data, we can use the enthalpy data (if available)
NOTE:
If Cp is constant, i.e., NOT dependent on T, then:
$!
$!
~ Cvl dt = Cvl T•$ = Cvl (T. − T& ) = Cvl ∆T
$" "
If Cp is a function of temperature the integral becomes more complex and integral tables
may be needed.
Given a stream of CO(g), flowing at 100 mol/min at a temperature of 350 K and pressure of 1 bar, calculate
the enthalpy.
Solution:
∴ 𝐻̇ = −11.053 × 10% 𝑘𝐽. 𝑚𝑖𝑛"# + 𝑛̇ 𝐶X4 ∆𝑇
𝐻̇ = −11.053 × 10% 𝑘𝐽. 𝑚𝑖𝑛"# + [100𝑚𝑜𝑙. 𝑚𝑖𝑛"# × 29.14𝐽. 𝐾 "# 𝑚𝑜𝑙"# × (350 − 298.15)𝐾]
𝐻̇ = −11.053 × 10% 𝑘𝐽. 𝑚𝑖𝑛"# + 151 090.9𝐽. 𝑚𝑖𝑛"#
𝐻̇ = −11.053 × 10% 𝑘𝐽. 𝑚𝑖𝑛"# + 151.0909𝑘𝐽. 𝑚𝑖𝑛"#
𝐻̇ = −10.902 × 10% 𝑘𝐽. 𝑚𝑖𝑛"#
NOTE:
From previous examples, we see that H = -11.053 x 103 kJ/min when T = 298.15K, and when
T = 350 K, H = -10.902 x 103 kJ/min. Therefore, energy is taken in when the reaction
proceeds at a higher temperature of 350K
(More negative to less negative à energy taken in)
Logically, to move from -10°C to -5°C, energy is taken in as we are heating the substances.
Therefore, a change from a negative larger magnitude value to a negative smaller
magnitude value, energy is absorbed or taken in.
Note:
- We cannot measure specific heat directly; we measure enthalpy relative to some basis.
- In practice, enthalpies are measured, and a polynomial model is fitted to the data in order
to obtain specific heat. The model has different constants depending on which units of
temperature, say, were used to make the enthalpy measurements
The data is often compiled into three main parts for water (steam):
- Saturated Water and Steam;
- Superheated and Supercritical Steam; and
- Further Properties of Water and Steam
While colloquially called steam tables, most compilations also include a section on the
properties of other materials such as Mercury, Ammonia and CFCs.
Solution:
(a) 909 kJ/kg
(b) 2743 kJ/kg
(c) 0.06665 m3/kg
(d) 0.0010030 m3/kg
(e) 3052 kJ/kg
(f) 2421 kJ/kg
(g) 572 kJ/kg
(h) 2857 kJ/kg
(i) 2974.5 kJ/kg
(j) 130°C and 2.7 bar
Kevin Harding (kevin.harding@wits.ac.za) 122
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – Steam Tables
We have saturated steam at 10 bar with 𝑉” = 0.02m3/kg. How much liquid and vapor do we have?
Solution:
For saturated steam:
P = 10 bar; Ts = 179.9°C; 𝑉”“ = 0.1944 m3/kg
But the specific volume of the system is only 0.02 m3/kg. Therefore, the system is not only steam; there
must be some liquid present in equilibrium with the steam.
By interpolation:
(0.1128 × 10%( − 0.1114 × 10%( )𝑚, t𝑉”d − 0.01114 × 10%( u𝑚,
𝑘𝑔 𝑘𝑔
=
(10.03 − 7.920)𝑏𝑎𝑟 (10 − 7.920)𝑏𝑎𝑟
(0.1128 × 10%( − 0.1114 × 10%( )𝑚, /𝑘𝑔 × (10 − 7.920)
𝑉”d = + 0.01114 × 10%( 𝑚, /𝑘𝑔
(10.03 − 7.920)
𝑉”d = 1.125 × 10%( 𝑚, /𝑘𝑔
Assuming 1 kg of steam:
𝑀9 𝑉m9 = 𝑀d 𝑉md + (𝑀9 − 𝑀d ) 𝑉m“
𝑀9 𝑉m9 = 𝑀d 𝑉md + 𝑀9 𝑉m“ − 𝑀d 𝑉m“
𝑀d t𝑉md − 𝑉m“ u = 𝑀9 t𝑉m9 − 𝑉m“ u
Steam quality derivation
(0.02 − 0.1944)𝑚, /𝑘𝑔
𝑀d = 1𝑘𝑔
(1.125 × 10%, − 0.1944)𝑚, /𝑘𝑔
𝑀d = 0.9024𝑘𝑔
Therefore:
Mg = (1 – 0.9024) kg AND Mg = 0.0976 kg
Find the enthalpy of steam at T = 175°C and P = 5 bar. Interpolate between saturation point and 200°C.
Solution:
At saturation: T = 151.8°C; P = 5 bar; hg = 2749 kJ/kg
Kevin Harding (kevin.harding@wits.ac.za) 123
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – Steam Tables
Therefore:
2857 − 2749
ℎ“ = 0 7 . (175 − 151.8) + 2749
200 − 151.8
hg = 2801 kJ/kg
Steam quality
The quality of steam is defined as the mass fraction (or percentage) of vapour in the mixture.
G
x = G” 4.17
~
And, by definition:
V = mV$ 4.19
Divide by mT:
" G~
o " • G•
o " ” G”
o
= + 4.21
G~ G~ G~
" "
$ = o• G • + o” G ”
V 4.22
G ~ G ~
We obtain:
G~ G G
− x = G~ − G” 4.23
G~ ~ ~
G~ # G”
1−x= 4.24
G~
G
1 − x = G• 4.25
~
And re-arranging:
$=V
V $u − x. V
$u + x. V
$v 11
$−V
V $u = x. /V
$v − V
$u 1 12
o" #o"•
x= " ” #o"•
13
o
In Out
When considering energy this can be in the form of energy due to material flowing through
the system, work or added heat.
Min Mou
t
where:
M = Mass (kg or kg/s)
W = Work (J or J/s)
Q = Heat (J or J/s)
Rearranging:
0 = Energy from flow=, − Energy from flow>?@ + Q̇ + Ẇ 4.29
But this is only true for steady state – if not at steady state:
(Change in Energy) = (Energy from flow=, − Energy from flow>?@ ) + Q̇ + Ẇ 4.30
Energy is going to be defined as the internal energy (U) of the system/flow terms:
D]
= (U=, − U>?@ ) + Q̇ + Ẇ 4.31
D@
D]
= ((H − PV)=, − (H − PV)>?@ ) + Q̇ + Ẇ 4.33
D@
Rearranging:
D]
= (H=, − H>?@ ) + Q̇ + Ẇ + Ẇqo 4.34
D@
D]
= ∑,=I& H= + Q̇ + Ẇ + Ẇqo 4.35
D@
What about kinetic and potential energy? These are forms of energy due to movement
(flow terms) or change with time (‘/dt’ terms).
Db]- wI - w– c
= ∑,=I&/H= + EM + El 1 + Q̇ + Ẇ + Ẇqo 4.36
D@
Or:
—!
Dx]- HY - yZz
! |!
= ∑,=I& Ṁ { %H
$+ + φ& + Q̇ + Ẇ + Ẇqo 4.37
D@ . =
Typically, the kinetic and potential energy terms are ignored as they do not influence
temperature and work terms are written as a single term, so:
D]
$ 1 + Q̇ + Ẇ + Ẇqo
= ∑,=I&/MH 4.38
D@ =
Or as a molar basis:
D]
= ∑,=I& %NH& + Q̇ + Ẇ + Ẇqo 4.39
D@ − =
Adiabatic process
D]
$ 1 + Ẇ + Ẇqo
= ∑,=I&/MH 4.43
D@ =
Sign Conventions
[Q]: What does positive or negative Q represent?
Heat Transfer
[A]: Positive Q is heat transferred from hot to cold, i.e., out of the system.
Therefore, definition of heat transfer is opposite in Energy Balance calculations to Heat Transfer
calculations.
BEWARE: When looking at both types of approaches in a single problem ensure the correct sign for Q.
Example 4.49:
Using Simpson’s Rule, determine the enthalpy from 20°C to 200°C given the data below.
Solution:
Note:
h = constant value of 30
7 points à 6 intervals (even number)
Therefore, we can use Simpson’s
Now:
(\\
𝐻= Q 𝐶𝑝(𝑥)
(\
Š F
Simpson’s Rule: ∫W 𝑓(𝑥) ≈ , [𝑓(𝑥& ) + 4𝑓(𝑥( ) + 2𝑓(𝑥, ) + ⋯ + 2𝑓(𝑥/%( ) + 4. 𝑓(𝑥/%& ) + 𝑓(𝑥/ )]
Therefore:
(\\
ℎ
𝐻= Q 𝐶𝑝(𝑥) ≈ [𝐶𝑝& + 4𝐶𝑝( + 2𝐶𝑝, + 4𝐶𝑝- + 2𝐶𝑝s + 4𝐶𝑝r + 𝐶𝑝˜ ]
(\ 3
𝐻 ≈ 5305.6𝐽/𝑚𝑜𝑙
NOTE:
If we were asked to find the integral up to 230°C, this would have been an odd number of intervals and
Simpson’s would not work. Therefore, would need to use Simpson’s as above AND then add the last
portion by another technique, e.g., Trapezoidal Rule.
#&
Rln(P/ ) 1
Tj = p + r
∆H|3l T/
where:
TB = the normal boiling point, K
R = the ideal gas constant, 8.314 J/K.mol
P0 = the vapor pressure at a given temperature, atm
DHvap = the heat of vaporization of the liquid, J/mol
T0 = the given temperature, K
O
HHV (BTU/lb) = 14.544C + 62.028. JH − K + 4050S
8
Where C, H, O and S are the weight fractions of Carbon, Hydrogen, Oxygen and Sulphur
respectively. These are typically obtained from flue gas analyses.
AH
Lh VO gMh
P = P/ J1 − K ≈ P/ . exp J− K
T/ RT/
where:
P = Atmospheric pressure (Pa)
1
Lowry, HH (ed), 1945. Chemistry of Coal Utilization, Ch. 4, Wiley: New York
Kevin Harding (kevin.harding@wits.ac.za) 131
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – Additional Reading
Therefore:
(}.7/886).(/./.7}855)
0.0052h (7.%&559).(/.//86) (9.80665). (0.0289644). h
P = 101325 J1 − K ≈ P/ . exp p− r
298.15 (8.31447). (298.15)
Additional Reading
Balmer, RT, 2011. Modern Engineering Thermodynamics, Academic Press.
Elliott, JR, Lira, CT, 2012. Introductory Chemical Engineering Thermodynamics, 2nd ed.,
Pearson International Edition.
Felder RM, Rousseau RW, 2000. Elementary Principles of Chemical Processes, 3rd ed.,
Wiley.
Himmelblau, DM, 2012. Basic Principles and Calculations in Chemical Engineering, 8th
ed., Prentice-Hall International.
Koretsky, MD, 2013. Engineering and Chemical thermodynamics, 2nd ed., Wiley.
Reklaitis, GV, 1983. Introduction to Material and Energy Balances, Wiley.
Rogers, G, Mayhew, Y, 1992. Engineering Thermodynamics – Work and Heat Transfer,
4th ed., Pearson – Prentice Hall.
Rogers, G, Mayhew, Y, 1995. Thermodynamic and Transport Properties of Fluids (SI
Units), 5th ed., Blackwell Publishing.
Sandler, SI, 2006. Chemical, biochemical, and Engineering Thermodynamics, 4th ed.,
Wiley.
Skogestad, S, 2009. Chemical and Energy Process Engineering, CRC Press, Taylor & Francis
Group.
Problems
Example 4.50: Understanding energy b) Assuming that all the energy is converted to kinetic energy, what is the
Solution: water’s speed at the bottom of the waterfall?
Explain your understanding of energy: a) The reference point is taken as the floor. c) The terminal velocity of water is given as roughly 10 m/s (depending on the
a) What is energy (giving examples)? What are the different units of energy? size of the droplets). How does this compare to the answer for (b)?
b) What is power? How does it differ from energy? b) Assuming: d) Assuming that the water droplet is only travelling at 10 m/s when it reaches
- Mass of piston can be neglected the bottom of the waterfall, what happened to the remaining energy?
c) What is the difference between energy and heat?
- No heat transfer between system and surroundings e) If the temperature of the water at the top of the waterfall is 25°C, calculate
d) List as many types of energy as possible. No work done on or by the system
- the temperature of the water at the bottom of the waterfall. Assumption: No
e) For each answer listed above, give an example where it would be important to water evaporates.
understand this energy. Calculating potential energy: 𝐸k = 𝑚𝑔𝑧 f) If water evaporates, would the temperature in (e) be higher or lower?
𝐸†k = 𝑔𝑧 Explain.
Solution:
a) Energy is the ability to perform work or release heat. Energy is a conserved Height equals vertical distance from reference to bottom of piston, plus half the Solution:
quantity. The S.I. unit for energy is the joule (J), whose equivalent unit is the hight of the gas: a) 𝐸†k A>G = 𝑔ℎA>G = 9.8 𝑚𝑠 &" × 270 𝑚 = 2058 𝑚" 𝑠 &" 𝑜𝑟 𝑁𝑚/𝑘𝑔 𝑜𝑟 𝐽/𝑘𝑔
Newton Meter (Nm = kgm2s-2). Other units of energy include: the calorie (cal), 𝐸†k = 9.81 𝑚𝑠 &" × (2 + 0.5(0.1))𝑚 #
b) 𝐸†k A>G = 𝐸†l 9>AA>= = # 𝑣9>AA>= "
British thermal unit (Btu) and the foot-pound (ft lbf). 𝐸†k = 20.11𝑘𝑔. 𝑚" . 𝑠 &" . 𝑘𝑔&# = 20.11 𝑁. 𝑚. 𝑘𝑔&# = 20.11 𝐽. 𝑘𝑔&#
b) Power is the rate at which work is performed. It is a measure of the amount of
energy consumed or transmitted per unit time. c) 𝐸†k = 𝑔ℎ 𝑣9>AA>= = ¨2 × 2058 𝑚" 𝑠 &" = 64.76 𝑚𝑠 &#
c) Heat is the disorganized transfer of energy that results when systems of
different temperatures are contacted. In contrast to heat, work is organised Height equals vertical distance from reference to bottom of piston, plus half the c) It is significantly smaller than the calculated velocity.
energy transfer between systems. Heat is the energy transfer between systems, hight of the gas d) The remaining energy is lost to other forms, mainly evaporation
whereas energy is measured within the system boundary. 𝐸†k = 9.81 𝑚𝑠 &" × (0 + 0.5(0.1))𝑚 = 𝐸†k = 0.49 𝐽. 𝑘𝑔&# e) Kinetic energy is converted to internal energy
d) Kinetic energy (Ek), rotational energy (ER), potential energy (Ep), electric
energy (Eel), surface energy (Es), internal energy (U), solar energy. d) 𝐸†k = 𝑔ℎ From an energy balance that includes kinetic and potential energy:
+
e) „ + 𝐸†J + 𝐸†G ~
𝑑}𝑈
„F + 𝐸†J + 𝐸†G ~ + 𝑄 + 𝑊
= °}𝐻
• Kinetic energy: wind turbines Height equals vertical distance from reference to bottom of piston, plus half the 𝑑𝑡
FS#
• Potential energy: hydro-electric power generation hight of the gas, which is half of the original height
$.#
• Internal energy: energy balance in chemical processing plant, nuclear 𝐸†k = 9.81 𝑚𝑠 &"
× -0 + 0.5( ). 𝑚 = 𝐸†k = 0.24 𝐽. 𝑘𝑔 &#
If we assume no heat is leaving the water (all to temperature change), and for a closed
"
energy generation system (water fixed volume of water falling), with no work, this becomes:
• Electric energy: electricity generation *Compressing the piston will also increase the internal energy of the system, „ + 𝐸†J + 𝐸†G ~
𝑑}𝑈
= 0
• Rotational energy: turbines which we have ignored here. 𝑑𝑡
• Solar energy: solar panels for energy generation
e) Energy is typically a difference and not an absolute value. The values above This can be simplified to:
∆}𝑈„ + 𝐸†J + 𝐸†G~ = 0
are relative to the reference. If a different reference point was chosen, the
„ + 𝐸†J + 𝐸†G~ − }𝑈
}𝑈 „ + 𝐸†J + 𝐸†G ~ = 0
Example 4.51: Potential energy of a gas in a piston values above would have been different. However, the differences in the values A>G 9>AA>=
above would have been the same, even if the reference point was different. [See „ + 𝐸†G ~ − }𝑈
}𝑈 „ + 𝐸†J ~ = 0
A>G 9>AA>=
A piston (diameter 5 cm; height 10 cm) containing a certain volume of air is placed how this works for other forms of energy later in the section.]
2 m above the ground. Rearranging:
a) Define a suitable reference point. f) A reference is a vertical height above the ground. In the calculations, we used ∆𝑈9>AA>=→A>G = 𝐸†G A>G − 𝐸†J 9>AA>=
b) What is the potential energy of the gas in the piston? the term ‘z’, and not ‘h’ to calculate potential energy. ‘h’ is the distance from 𝐶C . ∆𝑇 = 𝐸†G A>G − 𝐸†J 9>AA>=
the ground, while ‘z’ is the distance from the reference point. 1 " 1
c) The piston is moved and placed on the floor. What is the potential energy of 𝐸†G A>G − 𝐸†J 9>AA>= 𝑔ℎA>G − 2 𝑣 9>AA>= 9.81.210A>G − 10"
A system is a defined mass or volume of material that we will be investigating. 2 9>AA>=
the gas at its new position? ∆𝑇 = = = = 0.48 ℃
𝐶C 𝐶C 4180
d) While still on the ground, what is the potential energy of the gas if the piston Therefore, the temperature at the bottom is 25.48°C
is compressed to halve its original height?* Example 4.52: Water from a waterfall
e) Discuss the importance of defining a reference point. f) Since energy is required for evaporation, it would cool the water.
Water falls from a 210 m waterfall.
f) How does a system differ from and a reference point?
a) What is the potential energy of the water at the top of the waterfall?
Example 4.53: Nitrogen cooled in a closed vessel The energy balance assuming no kinetic or potential energy: By interpolation: Specific volume = 0.98436 m3/kg
+
𝑑𝑈
= ° 𝐻F + 𝑄 + 𝑊
Nitrogen is stored in a 10 L flask at an initial pressure of 5.921 atm. Calculate the 𝑑𝑡
FS#
heat that is transferred when the nitrogen is cooled from 90 to 30°C. For a closed system:
Given: 𝑑𝑈
= 𝑄+𝑊
Cp [kJ/mol. °C] = 0.0290 + 0.2199 × 10-5T + 0.5723 × 10-8T2 – 2.871 × 10-12T3 𝑑𝑡
Example 4.56: Labelling enthalpy-pressure diagram for water
Therefore:
Solution: 𝑑𝑈 Given the enthalpy-temperature diagram for water below, label all unknowns.
𝑄= −𝑊
From an energy balance that assumes no kinetic or potential energy: 𝑑𝑡
+
𝑑𝑈
= °𝐻 +𝑄 +𝑊 Substituting values: c b
𝑑𝑡
FS# 𝐵𝑇𝑈 0.001285𝐵𝑇𝑈
𝑄 = 8 𝑙𝑏. 4 − 2334 𝑓𝑡. 𝑙𝑏- . = 29 𝐵𝑇𝑈 g (from e à c) h (from c à e)
𝑙𝑏= 1𝑓𝑡. 𝑙𝑏-
For a closed system with no work added, this becomes: e
𝑑𝑈 a
=𝑄 Example 4.55: Specific volume from enthalpy calculation d
Enthalpy
𝑑𝑡
A stream of water (15 kg/hr) is heated from 25°C to T2 at a constant pressure of 2 bar i (from là d) j (from d à l)
Rearranging and rewriting for the definition of U:
such that the final enthalpy is 335630.9 kJ/hr. Calculate the specific volume of the
𝑄 = 𝑛 ž 𝐶C 𝑑𝑇 l
water at the new temperature.
k f
Cv = Cp – R Given:
= (0.0290 + 0.2199 × 10-5T + 0.5723 × 10-8T2 – 2.871 × 10-12T3) – (8.314 × 10-3) Cp [kJ/kg] = (4.2 – 2 × 10-3) T[K]
Temperature
= (0.02069 + 0.2199 × 10-5T + 0.5723 × 10-8T2 – 2.871 × 10-12T3) kJ/mol.°C Hf0 = -285820 kJ/kmol
Molar mass of water: 18 kg/kmol
Therefore: Solution:
5$
Solution: a) Liquid water
𝑄 = 𝑛 ž (0.02069 + 0.2199 × 10&* 𝑇 + 0.5723 × 10&) 𝑇 " – 2.871 × 10&#" 𝑇 5 )𝑑𝑇 ?@ b) Superheated steam
&")*)"$?A3B, Jn
„-° =
Conversion: 𝐻 = −15878.9
'$
#)
?C Jo c) Condensation point
0.2199 × 10&* ?A3B,
𝑄 = 𝑛 ²0.02069(30 − 90) + (30" − 90") d) Freezing point
2
&)
W"
e) Boiling point
0.5723 × 10 „-° + 𝑚̇ × ž 𝐶G 𝑑𝑇
+ (305 − 905) 𝐻 = 𝑚̇𝐻 f) Solid
3
2.871 × 10&#"
W#
g) Heat of vaporisation
− (307 − 907 )³ 𝑘𝐽 𝑘𝑔 𝑘𝐽
4 335630.9 = 715 × −15878.9 ; h) Heat of condensation
ℎ𝑟 ℎ𝑟 𝑘𝑔
𝑄 = 𝑛[−1.24921 𝑘𝐽/𝑚𝑜𝑙] 𝑘𝑔 W"
𝑘𝐽 i) Heat of fusion
+ 60 × ž (4.2 − 2 × 10&5)𝑇 ( ) j) Heat of solidification
ℎ𝑟 "') 𝑘𝑔
From the ideal gas law: 𝑘𝐽 𝑘𝑔 𝑘𝐽 k) Absolute zero temperature
𝑃𝑉 (5.921𝑎𝑡𝑚)(10.0𝐿) 335630.9 ℎ𝑟 − 715 ℎ𝑟 × −15878.9 𝑘𝑔 ; 4.2 − 2 × 10&5 𝑘𝐽
𝑛= = = 1.9877 𝑚𝑜𝑙 = ¤Ž • × (𝑇" − 298)"¥ l) Melting point
𝑅𝑇 𝑎𝑡𝑚 15 2 ℎ𝑟
´0.08206𝐿. 𝑚𝑜𝑙. 𝐾µ (363𝐾)
4.2 − 2 × 10&5
38254.28 = Ž • × (𝑇" − 298)" Example 4.57: Steam table calculations involving interpolation
2
And: "
18225 = (𝑇" − 298)
𝑄 = 1.9877[−1.24921] = −2.483 𝑘𝐽 Using the steam tables, calculate the following:
135 = 𝑇" − 298
𝑇" = 433𝐾 1) Enthalpy of saturated water at 98 bar
Example 4.54: Heat loss in a closed system with work 𝑇" = 160 °𝐶 2) Enthalpy of saturated steam at 131°C
3) Enthalpy of saturated steam at 8 bar and 240°C
Calculate the total heat for a closed stationary system of 8lbs, when 2334 ft.lbf of work From steam table, at 2 bar (superheated steam): 4) Internal energy of saturated steam at 0.05 bar and 190°C
is done on it. The internal energy of the system increases by 4 Btu/lbm. 5) Enthalpy of vaporisation for saturated system at a pressure of 178 bar
150°C 160°C 200°C
Given: 1 ft.lbf = 0.001285 BTU
0.9602 1.081 [kJ/kg]
Solution:
Solution:
A liquid at 180°C and 1 002.7 kPa has an internal energy of 762.0 kJ/kg and a specific Including the unit conversion:
volume of 1.128 cm³/g. H1 = 762 + (1 002.7 × [(1.28/1003).1000]) = 763.13 kJ/kg
5 Fluid Flow
5.1 Introduction
In this section we will look at what affects fluid flow in a pipe. This may include fluid
density, fluid viscosity, pipe diameter, pipe roughness (materials of construction), pressure,
pipe height and other things.
Standard atmosphere is given as: 101.325 Pa (or 101.325 kPa = 1 013.25 hPa)
Absolute pressure is the pressure referenced against a perfect vacuum. It is equal to gauge
pressure plus atmospheric pressure.
Fluid
Solution:
Step 1:
Divide the manometer into various sections: Vertically into A and B and horizontally into
point 0, 1, 2, 3 and 4 as shown below.
4 Fluid
3
2
h
1
0
A B
Step 2:
Calculate the pressure on Side A, at point 0.
This equals the pressure of the manometer fluid from 0 à 1 plus the pressure from the fluid
from 1 à2 plus from 2 à 3 plus the pressure due to the flow of the fluid at 4.
where:
man refers to the density of the fluid in the manometer
fluid refers to the density of the fluid in the pipe
Step 3:
Calculate the pressure on Side B, at point 0.
Step 4:
The fluid through the manometer is NOT moving, therefore pressure at 0 is not changing
and the pressure is equal
Rearranging:
PN,5 − Pj,5 = ρG3, gh&→. − ρ'[?=D gh&→. 5.5
Therefore:
−∆P = gh&→. (ρG3, − ρ'[?=D ) 5.6
Or simply:
−∆P = gh∆ρ 5.7
5.3 Buoyancy
5.3.1 Introduction
Simply: buoyancy is the upward force exerted by a fluid (liquid or gas) on an object.
In a column of fluid, pressure increases with depth (P = rgh). This means that a column of
fluid experiences greater pressure at the bottom of the column than at the top. This
difference results in a nett force that tends to accelerate an object upwards. The magnitude
of that force is proportional to the difference in the pressure between the top and the
bottom of the column, and is also equivalent to the weight of the fluid that would otherwise
occupy the column, i.e. the displaced fluid. For this reason, an object whose density is
greater than that of the fluid in which it is submerged tends to sink. If the object is either
less dense than the liquid or is shaped appropriately (as in a boat), the force can keep the
object afloat.
Or:
FA = ρ> . V> . g 5.10
where:
Fg = Gravitational force
ro = density of the object
Vd = Volume of the object
Since the Archimedes principle states the buoyancy force equals the volume of displaced
water, the Bouyancy force can be written as:
F2 = ρ' . V~ . g 5.11
where:
Fb = Buoyancy force
rf = density of the fluid displaced
VD = Volume displaced
If the object is stationary, this means the forces are equal, therefore:
FA = F2 5.12
ρ> . V> . g = ρ' . V~ . g 5.13
Therefore:
ρ> . V> = ρ' . V~ 5.14
And:
•Q .oQ
V~ = 5.15
•z
Example 5.1:
Solution:
𝐹5 = 𝜌5 𝑉5 𝑔
𝐹b = 𝐹š = 𝜌b 𝑉b 𝑔
Since the iceberg floats, the forces are equal (Fi = Fw) and:
𝜌5 𝑉5 𝑔 = 𝜌b 𝑉b 𝑔
Therefore:
𝜌5 𝑉5 = 𝜌b 𝑉b
And:
›m „m ›i „i
𝑉5 = ›i
𝑉b = ›m
Example 5.2:
What maximum load (M) can a helium filled balloon with radius, R = 20 m and mass (including support
cables and basket), m = 196kg carry?
Solution:
𝐹W5T = 𝜌W5T 𝑉W5T 𝑔 Buoyancy force
𝐹Zi = 𝜌Zi 𝑉Zi 𝑔
Now:
𝑉 = 4𝜋𝑅, /3
From Archimedes’:
𝐹W5T = 𝐹Zi + 𝑚𝑔 + 𝑀𝑔
𝐹W5T 𝐹Zi
𝑀= − −𝑚
𝑔 𝑔
Substituting values:
4
𝑀 = 0 𝜋7 . 12, . (1.25 − 0.160) − 196
3
𝑀 = 7690 𝑘𝑔
This mass would result in the balloon forces balancing and thus not moving. Anything lighter would mean
the balloon rises.
The principle can be applied to hydraulic presses. Since the pressure in the fluid is constant,
a pressure exerted at one end of the press results in an equal pressure at the other end.
However, if the two areas are of different sizes, the resulting force at the end with the
larger area is greater.
A
B
A force exerted at point A results in a pressure throughout the fluid. This pressure exerts
a force upwards at point B. Since the area at B is larger than A, the force upwards is larger.
Kevin Harding (kevin.harding@wits.ac.za) 141
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Fluid Flow – Terminal Velocity
.GA
V@ = ‹•Nk 5.16
}
where:
Vt = terminal velocity (m/s)
m = mas of object (kg)
g = acceleration due to gravity (m/s2)
r = density of object (kg/m3)
A = projected area (m2)
Cd = coefficient of drag
.€œ
CD = 5.17
•• ! N
where:
Cd = coefficient of drag
Fd = drag force (N)
v = speed of the object relative to the fluid (m/s)
A = reference area (m2)
Laminar Flow
Turbulent Flow:
Particles move and knock against each other. This movement results in a flat velocity
profile. [Chaotic eddies, vortices, and other flow instabilities]
Turbulent Flow
Definition:
Reynolds Number is defined by the formula:
•|~
Re = 5.18
‚
Where:
r = density (kg/m3)
v = velocity (m/s)
D = pipe diameter (m)
µ = viscosity (kg.m-1.s-1) (although commonly given in cP)
If:
- Re < 2300: Laminar Flow
- Re > 104: Turbulent flow
- Re between 2300 and 104: barrier between laminar and turbulent flow.
5.6.3 Viscosity
Viscosity is a measure of the resistance of a fluid, i.e., the ‘thickness’ of the fluid. A fluid
with a high viscosity will not splash as much as one with a lower viscosity and may be able
to form peaks. A fluid with a low viscosity flows more easily.
As temperature increases, the viscosity of liquids decreases, e.g., syrup in cold temperatures
does not flow out a bottle as easily. However, as temperature increases, the viscosity of
gases increases.
Example:
Using the attached nomograph, calculate the viscosity of Helium at:
a.) 20°C b.) 300°C
Solution:
For Helium:
x: 10.9 y: 20.5
Figure 5.2: Viscosity for various chemicals as a function of temperature (unknown source)
Kinematic viscosity (n) is the ratio of dynamic viscosity to the fluid density.
n = µ/r
Units: 1 cSt = 1 mm2.s-1 = 10-6m2.s-1
7‚On &.7‚On
∆P = = 5.19
TCn TDn
where:
DP = pressure (Pa)
µ = viscosity (Pa.s)
L = length (m)
Q = volumetric flow (m3/s)
r = radius (m)
D = diameter (m)
where:
DE = change in energy of the system
E = the flow of energy into and out of the system
Q = Heat
W = All forms of work added to the system
So:
Kevin Harding (kevin.harding@wits.ac.za) 149
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Fluid Flow – Bernoulli’s Equation
And:
∆U = U=, − U>?@ + Q̇ + Ẇ 5.26
Therefore:
D]
= /H=, − W'[>d,=, 1 − /H>?@ − W'[>d,>?@ 1 + Q̇ + Ẇ 5.27
D@
Rearranging:
D
/U + El + EM 1 = /Ḣ=, + El,=, + EM,=, 1 − /Ḣ>?@ + El,>?@ + EM,>?@ 1 + Q̇ + Ẇ −
D@
/Ẇ'[>d,=, + Ẇ'[>d,>?@ 1 5.29
D
/U + El + EM 1 = /Ḣ=, + El,=, + EM,=, 1 − /Ḣ>?@ + El,>?@ + EM,>?@ 1 + Q̇ + Ẇ 5.30
D@
Where:
W = all forms of work
Therefore:
D
/U + El + EM 1 = /Ḣ=, + El,=, + EM,=, 1 − /Ḣ>?@ + El,>?@ + EM,>?@ 1 + Q̇ + Ẇ + Ẇƒ@4BC
D@
5.31
D]
= /Ḣ=, − Ḣ>?@ 1 + Q̇ + Ẇ 5.32
D@
D]
= ∑,=I& Ḣ= + Q̇ + Ẇ 5.33
D@
where:
i = stream (streams flowing in are positive; streams flowing out are negative)
n = number of streams in the system
At Steady State:
D
=0 5.34
D@
Therefore:
D]
= ∑,=I& Ḣ= + Q̇ + Ẇ 5.35
D@
D]
= ∑,=I& Ḣ= + Q̇ + Ẇ 5.37
D@
0 = Q̇ + Ẇ 5.38
NOTE: The term ‘momentum balance’ is sometimes used in the concepts we are dealing
with here but is a different idea which should not be confused with what we are doing
here.
D
/U + El + EM 1 = /Ḣ=, + El,=, + EM,=, 1 − /Ḣ>?@ + El,>?@ + EM,>?@ 1 + Q̇ + Ẇ + Ẇƒ@4BC
D@
5.40
For thermal systems, Ek and Ep are approximately zero as they don’t affect the temperature
of the system. However, these terms do affect the velocity.
In a similar way, enthalpy (H) does not affect the velocity (nor does Q, Welec, Wsound and
other forms of work), even though they can/do affect the temperature.
D
/U + El + EM 1 = /El,=, + EM,=, 1 − /El,>?@ + EM,>?@ 1 + ẆF + Ẇ>@4BC 5.42
D@
At Steady State:
D
/U + El + EM 1 = 0 5.43
D@
Therefore:
D
/U + El + EM 1 = /El,=, + EM,=, 1 − /El,>?@ + EM,>?@ 1 + ẆF + Ẇ>@4BC 5.44
D@
Other Work:
1) Friction/Friction loss (F)
Can also be written as: F, F$, f, fl, ff
q,>?@
dP
ẆFCB| = m ~
q,=, ρ
With:
El = mgz [3]
&
EM = (mv . ) [4]
.
where:
m = mass (kg)
g = gravitational constant (m2/s2)
z = height (sometimes written as h) (m)
v = velocity (m/s)
NOTE:
- Friction and energy change due to pressure change is negative since energy is lost;
- For laminar flow, velocity at the centre of the pipe is greatest;
- Friction reduces the flow at the pipe walls;
- Re indicates whether flow is laminar or turbulent ;
- Flow is almost never in a straight line; it is more circular motion and eddies are formed; and
- If it is not stated or cannot otherwise be determined, we assume turbulent flow.
Final form of the energy balance: (substitute [3] and [4] into [2] and re-arrange)
& & q,>?@ Dq
Jmgz + . (mv . )K + F = Jmgz + . (mv . )K − m ∫q,=, + ẆF [5]
•
𝐨𝐮𝐭 𝐢𝐧
Units: J or J/s
|! ‰ F |! ‰ q,>?@ Dq i̇k
%gz + & + G = %gz + & − ∫q,=, + [6]
. >?@ . =, • G
q,>?@
dP 1 q,>?@ P>?@ − P=, P>?@ P=,
~ = ~ dP = = −
q,=, ρ ρ q,=, ρ ρ ρ
Rearranging:
P v. ∆P' P v. ẆF
p + gz + r + = p + gz + r +
ρ 2 >?@ ρ ρ 2 =, m
[7]
Units: J/kg
Units: Pa
Defining:
𝑊̇Œ . 𝜌
∆𝑃u‹Š; =
𝑚
ρv . ρv .
pP + ρgz + r + ∆P' = pP + ρgz + r + ∆P[='@
2 >?@ 2 =,
•|!
∆ %P + ρgz + & = ∆P[='@ − ∆P' [9]
.
ρv . ρv .
pP + ρgz + r = pP + ρgz + r
2 >?@ 2 =,
•• !
𝑃 + 𝜌𝑔𝑧 + = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 [10]
.
Where:
ik .•
∆P[='@ = and
G
F. ρ
∆P' =
m
Therefore:
ρv . F. ρ ρv . WF . ρ
pP + ρgz + r + = pP + ρgz + r +
2 >?@ m 2 =, m
Dividing by 𝜌𝑔:
q |! F q |! i
%•A + z + .A& + GA = %•A + z + .A& + GAk [11]
>?@ =,
Units: m
Summary:
J/kg.s
•|! •|!
%P + ρgz + & + ∆P' = %P + ρgz + & + ∆P[='@ Pa
. >?@ . =,
•|!
P + ρgz + = constant (ONLY when no friction, turbulent flow, no pumping)
.
q |! F q |! i
%•A + z + .A& + GA = %•A + z + .A& + GAk m
>?@ =,
q |! ∆qz q |! i
%•A + z + .A& + = %•A + z + .A& + GAk m
>?@ •A =,
5.7.5 Siphoning
A siphon is a device which typically involves a tube connecting two containers. Liquid in
the container of higher height will flow through the tube into the other container.
What is the linear velocity of a liquid at a tank surface when liquid is being removed in a thin pipe?
Mtank = Mpipe
Qtank.r = Qpipe. r
Where
Q = volumetric flow
Now:
Q = Cross sectional area × linear velocity
XSAtank.vtank = XSApipe.vpipe
Rearranging:
vtank = (XSApipe/ XSAtank).vpipe )
If we assume XSAtank >> XSApipe
(XSApipe/ XSAtank). à 0
Therefore: vtank = 0
4m
1m
2m
3
P1: Patm (this is why we needed to know it was at the coast – else it needed to be given)
r: density of water = 1000 kg/m3
g = 9.81 m2/s2
z = 0 (set a reference height at this point)
v1 = 0
Now at point 3:
56 )
Again, 𝑃 + 𝜌𝑔𝑧 + $
= 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
P3 = Patm
r: density of water = 1000 kg/m3
g = 9.81 m/s2
z = -1m
v3 = unknown
Constant = Value from point 1
At point 2:
56 )
Again, 𝑃 + 𝜌𝑔𝑧 + = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
$
ρv . α ρv . α
pP + ρgz + r + ∆P' = pP + ρgz + r + ∆P[='@
2 >?@ 2 =,
where:
f = friction factor (Fanning Friction factor) – Obtained from a graph.
L = Length (m)
D = Diameter (m)
r = density (kg/m3)
v = velocity (m/s)
Unit analysis:
𝑁 [𝑚] 𝑘𝑔 𝑚 (
= [𝑓] ? @{ |
𝑚 ( [𝑚] 𝑚, 𝑠
𝑘𝑔 𝑘𝑔 1
(
= [𝑓] ? @ ? ( @
𝑚. 𝑠 𝑚 𝑠
Therefore f à dimensionless
NOTE: Depending on the graph used (type of friction factor), the formula will change; e.g.
Fanning friction factor vs. Moody and others.
Where:
O
η = 4f ~
5.8.1 Finding f
To find the friction (f), we use a Moody friction factor graph. This gives the friction factor as
a function of Reynolds number, for different roughness values of pipes.
Where:
•|~
Re = [dimensionless]
‚
e = surface roughness [Units: length; typically, in mm]
D = pipe diameter [Units: length (same as for e)]
e will either be given (tests, exams etc.) or can be obtained from Perry or other texts, such as
the one in Figure 5.5.
Therefore:
16 L v .
∆P' = 4 ρ
Re D 2
16 L v .
∆P' = 4 ρ
ρvD/µ D 2
Lµv
∆P' = 32 .
D
O‚
∆P' = %32 ~!& v ∆P' ∝ v
Example:
Water flows through a commercial steel pipe at 5m/s. If the pipe is 10 cm in diameter,
what is the pressure drop due to friction over 100m of pipe?
Solution:
𝐿 𝑣$
∆𝑃7 = 4𝑓 𝜌
𝐷 2
e = 0.00015 ft
e = 0.004572 cm
Re = rvD/µ
= (1000)(5)(0.1)/(1 × 10-3)
= 5 × 105
𝐿 𝑣$
∆𝑃7 = 4𝑓 𝜌
𝐷 2
= 4(0.0042)(100/0.1)(1000)(5^2)/2
= 210000 Pa
= 210 kPa
The friction for a pipe (length and bends/valves/constrictions) can be calculated from the
general equation:
L v.
∆P' = 4f ρ
D 2
Oy |!
Where each term %4f ρ . & is pipe length and each bend, valve, and other fitting in the
~
system.
Why is this not factorized to pull the 4f and D terms to the left?
Reason:
Diameter could change.
Depending on the method to obtain the friction factor: 4fL/D = h
If the method of obtaining the friction factor does not involve using h, we can further factorize
this to:
v. LB LB LB LB
∆P' = ρ . 4f dJ K + J K + … + J K +J K e
2 D & D . D ,#& D ,
LB LB LB LB
∆P' = 2f. ρ. v . dJ K + J K + … + J K +J K e
D & D . D ,#& D ,
Le/D values for fittings are given in tables. Different researchers may have different values
(or formula) to calculate these.
Example:
Water at 90oF flows at 1000 lb/min through the system shown below. What is the friction term 𝐹h for the system?
Given:
e = 0.00015 ft or 0.00005 m
equivalent lengths:
elbow 30D
valve 12D
Pipe ID 2’’
Pipe ID 4’’ 100 ft
50 ft
10 ft
50 ft
Flow
25 ft 25 ft 50 ft
Solution:
𝐿8 𝑣 $
∆𝑃7 = 4𝑓 𝜌
𝐷 2
2 inch section:
1) Pipe length: 50 + 10 + 25+ 25 = 110ft
2) Pipe bends: 2
3) Number of valves: 1
4 inch section:
1) Pipe length: 50 + 50+ 100 = 200ft
𝐿8 𝑣 $ 𝐿8 𝑣 $ 𝐿8 𝑣 $
∆𝑃7 = 4𝑓 𝜌 + 4𝑓 𝜌 + 4𝑓 𝜌 + … . 𝑓𝑜𝑟 𝑒𝑎𝑐ℎ 𝑜𝑓 𝑡ℎ𝑒 𝑠𝑒𝑐𝑡𝑖𝑜𝑛𝑠 𝑎𝑠 𝑎𝑏𝑜𝑣𝑒.
𝐷 2 𝐷 2 𝐷 2
6)
Fanning: ∆𝑃7 = 4𝑓 @+ 𝜌
/
$
6)
Darcy: ∆𝑃7 = 𝑓@
/+
@
𝜌
$
(sometimes also called the Darcy–Weisbach friction factor or Moody friction factor)
Therefore:
𝑓@ = 4𝑓
The Darcy friction factor is simply four times larger than the Fanning friction factor.
Most Moody charts indicate on them which of the friction factor equations need to be used.
Example:
a.) Compare Fanning friction factor values obtained using the Moody chart, the Moody approximation and the
Haaland approximation for a smooth pipes at Reynolds numbers of 1 × 105, 1 × 106, 1 × 107 and 1 × 108
b.) Assuming the value obtained from the chart is correct (not a good assumption), calculate the absolute
percentage error in each approximation.
c.) Which approximation method is “better”?
Solution:
a.)
Re Chart Moody approx. Haaland approx.
1 × 105 0.0045 0.00434 0.00446
1 × 106 0.0028 0.00275 0.00290
1 × 107 0.002 0.00201 0.00203
1 × 108 0.0015 0.00167 0.00150
b.)
Re Moody approx. Haaland approx.
1 × 105 3.75 0.98
1 × 106 1.82 3.34
1 × 107 0.66 1.61
1 × 108 10.25 0.31
Average 4.12 1.56
c.) If the definition of “better” is ‘more accurate’, the Haaland method should be used. If one wishes to have a faster
method, the Moody or chart may be more appropriate.
Additional Reading
Coulson, JM, Richardson, JF, Backhurst, JR, Harker, JH, 1999. Coulson & Richardson’s Chemical
Engineering, Volume 1, 6th Edition, Elsevier.
Munson, BR, Okiishi, TH, Huebsch, WW, Rothmayer, AP, 2012. Fluid Mechanics, 7th Edition
(SI Version), Wiley.Perry
Perry, RH, and Green, DW, 1984. Perry's Chemical Engineers Handbook, 6th Edition.
Skogestad, S, 2009. Chemical and Energy Process Engineering, CRC Press.
Welty, JR, Wicks, CE, Wilson, RE, Rorrer, GL. 2008. Fundamentals of Momentum, Heat, and
Mass Transfer, 5th Edition, Wiley.
Young, DF, Munson, BR, Okiishi, TH, Huebsch, WW, 2012. Introduction to Fluid Mechanics,
5th Edition (SI Version), Wiley.
Problems
5.1 For the Figures 1 à 4 below, answer the following questions:
a.) The diagram in Fig. 1 shows a water pipeline. For the section 1 to 2, calculate the pressure difference ΔP and the friction loss yFh{.
b.) Calculate ΔP and Fh for the system shown in Fig. 2.
c.) Repeat (b) for the case where the liquid flowing in the pipe is an organic solvent of density 1200 kg.m-3.
d.) Calculate ΔP and Fh for the system shown in Fig. 3.
e.) What would the manometer reading be in the system of Fig. 3 if the flow in the pipe were downwards? Assume the value of Fh to be unchanged.
f.) Determine the friction loss Fh for the liquid flowing in the system shown in Fig. 4.
37 kPa 10 kPa
Water
1 2 10cm Hg
Fig. 1 Fig. 2
Water 10 cm
3m
0.5 m
15cm Hg
Fig. 3 Fig. 4
5.2 Water flows into a pipe which splits into two at a T-piece. Under what configuration/conditions will the flow split equally into the two arms? List at least 4 things.
5.3 Brine of s.g. 1.18 is to be pumped at 110 kg.min-1 through 5 cm piping in the system shown below. The friction term Fh is calculated as 0.418 kJ.kg-1. What is the power required to
do the pumping?
9m
1m
5.4 The figure below shows a water siphon. The siphon tube is initially filled with water before the lower end is opened to atmosphere.
1.5 m
20 mm siphon tube
2m
5.5 Brine of specific gravity 1.18 and viscosity 2.5 cp is to be pumped through 100 m of smooth copper tubing of 0.995” ID at 110 kg/min. What is the energy loss due to friction
measured in:
a) J/kg b) head (m) of liquid c) kPa
5.6 Water at 90oF flows at 1000 lb/min through the system shown below. What is the friction term ∆𝑃7 for the system?
Given:
e = 0.00015 ft or 0.00005 m
equivalent lengths: elbow 30D, valve 12D
100 ft
Pipe ID 2”
Pipe ID 4”
50 ft
50 ft
Flow
10 ft
25 ft 25 ft 50 ft
5.7 A pump takes water at 15°C from a large reservoir and delivers it to the bottom of an open elevated tank 8 m above the reservoir surface through an 8 cm ID pipe. The inlet to the
pump is located 3 m below the water surface, and the water level in the tank is constant at 11 m above the reservoir surface. The pump delivers 400 L/min. Assume a total
equivalent length for the pipeline of 25 m. If the pump and motor have an overall efficiency of 55%, determine the electrical power required to do the pumping.
5.8 Calculate the power delivered to an electrical generator by the water turbine shown in the figure below if the efficiency of the turbine is 80%. The turbine inlet has an ID of 20 cm
and the outlet has an ID of 30 cm. (Neglect friction.)
50 m
10 cm ID
10 m
5.9 The figure below shows an elevated tank connected to a pipeline. The system contains water at 180oF. What must be the height of the water surface in the tank to produce a flow
of 100 gal/min through the line? The water sprays into the air at the right hand end of the 2” pipe.
7.48 gal = 1 ft3 μ = 0.347 cp ρ = 60.6 lb.ft-3 elbow = 20D
Pipe ID 2”
H Pipe ID 4”
9’ 8
”
20’
50’
10
’
25’
5.10 An edible oils factory (as shown below in aerial view – not to scale) currently stores oil in Tank A before moving it to a second storage area (Tank B) at 250kg/min. Line A-B, which
has a 60 mm portion including two 45° bends, as well as a 40 mm portion needs to be modified in order to make way for administration buildings. As such, Tank C will be
constructed and the oil will flow, at the same flowrate, in a more direct route to Tank B. Valve V1, which will remain in place, will be closed so as to block flow when all old piping
and tanks are removed.
a.) Calculate the energy loss due to friction (in J/kg) of the current system (Tank A to Tank B).
b.) Compare this to the energy loss of the new system (Tank C to Tank B). Explain why the friction is higher or lower in this pipe.
c.) How will the following two suggested modifications to the system affect the flow to tank B:
i) The oil is heated up before being pumped to Tank B.
ii) The pipe diameter is increased
Additional information:
Take e = 0.00005
Equivalent lengths:
Elbow: 30D
Valve: 12D
45° bends: 15D
Straight through Tee: 20D
Right angle Tee: 70D
Sudden contraction: 22D
25 m 65 m
30 m
Tank C
V2 V3
20 m
V1
15 m
Flow
Tank B
20 m
35 m 35 m
Pipe ID 40 mm
Tank A
Pipe ID 60 mm
Flow
6 Pumping
6.1 Introduction
A pump is a device used to move fluids from one location to another
Pressure: Pressure at the outlet of the pump is higher than the inlet.
Velocity:
From the conservation of mass:
Mass in = mass out
Min = Mout
Therefore:
Qin. r = Qout.r where Q = volumetric flowrate
For an incompressible fluid, r is constant and:
Qin = Qout
Now:
Q = XSA × linear velocity
XSAin.vin = XSAout.vout
If the diameter of the pipe in = diameter of the pipe out (XSAin = XSAout) we are left with:
vin = vout
where:
WF . ρ
∆P[='@ =
m
6.2 Pumps
There are four basic types of pumps:
- Forward Action
- Reverse Action
Flow
Time
P v. F P v. WF
p + z+ r + =p + z+ r +
ρg 2g >?@ mg ρg 2g =, mg
Units: m (length)
|!
velocity head [a]
.A
q
pressure head [b]
•A
z potential head [c]
F
head loss due to friction [d]
GA
As flow rate increases, head decreases. The type of impellor fitted to the pump affects the
curve.
H (m)
Q (m3/s)
If there is zero flow, i.e., no pumping and Ws = 0, there can still be a head.
Velocity = 0
Friction = zero
Pressure change = zero
Therefore, the head is the potential head.
Kevin Harding (kevin.harding@wits.ac.za) 178
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Pumping – Geometrically Similar Pumps
Example:
Water is to be pumped in the system shown below:
2m
7m
1m
1.5 m
Table 1 gives the head loss due to friction in the pipeline at various flow rates. Table 2 gives Q-H data for a certain
10 cm diameter pump running at 3000 rpm.
a) What flow rate will be obtained using this pump in the system?
The manufacturers curve describes what is available, while the system curve describes what
is wanted.
𝑁𝑄&/.
𝑁Œ =
(𝑔𝐻)%/5
This gives rise to the Affinity Laws which give correlations between geometrically similar
pumps:
where:
1: Pump 1
2: Pump 2
Q: Volumetric flowrate (m3/s)
N: Impellor speed (rpm)
D: Impellor diameter (m)
H: Head (m)
P: Pump power (kW)
These laws assume that the pump efficiency remains constant. The laws work well for
constant diameter pumps but are less accurate for constant speed pumps.
Example:
A pump impellor is changed while the pump speed is kept constant.
Calculate the final flow rate and head. The initial flow was 100 gpm
at a head of 100 ft and 5 bph. The initial and final impellor
diameters were 8 and 6 inches respectively
Solution:
The final flow rate:
q2 = q1(D2/D1) = (100 gpm)(6 inches/8 inches) = 75 gpm
Example:
The speed of a pump is changed while the impeller size is kept
constant. Calculate the final flow rate and head. The initial flow
was 100 gpm at a head of 100 ft and 5 bph. The initial and final
speeds were 1750 rpm and 3500 rpm respectively
Solution:
The final flow rate:
q2 = q1(N2/N1) = (100 gpm)(3500 rpm/1750 rpm) = 200 gpm
System Curve
H (m)
Q (m3/s)
Doubling the head of the pump curves does not result in a doubling the head for the system.
For an infinite number of pumps in series, the flow rate tends to a constant value,
irrespective of the system curve.
System Curve
H (m)
Q (m3/s)
Doubling the flow of the pump curves does not result in a doubling the flow for the system.
For an infinite number of pumps in parallel, the head tends to a constant value, irrespective
of the system curve.
When two different pumps are added in series or in parallel, the resulting pump curve is as
below:
Pump curve 2
H (m)
Pump curve 1
3
Q (m /s)
Figure 6.11: Pump curve for different sized pumps; added in series
Pump curve 2
NEW Pump curve
H (m)
Pump curve 1
3
Q (m /s)
Figure 6.12: Pump curve for different sized pumps; added in parallel
Example:
Water is to be pumped in the system shown below:
2m
7m
1m
1.5 m
Table 1 gives the head loss due to friction in the pipeline at various flow rates. Table 2 gives Q-H data for a certain 10cm
diameter pump running at 3000 rpm.
a) What flow rate will be obtained using this pump in the system? (as from previous example)
b) At what speed would the pump need to run to give half this flow rate?
c) What flow rate will be obtained if the pump is replaced by two pumps in series if:
a. Both pumps are 10 cm in diameter running at 3000 rpm;
b. Both pumps are geometrically similar to the pump on Table 2, but are both 15cm in diameter and
run at 2500 rpm?
c. One pump is 10 cm in diameter running at 3000 rpm and one is 15 cm in diameter and runs at
2500 rpm
d) What flow rate would be obtained using two 10cm/3000 rpm pumps in parallel?
Table 2: Head loss for a 10cm diameter pump running at 3000 rpm
Q (m3/s) H (m)
0 25
0.01 24
0.02 22
0.04 18
0.06 12
0.08 2
Example:
Water is to be pumped through a smooth pipe from one reservoir to another on top of a hill as shown. What is the maximum
flow rate which can be obtained given the pump curve below (see solution) if the pipe ID is 16 cm and the equivalent length of
the pipe is 100 m?
40
m
Solution: We need to draw the system curve onto the pump curve and obtain intersection:
From the mechanical energy balance:
𝑃 𝑣$ ∆𝑃?+7H 𝑃 𝑣$ 𝑊=
o + 𝑧+ s + =o + 𝑧+ s +
𝜌𝑔 2𝑔 FGH 𝜌𝑔 𝜌𝑔 2𝑔 +, 𝑚𝑔
Assuming velocity at inlet & outlet are equal (same diameter pipe), pressures are equal & the reference height is the pump:
Selecting various flows, manipulating as needed and determining the head from this: System curve
Flow Re Fanning FF* Hpump
l/sec m/s - - m
0 0 0 - 40
2 0.099 15915.49431 0.00682258 40.009
4 0.199 31830.98862 0.00574756 40.029
6 0.298 47746.48293 0.00523253 40.059
8 0.398 63661.97724 0.00490797 40.099
10 0.497 79577.47155 0.00467657 40.147
12 0.597 95492.96586 0.00449944 40.204
14 0.696 111408.4602 0.00435742 40.269
16 0.796 127323.9545 0.00423976 40.342
40
30
Meters
20
10
0
0 5 10 15 20
Litres/second
Kevin Harding
Therefore, (kevin.harding@wits.ac.za)
the flowrate is approx. equal to 6 litres per second 186
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Pumping – Cavitation
6.8 Cavitation
In a pump, if you no longer have a liquid due to phase change, i.e., liquid changes to a gas,
you will no longer have flow. Pvap is the lowest pressure (at a particular temperature) at
which flow is possible.
If pressure is less than the vapour pressure (Pvap) at pump temperature, vaporisation will
occur. This results in bubble formation which is known as cavitation. Gasses are
compressible so pumping becomes more difficult. Further, this results in mechanical
damage to the pump.
When the bubbles formed from cavitation collapse, there is a high localized mechanical
energy release which causes impellor corrosion and pitting.
E F
C
A B D
Velocity
A B C D E F
Because of the conservation of mass, and incompressible fluid, the velocity along the system
remains constant.
Pressure
A B C D E F
NOTE: Pressure decreases along the length of a pipe. This is due to friction. If there is no
friction; there is no pressure drop along a pipe of constant height.
Patm + rgh
Pressure
Patm
A B C D E F
Patm + rgh
NPSHr
Pressure
NPSHa
Patm
Pvap
A B C D E F
where:
vi = velocity at the pump inlet (m/s)
Pin = pressure at the pump inlet (Pa)
Pvap = vapour pressure (Pa)
NOW:
From the mechanical energy balance for the section of pipe between the system inlet (in)
and the pump inlet (out):
P v. F P v. WF
p + z+ r + =p + z+ r +
ρg 2g >?@ mg ρg 2g =, mg
Rearranging and knowing that Ws = 0; since we are looking at the section before the pump:
P v. P v.
p + z + r = p + z + r + h'
ρg 2g =, ρg 2g >?@
Assuming the reference height is the pump AND the velocity of the incoming water is zero
(system is from an open tank)
Therefore:
vin = 0
Pin = Patm
P v. P v.
p + z + r = p + z + r + h'
ρg 2g =, ρg 2g >?@
P3@G P v.
+ z=, = p + r + h'
ρg ρg 2g >?@
Therefore re-defining ‘in’ as the pump inlet and ‘tank’ as the inlet to the system:
P3@G P v.
+ z@3,M = p + r + h'
ρg ρg 2g =,
Rearranging:
qK* |K* ! qžS€
+ = + z@3,M − h' [2]
•A .A •A
Therefore:
P|3l P3@G
NPSH3 + = + z@3,M − h'
ρg ρg
P3@G P|3l
NPSH3 = J − K + z@3,M − h'
ρg ρg
Note:
The reference height was the pump. It the tank surface is below the pump, z is negative.
Additional Reading
Coulson, JM, Richardson, JF, Backhurst, JR, Harker, JH, 1999. Coulson & Richardson’s
Chemical Engineering, Volume 1, 6th Edition, Elsevier.
Perry, RH, and Green, DW, 1984. Perry's Chemical Engineers Handbook, 6th Edition.
Welty, JR, Wicks, CE, Wilson, RE, Rorrer, GL. 2008. Fundamentals of Momentum, Heat,
and Mass Transfer, 5th Edition, Wiley.
Problems
6.1 A tank farm has several tanks for the storage of diesel, petrol and other fuels. Large road tankers deliver the fuel which is then stored until either smaller tankers, or commercial
farmers collect the fuel. In order to reduce the number of pumps on site, it is required that you design a system of pipes whereby, by only closing a few valves, fuel can be pumped
either from the road the road tanker to tank, or from the tank to road tanker, using a single pump. NOTE: A pump can only pump in one direction.
6.2 Water is to be pumped through the system shown, using a 4 ½“ centrifugal pump connected to a 2880 rpm motor. Pump characteristics for this pump and for a geometrically
similar 5” pump are provided. What flow rate will result and what will be the electricity consumption of a 95% efficient motor? (Obtain the characteristic in two ways and
compare.)
Pipe ID 2” length
100 ft
10 ft
12 ft
6.3 Water is to be pumped from a farm dam to a swimming pool beside the farmhouse, which is on top of a hill as shown. The farmer has available one 2880 rpm motor, one 3450 rpm
motor, one 4 ½“ pump and one 5” pump. What is the maximum flow rate which he can obtain?
Pipe ID 3” length
250 ft
120 ft
6.4 A pump takes water at 15°C from a large reservoir and delivers it to the bottom of an open elevated tank 8 m above the reservoir surface through an 8 cm ID pipe. The inlet to the
pump is located 3 m below the water surface, and the water level in the tank is constant at 11 m above the reservoir surface. The pump delivers 400 L/min. Assume a total
equivalent length for the pipeline of 25 m. If the pump and motor have an overall efficiency of 55%, determine the electrical power required to do the pumping. Assume a Fanning
friction factor of 0.006. [Include a diagram of the system in your solution.]
6.5 Explain your understanding of a pump manufacturer’s curve. Show using a diagram and explain for each of the available head, efficiency and power curves/lines.
a.) Show in a diagram and explain how 2 identical pumps in parallel change the overall pump curve.
b.) Show in a diagram and explain how 2 different pumps in series change the overall pump curve.
c.) What is NPSH? Include in your explanation the difference between NPSH available vs. NPSH required.
6.6 2.27 m3/h of water at 320K is pumped in a 40 mm ID pipe a distance of 150 m in a horizontal direction and then up through a vertical height of 10 m. In the pipe there is a control
valve for which the friction loss may be taken as equivalent to 200 pipe diameters and also other pipe fittings equivalent to 60 pipe diameters. The pipe also has a heat exchanger
across which there is an additional loss in head of 1.5 m of water.
a.) Ignoring the velocity head, calculate the total head loss (in metres) in the system.
b.) From the calculation in (a), what power (W) must be supplied to the pump if it is 60 per cent efficient in order to deliver the required flow?
c.) What is the percentage error in the head calculation (as in part a above) due to the assumption of ignoring the velocity head?
d.) The mass flowrate needs to double in the system. What do you expect to happen to the head in the pipe? Explain in terms of pressure, velocity, frictional and total head
loss.
Additional information:
At 320K, water has a viscosity of 0.65 x 10-3 Ns/m2
The Fanning friction factor which can be used in this system is 0.008.
The roughness of the pipe is 0.2 m
6.7 A junior engineer has been asked to increase both the flowrate and head of a certain line by adding additional pumps to the current system. In a warehouse, he finds two additional
pumps of the same size as the currently installed one. The pump curve (for a single pump) and system curve data are given below.
a.) He decides to add to the current system a second pump in parallel. What is the new flowrate and head? Show graphically.
b.) It is suggested that he rather adds the second pump in series. What is the flowrate and head of this system? Show graphically.
c.) Give one reason why it might be better to add the pumps in parallel rather than series.
d.) The engineer decides to use all available pumps and suggests the setups below. Rank these in order of highest to lowest flowrate as well as the highest to lowest head.
Design 1)
Design 2)
Design 3)
Design 4)
6.8 Water at 20°C is pumped from Tank A to tank B as shown in Figure 1 (not to scale) through a 10 cm pipe. On its way from A to B, the water passes through a gate valve, pump and
gate valve respectively. All bends in the pipe are 90° elbows. Roughness in the pipe is given as 248.8 µm.
a.) Assuming diameter of tank A is much larger than the diameter of the pipe leaving it, prove that the velocity in the tank is approximately zero.
b.) Instead of using a Moody chart, friction can be calculated using the Haaland equation (turbulent region only). Calculate the friction factor assuming water is flowing at
5000 kg/hr.
',
# B.E 8 -
Haaland: = −3.6𝑙𝑜𝑔#D 5 M8 + H%.N@I 6
L7
c.) Using the friction factor above (i.e. assume it is constant at all flow rates), derive an expression for the system curve (Q-H) of the tank setup given. Ignore the effects of
sudden enlargement and contraction.
d.) Given the attached pump curve, what is the head and flow of the system with the pump supplied? Draw on the given diagram to indicate how you obtained this.
e.) If the pump is replaced by 2 identical pumps in series, what is the new flowrate and head achieved? Include a new pump curve on the given diagram.
f.) If the pump is replaced by 2 identical pumps in parallel, what is the new flowrate and head achieved? Include a new pump curve on the given diagram.
Tank A
Pump curve
0.25
15 m
0.2
10 m 5m 5m
0.15
5m
H (m)
0.1
25 m
0.05 35 m
0 Tank B
0 0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 0.0045
Q (m3/s)
Aerial view of Tank A and B to be used in Question 5. (NOTE: There is no change in horizontal height from A to B.)
Solution:
a.) Mass in = Mass out
[(Volumetric flow rate). Density]in = [(Volumetric flow rate). Density]out
Since density is constant: [Volumetric flow rate]in = [Volumetric flow rate]out
[Linear velocity × Cross sectional area]in = [Linear velocity × Cross sectional area]out
[Linear velocity]in = [Linear velocity × Cross sectional area]out/[Cross sectional area]in
But since the tank diameter is MUCH larger than the pipe diameter, XSAout/XSAin à 0
And: [Linear velocity]in à 0
b.) Calculating required values:
Volumetric flowrate: 5m3/hr = 0.0013889m3/s
XSA pipe: 0007854 m2
Velocity in pipe: 0.17683 m/s
Viscosity: 1 cP = 0.001 Pa.s
Reynolds number (Re)= rvD/µ = 17683.883
Therefore (from Haaland): ff = 0.0100
/8 6 )
c.) 𝐻 = 2𝑓 @ O
NOTE: No change in height, no pumping, no change in pipe diameter.
#DD 6)
𝐻 = 2.0.001 H D.# + 2 ∗ 7 + 4 ∗ 20I E.P#
𝐻 = 2.303 𝑉 $
Dividing by A2
𝐻 = 36157.53 𝑄$
d.) Q = 0.0015; H = 0.08 (see diagram)
e.) Q = 0.0017; H = 0.15 (see diagram)
f.) Q = 0.002; H = 0.13 (see diagram)
0.6
0.5
0.4
System curve
H (m) 0.3
Pump curve
Pumps in series
0.2 Pumps in parallel
0.1
0
0 0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 0.0045
Q (m3/s)
7 Flow Measurement
7.1 Introduction
Measurement of velocity is typically done using:
i) Pitot Tube;
ii) Flow through a constriction (orifice meter, Venturi meter, nozzle);
iii) Flow through a variable area meter (rotameter);
iv) Weirs (or notches); or
v) Other flow measurement methods.
q,>?@
v. F v. dP WF
pgz + r + = pgz + r − ~ +
2 >?@ m 2 =, q,=, ρ m
Typically, the unknown value is the velocity – which is not easy to measure when inside a
metal pipe. However, if we know the pressure and height terms, we can calculate this.
If we assume there is no friction (sometimes not a good assumption) AND no work, the
energy balance becomes:
v. v. P=, P>?@
pgz + r = pgz + r + −
2 >?@ 2 =, ρ ρ
Now:
From a mass balance; mass in = mass out
Min = Mout
ρv=, A=, = ρv>?@ A>?@
1
2. 2 ρ ( P=, − P>?@ )3
V>?@ . = .
A
1 − % A>?@ &
=,
!
( qK* #qQRS )
V>?@ = ! [3]
L
&#Y QRS Z
LK*
1 0 2 3
U-Tube Manometer
where:
1: Point before the orifice plate where the pressure and density (r or v) equals the
original pressure and density, i.e., before any influence of the orifice meter
0: Point of the orifice plate
2: Vena contracta. Point where the flow diameter is at its minimum
3: Pount after the orifice plate where the pressure and density (r or v) equals the
original pressure and density, i.e., after the fluid has returned to its original
flow patterns
!
( q" #q! )
V. = L ! [4]
&#P ! Q
L"
Now: The diameter of the vena contracta is not measured easily. Therefore, define the co-
efficient of contraction:
A.
Ck =
A/
where:
Cc: co-efficient of contraction
A0: Area of the orifice plate
A2: Vena contracta. Area where the flow diameter is at its minimum
Therefore: A2 = CcA0
The flow before, inside and after the orifice plate may be laminar and/or turbulent.
Therefore, we need to re-introduce the alpha (a) term to account for turbulence.
2α.
( )
ρ P& − P.
G = ṁ = ρCE A/ ¡
α CE A/ .
1 − α. % A &
& &
[7]
& &
Reminder: • ≡ γ and : ’ ≡ ρ
CE A/ 2α. . γ( P& − P. )
G = ṁ =
γ ¤ α CE A/ .
1 − α. % A &
& &
In order to account for all the unknown terms and the unknown co-efficient of contraction,
new term is defined: Co-efficient of Discharge (CD)
This is effectively a fudge factor to take in to account the unknowns, turbulent and laminar
flow, frictional losses, and alpha (a) terms.
[8]
CD is obtained from a graph.
Additional Notes:
1) If the orifice meter diameter is much smaller than the diameter of the pipe:
D0 << D1
\A0 << A1
A/ .
1−J K → 1
A&
G = ṁ = C~ A/ i2ρ( P& − P. )
But:
∆P = ρgh
\G = ṁ = C~ A/ ρi2gh
NOTE: The height here is the equivalent height of the fluid in the pipe. This can be
converted to a height of mercury (or other material) as required.
Example:
Water flows through an orifice meter at a rate of 300 cm3/s. Calculate the height of a mercury manometer
(pressure drop) given a pipe diameter of 75 mm and an orifice diameter of 25mm as in the figure below.
dorifice = 25 mm
Water
3
300 cm /s
dpipe = 75 mm
h=?
𝜋
𝐴𝑟𝑒𝑎 𝑜𝑓 𝑜𝑟𝑖𝑓𝑖𝑐𝑒 = . (25)$ = 491 𝑚𝑚$ = 4.91 𝑥 10": 𝑚$
4
𝑐𝑚%
𝐹𝑙𝑜𝑤 𝑜𝑓 𝑤𝑎𝑡𝑒𝑟 = 300 = 3.0 𝑥 10": 𝑚% /𝑠
𝑠
3.0 𝑥 10": 𝑚% /𝑠
∴ 𝐿𝑖𝑛𝑒𝑎𝑟 𝑣𝑒𝑙𝑜𝑐𝑖𝑡𝑦 = = 0.61 𝑚/𝑠
4.91 𝑥 10": 𝑚$
𝑑FI+7+J8 0.025 1
= =
𝑑4+48 0.075 3
2
( )
𝜌 𝑃# − 𝑃$
𝐺 = 𝑚̇ = 𝜌𝐶@ 𝐴D „ $
𝐴
1 − H𝐴D I
#
𝐴D $ 25 $
1−O P = 1 − O P = 0.994 ≈ 1
𝐴# 75
\𝐺 = 𝑚̇ = 𝐶@ 𝐴D 𝜌u2𝑔ℎ
#D./ K0 #DDD0O
𝑚̇ = 3.0 𝑥 𝑥 = 0.30𝑘𝑔/𝑠
= K0
":
0.30 = (0.61 𝑥 4.91 𝑥 10 𝑥 1000)√2.9.81. ℎ
√ℎ = 0.266
ℎ = 0.051𝑚 𝐻$ 𝑂 = 51𝑚𝑚 𝐻$ 𝑂
Converting to mmHg
Pressure of water = Pressure of Mercury
𝜌𝑔ℎR>H8I = 𝜌𝑔ℎSO
1000. 9.81. ℎR>H8I = 13 600.9.81. ℎSO
1000. ℎR>H8I = 13 600. ℎSO
ℎSO = 3.75 𝑚𝑚𝐻𝑔
Kevin Harding (kevin.harding@wits.ac.za) 205
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Flow Measurement – Flow through a Constriction
Example:
Sulphuric acid flows through a pipe of diameter 50 mm. An orifice plate of 10 mm is used determine the flow. If a mercury
manometer gives a reading of 0.1 m, what is the mass flowrate through the pipe?
dorifice = 10 mm
H2SO4
dpipe = 50 mm
h = 0.1m
Now:
𝜌S$YZ: = 1300 𝑘𝑔/𝑚%
0.01 $
𝐴D = 𝜋 O P = 7.85398 𝑥 10"[ 𝑚$
2
𝑃# − 𝑃$ = ∆𝑃 = ∆𝜌𝑔ℎ = (13600 − 1300). 9.81.0.1 = 12066.3 𝑃𝑎
Cd à Guess value of 0.61
Into [1]:
2 2
𝐺 = 𝑚̇ = 𝜌𝐶@ 𝐴D ˆ ( 𝑃# − 𝑃$ ) = 1300. (0.61). 7.8539 𝑥 10"[ . ˆ ( 12066.3) = 0.0622√18.563 = 0.268345𝑘𝑔/𝑠
𝜌 1300
Example:
Sulphuric acid flows through a pipe of diameter 50 mm. An orifice plate of 10 mm is used determine the flow. If
we assume the ratio of (A0/A1)2 is negligible and a mercury manometer gives a reading of 0.1 , the mass flow rate
can be shown to be 0.299139 kg/s. What is the relative error of this reading due to the simplification assumption?
dorifice = 10 mm
H2SO4
dpipe = 50 mm
h = 0.1m
2
( )
𝜌 𝑃# − 𝑃$ "[ ˆ
2 ( 12066.3)
𝐺 = 𝑚̇ = 𝜌𝐶@ 𝐴D „ $ = 1300. (0.68). 7.8539 𝑥 10 . 1300 0.9984 = 0.0694√18.5932
𝐴
1 − H𝐴D I
#
= 0.299378𝑘𝑔/𝑠
A1 (total area)
A2 (annulus area)
Afloat (float area)
NOTE: Taper is
exaggerated
i) Tapered tube with smallest diameter at the bottom with free-flowing float inside
ii) When fluid is flowing, float will rise until force upward balances with gravity
iii) Pressure difference across float is equal to its weight divided by maximum cross-
sectional area in horizontal plane
iv) Area for flow is the annulus between the float and the wall of the tube
v) Consider as orifice meter with variable area
vi) Tube must be tapered otherwise the force upward would be constant (and greater to
gravity) and the float would exit the device
vii) Typically, a rotameter is used instead of an orifice meter since pressure drop is not as
high
From orifice plates the mass flow rate was given by:
2
( )
ρ P& − P.
G = ṁ = ρC~ A/ ¡
A .
1 − %A / &
&
2ρ( P& − P. )
G = ṁ = C~ A/
¤ A .
1 − % A/ &
&
Using a similar derivation as for orifice meters, the relationship as below is obtained for
Rotameters:
Kevin Harding (kevin.harding@wits.ac.za) 208
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Flow Measurement – Flow Through a Variable Area Meter
2ρ( P& − P. )
G = ṁ = C~ A.
¤ A .
1 − % A. &
&
2ρ(−∆P)
G = ṁ = C~ A.
¤ A .
1 − % A. &
&
[1]
Now:
Pressure drop across the float is given by:
where:
Vfloat = volume of float
rfloat = density of the float
rfluid = density of the fluid
Afloat = area of the float
[3]
where:
A1 = area of the pipe at the height of the float
A2 = area of the annulus
CD is dependent on the shape of the float (read off a graph)
Note: If the diameter of the annulus is required, the wetted diameter is to be used.
D2* = D1 – Dfloat
where:
D2* = wetted diameter of the annulus (m)
Dh = 4.Area/Wetted Perimeter
The wetted perimeter is the perimeter that is in contact with the fluid
Method 1:
Because A2 is unknown, CD is also unknown. Therefore, we need to iterate:
Iteration 1:
Guess CD and assume A2 >> A1
Therefore:
2ρcdefg Vcdhij (ρcdhij − ρcdefg )g
G = ṁ = CA A$ → ṁ = CA A$ u2ρcdefg Vcdhij(ρcdhij − ρcdefg )g
’ A $
Acdhij o1 − H $ I s
A#
Solve for A2,1
Iteration 2:
Using A2,1, calculate the velocity through the annulus, the Reynolds number and CD
$k;<=>? l;<@AB (k;<@AB "k;<=>?)m
Returning to the full form of the equation: G = ṁ = CA A$ ˆ C )
n;<@ABo#"WC)X p
'
Assume A2 in the square root is A2,1 and solve for the other A2 (A2,2)
Repeat until BOTH A2,n = A2,n-1 AND CD is constant.
Method 2:
Again CD is unknown. However, this time instead of iterating for BOTH CD and A2, solve for A2 as below and ONLY
iterate for CD as before.
2ρcdefg Vcdhij (ρcdhij − ρcdefg )g
G = ṁ = CA A$
’ A $
Acdhij o1 − HA$ I s
#
Example:
A rotameter has a tube length of 0.3 m and internal diameters of 25 mm (top) and 20 mm
(bottom). The float has a co-efficient of discharge of 0.7, a diameter of 20 mm, density of
4800 kg/m3 and a volume of 6.0 cm3.
What is the mass flow rate through the rotameter if the float is halfway up the tube?
Solution:
Cross sectional area (top of tube): π.(0.025m/2)2 = 4.91 x 10-4 m2
Cross sectional area (bottom of tube): π.(0.020m/2)2 = 3.14 x 10-4 m2
Area of float (Af): 3.14 x 10-4 m2 (same as bottom of tube)
Volume of float (Vf): 6 cm = 6 x 10 m3
3 -6
When the float is half way up the tube (A1): π(0.0225m/2)2 = 3.98 x 10-4m2
Area of annulus (A2): A1 – Af = 0.84 x 10-4m2
Mass flowrate:
2𝜌7?G+& 𝑉7?F>H y𝜌7?F>H − 𝜌7?G+& {𝑔
𝐺 = 𝑚̇ = 𝐶@ 𝐴$
’ 𝐴 $
𝐴7?F>H o1 − H𝐴$ I s
#
Additional Reading
Coulson, JM, Richardson, JF, Backhurst, JR, Harker, JH, 1999. Coulson & Richardson’s
Chemical Engineering, Volume 1, 6th Edition, Elsevier.
Crabtree, M, 2000. Mick Crabtree’s Flow Handbook, 2nd Edition, Crown Publications.
Perry, RH, and Green, DW, 1984. Perry's Chemical Engineers Handbook, 6th Edition.
Welty, JR, Wicks, CE, Wilson, RE, Rorrer, GL. 2008. Fundamentals of Momentum, Heat,
and Mass Transfer, 5th Edition, Wiley.
Problems
7.1 Explain the principle of the rotameter.
a.) Draw a rotameter showing the important features.
b.) Why is the top of the rotameter wider than the bottom?
c.) Is a rotameter placed vertically or horizontally? Could it be placed at 45 degrees?
d.) What is the float? What different floats are available?
e.) What is CD?
f.) Find a graph relating float type/shape to CD.
g.) What is the hydraulic mean diameter?
7.2 Oil flows through a 5” steel pipe at a rate of 300 US gal per minute. In the pipeline is a 3.5” standard sharp-edged orifice plate, to which is attached a mercury manometer. At the
flow temperature, the oil has a specific gravity of 0.87 and a viscosity of 15 cp. What will be the reading on the manometer?
I.D. of pipe = 5.047” US gal = 0.1338 ft3
7.3 Brine of s.g. 1.2 flows through a standard 3” pipe at 185 US gal/min. In order to measure the flow, it is desired to install an orifice meter in the line. The pressure difference is to be
measured on a simple mercury manometer. If the maximum reading is to be 40 cm, what size orifice should be installed? Estimate the accuracy to which it will be possible to
measure the flow by this method.
7.4 Oil of viscosity 11 cp and s.g. 0.9 flows through a 1” orifice installed in a 3” pipe. The pressure difference across the orifice is measured as 1.6” of oil. What is the flow rate in the
pipe?
7.5 A rotameter has a tube 12” long, which has an internal diameter of 1” at the top and 0.75” at the bottom. The float (of type B) has a diameter of 0.75”, a volume 0.4 in3 and a specific
gravity of 4.80. What will be the height of the float when there is a flow of 90 gallons of water per hour through the rotameter?
7.6 Water flows through a 30 cm rotameter (Type A) with an internal diameter of 25 mm at the top and 20 mm at the bottom. The float has a diameter of 20mm, a volume of 6.5 cm3
and a specific gravity of 5.
f.) Assuming that A2/A1 is negligible, what is the area of the annulus when there is a flow of 410 l/hr water?
g.) Recalculate the area of the annulus assuming A2/A1 is not negligible. Assume a constant CD.
h.) Compare and discuss the results of (a) and (b).
i.) Is the assumption of a constant CD valid?
7.7 A rotameter has a scale ranging from 0.014 m3/min to 0.14 m3/min. It is intended to use this meter for metering a gas of density 1.8 kg/m3 within a flow range of 0.0287 m3/min to
0.287 m3/min. A new float of the same dimensions is to replace the old one (density = 1905 kg/m3) in order for this system to work.
Given:
2ρcdefg Vcdhij (ρcdhij − ρcdefg )g
G = ṁ = CA A$
’ A $
Acdhij o1 − HA$ I s
#
a.) What density must the new float have? Both the floats can be assumed to have the same volume and shape. Further, it can be assumed that the coefficient of discharge is
constant.
b.) Explain (and include a diagram for) the different area terms in the equation above.
c.) How does a rotameter differ from an orifice place? How is Cd different in each?
Solution:
Simplifying:
D.[
𝐺 = 𝑚̇ = 𝑘u𝜌7?F>H − 𝜌7?G+& = 𝑘y𝜌7?F>H − 𝜌7?G+& {
Where k = all values EXCEPT the square root of density float minus fluid
Using the ratios of the new mass flowrate to the old flowrate and the equation above:
,.J ,.J
K̇"+H D.D$P W 0q5DE#G% "5DE$!Fr X W q5DE#G% "5DE$!F r X
K̇#EF
=
D.D#:
=2= ,.J
"+H
= ,.J
"+H
because K is constant
W 0q5DE#G% "5DE$!Fr X W q5DE#G% "5DE$!Fr X
#EF #EF
,.J
W q5DE#G% "#.Pr X
And: 2 = ( (#ED["#.P),.J)
"+H
#EF
Coefficient of Discharge (CD): Based on diameter of pipe vs. diameter of orifice plate, based on Re, more variable over Re than rotameter, at high Re CD is constant (single)
value
7.8 Water flows through a 30 cm rotameter (Type A) with an internal diameter of 25 mm at the top and 20 mm at the bottom. The float has a diameter of 20 mm, a volume of 6.5 cm3
and a specific gravity of 4.8.
a.) Assuming that A2/A1 is negligible, what is the area of the annulus when there is a flow of 390 l/hr water? Assume CD = 0.97.
b.) Recalculate the area of the annulus assuming A2/A1 is not negligible. Assume a constant CD.
c.) Compare and discuss the results of (a) and (b).
Solution:
Solving: A2 = 8.99 × 10-5 m2 = 89.9 mm2 (i.e. VERY small and the float is hardly raised)
b.) Iterating, using the full form of the equation:
Iteration 1: 79.29 mm2
Iteration 2: 80.35 mm2
Iteration 3: 80.24 mm2
Iteration 4: 80.25 mm2
Iteration 5: 80.25 mm2
c.) Numbers are different. We had to reintroduce the A2/A1 term.