Chemical Engineering With Computer Applications
Chemical Engineering With Computer Applications
Hussein K. Abdel-Aal
MATLAB® is a trademark of The MathWorks, Inc. and is used with permission. The MathWorks does not
warrant the accuracy of the text or exercises in this book. This book’s use or discussion of MATLAB® soft-
ware or related products does not constitute endorsement or sponsorship by The MathWorks of a particular
pedagogical approach or particular use of the MATLAB® software.
CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742
This book contains information obtained from authentic and highly regarded sources. Reasonable efforts
have been made to publish reliable data and information, but the author and publisher cannot assume
responsibility for the validity of all materials or the consequences of their use. The authors and publishers
have attempted to trace the copyright holders of all material reproduced in this publication and apologize to
copyright holders if permission to publish in this form has not been obtained. If any copyright material has
not been acknowledged please write and let us know so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced, transmit-
ted, or utilized in any form by any electronic, mechanical, or other means, now known or hereafter invented,
including photocopying, microfilming, and recording, or in any information storage or retrieval system,
without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copyright.
com (http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood
Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that provides licenses and
registration for a variety of users. For organizations that have been granted a photocopy license by the CCC,
a separate system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are used
only for identification and explanation without intent to infringe.
v
vi Contents
2.9.2 Thermodynamics........................................................26
2.9.3 Some General Tips .....................................................26
2.9.3.1 Parts-Per Notation......................................26
2.9.3.2 Frequently Asked Questions ....................... 27
End-of-Chapter Solved Examples ...................................................... 27
Reference ............................................................................................ 32
Section ii Fundamentals and Problem-Solving
Profile
xiii
xiv Preface
The material in each chapter in Section III, excluding Chapter 10, encompasses two
sections: one covers the theoretical principles in a condensed format, and the other
presents a set of numerically solved problems using interactive numerical software
packages.
Chapter 10, on the other hand, introduces a number of case studies that handle
practical technical problems. Many of these case studies are compiled by the author
in the field of applied research, in both academy and industry.
This primer in chemical engineering would serve as a text for students majoring
in chemical engineering, applied chemistry, and biochemical process engineering.
This book is designed to be used by students in a variety of ways, dependent on
the level of enrollment. Possible recommendations are as follows:
This book takes a highly pragmatic approach to presenting the principles and appli-
cations of chemical engineering and offering a valuable, easily accessible guide to
solve problems using computers. MATLAB and spreadsheeting are applied in solv-
ing many examples.
Problems assigned for each chapter and found in Appendix A, ranging from simple
to difficult, allow readers to gradually build their skills and be able to tackle a broad
range of problems in chemical engineering. In addition, the numerous real examples
throughout this book (more than 70) include computer or hand solutions, or in many
cases both. Many solved example problems reinforce the concepts covered. The text is
written with a view that solving problems is an essential part of the learning process.
It gives a taste of the kinds of problems with which chemical engineers grapple. It puts
a wealth of solved practical problems at your fingertips.
One of the main unique features of this book is that it is a convenient source in
chemical engineering that can be used at ease, without the need for many compre-
hensive texts. It is an integrated form of a single source.
In conclusion, one can say that two things are important: “Science…means un-
resting endeavor and continually progressing development toward an aim which the
poetic intuition may apprehend, but the intellect can never fully grasp.” Chemical
engineering is the discipline that uses a molecular understanding of matter to
produce—from raw materials—all the synthetic materials humankind is using.
xvii
Author
Professor Hussein K. Abdel-Aal is emeritus professor of chemical engineering and
petroleum refining, NRC, Cairo, Egypt. He received his BS in chemical engineering
in 1956 from Alexandria University and MS and PhD in chemical engineering in
1962 and 1965, respectively, from Texas A&M University, College Station, Texas.
Professor Abdel-Aal joined the Department of Chemical Engineering at
KFUPM, Dhahran, Saudi Arabia (1971–1985). He was the head of department for
the period 1972–1974. He was also visiting professor with the Chemical Engineering
Department at Texas A&M (1980–1981). In 1985–1988, he assumed the responsi-
bility of the head of the Solar Energy Department, NRC, Cairo, before rejoining
KFUPM for the period 1988–1998.
Professor Abdel-Aal conducted and coordinated projects involving a wide range
of process development, feasibility studies, industrial research problems, and con-
tinuing education programs for many organizations, including Suez Oil Processing
Company, Petromin in Riyadh (Saudi Arabia), Arab Petroleum Investment Corp in
Dhahran (Saudi Arabia), Hagler-Baily & Company, Washington DC, Mobil/Esso Oil
Companies in Libya, and Kuwait oil companies.
Professor Abdel-Aal has contributed to more than 90 technical papers, is the edi-
tor of Petroleum Economics & Engineering, third edition (2014), and is the main
author of the textbook Petroleum and Gas Field Processing, second edition, (2016).
Both books were published by Taylor & Francis Group/CRC Press.
Professor Abdel-Aal is listed in Who’s Who in the World (1982) and is a member
of AICHE, Sigma Si, Phi Lambda Upsilon. He is a fellow and founding member of
the board of directors of the International Association of Hydrogen Energy.
xix
Introduction: Initial Thoughts
I.1 OBJECTIVES
For a student, taking the first course in their major is an exciting but scary step
into the unknown. Effectively introducing chemical engineering education to under-
graduate students is usually difficult, especially when teaching introductory courses.
A prime pitfall of the larger texts is that students are disinclined to read several pages
of text to solve a problem.
Experience in teaching the first course has shown that there is a need for such a
handy, collective, and summarized source of basic principles in chemical engineer-
ing, especially when solving problems. The methodology of solving problems and
their solutions should be presented hand in hand with the fundamentals, as discussed
in this text.
This primer is to be considered a survival companion for students of chemical
engineering. As you navigate through the text, you can maneuver smoothly from
chapter to chapter. The two main themes embodied in this text are as follows:
• Theme No. 1: To provide students with some guide steps on how to identify
the type of problems to solve and how to target a solution. Focus will be on
the use of systematic algorithms that employ numerical methods to solve
different chemical engineering problems by describing and transforming
the information.
• Theme No. 2: To present the basic principles and techniques in chemical
engineering and the underlying unit operations and chemical processes in a
concise and nontraditional format with numerous applications.
Let us first acquaint our readers with the field of chemical engineering. All engi-
neers employ mathematics, physics, and the engineering art to overcome technical
problems safely and economically. Yet it is the chemical engineer alone that draws
upon the vast and powerful science of chemistry to solve a wide range of problems.
Chemical engineers use chemistry and engineering to turn raw materials into usable
products, such as petrochemicals, plastics, and pharmaceutical ingredients. They
apply physical science together with mathematics and computers to processes that
involve physical changes and/or chemical reactions in order to produce more valu-
able products.
Every scientific discipline, including chemical engineering, relies on the use of
computers and software. The capabilities provided by computers, such as fast cal-
culations, large storage capacity, and software, permit engineers to solve different
problems.
xxi
xxii Introduction: Initial Thoughts
The theoretical component involves the use of systematic algorithms to solve dif-
ferent scientific problems by describing and transforming the information, while the
practical component involves the implementation of the computational hardware
and software.
Traditionally, in solving scientific problems, researchers use either of the following:
• Theoretical approach
• Laboratory experiment
In the theoretical approach, scientists rely upon the use of some classified models,
in accordance with how they are derived. In general, three models are well known:
Computational science, on the other hand, is different from these two traditional
approaches, the theoretical and the practical, as illustrated in Figure I.1. The scientific
Architecture
(computing environment)
n
tio
Ex
uta
pe
ri
mp
me
Co
Model
nt
Theory
Application Algorithm
(science) (mathematical model)
The physical operations (known as unit operations and shown in the earlier exam-
ple) are fluid flow, heat transfer, and distillation. They are basically based on these
xxiv Introduction: Initial Thoughts
Crude
2 atm 43°C Reflux
oil Reflux
drum
modes of transfer. Unit operations deal chiefly with the transfer of energy and the
transfer, separation, and conditioning of materials by physical means.
At this stage, two basic questions arise:
The answer to the first question deals with theory of transfer or transport, as explained,
within the boundaries of our system. In answer to the second question, the combined
effect of momentum, heat, and mass (MHM) is responsible for the physical changes
that take place in the distillation column to produce the finished products.
Further, assume that the gasoline exiting the distillation column is introduced into
what is known as a reforming unit in order to obtain a higher-grade gasoline. This
reforming process represents a typical example of a chemical conversion or chemi-
cal reaction—known as unit process—where hydrocarbons undergo molecular
changes and rearrangement leading to high-octane gasoline. Unit processes involve
primarily the conversion of materials by means of chemical reactions. Again, it
should be pointed out that the three modes of transfer, MHM, take place for opera-
tions involving chemical reactions, or chemical changes, as well.
Introduction: Initial Thoughts xxv
Now, having presented this example, one should elaborate the concepts of unit
operations, unit processes, and transport phenomena. This will pave the way to the
students when they come to the heart of the materials embodied in the text.
All processing schemes can be composed from and decomposed into a series of indi-
vidual or unit steps.
xxvi Introduction: Initial Thoughts
These unit operations cut across widely different processing applications, including
the manufacture of chemicals, fuels, pharmaceuticals, pulp and paper, and processed
foods. The unit operation approach serves as a very powerful form of morphological
analysis, which systematizes process design and greatly reduces both the number of
concepts that must be taught and the number of possibilities that should be consid-
ered in synthesizing a particular process, especially in plant design.
Most unit operations are based mechanistically upon the fundamental transport
processes of mass transfer, heat transfer, and fluid flow (momentum transfer), as
shown in the following table:
It is important to point out that the study of transport phenomena provides a unifying
and powerful basis to understand different unit operations. The subject of transport
phenomena includes three closely related topics as shown earlier; of interest is mass
transfer, which is concerned with the transport of mass of various chemical spe-
cies. They frequently occur simultaneously in industrial applications like distilla-
tion, where the transport of mass occurs by diffusion for various chemical species.
The basic equations that describe the three modes of transport phenomena are also
closely related and can often be analyzed by analogy:
¶rux
tzx = - n
¶z
Introduction: Initial Thoughts xxvii
where
τzx is the flux of x-directed momentum in the z direction
ν is μ/ρ, the momentum diffusivity
z is the distance of transport or diffusion
ρ is the density
μ is the viscosity
Newton’s law is the simplest relationship between the flux of momentum and the
velocity gradient.
The net flux of energy through a system, on the other hand, equals the conductiv-
ity, k times the rate of change of temperature with respect to position dT/dx.
where A is the surface area and q is the heat flow per unit time.
Finally, for mass transfer, the diffusion flux, J, from a higher concentration to a
lower concentration is proportional to the gradient of the concentration, dc/dx, of the
substance and the diffusivity, D, of the substance in the medium.
Based on this presentation, the following conclusions are drawn:
However, it is believed that there is no real conflict between the two approaches: the
unit operations and the transport phenomena. Most important is the development of
the skills that lead to the creative engineering use of the fundamentals or a combina-
tion of both, as demonstrated and presented in our book.
REFERENCE
American Institute of Chemical Engineers. AIChE Constitution. http://www.aiche.org/About/
WhoWeAre/Governance/Constitution.aspx, 2003. Retrieved August 13, 2011.
Section I
An Insight into Chemical
Engineering
Much of the background needed for an understanding of the materials embodied in
the text is basic science and engineering. To analyze and follow-up quantitatively
and to optimize the underlying processes in chemical engineering a profound knowl-
edge of theoretical background is needed. Section I, consisting of Chapters 1 and 2,
highlights a technical introduction to the basic principles and calculation techniques
introduced in the text.
Solved examples given in the chapters are straight forward direct applications.
They are hand-solved. It is meant to be a kind of refresher for a first course in chemi-
cal engineering curriculum. Some are worked out by Excel. Application problems
for Section I are included in Appendix A.
1 Introductory Concepts
As indicated by the title of this chapter, the essentials of some chemical engineering
aspects are introduced. Most of the common systems of units used in calculations are
presented first. Temperature (T) and pressure (P) scales are discussed in detail next,
since both T and P are key process variables. For example, to design a distillation
column, information is needed about the boiling temperatures of the components
that should be separated. Physical gas laws that are frequently used in the calcula-
tions are described as well.
• Length L
• Time θ
• Mass M
• Temperature T
To express any physical property, you need two items: a value and a unit. For exam-
ple, we say the room temperature is 25°C. Another example is the dimensions of a
distillation column: 30 ft in length and 4.5 ft in diameter.
The most common systems of units are as given next:
In addition, we have what we call derived units, for example, velocity = ft/s = L/θ.
The International System of Units (abbreviated SI from the French Le Système inter-
national d’unités) is the modern form of the metric system. The standards, published
in 1960, are based on the meter–kilogram–second system rather than the centimeter–
gram–second (cgs) system, which, in turn, had several variants. It comprises a coherent
system of units of measurement built around seven base units, as shown in Figure 1.1.
The system has been nearly globally adopted. In the United States, metric units
are not commonly used outside of science, medicine, and the government. However,
* This section is modified after International Bureau of Weights and Measures (n.d.; 2016).
3
4 Chemical Engineering Primer with Computer Applications
A s
mole m
cd kg
FIGURE 1.1 The seven base units of International System of Units: Kelvin (temperature),
second (time), meter (length), kilogram (mass), candela (luminous intensity), mole (amount
of substance), and ampere (electric current). (With kind permission from Springer: Units
of Measurement: Past, Present and Future, International System of Units, 2009, 16,
Gupta, S.V.)
U.S. customary units are officially defined in terms of SI units. The United Kingdom
has officially adopted a partial metrication policy, with no intention of replacing
imperial units entirely.
The basic laws that we use can be easily memorized using the following math-
ematical relationships:
Newton
N
Joule
J
Watt
W
Introductory Concepts 5
It is worth mentioning that many of the derived units in the SI system are given
special names with corresponding symbols. For example, force, which has the units
kg · m/s2, has been named Newton and given the symbol N. Similarly, energy is
defined as Newton ∗ meter or N · m, having the units m2 · kg/s2. Hence, it is named
Joule and given the symbol J, and the corresponding unit of power is joule per sec-
ond named Watt.
Now, let us examine the systems of units when it comes to using the conversion
factor, known as gc. Basically, we start with the Newton’s law:
F = C×m ×a (1.1)
where
F is the force
m is the mass
a is the acceleration
C is a constant, whose numerical value and units depend on those selected for the
units F, m, and a
• The cgs system: The unit of force is dyne where 1 g is accelerated at 1 cm/s2.
Equation 1.1 is rewritten as
1. Temperature
2. Pressure
3. Flow
4. Level
We will limit our talk to the first two process variable, T and P.
1. Fahrenheit
2. Celsius
3. Kelvin
Fahrenheit temperature scale is a scale based on 32 for the freezing point of water
and 212 for the boiling point of water, the interval between the two being divided
into 180 parts.
Celsius temperature scale, also called centigrade temperature scale, is the scale
based on 0 for the freezing point of water and 100 for the boiling point of water. The
following formula can be used to convert a temperature from its representation on
the Fahrenheit (F) scale to the Celsius (C) value: C = 5/9(F − 32). The scale is illus-
trated in Figure 1.2.
Kelvin temperature scale is the base unit of thermodynamic temperature mea-
surement in the SI of measurement. Such a scale has as its zero point (absolute zero),
the theoretical temperature at which the molecules of a substance have the lowest
energy, as shown in Figure 1.3. Many physical laws and formulas can be expressed
more simply when an absolute temperature scale is used; accordingly, the Kelvin
scale has been adopted as the international standard for scientific temperature mea-
surement. The Kelvin scale is related to the Celsius scale. The difference between
the freezing and boiling points of water is 100° in each so that the Kelvin has the
same magnitude as the degree Celsius.
The following equation is used for temperature intervals rather than specific tem-
peratures: 1 K = 1°C = 9/5°F = 9/5°R.
Introductory Concepts 7
Water
212°F 100°C
boils
50°C
98.6°F 37°C
Water
32°F 0°C
freezes
–40°F –40°C
–50°C
Dry ice
–108°F –78°C
(solid CO2)
–100°C
–150°C
Liquid –191°C
–312°F air –200°C
–250°C
Absolute –273°C
–459°F
zero
Fahrenheit Celsius
(centigrade)
• 1.0 atm
• 14.7 psia
• 760 mm Hg
• 29.92 in. Hg
• 33.91 ft water
• 100 × 103 Pa (pascals) or N (Newton) per sq. meter (N/m2) = 100 kPa
8 Chemical Engineering Primer with Computer Applications
400 K
Water
100°C 373 K
boils
350 K
50°C
300 K
0°C Water 273 K
freezes
250 K
–40°C
–50°C
Dry ice 200 K
–78°C
(solid CO2) 195 K
–100°C
150 K
–150°C
–191°C 100 K
–200°C Liquid 82 K
air
50 K
–250°C
Absolute
–273°C 0K
zero
Celsius Kelvin
(centigrade) (absolute)
Gauge pressure
Local
atmospheric Absolute pressure
pressure
Absolute zero
(complete vacuum)
If pressure is measured by means of height of the column of liquid, other than mer-
cury, the following equation is used to relate height to pressure:
Pa = Pb + rgh
where
Pa is the pressure at the bottom of the column
Pb is the pressure at the top of the column
ρgh are the density of fluid, acceleration of gravity, and height of fluid column,
respectively
V1 T1
= (1.7)
V2 T2
PT = Sp i = p1 + p2 + p3 + + p n (1.8)
pA pA
XA = = (1.9)
PT Sp i
10 Chemical Engineering Primer with Computer Applications
The most useful form in which the laws of a mixture of ideal gases can be put is as
follows:
PV = nRT (1.10)
where
R is the gas constant
n is the number of moles
T is the temperature
P is the pressure; both are in absolute units
PV
R=
nT
where
P is the pressure
V is the volume
n is the chemical amount
T is the temperature
As pressure is defined as force per unit area, the gas equation can also be written as
Force/Area × Volume
R=
Amount × Temperature
Force/(Length)2 ´ (Length)3
R=
Amount ´ Temperature
Introductory Concepts 11
Since
Work
R=
Amount ´ Temperature
The physical significance of R is work per degree per mole. It may be expressed in
any set of units representing work or energy (such as joules), other units representing
degrees of temperature (such as degrees Celsius or Fahrenheit), and any system of
units designating a mole or a similar pure number that allows an equation of mac-
roscopic mass and fundamental particle numbers in a system, such as an ideal gas.
Values of R are found in Appendix B.
How much gas is present could be specified by giving the mass instead of the
chemical amount of gas. Therefore, an alternative form of the ideal gas law may
be useful. The chemical amount (n) (in moles) is equal to the mass (m) (in grams)
divided by the molar mass (M) (in grams per mole):
m
n=
M
By replacing n with m/M and subsequently introducing the density ρ = m/V, we get
m
PV = RT
M
R
P=r T
M
The ideal gas law is the equation of state of a hypothetical ideal gas. It is a good
approximation to the behavior of many gases under many conditions, although it has
several limitations, as shown in Figure 1.5.
Nevertheless, this form of the ideal gas law is very useful because it links pres-
sure, density, and temperature in a unique formula independent of the quantity of
the considered gas. The behavior of real gases usually agrees with the predictions
of the ideal gas equation to within ±5% at normal temperatures and pressures. At
low temperatures or high pressures, real gases deviate significantly from ideal gas
behavior.
One of the most useful equations to predict the behavior of real gases was devel-
oped by Johannes van der Waals (1837–1923). The complete van der Waals equation
is therefore written as follows:
ìï é n 2a ù üï
íP + ê 2 ú ý{V - nb} = nRT (1.11)
îï ë V û þï
12 Chemical Engineering Primer with Computer Applications
7.9
0 5 10 15 20
P (atm)
From the definition of partial pressure as given before, one can say for a number of
components of a gas mixture:
Since
p1 n
V1 = V2 and T1 = T2 ; = 1
p2 n 2
pi n
= i , or p i = PT y i (1.12)
PT n T
where yi is the mole fraction of the ith component in the total mixture of n compo-
nents. The next relationship provides a way to determine the volume-based concen-
tration of any individual gaseous component
Ci
p i = PT (1.13)
106
component present in the solution. Therefore, the individual vapor pressure for each
component is given by the next equation:
p i = Pi*x i (1.14)
where
pi is the partial pressure of the component i in the mixture (in the solution)
Pi* is the vapor pressure of the pure component i
xi is the mole fraction of the component i in the mixture (in the solution)
P*x i = PT y i
or
P* Yi
= = K i (applied only for ideal mixtures)
PT x i
Therefore,
y i = K ix i (1.15)
x i = Hip (1.16)
where
xi is the solubility of a gas at a fixed temperature in a particular solvent (in units
of M or mL gas/L)
Hi is Henry’s law constant (often in units of M/atm)
p is the partial pressure of the gas (often in units of atm)
Henry’s law states that at a constant temperature, the amount of a given gas that
dissolves in a given type and volume of liquid is directly proportional to the partial
pressure of that gas in equilibrium with that liquid.
Example 1.1
[h] = QL-2q-1T -1
= (369)(0.252)/(0.3408)2 (1/18
. ) kcal/(m)2 (°C)(h)
Example 1.2
1 2 3 4
Comp % MWt (2) (3)
O2 0.21 32 6.72
N2 0.78 28 21.84
Ar 0.01 20 0.2
Avg. MW 28.76
Introductory Concepts 15
Example 1.3
Calculate the total pressure and the composition of the following gas mixture,
given the partial pressure of each component:
Solution
Substitution using Equation 1.8 gives the following:
PT = å pi = 75 + 50 + 595 + 26 = 746 mm Hg
75
XCO2 = = 0.101
746
Similarly, calculations are carried out to find the mole fractions of CO, N2, and O2:
0.067, 0.797, and 0.035, respectively.
Example 1.4
Five hundred pounds of water is flowing a heat exchanger pipe at the rate of
20.0 ft/s. Calculate the kinetic energy of the flowing water in (ft)(lbf ).
Solution
1 2 1 é 1 ù
K.E. = mv = éë(500 lbm )(20 ft/s)2 ùû ê 2 ú
2 2 ë 32.174 (ft)(lbm )/(s )(lbf ) û
Example 1.5
Calculate the volume of 1 mol of air at 20°C on top of a mountain (altitude 4.2 km),
where the air pressure is approximately 60 kPa.
Solution
Apply Boyle’s and Charles’ laws as successive correction factors to the standard
sea-level pressure of 101.3 kPa:
The standard molar volume 22.4 L/mole is applicable at standard temperature and
pressure STP.
16 Chemical Engineering Primer with Computer Applications
Example 1.6: To Check the Validity of the van der Waals equation
Compare the experimental PV data for propane gas to the predicted values as
shown next:
6.0 P experimental
5.0
4.0
3.0
2.0
1.0
0.0
0 10 20 30 40 50 60
Conclusion: The PV values for both experimental and predicted ones match
closely.
Example 1.7
Solution
pi = PTyi.
PT is the sum of pi = pCO2 + pCO.
First: At 30 atm,
p2CO
Kp =
pCO2
(302 )(0.83)2
= = 1215
. atm.
(30)(0.17)
Example 1.8
Compute the work in joules if a steady force of 3 lbf is needed to push 5 lbm in a
distance of 100 ft.
Solution
= 406.7 J
Example 1.9
Calculate the molar volume (ft3/lb · mole) and the density (lb/ft3) of methane gas
at 60°F (519.67°R) and a pressure of 3.5 atm. Molecular weight (M) for methane
is 16.04.
Solution
Methane density is
= 0.148 lb /ft 3
18 Chemical Engineering Primer with Computer Applications
Example 1.10
Example 1.11
Find the absolute pressure (mm Hg) of a gas under draft of 5.0 in water. The baro-
metric pressure reads 745.0 mm Hg.
Solution
= 9.338 mm Hg
Therefore, the absolute pressure in the duct is 745.0 + (−9.338) = 735.66 mm Hg.
REFERENCES
BIPM, n.d. SI Brochure, Section 2.1.1.5, http://www1.bipm.org/en/si/si_brochure/chapter2/
2-1/2-1-1/kelvin.html.
International Bureau of Weights and Measures, n.d., Brief History of the SI, http://www.bipm.
org/en/si/history-si/, retrieved April 21, 2009.
International Bureau of Weights and Measures,” Wikipedia, last modified date Aug. 2, 2016,
https://en.wikipedia.org/wiki/International_Bureau_of_Weights_and_Measures.
2 Basic Principles
and Introduction
to Calculations
In this chapter, a review of some of the well-known physical and chemical laws is
given in continuation to Chapter 1. In addition, the principles underlying chemical
engineering problems are presented, showing readers how to apply such principles
to the different topics covered in the text.
Definitions of many of the concepts and terms that underlie much of the chemical
and process industries are summarized, such as unit operations and unit processes,
to name a few. Design variables for sizing piping and process equipment are briefly
presented. Systems of units, different temperature scales, and pressure scales are
described.
Much of the materials given in this chapter would add basic knowledge to the
students and enhance the ability in problem solving.
Section I, consisting of Chapters 1 and 2, would pave the way for the readers,
especially when it comes to topics on unit operations and plant design. It may be
considered as a supplement to the chapters found in Section III.
1. In many chemical processes, each one can be broken down into a series of
steps or operations that are repeated in other processes.
2. These individual steps or operations are based on the same scientific prin-
ciples and guidelines.
19
20 Chemical Engineering Primer with Computer Applications
Unit operations deal mainly with the transfer and the change of both materials
and energy primarily by physical means, arranged as needed by a chemical process.
The following is a partial list of some important unit operations:
• Fluid flow: It deals with the governing principles for flow and transporta-
tion of fluids.
• Heat transfer: It deals with the underlying principles of heat transfer by
different modes.
• Diffusional mass transfer unit operations include distillation, absorption,
extraction, and drying. Separation in these unit operations is accomplished
by the transfer of molecules from one phase to the other by diffusion.
Chemical
engineering
1. The design and the operation of plants in chemical industry and related
processes. This group is known as chemical process engineers.
2. The development of new chemicals or products ranging from foods and
beverages to cosmetics to cleaners to pharmaceutical ingredients, among
many other products. This group is known as chemical product engineers.
Purely mechanical changes are usually not considered part of the chemical processes.
A chemical plant is schematically shown in Figure 2.1.
As a result of the transformation of raw materials in the presence of energy-
producing finished products, chemical processes are analyzed and studied using the
following chemical engineering tools:
R. M.
Chemical plant Product
Energy
FIGURE 2.1 Chemical plant with input of raw materials and energy source.
22 Chemical Engineering Primer with Computer Applications
4. Unit operations
5. Instrumentation and control
6. Economics
These six subjects, known as the backbone courses in a typical BS degree curriculum,
are closely dependent on each other, but their principles can be studied individually.
–4.5
–5.0
Straight line
ln[NO2]
–5.5
Data
–6.0
In other words, given the value of ΔG, calculate k (the chemical reaction
equilibrium constant), hence the composition of phases at equilibrium).
3. Kinetics: This law is used to determine the rate at which a chemical com-
pound reacts, that is, finding the rate of reaction, k. It is needed in the design
of chemical reactors.
4. Kinetics versus thermodynamics.
As stated earlier, the maximum extent to which a chemical reaction can proceed may
be calculated by thermodynamics, which occurs at equilibrium. It follows that the
net rate of a chemical reaction must be zero at this equilibrium point. Thus, a plot of
reaction rate (r) versus time would always approach zero, as the time approaches ∞,
as illustrated in Figure 2.2 for the thermal decomposition of NO2 to NO + O2.
Concentration of products
k= .
Concentration of reactants
å k x -1 = f (T
i i BP ),
V
y
F
xF
L
x
This equation represents the mathematical function for the physical statement of
the BP. Therefore, the solution for the BP becomes a trial and error problem by solv-
ing the function, since K i is a function of T or K i = f(Ti).
When the two laws are merged together, we get what is known as the law of con-
servation of mass–energy: the total amount of mass and energy in the universe is
constant, where mass and energy can interconvert.
For heat balance calculations for a system, the energy or heat input in the form of
enthalpy processed by materials in, plus any heat added, should be equal to heat out.
26 Chemical Engineering Primer with Computer Applications
2.9.2 thermodynamIcs
The four laws of thermodynamics define the fundamental physical quantities (tem-
perature, energy, and entropy) that characterize thermodynamic systems. The laws
describe how these quantities behave under various circumstances and forbid certain
phenomena (such as perpetual motion):
Example 2.1
Sodium sulfite is used to remove dissolved oxygen from boiler-feed water to avoid
corrosion, as given by the equation:
2 Na2SO3 + O2 → 2Na2SO4
Example 2.2
p2CO
Kp =
pCO2
(302 )(0.83)2
= = 1215
. atm.
(30)(0.17)
2
Second: At 20 atm, KP = 1215 . = pCO /pCO2 = {20 - pCO2 }/pCO2.
Solving for pCO2, we obtain pCO2 = 2.55 atm.
Percent CO2 = [2.55/20] × 100 = 12.75.
Example 2.3
Calculate the volume of 1 mol of air at 20°C on top of a mountain (altitude 4.2 km),
where the air pressure is approximately 60 kPa.
Solution
Apply Boyle’s and Charles’ laws as successive correction factors to the standard
sea-level pressure of 101.3 kPa:
Basic Principles and Introduction to Calculations 29
293 K 1013
. kPa
V2 = (22.4 L/mole)
273 K 60 kPa
= 41L
The standard molar volume 22.4 L/mole is a value that it is valid only at STP.
Example 2.4
For the industrial production of CO2, fuel oil is injected into a furnace to be burned
with exactly 1.5 times the theoretical amount required for complete combustion.
Calculate the exact composition of the stack gases, and determine the CO2 purity
of the product to be sold. For simplicity, assume the fuel oil to be C12H26.
Solution
Products kg · mole %
CO2 (100)(12) = 1,200 1,200/13,732 = 8.7
H2O (100)(13) = 1,300 9.5
O2 925 6.7
N2 10,307 75.1
Total 13,732 100
Example 2.5
W kg/h water
Solution
This is a physical separation process. We have two unknowns to calculate by using
total MB (TMB) and component MB (CMB).
TMB:
1000 = W + B (2.1)
CMB on solids:
Solving the two equations in two unknowns, we obtain W = 878 kg/h and
B = 122 kg/h.
Example 2.6
Example 2.7
Calculate the kinetic energy of a missile moving at 12,000 mile/h above the earth,
with an acceleration due to gravity of 30 ft/s2.
Solution
We arbitrarily take 1 lb mass as a basis.
1 2
K.E. = mv
2
1
= (1 lbm )[1 lbf × (s2 ) / 32.2 ft × lbm ][(12, 000 miles)2 /h2 ]
2
´ [(5280 ft )2 /mile2 ][h2 /(3600 s)2
= 4.81 ´106 (ft)(lbf )
Basic Principles and Introduction to Calculations 31
Example 2.8
The well-known equation to calculate the pressure drop due to friction losses for
fluids in pipes is given by ΔP = [2fLρv2/D], where ΔP is the pressure drop, v is the
velocity, L is the length of pipe, ρ is the density of fluid, and D is the diameter of pipe.
Find out if this equation is dimensionally consistent. Support your answer.
Solution
DP = [2fLrv2 /D] = (ft /ft )(lbm /ft 3 )(ft /s)2[lbf × (s2 ) /lbm × ft ]
= lbf /ft 2
Example 2.9
Assuming that 1 kg of CO2 is released per day per person, calculate the amount
of LiOH required.
Solution
Basis 1 kg CO2
[1kg CO2 ][1kg × mole CO2 /44.01kg CO2 ][2 kg × mole LiOH/1kg × mole CO2 ]
Example 2.10
A solvent recovery system delivers an inert gas, saturated with benzene vapor
(C6H6). This gas is at 70°F and 750 mm pressure. The gas is compressed to 5 atm
and cooled to 70°F, after compression.
Calculate how many pounds of benzene are condensed per 1000 ft3 of the
feed gas. The vapor pressure of benzene at 70°F is 75 mm.
Feed gas
Benzene
32 Chemical Engineering Primer with Computer Applications
Solution
The basis is 1000 ft3 of feed gases consisting of benzene plus inert. Applying
volume % = pressure %, where the total pressure is 750 mm, while the partial
pressure of benzene is 75 mm, the volume % of the inert = (750 − 75)/750.
The volume of inert gas in the feed = (1000)(750 − 75)/750 = 900 ft3.
Then the volume of benzene = 100 ft3.
Convert the volumes to moles: moles of inert = 2.3 and moles of benzene =
0.2556.
The ratio of benzene to inert, before compression, is 0.2556/2.3 = 0.1111.
Now after compression and cooling, a portion of benzene will condense and
be separated. Ratio will be 75/[(760 × 5) − 5] = 0.0201. Hence (0.1111 − 0.0291) =
0.0910 mole of benzene is condensed per mole of inert.
For 2.3 mol of inert, benzene separated is 0.0910 × 2.3 × 78.1 = 16.35 lb.
REFERENCE
Schwartz, S.E. and Warneck, P. Units for use in atmospheric chemistry. Pure &
Applied Chemistry 67: 1377–1406, 1995. Retrieved on March 9, 2011, http://
www.iupac.org/publications/pac/1995/pdf/6708x1377.pdf.
Section II
Fundamentals and
Problem-Solving Profile
Chapter 3 Numerical Methods and Chemical Engineering Computations
Chapter 4 The Approach to Solve Problems
Section II encompasses Chapters 3 and 4. Chapter 3, “Numerical Methods and
Chemical Engineering Computations,” is devoted to the introduction of numerical
methods, which are often divided into two main categories:
One of the main objectives of Chapters 3 and 4 (Section II) is to devise, evaluate, and
apply numerical techniques for employing computers to solve problems in chemi-
cal engineering. In Chapter 3, numerical methods are introduced with emphasis on
elementary applications. Definitions of basic concepts are given; areas and domain
of numerical methods are identified; and basic applications for many examples are
presented.
3.1 INTRODUCTION
Computer scientists, programmers, and technology support staffs take a very dif-
ferent view of the computer from other users. They must know the details of how
operating systems work and how to code various scripts that control function. They
must be able to control the low-level details that a user simply assumes.
The focus of numerical methods is simply translating engineering problems into
algorithms and implementing them using computers. Numerical method is a branch
of applied mathematics, concerned with methods for solving complicated equa-
tions using arithmetic operations, often so complex that they require a computer, to
approximate the analysis process. Numerical analysis naturally finds applications in
all fields of engineering and the physical sciences. It is an important underpinning
for techniques used in computational science in general. In other words, using com-
puters to solve problems is automatically tied up to numerical techniques.
In the development of numerical methods, simplifications need to be made to
progress toward a solution. For example, general functions may need to be approxi-
mated by polynomials, and computers, on the other hand, cannot generally represent
numbers exactly. Therefore, numerical methods do not usually give the exact answer
for a given problem at hand. However, it could be stated that the overall goal of the
field of numerical analysis is the design and analysis of techniques to give approxi-
mate but accurate solutions to problems.
35
36 Chemical Engineering Primer with Computer Applications
3.2.1 algorIthms
According to Blass and Gurevich (2003), an algorithm is defined as an effec-
tive method expressed as a finite list of well-defined instructions for calculating a
function. Starting from an initial state and initial input, the instructions describe a
computation that, when executed, proceeds through a finite number of well-defined
successive states, eventually producing output and terminating at a final ending state.
In computer systems, an algorithm is basically an instance of logic written in
software by software developers to be effective for the intended target computer(s)
for the target machines to produce output from given input (perhaps null).
Algorithms resemble recipes, which tell you how to accomplish a task by per-
forming a number of steps. For example, to bake a cake the steps are as follows:
preheat the oven; mix flour, sugar, and eggs thoroughly; pour into a baking pan; and
so forth.
Algorithms were originally born as part of mathematics. The word algorithm
comes from the Arabic writer Muhammad ibn Mūsā al-Khwārizmi, but currently the
word is strongly associated with computer science.
components and f(x) a vector with m components. The most important case
is where n = m.
3. Trial and error procedures: This refers to solution of equations using the
iterative convergence methods (Newton’s method).
4. Problems that require numerical differentiation/integration.
5. Optimization: This refers to minimizing or maximizing a real-valued func-
tion f(x). The permitted values for x = (x1,…,xn) can be either constrained
or unconstrained. The linear programming problem is a well-known and
important case; f(x) is linear, and there are linear equality and/or inequality
constraints on x.
6. Use computable functions p(x) to approximate the values of functions f(x)
that are not easily computable or use approximations to simplify dealing
with such functions. The most popular types of computable functions p(x)
are polynomials, rational functions, and others.
7. Regression and correlation of a given set of data.
Algorithm
• Set max to 0.
• For each number x in the list L, compare it to max. If x is larger, set max
to x.
• max is now set to the largest number in the list.
• end
• max
• else
• disp(‘insert a correct set of numbers!’)
• end
q = integral (fun,0,Inf)
q = 1.9475
Solve the following system of equations for x, y, and z by using the matrix method:
10x + 3y + 5z = 8
8x + 2y + 2z = 7
3x + 4y - z = 6
Solution
A = [10, 3, 5; 8, 2, 2; 3, 4, - 1];
B = [8; 7; 6];
X = A\B;
x = 0.7917
y = 0.7917
z = -0.4583
40 Chemical Engineering Primer with Computer Applications
Example 3.6
The experimental values calculated for the heat capacity of ammonia from 0 to
500 are as follows:
Cp = a + bT + CT2 + DT3
where T is in °C
(b) Calculate the amount of heat Q required to increase the temperature of
150 mol/h of ammonia vapor from 0°C to 200°C if you know that
Tout
Q=n
ò Cp dt
Tin
Solution
T =[0,18,25,100,200,300,400,500]
Cp =[8.371,8.472,8.5149
, .035,9.82410
, .60611
, .347,12.045]
P = Polyfit(T,Cp,3)
n = 150; syms t
Q = n*int(Cpf,0.200)
2.7180e+005
é5 1 8ù éx ù é 46ù
ê ú ê ú ê ú
with A = ê 4 -2 0ú , X = êy ú , B = ê12 ú
êë6 7 4úû êë z úû êë50 úû
If A−1 (the inverse of A) exists, we can multiply both sides by A−1 to obtain X =
A−1B. To solve this system of linear equations in Excel, execute the following steps:
Solution
log(initial con + min con)
n=
log[1 + lost% / 100]
min _ con = 50
lost = 0.01
n = 59
Polonium has a lifetime of 140 days. Starting with 10 g today, calculate how much
is left after 250 days?
Solution
Amount left = initial amount (0.5)time half time
≫ initial _ amount = 10;
half _ life = 140; time = 250 ;
Example 3.10
Solve the following system of equations for x, y, and z by using the matrix method
10x + 3y + 5z = 8
8x + 2y + 2z = 7
3x + 4y - z = 6
Solution
x = 0.7917
y = 0.7917
z = -0.4583
Example 3.11
Xwater
= 0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000
0.7000 0.8000 0.9000 1.0000
Xethanol
= 1.0000 0.9000 0.8000 0.7000 0.6000 0.5000 0.4000
0.3000 0.2000 0.1000 0
Answer : rav
= 780 802 824 846 868 890 912
934 956 978 1000
REFERENCES
Blass, A. and Gurevich, Y. Algorithms: A quest for absolute definitions. Bulletin of
European Association for Theoretical Computer Science 81. http://research.microsoft.
com/~gurevich/Opera/164.pdf, 2003. Retrieved July 2012.
Heilbron, J.L. The Dilemmas of an Upright Man. Max Planck and the Fortunes of German
Science, With a New Afterword, Harvard University Press, 2000.
4 The Approach to Solve
Problems by Computers
An experiment is a question which science poses to Nature and a measurement is the
recording of Nature’s answer.
(Planck, 1903)
The main objective of this chapter is to learn how to set up and validate mathemati-
cal models in order to solve chemical engineering problems. The implementation of
models in a structured programming language such as MATLAB® or using spread-
sheets is presented.
Model process development and mathematical formulation is detailed in this
chapter. Applications covering numerous examples are also included.
4.1 INTRODUCTION
Learning how to solve problems is an important part of developing competency in
science and engineering. It is worth noting that most engineering problems are based
upon one of three underlying principles:
45
46 Chemical Engineering Primer with Computer Applications
Process Design
calculations calculations
include two main and
categories computation
Physical Chemical
operations: reactions:
e.g., evaporat., use mole
distillat. balance
Physical Chemical
operations processes
• Find out, for the problem as identified in the first procedure, if you need
additional help from the sources available at your fingertip (found in the
text).
4.2 METHODOLOGY
One of our main objectives in this primer is to demonstrate how to solve chemical
engineering problems that require numerical methods by using standard algorithms,
such as MATLAB or spreadsheets. To say it in simple words, it is to device and
evaluate numerical techniques for employing computers to solve problems in chemi-
cal engineering.
Problem solving using computers could be handled by using
A spreadsheet such as Excel is a program that let you analyze moderately large
amounts of data by placing each data point in a cell and then perform the same
operation on groups of cells at once. One of the advantages of spreadsheets is that
data input and manipulation are relatively intuitive and hence easier than doing the
same tasks in MATLAB (Towler and Sinnott, 2013).
The Approach to Solve Problems by Computers 47
Details on using Excel showing how to input and manipulate data and perform
operations and others are fully explained in many references found in the open
literature.
Spreadsheet software has become indispensable tools in solving chemical engi-
neering problems because of the availability of personal computers, ease of use, and
adaptability to many types of problems. Handheld calculations are encouraged as
well to learn a numerical technique in solving a problem.
The following basic four steps are normally involved in solving problems by
computers:
1. Go to refresher in
2. Formulation of the
Start Text for Ch.E.
Problem: Given, Find?
Fundamentals
3. Selection of an
Algorithm to solve your
problem
4. Computation
End
Formulation step
may result in:
Computer
Input Output: Solution
These steps are presented by Figure 4.1, which illustrates different options of numer-
ical methods involved in solving problems.
Ax = b Linear
X= A–1b Other
Conceptual model
for a physical event or a Mathematical ODE
stated principle model
Nonlinear
4.4 APPLICATIONS
To demonstrate the earlier procedure, the problem of flash distillation is considered.
Flash calculations are very common, perhaps one of the most common chemical
engineering calculations. They are a key component of simulation packages like
Hysys and Aspen.
F
xF
L
x
Now, mathematical formulation is pursued through these two equivalent definitions for
the BP temperature. But let us first present the following fundamental relationships:
p i = Pi0 x i (4.1)
• Dalton’s law: In the vapor phase, where the vapor is in equilibrium with the
liquid, the partial pressure of component, i, is
p i = PT y i (4.2)
åp = P i T
The Approach to Solve Problems by Computers 51
where
pi is the partial pressure of component i
Pi0 is the vapor pressure of pure component i
PT is the total pressure
y i Pi0
= = Ki (4.3)
x i PT
where
A, B, and C are Antoine coefficients that vary from substance to substance
P is the vapor pressure of the pure component
f (Tassu ) = éê
ë å k x ùûú -1 goes to zero at T
i i assu = TBP (4.4)
At the BP, the sum of partial pressure of the components should be equal to the
total pressure on the system, P T, or f (Tassu ) = {[ å p i ] - PT} goes to zero at Tassu = T BP.
YES
Calculate:
YES
The solution of these nonlinear algebraic equations for the BP, Equations 4.4 and 4.5,
using MATLAB, is presented in the following examples.
Similarly, the DP temperature is physically defined as the temperature at which
the first liquid drop would form when the temperature of a mixture of vapors is
slowly decreased (cooled) at a specified constant pressure.
Example 4.1
and
dKi
f ¢(100) = å xi = 0 .0 2
dT
ì f (T ) ü
T2 = T1 – í 1 ý
î f ¢(T1) þ
æ 1 ö
= 100 – ç ÷
è 0.02 ø
= 50°F
B C D E F G H I J K L
T1 Component Xi Ki XiKi [dKi/dT]Tn (xi)(dki/dT)
100 1 1/3 1 1/3 0.01 1/300
2 1/3 2 2/3 0.02 1/150
3 1/3 3 1 0.03 1/100
∑ 6 1 6 2 0.06 0.02
f(100) = [∑KiXi] – 1 1
f΄(100) = ∑ XidKi/dT 0.02
0
T2 = T1–[f(T1)/f΄(T1)] 50 F
54 Chemical Engineering Primer with Computer Applications
Example 4.2
A = [13.8858, 14.0045]
B = [2788.51, 3279.47]
C = [220.79, 213.201]
T = 100
Y = [0.5,1- 0.5]
= 180.0377 34.2488
The total pressure and the liquid mole fraction are calculated:
P = 1/sum (y./Psat )
x = P*y./Psat
P = 57.5498
x = 0.1598 0.8402
Example 4.3
35 kg B
Column 1
Column 2
n1 n5
50 kg T
0.673 kg B/kg 0.059 kg B/kg
15 kg X
0.306 kg T/kg 0.926 kg T/kg
0.021 kg X/kg 0.015 kg X/kg
n6 kg B
n2 kg B n7 kg T: 10% of T in Column 1
n3 kg T n8 kg X: 90% of X in Column 1
n4 kg X
Column 1 Balances
Column 2 Balances
B: n2 = 0.059n5 + n6 (4.3.4)
T: n3 = 0.926n5 + n7 (4.3.5)
X: n4 = 0.015n5 + n8 (4.3.6)
Example 4.4
The solution procedure is to choose a temperature, evaluate PA* and PB* for that tem-
perature from the Antoine equation, evaluate f(Tbp) from the earlier equation, and
repeat the calculations until a temperature is found for which f(Tbp) is sufficiently
close to 0.
Solve using initial guess (100°C):
1203.5311 1346.773
6.89272- 6.95805-
0.40 × 10 T + 219.888 + 0.6 × 10 T + 219.693 - 760 = 0
T = fzero(@(T)040*10^89,272 - 12035
, ,311/(T+219,888))
T = 95.1460
pA = 040*10^89,272 12035
, ,311/(T+219,888))
pA = 472.5616
pB = 06*10^95,805 1,346,773/(T+219,693))
pB = 287.4384
yA = pA/760
yA = 0.6218
yb = pB/760
yb = 0.3782
y APV yBP
+ + = 1
p*A (Tdp ) pB* (Tdp )
Example 4.5
For the mixer shown in the following, write a code to find the values of streams
A, B, and C.
W = 100 kg/h
100% benzene
A=? B=?
50% xylene 30% xylene
20% toluene 30% toluene
30% benzene 40% benzene
C=?
40% xylene
20% toluene
40% benzene
Solution
By making a component-material balance on each component within the mixer,
you can reach to a system of three equations that can be solved by using the com-
mand to find the unknowns A, B, and C.
Type the following command:
[A,B,C] = solve(‘.5* A + .3*B = .4*C’, ‘.2*A + .3*B = .2*C’,‘.3*A +
.4*B + 100 = .4*C’). The results will be the following:
A = 600.
B = 200.
C = 900.
A = 600.
B = 200.
C = 900.
W = 100.
58 Chemical Engineering Primer with Computer Applications
Example 4.6
For the following distillation column, calculate the values of F1, F3, and F4.
F3 = ?
50% A
F1 = ? 30% B
20% A 20% D
30% B
50% D
F2 = 500
50% A
50% B
F4 = ?
20% A
40% B
40% D
For the vapor liquid separater shown in the figure, write a program to calculate the values of
XA, XB, YA, YB, L, and V
If you know: V=?
YA = ?
XA + XB = 1
YA + YB = 1 YB = ? F = 100 kg/h
Vapor
YA = KA*XA = 2XA 40% A
YB = KB*XB = 2XB
Liquid 60% B
L=?
XA = ?
XB = ?
Solution:
A=[1,1,0,0;0,0,1,1;–2,0,1,0;0,–.5,0,1];
B=[1;1;0;0]; X=A\B;
xa=X(1),xb=X(2),ya=X(3),yb=X(4)
a=[xa,ya;xb,yb];
b=[.4*100;.6*100];
x=a\b;
L=x(1),V=x(2)
The Approach to Solve Problems by Computers 59
Plot the Txy diagram for ethanol–water system, knowing the vapor pressure for
three components is calculated by
Solution:
Xe = 0:.1:1;
Xw = 1–Xe;
for m = 1:11
for T = 273.15:.01:450;
375
370
Temperature (K)
365
360
355
350
345
340
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x, y for ethanol
60 Chemical Engineering Primer with Computer Applications
The main five equations relating the seven variables x, T, p*a , p*b , pa, pb, and y are
given in the following.
Equations 4.9.1 and 4.9.2 are the Antoine equations for benzene (a) and toluene
(b), respectively. These give the vapor pressure pi* in mm Hg of the pure liquids,
a or b versus temperature T in °C. Here, xa is the mole fraction of benzene in the
liquid, ya is its mole fraction in the vapor, and pi is the partial pressure of i for a or
b. Equations 4.9.3 and 4.9.4 come from assuming an ideal liquid mixture, which is
approximately true for benzene and toluene.
Find
The goal is to create a Txy phase diagram for mixtures of benzene and toluene,
where T is the temperature, x is the mole fraction of benzene in the liquid, and y
is the mole fraction of benzene in the vapor. A horizontal line drawn for a given T
gives the compositions of liquid and vapor in equilibrium at that T. Such diagrams
are very useful for distillation calculations
Approach: Using these equations, you will create a Txy graph for benzene–toluene
mixtures at 1 atm. This graph will show the mole fraction ya of benzene in the
vapor corresponding to equilibrium with liquid of mole fraction xa at temperature
T required to give a total pressure of 760 mm Hg:
1090
log10 p*a = 6.814 - (4.9.1)
197.1+ T
1457
log10 p*b = 7.136 - (4.9.2)
231.8 + T
pa = xp*a (4.9.3)
pa
y= (4.9.5)
pa + pb
Proposed Procedure
FIGURE 4.4 Calculation of the equilibrium data for Benzene/Toluene mixture using Excel.
(Courtesy of Clarkson University, Engineering spreadsheets Using Excel II. Goalseek and
Graphing, Potsdam, NY, n.d., http://clarkson.edu/~wwilcox/ES100/xl-tut2.pdf, accessed
October 14, 2014.)
62 Chemical Engineering Primer with Computer Applications
115
Calculated Txy diagram for benzene and toluene at 1 atm
110
105
Vapor
Temperature (°C)
100
Vapor + liquid
95
90 Liquid
85
80
75
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Mole fraction benzene
Example 4.10
Given
Find: Calculate the vapor and liquid compositions in equilibrium (in mole frac-
tions, y and x) for the benzene–toluene system, and plot the Txy diagram.
TABLE 4.1
Temperature–Vapor Pressure Data for Benzene–Toluene
at 101.32 kPa
Vapor Pressure
Temperature Benzene Toluene
K °C kPa mm Hg kPa mm Hg
353.3 80.1 101.32 760 — —
358.2 85.0 116.9 877 46.0 345
363.2 90.0 135.5 1016 54.0 405
368.2 95.0 155.7 1168 63.3 475
373.2 100.0 179.2 1344 74.3 557
378.2 105.0 204.2 1532 86.0 645
383.3 110.6 240.0 1800 101.32 760
The Approach to Solve Problems by Computers 63
Solution
Step 1
Write Raoult’s law for a binary mixture of A and B:
PT - PB*
xA =
PA* - PB*
æ P* ö
y A = ç A ÷ xA
ç PT ÷
è ø
Step 2
The vapor-liquid-equilibrium (VLE) data can be calculated at each temperature by
substituting for total pressure (101.32 kPa) and the appropriate vapor pressures.
For example, at 85.0°C,
101.32 - 46.0
xA = = 0.7803
116.9 - 46.0
æ 116.9 ö
yA = ç ÷ (0.7803)
è 101.32 ø
Step 3
Repeating for other temperatures using Excel, the following results are obtained
110
105
y
Temperature (°C)
100
95 x
90
85
80
0.0 0.2 0.4 0.6 0.8 1.0
Mole fraction of benzene in liquid (x) and vapor (y)
Txy diagram for benzene-toluene system at 1 atm
Example 4.11
Plot Pxy diagram for the binary system of acetonitrile (1)/nitromethane (2). Vapor
pressures for the pure species are given by the following equations:
For acetonitrile (1), P1o = exp(14.2724 - 2945.47/(T + 224)).
For nitromethane (2), P2o = exp(14.2043 - 2972.64/(T + 209)),
in which
T = 75°C
P1 and P2 in kPa
Pio
Ki =
Pt
At bubble point, å yi = å Ki ´ xi = 1
Solution
Write the following code:
X1 = 0:.1:1;
X2 = 1−X1;
T = 75;
P1 = exp(14.2724 - 2945.47/(T + 224));% Vapor pressure of
acetonitrile
P2 = exp(14.2043 - 2972.64/(T + 209));% Vapor pressure of
nitromethane for m = 1:11
for Pt = .1:.01:120;
K1 = P1/Pt; % acetonitrile
K2 = P2/Pt; % nitromethane
The Approach to Solve Problems by Computers 65
90
80
70
Pt (kPa)
60
50
40
30
20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x, y for acetonitrile
REFERENCES
Clarkson University, n.d., Engineering Spreadsheets Using Excel II. Goalseek and Graphing,
Potsdam, NY. http://clarkson.edu/~wwilcox/ES100/xl-tut2.pdf (accessed October 14,
2014).
Planck, M. Treatise on Thermodynamics, A. Ogg (Trans.), London: Longmans, Green & Co.
https://archive.org/details/treatiseonthermo00planrich, 1903.
Towler, G. and Sinnott, R.K. Chemical Engineering Design: Principles, Practice and
Economics of Plant and Process Design, New York: Elsevier, 2013.
Section III
Backbone Materials
Fluid
Fluid Basic definitions
statics
mechanics is of fluids
classified Fluid
into: dynamic
1. Compressible and
noncompressible.
2. Pressure concept
Classification of
3. Newton’s second law
fluid flow states that: ∑ F = m . a
4. Poise is viscosity unit
5. Kinematic viscosity:
v = µ/ρ , Stoke is the
unit
Ideal and Steady and Uniform and Laminar and
real fluid unsteady nonuniform turbulent
flow state flow
3. Absolute pressure = gauge pressure (taken with the proper sign) + atmospheric pressure.
4. Fluid statics is concerned with the static properties and behavior of fluids
Bernoulli’s equation
1. The diameter of a pipe is a function of the flow rate of the fluid: D = f (Q), while
2. The thickness of a pipe (t) is a function of the working pressure inside the pipe:
t = f (p)
5.1 INTRODUCTION
Before a formal definition of fluid is given, it is known that the three common states
of matter are solid, liquid, and gas. A fluid is either a liquid or a gas. A fluid is defined
as a substance that deforms continuously under the application of a shear stress.
The definition of a stress would illustrate some basic principles. A stress is defined
as a force per unit area, acting on an infinitesimal surface element. Stresses have
69
70 Chemical Engineering Primer with Computer Applications
both magnitude (force per unit area) and direction, and the direction is relative to the
surface on which the stress acts. Stresses are classified into two main types: normal
stresses and tangential stresses. Pressure is an example of a normal stress and acts
inward, toward the surface, and perpendicular to the surface; on the other hand, a
shear stress is an example of a tangential stress, that is, it acts along the surface,
parallel to the surface. Friction due to fluid viscosity is the primary source of shear
stresses in a fluid.
Fluid mechanics has two branches, as shown next:
Fluid mechanics
Fluid Fluid
statics dynamics
Fluid statics treats fluid in the equilibrium state (no motion), while fluid dynam-
ics treats fluids when portion of the fluid are in motion (concerned with the relation
between the fluid velocity and the forces acting on it).
åF = m ×a
Fluid Flow and Transport of Fluids 71
y τ = F/A u
F du
=τ=η
du/dy A dy
This law is applied for fluid statics and fluid dynamics, while for fluid statics,
å F = 0.0
5.2.4 nature of fluIds
As explained earlier, a fluid undergoes continuous deformation when subjected to a
shear stress. The resistance offered by a real fluid to such deformation is called its
consistency. Now, if the static pressure and temperature are fixed, the consistency of
the fluid is constant. This is defined as a Newtonian fluid (simple fluids and gases), and
consistency is called viscosity. If, however, consistency is a function of shear stress as
well as of temperature, the fluid is called nonnewtonian fluid (see Figure 5.1).
The unit of viscosity (i.e., absolute viscosity) in the cgs system is defined as the poise
(p) = 1 g/(cm)(s). Viscosities are usually tabulated in centipoises: 1 (cp) = 0.01 poise.
Kinematic viscosity of a fluid of density lb/cu × ft and viscosity µ lb/(ft)(s) is given
by
m 2
v= ft /s
r
1 stoke = 1 cm 2 /s
This type of flow is also called potential flow or irrotational flow. This can exist at a
distance not far from a solid boundary (outside the boundary layer).
This type exists inside the boundary layer where the fluid is affected by the presence
of solid boundaries.
• Steady state: In this type of flow, the conditions are independent of time
(invariant with time).
• Unsteady state: The conditions are dependent to time (change with time).
Velocity profile
Inviscid
U(x)
region
Viscous
region
u(x,y)
Boundary
layer
thickness
Viscosity 3
velocity gradient =
wall shear stress
No-slip at a
solid surface
x
Velocity of the fluid is
zero at the plate
By plotting (y) versus (du/dy), the figure shows that the velocity gradient (rate of
shear stress) is maximum at the wall and decreases as the distance increases to reach
the minimum value at the free stream velocity.
gauge pressures. Absolute pressure, on the other hand, is the true total pressure
and is given by
The term static head generally means the pressure of a fluid exerted by the head of a
fluid above the point in question.
Newton’s law, which is given by the equation force = mass × acceleration, is used
to determine the pressure as shown.
For constant density fluids, the static head is given by
hrg
ph =
gc
where
h is the head of liquid above the point, ft
ρ is the liquid density, lb/cu · ft
g is the local acceleration due to gravity, ft/s2
gc is the dimensional constant, 32.17 (lb)(ft)/(lb · force)(s2)
The total amount of energy entering the system is the sum of all of the different types
entering the system. Energy may be transferred between a system and a boundary
by two modes:
1. Heat: Q
2. Work: W
or
where WT is the sum of Ws + ΔPV. This term, ΔPV, is called the flow work that is the
work done on the fluid at the inlet minus work done by the fluid at the outlet.
If we define
DU + DPV = DH ,
DP
DK.E. + DP.E. + + (DU - Q) = -Ws (5.3)
r
This term, ΔU – Q, represents the friction losses in our system. It means that this
energy appears as friction losses, F. Equation 5.3 can be rewritten as
DP
DK.E. + DP.E. + + F + Ws = 0 (5.4)
r
W Q
Vin2 g V2out g
min uin+ Pin vin+ + Z mout uout + Pout vout + + Z
2gc gc in 2gc gc out
System boundary
p1 g u2 p g u2
+ Z1 + 1 = 2 + Z 2 + 2
r g c 2g c r g c 2g c
This equation is the simplified form of Bernoulli’s equation and is used for friction-
less flow; no work is applied on the fluid and no heat is added or lost from the fluid.
Each term in Bernoulli’s equation represents a form of energy possessed by a
moving fluid (potential-, kinetic-, and pressure-related energies). In essence, the
equation physically represents a balance of the KE, PE, and PV energies so that if
one form of energy increases, one or more of the others will decrease to compensate
and vice versa. It is to be noted that each term of Bernoulli’s equation has a unit of
energy per unit mass (J/kg or lbf · ft/lbm).
78 Chemical Engineering Primer with Computer Applications
The units for all the different forms of energy in this equation are measured in units
of distance; these terms are sometimes referred to as heads (pressure head, velocity
head, and elevation head). Each of the energies possessed by a fluid can be expressed
in terms of head.
It has been established that the elevation head represents the PE of a fluid due to
its elevation above a reference level. The velocity head represents the KE of the fluid.
It is the height in feet that a flowing fluid would rise in a column if all of its KE were
converted to PE. The pressure head represents the flow energy of a column of fluid
whose weight is equivalent to the pressure of the fluid.
The sum of the elevation head, velocity head, and pressure head of a fluid is called
the total head.
Thus, Bernoulli’s equation states that the total head of the fluid is constant.
Bernoulli’s equation can be modified to take into account the friction losses in the
fluid flow to be in the form:
V2 V2
åF
p1 g p g
+ Z1 + 1 = 2 + Z 2 + 2 +
r gc 2a1gc r g c 2a 2 g c
where F is the friction loss in the piping system between point 1 and point 2 includ-
ing both the skin and form friction.
If the flow line contains a device that adds work to the fluid, Bernoulli’s equation
will take the form
V2 V2
åF
p1 g p g
+ Z1 + 1 + hW = 2 + Z 2 + 2 +
r gc 2a1gc r g c 2a 2 g c
where
W is the mechanical work done by the pump per unit mass of fluid
η is the pump efficiency
W - friction losses in the pump
h=
W
The relation between the Δpf due to friction and f is given as follows:
éæ DL ö æ u2 ö ù
Dpf = 4fr êç ÷ ç ÷ú
êëè D ø è 2 ø úû
or
Dpf æ DL ö æ u2 ö
Ff = = 4f ç ÷ç ÷
r è D øè 2 ø
16
f=
N Re
This relationship holds up to a Reynolds (Re) number of 2100. Beyond that and for
turbulent flow and for design purposes, Moody chart, shown in Figure 5.4, is used to
predict the value of f and, hence, the frictional pressure drop of round pipes.
Moody diagram
0.1
0.09
0.08 Transition region
0.07 0.05
0.04
0.06 0.03
0.05 0.02
0.015
d
Relative pipe roughness ε
0.04 0.01
Friction factor
0.03 0.005
Material ε (mm)
Concrete, coarse 0.25
Concrete, new smooth 0.025
Drawn rubing 0.0025
Glass, plastic, perspex 0.0025
Iron, cast 0.15
Sewers, old 3.0
Steel, mortar lined 0.1
Steel, rusted 0.5
Steel, structural or forged 0.025
Water mains, old 1.0
5.7.2 PiPing
5.7.2.1 Sizing of Pipelines
Choosing a pipeline size and determining wall thickness represent a major target in the
design of piping system. In this respect, one should be aware of two fundamental concepts
1. The diameter of a pipe is a function of the flow rate of the fluid: D = f(Q)
2. The thickness of a pipe is a function of the working pressure inside the pipe:
t = f(p)
By sizing, we mean to determine the pipe diameter first. An engineer in charge must
specify the diameter of the pipe that will be used in a given piping system. Normally,
the economic factor must be considered in determining the optimum pipe diameter.
To calculate the pipe diameter for noncompressible fluids, one can apply the fol-
lowing well-known equation:
Pipe diameter, d, is readily calculated from this equation for a specified flow rate Q
(bbl/h) and for an assumed fluid velocity u (ft/s).
The first alternative involves a higher capital cost with lower running costs; the sec-
ond involves a lower capital cost with higher running costs specifically because of
the need for more pumps.
So it is necessary to arrive at an economic balance between the two alternatives.
Unfortunately, there are no hard and fast rules or formulas to use; every case is different.
Fluid Flow and Transport of Fluids 81
Total costs
Annual cost, dollars/(year)
Annual variable
cost for pumping
Optimum economic
pipe diameter
o Pipe diameter
The costs of actual pumping equipment undoubtedly must be considered, but the area
in which the pipes will run is also important. For instance, to obtain the same pumping
effort in the desert as opposed to a populated area could involve much higher costs in the
form of providing outside services and even creating a small, self-contained township.
In the flow of oil in pipes, the fixed charges are the cost of the pipe, all fittings,
and installation. All these fixed costs can be related to pipe size to give an approxi-
mate mathematical expression for the sum of the fixed charges.
In the same way, direct costs, or variable costs, comprising mostly the costs of
power for pressure drop plus costs of minor items such as repairs and maintenance,
can be related to pipe size. For a given flow, the power cost decreases as the pipe size
increases. Thus, direct costs decrease with pipe size. And total costs, which include
fixed charges, reach a minimum at some optimum pipe size. The ultimate solution lead-
ing to the optimum economic diameter is found from the graph shown in Figure 5.5.
• Pick up a nominal pipe size with the specified schedule no. that gives Di for
our flow equal or slightly greater than Di obtained before.
• As a final check, use the following equation to calculate the safe working
pressure:
Schedule No = 1000 P/S = 2000(tm /Dav); solving this relationship to obtain
ét ù
Ps = 2Ss ê m ú
ë Dav û
where
P is the operating pressure
Ps is the safe working pressure
Ss is the safe working fiber stress
S is the tensile strength, the greatest longitudinal stress a material can
bear without tearing apart
tm is the minimum thickness of pipe
Dav is the average diameter of Di and Do
There are many factors that affect the pipe wall thickness, which include
5.7.2.4 Relationship between the Pipe Diameter and the Pressure Drops
Two scenarios could be followed: either to assume a value for the velocity u and cal-
culate ΔP or to consider an allowable value for the pressure drop ΔP and calculate
the corresponding u.
The first scenario is illustrated by the following block diagram:
(∆P)cal is ≤ (∆P)spec
Solution
Fluid Flow and Transport of Fluids 83
To determine the fluid velocity in a pipe, the rule-of-thumb economic velocity for
turbulent flow is used, as given next:
• Initial pressure
• Flow characteristics and fluid density
• Ground elevation
• Delivery pressure
5.7.3 PumPs
5.7.3.1 Introduction
A fluid moves through a pipe or a conduit by increasing the pressure of the fluid
using a pump, which supplies the driving force for flow. In doing so, power must be
provided to the pump. There are six basic means that cause the transfer of fluid flow:
gravity, displacement, centrifugal force, electromagnetic force, transfer of momen-
tum, and mechanical impulse. Excluding gravity, centrifugal force is the means most
commonly used today in pumping fluids.
Centrifugal force is applied by using centrifugal pump or compressor, where the
basic function of each is the same. That is to produce KE by the action of centrifugal
force and then converting this KE into potential energy (PE) by efficient reduction of
the velocity of the flowing fluid.
Fluid flow in pipes applying centrifugal devices has in general the following basic
advantages and features:
A concise summary for the comparison between different types of pump is found in
Table 5.1.
Pumps are used for many different applications. Understanding which pump type
one needs for this application is very important. For the oil and gas industry, some
basic features are listed next:
1. Pumps should handle the fluids with low shear and least damage to droplet
sizes causing no emulsions for the effective separation of water from oil.
2. Pumps should be self-priming and experience no gas locking.
3. The requirement of having low net positive suction head is an advantage.
This is advantageous for vessel-emptying applications such as closed-drain
drums or flare knockout drums or any applications encountering high-vapor
pressure liquids.
4. Pumps should handle multiphase fluids.
TABLE 5.1
Comparison between Types of Pumps
Type of Pump Features
Centrifugal Most common; high capacity; discharge lines can be shut off (safe).
Handle liquids with solids.
Reciprocating Low capacity and high head; can handle viscous fluids; used to
discharge bitumen (asphalt) in vacuum distillation columns.
Rotary positive displacement Combination of rotary motion and positive displacement; used in gas
pumps, screw pumps, and metering pumps.
Air displacement Nonmechanical, airlift type; used for acid eggs and jet pumps
Fluid Flow and Transport of Fluids 85
• Solar pumps: Powered by the sun. They can be positive displacement pumps
or centrifugal pumps.
• Water pumps: A type of equipment used to move water through a piping
system. They rely upon principles of displacement, gravity, suction, and
vacuums to move water. They can be either positive displacement or cen-
trifugal pumps.
• Well pumps: Designed to draw water to the surface from an underground
water source. Depending on the well depth and configuration, these pumps
can be jet pumps, centrifugal pumps, or submersible pumps.
[H × × Q] [P × Q]
HP (hydraulic horsepower) = =
3960 1714
HP(hydraulic)
Brake HP (actual) =
a
where
H is the head in ft
is the specific gravity
Q is the flow rate (gpm)
P is the gauge pressure in lb/in.2
α is the pump efficiency; 60% is used for centrifugal pumps
Assume frictionless flow in a long, horizontal, conical pipe. The diameter is 2.0 ft
at one end and 4.0 ft at the other. The pressure head at the smaller end is 16 ft of
water. If water flows through this cone at a rate of 125.6 ft3/s, find the pressure
head at the larger end.
Fluid Flow and Transport of Fluids 87
Solution
By applying Bernoulli’s equation between the two ends:
p1 g c u2 p g u2
+ Z1 + 1 = 2 c + Z 2 + 2
r g 2g r g 2g
V o 125.6 ft 3 /s
V1 = = = 40 ft /s
A1 p(1)2 ft 2
V o 125.6 ft 3 /s
V2 = = = 10 ft /s
A2 p(2)2 ft 2
(40)2 p g (10)2
16 + 0.0 + = 2 c + 0 .0 +
2´ 32.17 r g 2´ 32.17
p2 g c
\ Outlet pressure head = = 39.9 ft
r g
Example 5.2
Water is pumped from a large reservoir to a point 65 ft higher than the reser-
voir. How many feet of head must be added by the pump if 8000 lbm/h flows
through a 6 in. pipe and the frictional head loss is 2 ft? The density of the fluid is
62.4 lbm/ft3 and the pump efficiency is 60%. Assume the kinetic energy correction
factor equals 1.
Solution
V2 V2
åF
p1 g p g
+ Z1 + 1 + hW = 2 + Z 2 + 2 +
r g c 2a1g c r g c 2a2g c
p1 g c V2 p g V2
+ Z1 + 1 + hHp = 2 c + Z 2 + 2 +Hf
r g 2a1g r g 2a2g
where
Hp is the pump head
Hf is the frictional head loss
88 Chemical Engineering Primer with Computer Applications
To use the modified form of Bernoulli’s equation, reference points are chosen
at the surface of the reservoir (point 1) and at the outlet of the pipe (point 2).
The pressure at the surface of the reservoir is the same as the pressure at the
exit of the pipe, that is, atmospheric pressure. The velocity at point 1 will be
essentially zero.
Manual solution is presented first, followed by solving the problem by Excel.
((0.178)2 - 0.0)
0.6 ´ Hp = 65 + + 0 .0 + 2
2 ´ 32.17
Hp = 111.66 ft
Fluid Flow and Transport of Fluids 89
Example 5.3
Calculate the diameter of a pipeline handling 10,000 bbl of oil per hour, assuming
that the velocity of flow is about 5 ft/s.
Solution
Q = u(p / 4)d2
Example 5.4
Find the HP for a pump that is handling 500 gpm of oil against a 1000 ft pipeline
with 19 ft equivalent to pipe fittings and valves. Assume friction losses account to
20% of the total head and the specific gravity of the oil is 0.8.
Solution
Total equivalent H = 1000 + 19 + 0.2(1019) = 1223 ft
HP = [(1223)(500)(0.8)]/3960
= 123.5
Brake HB = 123.5/ 0.6
= 206
Example 5.5
1. Manual solution
The pressures at two sections of a horizontal pipe, Sections 1 and 2, are
0.3 and 0.6 kgf/cm2, and the diameters are 7.5 and 15 cm, respectively.
Determine the direction of flow if water flows at a rate of 8.5 kg/s. State
your assumptions.
Data:
Section 1 Section 2
Equations
Equation of continuity:
Bernoulli’s equation
For the flow direction from 1 to 2,
p1 v12 p v2
+ + z1 = 2 + 2 + z2 + h + w - q
r1g 2g r2g 2g
Section 1 Section 2
-3 -3
V1 = 8.5 ´ 10 /(pD / 4) = 1.924 m/s
2
1 V1 = 8.5 ´ 10 /(pD22 / 4) = 0.481 m/s
P1 = 0.3 kgf/cm2 = 2.9436 × 104 N/m2 P2 = 0.6 kgf/cm2 = 5.8872 × 104 N/m2
29.436 + 1851
. = 58.872 + 0.116 + h + w - q
In our problem the work done by fluid (w) and pump work on fluid (q)
are zero.
So to balance this equation, the quantity h has to have negative val-
ues. This is not possible. This equation will be a correct one if the flow
is from 2 to 1.
Therefore, the flow direction is from the end at which pressure is 0.6 kgf/cm2
and diameter is 15 mm to the end at which pressure is 0.3 kgf/cm2 and
diameter is 7.5 mm.
2. Excel solution
Given
P1 = 0.3 kgf/cm2 P1 = 29,420.4178 N/m2
P2 = 0.6 kgf/cm2 P2 = 58,840.8355 N/m2
D1= 7.5 cm V1 = 0.0001925 m/s
D2 = 15 cm V2 = 0.0000481 m/s
Water flow rate = 8.5 kg/s Volumetric flow rate = 0.0085 m3/s
p1 v12 p v2
+ + z1 = 2 + 2 + z2 + h + W - q
r1g 2g r2g 2g
Fluid Flow and Transport of Fluids 91
h + w − q = −29.42041775253
w=0
q=0
h = −29.420
The previous equation will be a correct one if the flow is from 2 to 1.
Therefore the flow direction is from the end at which pressure
is 0.6 kgf/cm2 and diameter is 15 mm to the end at which pressure is
0.3 kgf/cm2 and diameter is 7.5 mm.
Example 5.6
where
Vt is the terminal velocity in m/s
g is the acceleration due to gravity
pp is the particle density
Dp is the diameter of the spherical particle in m
CD is the dimensionless drag coefficient
Velocity(1) = vt;
for m = 1:1:20
Re = (Dp*vt*p)/mu;
if Re < 0.1
CD = 24/Re;
elseif Re ≥ 0.1 & Re ≤ 1000
CD = 24*(1 + 014*Re^7)/Re;
elseif Re>1000 & Re ≤ 350,000
CD = 0.44;
elseif Re>350,000
CD = 0.19−8*10^4/Re;
end
vt = ((4*g*(pp‐p)*Dp)/(3*CD*p))^5;
Velocity(m + 1) = vt;
if abs (Velocity(m + 1)−Velocity(m))<.0001
break end
end
Velocity’
ans =
1.0000
0.0538
0.0251
0.0190
0.0170
0.0163
Example 5.7
Calculate the minimum area and diameter of a thickener with a circular basin to
treat
0.1 m3/s of a slurry of a solid concentration of 150 kg/m3. The results of batch set-
tling tests are as follows:
160
140
120
100
80
u
60
40
20
0
0 150 300 450 600 750 900 1050 1200 1350 1500
G
0.02
0.018
0.016
0.014
0.012
0.01
0.008
0.006
0.004
0.002
0
0 150 300 450 600 750 900 1050 1200 1350 1500
G = 0.0154 kg/m2s
REFERENCES
FAO Corporation Document Repository, Water lifting devices, pump classifications. Fao.org.,
http://www.fao.org/docrep/010/ah810e/AH810E05.htm\l5.3.1, retrieved on May 25, 2011.
Hydraulic Institute, http://www.pumps.org/, retrieved on May 25, 2011.
6 Heat Transmission
Convection: A wall
Conduction: Through wall
Convection: wall B
Heat Trans Mechanism in HT.
Exchangers
6.1 INTRODUCTION
Heat transfer, as a topic, is concerned with two main items:
1. Temperature
2. Flow of heat
95
96 Chemical Engineering Primer with Computer Applications
Heat transfer always occurs from a region of high temperature to another region
of lower temperature. Heat transfer changes the internal energy of both systems
involved according to the first law of thermodynamics. Heat transfer is a process
function (or path function), as opposed to functions of state; therefore, the amount
of heat transferred in a thermodynamic process that changes the state of a system
depends on how that process occurs and not only on the net difference between the
initial and final states of the process.
In many practical problems, these three mechanisms combine to generate the total
energy flow.
Fourier’s law is the fundamental differential equation for heat transfer by conduction:
dq dT
= Q = - kA (6.1)
dq dx
where
dq/dθ is the rate of flow of heat
Q is the quantity per unit time
A is the area at right angles to the direction of heat flow
−dT/dx is the rate change of temperature with a distance (temperature gradient)
k is the thermal conductivity, which is a function of the material through which
heat flows
The equation for the case of one-dimensional steady-state heat conduction, Figure
6.1, is obtained as follows:
Q = kA avg DT / Dx (6.3)
There are several ways to correlate the geometry, physical properties, and tempera-
ture difference of an object with the rate of heat transfer through the object. In heat
transfer by conduction, the most common means of correlation is through Fourier’s
law of conduction. The law, stated by the equation, is applicable most often in objects
98 Chemical Engineering Primer with Computer Applications
T1 T2 < T1
dT
dx
of rectangular or cylindrical form (e.g., pipes and cylinders). Equation 6.3 is used for
the case of those of rectangular form, while Equation 6.4 is designated for those of
cylindrical form:
æ DT ö
Q = - kA ç ÷ (6.4)
è Dr ø
where
Q is the rate of heat transfer (Btu/h)
Q/A is the heat flux (Btu/h · ft2)
A is the cross-sectional area of heat transfer (ft2)
x is the thickness of slab (ft)
r is the thickness of cylindrical wall (ft)
T is the temperature difference (°F)
k is the thermal conductivity of slab (Btu/h · ft · °F)
T1 - T2
Q = kA lm
r2 - r2
where Alm refers to the log mean area = (A2 − A1)/ln (A1/A2).
Q = [T1 - T4 ]/{R A + R B + R C}
Example 6.1
The heat transfer rate through a section of insulating material (as shown in Figure
6.2) that measures 1 ft 2 in cross-sectional area is equal to 1000 Btu/h. The thick-
ness is 1 in. and the thermal conductivity is 0.12 Btu/h · ft · °F.
Solution
= kA æ DT ö
(a) Q ç X÷
èD ø
Solving for ΔT,
æ DX ö
DT = Q ç kA ÷
è ø
1000 (Btu/h)(1/12 ft )
DT =
0.12 Btu/h × ft × °F
100 Chemical Engineering Primer with Computer Applications
Insulating
material
.
Q = 1000 Btu/h
A = 1 ft2
k = 0.12 Btu/h ft °F
∆x
1 in.
q = 1000
k = 0.12
a=1
800
60
700
50
Temperature (°F)
600
500 40
400 30
300 20
200 10
100
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08
Distance (deltaX)
Columns 1 through 6
8.000000000000000 7.166666666666666 6.333333333333332
5.500000000000000 4.666666666666666 3.833333333333333
Columns 7 through 9
3.000000000000000 2.166666666666666 1.333333333333333
≫ plot(x,temperature)
>
Example 6.2
= 0.8(10 / 0.333)
= 24 Btu/h × ft 2
Ceiling
Radiator
Floor
Solid surface
Natural convection Forced convection
of a heat exchanger to the bulk of a fluid being pumped through the heat exchanger
is an example of forced convection.
Heat transfer by convection is more difficult to analyze than heat transfer by con-
duction because no single property of the heat transfer medium, such as thermal
conductivity, can be defined to describe the mechanism. Heat transfer by convection
varies from situation to situation (upon the fluid flow conditions), and it is frequently
coupled with the mode of fluid flow.
In practice, the analysis of heat transfer by convection is treated empirically (by
direct observation). Convection heat transfer is treated empirically because of the
factors that affect the stagnant film thickness such as fluid velocity, fluid viscosity,
heat flux, surface roughness, and the type of flow (single phase/double phase).
Q = hADT
where
Q is the rate of heat transfer (Btu/h)
h is the convective heat transfer coefficient (Btu/ft2 · °F)
A is the surface area for heat transfer (ft2)
ΔT is the temperature difference (°F)
Example 6.3
An uninsulated steam line, 22 ft long, crosses a room. The outer diameter of the
steam line is 18 in. and the outer surface temperature is 280°F. The convective
heat transfer coefficient for the air is 18 Btu/h · ft 2 · °F. Calculate the heat transfer
rate from the pipe into the room if the room temperature is 72°F.
Solution
Q = hADT
= h(2prL)DT
( )
= 18 Btu/h × ft 2 × °F ( 2(3.14)(0.75 ft)(22 ft) ) (280 °F - 72°F)
5
= 3.88 ´ 10 Btu / h
6.5.2 discussion
Recalling the basic equation for the quantity of heat transfer per unit time:
Q = U o ADT
Fluid A
Fluid B
TA
.
. Q
Q T1 T2 TA T1 T2 TB
h1 h2
TB 1 ∆x 1
k
h1 A kA h2 A
∆x
. TA – TB
Q=
1/h1A + ∆x/kA + 1/h2A
.
Q = UoA∆Toverall
1
Uo =
1/h1 + ∆x/k + 1/h2
FIGURE 6.5 Relationships between U and individual conduction and convection terms.
T3 T4
∆r
T2
T1
A1
Tube
A2
FIGURE 6.6 Example for the illustration of the concept of Uo cylindrical geometry.
Thus, there are three processes involved. Each has an associated heat transfer coef-
ficient, cross-sectional area for heat transfer, and temperature difference. The basic
relationships for these three processes can be expressed as follows:
= h1A1 (T1 - T2 )
Q
k
Q = A lm (T2 - T3 )
Dr
Q = h 2 A 2 (T3 - T4 )
If the basis relationship for each process is solved for its associated temperature dif-
ference and substituted into the expression for ΔT, the following relationship results:
æ 1 Dr 1 ö
DTo = Q ç + + ÷
h A
è 1 1 kA lm h 2A2 ø
Q æ A o DrA o Ao ö
DTo = + +
A o çè h1A1 kA lm h 2 A 2 ÷ø
1
Q = A o DTo
æ Ao DrA o Ao ö
çhA + + ÷
è 1 1 kA lm h 2 A 2 ø
or
Q = UADT
where
1
U=
æ Ao DrA o Ao ö
ç h A + kA + h A ÷
è 1 1 lm 2 2 ø
Example 6.4
Calculate the heat rate per foot of a tube from a condenser under the follow-
ing conditions. T = 232°F and the outer diameter of the copper condenser tube
is 0.75 in. with a wall thickness of 0.1 in. Assume the inner convective heat
transfer coefficient is 2000 Btu/h · ft 2 · °F and the thermal conductivity of cop-
per is 200 Btu/h · ft · °F. Also, the outer convective heat transfer coefficient is
1500 Btu/h · ft · °F.
Solution
1
Uo =
(1/h1) + (Dr/k) + (1/h2 )
1
=
(1/ 2000) + (0.1 in./200)(1 ft/12 in.) + (1/1500)
= 827.6 Btu/h × ft 2 × °F
106 Chemical Engineering Primer with Computer Applications
= U A DT
Q o o lm
Q U A DT
= o o lm
L L
= Uo 2prDTlm
æ 1 ft ö
( )
= 827.6 Btu/h × ft 2 × °F (2p)(0.375 in.) ç ÷ (232°F)
è 12 in. ø
= 37, 700 Btu/h × ft
Thermal
radiation X-ray
log λ (m)
Ultra
Infrared
violet
Radio waves γ rays
Visible light
Q = FAT 4 (6.5)
where
Q is the heat transfer rate (Btu/hr)
Φ is the Stefan–Boltzmann constant (0.174 Btu/hr · ft2 · °R4)
A is the surface area (ft2)
T is the temperature (°R)
• Flow arrangement
• Construction
• Degree of compactness
We will elaborate on the classification by flow arrangement. Three types are known:
1. Parallel flow: Both fluids run side by side in the same direction.
2. Countercurrent flow: Similar to parallel flow, except the streams go in the
opposite side to each other.
3. Cross flow: Here, we have the two-flow stream normal to one another.
TA > TB
q
Convection: A wall
Conduction: Through wall
Convection: wall B Fluid A Fluid B
Inlet Inlet
Exit
Th, i Th, i
T
T
Tc, o
Cold Inlet
Cold Tc, i
Tc, i
(a) Distance along heat exchanger, x (b) Distance, x
FIGURE 6.9 Temperature distribution in both flow arrangements. (a) Parallel flow and
(b) counterflow.
Secondary
fluid out
Steam in
Shell
Channel
(end box or header) U-tube bundle
Pass partitions
Condensate
Secondary
out
U-tube fluid in
(T2 - t1 ) - (T1 - t 2 )
DTlm = (6.6)
ln ( (T2 - t1 ) / (T1 - t 2 ) )
Heat Transmission 109
Tube Shell
outlet inlet Baffles
Shell Tube
outlet inlet
Shell and tube
where
T1 is the temperature of hot fluid stream at the entrance
T2 is the temperature of hot fluid stream at the exit
t1 is the temperature of cold fluid stream at the entrance
t2 is the temperature of cold fluid stream at the exit
Q = UADT
Given: Q, U, ΔT
Find: A
Now, it is not that simple to find the area of a heat exchanger.
The thermal design of a shell and tube HE is briefly described by the following
procedure:
In addition, to carry out your design, it would require the following input:
Example 6.5
DT2 - DT1
DTlm =
ln ( DT2 /DT1 )
(280°F - 190°F)
=
ln ( 280°F/190°F )
= 232°F
Example 6.6
Consider a rectangular plate with some heat sources on each edge. Assuming
these heat sources remain constant, the plate will eventually reach an equilib-
rium temperature. We will assume a constant heat source, while the ambient
air is staying at a constant temperature. This would make one side of the cube
to gain a very high temperature and the other five sides a relatively low tem-
perature (room temperature). We are, of course, ignoring any fans or buildup of
heat within the enclosure. The 2D analog to this problem would be some cross
sections cut out from this cube. Consider discretizing our rectangular plate with
a rectangular grid (Figure 6.12). To determine the equilibrium, we may assume
that if the plate is at equilibrium and xi is a grid point not on the boundary,
then the temperature at xi is given by the average of the temperatures of the
four closest grid points to xi. This creates a linear system for xi, which we can
then solve.
This problem corresponds to solving a partial differential equation using finite
differences in two dimensions using MATLAB.
Heat Transmission 111
20°C
x3 x4
20°C 20°C
x1 x2
50°C
FIGURE 6.12 Rectangular plate with different heat sources on each edge.
Hint: The temperature at any point will be given by the average of the four tem-
peratures around it.
Solution
Given:
The temperature at each point (say x1, x2, x3, and x4) is the average of the adjacent
four points:
50+20+x2+x3
x1 =
4
4x1 x2 x3 0x4 = 70
x1 + 4x2 + 0x3 x4 = 70
x1 + 0x2 + 4x3 x4 = 40
0x1 x2 x3 + 4x4 = 40
X = A\B
X = 31.2500
= 31.2500
= 23.7500
= 23.7500
112 Chemical Engineering Primer with Computer Applications
Example 6.7
Water at the rate of 30,000 Ibm/h is heated from 100°F to 130°F in a shell and
tube heat exchanger. On the shell side, one pass is used with water entering as the
working fluid that is heated at a flow rate of 15,000 Ibm/h at 200°F.
Given: Uo = 250 Btu/h · ft 2 · °F and the average water velocity in the ¾ in. 1D tubes
is 1.2 ft/s.
Find: The total area, A; the number of tubes, n; and the length of the tube, L,
assuming countercurrent flow and one tube pass is used
Solution (Manual)
The exit temperature of the hot water is unknown and to be calculated by
Q = mccpDt c = mhcpDt h
= UADTlm
9 ´ 105
A= = 67.3 ft 2
250 ´ 53.6
The cross sectional area of flow, S, in ft 2, is calculated, given velocity in ft/h and
density in Ibm/ft3 and mc = 30.000 lbm/h:
Therefore,
30, 000
S= = 0.111ft 2
62.4 ´ 12
. ´ 3600
where S is equal to the number of tubes (n) multiplied by the flow area per tube,
that is,
S = n × (π/4)(d2)
Therefore, n = 36.
Excel Solution
Calculations
Estimate a relation between the variation of the number of tubes and the fluid velocity.
V n
0.2 217
0.5 87
0.8 54
1.1 40
1.2 36
1.5 29
1.8 24
150 R2 = 0.9903
100
50
0
0.2 0.5 0.8 1.1 1.2 1.5 1.8
Fluid velocity ft/s
Example 6.8
Potential (V) 1 2 3 4 5 6
Resistance (Ω) 15.47 15.63 15.91 16.32 16.83 17.48
Rw
= 0 .004 Tw - 0.092,
R0
114 Chemical Engineering Primer with Computer Applications
where
R0 is the resistance of the wire (equal to 14 Ω) at 273 K
Tw is in K
Estimate the bed temperature and the value of the heat transfer coefficient between
the surface and the bed.
Solution
The heat generation rate by electrical heating is (V2)/R.
The rate of heat dissipation is hA(Tw − TB),
where Tw and TB are the wire and bed temperatures, respectively.
At equilibrium,
But
Rw
= 0.004Tw - 0.092
R0
so that
æ Rw ö
Tw = 250 ç ÷ + 23
è R0 ø
Thus,
æ Rw ö
(V 2 ) = ç 250hARw - (hARw(TB - 23)),
è R0 ÷ø
where Rw is a mean value of Rw noting that the mean cannot be used inside the
bracket in the equation for Tw.
Thus, a plot of (V2) against Rw should yield a line of slope = 250hARw/R0 as given
in the next diagram, from which the value of the slope is found to be 17.4.
Hence,
17.4 ´ 14
h= = 49.2 W/m2K
250 ´ 12 ´ 10 -4 ´ 16.5
The bed temperature is found by the intercept at (V2) = 0, that is, when
Rw = 15.4 W
Thus,
TB = 250(15.4/14) + 23 = 298 K
Area = 12 cm2
R0 at 273 K = 14 Ω
Heat Transmission 115
Rw
= 0 .004Tw - 0.092
R0
so that
Rw
Tw = 250 + 23
R0
Thus,
æ Rw ö
(V 2 ) = ç 250hARw - (hARw(TB - 23))
è R0 ÷ø
40
35
30
25
y = 17.388x – 267.79
20
V2
R2 = 0.9999
15
10
0
15 15.5 16 16.5 17 17.5 18
Resistance, Rw
where Rw is a mean value of Rw noting that the mean cannot be used inside the
bracket in the equation for Tw.
Thus, a plot of (V2) against Rw should yield a line of slope = 250hARw /R0. This
is shown in the figure from which the value of the slope is 17.4
h = 49.21 W/m2K
116 Chemical Engineering Primer with Computer Applications
The bed temperature is found by the intercept at (V2) = 0 that is when Rw = 15.4 Ω
Rw = 15.4 Ω
TB = 298 K
Example 6.9
Write a program to estimate the physical properties for water in the range of tem-
peratures from 273 to 323 K.
The density for water in the range of
Solution
T = 273:.5:323;
%The density
P = 120,092 − 10,056* T + 0001084* T2;
%The conductivity
K = 034 + 9278* 10 −4* T;
%The specific heat
CP = 0015539* (T − 3082)2 + 4180.9;
A = [T′,P′,K′,CP′]
1.E+04
Nusselt number
1.E+03
Pr = 1000
1.E+02 Pr = 100
Pr = 10
Pr = 1
1.E+01 Pr = 0.05
Pr = 0.005
Your value
1.E+00
1.E+02 1.E+03 1.E+04 1.E+05 1.E+06
Reynolds number
REFERENCE
Ribando, R.J. An excel/visual basic for applications (VBA) primer, Computers in
Education Journal, VIII(2): 38–43, April–June 1998 [A version of this article
updated for Excel 2007].
7 Two-Phase Gas–Liquid
Diffusional Operations
Distillation and Absorption
Liquid–gas contacting systems, in general, are utilized for transferring mass, heat,
and momentum between the phases, subject to constraints of physical and chemical
equilibrium.
Emphasis is placed in this chapter on the diffusional transfer of mass between
phases, known as interphase operations. The transfer of mass from one phase to
another is involved in typical liquid–gas mass transfer operations described as
follows:
Distillation and absorption, as two main unit operations, are covered in this chapter.
119
120 Chemical Engineering Primer with Computer Applications
Binary
Atm P
Based on number of Complex
components in feed Multi
Based on
Hi P design pressure
Single stage
Based on nNumber
Vac. of stages:
Multistage
Two streams:
simple Based on Classification,
stream concepts and tech.
Side streams: no.
comple aspects
Interphase operation?
Transfer of mass from one
Change in phase phase to another in distillation
Separation using distillation and absorption operations is an
and absorption requires interphase operation.
Distillation
change in phase.
Design aspects
N = number of trays, D = diameter
In column design:
N = f (degree of
Design of columns based on
separation)
equilibrium-stage concept involves 3
Major steps: D = f (capacity)
1. Collect thermodynamic V/L
equilibrium data
2. Calculate NT (theor.no. of stages) Two physical models are used
in the design analysis of
3. Find efficiency and convert NT to
distillation columns:
NA(actual no. of stages)
Two key factors to be 1. Equilibrium stage model:
considered in column (finite stage contactors)
design are:
Tray efficiency: It 2. Continuous contact model:
1. Operating pressure
describes the extent to (Differential stage contactors)
2. Reflux ratio
which performance of
NA duplicates
performance of NT
(y i / x i ) K i
a= =
(y j / x j ) K j
where
α is the relative volatility of the more volatile component i to the less
volatile j
yi is the vapor–liquid equilibrium concentration of component i in
the vapor phase
xi is the vapor–liquid equilibrium concentration of component i in
the liquid phase
yj is the vapor–liquid equilibrium concentration of component j in
the vapor phase
xj is the vapor–liquid equilibrium concentration of component j in
the liquid phase
yi/xi = K is commonly called the K value or vapor–liquid distribu-
tion ratio of a component i
Fractional
Single
distillation
stage
(rectification)
Simple2 Simple
Flash1 steam
(batch)
(equilibrium) distillation
differential
1 Vapors are kept in intimate contact with the liquid.
2 Vapors are withdrawn as quickly as they are formed to be condensed by a condenser.
Tray n
Ln Vn+1
Two-Phase Gas–Liquid Diffusional Operations 123
N = f (separation)
D = f (capacity )
DH + DK.E. + DP.E. = Q - Ws
DH = Q
124 Chemical Engineering Primer with Computer Applications
This means that the change in enthalpy (Hout – Hin) of 1 lb of mass flowing is equal
to the heat exchanged by the system, per 1 lb of mass flowing.
A heat balance around the distillation column gives
Fh F + Q R = DH D + BH B + Q C
Q R - Q C = [ DH D + BH B ] - FH F
or
Q net = DH
where
∆H is the enthalpy change
ΔK.E. is the change in kinetic energy
ΔP.E. is the potential energy change
Q is the heat added to the system
Ws is the shaft work done by the system.
QC
D
F
QR
Condenser
Reflux
Cooling water drum
Reflux
Pump
Overhead
Distillation column
product
Feed
Distillation plate
Vapor
Reboiler
Steam
Condensate
Bottoms
Bottoms product
liquid
Upflowing vapor
Downflowing liquid
Columns
classification
Function
Number of
Flash Stage components
(equilibrium)
Batch
Binary Multi
Number of
streams Design
pressure
Two Side
streams streams
Atmospheric Under P
Vacuum
Column
internals
Bubble cap
Sieve
Total condenser y1 xD
Vapor
Steam
Bottoms
the upflowing vapors, thereby increasing the efficacy of the distillation tower. The
reflux ratio is an important parameter in column operation. It is normally defined as
the ratio of reflux to distillate (L/D).
The cooling media typically used are air, water, and refrigerant.
Air cooling is the least expensive. With cooling water, process temperatures of
95°F–125°F are possible. Below 95°F, mechanical refrigeration must be used to
achieve the desired condensing temperature. It is the most expensive method.
128 Chemical Engineering Primer with Computer Applications
TABLE 7.1
Types and Features of Distillation Operations
Atmospheric Vacuum
Operation Features Distillation Distillation Pressure Distillation
Application Fractionation of crude Fractionation of Fractionation and/or
oils heavy residues separation of light
(fuel oil) hydrocarbons
Justification Always work near To avoid thermal To allow condensation of
atmospheric pressure decomposition the overhead stream
using cooling water
Extra equipment Steam jet ejectors Stronger thickness for the
(as compared with and condensers vessel shell
atmospheric distillation) to produce and
maintain vacuum
Extra design features Larger diameter Increased number of trays
(as compared with because of higher (N) because separation
atmospheric distillation) vapor flow rate becomes more difficult;
increased reflux ratio
å yi
å xi = = 1.0
Ki
Thus, the DP is a function of the Ki values, that is, DP = f(å yi /Ki ) and [(å yi /Ki ) - 1]
goes to zero at the DP.
Since Ki is a function of the operating pressure, P, we can say that
æ å yi ö
f(P) = ç ÷ -1
è Ki ø
Two-Phase Gas–Liquid Diffusional Operations 129
The method of false position (FP) will be used to perform the iterative calculations.
Now, we can say that
ì (Pk +1 – P(k ) ) ü
Pk + 2 = Pk -1 – f(Pk +1) ´ í k +1 k ý
î[f(P ) – f(P )] þ
To initialize calculations by the FP, two values are assumed for P: 100 and 190 psia.
As an example, calculations at 100 psia are done:
é æ å y i ö ù é æ 0 .9 ö æ 0 . 1 ö ù
f(P) = êç ÷ - 1ú = êç . ÷ + ç 0.675 ÷ ú – 1.0 = -0.395
ëè Ki ø û ëè 197 ø è øû
where, 1.97 and 0.675 are the k values for propylene and 1-butene at 100 psia,
respectively.
Similar iteration for next trial at P = 190 psia would lead to better values of P.
The very final result will give an operating pressure of about 186 psia for the con-
denser outlet. At this pressure, the k values of propylene and 1-butene are 1.178
and 0.424, respectively.
The liquid composition consists of 0.764 propylene and 0.236 1-butene.
A distillation column can only accomplish a desired separation between the two
limits of minimum reflux (MR) and minimum number of trays, as depicted in
Figure. 7.4 Neither of these situations represents real operations.
130 Chemical Engineering Primer with Computer Applications
Condenser
Reflux drum
Enriching Reflux
(rectification)
section
Distillate
Feed
Stripping
section
Heat in
Reboiler
Heat out
Bottoms
Consideration of these two key factors is what we may call the process design
stage of distillation column. On the other hand, the purpose of the mechanical design
is to select the tower internals and column diameter and height. Some of the factors
involved in design calculations include feed load size and properties and the type of
distillation column utilized.
unknown variables, and a design solution procedure is carried out using the
values guessed. The results are compared with the column specifications,
and if they do not match, a new trial is made and the solution is repeated.
1. To establish a standard for measuring the 1. In this model, equilibrium is never allowed to
performance, the ideal stage concept is be achieved between phases at any point.
introduced. The approach to equilibrium Should equilibrium occur anywhere in the
realized in any stage is defined as stage system, this will be equivalent to an ∞ number
efficiency η of trays.
2. The design of tray column would involve three 2. In the design of packed column, the efficiency
basic steps: η, is function of the transport rates of the
(1) VLE data species and the construction of the separating
(2) Calculation of Ntheor device. Ntheor is calculated by a differential
(3) Finding the value of η to calculateNact balance on the column
Methods applied to calculate the number of trays for distillation columns are usually
grouped into the following categories:
7.1.5.2 Assumptions
• Constant molar overflow.
• Molar heat of vaporization of the components is roughly the same.
7.1.5.3 Procedure
1. The equation of the operating line for the rectification section is derived by
making a material balance on two passing streams on a tray, above the feed.
æ L ö æDö
y = ç ÷ x+ ç ÷ x D
èVø èVø
é R ù é 1 ù
=ê ú x + ê (R +1) ú x D
ë ( R +1) û ë û
where
y is the mole fraction of the more volatile component in the vapor phase
x is the mole fraction of the more volatile component in the liquid phase
xD is the mole fraction of the more volatile component in the distillate
D is the distillate product
L is the flow rate of liquid stream in the rectifying section
V is the flow rate of vapor stream in the rectifying section
R is the reflux ratio L/D
1.0
0.8 6
5
0.6
Vapor (y)
0.4 3
2
0.2
1
xF
0.0
0.0 xb 0.2 0.4 0.6 0.8 xD 1.0
Liquid (x)
2. The equation of the operating line for the stripping section is derived using
a similar approach for a tray below the feed.
æLö æBö
y = ç ÷ x – ç ÷ xB
èVø èVø
where
y is the mole fraction of the more volatile component in the vapor phase
x is the mole fraction of the more volatile component in the liquid phase
xB is the mole fraction of the more volatile component in the bottom
product
B is the bottom product
L is the flow rate of liquid stream in the rectifying section
V is the flow rate of vapor stream in the rectifying section
The operating line for the stripping section is then drawn in a similar man-
ner. Here, the starting point is the desired bottom product composition.
A vertical line is drawn from this point to the diagonal line, and a line of
slope L/V is drawn.
3. Given the operating lines for both the stripping and the rectification sec-
tions, the graphical construction of the steps is applied. In this way, we
count the number of steps as shown in the diagram. This particular example
shows that the seven theoretical stages are required to achieve the desired
134 Chemical Engineering Primer with Computer Applications
1.0
q>1
q=1
0.8 0<q<1
0.6
Vapor (y)
q=0
0.4
q<0
ZF
0.2
0.0
0.0 0.2 0.4 0.6 0.8 1.0
Liquid (x)
separation. The required number of trays (as opposed to stages) is one less
than the number of stages since the graphical construction includes the con-
tribution of the reboiler in carrying out the separation.
The actual number of trays required is
N theorit
N actual =
h
0.0
0.0 1.0
Mole fraction of lighter component in liquid phase
Simply enter your parameters at the top of the spreadsheet, including the
The spreadsheet will automatically calculate flow rates throughout the column,
the number of theoretical plates, the feed plate position, the MR ratio (from the
Underwood equation), and the minimum number of theoretical plates (from the
Fenske equation). (Microsoft Office Excel, 2007).
Concept
Thiel and Geddes method is used to calculate the number of trays. It involves
the simultaneous solution of equilibrium relationships (VLE) and the operating
136 Chemical Engineering Primer with Computer Applications
line, where the operating line is used to compute the composition of one of the
two streams passing each other for two consecutive plates, while the equilibrium
relationship is used to compute the composition of either the vapor or liquid
(in equilibrium) on the same plate.
This example is an oversimplified one for illustrative purposes.
Given
é R ù é 1 ù
yn+1 = ê ú Xn + ê ú XD. (7.2.1)
ë R+1û ë R + 1û
where R is the reflux ratio (RR) = L/V and XD is the composition of the
overhead product.
2. Equilibrium data
Yn
Xn = (7.2.2)
[Yn + a(1- Yn )]
A liquid mixture at the boiling point that consists of 70 mole% benzene and
30 mole% toluene is fed to a stripping column. Pressure is taken 1 atm. Feed rate
is 400 kg mole/h. Stripping operation is carried out to achieve a bottom product
W = 60 kg · mol/h that contains no more than 2 mole% benzene.
Solve the problem using Excel, in order to determine the number of the-
oretical trays, N, required to obtain the desired specifications of the bot-
tom product W. Use αAB relative volatility for benzene/toluene, where
a AB = KBenzene / K Toluene = PB0 / PT0 = 2.45
Solution
The number of trays required to reach a bottom product, exiting the stripping col-
umn, is found to be around 11 trays, as seen next; that corresponds to 1.9 mole%
benzene.
138 Chemical Engineering Primer with Computer Applications
The design trade-off between reflux and stages is the standard economic optimi-
zation problem chemical engineers always face: balancing the capital costs for
the number of trays versus the operating cost for the reflux to be recirculated.
A good design will operate near a cost optimum reflux ratio, as illustrated in
Figure 7.7.
The TR condition represents operation with no product removal. The overhead
vapor is condensed totally and returned as reflux. Consequently, the reflux ratio
(L/D) is infinite. The operating lines coincide with the diagonal line.
The aforementioned conditions of TR and MR can be further visualized by the
following approach (with the aid of the opposite diagrams):
B
L=0 B=0
MR TR
The MR condition, on the other hand, represents the theoretical opposite of TR—
an infinite number of ideal separation stages. In this case, the intersection of the
operating lines lies on the equilibrium curve itself.
This intersection is called a pinch point. A simple column will have two pinch
points (because there are two operating lines). The points change when the operat-
ing lines do. A pinch at the intersection of the feed line and the equilibrium curve
indicates that the column is operating at MR.
150
90
Steam and cooling–water
costs (1)
60
One can replace the true equilibrium curve with an effective equilibrium curve
(EEC) located between the true curve and the operating lines as shown in Figure 7.8.
In this case, the number of stages counted will be the true or actual trays.
Murphree efficiency is as indicated in the following equation:
[ Yn - Yn +1 ]
EEC =
[ Yn* - Yn +1 ]
Z1
=
Z2
140 Chemical Engineering Primer with Computer Applications
Equilibrium Effective
curve equilibrium curve
Y*n
Yn
Y*n
Yn n
z1 z2 Yn+1
Yn+1 n+1
Liquid composition, x
The results of the correlation for fractionators and absorbers are graphically pre-
sented by Peters and Timmerhaus. Approximate estimation of efficiency is done by
a number of methods as described by Van Winkle.
1. To allow the contact of the V and L in such a manner, where the interfacial
area and time of contact are at their maximum for the stage to approach
equilibrium
2. To handle the required quantities of V and L without excessive pressure drop
The approach to equilibrium, on the other hand, is a rate factor, intimately tied up
into the rate of mass and heat transfer. A great number of possible variables con-
tribute to these rate processes, which makes the evaluation of efficiency extremely
difficult.
Therefore, once NT (Ntheor) is calculated, estimation of the efficiency is carried out
to find NA (Nact) rather than taking the rate processes to evaluate the efficiency as
indicated in Figure 7.9a.
In other words, to determine NA, one can pursue one of the two options:
Calculation of NT
Go to
Obtain
data on
efficiency
Calculate NA
= NT/E
(a)
Face complex
problems of HT.
Head on approach and M. transfer Find
Start on the trays rates?
NT Find E NA
(b)
FIGURE 7.9 (a) Approach to find efficiency. (b) Head on approach to find the effi-
ciency (N T, NA, and E are theoretical number of trays, actual number, and tray efficiency,
respectively).
142 Chemical Engineering Primer with Computer Applications
where
N is the minimum number of theoretical plates required at TR (of which the
reboiler is one)
Xd is the mole fraction of the more volatile component in the overhead distillate
Xb is the mole fraction of the more volatile component in the bottom products
αavg is the average relative volatility of the more volatile component to the less
volatile component
1. First equation
where
zF is the mole fraction of component n, in the feed
q is the thermal condition of feed (heat required to vaporize 1 mol of
feed divided by molar latent heat of feed)
Two-Phase Gas–Liquid Diffusional Operations 143
2. Second equation
(a1 )(x D1 ) (a 2 )(x D 2 )
Rm + 1 = + +
[(a1 ) - q] [(a 2 ) - q]
The first equation is solved for the value of the θ, which lies between the relative
volatilities of the key components.
Having obtained θ, we calculate Rm from the second equation.
Infinite
1.0
stages
0.8
0.6
N – Nmin
N+1
0.4
0.2
Min 0
stages 0 0.2 0.4 0.6 0.8 1.0
Min R – Rmin Total
reflux reflux
R+1
[ L - V ]
Vm = K v
L
[r L - r v ]
G m = Vmrv = K v
rL
where
Vm is the maximum allowable vapor velocity, based on cross-sectional area of
empty tower, ft/s
Kv is the an empirical constant, ft/s. It is a major function of tray spacing
ρL and ρv are the liquid and vapor densities, lb/ft3, respectively
Gm is the maximum allowable mass velocity of vapor, lb/(s) · (ft2)
Solving the equation for the cross-sectional area A, the tower diameter is readily
obtained.
here, methods used in practice to produce and hence separate two-phase mix-
tures are as follows:
For flashing to take place, the feed has to be two-phase mixture, that is, it satisfies
the following:
TBP < Tf < TDP (as indicated in the VLE diagram, Figure 7.11).
Or,
The sum of (ZiK i) for all components is greater than 1 and the sum of [zi/K i] for
all components is less than 1, where TBP, Tf, and TDP are the bubble point of the feed
mixture, flash temperature, and dew point of the feed mixture, respectively.
Zi and K i are the feed composition and equilibrium constant, respectively, for
component i.
Vapor
T.D.P.
Tf
Va
Temperature
po
r+
T.B.P. liq
uid
Liquid
X, feed
Concentration (mole fraction)
1. Simple form
Zi
Xi =
1 - (V/F )(1 - K i )
For given conditions of P and T, the solution of the equation to find the value
of Xi is obtained by trial and error assuming a value for V/F (take F unity),
until the sum of [Xi] = 1 is satisfied.
2. Functional form
i =c
é ù
å êë1 - g(1 - K ) úû - 1
Zi
f (g ) =
i =1 i
where g = V/F. The desired root to this function, gr = (V/F)r, is the value
that makes the function f(g) go to zero or the sum of [Xi] = 1. The solution
of this flash equation is carried out using computers. Details on using func-
tional technique are fully described in many references.
V, Yi
F, Zi
L, Xi
Yi
Ki =
Xi
where
Yi is the mole fraction of component i in the vapor phase
Xi is the mole fraction of component i in the liquid phase
K i is the equilibrium constant of component i for a given T and P, is considered
the key concept used in the computation of phase behavior of hydrocarbon
mixtures in oil and gas streams
Stripping (desorption) on the other hand, is the transfer of gas, dissolved in a liquid,
into a gas stream. It is the reverse of absorption. The fluid to be processed is liquid;
148 Chemical Engineering Primer with Computer Applications
the aim is to liberate a dissolved gas from the fluid. Stripping can be carried out by
more than one option:
Design methods for absorbers are presented. This includes graphical technique,
numerical (analytical) approach, and the unit concept method.
Contact between
gas and liquid
Stripping is the phases is brought A Typical example is
removal of a solute washing of ammonia
by two means: from a mixture of
from a liquid by
(ammonia+air) by
using heat or inert
liquid water
Gravity Mech.
Solution of
Absorption Design aspects design problems
tskes place by is reached by:
using:
Graphical
methods
Plate Packed Spray
columns towers columns
Analytical
methods
Lean oil
Lean gas
Rich gas Rich oil
1. Graphical methods
2. Analytical methods using absorption factors
3. Transfer unit concept
The following parameters are the most important to consider in the design of
absorbers/strippers:
In absorption problems, there are two carriers and a solute which is partitioned
between them. An example is the use of sulfuric acid to absorb humidity (water)
from moist air. The two carriers are: air and acid, where water is the solute. This is
represented by
Moist
air Water Acid
Equilibrium
(Tower bottom) line
y1
(Top)
Mole fraction, y
y2
Operating
line
Equilibrium
line
y2 Operating
(Top) line
y1 (Tower bottom)
0 0
0 x2 x1 0 x1 x2
FIGURE 7.13 Illustrations for the operating line and equilibrium curve for absorbers/
strippers. (a) Case of absorption. (b) Case of stripping.
Two-Phase Gas–Liquid Diffusional Operations 151
1. The gas stream flow rate G (mol/h), at any point in the tower consists of two
components:
a. Diffusing solute, A, composition y, or Y = y/(1 − y) and GA = G · y; and
b. Nondiffusing (inert gas), composition is (1 − y); Vs = V(1 − y) and (Vs +
VA) = [V(1 − y) + V · y] = V
2. The liquid stream, L mole/h · ft2, at any point, is composed of two species:
a. Nonvolatile liquid, Ls
b. The absorbed solute of the gas, in the liquid phase, having mole frac-
tion: X = x/(1 − x) and Ls = L(1 − x), LA = Lx
Since Vs and Ls are almost constant in quantity, as they pass through the tower, it is
convenient to express the material balance in terms of these flow rates.
Solute balance around the boundary shown in the diagram yields
Vs (Y1 - Y) = L s (X1 - X)
This represents an equation of straight line (operating line):
Y1 – Y =[Ls/Vs] X2 – X , or
where
Y is the moles of solute/moles of inert gas
X is the moles of solute/moles of inert liquid
Finally, we obtain
L [y b - ya ]
=
V [x b - xa ]
Y X
Example 7.4
Moist air that enters is to be dried in a countercurrent drier using sulfuric acid to
absorb the moisture. Conditions are as given next.
Given
Find
Solution
(a) As shown earlier, the slope of the operating line EF is
Y/X = ratio of: H2SO4 Acid to Air treated
= 2.86 × 10 −3, mole acid to mole dry air
= 2.86 mole acid for 1000 moles dry air
(b) Line FH has the lowest slope through F that can provide a value for Yb =
2.8 × 10 −2.
0.032
Yb E H C
ium
ilibr
Equ r ve
0.024 cu
B
Moles dry air
Moles H2O
D
P
0.016 N O
L
J M a b
H K Air, H2O Air, H2O
I
Y=
A Xa Xb
0
0 4 8 12 16 20
Moles H2O
X=
Moles H2SO4
This line is a tangent to the equilibrium curve at point D. It has a slope of 1.64 ×
10 −3 mol acid per mole air. Although this is the smallest amount of acid to be used,
it requires an infinite number of stages as shown in the diagram (pinch point).
Example 7.5
Find
(a) The minimum oil rate (lb/h) required to recover 90% of the entering benzene
(b) The number of theoretical stages needed if an oil rate that is 1.5 times the
minimum is used.
0.05
0.04
Ya, lb . mole benzene/lb . mole air
0.03
0.02
0.01
(X0, Y1)
0
0 0.1 0.2 0.3 0.4 0.5 0.6
Xa, lb . mole benzene/lb . mole absorbent oil
Solution
First. Calculation of the equilibrium data, Ya /Xa, for benzene/air is done.
For benzene,
æP ö
Ya = ç a ÷ xa
èPø
0 .0 5
X 0 = 0 .0 and YN+1 = = 0.0526
0.95
Operating oil rate = (1.5)(Oil rate)m = 15,600 lb/h, and the value of XN = 0.345
The number of stages are figured out to be about five stages for the given oil
rate.
A detailed discussion for the design methods of plate and packed absorbers is
beyond the scope of this chapter. Reference to Perry’s Handbook (50th edition),
pp. 14-14–14-24 is recommended. Heat effects in gas absorption, multicompo-
nent systems, and absorption with chemical reaction are fully discussed as well
on pp. 14-24–14-31.
[y b - ya ]
N=
[ y a - y*a ]
[ xa - x b ]
N=
[ x b - x*b ]
Example 7.6
Calculate:
Solution
(a) Material balance based on 1 mole of air and using mole ratio units is
4
Yb = = 0.04167
96
156 Chemical Engineering Primer with Computer Applications
= 14 trays
Z = (H G )(N G )
where
HG is the height of a transfer unit based on the gas film. The units of H are in
meter
NG is the number of transfer units
When the solutions in gas absorption are dilute (concentrations < 10%), the following
equation is applicable:
é (1 - y)LM ù é (y1 - y 2 ) ù
Z = HGNG = HG ê ú ê ú
ë (1 - y) û av ë (y - y i ) û
where
y1 and y2 are the mole fractions of component A in gas inlet and gas outlet,
respectively
y and yi are the mole fraction in bulk of gas phase and at gas side of interface,
respectively
Sweet NGL
to pipeline
shipping
Cooler
Chemical and physical solvents are the most commonly used absorption pro-
cesses for acid gas treatment. In the chemical process, mono- and diethanolamine
are used to react with the acid gases reversibly and form products that can be gener-
ated by a change in temperature as given in Figure 7.16 The physical process, on
the other hand, uses a solvent to absorb the sour gases (Sulfinol process) from the
natural gas. The absorbed impurities are then separated from the rich solvent by
steam stripping.
Another important application for gas absorption is the separation of the com-
ponents of natural gas liquid, known as C2+, ethane and heavier. This is done using
heavy hydrocarbon cut to absorb these components, which are generated from the
rich solvent by distillation, as shown in Figure 7.17.
Example 7.7
For the following distillation column, write a code to find the value of stream B and
the compositions of stream B.
158 Chemical Engineering Primer with Computer Applications
MATLAB®
Results
Two-Phase Gas–Liquid Diffusional Operations 159
Example 7.8
Xylene, styrene, toluene, and benzene are to be separated with the array of distilla-
tion columns that are shown later. Write a program to calculate the amount of the
streams D, B, D1, B1, D2, and B2 and also to calculate the composition of streams
D and B.
7% xylene
D1
4% styrene
54% toluene
35% benzene
#2
18% xylene
B1 24% styrene
#1 42% toluene
15% xylene 16% benzene
25% styrene
40% toluene
15% xylene
20% benzene D2 10% styrene
54% toluene
F = 70 moles/min
21% benzene
#3
B
B2 24% xylene
65% styrene
10% toluene
1% benzene
Solution
Making material balances on individual components on the overall separation
train yields the following equation set:
Xylene: 0.07D1 + 0.18B1 + 0.15D2 + 0.24B2 = 0.15 × 70
Styrene: 0.04D1 + 0.24B1 + 0.10D2 + 0.65B2 = 0.25 × 70
Toluene: 0.54D1 + 0.42B1 + 0.54D2 + 0.10B2 = 0.40 × 70
Benzene: 0.35D1 + 0.16B1 + 0.21D2 + 0.01B2 = 0.20 × 70
Overall material balances and individual component balances on column 2 can
be used to determine the molar flow rate and mole fractions from the equation of
stream D.
Molar flow rates: D = D1 + B1
Xylene: XDxD = 0.07D1 + 0.18B1
Styrene: XDsD = 0.04D1 + 0.24B1
Toluene: XDtD = 0.54D1 + 0.42B1
Benzene: XDbD = 0.35D1 + 0.16B1
160 Chemical Engineering Primer with Computer Applications
where
XDx is the mole fraction of xylene
XDs is the mole fraction of styrene
XDt is the mole fraction of toluene
XDb is the mole fraction of benzene
Similarly, overall balances and individual component balances on column 3 can
be used to determine the molar flow rate and mole fractions of stream B from the
equation set.
Molar flow rates: B = D2 + B2
Xylene: XBxB = 0.15D2 + 0.24B2
Styrene: XBsB = 0.10D2 + 0.65B2
Toluene: XBtB = 0.54D2 + 0.10B2
Benzene: XBbB = 0.21D2 + 0.01B2
where F, D, B, D1, B1, D2, and B2 are the molar flow rates in mole/min.
Solution
clear all and clc
A =[0.07,.
0 18,.0 15,.
0 24; 0.040
,.240,.10,.
0 65;0.540,.420,.540
,0.1;0.35,.
0 160 0 01];
,.21,.
B=[0.15*70; 0.25*70; 0.4*70; 0.2*70];
X = A\B; D1 = X(1),B1 = X(2),D2= X(3),B2= X(4), D = D1+B1
B = D2+B2
XDx =(.07*D1+.18*B1)/D XDs =(.04*D1+.24*B1)/D
XDt =(.54*D1+.42*B1)/D XDb =(.35*D1+.16*B1)/D
XBx =(.15*D2+.24*B2)/B XBs =(.1*D2+.65*B2)/B
XBt =(.54*D2+.1*B2)/B XBb =(.21*D2+.01*B2)/B
Example 7.9
For the following separation system, we know the inlet mass flow rate (in kg/h) and
the mass fractions of each species in the inlet (stream 1) and each outlet (streams
2, 4, and 5). Calculate the unknown mass flow rates of each outlet stream.
Two-Phase Gas–Liquid Diffusional Operations 161
Solution
If we define the unknowns as x1 = F1, x2 = F2, and x3 = F3 and set up the
mass balances for the following:
F1 = ?
A = 4%
B = 93%
C = 3%
F = 10 F2 = ?
A = 20% A = 54%
B = 60% B = 24%
C = 20% C = 22%
F3 = ?
A = 26%
é 1 1 1 ù é x1ù é10 ù
ê ú ê ú ê ú
ê0.04 0.54 0.26ú êx2ú = ê 2 ú
êë 0.93 0.24 0 úû êë x3úû êë 6 úû
Example 7.10
a) Determine the mass transfer coefficient through the gas film in kg/m2s
(kN/m2).
b) Determine the mass transfer coefficient through the gas film in kg/m2s
(kN/m2) if the flow rate increases from 2 to 10 kg/m2s.
Solution
Given
Flow rate = 2.00 kg/m2s
Reynolds no. = 5100
Friction factor R/pu2 = 0.020
Viscosty = 0.0180 mN · s/m2
Density = 1.154 kg/m3
Diffusion coefficient = 0.000020 m2/s
(a)
[(hd/u)*(PBm/P)*(u/pD)exp0.56]=jd
(u/pD)exp0.56 = 0.88
(hd/u)*(PBm/P)= 0.02261364
u = G/p = 1.73310225
KG = [(hd/R*T)*(PBm/P)] = 1.5319E-05
KG = 0.00026392 kg/m2s (kN/m2)
KG = 2.70 × 10–4 kg/m2s (kN/m2)
(b)
The previous steps were repeated for different ammonia flow rates and
the mass transfer coefficient was calculated in each case.
The results were shown in the following table and were plotted in
the figure.
Flow Rate kg
2 0.0002
4 0.0005
6 0.0008
8 0.00107
10 0.00134
Two-Phase Gas–Liquid Diffusional Operations 163
Example 7.11
Mole fraction of heptane in vapor 0.96 0.91 0.83 0.74 0.65 0.50 0.37 0.24
Mole fraction of heptane in liquid 0.92 0.82 0.69 0.57 0.46 0.32 0.22 0.13
Determine the minimum reflux ratio required. What diameter column would be
required if the reflux used were twice the minimum possible?
164 Chemical Engineering Primer with Computer Applications
heptane in vapour 0.96 0.91 0.83 0.74 0.65 0.5 0.37 0.24 0
heptane in liquid 0.92 0.82 0.69 0.57 0.46 0.32 0.22 0.13 0
Xd = 0.695
100 = (D + W)
and: 100 × 0.695 = 0.99D + 0.01W since 99 per cent n-octane is required.
Hence: D = 69.9 and W = 30.1
and: Ln = 2.16D = 151 and Vn = Ln + D = 221
The mean molecular mass of the feed = (0.695 × 100) + (0.305 × 114) =
104.3 kg/kmole.
Thus: feed rate = (1.25/104.3) = 0.0120 kmole/s
The vapor flow at the top of the column = (221/100) × 0.0120 = 0.0265 kmole/s.
The vapor density at the top of the column = (1/22.4) (273/372) =01
Hence, the volumetric vapor flow = (0.0265/0.0328) = 0.808 m3/s
If the vapor velocity = 0.6 m/s, the are a required = (0.808/0.6) = 1.35 m2 equivalent
to a column diameter of [(4 × 1.35)/π]0.5 = 1.31 m.
REFERENCE
Microsoft Office Excel 2007, Excel 2007 Win32 English VUP CD, http://excelcalculations.
blogspot.com/?, Mindspark Interactive Network, Inc., 2016.
8 Reaction Kinetics,
Chemical Reactors,
and Thermodynamics
Chapter 8 consists of three parts: Part I “Fundamentals of Reaction Kinetics,”
Part II “Types and Design of Chemical Reactors,” and Part III “Introduction to
Thermodynamics.”
In Part I, the two forms of the rate law for chemical kinetics are presented: the
differential rate law and the integrated rate law. Most chemical reactions obey one
of the three differential rate laws:
The main factors that influence the reaction rate are discussed. For each inte-
grated rate law, the characteristic plots are shown and presented in table form as
well.
A summary of the types of laboratory and industrial catalytic reactors is pre-
sented in Part II, along with the equations underlying the design of a chemical
reactor.
Part III, on the other hand, introduces thermodynamics. Thermodynamics applies
to a wide variety of topics in science and engineering, especially physical chem-
istry, chemical engineering, thermal power generation, and industrial applications
of steam and combustion turbines. It is defined as the subject that relates heat to
forces acting between contiguous parts of bodies and the relation of heat to electrical
energy.
165
166 Chemical Engineering Primer with Computer Applications
A simple rule of thumb: for every 10°C increase, the reaction rate doubles.
(1) Based on models and types (2) Based on heating and cooling
Two
fundamental
questions?
One concludes that thermodynamics gave us a partial answer to our problem. Not
only do we have to know whether a reaction is thermodynamically favored, but we also
have to know whether the reaction can or will proceed at a finite rate. This completes the
other part of the chemistry story by the study of the rate of reactions or chemical kinetics.
8 Fe + S8 ® 8 FeS
2H 2 O ® 2H 2 + O 2
4. Single displacement: This is when one element trades places with another
element in a compound.
One example of a single displacement reaction is when magnesium
replaces hydrogen in water to make magnesium hydroxide and hydrogen gas:
Mg + 2H 2O ® Mg(OH)2 + H 2
168 Chemical Engineering Primer with Computer Applications
5. Double displacement: This is when the anions and cations of two different
molecules switch places, forming two entirely different compounds.
One example of a double displacement reaction is the reaction of lead
nitrate with potassium iodide to form lead iodide and potassium nitrate:
The differential rate law relates the rate of reaction to the concentrations of the vari-
ous species in the system. Differential rate laws can take on many different forms,
Reaction Kinetics, Chemical Reactors, and Thermodynamics 169
Each rate law contains a constant, k, called the rate constant. The units for the rate
constant depend upon the rate law, because the rate always has units of mole/L ⋅ s and
the concentration always has units of mole/L.
It should be stated that the rate constant, k, is dependent on the temperature of
which the reaction takes place. This can be seen through the Arrhenius equation
shown as follows:
k = Ae - E / RT
r=k (8.1)
r = k[ A] (8.2)
r = k[ A]2 (8.3)
A®B
d[ A]
r=- (8.4)
dt
r = k[ A], (8.2)
d[ A]
r =- = k[ A] (8.5)
dt
This equation is a differential equation that relates the rate of change in a concentra-
tion to the concentration. By multiplying both sides of Equation 8.5 by dt/[A], one
obtains the following equation:
d[ A]
= - k dt (8.6)
[ A]
Reaction Kinetics, Chemical Reactors, and Thermodynamics 171
Integration of this equation between the time limits, t = o and t = t, gives the cor-
responding integrated rate law (Equation 8.7), which relates the concentration to
time:
[ A] = [ A]0 e - kt (8.7)
Slope of
Reaction Differential Characteristic Kinetic Units of Rate
Order Rate Law Integrated Rate Law Kinetic Plot Plot Constant
d éA ù
Zero - ë û =k [A] = [A]0 − kt [A] vs. t −k mole/L · s
dt
d éA ù
First - ë û = k éë A ùû [A] = [A]0 e−kt ln [A] vs. t −k s−1
dt
d éA ù éë A ùû 0
- ë û = k éë A ùû
2
Second éë A ùû = 1/[A] vs. t k L/mole · s
dt 1 + kt éë A ùû 0
The series of three graphs shown later illustrates the use of the characteristic kinetic
plots. The graph on the left shows [A] versus t plots for a zero-order (red line), first-
order (green line), and second-order (blue line) reaction. The graph in the middle
shows ln [A] versus t plots for each reaction order, and the graph on the right shows
1/[A] versus t plots for each reaction order.
172 Chemical Engineering Primer with Computer Applications
ln[concentration L/mole]
Concentration (mole/L)
Red
–4.00
50.00
Green
–6.00
0.00
0.00 5.00 10.00 0.00 5.00 10.00
Time (s) Time (s)
1000.00
1/concentration vs. time
1/concentration (L/mole)
500.00
Blue
0.00
0.00 5.00 10.00
x
On the other hand, there are two main basic vessel types (in accordance to their
shapes):
1. A tank reactor
2. A pipe or tubular reactor (laminar flow reactor)
Both types can be used as continuous reactors or batch reactors, and either may
accommodate one or more solids (reagents, catalyst, or inert materials), but the
reagents and products are typically fluids.
A chemical reactor, typically a tubular reactor, could be a packed bed. The pack-
ing inside the bed may have a catalyst to catalyze the chemical reaction. A chemical
reactor may also be a fluidized bed.
Furthermore, catalytic reactors require separate treatment, whether they are
batch, CST, or PF reactors, as the many assumptions of the simpler models are not
valid.
C0 C ΔC 0
(b)
(a) C0 C ΔC large
(d)
C0 C ΔC large C
+cat ΔC large
(c)
C0 C ΔC large
C0
+ cat
(e) (f)
FIGURE 8.1 Types of laboratory catalytic reactors: (a) pulse, (b) differential, (c) integral,
(d) CSTR, (e) recycle, and (f) transport.
174 Chemical Engineering Primer with Computer Applications
Motor
Feed
Cooling jacket
Baffle
Agitator
Mixed product
With a solid phase catalyst and fluid phase reagents, this is proportional to the exposed
area, efficiency of diffusion of reagents in and products out, and efficacy of mixing.
The behavior of the catalyst is also a consideration. Particularly in high-temper-
ature petrochemical processes, catalysts are deactivated by sintering, coking, and
similar processes.
A common example of a catalytic reactor is the catalytic converter following an
engine. However, most petrochemical reactors are catalytic and are responsible for
most of industrial chemical production in the world, with extremely high-volume
examples such as sulfuric acid, ammonia, and reformate/benzene, toluene, ethylben-
zene, and xylene (BTEX).
A catalyst performance depends directly on the following four factors:
1. Activity
2. Selectivity
3. Life and stability
4. Economics and process costs
Reaction Kinetics, Chemical Reactors, and Thermodynamics 175
C0 C0
(a) (b)
C Liquid + catalyst
C0 C
C0 Liquid + catalyst
(c) (d)
C0 Liquid Catalyst C
C + Liquid C0 Catalyst
(e) (f )
C + catalyst
C
Catalyst
C0 + catalyst C0 + catalyst
(g) (h)
FIGURE 8.3 Heterogeneous catalysis reactor types: (a) fixed bed, (b) batch fluid bed,
(c) slurry, (d) catalytic gauze, (e) trickle bed, (f) moving bed, (g) continuous fluid bed, and
(h) transport line.
Products and
unreacted materials
to separation
Catalyst on
support
Diffuser
Reactants
FIGURE 8.4 Fixed-bed reactor (e.g., this is used in the manufacture of ammonia).
Definition of θm: It is the average of time periods during which reaction mixtures
stay in the reactor and is given by the equation
dN A
qm =
ò V(-r ) ,
A
between the limits N A 0 and N Af (8.8)
where
NA is the number of moles of reactant A
V is the reactor volume
–rA is the rate of disappearance of A
Reaction Kinetics, Chemical Reactors, and Thermodynamics 177
dx A
qm = N A 0
ò V(-r ) ,
A
between the limits 0 and x Af (8.9)
Under steady-flow conditions with constant fluid density, θm for a continuous flow
reactor may be calculated by using
V (8.10)
qm =
q
V V V 1
f= = = = (8.11)
(F /) (FA 0 /CA 0 ) q 0 S
φ = CA0 dxA/(–rA)
= V/(F/ℓ) =
V/(FA0/CA0) = V/q0 = 1/S
Third:
find V = [FA0/CA0]
178 Chemical Engineering Primer with Computer Applications
where
F and FA0 are the flow rates in moles per unit time of the total feed and of reactant
A in the feed, respectively
is the molar density in moles per unit volume of the feed
CA0 is the concentration of A in the feed
q0 is the volumetric flow rate of the feed at entering conditions
S is the space velocity
To summarize this, the required capacity of the reactor V is computed from the
specified feed rate and feed concentration of component A and with the knowledge
of space time, ϕ.
Po
Mass A
Gas Gas
T1 T2 Insulation
Gas at
Copper partition
Mg + PoA = PA pressure, P
Over time, T1 T2
where υ is the volume per unit mass, (1/), that is, if we know υ and T is the tem-
perature, we know P is the pressure, etc. Any of these is equivalent to the equation
f(P, u, T) = 0, which is known as an equation of state. The equation of state for an ideal
gas, which is a very good approximation to real gases at conditions that are typically
of interest for aerospace applications, is
Pu = RT,
where
u is the volume per mol of gas
R is the universal gas constant
8.3.2.1 Entropy
The idea of entropy comes from a principle of thermodynamics dealing with energy.
It usually refers to the idea that everything in the universe eventually moves from
order to disorder and entropy is the measurement of that change. The following is a
list of definitions of entropy from a collection of textbooks:
We see evidence that the universe tends toward highest entropy many places in our
lives. A campfire is an example of entropy. The solid wood burns and becomes ash,
smoke, and gases, all of which are more disordered than the solid fuel. For isolated
systems, entropy never decreases. This fact has several important consequences in
science: it prohibits perpetual motion machines, it implies that the arrow of entropy
has the same direction as the arrow of time, and so on. Increases in entropy cor-
responds to irreversible changes in a system, because some energy is expended as
waste heat, limiting the amount of work a system can do.
r1
¾¾
aA + bB ¬ ¾ ® cC
¾ r2
r1 = k1[A]a [B]b
r2 = k 2 [C]c
Now, at equilibrium
r1 = r2
Therefore,
and we get
k1 [C]c [ Products]
= Kp = =
k2 [A]a [B]b [ Reactants]
Given the free energy of formation of the reaction, ∆G, we can write
DG = -RT ln K p
where Kp is the equilibrium constant for the reaction at the given conditions of T and
P, expressed as the concentrations of the participating species in terms of the partial
pressures.
Reaction Kinetics, Chemical Reactors, and Thermodynamics 181
DG = DH - TDS
where ΔH and ΔS are the increase in total heat content and entropy in the reaction,
respectively.
Example 8.1
Calculate the entropy of the surroundings for the following two reactions:
DH = -2045 kJ
DH = +44 kJ
Solution
The change in entropy of the surroundings after a chemical reaction at constant
pressure and temperature can be expressed by the formula
DH
DSsurr = -
T
where
ΔSsurr is the change in entropy of the surroundings
−ΔH is the heat of reaction
T is the absolute temperature in Kelvin
Reaction (a)
DH
DSsurr = -
T
-2045 kJ
DSsurr = -
25 + 273
2045 kJ
DSsurr =
298 K
DSsurr = 6.86 kJ/K or 6860 J/K
Note that the increase in the surrounding entropy since the reaction was
exothermic.
182 Chemical Engineering Primer with Computer Applications
Reaction (b)
DH
DSsurr = -
T
+44 kJ
DSsurr = -
298 K
This reaction needed energy from the surroundings to proceed and reduced the
entropy of the surroundings.
Conclusions: The change in entropy of the surroundings of reactions 1 and 2 was
6860 and −150 J/K, respectively.
Example 8.2
We can calculate ΔrG° for this reaction from the tables of free energies of forma-
tion. We find that ΔrG° for this reaction is very large and of negative value, which
means that the reaction is very feasible. We conclude, right away, that the equilib-
rium constant, K, for the reaction is very large:
DG = -RT lnK
However, we can mix hydrogen gas and oxygen gas together in a container, even
in their correct stoichiometric proportions, and they will stay there for centuries,
perhaps even forever, without reacting. The problem is not that the reactants do
not want to form the products, they do, but they cannot find a pathway to get from
reactants to products. In other words, thermodynamics is not the only main focus
in chemistry. Not only do we have to know whether a reaction is thermodynami-
cally favored, we also have to know whether the reaction can or will proceed at a
finite rate. The study of the rate of reactions is called chemical kinetics.
Example 8.3
CH3OH ® HCHO + H2
Reaction Kinetics, Chemical Reactors, and Thermodynamics 183
The fresh feed to the process was 0.5 kmole/h of O2 and an excess methanol. All
of the O2 reacts in the reactor. Formaldehyde and water are removed from the
product stream first, after which H2 is removed from the recycled methanol. The
recycle flow rate of methanol was 1 kmole/h. The ratio of methanol reacting by
decomposition to that by oxidation was 3. Draw the flow diagram and then cal-
culate the per pass conversion of methanol in the reactor and the fresh feed rate
of methanol.
Calculations:
The flow diagram for the aforementioned process is given here:
Hydrogen
Recycle methanol
1 kmole/h Hydrogen
separator
The numbers of moles taking part in the reactions are given as follows:
From reaction (8.3.1): In 1 h, 0.5 mole of O2 reacts with 1 mole of CH3OH.
Therefore, for complete conversion of oxygen, 0.5 moles/h of O2 needs
1 mole/h of CH3OH.
Obtaining the ratio of methanol conversion, we then find that the flow rate in
methanol decomposition is three times greater than that in methanol oxidation,
that is, 3 moles/h of CH3OH.
Therefore,
Total methanol entering the reactor = methanol reacted by oxidation + metha-
nol reacted by decomposition + recycled methanol = 1 + 3 + 1 = 5 kmoles/h
Conversion per pass = 100 × (total methanol reacted/total methanol entering
the reactor) = 100 × 4/5 = 80%
Fresh methanol rate = total methanol entering the reactor − recycled metha-
nol = 5 − 1 = 4 kmoles/h (Subbu, 2014)
184 Chemical Engineering Primer with Computer Applications
CO + 2H2 Û CH3OH,
PCH3OH y 1
Kp = = CH3OH
pCOpH22 yCOyH22 P2
When the carbon monoxide and hydrogen enter the reactor in a 1: 2 ratio, find the
equilibrium conversion when the pressure is 50 atm and Kp = 0.0016.
Feed Product
Equilibrium
1 mole CO reactor 1–x mole CO
2 moles H2 2–2x mole H2
x mole CH3OH
A conceptual model is formed as follows, where x is the fraction reacted from the
carbon monoxide at equilibrium.
Hence,
x (3 - 2x )
2
1
Kp = 0.0016 =
(1- x ) (1- x )
2
4 ´ 502
Example 8.5
Methane is used to produce acetylene as the main product by the following ther-
mal reaction:
Given
Calculate the product component flow rates and the yield of moles of C2H2/mole
of CH4. Also, calculate the yield of H2 per mole of H2 introduced to the reactor in
the form of CH4.
N2 mole C2H2
N4 mole H2
(Products/s)
Solution
The conversion rate of methane is 60%; hence, N1 = 10(1 − 0.6) = 4 moles CH4
(unreacted).
Carbon balance, in = out
or
2N2 + N3 = 6 (8.5.1)
or
Note that the moles of hydrogen in this balance are meant to be atoms of hydrogen.
Now, we have two equations and three unknowns (N2, N3, N4). We use the
enthalpy balance next to come up with the third relationship, in order to have
three equations in three unknowns.
186 Chemical Engineering Primer with Computer Applications
Q = DH = éê
ë ånH ùûú
i i
products
- éê
ë ånHùûú
i
reactants
Rearranging and substituting for the enthalpies, Hs with the corresponding values,
we obtain
Now, we can enter these three equations, 8.5.1 through 8.5.3, into the MATLAB
matrix form: A = xb. Solving for the value of
Yield of hydrogen = 9.5017 / éë10 moles methane ´ 4 (hydrogen atoms per mole ) ùû
= 0.2375 mole /mole of methane
REFERENCES
Dixon, J., Espenson, J. H., Kinetics and mechanism of oxygen atom abstraction from a
dioxo-rhenium(VII) complex, Inorganic Chemistry 41: 4727, 2002.
Myers, R. The Basics of Chemistry, London, U.K.: Greenwood Publishing Group, 1951.
Subbu, M. Chemical Engineering Learning Resources, 2014. http://www.msubbu.in.
9 Chemical Plant Design
and Process Economics
Safety and
Environmental
Health loss
problems
prevention
Process design:
Design procedure: A–Z chart
1. A written description of
Feasibility study
the process
Process design 2. Make a qualitative BD
Economic evaluation
3. Perform MB and EB
Plant erection
Startup
Production
The successful convergence of theory and practice has a great impact on developing
courses in plant design. This chapter will highlight the different aspects involved in
the design of chemical plants. Since the topic is too involved, we attempt to provide
a balanced presentation to the strategy of process synthesis in plant design con-
sistent with industrial practice and within the economic frame of feasible projects.
The dollar sign must always be kept in mind when carrying out the design of a
chemical plant.
This chapter covers the basic fundamentals in plant design that the reader should
acquire, from the inception of an idea all the way to the construction of a plant.
Process design is a cornerstone in this presentation. General design considerations,
187
188 Chemical Engineering Primer with Computer Applications
9.1 INTRODUCTION
A chemical plant is one that belongs to the chemical industry. Most processes in
the chemical industry involve a chemical change, which could imply two distinctive
groups of industry:
1. One The first group of industry involves a chemical change, for example,
in the case of chemical reaction of nitrogen with hydrogen that produces
ammonia.
2. The second group of industry involves physicochemical change, involved in
the separation of a given feed to components or fractions, for example, the
distillation of hydrocarbons and water desalination. These industries are
normally recognized as physical separation operations that utilize strictly
unit operations.
Purely mechanical changes or separations are not considered part of the
chemical processing.
Plant design for processes of the first group utilizes, in its initial stage, the basic
chemical stoichiometric equations describing the chemical reactions, in order to
do material balance (MB) and heat balance (HB) calculations. However, for indus-
tries of the second group, calculations are usually based on setting up the total
MB and component MB, for example, the solution of binary-distillation problems
involved in the setting up of two equations in two unknowns, as was presented in
Chapter 6.
Input
Heat transfer
Fluid mechanics
Diffusional (mass
transfer)
Kinetics and reactor
design
Instrumentation
and procees control
Output
Thermodynamics PD
A profitable
working plant
Fundamentals or process
of
economics
Design considerations:-
1. Plant location
2. Plant layout
3. Others
In a plant design course, students are basically well prepared to launch a design
project, once they cover the following additional topics: economics, materials selec-
tion, computer-aided design (CAD), and optimization technique, and others. Their
personal information has an impact in handling a design project as well.
Material
balance
No
No
Economic Complete
Procurement Development
Detailed eval. process
based on
eng design preconst design
quotations
cost est.
a collective and comprehensive flowchart that takes the student through the whole
plant to be: from feasibility study to commercial production via the intermediate
steps of process design followed by economic evaluation.
In general, the following three cases of design projects are common to handle:
Full size
Lab Bench Pilot commercial
scale scales plant production
(i) Definition: It is the amount of overdesign to account for changes in operating performance
of an equipment.
(ii) Applications: 1. To safeguard against fouling in: (a) heat exchangers and (b) furnace pipes
2. For potential increase in capacity requirements
(iii) Calculations for Ht Ex:
1. To find out if a given Ht Ex is fouled or not?
RT = Total resistance = Rfi + Rw + Rf0 + Rsi + Rso, define Rsi + Rso = Rst
Steps:
Calculate Uo = Q/A ∆Tm
Find the values of Rfi, Rfo, and Rw
Calculate Rst, using Uo = 1/[Rfi + Rfo + Rw+ Rst]
Compare the value of Rst, calculated, with the value for clean
Ht Ex....If Rst is >> 0.001, then the tube are fouled.
• Health and safety hazards: One should consider the toxicity of materials
and frequency of exposure, fire, and explosion hazards.
• Loss prevention: HAZOP study.
• Environmental protection and control include air pollution, water pollution,
solid wastes, thermal effects, and noise effects.
• For plant location: Both primary and specific factors are to be considered.
• For plant layout: Optimum arrangement of equipment within a given area
is a strategic factor.
• For plant operation and control: The designer should be aware of
• Instrumentation
• Maintenance
• Utilities
• Structural design
• Storage
• Materials handling, pipes, and pumps
• Patents aspects
Chemical Plant Design and Process Economics 193
Obviously, the estimation of the fixed capital investment (FCI) would enable the
designer to answer the following questions:
(1) (2)
Cost Profitability
estimation analysis
T. capital T. product
investment cost (TPC)
Direct Indirect†
1. Land
2. Process equipment
3. Buildings
4. Storage facilities
5. Utilities
6. Auxiliaries: workshops, laboratories, others
Apply
Calculate economic
income indicators
FIGURE 9.6 Economic aspects in plant design. *Working capital is the money needed for raw
materials and cash at start up and is equal to 10%–20% of T.C.I. †Indirect fixed capital is spent
in engineering expenses, supervision for construction and construction tools and equipment.
194 Chemical Engineering Primer with Computer Applications
Lang factors are recommended in this regard. Once the costs of delivered equipment
are estimated, factors are applied as shown next:
Factor to Estimate
Classification of the Plant Function FCI TCI
Solid processing operation 3.9 4.6
Solid/fluid processing operation 4.1 4.9
Fluid processing operation 4.8 5.7
The six-tenth factor is another handy method to determine the cost of an equipment
of a given capacity (A), if the cost and the capacity of a similar one (B) are known:
é i ù
P = A d (SFF ) and A d = P ê ú
ë (1 + i) - 1 û
n
2. If “A” deposited annually, on the other hand, is the amount of the annual
capital recovery of the equipment, A r, then the amount of money generated
by the sinking fund will be the future value for the cost of this equipment:
F = Ar(SFF)
To make it clearer, the Ad, defined as the annual depreciation cost, when deposited in the
sinking fund, will generate the P of the equipment. This is simply to recover the original
cost paid in the equipment. On the other hand, when Ar(>Ad), defined as the annual
capital recovery, is deposited in the sinking fund, then the amount generated will be F.
In this case, F will be equal to P + the interest on P, as if it is deposited in a bank
for n year or simply
F = P(1 + i)n
Bank Industry
Fb = P(1+i)n (1+i)n – 1
Fp = A
i
Let:
Fb = Fp
(1+i)n i
Ar = P
(1+i)n – 1 Sinking fund F (future value)
i
Ad = P Sinking fund P (present worth)
(1+i)n – 1
Annual profit
ROI = 100
Capital investment
For projects, where the cash flow extends over a number of years, the average rate of
return is calculated using an average value for the profit, by dividing the sum of the
annual profits by the useful lifetime:
é
åAnnual profits/n ùú
n
ê 1= y
ROI = ê ú 100
ê (capital investment) ú
ë û
196 Chemical Engineering Primer with Computer Applications
9.5.4.3 Discounted Cash Flow Rate of Return and Present Value Index
If we have an oil asset (oil well, surface treatment facilities, a refining unit, etc.) with
an initial capital investment “P,” generating annual cash flow over a lifetime n, then
the DCFR is defined as the rate of return or interest rate that can be applied to yearly
cash flow so that the sum of their present value equals P.
From the computational point of view, DCFR cannot be expressed by an equation
or formula, similar to the previous methods. A three-step procedure involving trial
and error is required to solve such problems, as illustrated in the solved example in
Example 9.3.
Example 9.1
In a gas plant, a stream of natural gas is processed in order to recover the valu-
able hydrocarbons C2 and C3/C4 and to sell the bulk of methane (C1) as a fuel gas,
sales gas.
Feed: n.g.
Gas plant C2………… E
F
C3/C4……. P
Chemical Plant Design and Process Economics 197
The feed rate is 2.5 × 106 g ⋅ mole/h. The feed gas composition is as follows:
Component Mole%
C1 83.0
C2 12.0
C3 4.0
n‐C4 1.0
Solution
Mole balance will give the following
Given
y1 + y 2 = 1.0 (9.1.5)
We have six equations and seven unknowns; F is given. Since the value of C2 in
the sales gas, y1, is specified to be 0.0, the solution is obtained readily as follows:
Answer
S = 2.07 ´ 106
E = 0.31´ 106
P = 0.12 ´ 106
Example 9.2
mass flow rate (kg/h) and the mass fractions of each species in the inlet (stream 1),
F, and each outlet are known.
F1 = ?
Solution
A = 4%
If we define the unknowns as
B = 93%
x1 = F1, x2 = F2, x3 = F3
and set up the mass balances for C = 3%
F = 10 F2 = ?
1. The total mass flow rate
x1 + x2 + x3 = 10 A = 20% A = 54%
B = 60% B = 24%
2. The mass balance on species 1
C = 20% C = 22%
0.04x1 + 0.54x2 + 0.26x3 = 0.2*10
3. The mass balance on species 2
0.93x1 + 0.24x2 = 0.6*10
These three equations can be written in matrix form
1 1 1 x1 10
F3 = ?
0.04 0.54 0.26 x2 = 2
A = 26%
0.93 0.24 0 x3 6 B = 0%
C = 74%
To find the values of unknown flowrates write the code:
A = [1,1,1;.04,.54,.26;.93,.240
, ]; B = [10;.2*10;.6*10];
F2 = 2.4330
F3 = 1.743
Example 9.3
The capital cost of a heat exchanger, P, is $10,000 with a lifetime of 10 years. If the
money can be invested at 6% (annual interest rate), calculate
Solution
(c) The difference between the two values is = 1359–759 = $600; the
interest on the capital, if invested in a sinking fund, gives 600(13.18) =
$7911
The Ar, on the other hand, if invested in sinking fund, should give
Example 9.4
Assume an oil company is offered a lease of oil wells that would require a total
capital investment of $110,000 for equipment used for production. This capital
includes $10,000 working money, $90,000 depreciable investment, and $10,000
salvage value for a lifetime of 5 years.
Cash flow to project (after taxes) gained by selling the oil is as given next.
Based on calculating the DCFR, a decision has to be made: should this project be
accepted?
Solution
Our target is to set the following equity: by the end of 5 years, the future worth of
the cash flow recovered from oil sales should break even with the future worth of
the capital investment, had it been deposited for compound interest in a bank at
an interest rate i. This amounts to say that
Fo = Fb
where
Fb = 110,000(1 + i)5 for banking
Fo = Si5=1Fi for oil investment that represents the cash flow to the project, com-
pounded on the basis of end-of-year income; hence
Notice that the $20,000 represents the sum of WC and salvage value; both are
released by the end of the fifth year. Setting up FB = Fo, we have one equation
involving i as the only unknown, which could be calculated by trial and error.
The value of i is found to be 0.207, that is, the DCFR = 20.7%.
Example 9.5
Calculate the NPV of the cash flow for the oil lease described in the example
earlier, if the money is worth 15%.
Solution
At i = 0.15, the annual cash flow is discounted. The present value of the sum of
the cash flows = $127,000.
That is, the oil lease can generate $17,000 (evaluated at today’s dollar value) more
than and above the totally recovered capital investment. Solution is illustrated in
Table 9.1:
TABLE 9.1
DCFR: Results for Solved Example 9.4
i = 5% i = 20% i = 25% i = 20.7%
Year Cash Present Present Present Present
(y) Flow dy Value ($) dy Value ($) dy Value ($) dy Value ($)
0 110,000
1 30,000 0.8696 26,088 0.8333 24,999 0.8000 24,000 0.8290 24,870
2 31,000 0.7561 23,439 0.6944 21,526 0.6400 19,840 0.6870 21,297
3 36,000 0.6575 23,670 0.5787 20,833 0.5120 18,432 0.5700 20,520
4 40,000 0.5718 22,872 0.4823 19,292 0.4096 16,384 0.4720 18,880
5 43,000 0.4971 21,375 0.4019 17,282 0.3277 14,091 0.3910 16,813
20,000
Total 117,444 103,932 92,747 102,380
P.V.I 1.07 0.94 0.84 0.93
N.P.V $7,444.40 −$6,067.70 −$17,252.90 −$7,620.00
10 Case Studies
10.1 INTRODUCTION
This chapter is best introduced by recalling the definition that chemical engineering
is the field that applies the core scientific disciplines of chemistry, physics, biology,
and mathematics to transform raw materials or chemicals into more useful or valu-
able forms, invariably in processes that involve chemical change.
Case Studies in Chemical Engineering is one of a set of resources published
by the chemical company ECI between 1971 and 1982. The booklets were the
result of a series of workshops at which contributors from industry and schools
collaborated.
The booklet is perhaps only of historical interest, but it does show how the
principles of chemical engineering were applied to improve the existing processes
and devise new ones in the days before the principles of green chemistry had been
enumerated.
CACHE (Computer Aids for Chemical Engineering; CACHE Corporation, 2016)
Process Design Case Studies represents an outstanding resource of many case
studies. These are covered in 11 different volumes that span “Separation System for
Recovery of Ethylene and Light Products from Naphtha Pyrolysis Gas Steam” in
Volume 1, all the way to “Conceptual Design of an Aromatics Plant from Shale Gas”
in Volume 2.
In Chapter 10, we attempt to introduce a number of case studies that touch on
practical technical problems, some of which the author has personally experienced
in his research work, both in academic institutions and in industry. Many of the
principles and fundamentals presented in the preceding chapters of our book are
introduced and illustrated as given in these case studies. Case studies 5 and 6
are taken from the open literature. The list of case studies presented includes the
following:
201
202 Chemical Engineering Primer with Computer Applications
10.2.3 aPProacH
The possible alternatives are
Step 1: Since ammonia is a vapor under normal conditions, one might store it
in large gasholders. This idea is ruled out because of the excessive volume
of ammonia at standard T and P, as shown next.
Reject brines
Soda ash and
CO2 magnesium
NH3 chloride products
Chemical reactions to
convert saline water Ammonium
NH3
into valuable products chloride and partially
synthesis
desalted water
H2 Raw mat. Ammonia-based
Solar energy
Water fertilizers
electrolysis Source
of energy
Desert development
Step 2: Let us calculate the vapor pressure of liquid ammonia at low temperatures.
Let us assume that ammonia can be stored as a liquid at its normal boiling
point. Then, we have to take into consideration the heat losses from ammo-
nia tank to the atmospheric conditions.
Step 3: The primary resistance to heat transfer would be offered by the tank–
air interface, with h equal to about 10 Btu/h · ft2 · °F and k, of an insulation,
is assumed to be equal to 0.15 Btu/h · ft · °F.
The surface area of the tank exposed to the atmosphere is
A = (5p / 4)(4V/p)2 / 3
where
Therefore,
A = 8.5 × 102 ft2
The heat transfer rate is estimated as a function of the insulation thickness
X (inches):
10.2.4 conclusions
It is observed that the heat transfer Q to the ammonia storage tank would lead to
excessive loss of ammonia. Even with a foot of insulation, about 105 Btu/h will be
transferred into the tank. This will result in the boiling of about 200 lb of ammonia
per h (1 Btu vaporizes 590 lb ammonia).
It is therefore suggested to use 6 in. insulation that reduces the heat loss to
Safety precautions have to be applied and start-up purge plan has to be considered
to eliminate hazards of explosion involved during filling. In addition, the storage
capacity of 60,000 tons could be divided into two units: one for loading and the other
for unloading.
10.2.5 discussions
1. Consider other possibility to store ammonia in a solid compound, from
which ammonia is released as demanded; is it a feasible proposal?
Could large underground cavern be considered as another option?
2. An optimization model could be carried out to calculate the optimum insu-
lation thickness for the storage tank. This could lead to a substantial savings
when choosing a small insulation thickness.
6
1 2 Sludge
Taking the treatment 7
wastewater Screening 5
Back to
away 4 Final the river
3
Secondary treatment
Primary
treatment treatment
FIGURE 10.2 Sewage treatment process. (After Thames Water, Facts and Figures, http://www.thameswater.co.uk/about-us/4625.htm, retrieved on
May 31, 2014.)
205
206 Chemical Engineering Primer with Computer Applications
Water
Dewatered sludge
Heat input
Power Incineration
generation Resd
Exhaust steam
Electricity
FIGURE 10.3 Proposed scheme for treating sewage sludge to produce water and electricity.
A
B
F H
C
G
D
E
heat exchanger. The steam cools and condenses against the heat exchanger tubes,
thereby heating the feedwater as described earlier. By that time, the sewer sludge
has been subjected to boiling, a process that removes all liquid, to be captured as
water vapor upon flashing inside the chambers.
The water is collected through troughs located beneath the heat exchangers (A).
Further treatment and ultrarefining steps are carried out for the produced water,
making it suitable for drinking.
Dewatered sludge is directed to an incinerator, which produces high-pressure
steam that fuels a generator. The generator creates electricity that is used as a source
of power. There is even a little extra electricity left that can be transferred into the
power grid. Exhaust steam, on the other hand, is used for heating purpose and for
steam-jet ejector to create a vacuum.
Based on daily basis, the results of S200 Omni Processor are reported as follows:
1. Plant input
a. A quantity of 10–12 tons of dry fuel
b. It could be
i. Any kind of biomass from sewage sludge
ii. Animal by-products
iii. Agricultural by-products
iv. Garbage (glass- and metal-free)
The fuel does not have to be dry as such. It can contain an appre-
ciable amount of water: 12 tons of fuel can hold up to 70 tons of
water.
Arrangement can be made to have a feed input that contains a mix-
ture of
A. Sewage feed, having more than 90% water
B. Sludge or other solid material, having 30%–90% water
2. Plant output
Three items are produced:
a. The power output is 300 kW.
b. The water output has a total volume of 70,000 L (water produced will
be clear and sterile, but with residual odor; ultrapurification equipment
is required to make it odor-free).
c. About 10%–20% of the dry material input will come out as ash.
10.3.5 conclusions
The proposed scheme presented in our chapter consists of two stages. However,
the number of flashing units (chambers) could multiple. This number is a function
of the water content and the type of sludge. The advantage of applying vacuum in
the first unit is obvious to reduce the temperature in the next units causing flashing
of the entering feed. In addition, reduction of vacuum inside the chambers leads
to lower operating pressure, and hence lower operating temperatures as seen in
Figure 10.5.
208 Chemical Engineering Primer with Computer Applications
10.3.6 discussions
1. Derive the flash evaporation equation that applies for a single component
using a heat balance to predict how much is vaporized from the feed to the
evaporator. In this case, it will be assumed that the feed contains the solid
material (A) suspended in a bulk of water, similar to seawater with soluble
sodium chloride.
2. Elaborate on the incineration process. What are the basic design parameters
to consider in building an incineration furnace?
Solar energy
Photovoltaic array
Sea water
Brine
Desalination Preferential salt separator
Electrolysis cells
Water Magnesium
electrolysis Magnesium
removed by
chloride
Anode Chlorine gas
vacuum
Magnesium
metal
H2
Electrolyte
Magnesium
10.4.4 discussions
1. Draw a block diagram for the aforementioned process, indicating the type
of operation in each step: physical (unit operation) or chemical conversion
(unit process).
2. Identify the chemical reactions that take place.
3. Compare between the proposed method of extracting Mg chloride from
seawater (PSS) and the one in current practice by Dow Chemical Company.
10.5.2 metHodology
In its general term, desalting refers to a water treatment process that removes salt
from saline water resources. It is simply represented as follows:
Case Studies 211
Conventional
Source of energy desalting Brine
process
Fresh water
The proposed chemical desalting process, on the other hand, when compared to the
conventional desalting process, differs on the following grounds:
NH3 CO2
Chemical
Brine (sat.) desalting Soda ash
Fertile water
(contains NH4Cl)
It could be safely stated that the chemical desalting process complements seawater
desalination.
• Primary reaction: The reaction between CO2 and NH3 is described as follows:
Fresh water
Seawater
M S F pant
Brine (30%)
NH3 Mg(OH)2
Absorber Double MgCl2
NH4Cl decomposition
NH4OH
Separation
CO2 Formation of
chemical NH3
products regeneration
Soda ash
(Na2CO3)
NH3 recycle
Reaction 3:
2NH 4OH + 2CO2 + 2NaCl ® 2NaHCO3 + 2NH 4Cl (use 2 moles of NH 4OH)
NH 3 ( recycle )
TABLE 10.1
Consumption–Production Analysis for the Chemical Desalting Process
Reactants Products
Reactions NH3 H2O CO2 MgCl2 NaCl NH4OH Mg(OH)2 NaHCO3 Na2CO3 NH4Cl CO2 W
No. 1 −4 −4 +4
No. 2 −1 −2 +1 +2
No. 3 −2 −2 −2 +2 +2
No. 4 −2 +2 +1
Net −4 −4 −1 −1 −2 0 +1 0 +2 +4 *
Chemical Engineering Primer with Computer Applications
Case Studies 215
10.5.4 conclusion
The proposed process offers a scheme that provides three products as compared to
one product only by the Solvay process. Magnesium chloride is an important product
for the manufacture of magnesium metal. As a matter of fact, one can claim that the
proposed process could compete with the well-known Dow process for the extrac-
tion of magnesium chloride from seawater.
The process synthesis suggested in this paper focuses on the use of ammonia
as a recycle reagent. Regeneration of ammonia is accomplished in the absence of
Ca(OH)2 used in the Solvay process as indicated in the earlier comparison.
The production of fertile water (partially–desalted water) containing NH4Cl
could be a feasible choice; NH4Cl could be sold as a solution to be used in fertilizers
for rice crops.
Total balance F = D + B
F Concentrator/
evaporator
B
216 Chemical Engineering Primer with Computer Applications
The main objective in this case study is to figure out how much chemical
products we produce for a given feed input of brine. Also, it is important
to calculate how much ammonia and carbon dioxide we consume in the
process.
3. Comment on the proposed process as compared to the traditional water
desalination.
Stream A Stream B
1.5 L/min 2 L/min
25 g/L 0 g/L
Stream C
3.5 L/min
x1 g/L
environment. Here, it is assumed that as the two streams, A and B, enter into the
tank, instantaneously the chloride concentration in the tank changes to the exit
concentration, x1.
A material balance on sodium chloride is carried out for the tank system as
follows:
dx1
= (25 g/L)(1.5 L/min) + (0 g/L)(2 L/min ) - (x1 g/L)(3.5 L/min) (10.6.2)
dt
dx1
+3.5x1 = 37.5 (10.6.3)
dt
1 dx1 dx
+ x1 = 0.286 1 + x1 = 10.7 (10.6.4)
3.5 dt dt
X1 = 10.7 + C e - t/0.286
For the initial boundary conditions, at t = 0, the salt concentration in the tank was
given as 750 g/15 L (50 g/L) or at t = 0, x1(0) = 50 g/L.
Therefore,
Checking:
As t increases, e−t → 0, and X1 → a, value = 10.7
And at t = 0, e−t = 1, and x1 → a, value = 50
10.6.4 discussions
• Using Equations 10.6.3 and 10.6.4, numerically determine how the concen-
tration of the discharged salt changes with time.
• Plot a graph for the relationship that represents your solution.
• How long did it take to achieve the minimum required local specifications?
At steady state, what is the concentration of salt being discharged from this
local soap factory? (Kaw and Kalu, 2011)
10.7.2 oBjective
In comparing the traditional versus the reactive distillation methods, such as in the
production of methyl acetate, for example (Figure 10.9), the latter methods have
the following advantages:
• Improve selectivity
• Reduce by-products
• Prevent pollution
• Reduce energy use
• Handle difficult separations
• Avoid separating reactants
Case Studies 219
Condenser
Product 1
Rectification
Feed 1 section
Catalyst
Reaction zone
(extractive distillation zone)
Feed 2
Stripping
section
Evaporator
Product 2
• Eliminate/reduce solvents
• Enhance overall rates
• Beat low equilibrium constants
10.7.3 discussion
The introduction of a separation process side by side in the reaction zone, or vice
versa, could lead to many difficulties in reactive distillation process. Give examples
of some of these problems. On the other hand, some of the merits that apply in the
case of reactive distillation could be applied to other aspects in many unit operations.
Give examples of similar cases, especially in surface oil operations that involve crude
oil and gas field treatment and processing (see texts by the author in bibliography).
REFERENCES
CACHE Corporation, CACHE Process Design Case Studies, The CACHE Corporation, 2016.
http://cache.org/super-store/cache-process-design-case-studies.
The Essential Chemical Industry (ECI), Case Studies in Chemical Engineering, Imperial
Chemical Industries (ICI), London, U.K., 1970–1979.
Kaw, A. and Kalu, E.E., Numerical Methods with Applications: Abridged, 2nd edn.
Tampa, FL: University of South Florida, March 2011.
Appendix A: Application
Problems
A.1 SECTION I (CHAPTERS 1 AND 2)
This set of problems is aimed to go hand in hand with the materials covered in
Chapters 1 and 2. The solution of some simple problems could be done either manu-
ally or using MATLAB® or Excel. It involves basic conversion of units and other
fundamental calculations.
where T is in °C
It is required to have the heat capacity expression in terms of °R.
2. A tank contains CO2 is used to supply a system with CO2. The pressure
gauge reads 50 psi. Determine the absolute pressure in the tank, in psi, if
the barometer reads 28.0 in. · Hg.
3. Calculate the gauge pressure in lbf/in.2 at a depth of 4.0 miles below the sea
surface. Assume the water temperature is 60°F and the specific gravity is
1.042.
4. Calculate the diameter of a pipeline handling 12,000 bbl of oil/h assuming
the velocity of flow in the pipe line is 5 ft/s (1 ft3 = 0.413 barrel).
5. Polonium has a half-life of 140 days. Starting with 50 g today, how much is
left after 200 days?
6. Find the volume of a tubular reactor, with a residence time of the reaction
being 60 s. The input feed flow rate to the reactor is 400 ft3/min.
7. What is the weight of air in a room having a volume of 2000 ft3 at T = 80°F
and P = 750 mm Hg.
8. In one of your experiments in unit operations laboratory, a reading of
400 cm3/min was measured. How do you covert it to in.3/h?
9. In experiment in problem 8, the pressure gauge reads 12 psi and the
atmospheric pressure is 14 psi. What is the absolute pressure in this
experiment?
221
222 Appendix A: Application Problems
10. What will be the concentration of O2 dissolved in water at 298 K when the
solution is in equilibrium with air at 1 atm total pressure?
Hint: Apply Henry’s law, pA = HxA,
where
H is the Henry’s constant = 4.38 × 104 atm/mol fraction
pA is the partial pressure of O2 in air, which is taken as 0.21 atm.
11. Find the density of methane gas at T = 70°F and P = 760 mm Hg, by using
the gas equation.
12. An orifice meter is used to measure flow rates in a pipeline. Given the fol-
lowing relationship,
u = C [ DP/ ]0.5,
where
u is the fluid velocity
ΔP is the pressure drop
ℓ is the density
C is a constant
16. A chemical plant located at the base of a mountain (about 300 ft) has to
dispose of its CO2 by-product. An engineer proposed a unique method for
absorbing the CO2 in water; water is then lifted up the mountain. The CO2
was then vented to the atmosphere and water is recirculated as shown in
the diagram. Assume that CO2 is to be absorbed in water at a CO2 partial
pressure, pCO2 , of 10 psi above that required to lift the water up to the top of
the mountain (300 ft). Calculate the amount of water required to dispose of
5000 ft3 of CO2 at STP.
Hint: (a) Calculate pCO2 to take care of the pressure head for uplifting.
(b) Apply Henry’s law to find the X CO2 absorbed in water: p = Hx.
CO2 venting
Top of the mountain
Water
CO2 absorbed
in water
Chemical
plant
V, yi
Flash
Distill
F, Zi
L, xi
18. A stream of air is passing through a tube. Its flow rate, F, is unknown. In
order to determine the value of F, another gas, CO2, is introduced in the
tube at the entrance point, as shown in the diagram at a rate = 1 ft3/s. The
composition of the effluent stream (air + CO2) that exits the tube was mea-
sured and found to have 0.5% CO2 by volume. Calculate F, using compo-
nent MB on CO2.
CO2
Air + CO2
Tube
Air, feed
QC
D, hD, XD
F, hf
QR
B, hB, XB
22. Applying the Newton’s law of cooling, dT/dt = −k(T−30), where the tem-
perature of moving air is 30°C, T is the temperature of the substance at time
t minutes, and k is a constant, find the time necessary to cool a substance
from 100°C to 40°C, if the same substance cools from 100°C to 70°C in
15 min time.
Appendix A: Application Problems 225
25. A pilot plant distillation column is used to strip alcohol from a feed stream
of ethyl alcohol and water, containing 1.5 mole% alcohol. The following
conditions apply:
where the distillate contains 87% alcohol and the bottom product contains
0.05 % alcohol. Calculate all stream flow rates.
26. At 1000°C and a total pressure of 30 atm, the equilibrium in the reaction
is such that 17 mole% of the gas is CO2. What percent will be CO2 if the
total pressure is 20 atm instead of 30 atm?
27. Draw a computer block diagram for solving flash separation problem by
trial and error.
28. Find the value of V/F for a feed consisting of 10, 20, 30, and 40 mole% of
propane, butane, pentane, and hexane, respectively, to be flashed at 200°F
and 100 psia. The corresponding values of K i for the hydrocarbons are 4.5,
1.8, 0.75, and 0.37 (in the same order).
29. Calculate the temperature and composition of a liquid in equilibrium with
a gas mixture containing 10.0 mole% benzene, 10.0 mole% toluene, and the
balance nitrogen (which may be considered noncondensable) at 1 atm. Is the
calculated temperature a bubble-point or dew-point temperature?
30. Under limited conditions, volume balance can be applied to liquids in con-
junction with a mass balance. Non–polar mutually–soluble liquids in gen-
eral are blended so that the volume is approximately equal to the sum of the
226 Appendix A: Application Problems
liquid volumes of the components. Find the yields of light and heavy cuts
that can be produced from 100 gal of gas oil under the given conditions:
L, ρ = 0.8
H, ρ = 0.93
A ® B.
Conductivity
K = 0.34 + 9.278 * 10 -4 T
34. Make a table of (x vs. y) values for the three components such as ethanol,
water, and benzene, at T = 373.15 K. The vapor pressure of these three
components is calculated by (Note: Take 11 points for each component from
0 to 1)
æ 18.5242 - 3578.91 ö
Ethanol Peo = exp ç ÷
è T - 50.5 ø
Appendix A: Application Problems 227
æ 18.3036 - 3816.44 ö
Water Pwo = exp ç ÷
è T - 46.13 ø
æ 15.9008 - 2788.51 ö
Benzene Pbo = exp ç ÷
è T - 52.36 ø
where
K i = Pio /Pt
Pt = 760
yi = K i ´ xi
N 2 + 3H 2 ® 2NH 3
The total feed (T moles) is made up of fresh feed (F moles) + recycle stream
(R moles). The gross product (G moles) is 12% NH3 moles, which goes to a
separator, removing all NH3 and allowing recycled stream to return to the
convertor. For F = 100, calculate G and R.
F
Feed
20% C2 A
40% C3
Xp = ?
40% C4
B
XB = ?
Propane fraction
8.4% C3
91.6% C4
228 Appendix A: Application Problems
37. Equilibrium constants (k values) are needed for our distillation calcula-
tions, as was explained earlier. Write a program to calculate the vapor
pressure for the following system and to calculate and tabulate its k
values.
æ 14.2724 - 2945.47 ö
P1 = exp ç ÷
è T + 224 ø
æ 14.2043 - 2972.64 ö
P2 = exp ç ÷
è T + 2 09 ø
Pt = 101.3
Pi
Ki =
Pt
T1 = 1000°F q
3 in.
T2 = 200°F
Mole fractions of hexane in liquid x 1.00 0.69 0.40 0.192 0.045 0.00
Mole fractions of hexane in vapor y 1.00 0.932 0.78 0.538 0.1775 0.00
where
x is the mole fraction of methanol in liquid
y is the mole fraction of methanol in vapor
Insulation D1 D2 D3
Pipe
1
U=
æD ö æ D3 ö
D3 × ln ç 2 ÷ D3 × ln ç D ÷
D è D1 ø+ è 2ø + 1
+
D1 × h in 2 × k PIPE 2 × k INSULATTION h AIR
where
k PIPE and k INSULATION are the thermal conductivities of the pipe and
insulation
hin is the heat transfer coefficient for fluid flowing in pipe
hAIR is the heat transfer coefficient due to air flowing outside the pipe
Appendix A: Application Problems 231
The first two terms of the denominator in this equation are generally smaller
compared to the remaining terms and can be neglected. For this example,
the first term due to pipe fluid is ignored.
47. An evaporator is fed with 10,000 kg/h of a solution containing 1% solute
by weight. It is to be concentrated to 1.5% solute by weight. The feed is
at a temperature of 37°C. The water is evaporated by heating with steam
available at a pressure of 1.34 atm absolute, corresponding to a temperature
of 108.3°C. The operating pressure in the vapor space is 1 atm absolute.
Boiling point elevation and other effects can be neglected. The conden-
sate leaves at the condensing temperature. All the physical properties of the
solution may be taken to be the same as that of water. What is the quantity
of steam required per hour?
Data
Enthalpy of feed = 38.1 kcal/kg
Enthalpy of solution inside the evaporator (at 100°C) = 98 kcal/kg
Enthalpy of vapor at 100°C = 644 kcal/kg
Latent heat of vaporization of steam = 540 kcal/kg
48. Methanol is produced by the following reaction:
CO + 2H 2 ® CH 3OH
F2
Reactor Separator
F1 33% CO, 66% H2, 1% impurities F3
æ 14.2724 - 2945.47 ö
P1 = exp ç ÷
è T + 224 ø
æ 14.2043 - 2972.64 ö
P2 = exp ç ÷
è T + 209 ø
Pt = 101.3
Pi
Ki =
Pt
235
236 Appendix B: Mathematical Methods and Review Notes
B.2.1 introduction
Two approaches can be pursued to seek a solution: analytical and numerical.
An analytical solution to an equation or system is a solution that can be arrived at
exactly using some mathematical tools. For example, consider the function y = ln (x),
graphed in the following.
1 2 3 4 5
–1
–2
–3
The root of this function is, by convention, when y = 0 or when this function
crosses the x-axis. Hence, the root will occur when ln (x) = 0 ≥ x = e0 = 1.
The answer x = 1 is an analytical solution because through the use of algebra,
we were able to come up with an exact answer. On the other hand, attempting to
solve other equations could be hard to achieve. In such cases, it is necessary to
seek a numerical solution, in which guesses are made until the answer is close
enough.
Bisection
Bracketing
methods Regula falsi
Interpolation
Newton’s
Open
Secant
methods
Interpolation
Combination Brent’s
1. Bisection method
2. Newton–Raphson method (Tangent method)
3. Secant method
Proof:
Use theory
x = cos x
All roots of x – cos x = 0 lie in the interval [–1,1] |x|=|cos x|≤1
Difficult cases:
Bisection method
a+b
Given a bracketing interval [a,b], compute x = and sign( f (x));
2
repeat using [a,x] or [x,b] as new bracketing interval.
function x=bisection(f,a,b,tol)
sfb = sign(f(b)); >> f = @(x) x–cos(x);
width = b–a; >> bisection(f,0.7,0.8,1e–3);
disp(‘ a b sfx’) a b sfx
while width > tol 0.70000000 0.80000000 1
width = width/2; 0.70000000 0.75000000 –1
x = a + width; 0.72500000 0.75000000 –1
sfx = sign(f(x)); 0.73750000 0.75000000 1
disp(sprintf(‘%0.8f %0.8f %2.0f’, [a b sfx])) 0.73750000 0.74375000 1
if sfx == 0, a = x; b = x ; return 0.73750000 0.74062500 –1
elseif sfx == sfb, b = x; 0.73906250 0.74062500 1
else, a = x; end
end
• Guaranteed convergence
• Predictable convergence rate
• Rigorous error bound
240 Appendix B: Mathematical Methods and Review Notes
Disadvantages
f(xk)
xk+1 = xk –
f΄(xk)
x2 x1 x0
function x=newtonraphson(f,df,x,nk)
disp(‘k x_k f(x_k) f’’(x_k) dx’)
for k = 0 :nk
dx = df(x)\f(x);
disp(sprintf(‘%d %0.12f %9.2e %1.5f %15.12f’,[k,x,f(x),df(x),dx]))
x = x – dx;
end
Iteration formula
f (xk )
xk+1 = xk –
f (xk ) – f (xk –1)
xk – xk –1 x2 x1 x0
function xx = secant(f,xx,nk)
disp(‘k x_k f(x_k)’)
ff = [f(xx(1)), f(xx(2))];
>> f = @(x) x–cos(x);
h = 10*sqrt(eps); >> secant(f,[0.7 0.8] , 6);
for k = 0 : nk
disp(sprintf(‘%d %17.14f %14.5e’,... k x_k f(x_k)
[k,xx(1),ff(1)])) 0 0.70000000000000 –6.48422e–02
if abs(diff(xx)) > h 1 0.80000000000000 1.03293e–01
df = diff(ff)/diff(xx); 2 0.73856544025090 –8.69665e–04
else 3 0.73907836214467 –1.13321e–05
df = (f(xx(2)+h)–ff(2))/h; 4 0.73908513399236 1.30073e–09
end 5 0.73908513321516 –1.77636e–15
xx = [xx(2), xx(2)–ff(2)/df]; % update xx 6 0.73908513321516 0.00000e+00
ff = [ff(2), f(xx(2))]; % update ff
end
Appendix B: Mathematical Methods and Review Notes 241
Formulation step
may result in:
REFERENCE
Jensen, W.B. The universal gas constant R, Journal of Chemical Education 80(7): 731,
July 2003.
Weisstein, E.W., Constant, Wolfram MathWorld. http://mathworld.wolfram.com/Constant.
html. Retrieved April 13, 2011.
Appendix C:
Conversion Tables
Common Engineering Design Conversion Factors
Given Equivalent Value Unit of Measure
Length [L]
Foot (ft) 0.304800 Meter (m)
Inch (in.) 25.4000 Millimeter (mm)
Mile (mi) 1.609344 Kilometer (km)
Area [L]2
ft2 0.092903 m2
in.2 645.16 mm2
in.2 6.45160 cm2
Volume [L]3 and capacity
in.3 16.3871 cm3
ft3 0.028317 m3
ft3 7.4805 Gallon
ft3 28.3168 Liter (L)
Gallon 3.785412 Liter
Energy, work, or heat [M] [L]2 [t]−2
Btu 1.05435 kJ
Btu 0.251996 kcal
Calories (cal) 4.184* Joules (J)
ft · lbf 1.355818 J
ft · lbf 0.138255 kgf · m
hp · h 2.6845 MJ
kWh 3.600 MJ
m · kgf 9.80665* J
N·m 1 J
Flow rate [L]3 [t]−1
ft3/min 7.4805 gal/min
ft3/min 0.471934 l/s
gal/min 0.063090 l/s
Force or weight [M] [L] [t]−2
kgf 9.80665* Newton (N)
lbf 4.44822 N
lbf 0.453592 kgf
(Continued)
243
244 Appendix C: Conversion Tables
Fracture toughness
ksi · sqrt(in.) 1.098800 MPa sqr(m)
Heat content
Btu/lbm 0.555556 cal/g
Btu/lbm 2.324444 J/g
Btu/ft3 0.037234 MJ/m3
Heat flux
Btu/h · ft2 7.5346E−5 cal/s · cm2
Btu/h · ft2 3.1525 W/m2
cal/s · cm2 4.184* W/cm2
Mass density [M] [L]−3
lbm/in.3 27.68 g/cm3
lbm/ft3 16.0184 kg/m3
Power [M] [L]2 [t]−3
Btu/h 0.292875 Watt (W)
ft · lbf/s 1.355818 W
Horsepower (hp) 745.6999 W
Horsepower 550* ft · lbf/s
Pressure (fluid) [M] [L]−1 [t]−2
Atmosphere (atm) 14.696 lbf/in.2
atm 1.01325E5* Pascal (Pa)
lbf/ft2 47.88026 Pa
lbf/in.2 27.6807 in. H2O at 39.2°F
Stress [M] [L]−1 [t] −2
kgf/cm2 9.80665E−2* MPa
ksi 6.89476 MPa
N/mm2 1 MPa
kgf/mm2 1.42231 ksi
Specific heat
Btu/lbm · °F 1 cal/g · °C
Temperature
Fahrenheit (°F − 32)/1.8 Celsius
Fahrenheit °F + 459.67 Rankine
Celsius °C + 273.16 Kelvin
Rankine R/1.8 Kelvin
Thermal conductivity
Btu · ft/h · ft2 · °F 14.8816 cal · cm/h · cm2 · °C
Appendix C: Conversion Tables 245
1 Quad = 1015 Btu (world energy usage is about 300 Quads/year, and United
States is about 100 Quads/year in 1996).
1 therm = 100, 000 Btu .
1000 kWh = 3.41 million Btu .
247
248 Bibliography
Griskey, R.G., Transport Phenomena and Unit Operations: A Combined Approach (paper-
back), New York: John Wiley, 2006.
Harj, B., Cambridge Reactor Design Limited, Cambridge, U.K.: Cambridge University.
Kunesh, J.G., Lahm, L., and Yanagi, T., Commercial scale experiments that provide insight
on packed tower distributors, Industrial and Engineering Chemistry Research 26(9):
1845, 1987.
Long, C. and Sayma, N., Heat Transfer: Exercises, 1st edn., Brighton, UK: Ventus Publishing
ApS, 2010. http://www.leka.lt/sites/default/files/dokumentai/heat-transfer.pdf. Retrieved
June 12, 2013.
Mellon, C., MPEC problem formulations in chemical engineering applications, http://numero.
cheme.cmu.edu/uploads/rp.pdf, Chemical Engineering Department. Retrieved on June
16, 2007.
Metcalf & Eddy, Wastewater Engineering: Treatment and Reuse, 4th edn., New York:
McGraw-Hill, 2003, p. 1449.
Misterguch, The six types of chemical reactions, 2012, http://misterguch.brinkster.net/
6typesofchemicalrxn.html.
Moore, F. and Rukovena, F., Random packing, vapor and liquid distribution: Liquid and gas
distribution in commercial packed towers, Chemical Plants and Processing 11–15,
August 1987.
Needham, J., Science and Civilisation in China, Cambridge, U.K.: Cambridge University
Press, 1980.
Nonweiler, T.R., Computational Mathematics: An Introduction to Numerical Approximation,
New York: John Wiley & Sons, 1986.
Othmer, D.F., Distillation—Some steps in its development, in W.F. Furter (ed.), A Century of
Chemical Engineering; Separation of liquid–liquid mixtures (solutions), IUPAC, 1982.
Palm, W.J., Introduction to MATLAB (R) for Engineers, 3rd edn., New York: McGraw-Hill,
2011.
Perry, R.H. and Green, D.W., Perry’s Chemical Engineers’ Handbook, 6th edn., New York:
McGraw-Hill, 1984.
Perry, R.H. and Green, D.W. (Eds.), Perry’s Chemical Engineers’ Handbook, 7th edn.,
New York: McGraw-Hill, 1997.
Plawsky, J.L., Transport Phenomena Fundamentals, New York: CRC Press, April 2001,
pp. 1–3.
Press, W.H. et al., Numerical Recipes in C: The Art of Scientific Computing, 3rd edn.,
Cambridge, U.K.: Cambridge University Press, 2007.
Rice, R.G. and Do, D.D., Applied Mathematics and Modeling for Chemical Engineers,
New York: John Wiley & Sons, Inc., 1995.
Roy, G.K., Solved Examples in Chemical Engineering, Delhi, India: Khanna Publishers,
2010.
Schaschke, C.J., Solved Practical Problems in Fluid Mechanics, Dundee, Scotland: Abertay
University.
Seader, J.D. and Henley, E.J., Separation Process Principles, New York: Wiley, 1998.
Shalabi, M.A., Harji, B.H., and Kenney, C.N., Journal: Cheminform 27(28), 2010.
Spiegel, L., A new method to assess liquid distributor quality, Chemical Engineering and
Processing 45(11): 1011, 2006.
Study Guides and Strategies, Welcome to Problem Solving with Algorithms and Data
Structures, n.d., http://www.studygs.net/scimethod.htm, retrieved on December 4,
2014.
Subramanian, M., Process Calculations—Solved Problems in Chemical Engineering, n.d.,
http://www.msubbu.in/sp/pc/, retrieved on September 13, 2014.
Bibliography 249
251
252 Index