Overall PDF

See discussions, stats, and author profiles for this publication at: https://www.researchgate.
net/publication/267811360
TEACHING MODULE: MATERIAL BALANCES, ENERGY BALANCES AND HEAT

TRANSFER
Book · September 2010

DOI: 10.13140/2.1.4401.5049
CITATIONS READS
0 14,682
3 authors, including:
Norazian Mohamed Noor Masitah Hasan

Universiti Malaysia Perlis Universiti Malaysia Perlis
64 PUBLICATIONS 281 CITATIONS 7 PUBLICATIONS 340 CITATIONS
SEE PROFILE SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Variation of particulate matter (PM10) in Malaysia View project
Investigation of adsorption mechanism of elephant apple (Dilelnia Indica) fruit as an adsorbent for the removal of emerging pharmaceutical contaminants View project
All content following this page was uploaded by Norazian Mohamed Noor on 06 November 2014.
The user has requested enhancement of the downloaded file.

CHAPTER 1
INTRODUCTION TO ENGINEERING CALCULATIONS
1.1 Introduction
What do chemical engineers do? Although their backgrounds and professional skills are
similar, chemical engineers work in a wide variety of industries, in addition to chemicals
and petroleum, such as:
Biotechnology Lime and cement

Consulting Man-made fibers
Drugs and pharmaceuticals Metallurgical and metal products
Fats and oils Paints, varnishes, and pigments
Fertilizer and agricultural chemicals Pesticides and herbicides
Environment Waste water treatment
All the industries as mentioned above are involving numerous of chemical process unit.
Therefore, chemical engineers play an important role on design, operation, control,
troubleshooting, research and management in the chemical process. Chemical process is a
combination of process equipment designed to efficiently convert raw materials into
finished or intermediate products. Figure 1 shows the example of chemical processes
converting the raw material into desired product.
Figure 1.1: Chemical process
1.2 Units and Dimensions
“What are units and dimensions and how do they differ?”
Dimensions are basic concepts of measurement such as length (L), mass (M), time (t),
temperature (T), amount of substance (n) and so on. Besides, units are the mean of
expressing the dimensions as feet or centimeters for length, or hours or seconds for time.
By attaching units to all numbers that are not fundamentally dimensionless, you are able
to easy interrelating the physical meaning to the numbers use. Moreover, a logical
1
approach to the problem rather than remembering a formula and plugging numbers could
also help the chemical engineers in engineering calculation. SI units are universally
accepted for engineering calculation. Thus, American engineering system (based on
British standards) is still used extensively in the U.S.
Example 1.1
What are the dimensions of mass flux (mass flow rate per unit area perpendicular to the
flow)?
1 dm M
G= dimensions are
A dt L2t
The rules for handling units are essentially quite simple by addition, subtraction or
equality.
¾ Values could be added if UNITS are the same.

¾ Values cannot be added if DIMENSIONS are different.
Example 1.2
(i) 6 ft + 10 0C =???
* Different dimensions: length, temperature -- cannot be added
* Same dimension: length, different units -- can add
2
Table 1.1: SI Units
Symbol Definition of
Physical Quantity Name of Unit
for Unit* Unit
Basic SI Units
Length metre, meter m
Mass kilogramme, kilogram kg
Time second s
Temperature kelvin K
Amount of
mole mol
substance
Derived Sl Units
Energy joule J kg.m2.s-2
Force newton N kg.m.s-2 = J.m-1
Power watt W kg.m2.s-3 = J.s-1
Density kilogram per cubic meter kg.m-3
Velocity meter per second rn.s-1
Acceleration meter per second squared rn.s-2
newton per square meter,
Pressure N.m-2, Pa
pascal
Heat Capacity joule per (kilogram -
J.kg-1 ,K-1
kelvin)
Alternative Units
Time minute, hour, day, year min, h, d, y
Temperature degree Celsius °C
Mass tonne, ton (Mg), gram t, g
Volume litre, liter (dm3) L
Table 1.2: American Engineering System Units
Physical Quantity Name of Unit Symbol

Basic Units
ft
Length feet
Mass pound (mass) lbm
Force pound (force) lbf
Time second, hour s, hr
Temperature degree Rankine °R
Derived Units
Energy British thermal unit, foot pound (force) Btu, (ft)(lbf)
Power horsepower hp
Density pound (mass) per cubic foot lbm/ft3
Velocity feet per second ft/s
Acceleration feet per second squared ft/s2
Pressure pound (force) per square inch lbf/in2
Heat capacity Btu per pound (mass) per degree F Btu/lbm.0F
3
1.3 Conversion of Units
Conversion factors are statements of equivalent values of different units in the same
system or between systems of units. The concept is to multiply any number and its
associated units with dimensionless ratios termed conversion factors to arrive at desired
answer and its associated units. The factors for conversion units are show in table 1.3.
Table 1.3: Factors for unit conversions
Example 1.3
Convert an acceleration of 1 cm/s2 to its equivalent in km/yr2.
1cm 1m 1km 3600 2 s 2 24 2 hr 2 365 2 day 2

× × × × × 2 2
s 2 100cm 1000m 1hr 2 1day 2 1 yr
= 9.95 × 10 9 km / yr 2
4
Do It Yourself:
Convert 400 in3/day to cm3/min.

(Answer: 4.56 cm3/min)
1.4 Processes and Process Variables
A process is any operation or series of operations by which a particular objective is

accomplished. Those mentioned operations are involving a physical or chemical change
in a substance or mixture of substances. The material that enters a process is referred to
as the input or feed, and that which leaves is the output or product. Therefore, several
process variables are associated through input or output of a process stream.
1.3.1 Instructional Objectives
The objectives in studying this section are to be able to:
1. Draw a simple block flow diagram representing a process, showing input and
output streams, and essential process variables.
2. Calculate the quantities of mass (or mass flow rate), volume (or volumetric flow
rate), and moles (or molar flow rates) from a knowledge of the third quantity for
any species of known density and molar mass.
3. Explain:
(a) The meaning of gram-mole, lb-mole, mol and kmol;
(b) At least two methods of measuring temperature and at least two for
measuring fluid pressure;
(c) The meaning of the terms absolute pressure and gauge pressure.
4. Convert a pressure expressed as a head of a fluid to the equivalent pressure

expressed as a force per unit area, and vice versa
5. Convert a manometer reading into a pressure difference for an open end

manometer, a sealed end manometer and a differential manometer.
6. Convert among temperatures expressed in K, °C, °F and °R.
5
1.3.2 Process
A process is any operation or series of operations that cause a physical or chemical

change in a substance or mixture of substances. Figure 1.2 shows a process stream with
several examples of process variables of input and output. Meanwhile, the details about
process variables are stated in table 1.4.
Figure 1.2: Process streams with various process variables.
Table 1.4: Process variables
6
1.3.3: Mass and Volume
The density of a substance is the mass per unit volume of the substance (kg/m3, g/m3,
lb/ft3, etc). The specific volume of a substance is the volume occupied by a unit mass of
the substance; it is the inverse of density. Densities of pure solids and liquids are
essentially independent of pressure and vary relatively slightly with temperature.
Densities of many pure compounds, solutions and mixtures can be found in standard
references. The density of a substance can be used as a conversion factor to relate the
mass and the volume of a quantity of the substance.
Example 1.4
The density of carbon tetrachloride is 1.595 g/cm3, the mass of 35 cm3 of CCl4 is
therefore,
35 cm3 1.595 g = 55.825 g

cm3
and the volume of 9.3 lbm of CCl4 is
9.3 lbm 454 g cm3

= 2647.47 cm3
1 lbm 1.595 g
The specific gravity of a substance is the ratio of the density ρ of the substance to the
density ρref of a reference substance at a specific condition:
ρ (1.1)
SG =
ρ ref
The reference most commonly used for solids and liquids is water at 4 oC, which has the
following density:
ρH2O (4 °C) = 1 g/cm3

= 1000 kg/m3
= 62.43 lbm/ft3
If you are given the specific gravity of a substance, multiply it by the reference density in
any units to get the density of a substance in the same units. Special density units called
degrees Baumé (°Bé), degrees API (°API) and degrees Twaddell (°Tw).
7
Example 1.5
Calculate the density of mercury in lb/ft3 from a tabulated specific gravity, and calculate
the volume in ft3 occupied by 215 kg of mercury.
(Given that the specific gravity of mercury at 20 oC as 13.546)
ρ Hg = S .G × ρ ref
= 13.546 × 62.43 lb / ft 3
= 845.67 lb / ft 3
Volume = 215 kg lbm ft3 = 0.56 ft3
0.454 kg 845.67 lb
Do It Yourself:
A liquid has a specific gravity of 0.5.

i. What is its density in g/cm3? (Answer: 0.5 g/cm3)
ii. What is its specific volume in cm3/g? (Answer: 2 cm3/g)
iii. What is its density in lb/ft3? (Answer: 0.032 lb/ft3)
iv. What is the mass of 3 cm3 of this liquid? (Answer: 1.5 g)
v. What volume is occupied by18g? (Answer: 36 cm3)
1.3.4: Flow rate
Most processes involve the movement of material from one point to another. The rate at
which a material is transported through a process line is the flow rate of that material.
The flow rate of a process stream can be expressed as a mass flow rate (mass/time) or as
a volumetric flow rate (volume/time) as show in figure 1.3 as below.
Figure 1.3: Flow rate
8
Example 1.6
The volumetric flow rate of CCl4 in a 1.0-cm-diameter pipe is 100 cm3/min.

(Given that the molecular weight of CCl4 = 153.838 and density of ρ CCl4 = 1.595 g/cm3.)
i. What is the mass flow rate?

• •
3 3
m = V ρ = 100 cm /min x 1.595 g/cm = 159.5 g/min
ii. What is the molar flow rate?
= 159.5 g-CCl4/min x g mol-CCl4/153.838 g-CCl4
= 1.034 g mol-CCl4/min
iii. What is the linear velocity of CCl4?
m 159.5 g / min
v= = = 127.32cm / min
Aρ π (1cm) 2
×1.595 g / cm 3
4
Do It Yourself:
The mass flow rate of n-hexane (ρ=0.659 g/cm3) in a pipe is 6.59 g/s.
(a) What is the volumetric flow rate of the hexane?

(b) What is the linear velocity of hexane in the pipe with internal diameter of 5 cm?
(Answer: (a) 10 cm3/s; (b) 0.509 cm/s)
1.3.5: Chemical Composition
Most materials encountered in nature and in chemical process systems are mixtures of
various species. The physical properties of a mixture depend strongly on the mixture
composition. In this section we will review different ways to express mixture
compositions and also outline the methods of estimating physical properties of a mixture
from the properties of the pure components.
1.3.5(a): Moles and Molecular Weight
The atomic weight of an element is the mass of an atom on a scale that assigns. The
atomic weights of all the elements in their naturally isotopic proportions are listed in the
table given. The molecular weight of a compound is the sum of the atomic weights of the
atoms that constitute a molecule of the compound: For example atomic oxygen (O) has
an atomic weight of approximately 16 and therefore molecular oxygen (O2) has a
molecular weight of approximately 32. A gram-mole (g-mole or mol in SI units) of a
9
species is the amount of that species whose mass in grams is numerically equal to its
molecular weight. Hence, one g-mole of any species contains approximately 6.02x1023
(Avogadro’s number) molecules of that species.
Moreover, there are other types of moles such as kg-moles, lb-moles, and ton-moles. For
example: Carbon monoxide (CO) has a molecular weight of 28; 1 mol of CO therefore
contains 28g, 1 lb-mole contains 28 lbm, 1 ton-mole contains 28 tons and so on. Besides,
the same factors used to convert masses from a unit to another can also be used to convert
the equivalent molar units: there is 454 g/lb for example, and therefore there is 454
mol/lb-mole, regardless of the substance involved.
100 g CO2 1 mol CO2

44.01 g CO2
Example 1.7
How many of each of the following are contained in 100g of CO2 (M=44.01)?
i. Mol CO2
= 2.273 mol CO2
ii. lb-moles CO2
2.273 mol CO2 1 lb-mol = 5.011 x 10-3 lb-mole CO2

453.6 mol
Each molecules of CO2 contains one atom of C, one molecule of O2 or two atoms of O.
Therefore, each 6.02x1023 molecules of CO2 (1mol) contains 1 mol C, 1 mol O2, or 2 mol
O. Thus,
iii. mol C
2.273 mol CO2 1 mol C = 2.273 mol C

1 mol CO2
10
iv. mol O
2.273 mol CO2 2 mol O = 4.546 mol O

1 mol CO2
v. mol O2
2.273 mol CO2 1 mol O2 = 2.273 mol O2

1 mol CO2
vi. gram O
4.546 mol O 16 g O = 72.7 g O

1 mol O
vii. gram O2
2.273 mol O2 32 g O2 = 72.7 g O2

1 mol O2
viii. molecules of CO2
2.273 mol CO2 6.02 x 1023 molecules = 1.37 x 1024 molecules

1 mol
Do It Yourself:
Calcium carbonate is a naturally occurring white solid used in the manufacture of lime
and cement. Calculate the number of lb mol of calcium carbonate in:
(a) 50 g mol of CaCO3.

(b) 150 kg of CaCO3.
(c) 100 lb of CaCO3.
(Answer: (a) 0.11 lb mol; (b) 3.3 lb mol; (c) 1 lb mol)
1.3.5 (b): Mass and Mole Fractions
Process streams occasionally contain one substance, but more often they consist of
mixtures of liquids or gases, or solution s of one or more solutes in a liquid solvent. The
following terms can be used to define the composition of a mixture of substances,
including a species A.
11
Example 1.8
A solution contains 15% A by mass (xA = 0.15) and 20 mole % B (yB = 0.20)
a) Calculate the mass of A in 175 kg of the solution.
175 kg solution 0.15 kg A = 26.25 kg A

kg solution
b) Calculate the mass flow rate of A in a stream of solution flowing at a rate of 53 lb/h.
53 lb 0.15 lb A = 7.95 lb A/ h
h lb
c) Calculate the molar flow rate of B in a stream flowing at a rate of 1000 mol/min.
1000 mol 0.2 mol B = 200 mol B/min

min mol
d) Calculate the total solution flow rate that corresponds to a molar flow rate of
28 kmol B/s.
28 k mol B 1 k mol solution = 140 kmol solution/s

s 0.2 k mol B
e) Calculate the mass of the solution that contains 300 lb of A.
300 lb A 1 lb solution = 2000 lb solution

0.15 lb A
12
1.3.5 (c): Concentration
The mass concentration of a component of a mixture or solution is the mass of this

component per unit volume of the mixture.
The molar concentration of a component is the number of moles of the component per
unit volume of the mixture.
The molarity of a solution is the value of the molar concentration of the solute expressed
in g-moles solute / liter solution.
Example 1.9
A 0.5 molar aqueous solution of sulfuric acid flows into a process unit at a rate of 1.25
m3/min. The specific gravity of the solution is 1.03 and molecular weight of sulfuric acid
is 98.08. Calculate:
a) The mass concentration of H2SO4 in kg/m3

0.5 mol H 2 SO4 98 g 1kg 1000 L
= × × × = 49 kg H 2 SO4 / m 3
L mol 1000 g 1m 3
b) The mass flow rate of H2SO4 in kg/s

49 kg H 2 SO4 1.25 m 3 1 min
= × × = 1.02 kg H 2 SO4 / s
m3 min 60 s
c) The mass fraction of H2SO4
The mass fraction of H2SO4 equals the ratio of the mass flow rate of H2SO4 to the total
mass flow rate, which can be calculated from the total volumetric flow rate and the
solution density.
1000kg
ρ Solution =1.03 × ( 3
) = 1030kg / m 3
m
13
kg 1.25m 3 solution 1030kg 1 min
m solution ( s ) = min
× 3 ×
m solution 60 s
= 21.46kg / s
m H 2 SO4 1 kg H 2 SO4 / s
x H 2 SO4 = = = 0.048 kg H 2 SO4 / kg solution
m solution 21.46 kg solution / s
1.3.6: Pressure
A pressure is the ratio of a force to the area on which the force acts. Pressure units are
force units divided by area units such as N/m2 or Pascal (Pa), dynes/cm2, and lbf/in2 or
psi.
Hydrostatic pressure = pressure at the base of a fluid column
P=P +ρgh
0
P-P = ρ g h
0
If P is atmospheric pressure,
0
then P-P = ρ g h is called the gauge
0
pressure, and P is the absolute pressure.
P =P +P
abs gauge atm
Fluid Pressure Measurement
Most common pressure measuring devices are stated in Figure 1.4. Bourdon gauge
manometers can show measurement nearly perfect vacuums to about 700 atm.
Meanwhile, manometers only can measure pressures below about 3 atm.
14
Figure 1.4: Pressure measurement device.
Manometer principle is showing in figure 1.5. The fluid pressure must be the same at any
two points at the same height in a continuous fluid.
Figure 1.5: Manometer principle
General manometer equation:
P1 + ρ1gd1 = P2 + ρ2gd2 + ρfgh
Differential manometer equation:
P1 – P2 = (ρf –ρ) gh, since ρ1 = ρ2
15
1.3.7: Temperature
Temperature is a measurement of the average kinetic energy possessed by the substance

molecules. It must be determined indirectly by measuring some temperature-dependent
physical properties of another substance. The temperature measuring devices are:
a) Resistance thermometer (by means of electrical resistance of a conductor)

b) Thermocouple (by voltage at the junction of two dissimilar metals)
c) Pyrometer (by spectra of emitted radiation)
d) Thermometer (by volume of a fixed mass of fluid)
The temperature conversions are:
T(K) = T(0C) + 273.15; T(0R) = T(0F) + 459.67;
T(0R) = 1.8 T (K); T(0F) = 1.8 T(0C) + 32
Problems:
1. Convert the following to the desired units:

(a) 60 mi/hr to m/s (b) 30 N/m2 to lbf/ft2
(c) 16.3 J to Btu (d) 4.21 kW to J/s
2. Change the following to the desired units:

(a) 235 g to pounds. (b) 610 L to cubic feet.
(c) 30 g/L to pounds/cubic feet. (d) 14.7 lb/in2 to kg/cm2
3. Convert the following quantities to the ones designated:

(a) 42 ft2/hr to cm2/s.
(b) 25 psig to psia.
(c) 100 Btu to hp-hr.
4. The specific gravity of a fuel oil is 0.82.

(a) What is the density of oil in lb/ft3?
5. The density of a liquid is 1500 kg/m3 at 20 °C.

(a) What is the specific gravity 20°C/4°C of this material?
(b) What volume (ft3) does 140 lbm of this material occupy at 20°C?
6. Silver nitrate (lunar caustic) is a white crystalline salt, used in marking inks, medicine
and chemical analysis. How many kilograms of silver nitrate (AgNO3) are there in:
(a) 13.0 lb mol AgNO3.
(b) 55.0 g mol AgNO3
16
7. Complete the table below with the proper equivalent temperatures.
°C °F K °R
- 40.0
77.0
698
69.8
17
CHAPTER 2
FUNDAMENTALS OF MATERIAL BALANCES
2.1 Introduction
Material balances are important first step when designing a new process or analyzing an
existing one. They are almost always prerequisite to all other calculations in the solution
of process engineering problems. Material balances are nothing more than the application
of the law of conservation of mass, which states that mass can neither be created nor
destroyed. Thus, as an example, you cannot specify an input to a reactor of one ton of
naphtha and an output of two tons of gasoline or gases or anything else. One ton of total
material input will only give one ton of total output, i.e. total mass of input = total mass
of output.
A material balance is an accounting for material. Thus, material balances are often
compared to the balancing of current accounts. They are used in industry to calculate
mass flow rates of different process streams entering or leaving chemical or physical
processes.
2.2 Process Classification
A. Based on manner of delivering process streams
Batch process
The feed is charged into a vessel at the beginning of the process and the vessel contents
are removed some time later. No mass transfer across the boundary between the time the
feed is charged and the time the product is removed.
Example: Rapidly add reactants to a tank and remove the products and unconsumed
reactants sometime later when the system has come to equilibrium.
Continuous process
Inputs and outputs flow continuously through the duration of the process.
Example: Pump a mixture of liquids into a distillation column at constant rate and
steadily withdraw product streams from the top and bottom of the column.
Semi batch process
Any process that is neither batch nor continuous.
Example: Allow the contents of pressurized gas container to escape to the atmosphere;
slowly blend several liquids in a tank from which nothing is being withdrawn.
18
B. According to variation of the process variables with time
Steady state
The values of all process variables (i.e. all temperatures, pressures, compositions, flow
rates) do not change with time. For example, batch or semi batch processes.
Unsteady state or transient
Any of the process variables change with time. Continuous processes may be either
steady-state or transient.
2.3 The General Balance Equation
Figure 2.1: Process balance
A balance on a conserved quantity (total mass, mass of a particular species, energy, and
momentum) in a system as show in figure 2.1 can be written in the following general way:
Input + generation - output - consumption = accumulation [2.1]
where,
Input - enters through system boundary
Generation - produced within system
Output - leaves through system boundary
Consumption - consumed within system
Accumulation - builds up within system
Example 2.1
Each year 50000 people move into a city, 75000 people move out, 22000 are born, and
19000 die. Write a balance on the population of the city.
Input + generation - output - consumption = accumulation
50000 + 22000 – 75000 – 19000 = - 22000
Therefore, the city’s population decreases by 22,000 people each year.
19
Two types of balances:
A. Differential balance
¾ Balance that indicates what is happening in a system at an instant in time.

¾ Each term of the balance equation is a rate (rate of input, rate of generation, etc)
and has units of the balanced quantity unit divided by the time unit (people/yr,
barrels/day).
¾ This is the type of balance usually applied to a continuous process.
B. Integral balance
¾ Balance that describes what happens between two instants of time.

¾ Each term of the equation is an amount of the balanced quantity and has the
corresponding unit (people, barrels).
¾ This type of balance is usually applied to a batch process, with the two instants of
time being the moment after the input takes place and the moment before the
product is withdraw.
The following rules can be used to simplify the material balance equation:
¾ If the balanced quantity is total mass, generation = 0 consumption = 0:
Input - output = accumulation
¾ If the balanced substance is non-reactive species, generation = consumption = 0:
Input - output = accumulation
¾ If a system is at steady state, accumulation = 0:
Input + generation = output + consumption
2.4 Balance on Continuous Steady-State Process
For continuous processes at steady state, the accumulation term in the general balance
equation, (Eq 2.1), equals zero and the equation simplifies as below:
Input + generation = output + consumption [2.2]
If the equation is on a non-reactive species or on total mass, the generation and

consumption terms equal zero and the equation reduces to the equation as below:
Input = Output [2.3]
20
Example 2.2
One thousand kilograms per hour of mixture of benzene (B) and toluene (T) containing
50% benzene by mass is separated by distillation into two fractions. The mass flow rate
of benzene in the top stream is 450 kg B/h and that of toluene in the bottom stream is 475
kg T/h. The operation is at steady state. Write balance on benzene and toluene to
calculate the unknown component flow rates in the output streams.
The process can be depicted schematically as follows:
450 kg B/h
m1 (kg T/h)
500 kg B/h
500 kg T/h
m2 (kg B/h)
475 kg T/h
Since the process is at steady state there can be no buildup of anything in the system, so
the accumulation term equals zero in all material balances. In addition, since no chemical
reactions occur, there can be no nonzero generation or consumption terms. For all
balances, Eq 2.3 are apply here:
Input = Output
For Benzene balance:

500 kg B/h = 450 kg B/h + m2
m2 = 50 kg B/h
For Toluene balance:

500 kg T/h = m2 + 475 kg T/h
m1 = 25 kg T/h
Hence,
Total Mass Balance:

1000 kg/h = 450 kg/h + m1kg/h + m2 kg/h + 475 kg/h
1000 kg/h = 1000 kg/h
21
2.5 Integral Balance on Batch Process
Ammonia is produced from nitrogen and hydrogen in a batch reactor. At time t=0 there
are n0 mol of NH3 in the reactor, and at later time tf the reaction terminates and the
contents of the reactor, which include nf mol of ammonia are withdraw. Between t0 and tf
no ammonia enters or leaves through the reactor boundaries, so general balance equation
(1) is simply generation = accumulation. Moreover, the quantity of ammonia that builds
up (accumulates) in the reactor between t0 and tf is simply nf – n0, the final amount minus
the initial amount. The same reasoning may be applied to any substance participating in a
batch process to obtain the following equation:
Accumulation = final output – initial input (by definition)

= generation – consumption
Therefore, equating these two expressions for the accumulation yields:
Initial input + generation = final output + consumption [2.4]
Example 2.3: Balance on a Batch Mixing Process
Two methanol water mixtures are contained in separate flasks. The first mixture contains
40 wt% methanol, and the second contains 70 wt% methanol. If 200g of the first mixture
is combined with 150g of the second, what are the mass and composition of the product?
200g
0.4g CH3OH/g m (g)
0.6g H2O/g
x (g CH3OH/g)
(1 –x )(g H2O/g)
150g
0.7g CH3OH/g
0.3g H2O/g
Observe that the input and output ‘streams’ shown on the chart denote the initial and final
states for this batch process. Since no reactions are involved, the generation and
consumption terms of Eq (4) can be omitted so that all balance have the simple form
‘input = output’.
Total Mass Balance

200g + 150g = m
m = 350g
22
Methanol Balance:
200 g × 0.4 gCH 3 OH 150 g × 0.7 gCH 3 OH m ( g ) × x ( gCH 3OH )

( )+( )=( )
g g (g)
x = 0.529 g CH3OH/g
2.6 Integral Balance on Semi-batch Process
Integral balances can also be written for semi-batch and continuous processes. The
procedure is to write a differential balance on the system and then to integrate it between
two instants of time. In most cases the required calculations are more complex than those
we have seen so far. However, some problems of this type are relatively straight forward
solutions.
Example 2.4
Air is bubbles through a drum of liquid hexane at a rate of 0.1 kmol/min. The gas stream
leaving the drum contains 10 mole% hexane vapors. Air may be considered insoluble in
liquid hexane. Use an integral balance to estimate the time required to vaporize 10m3 of
the liquid.
0.1 kmol /min n
0.1 kmol C6H14/kmol

0.9 kmol air/kmol
We begin with a differential balance on air. Since we assume that air neither dissolves in
the liquid nor reacts with hexane in the process unit, the balance reduces to input =
output:
0.1 kmol air 0.9 kmol air n kmol

= ×
min kmol min
n = 0.111 kmol / min
We next write an integral hexane balance, proceeding from time t=0 to tf (min), the time
we seek to calculate. The balance has the form accumulation = - output. The
accumulation term, which is the total change in the moles of liquid hexane in the system
during time tf, must be negative since hexane is being lost from the system. Since the
23
total number of moles of hexane evaporated occupied a liquid volume of 10 cubic meters
and the specific gravity of liquid hexane is 0.659, the accumulation term equals:
0.659 kg 1 kmol
Δn = −10 m 3 × × = −76.45 kmol C 6 H 14
m3 86.2 kg
The output term in the balance is the rate at which hexane is leaving the system [0.1n
(kmol C6H14/min)] times the total process time, tf (min). The balance (accumulation = -
output) is therefore;
- 76.45 kmol C6H14 = -0.1n tf

n = 0.111 kmol /min
tf = 6887 min
2.7 General procedure for material balance calculations
1) Choose a basis of calculation an amount or flow rate of one of the process streams.
¾ If no stream amount or flow rate is specified in the problem statement, take as a
basis an arbitrary amount or flow rate of a stream with a known composition.
2) Draw the flowchart of the process, using boxes or other symbols to represent the
process units, and lines with arrows to represent input and output streams.
3) Fully label the chart when it is first drawn with values of known process variables
and symbols for unknown variables being written for each input and output stream.
4) Do the degree-of-freedom analysis.
5) If the problem is solvable, the starting balance should be an equation with minimum
number of unknowns, preferably with only one unknown.
6) After the one unknown in step 4 has been calculated, use that value to give an
equation with one variable for another unknown.
7) As each unknown variable is determined, its value is filled so that the chart provides
a continuous record of where the solution stands and what must still be done.
Notes:
¾ The maximum numbers of independent equations that can be derived by
writing balances on a non-reactive system equal the number of chemical
species in the input and output streams.
¾ The additional equations can be written from the fact that the component mass or
mole fractions of a stream adds up to 1.0 (Σ xi = 1.0).
(You can use these constrains to reduce the number of unknowns in the flow
chart (step 2 above) to a minimum.)
24
¾ If you are given that the mass of stream 1 is half that of stream 2, label the masses
of these streams m and 2m rather than m1 and m2; if you know that there is three
times as much nitrogen (by mass) in a stream as oxygen, label the mass fractions
of O2 and N2 y and 3y rather than y1 and y2.
(This can also reduce the number of unknowns in step 2 above, or can
provide additional equation needed for solution.)
¾ Balances are not normally written on volumetric quantities. If a volumetric flow
rate is given, it is generally useful to calculate the mass or molar flow rate for the
balance.
Example 2.5
An experiment on the growth rate of certain organisms requires an environment of humid

air enriched in oxygen. Three input streams are fed into an evaporation chamber to
produce an output stream with the desired composition.
A: Liquid water, fed at a rate of 20 cm3/min

B: Air (21 mole% O2, the balance N2)
C: Pure oxygen, with a molar flow rate one fifth of the molar flow rate of stream B.
The output gas is analyzed and is found to contain 1.5 mole% of water. Draw and label a
flowchart of the process, and calculate all unknown stream variables
0.2 n1 (mol
n3
n1 (mol air/min) 0.015 mol H2O/mol

y mol O2/mol
(0.985 – y) (mol N2/mol)
0.21g mol O2/mol
0.79 mol N2/mol
20 cm3 H2O/min
n2 (mol H2O/min)
Notes on the labeling:
1. Since the one known flow rate (20 cm3 H2O/min) is given on a per minute basis, it
is most convenient to label all stream flow rates on this basis.
2. Once the variable name (n1) is chosen for the air flow rate, the given information
about the ratio of the air and O2 flow rates may be used to label the O2 flow rate
0.2n1.
3. The mole fractions of the components of any stream must add up to 1. Since the
mole fraction of H2O in the outlet stream is known to be 0.015, once the mole
25
fraction of O2 is labeled y, that of N2 must be 1 – (y + 0.015) = (0.985 – y) (mol
N2/mol).
The quantity n2 can be calculated from the given volumetric flow rate and the density of
liquid water:
20 cm 3 H 2 O 1 g 1 mol
n2 = × 3× = 1.11 mol H 2 O / min
min cm 18.02 g
The three remaining unknowns (n1; n2 and y) can be determined from balance, all of
which have the simple form input = output for this non-reactive steady state process. The
balances are easily written by referring to the flow chart.
Water balance:
n2 (mol H2O/min) = n3 (0.015 mol H2O/min)
n3 = 74 mol/min
Total Mole Balance:
0.2n1 + n1 + n2 = n3
n1 = 60.74 mol/min
Nitrogen Balance:
n1 mol 0.79 mol N 2 n3 mol (0.985 − y ) mol N 2

× = ×
min mol min mol
∴ 0.79 n1 = n3 (0.985 − y )
y = 0.337 mol O2 / mol
Do It Yourself:
In the steady state flow process shown in the figure below, a number of values of m
(mass fraction) are not given. Mary says that nevertheless the problem has a unique
solution for the unknown values of m. Kelly says that 4 values of m are missing, that you
can write 3 component material balances, and that you can use 3 relations for ∑mi =1,
one for each stream, a total of 6 equations, so that a unique solution is not possible. Who
is right?
26
F = 10kg P = 16kg
m1 = 0.175
m1 = 0.01 m2 =?
m2 =? m3 =?
A = 6kg m1 = 0.30
m2 =?
m3 = 0.20
(Answer: Stream F, m2 = 0.90; stream A, m2 = 0.5; stream P, m2 = 0.75, m3 = 0.075)
2.7 Balancing a process
Consider the following process flowchart for a separation unit at steady state:
Number of unknowns = 7
Number of mass balance equations:
4 component balances: F x = D x + W x for i = A, B, C, and D

iF iD iW
1 overall balance: F = D + W
3 mass fraction constraints of streams: Σ x = 1.0, Σ x = 1.0, Σ x = 1.0
iF iD iW
Total number of equations = 8. They are not all independent, because the summing of the
four component balances and then using the mass fraction constraint of each stream
results in the overall balance. Thus the number of independent equations = 7.
Since number of unknowns = number of independent equations, the problem is solvable.
Then how will you solve the problem in the easiest way?
27
Example 2.6
An aqueous solution of sodium hydroxide contains 20% NaOH by mass. It is desired to

produce an 8.0% NaOH solution by diluting a stream of the 20% solution with a stream
of pure water. Calculate the ratios (liters H O needed / kg feed solution) and (kg product
2
solution / kg feed solution).
Solution:
Basis: 100 kg of 20% feed solution
Step 2:
Step 3 (in conjunction with Steps 1 and 2):
Degree of freedom analysis:

Number of unknowns needed to be solved = 2 (V in terms of m , and m )
1 1 2
Since it is a nonreactive process and involves two species,
Number of independent equations = 2
So, degree of freedom = 2-2 = 0. Thus the problem is solvable.
Step 4 (solution procedure):
Number of mass balance equations involved: 2
The balance equation that contains only one unknown is NaOH balance since it is not
present in pure water stream. So choose it.
NaOH mass balance: 100 x 0.2 = 0.08 m2

m2 = 250 kg
The remaining two balances (H2O and overall) contain 2 unknowns, and simplest of these
is overall balance, so choose it.
Overall mass balance: 100 + m1 = m2

So, m1 = m2 – 100
= 250 – 100 = 150 kg
28
Volume of diluents water, V = m /ρ = 150 kg x 1.0 liter/kg = 150 liter
1 1 w
Thus,
V1/100 kg = 150 liter/100 kg = 1.5 liters H2O/kg feed solution
m2/100 kg = 250 kg/100 kg = 2.5 kg product solution/kg feed solution
Do It Yourself:
A solution composed of 50% ethanol (EtOH), 10% methanol (MeOH), and 40% water
(H2O) is fed at the rate of 100 kg/hr into a separator that produces one stream at the rate
of 60 kg/hr with the composition of 80% EtOH, 15% MeOH, and 5% 5% H2O, and a second
stream of unknown composition.
(a) Sketch and completely label a flow chart of the process.
(b) Calculate the composition (in %) of the three compounds in the unknown
stream and its flow rate in kg/hr.
Answer: (a)
100 kg/hr 60 kg/hr
50 % EtOH 80 % EtOH
10 % MeOH 15 % MeOH
40 % H2O 5 % H2O
A kg/hr
x % EtOH
y % MeOH
z % H2O
Answer: (b) A = 40 kg/hr; x = 5%; y = 2.5%; z = 92.5%
2.8 Balances on Multiple Unit Processes
Plants in the process industries are comprised of many interconnected units. However,
you can easily apply the same techniques discussed in the previous sections to solve
material balance problems in such plants.
List and count up the number of variables whose values are unknown, making sure you
do not count the same variable more than once, and then list and count up the number of
independent balances that you can make, making sure that balances for one unit do not
render formerly independent balances for another unit into dependent balances.
If the number of independent equations equals the number of variables whose values are
unknown, at least for a set of linear equations, you can generally solve the equations for a
unique answer.
29
If you ignore all the internal streams and variables within a set of connected subsystem,
you can treat the overall system exactly as you treated a single system by drawing a
boundary about the entire set of subsystem.
A flow chart for a two unit process is show as below.
Feed 2
A
C E
B
Unit 1 D Unit 2
Feed 1 Product
Product Product Feed 3
Five boundaries drawn above portions of the process define systems on which balances
can be written. Boundary A encloses the entire process; the system defined by this
boundary has as inputs Feed Streams 1, 2 and 3 and Product Stream 1, 2 and 3.Balances
on this system are referred to as overall balances. The stream that connects Unit 1 and 2
is internal system are referred to this system and so would not either into overall system
balances.
Boundary B encloses a feed stream mixing point. Feed Streams 1 and 2 are inputs to this
system and the stream flowing to Unit 1 is an output. Boundary C encloses Unit 1 (one
input stream and two output streams). Boundary D encloses a stream splitting point (one
input stream and two output streams). Boundary E encloses Unit 2 (two input streams and
one output stream).
Example 2.7
A labeled flow chart of a continuous steady state process is shown below. Each stream
contains of two components; A and B, in different proportions. Three streams whose flow
rates and compositions are not known are labeled 1, 2 and 3.
30
40 kg/h 30 kg/h
0.9 kg A/kg 0.6 kg A/kg
0.1 kg B/kg 0.4 kg B/kg
100 kg/h m1 kg/h m2 kg/h m3 kg/h
0.5 kg A/kg x1 kg A/kg x2 kg A/kg x3 kg A/kg

0.5 kg B/kg 1 - x1 kg B/kg 1 – x2 kg B/kg 1 – x3 kg B/kg
30 kg/h
0.3 kg A/kg
0.7 kg B/kg
Solution:
Overall Mass Balance:
(100 + 30) kg/h = (40 + 30) kg/h + m3

m3 = 60 kg/h
Overall Mass Balance on A: (kg A/h)
(0.5)(100) + (0.3)(30) = (0.9)(40) + (0.6)(30) + x3 (60)

x3 = 0.0833 kg A/kg
Mass Balance on Unit 1:
100 = 40 + m1
m1 = 60 kg /h
Balance on A in Unit 1: (kg A/h)
(0.5)(100) = (0.9)(40) + x1 (60)

x3 = 0.233 kg A/kg
Mass Balance on Stream Mixing Point:
m1 + 30 = m2; m1 = 60 kg/h
m2 = 90 kg/h
31
Balance on A in Stream Mixing point: (kg A/h)
x1 m1 + (0.3)(30) = x2 m2; x1 = 0.233 kg A/kg

m1= 60 kg /h
m2= 90 kg /h
x2 = 0.255 kg A/kg
Do It Yourself:
A labeled flow chart of a continuous steady state double effect evaporator is shown below.
The system of evaporator is designed to reduce water from an incoming brine (NaCl +
H2O) stream from 60 wt % to 20 wt %. If the evaporator unit is receiving 12000 lb/hr of
NaCl (along with 60 wt % H2O), determine:
(a) The product rate (P2) of brine in lb/hr.

(Answer: P2 = 6000 kg/hr)
(b) The water removed from the brine in each evaporator (V1 V2) in lb/hr.
(Answer: V1= 3272.73 kg/hr; V2= 2727.27 kg/hr)
V1 V2
H2O 1.00 H2O 1.00
Brine
12000 lb/hr I II
NaCl 0.40
H2O 0.60
P1 P2
NaCl 0.55 NaCl 0.80
H2O 0.45 H2O 0.20
32
Problems:
1. Classify the following processes as (1) batch, (2) flow, (3) neither, or (4) both on
a time scale of one day:
(a) Oil storage tank at a refinery
(b) Flush tank on a toilet
(c) Catalytic converter on an automobile
(d) Gas furnace in a home
2. A continuous still is to be used to separate acetic acid, water, and benzene from
each other. On a trial run, the calculated data were as shown in the figure. Data
recording the benzene composition of the feed were not taken because of an
instrument defect. The problem is to calculate the benzene flow in the feed per
hour. How many independent material balance equations can be formulated for
this problem? How many variables whose values are unknown exist in the
problem?
3. A liquid adhesive, which is used to make laminated boards, consists of a polymer

dissolved in a solvent. The amount of polymer in the solution has to be carefully
controlled for this application. When the supplier of the adhesive receives an
order for 3000 kg of an adhesive solution containing 13 wt % polymer, all it has
on hand is (A) 500 kg of a 90 wt % solvent, (B) a very large quantity of a 20 wt %
solvent, and (C) pure solvent. Sketch the flow charts for the scaled processes and
calculate the weight of each of the three stocks that must be blended together to
fill the order.
4. An evaporator is fed continuously with 25,000 kg/h of a solution containing 10%

NaOH, 10% NaCl, and 80% H2O. During evaporation, water is boiled off, and
salt precipitates as crystals, which are settled and removed from the remaining
liquor. The concentrated liquor leaving the evaporator contains 50% NaOH, 2%
NaCl, and 48% H2O. Calculate:
a) The kilograms of water evaporated per hour
b) The kilograms of salt precipitated per hour
c) The kilograms of concentrated liquor produced per hour.
33
5. A waste solution of acid from titration process contains 33% of sulfuric acid
(H2SO4), 36% of nitric acid (HNO3) and 31% of water (H2O) by mass. 98% of
pure sulfuric acid and 78% of pure nitric acid have to be added through the waste
solution of acid to produce final mixed solution of acid which contains 40%
H2SO4 and 43% HNO3.
(a) Draw and completely label a flow chart of the process.
(b) Calculate the quantity of waste solution of acid and pure acid that’s
needed to be added together to produce 1500 kg of mixed solution of
acid.
6. A labeled flow chart of a continuous steady state triple effect evaporator is shown
below. The system of evaporator is designed to reduce water from an incoming
brine (NaCl + H2O) stream from 75 wt % to 3 wt %. If the evaporator unit is to
produce 14,670 lb/hr of NaCl (along with 3 wt % H2O), determine:
(a) The feed rate (F) of brine in lb/hr.
(b) The water removed from the brine in each evaporator (V1 V2 V3) in
lb/hr.
V1 V2 V3
H2O 1.00 H2O 1.00 H2O 1.00
Brine I II III
F lb/hr
NaCl 0.25
H2O 0.75
P1 P2 P3 14670 lb/hr
NaCl 0.33 NaCl 0.50 NaCl 0.97
H2O 0.67 H2O 0.50 H2O 0.03
34
CHAPTER 3
ENERGY AND ENERGY BALANCES
3.1 Forms of Energy
Three component of total energy of a system:
` Kinetic energy (Ek)

◦ energy due to the translational motion of the system as a whole relative to
some frame of reference (usually the earth’s surface) or to rotation of the
system about some axis.
` Potential energy (Ep)
◦ energy due to the position of the system in a potential field (such as a
gravitational or electromagnetic field).
` Internal energy (U)
◦ all energy possessed by a system other than kinetic and potential energy;
or
◦ Energy due to translation, rotation, vibration & electromagnetic
interactions of the molecules, atom and subatomic particle within the
system.
3.2 Transfer of Energy
` In closed system (i.e. no mass is transferred across the system boundaries while
the process is taking place), energy may be transferred between such a system and
its surroundings in two ways as heat or work.
Heat
◦ Energy that flows as a result of temperature difference between a system
and its surroundings.
◦ The direction of flow is always from a higher temperature to a low one.
◦ Heat is defined as positive when its transferred to the system from the
surroundings.
Work
◦ Energy that flows in response to any driving force other a temperature
difference, such as a force, a torque or a voltage.
◦ Work is defined as positive when it is done by the system on the
surroundings.
3.3 First Law of Thermodynamics
` Law of conservation of energy, which state that energy can neither be created nor
destroyed.
` General form of first law of thermodynamics
35
Inlet Energy + Heat - Outlet Energy – Work = Accumulation
where;
Inlet energy and outlet energy is summation/total of all energy such as potential,
kinetic and internal energy.
3.4 Kinetic Energy Equation (Ek)
` Kinetic energy, Ek (J) of an object of mass m (kg) moving with velocity u (m/s)
relative to the surface of the earth is
1
Ek = mu 2 [3.1]
2
` If the fluid enters a system with a mass flow rate m& (kg/s) and uniform velocity u
(m/s), the rate at which kinetic energy E& K (J/s) is transported into the system is
1 [3.2]
E& k = m& u 2
2
Example 3.1
Water flows into a 2 cm internal diameter (ID) pipe at a rate of 2.00 m³/h. Calculate E& k
for this stream in J/s.
Solution:
1
E& K = m& u 2
2
First, we calculate the velocity (u).
V&
Velocity (u) = Volumetric flowrate =u=
A(πr 2 )
Pipe cross sectional area
2.00m 3 1 100 2 cm 2 1h
= × × × = 1.77m / s
h π (1) cm
2 2 2
1 m 2
3600s
Then, we calculate the mass flow rate of the water (m& ).
2.00m 3 1000kg 1h
m& = V&ρ = × 3
× = 0.556kg / s
h m 3600s
36
Finally, calculate E& K
& 1 0.556kg / s (1.77) 2 m 2 1N

E K = m& u =
2
× 2
× = 0.870 N .m / s = 0.870 J / s
2 2 s 1kg.m / s
3.5 Potential Energy Equation (Ep)
` Gravitational potential energy, Ep of an object is given below
E p = mgz [3.3]
` If the fluid enters a system with a mass flow rate m& (kg/s) and an elevation
z relative to the potential energy reference plane.
E& p = m& gz [3.4]
` Normally we are interested in the change of potential energy during energy

balance calculation;
ΔE& p = E& p 2 − E& p1 = m& g ( z 2 − z1 ) [3.5]
Example 3.2
Crude pump oil is pumped at a rate of 15.0 kg/s from a point 220 meters below the
earth’s surface to a point 20 meters above the ground level. Calculate the attendant rate of
increase of potential energy.
Solution:
ΔE& p = E& p 2 − E& p1 = m& g ( z 2 − z1 )
20 m 15 kg 9.81m [20 − ( −220 ) ]m 1N
= × 2
× ×
Ground level s s 1kg .m / s 2
= 35300 N .m / s = 35300 J / s = 35300W
220 m
37
3.6 Energy Balances on Closed System
` Closed system is defined as no mass is transferred across the system boundaries

while the process is taking place.
` Energy balance for closed system;
Final System Energy – Initial System Energy = Net Energy Transferred to the System
where;
Initial energy system = U ( initial ) + EK (initial ) + EP ( initial )
Final energy system = U ( final ) + EK ( final ) + EP ( final )
Net energy transfer = Q −W
(U ( final ) − U ( initial ) ) + ( EK ( final ) − EK ( initial ) ) + ( EP ( final ) − EP ( initial ) ) = Q − W
or
ΔU + ΔEk + ΔE p = Q − W [3.6]
` When applying energy balance equation to a given process, the following point
must be aware;
1. The internal energy of a system depends almost entirely on the chemical
composition, state of aggregation (solid, liquid, or gas), and temperature of
the system materials. If no temperature changes, phase changes, or
chemical reactions occur in a closed system and if pressure changes are
less than a few atmospheres, then ∆ U ≈ 0.
2. If a system is not accelerating, then ∆ EK = 0. If a system is not rising or
falling, then ∆ EP = 0.
3. If a system and its surroundings are at the same temperature or the system
is perfectly insulated, then Q = 0. The process is then termed adiabatic.
4. Work done on or by a closed system is accomplished by movement of the
system boundary against a resisting force or the passage of an electrical
current or radiation across the system boundary. If there no moving parts
or electrical current at the system boundary, then W = 0.
3.7 Energy Balances on Open System

` In open system, mass is transferred across the system boundaries while the
process is taking place.
` Therefore work must be done on open system to push mass in and work is done
on the surrounding by mass that emerges from the systems.
` Both work terms must be include in the energy balance for open system.
38
` The net work done by an open system;
W& = W& s + W& fl
W& fl = W&out − W&in = PoutV&out − PinV&in
where;
W& S - shaft work (rate of work done by the process fluid on a moving part within
the system such as a pump rotor.)
W& fl - flow work (rate of work done by the fluid at the system outlet minus rate of
work done by the fluid at the system inlet.)
` ^ symbol is used to denote the specific property (property divided by mass or by

) )
mole) such as specific internal energy ( U kJ/kg), specific volume ( V m3/kg) and
so on.
` One important property for energy balance on open system is specific enthalpy
( Ĥ kJ/kg).
Hˆ = Uˆ + PVˆ [3.7]
` Sometimes, universal gas law constant can be used as a conversion factor to

evaluate specific enthalpy.
Example 3.3
The specific internal energy of helium at 300 K and 1 atm is 3800 J/mol, and the specific
molar volume at the same temperature and pressure is 24.63 L/mol. Calculate the specific
enthalpy of helium at this temperature and pressure, and the rate at which enthalpy is
transported by a stream of helium at 300 K and 1 atm with a molar flowrate of 250
kmol/h.
Solution:
Hˆ = Uˆ + PVˆ = 3800 J / mol + (1atm)(24.63L / mol ) = 3800 J / mol + 24.63L.atm / mol
To convert L.atm/mol into J/mol, we need the factor J/(L.atm). From the gas constant
table,
0.08206 L.atm/(mol.K) = 8.314 J/(mol.K)
By dividing these two, the desired factor will be obtained;
8.314 J /( mol.K )
= 101.3 J /( L.atm)
0.08206 L.atm /( mol.K )
39
So;
⎡ 24.63L.atm 101.3J ⎤
Hˆ = 3800 J / mol + ⎢ × = 6295 J / mol
⎣ mol L.atm ⎥⎦
If n = 250 kmol;
250kmol 6295 J 1000mol 1h

H& = nHˆ = × × × = 437150 J / s = 437kW
h mol 1kmol 3600s
Answer:
` Energy balances equation for open system;
ΔH& + ΔE& k + ΔE& p = Q& − W& s [3.8]
where;
ΔH& = ∑ m Hˆ
output
j j − ∑ m Hˆ
input
j j
stream stream
2 2
m ju j m ju j
ΔE& K = ∑
output 2
− ∑
input 2
stream stream
ΔE& P = ∑ m gz
output
j j − ∑ m gz
input
j j
stream stream
3.8 Reference States and State Properties
` It is not possible to know the absolute value of Û and Ĥ for a process material,
but we can determine the change in Δ Û and change in Δ Ĥ corresponding to a
specific change of state (temperature, pressure, phase).
` A convenient way to tabulate Δ Û and Δ Ĥ is to choose a temperature, pressure
and state of aggregation (i.e. phase) as a reference state.
` Since Ĥ cannot be known absolute, for convenience we may assign a value
Hˆ O = 0 to be a reference state. Then ΔHˆ 1 = Hˆ 1 − 0; ΔHˆ 2 = Hˆ 2 − 0; and so on.
` Some enthalpy tables give the reference states on which the listed values of Ĥ are
based and others do not.
` However, we do not have to know the reference state to calculate Δ Ĥ for the
transition from one state to another state since the values are based on the same
reference in the table.
40
` Ĥ and Û can be said as state properties (property whose change of value in any
process depend only on it initial and final states and do not depend on the path
take to reach the state.
Example 3.4
The following data are for methyl chloride:
State T(°F) P(psia) Vˆ ( ft 3 / lbm ) Hˆ ( Btu / lbm )

Liquid -40 6.878 0.01553 0.0
Vapor 0 18.9 4.969 196.23
Vapor 50 51.99 1.920 202.28
(a) What reference state was used to generate the given enthalpies?
(b) Calculate ΔĤ and Δ̂U for the transition of saturated methyl chloride vapor from
50°F to 0°F.
Solution:
(a) Liquid at -40 °F and 6.878 psia (the state at which Hˆ = 0 ).
(b) ΔHˆ = Hˆ final − Hˆ initial = (196.23 − 202.28) Btu / lbm = −6.05Btu / lbm
From equation 3.7,
ΔHˆ = ΔUˆ + ΔPVˆ

⎡ (93.914 − 99.821) psia. ft 3 / 1.987 Btu ⎤
∴ ΔU = ΔH − Δ(PV ) ) = −6.05Btu / lbm − ⎢ × ⎥
⎣ lbm 10.73 ft 3 . psia ⎦
ΔU = −4.96 Btu / lbm
3.9 Steam Table
The simplest phase diagrams are pressure-temperature diagrams of a single simple

substance, such as water. Figure 3.1 shows the water phase diagram. The axes correspond
to the pressure and temperature. The phase diagram shows, in pressure-temperature space,
41
the lines of equilibrium or phase boundaries between the three phases of solid, liquid,
and gas.
The phase boundary between liquid and gas does not continue indefinitely. Instead, it
terminates at a point on the phase diagram called the critical point. This reflects the fact
that, at extremely high temperatures and pressures, the liquid and gaseous phases become
indistinguishable, in what is known as a supercritical fluid. In water, the critical point
occurs at around Tc=647.096 K (1,164.773 °R), pc=22.064 MPa (3,200.1 psi) and
ρc=356 kg/m³.
Figure 3.1: Water Phase Diagram
Appendix 1 lists properties of saturated liquid water and saturated steam at temperature
from 0.01 °C to 100 °C. The following properties can be determined for each tabulated
temperature:
Column 2. The pressure
Column 3. The specific volumes, vˆ g / v g (m³/kg), of gas or saturated steam.
Column 4 – 5. The specific internal energies; u f for liquid/fluid, u g for gas.
42
Column 6-8. The specific enthalpies; h f for liquid/fluid, h fg for intermediate (liquid and
gas), hg for gas.
Column 9 – 11. The entropies. s f for liquid/fluid, s fg for intermediate (liquid and gas),
s g for gas.
Page 2 (Appendix 1) list the same properties as page 1(Appendix 1), except the pressure
is the first column and temperature the second and the table covers a much broader range
of temperatures and pressures. Appendix 2 is known as superheated steam table. It is used
for superheated properties of steam.
Example 3.5
(a) Determine the pressure, specific internal energy and specific enthalpy of saturated
steam at 330.8 °C.
(b) Show that water at 400 °C and 10 bar is superheated steam and determine its specific
volume, specific internal energy and specific enthalpy.
Solution:
(a) Data in page 1 (Appendix 1) does not go to 330.8 °C, so we take a look at Page 2.
p = 130bar; u g = 2496kJ/kg; hg = 2662kJ/kg
(b) The data in Appendix 1 does not cover temperature more than 374.15, therefore it is
superheated steam.
v = 0.3065; u g = 2957kJ/kg; hg = 3262kJ/kg
3.9 Energy Balance Tips
` When labeling flowchart, write down together the temperature, pressure and state
of aggregation of the process material.
` Normally (depend on the process description) for chemical process unit; shaft
work, kinetic and potential energy change tend to be negligible compared to heat
flows, internal energy and enthalpy changes.
` Then simplified energy balance become;
For closed system
Q = ΔU [3.9]
For open system
Q& = ΔH& [3.10]
43
Problems
1. Liquid methanol is pumped from a large storage tank through a 1 inch internal
pipe (ID) at the rate of 3.00 gal/min. At what rate in ft.lbf/s and hp is kinetic
energy being transported by the in the pipe methanol? (Given: Density of
methanol = 49.5 lbm/ft3).
2. Air at 300 °C and 130 kPa flows through a horizontal 7 cm ID pipe at velocity of
42.0 m/s. Calculate E& K (W ) , assuming ideal gas behavior.
3. If you pour 1 gallon of water on a yowling cat 10 ft below your bedroom window,
how much potential energy (ft.lbf) does the water lose?
4. Carbon monoxide (CO) at 120 K and 45 atm has a tabulated specific volume of
5.23 cm3/g and specific internal energy of 1813 J/mol. Calculate the specific
enthalpy of CO in this state.
5. Oxygen at 150 K and 41.64 atm has a tabulated specific volume of 4.684 cm³/g
and a specific internal energy of 1706 J/mol. Calculate the specific enthalpy of O2
in this state.
6. Values of the specific internal energy of a fuel gas at two conditions are listed
below.
State T(K) P(bar) Vˆ (L/mol) Û (kJ/mol)

[Fasa]
Liquid 320 0.505 0.0978 0.000
[Cecair]
Vapor 320 0.550 97.78 35.4
[Wap]
i) What reference state was used to generate the listed specific internal
energies?
ii) Calculate ΔÛ (kJ/mol) for a process in which a fuel gas vapor at 320 K is
condensed at constant pressure. Then, calculate ΔĤ (kJ/mol) for the same
process. Finally, calculate ΔĤ (kJ) for 25 mol of the fuel gas that undergo
the process.
44
7. Complete the following table.
T(OC) P(Bar) Vˆ (m3/kg) Û (kJ/kg) Phase

[Fasa]
88 ........... ................. ................ Saturated steam
.......... ........... ................. 1511 Water
.......... 32 ................. ................ Water
475 15 ................. ................ ........................................
.......... 197 ................. ................ Saturated steam
45
CHAPTER 4
ENERGY BALANCE ON NONREACTIVE SPECIES
4.1 Introduction
` Normally in chemical process unit, WS =0; Δ EP =0; Δ EK =0; Then energy balance
equation become:
Close System Open System
Q = ΔU Q = ΔH
` For this chapter, we will learn the procedure for evaluating ΔU and ΔH when
table Ĥ and Û are not available for all process species.
` Method to calculate Δ Ĥ and Δ Û associated with certain process such as:
1. Change in P, at constant T & constant state of aggregation.
2. Change in T, at constant T & constant state of aggregation.
3. Phase changes at constant T & constant P.
4. Mixing at constant T & constant P.
5. Chemical reaction at constant T & constant P.
4.2 Hypothetical Process Path
` State properties
o Properties that depend on the state of the species (primarily on its
temperature and state of aggregation, and to lesser extent on its pressure).
o Specific enthalpy ( Ĥ ) and specific internal energy ( Û ) are state
properties species.
o When a species passes from one state to another state, both Δ Ĥ and Δ Û
for the process are independent of the path taken from the first state to the
second state.
` We can construct a hypothetical process path which can consist of several step
based on our convenience, as long as we reach to the final state starting from their
initial state.
` For example, the enthalpy change (Δ Ĥ ) of solid phenol at 25 oC and 1 atm to
phenol vapor at 300 oC and 3 atm.
ΔĤ phenol = (vapor, 300˚C, 3 atm) – (solid, 25˚C, 1 atm)
o Cannot determine directly from enthalpy table – must use hypothetical
process path consist of several step.
o Check Table of Physical Properties Data Appendix 5: P= 1 atm; Tm=
42.5°C and Tb= 181.4°C
46
ΔHˆ = ΔHˆ 1 + ΔHˆ 2 + ΔHˆ 3 + ΔHˆ 4 + ΔHˆ 5 + ΔHˆ 6
4.3 Procedure Energy Balance Calculations
1. Perform all required material balance calculations.

2. Write the appropriate form of the energy balance (closed or open system) and
delete any of the terms that are either zero or negligible for the given process
system.
3. Choose a reference state – phase, temperature, and pressure – for each species
involved in the process.
4. Construct inlet-outlet table for specific internal energy (close system) or specific
enthalpy (close system)
{ For closed system, construct a Table with columns for initial and final
amounts of each species (mi or ni) and specific internal energies ( Û )
relative to the chosen reference states.
{ For an open system, construct a table with columns for inlet and outlet
stream component flow rates (mi or ni) and specific enthalpies ( Ĥ )
relative to the chosen references states.
5. Calculate all required values of Ĥ or Û and insert the values in the appropriate
places in the table. Then calculate Δ Ĥ or Δ Û for the system.
6. Calculate any work, kinetic energy, or potential energy terms that you have not
dropped from the energy balance.
7. Solve the energy balance for whichever variable is unknown (often Q).
47
Example of Inlet-Outlet Enthalpy Table:
References: Ac (liquid, 20˚C, 5atm); N2 (gas, 25˚C, 1atm)

Substance Inlet Outlet
n&in Ĥ in n&out Ĥ out
Ac (vapor) 66.9 3.35
Ac (l) - - 63.55 0
N2 33.1 33.1
4.4 Change in Pressure (P) at Constant Temperature (T) & Constant Phase
` Solid & Liquid

o Nearly independent of pressure.
ΔUˆ = 0
ΔHˆ = VˆΔP
` Ideal Gases
o Independent of pressure (unless undergo very large pressure changes).
Δ Uˆ = 0
Δ Hˆ = 0
4.5 Change in Temperature (T) at Constant Pressure (P) & Constant Phase
` Sensible heat is the heat that must be transferred to RAISE or LOWER the
temperature of substance or mixture of substance.
o C P - Heat capacity at constant pressure (given in Appendix 3 the form of
polynomial and the equation is the function of temperature).
o CV - Heat capacity at constant volume.
C p = Cv Liquid & Solid
C p = Cv + R Ideal Gas
` Specific internal energy change:

T2 Ideal gas : exact
ΔUˆ = Cv (T )dT Solid or liquid
∫ : good approximation
Nonideal gas : valid only if volume (V) is
T1
constant
48
` Specific enthalpy change:
T2
ΔHˆ = ∫ C p (T )dT Ideal gas : exact
Nonideal gas : exact only if pressure (P) is
T1
constant
T2 Solid & liquid

ΔHˆ = VˆΔP + ∫ C p (T )dT
T1
4.6 Heat Capacities, Cp
` Estimation of heat capacities, Cp
o Kopp’s rule is the simple empirical method for estimating Cp of solid or

liquid at 20 °C based on the summation of atomic heat capacities of the
molecular compound.
(CP )Ca ( OH ) 2 = (CPa )Ca + 2(CPa )O + 2(CPa ) H
= 26 + (2 × 17) + (2 × 9.6) = 79 J / mol.O C
` Estimation for heat capacities of mixtures
(C p ) mix (T ) = ∑ yi C pi (T )
CPi = Cp for ith component
yi = mass or moles fraction
Example 4.1
Calculate the heat required to raise 200 kg nitrous oxide (N2O) from 20 °C to 150 °C in a
constant volume vessel. The constant volume heat capacity of N2O in this temperature
range is given by this equation;
Cv (kJ / kg.o C ) = 0.855 + 9.42 ×10 −4 T
where T is °C.
Solution:
150o C
∫ (0.855 + 9.42 ×10 T )dT

−4
ΔU (kJ / kg ) =
20o C
150 o C
⎡ 9.42 × 10 − 4 T 2 ⎤
= [0.855T ]20 o C
o
150 C
+⎢ ⎥
⎣ 2 ⎦ 20 o C
⎡ 9.42 × 10 − 4 (1502 − 30 2 ) ⎤
= [0.855(150 − 20)] + ⎢ ⎥
⎣ 2 ⎦
= 111.15 + 10.491
= 121.56kJ / kg
49
Example 4.2
15 kmol/min of air is cooled from 430 °C to 100 °C. Calculate the required heat removal
rate using 1) heat capacities formulas from Heat Capacities Table (Appendix 3) and 2)
Specific Enthalpies Table (Appendix 4).
Solution:
1. Write the energy balance for open system:
0 0 0
ΔH& + ΔE& K + ΔE& p = Q& + Ws
Q& = ΔH& = n&air Hˆ air , out − n&air Hˆ air ,in = n& Hˆ
2. Using the heat capacities table (Appendix 3):

100o C
ΔHˆ (kJ / mol ) = ∫C p (T )dT
o
430 C
100o C
ΔHˆ (kJ / mol ) = ∫ [28.94 ×10 ]
−3
+ 0.4147 × 10 −5 T + 0.3191× 10 −8 T 2 − 1.965 × 10 −12 T 3 dT
o
430 C o
100 C o
100 C 100o C
⎡ 0.4147 × 10−5 T 2 ⎤ ⎡ 0.3191× 10−8 T 3 ⎤ ⎡1.965 ×10 −12 T 4 ⎤
[
ΔHˆ (kJ / mol ) = 28.94 × 10 −3 T 430oC + ⎢ ]
100o C
2
⎥ + ⎢
3
⎥ − ⎢
4
⎥
⎣ ⎦ 430oC ⎣ ⎦ 430oC ⎣ ⎦ 430oC
⎡ ⎤ ⎡ )⎤
[ ]
−5 −8
0 .4147 × 10 (100 2
− 430 2
) 0 .3191 × 10 (100 3
− 430 3
ΔHˆ ( kJ / mol ) = 28.94 × 10 − 3 T (100 − 430) + ⎢ ⎥ + ⎢ ⎥
⎣ 2 ⎦ ⎣ 3 ⎦
⎡1.965 × 10 −12 (100 4 − 430 4 ) ⎤
−⎢ ⎥
⎣ 4 ⎦
ΔHˆ (kJ / mol ) = −9.98kJ / mol
3. Using Tabulated Enthalpies (Appendix 4):
Read the value directly from the table according to the temperature desired:
Hˆ (100o C ) = 2.19kJ / mol

Hˆ (430o C ) = 12.179kJ / mol
ΔHˆ = (2.19 − 12.179)kJ / mol = −9.98kJ / mol
4. Calculate the heat removal rate.
15kmol 1000mol 1 min − 9.98kJ 1kW

Q = ΔH& = n& ΔHˆ = × × × × = 2495kW
min 1kmol 60 s mol 1kJ / s
50
Note: We can use Table in Appendix 4 if the gases are covered in the table. If not,
you have to calculate using the harder way that is heat capacities (Appendix 3).
Example 4.3
Calculate the heat required to bring 150 mol/h of a stream containing 50% C2H6 and 50%
C3H8 by volume from 0 OC to 400 OC. Determine a heat capacity for the mixture.
Solution:
From heat capacities table;

CP (C 2 H 6 ) = 49.37 × 10 −3 + 13.92 × 10 −5 T − 5.816 × 10−8 T 2 + 7.28 × 10−12 T 3
CP (C3 H 8 ) = 68.032 × 10 −3 + 22.59 × 10 −5 T − 13.11 × 10 −8 T 2 + 31.71 × 10−12 T 3
400o C
∫ 49.37 ×10
−3
(Cp) mix [kJ /( mol. C )] = 0.5
o
+ 13.92 × 10 −5 T − 5.816 × 10 −8 T 2 + 7.28 ×10 −12T 3
0o C
400o C
∫ 68.032 ×10
−3
+ 0.5 + 22.59 × 10 −5 T − 13.11× 10 −8 T 2 + 31.71×10 −12T 3
0o C
Try to integrate and calculate by yourself.
ΔHˆ = 36.19kJ / mol

150mol 36.19kJ 1h
Q& = ΔH& = n&ΔHˆ = × × = 1.51kW
h mol 3600s
Example 4.4
A stream of gas containing 10 % CH4 and 90 % air by volume is to be heated from 20 °C

to 300 °C. Calculate the required rate of heat input in kW if the flow rate of the gas is
2.00 x 103 liters (STP)/min.
Solution:
1. Draw the flow chart.
2. Change the flow rate in STP to mol:
2000 L( STP ) 1mol

n& = × = 89.3mol
min 22.4 L( STP )
CH4= 0.1(89.3 mol) =8.93 mol; Air = (89.3-8.93)mol =80.4mol
51
3. Write the energy balance;
0 0 0
Δ H& + Δ E& K + Δ E& p = Q& + W s
Q& = Δ H&
4. Construct inlet-outlet table.
References: CH4 (gas, 20˚C, 1atm); air (gas, 25˚C, 1atm) Why we choose air at 25°C
and 1 atm as our reference?
Substance Inlet Outlet This is because it is the
reference state for gases in
Table B.8, so we can still
n&in Ĥ in n&out Ĥ out calculate the enthalpy at
20°C.
CH4 8.93 0 8.93 Ĥ1
air 80.4 Ĥ 2 80.4 Ĥ 3
5. Calculate the enthalpies.

300o C
Hˆ 1 = ∫ (C )
p CH 4 dT
20o C
300 o C
∫ (34.31×10
−3
= + 5.469 ×10 −5 T +0.3661×10 −8 T 2 − 11.0 × 10 −12 T 3 )dT
20o C
= 12.1kJ / mol
To calculate H2 and H3, use Table in Appendix 4:

Why we calculate ΔĤ of air
Hˆ o
( air , 25 C )
= 0.0kJ / mol Hˆ o
( air , 25 C )
= 0.0kJ / mol using Table B.8 instead of
Table B.2?
It is because the air is
Hˆ ( air , 20 o C )
= −0.144kJ / mol Hˆ ( air , 300o C )
= 8.17 kJ / mol included in Table B.8. We can
still integrate and calculate
ΔHˆ 2 = (−0.144 − 0)kJ / mol ΔHˆ 2 = (8.17 − 0)kJ / mol
using the long way. You will
get the same answer.
Insert the value of enthalpies inside the inlet-outlet table.
References: CH4 (gas, 20˚C, 1atm); air (gas, 25˚C, 1atm)

CH4 8.93 0 8.93 12.1
air 80.4 -0.144 80.4 8.17
52
6. Calculate the heat.
Q& = ΔH& = ∑ n&i Hˆ i − ∑ n&i Hˆ i

out in
= [(8.93mol / min)(12.1kJ / mol ) + (80.4mol / min)(8.17 k ) / mol ]

− [(80.4mol / min)(−0.144kJ / mol ) + (8.93mol / min)(0kJ / mol )]
kJ 1 min
= 776.5 × = 12.94kW
min 60s
4.7 Phase Change Operations
Phase change such as melting and evaporation are usually accompanied by large changes
in internal energy and enthalpy.
` Latent heat
o Specific enthalpy change associated with the phase at constant
temperature and pressure.
` Heat of fusion or heat of melting, ΔĤm (T,P)

o Specific enthalpy different between solid and liquid forms of
species at T & P.
o Heat of solidification (liquid to solid) is negative value of heat of
fusion.
` Heat of vaporization, ΔĤv (T,P)

o Specific enthalpy different between liquid and vapor forms of
species at T & P
o Heat of condensation (vapor to liquid) is negative value of heat of
vaporization.
The latent heat of phase change may vary considerably with the temperature at which the
changes occurs but hardly varies with the pressure at the transition point.
4.8 Estimation of Heat of Vaporization
1. Trouton’s rule – accuracy between 30%.

ΔHˆ v (kJ / mol ) = 0.088Tb nonpolar liquid
ΔHˆ (kJ / mol ) = 0.109T
v b water or low MW alcohol
where;
Tb = Boiling point temperature
2. Chen’s equation – accuracy between 2%.

Tb [0.0331(Tb / Tc ) − 0.0327 + 0.0297 log10 Pc ]
ΔHˆ v (kJ / mol ) = 53
1.07 − (Tb / Tc )
where;
Tb = Normal boiling point temperature
Tc = Critical temperature
Pc = Critical Pressure
3. Clausius-Clapeyron equation - plot ln p* versus 1/T.
ΔHˆ v
In p * = − +B
RT
4. Chaperon equation
d ( In p* ) ΔHˆ v
=−
d (1 / T ) R
5. Watson correlation – estimate ΔĤv at T2 from known ΔĤv at T1.

0.38
⎛T −T ⎞
ΔHˆ v (T2 ) = ΔHˆ v (T1 )⎜⎜ c 2 ⎟⎟
⎝ Tc − T1 ⎠
where;
Tc = Critical temperature
4.8 Estimation of Heat of Fusion
ΔĤm (kJ/mol) = 0.0092 Tm (K) metallic elements

= 0.0025 Tm (K) inorganic compound
= 0.050 Tm (K) organic compound
Example 4.5
The normal boiling point of methanol is 337.9 K, and the critical temperature of this
substance is 513.2 K. Estimate the heat of vaporization of methanol at 300 °C .
Solution:
First, we calculate the standard heat of vaporization using Trouton’s rule:

ΔHˆ (kJ / mol ) ≈ 0.109T ( K )
V b
≈ (0.109(337.9))
≈ 36.83kJ / mol
Then, we calculate the latent heat using Watson’s correlation:
0.38 0.38
⎛T −T ⎞ ⎛ 513.2 − 473 ⎞
ΔHˆ V (T2 ) = ΔHˆ V (T1 )⎜⎜ c 2 ⎟⎟ = 36.83⎜ ⎟ = 21.0kJ / mol
⎝ Tc − T1 ⎠ ⎝ 513.2 − 337.9 ⎠
54
Problems:
1. Estimate the specific enthalpy of steam (kJ/kg) at 350 °C and 100 bar relative to
steam at 100 °C and 1 atm using:
(a) The steam tables
(b) Heat capacities table
2. Chlorine gas is to be heated from 100 °C and 1 atm to 200 °C.

(a) Calculate the heat input (kW) required to heat a stream of the gas flowing
at 5.0 kmol/s at constant pressure.
(b) Calculate the hat input (kJ) required to raise the temperature of 5.0 kmol
of chlorine in a closed rigid vessel from 100 °C and 1 atm to 200 °C.
3. A stream of water vapor flowing at a rate of 250 mol/h is brought from 600 °C
and 10 bar to 100 °C and 1 atm.
(a) Estimate the required cooling rate (kW) three ways: (i) from the steam
tables, (ii) using heat capacity data and (iii) using specific enthalpy data.
(b) Which of the answers in part (a) is most accurate and why?
4. A fuel gas containing 95 mole% methane and the balance ethane is burned
completely with 25% excess air. The stack gas leaves the furnace at 900 °C and is
cooled to 450 °C in a waste heat boiler, a heat exchanger in which heat lost by
cooling gases is used to produce steam from liquid water for heating , power
generation, or process applications.
(a) Taking as abasis of calculation 100 mol of the fuel gas fed to the furnace,
calculate the amount of heat (kJ) that must be transferred from the gas in
the waste heat boilerto accomplish the indicated cooling.
(b) How much saturated steam at 50 bar can be produced from boiler
feedwater at 40 °C for the same basis of calculation?
55
CHAPTER 5
ENERGY BALANCE OF REACTIVE SYSTEMS
5.1 Types of reaction
` Exothermic reaction: the product molecules have lower internal energies than the
reactants at the same T and P. ΔH = NEGATIVE, reaction produces heat.
` Endothermic reaction: the product molecules have higher internal energies than
the reactants. ΔH = POSITIVE, reaction consumes heat.
5.2 Heat of reaction
` ΔH depends on how the stoichiometric equation is written.
CH4 (g) + 2O2 (g) Î CO2(g) + 2H2O(l) ΔHr1 (25OC) = -890.3 kJ/mol Î for 1 CH4
2CH4 (g) + 4O2 (g) Î 2CO2 (g) + 4H2O (l) ΔHr2 (25OC) = -1780.6 kJ/mol Î for 2 CH4
` ΔH depends on the states of aggregation (gas, liquid, or solid)
CH4 (g) + 2O2 (g) Î CO2 (g) + 2H2O (l) ΔHr1 (25°C)= -890.3 kJ/mol
CH4 (g) + 2O2 (g) Î CO2 (g) + 2H2O (g) ΔHr2 (25°C)= -802.3 kJ/mol
5.3 Standard heat of reaction (ΔHro)
` heat of reaction when both reactants and products are at reference conditions
(usually 25 °C and 1 atm).
C4H10 (g) + 13/2O2 (g) Î 4CO2 (g) + 5H2O(l) ΔHr1 (25 °C)= -2878 kJ/mol
Example: For 2400 mol/s CO2 produced;
2400 ⎛ kJ ⎞
ΔH r 2 = × ⎜ − 2878 ⎟ = −1.73 × 10 kJ / s
6
4mol / s ⎝ mol ⎠
5.4 Reaction in a REACTOR (constant V)
ΔUr(T)=Uproducts – Ureactants
⎛ ⎞
⎜ ⎟
ΔU r (T ) = ΔH r (T ) − RT ⎜ ∑ vi − ∑ vi ⎟
⎜ gaseous gaseous ⎟
⎝ products reactants ⎠
56
5.5 Hess’s Law
` Look at this reaction.
C (s) + ½ O2 (g) Î CO (g): ΔHr = ?

(1) C (s) + O2 (g) Î CO2 (g): ΔHr1 = -393.51 kJ/mol
(2) CO (g) + ½ O2 (g) Î CO2 (g): ΔHr2 = -282.99 kJ/mol
ΔHr = ΔHr1 + (- ΔHr2) = (-393.51 + (-282.99)) = -110.52 kJ/mol
Example 5.1
The standard heats of the following combustion reactions have been determined
experimentally:
1. C2H6 + 7/2 O2 Î 2 CO +3 H2O: ΔHˆ r01 = −1559.8kJ / mol
2. C + O2 Î CO2: ΔHˆ r02 = −393.5kJ / mol
3. H2 + ½ O2 Î H2O: ΔHˆ r03 = −285.8kJ / mol
Use Hess’s law and the given heats of reaction to determine the standard heat of reaction
4. 2 C + 3 H2 Î C2H6
ΔHˆ r 4 = ?
0
Solution:
1) Rearrange the given equations:
(Eq. 2 x 2) 2C + 2O2 Î 2CO2 ΔHˆ r02 = (−393.5 × 2)kJ / mol = −787kJ / mol
(Eq. 3 x 3) 3H2 + 3/2 O2 Î 3H2O ΔHˆ r03 = (−285.8 × 3)kJ / mol = −857.4kJ / mol
(Eq. 1) C2H6 + 7/2 O2 Î 2 CO +3 H2O ΔHˆ r03 = −1559.8kJ / mol
2) Substrate or add these equations to form the desired equation:

3H2 + 3/2 O2 Î 3H2O ΔHˆ r03 = −857.4kJ / mol
+ 2C + 2O2 Î2CO2 ΔHˆ r02 = −787 kJ / mol
3H2 + 7/2 O2 + 2C Î 3H2O + 2CO2 ΔHˆ r0( 3− 2) = −1644.4kJ / mol

- C2H6 + 7/2 O2 Î 2 CO2 +3 H2O ΔHˆ r01 = −1559.8kJ / mol
3H2 + 2C ÎC2H6 ΔHˆ r04 = −84.6kJ / mol
57
5.6 Heat of combustion (ΔHc)
` The heat of the combustion of a substance with oxygen to yield specific products.
C2H5OH (l) + 3 O2 (g) Î 2 CO2 (g) + 3 H2O (l): ΔHc (25 oC, 1 atm) = -1366.9 kJ/mol
` The value of Hc for the substance is indicated in Table B.1.
⎛ ⎞
ΔH r = −∑ vi (ΔH co ) i =⎜⎜ ∑ vi (ΔH co ) i − ∑ vi (ΔH co ) i ⎟⎟
o
i ⎝ reactants products ⎠
Example 5.2:
Calculate the heat of reaction for the dehydrogenation of ethane:
C2H6 Î C2H4 + H2
Then, use Hess’s law to find the standard heat of reaction for the above equation, using
these reactions:
1. C2H6 + 7/2O2 Î2CO2 + 3H2O

2. C2H4 + 3 O2 Î2CO2 + 2H2O
3. H2 + ½ O2 ÎH2O
Compare the answer.
Solution:
1) Find standard heat of combustion for every substance. From Appendix 5;
(ΔH c0 ) C2 H 6 = −1559.9kJ / mol

(ΔH c0 ) C2 H 4 = −1410.99kJ / mol
(ΔH c0 ) H 2 = −285.84kJ / mol
2) Use the formula to calculate the standard heat of reaction:

⎛ ⎞
ΔH r = −∑ vi (ΔH co ) i =⎜⎜ ∑ vi (ΔH co ) i − ∑ vi (ΔH co ) i ⎟⎟ = [(−1559.9) − (−1410.99 − 285.84)kJ / mol )
o
i ⎝ reactants products ⎠
o
ΔH r = −136.9kJ / mol
58
3) Find the standard heat of reaction for each equation using standard heat of
combustion:
ΔH r0 ( Eq.1) = −1559.9kJ / mol

ΔH r0 ( Eq.2 ) = −1410.99kJ / mol
ΔH r0 ( Eq.3) = −285.84kJ / mol
4) Substrate or add these equations to form the desired equation:
(1) – (2) – (3) ΔH r0 ( Eq.1) = −1559.9kJ / mol

C2H6 + 7/2O2 Î2CO2 + 3H2O
ΔH r0 ( Eq.2 ) = −1410.99kJ / mol
- C2H4 + 3O2Î2 CO2 + 2 H2O
C2H6 + ½ O2ÎC2H4 + H2O ΔH r0 ( Eq.1−2 ) = −148.91kJ / mol
- H2 + ½ O2 ÎH2O ΔH r0 ( Eq.3) = −285.84kJ / mol

C2H6 ÎC2H4 + H2 ΔH r0 = −136.93kJ / mol
5.7 Heat of formation (ΔHf)
` Enthalpy change associated with the formation of 1 mole of a compound from its
elemental constituents (in nature) at a reference T and P.
N2 (g) + 2 H2 (g) + 3/2 O2 (g) Î NH4NO3 (c) ΔHro = -365.14 kJ/mol

6 C(s) + 3 H2 (g) Î C6H6 (l) ΔHro = 48.66 kJ/mol
Note: standard heat of formation of an elemental species is ZERO.
⎛ ⎞
ΔH r = ∑ vi ΔH ofi =⎜⎜ ∑ vi ΔH ofi − ∑ vi ΔH ofi ⎟⎟
o
i ⎝ products reactants ⎠
Example 5.3
Determine the standard heat of reaction for the combustion of liquid n-pentane, assuming
H2O(l) is a combustion product.
C5H12 (l) + 8O2 (g) Î5CO2 (g) + 6H2O (l)
Then, use Hess’s law to find the standard heat of reaction for the above equation, using
these reactions:
1. 5C (s) + 6H2 ÎC5H12
2. C (s) + O2 (g) ÎCO2
3. H2 (g) + ½O2 (g) ÎH2O (l)
Compare the answer.
59
Solution:
1) Find standard heat of formation for every substance. From table B.1;
(ΔH 0f )C 5 H 12 = −173.0kJ / mol
(ΔH 0f )CO2 = −393.5kJ / mol
(ΔH 0f ) H 2 O = −285.84kJ / mol
2) Use the formula to calculate the standard heat of reaction:
⎛ ⎞
ΔH r = −∑ vi (ΔH of )i =⎜⎜ ∑ vi (ΔH of )i − ∑ vi (ΔH of )i ⎟⎟ = [(5 × −393.5) + (6 × −285.84) − (−173.0)kJ / mol )
o
i ⎝ products reactants ⎠
o
ΔH r = −3509kJ / mol
3) Find the standard heat of reaction for each equation using standard heat of
formation and use Hess Law to calculate the desired equation:
1. 5C (s) + 6H2 ÎC5H12

2. C (s) + O2 (g) ÎCO2
3. H2 (g) + ½O2 (g) ÎH2O (l)
o
The desired reaction can be obtained by; 5 × (2) + 6 × (3) − (1) , and the ΔH r can be
obtained from Hess Law.
5.8 Energy Balances (General Procedures)
` Heat of reaction method:
60
` Heat of formation method:
Some notes (extent of reaction)

` For a single reaction, extent of reaction (ξ) can be calculated:
n A−out − n A−in
ξ=
vA
where A is any reactant or product known
Some notes (inlet-outlet enthalpy table)
Components nin Hin nout Hout

A
B
C
Some notes
` Latent heat Î heat transferred without change of T. It could be as heat of

vaporization (or condensation). There is a phase change.
ΔHˆ = ΔHˆ v
` Sensible heat Î heat transferred due to the T difference. There is a change of T,
but no phase change.
T2
ΔHˆ = ∫ C p dT
T
Since Cp = 1f (T), then don’t forget to integrate it. See Appendix 3.
61
STEP-by-STEP procedure of calculation:
y Material balance calculation of reactor.
y Choose reference states (usually 25 °C, 1 atm).
y Calculate extent of reaction.
y Prepare inlet-outlet enthalpy table.
y Calculate unknown component enthalpy.
y Calculate ΔH for the reactor.
ΔH = ξΔHˆ ro + ∑ nout Hˆ out −∑ nin Hˆ in (single reaction)

ΔH = ∑ ξ ΔHˆ
reaction
j
o
rj + ∑ nout Hˆ out −∑ nin Hˆ in (multiple reactions)
then
Q − W = ΔH + ΔEk + ΔE p can be evaluated
Example 5.4
Normal heptane is dehydrocyclized to toluene in a continuous vapor-phase reaction:
C7 H16 ÎC6H5CH3 + 4H2
Pure heptane at 400 °C is fed to the reactor. The reactor operates isothermally at 400 °C
and the reaction goes to completion.
(a) Taking basis of 1 mol of heptane fed, draw and label a flowchart.
(b) Taking elemental species [C(s), H2(g)] at 25 °C as references, prepare and
fill in an inlet-outlet enthalpy table.
(c) Calculate the required heat transfer.
Given: C7H16(g)= 137.44x10-3 + 40.85x10-5 T – 23.92x10-8 T2 + 57.66x10-12 T3.
Solution:
1) Draw the flowchart.
Basis: 1 mol heptane fed
62
2) Prepare inlet-outlet Table
References: C(s), H (g) at 25 °C

C7H16 1 Ĥ1 - -
C7H8 - - 1 Ĥ 2
H2 - - 4 Ĥ 3
3) Find all the enthalpies

400o C
Hˆ 1 = Hˆ 0f ( C H ) + ∫C p dT
7 16
o
25 C
400o C
∫ 137.44 ×10
−3
= −187.8kJ / mol + + 40.85 × 10 −5 T − 23.92 ×10 −8 T 2 + 57.66 ×10 −12 T 3 dT
25o C
= (−187.8 + 79.34)kJ / mol = −108.46kJ / mol

400o C
Hˆ 2 = Hˆ 0
f (C H ) + ∫ C dT p
7 8
25o C
400o C
∫ 94.18 ×10
−3
= +50.00kJ / mol + + 38.00 × 10 −5 T − 27.86 × 10 −8 T 2 + 80.33 ×10 −12 T 3 dT
25o C
= (+50.00 + 60.17)kJ / mol = 110.17 kJ / mol

H3 can be obtained directly from Appendix 4:
Hˆ 3 = Hˆ ( H o − Hˆ ( H o = (10.89 − 0)kJ / mol = 10.89kJ / mol

2 , 400 C) 2 , 25 C)
63
3) Insert the H values calculated inside the inlet-outlet table.
References: C(s), H (g) at 25 °C

C7H16 1 Ĥ1 = - - -
108.46
C7H8 - - 1 Ĥ 2 = 110.17
H2 - - 4 Ĥ 3 = 10.89
5) Calculate the standard heat of reaction.

⎛ ⎞
ΔH r = ∑ vi ΔH ofi =⎜⎜ ∑ vi ΔH ofi − ∑ vi ΔH ofi ⎟⎟
o
Hˆ r0 = [(4(0) + 50.00⎝− (−187.8)]kJ / mol ⎠

i products reactants
= 237.8kJ / mol
6) Calculate the overall enthalpy:
ΔH = ξΔHˆ ro + ∑ nout Hˆ out −∑ nin Hˆ in
nC7 H16 ,r 1
ξ& = = =1
vC7 H16 1
ΔH = ξΔHˆ ro + ∑ nout Hˆ out −∑ nin Hˆ in

= (1)(237.8kJ/mol) + [(1)(60.17) + (4)(10.89) - (1)(79.34)]kJ/mol
= 262.2kJ / mol
64
Problems:
1. Use Hess’s law to calculate the standard heat of the water-gas shift reaction;
CO(g) + H2O(v) ÎCO2(g) + H2(g)

from each of the two sets of data given here.
(a) CO(g) + H2O(l) ÎCO2(g) + H2(g): ΔHˆ r0 = +1226 Btu / lbmole
H2O(l) ÎH2O(v): ΔHˆ r0 = +18,935Btu / lbmole
(b) CO(g) + ½ o2(g) ÎCO2(g): ΔHˆ r0 = −121,740 Btu / lbmole

H2(g) + ½O2(g) Î H2O(v): ΔHˆ r0 = −104,040 Btu / lbmole
2. Formaldehyde maybe produced in the reaction between methanol and oxygen:
2CH3OH(l) + O2(g) Î 2HCHO(g) + 2H2O(l):
The standard heat of combustion of hydrogen is
H2(g) + ½ O2(g) Î H2O(l):
Use these heats of reaction and Hess’s Law to determine the standard heat of the
direct decomposition of methanol to form formaldehyde:
CH3OH(l) Î HCHO(g) + H2(g)
3. Trichloroethylene, a widely used degreasing solvent for machine parts, is

produced in a 2 steps reaction sequence. Ethylene is first chlorinated to yield
tetrachloroethane, which is dehydrochlorinated to form trichloroethylene.
C2H4(g) + 2Cl2(g) Î C2H2Cl4(l) + H2(g) ΔHˆ r0 = −385.76kJ / mol
C2H2Cl4(l) Î C2HCl3(l) + HCl(g)
The standard heat of formation of liquid trichloethylene is -276.2 kJ/mol.

(a) Use the given data and tabulated standard heats of formation of ethylene
and hydrogen chloride to calculate the standard heat of formation of
tetrachloroethane and the standard heat of the second reaction.
(b) Use Hess’s Law to calculate the standard heat of the rection
C2H4(g) + 2Cl2(g) Î C2HCl3(l) + H2(g) +HCl(g)
65
4. The standard heat of reaction for the oxidation of ammonia isgiven below:
4NH3(g) + 5O2(g) Î 4NO(g) + 6H2O(v)
500 mol NH3/s and 800 mol O2 at 25 °C are fed into a reactor in which the
ammonia is completely consumed. The product gas emerges at 300 °C. Calculate
the rate at which heat must be transferred to or from the reactor, assuming
operation at approximately 1 atm.
66
CHAPTER 6
BASIC CONCEPT OF HEAT TRANSFER
6.1 Introduction
o What is heat transfer?

Heat transfer is thermal energy in transit due to a temperature difference. In other
words, heat transfer is to predict the energy transfer between material bodies.
[In the simplest of terms, the discipline of heat transfer is concerned with only two
things: temperature, and the flow of heat. Temperature represents the amount of
thermal energy available, whereas heat flow represents the movement of thermal
energy from place to place.]
o How Science of heat transfer explain:

¾ How heat energy may be transfer?
¾ Predict the rate at which the exchange will take place?
¾ Different between thermodynamics and heat transfer ?
o Thermodynamics deal with system in equilibrium, it may used to predict of

energy required to change a system from one equilibrium state to another. But it
does not predict how fast a change will take place.
Example 6.1
Consider the cooling of a hot steel bar which is placed in a pail of water.
Thermodynamics may used to predict the final equilibrium temperature of steel bar-water
combination. However, thermodynamics will not tell:
¾ how long to reach the equilibrium or,

¾ what the temperature of steel bar after certain length of time?
67
Table 6.1: Comparison between heat transfer and thermodynamics
HEAT TRANSFER THERMODYNAMICS
• Heat transfer is a study which

• Thermodynamics deals with
predicts the energy transfer which
systems in equilibrium.
takes place between material
bodies. • Thermodynamics may used to
predict the amount of energy
• It is due to the temperature
required to change a system from
difference
one equilibrium state to another.
• Heat transfer explains how heat
• Thermodynamics may not used to
energy may be transferred.
predict how fast a change will take
place since the system is not in
• It also predicts the rate at which the
equilibrium during the process.
exchange will take place under
certain condition.
Three modes of heat transfer are:

(1) Conduction
(2) Convection
(3) Radiation
6.2 Conduction Heat Transfer
o Energy transfer from high temperature region to low temperature region. We said
that the energy is transferred by conduction. And, the heat transfer rate per unit
area is proportional to the normal temperature gradient.
q/A ~ dT/dx
where,
q/A = heat transfer rate (W/m2)
dT/dx = temperature gradient in the direction of the heat flow
When proportionality constant is inserted,
q = - kA dT/dx [6.1]
o The positive constant k = thermal conductivity of the material. The minus (-) sign
is inserted so that the second principle of thermodynamics will be satisfied i.e.
heat must flow downhill on the temperature scale as indicated in Figure 6.1.
68
Equation (8.0) is called Fourier’s Law of heat conduction. Above is defining
equation for the thermal conductivity and k has the unit of Watts per meter per
Celsius degree (W/m.0C), which the heat flow is expressed in watts.
Figure 6.1: Sketch showing direction of heat flow
Table 6.2: Lists typical values of the thermal conductivities of some metal (thermal
conductivity at 300 K(W/m K)).
Metal k
Copper, pure 396
Aluminium 238
Carbon steel, 1% C 42
Plastics 0.2 - 0.3
Air 0.026
69
Figure 6.2: One –dimensional plane wall
o Consider the one-dimensional plane wall shown in Figure above, if the system in
a steady state, i.e., if the temperature does not change with time, then only
integrate Equation (8.0) and substitute the appropriate values.
x T1
q
A ∫0
dx = ∫ kdT
T2
o Where temperature at the left face (x=0) is uniform at T1 and the temperature at
right face is uniform at T2. If k is independent of T, we obtain after integration :
q T1 − T2
= −k
A L
o Under the steady state conditions, where the distribution is linear, the temperature
gradient may be expressed as:
dT T2 − T1
=
dx L
o And the heat transfer rate:
q T −T
= −k 2 1
A L
or;
q T −T ΔT
=k 1 2 =k
A L L
70
o Since dT/dx = - q/k for the same q, if k is low (i.e: for an insulator), dT/dx will be
large .i.e. there will be a large temperature difference across the wall, and if k is
high (i.e. for a conductor), dT/dx will be small, or there will be a small
temperature difference across the wall.
Example 6.2
One face of a copper plate 3 cm thick is maintained at 400 °C, and the other face
maintained at 100 °C. How much heat is transferred through the plate? Given: k = 370
W/m.K.
From Fourier’s Law,
q dT
= −k
A dx
Integrating gives
q ΔT − ( 370 ) W / m.K (100 − 400) K

= −k =
A Δx 3 × 10−2 m
2
q/A = 3.7 MW/m
6.3 Convection Heat Transfer
o It is well known that hot plate of metal will cool faster when placed in front of a
fan than when exposed to still air. We say that, heat is convected away, and we
call the process, convection heat transfer. The velocity at which the air blows over
the hot plate obviously influence the heat transfer rate.
6.3.1 Mechanism of convection
o Consider the heated plate shown on Figure 8.3. Temperature of plate is Tw and
temperature of the fluid is T∞.The velocity of the flow will appear as shown in
the figure.
71
Figure 6.3: Convection heat transfer from a plate
o The velocity is being reduced to zero at the plate as a result of viscous action.
Since the velocity of fluid layer at the wall will be zero, the heat must be
transferred only by conduction. Thus we might compute the heat transfer, using:
q = -kA (dT/dx)
but with little changes. We use thermal conductivity of fluid and the fluid
temperature gradient at the wall.The temperature gradient is dependent on the rate
at which the fluid carries the heat away. High velocity produces a large
temperature gradient. Thus, the temperature gradient at the wall depends on the
flow field. To express the overall effect of convection, we use Newton’s law of
cooling:
q = hA (Tw-T∞) [6.2]
where,
h = convection heat-transfer coefficient (heat-transfer coefficient = film

conductance)
h = watt per square meter per Celsius degree (W/m² .°C) when the heat flow is in
Watt.
o Heat-transfer rate is related to overall temperature difference between the wall

and fluid and the surface area.
6.3.2 Free and Forced Convection
o If heated plate were exposed to ambient room air without an external source of
motion, a movement of the air would be experienced as a result of the density
gradients near the plate.
o We call this as free convection. When the mass motion of the fluid is caused by
an external device like a pump, compressor, blower or fan, the process is called
forced convection.
Example 6.3
72
Air at 20 °C blows over a hot plate 50 by 75 cm maintained at 250 °C. The convection
heat transfer coefficient is 25 W/m².°C. Calculate the heat transfer?
Solution:
From Newton’s Law of cooling;
q = hA(Tw − T∞ )
2 2
q = (25)W/m .K (0.5x0.75) m (250-20)K
q = 2.156 kW
Example 6.4
An electric current is passed through a wire 1 mm in diameter and 10 cm long. The wire
is submerged in liquid water at atmospheric pressure, and the current is increased until
the water boils. For this situation h = 5000 W/m².°C, and the water temperature will be
100 °C. How much electric power must be supplied to the wire to maintain the wire
surface at 114 °C.
Solution:
Total convection loss is

q = hA(Tw − T∞ )
The surface area of the wire

-3 -2
A = πdL = π (1x10 )(10x10 )
-4 2
= 3.142 x 10 m
Therefore, the heat transfer is

-4
q = (5000w/m².°C)(3.142 x 10 m²)(114-100)
= 21.99 W and this is equal to the electric power which must be applied.
6.4 Radiation Heat Transfer
o Heat may also transfer through regions where perfect vacuum exist. The
mechanism in this case is electromagnetic radiation which is propagated as a
result of a temperature difference. This is called thermal radiation.
o This is confirmed by our experienced everyday experience of being warmed by
the sun, which separated from the earth by approximately 1.5 x 1011 m of empty
space.
o Thermal radiation can of course transport through any ‘transparent’ medium such
as air. Ideal thermal radiator, or black body, will emit energy at a rate
73
q emitted = σAT4 [6.3]
where;
σ is the proportionality constant and is called the Stefan-Boltzmann constant with
-8 2 4
value of 5.669 x 10 W/m .K .
o This equation is called Stefan-Boltzmann law of thermal radiation and it apply to

only blackbodies. This equation is valid only for thermal radiation; other type of
electromagnetic radiation may not be treated so simply.
q emitted = σAT4 govern only radiation emitted by a blackbody.
The net radian exchange between two surface will be proportional to the difference in
4 4
absolute temperature to the fourth power. i.e: q net exchange/ A  (T1 – T2 )
4 4
q net exchange/ A = σ (T – T )
1 2
6.4.1 Blackbody
o It is a perfect emitter of radiation. At a particular temperature the blackbody

would emit the maximum amount of energy possible for that temperature.
o This value is known as the blackbody radiation. It would emit at every
wavelength of light as it must be able to absorb every wavelength to be sure of
absorbing all incoming radiation.
o The maximum wavelength emitted by a blackbody radiator is infinite. It also
emits a definite amount of energy at each wavelength for a particular temperature,
so standard black body radiation curves can be drawn for each temperature,
showing the energy radiated at each wavelength (Fig. 6.4).
74
Figure 6.4: Theoretical black body curve for 5000K
4
o Again, blackbody is a body which radiates energy according to the T law. Other
type of surface, such as glossy painted surface or a polished metal plate, do not
radiate as much energy as the blackbody.
4
o However, the total radiation emitted by this body still generally follows the T1
proportionality. To take account of the “gray” nature of such surfaces we
introduce another factor, called the emissivity,  which relates the radiation of the
“gray” surface to that of an ideal black surface. In addition, we must take into
account the fact that not all the radiation leaving one surface will reach the other
surface since electromagnetic radiation travels in straight line and some will lost
to the surrounding.
o Therefore two new factors in (Eq 8.2)take into account both situation, so that
4 4
q= FF σ A (T - T ) [6.4]
G 1 2
F = emissivity function
F = geometric “view factor” function
G
A simple radiation problem is encounter when we have a heat-transfer surface at

temperature T completely enclosed by much larger surface maintained at T . The net
1 2
radiant exchange in this case can be calculated with,
4 4
q =  Aσ (T - T )
1 1 1 2
where,
 = emissivity of material
1
Example 6.5
0 0
Two infinite black plates at 800 C and 300 C exchange heat by radiation. Calculate the
-8 2 4
heat transfer per unit area. Given: σ = 5.669 x 10 W/m .K
Solution:
4 4
q/A = σ (T1 – T2 )
-8 2 4 4
= 5.669 X 10 )W/m .K (10734 - 5734)K
2
= 69.03 kW/m
75
Problems
1. If 3 kW is conducted through a section of insulating material 0.6 m2 in cross-

section and 2.5 cm thick and the thermal conductivity may be taken as 0.2
W/m0C, compute the temperature difference across the material.
2. A temperature difference of 850C is impressed across a fiberglass layer of 13
cm thickness. The thermal conductivity of the fiberglass is 0.035 W/m.0C.
Compute the heat transferred through the material per hour per unit area.
3. Two perfectly black surfaces are constructed so that all the radiant energy
leaving a surface at 8000C reaches the other surface. The temperature of the
other surface is maintained at 2500C. Calculate the heat transfer between the
surfaces per hour and per unit area of the surface maintained at 8000C.
4. One side of a plane wall is maintained at 1000C, while the other side is
exposed to a convection having T = 100C and h=10 W/m2.0C. The wall has k
= 1.6 W/m.0C and is 40 cm thick. Calculate the heat transfer rate through the
wall.
5. A vertical square plate, 30cm on a side, is maintained at 500C and exposed to
room air at 200C. The surface emissivity is 0.8. Calculate the total heat lost by
both sides of the plate.
CHAPTER 7
HEAT EXCHANGERS
7.1 Introduction
A heat exchanger is a device that facilitates exchange of heat between two fluids that are
at different temperatures and separated by a solid wall. In other words, heat exchangers
are devices built for efficient heat transfer from one fluid to another and are widely used
in engineering processes. The specific applications of heat exchanger are space heating
and air conditioning, power production, waste heat recovery, food and chemical
processing, oil refining, and in vehicles.
76
7.2 Heat exchanger types
Heat exchangers are classified according to flow arrangement and type of construction.
Six (6) types of heat exchanger are:
(a) Concentric Tube Heat Exchangers

(b) Cross-Flow Heat Exchangers
(c) Cross-counter Flow (Coil) Heat Exchangers
(d) Shell-and-Tube Heat Exchangers
(e) Compact Heat Exchangers
(f) Plate-Type Heat Exchangers
(a) Concentric Tube Heat Exchanger

It also called double-pipe heat exchangers or co-axial flow heat exchangers
[Figure 7.1].
One fluid flows inside the tube and the other fluid flows inside the annulus.
In parallel-flow heat exchangers, the two fluids enter the exchanger at the same
end, and travel in parallel to exit at the other side. [Figure7.2a]
In counter-flow heat exchangers the fluids enter the exchanger from opposite ends.
Counter current design is more efficient, in that it can transfer more heat [Figure
7.2b].
Figure 7.1: Concentric tube heat exchangers
(a) Parallel flow (b) Counterflow
Figure 7.2: Concentric tube heat exchangers
77
(b) Cross Flow Heat Exchanger
In a cross-flow heat exchanger, the fluids travel roughly perpendicular to one
another through the exchanger.
In finned tubular heat exchangers, the fin-side fluid is unmixed because the fins
confine the flow to one direction. Example: automobile radiator.
In unfinned tubular heat exchangers, the fin-side fluid is mixed because the flow
in transverse direction is possible. [Figure 7.3]
The use of fins to improve the convection coefficient of fin-side fluid by
increasing the outside surface area.
(a) Finned with both fluids (b) Unfinned with one fluid mixed
unmixed and the other unmixed
Figure 7.3: Cross-flow heat exchangers
(c) Cross-counter Flow (Coil) Heat Exchangers
78
(a) Two-pass (b) Four-pass
Figure 7.4: Cross-counter Flow (Coil) Heat Exchangers
(d) Shell-and-Tube Heat Exchangers
Tubular heat exchangers consist of a tube bank enclosed by a shell. One fluid
flows inside the tubes and the other flows inside the shell.
Figure 8.5, show the simplest form of shell and tube heat exchanger which
involves single tube and shell passes.

Figure 7.5: Shell and tube heat exchanger with one-shell pass one-tube pass (1-1)
(cross-counterflow mode of operation)
Baffles are usually installed to increase the convection coefficient of the shell-side
fluid by inducing turbulence and a cross-flow velocity component.
79
Baffled heat exchanger with one shell pass and two tube passes and with two shell
passes and four tube passes are shown in Figures 7.6a and 7.6b, respectively.
Figure 7.6: Shell-and-tube heat exchanger. (a) One shell pass and two tube passes. (b)
Two shell passes and four tube passes
Shell-and-tube heat exchangers are classified according to the number of shell and
tube passes involved.
Heat exchangers in which all the tubes make one U-turn in the shell, for example,
are called one-shell-pass and two-tube-passes heat exchangers.
Heat exchanger that involves two passes in the shell and four passes in the tubes
is called a two-shell-passes and four-tube-passes heat exchanger [Figure 7.7]
80
(a) One-shell pass and two-tube passes
(b) Two-shell passes and four-tube passes
Figure 7.7: Multipass flow arrangement in shell-and-tube heat exchangers
(e) Compact Heat Exchanger
Used to achieve very large heat transfer area per unit volume.
Have dense arrays of finned tubes or plates.
Typically used when at least one of the fluids is a gas, characterized by small
convection coefficient.
The tubes may be flat or circular, and the fins may be plate or circular (Figure
7.8a, 7.8b and 7.8c).
Parallel plates may be finned or corrugated and may be used in single-pass or
multi-pass mode (Figure 7.8d and 7.8e).
81
Figure 7.8: Compact heat exchanger cores. (a) Fin-tube (flat tubes, plate fins; (b) Fin-tube
(circular tubes, plate fins); (c) Fin-tube (circular tubes, circular fins); (d) Plate-fin (single
pass); (e) Plate-fin (multi-pass)
(f) Plate-Type Heat Exchangers
• Gasketed plate exchanger, alternative to shell-and-tube exchangers for

applications at moderate temperature and pressure.
• Consists of many corrugated stainless steel sheets separated by polymer gaskets
and clamped in a steel frame.
• Inlet portals in the gaskets direct the hot and cold fluid to alternate spaces between
plates.
• Adjoining plates have different pattern or angle of corrugation. Corrugations
induce turbulence for improved heat transfer.
• Widely used in dairy and food processing
82
(a) General layout (b) Detail of plate design
Figure 7.9: Plate-Type Heat Exchangers
The Overall Heat Transfer Coefficient
• The heat transfer between the two fluids across the solid wall involves convection
of fluid films adjacent to the wall and conduction across the wall.
• The rate of heat transfer can be expressed by a single equation like Newton’s law
of cooling, with the overall heat transfer coefficient U incorporating convection
and conduction terms:
q = UAΔTm [7.1]
where,
ΔTm = mean temperature difference between
the two fluids along the exchanger length
• For the unfinned tubular heat exchanger, U can be calculated as follows:

1 1 1
= = [7.2]
UA U i Ai U o Ao
1 1 R ln( Do / Di ) R fo 1
= + fi + + + [7.3]
UA hi Ai Ai 2π k L Ao ho Ao
83
where Ai , Ao = inside and outside heat transfer areas
U i , U o = overall heat transfer coefficients based on
inside and outside surface areas
hi , ho = inside (tube-side) convection coefficient and
outside (shell-side) convection coefficient
R fi , R fo = fouling resistances at inside and outside surfaces
Di , Do = inside and outside diameters of the tube
k = thermal conductivity of the tube wall
L = tube length of heat exchanger
7.3 Heat Exchanger Analysis: Use of the Log Mean Temperature Difference
(LMTD)
To design or to predict the performance of a heat exchanger, it is essential to relate the

total heat transfer rate to quantities such as the inlet and outlet fluid temperatures, the
overall heat transfer coefficient, and the total surface area for heat transfer. A few steps
to design or predict the performance of a heat exchanger:
• Step 1
- Write down the overall energy balances between heat gain of cold fluid, heat
loss of hot fluid, and heat transfer across the wall separating the two fluids.
) )
q = m& c ΔH c = m& h ΔH h = UAΔTlm [7.4]
- If the fluids are not undergoing a phase change and constant specific heats are
assumed, the equation becomes:
q = m& hC p ,h (Th ,i − Th ,o ) = m& c C p ,c (Tc ,o − Tc ,i ) = UAΔTlm [7.5]
- Determine any unknown that can be directly calculated from the above
relations.
84
q = rate of heat transfer
m& c , m& h = mass flow rate of cold fluid and hot fluid
Tci , Tco , Tc = inlet, outlet, and mean temperature of cold fluid
Thi , Tho , Th = inlet, outlet, and mean temperature of hot fluid
) )
H c , H h = specific enthalpies of cold fluid and hot fluid
U = overall heat ransfer coefficient
A = heat transfer surface area
ΔTlm = log mean temperature difference (LMTD)
ΔT1 − ΔT2
= , where ΔT1 and ΔT2 are
ln ( ΔT1 ΔT2 )
temperature differences at the two ends of HE.
ΔT1
ΔT2
) )
q = m& c ΔH c = m& h ΔH h = UAΔTlm
Enthalpy change of cold Enthalpy change of hot

fluid fluid
Figure 7.10: Overall Energy Balances of the Hot and Cold Fluids of a Two-fluid Heat
Exchanger
• Step 2
- Enthalpies and LMTD depends on fluid temperature behavior:
(a) Parallel flow
For parallel flow with no phase change,
ΔHh = Cph(Thi-Tho), ΔHc = Cpc(Tco-Tci)
85
(Thi − Tci ) − (Tho − Tco )
ΔTlm =
ln[(Thi − Tci ) (Tho − Tco )]
Figure 8.11: Temperature scheme for parallel flow
(b) Counterflow
For counterflow with no phase change,
ΔHh = Cph(Thi-Tho), ΔHc = Cpc(Tco-Tci)
(Tho − Tci ) − (Thi − Tco )

ΔTlm =
ln[(Tho − Tci ) (Thi − Tco )]
Figure 8.12: Temperature scheme for counterflow
(c) Condensers
ΔHh = λh, ΔHc = Cpc(Tco-Tci)
Tco − Tci
ΔTlm =
ln[(Th − Tci ) (Th − Tco )]
86
Figure 8.13: Temperature scheme for condenser
(d) Evaporators
ΔHc = λc, ΔHh = Cph(Thi-Tho)
Thi − Tho
ΔTlm =
ln[(Thi − Tc ) (Tho − Tc )]
Figure 8.14: Temperature scheme for evaporator
where λ= latent heat of vaporization
- Miscellaneous
• Overall heat transfer coefficient (U), if not known, can be determined
Equation (10.3).
• Heat transfer area (A) is related to tube length as A = 2πrL.
• Step 3
- Determine any other unknowns from the overall energy balances (Eq. 8.4)
- Determine the HE effectiveness and number of transfer units.
87
• If the heat exchanger other than the double pipe is used, the heat transfer is
calculated by using a correction factor applied to the LMTD for a counter flow
double-pipe arrangement with the same hot and cold fluids temperatures.
• The heat-transfer equation becomes :
q = UAF ΔTm [7.6]
• Values of the correction factor F are plotted in Figures 7.14 to 7.17 for several
different types of heat exchangers.
• When phase changed is involved, as in condensation or boiling (evaporation), the
fluid normally remains at essentially constant temperature.
• For this condition, P and R becomes zero and F = 1.0 (for boiling and
condensation)
Figure 7.15: Correction factor plot for exchanger with one shell pass and two, four, or
any multiple of tube passes.
88
Figure 7.16: Correction factor plot for exchanger with two shell pass and four, eight, or
any multiple of tube passes.
89
Figure 7.17: Correction-factor plot for single-pass cross-flow exchanger, both fluids
unmixed.
90
Figure 7.18: Correction-factor plot for single-pass cross-flow exchanger, one fluid mixed,
the other unmixed.
91
Example 7.1
Water at the rate 68 kg/min is heated from 35 to 750C by an oil having specific heat of
1.9 kJ/kg.0C. The fluids are used in a counterflow double pipe heat exchanger and the oil
enter the exchanger at 1100C and leaves at 750C. The overall heat-transfer coefficient is
320 W/m2.0C. Calculate the heat exchanger area.
Solution:
The total heat transfer is determined from the energy adsorbed by the water:
.
q = m w Cw ΔTw
= (68)(4180)(75 − 35)
= 11.37 MJ / min
= 189.5kW
Since all the fluid temperature are known, the LMTD can be calculated by using the
temperature scheme in Figure 10.12,
(110 − 75) − (75 − 35)
ΔTm = = 37.440 C
ln(110 − 75) /(75 − 35)
Then, since q = UAΔTm
1.895 x105
A= = 15.82m 2
(320)(37.44)
Example 7.2
Instead of the double-pipe heat exchanger of Example 1 above, it is desired to use a shell-
and-tube exchanger with the making one shell pass and the oil making two tube passes.
Calculate the area required for this exchanger, assuming that the overall heat-transfer
coefficient remains at 320 W/m2.0C.
Solution:
To solve this problem we determine a correction factor from Figure x to be used with the
LMTF calculated on the basis of a counterflow exchanger. The parameters according to
the nomenclature of Figure 8.15 are:
T1=35 0C T2=75 0C t1=110 0C t2=75 0C
92
t2 − t1 75 − 110
P= = = 0.467
T1 − t1 35 − 110 t2 − t1 75 − 110
P= = = 0.467
T1 − T2 35 − 75 T1 − t1 35 − 110
R= = = 1.143
t2 − t1 75 − 110 T1 − T2 35 − 75
R= = = 1.143
1.895 x10 5 t2 − t1 75 − 110
A= = 19.53m 2
(320)(0.81)(37.44)
So the correction factor is
F = 0.81
And the heat transfer is
q = UAFΔTm
so that
1.895 x105
A= = 19.53m 2
(320)(0.81)(37.44)
Example 7.3
A cross flow heat exchanger is used to heat an oil in the tubes, c = 1.9 kJ/kg.0C from
150C to 850C. Blowing across the outside of the tubes is steam which enters at 1300C and
leaves at 1100C with mass flow of 5.2kg/sec. The overall heat transfer coefficient is 275
W/m2.0C and c for steam is 1.86 kJ/kg.0C. Calculate the surface area of the heat
exchanger.
Solution:
The total heat transfer may be obtained from an energy balance on the steam:
.
q = m s cs ΔTs = (5.2)(1.86)(130 − 110) = 193kW
We can solve for the area from Eq.(8.6). The value of ΔTm calculated is as if the
exchanger were counterflow double pipe.
Thus,
(130 − 85) − (110 − 15)
ΔTm = = 66.90 C
⎛ 130 − 85 ⎞
ln ⎜ ⎟
⎝ 110 − 15 ⎠
From Figure 8.18, t1 and t2 will represent the unmixed fluid (oil) and T1 and T2 will
represent the mixed fluid (the steam) so that
T1 = 1300C T2 = 1100C t1 = 15 0C t2 = 850C
93
and we calculate
130 − 110
R= = 0.286
85 − 15
85 − 15
P= = 0.609
130 − 15
so the area is calculated from

q
A=
UF ΔTm
7.4 Heat Exchanger Analysis: Use of the Effectiveness-NTU method
7.4.1 Introduction
The driving temperature across the heat transfer surface varies with position, but an
appropriate mean temperature can be defined. In most simple systems this is the log mean
temperature difference (LMTD). Sometimes direct knowledge of the LMTD is not
available and the NTU method is used.
• The LMTD approach to heat-exchanger analysis is useful when inlet and outlet
temperature are known or are easily determined.
• However, when the inlet or exit temperatures are to be evaluated for a given heat
exchanger, the analysis involves an interactive procedure because of the
logarithmic function in the LMTD.
• In these cases the analysis is performed more easily by utilizing a method based
on the effectiveness of the heat exchanger in transferring a given amount of heat.
• The heat exchanger effectiveness can be define as:
Actual rate of heat transfer

Effectiveness, ε =
Maximum possible rate of heat transfer
• The actual rate heat transfer, q may be computed by calculating either the energy
lost by the hot fluid or the energy gained by the cold fluid.
• For the parallel-flow exchanger:

. .
q = m h ch (Th1 − Th 2 ) = m c cc (Tc 2 − Tc1 )
• For the counter-flow exchanger :
94
. .
q = m h ch (Th1 − Th 2 ) = m c cc (Tc1 − Tc 2 )
• The maximum possible heat transfer, qmax is the rate of heat transfer that a heat
exchanger of infinite area would transfer with given inlet temperatures, flow rates,
and specific heat.
• qmax occurs when the fluid with minimum product of flow rate and specific heat
changes temperature to the entering temperature of the other fluid.
• Maximum possible heat transfer is expressed as

.
qmax = (m c) min (Thi − Tci )
• The minimun fluid may be either hot or cold fluid depending on mass-flowrates,
m and specific heats, c.
• For the parallel exchanger :

.
. .
m h ch (Th1 − Th 2 )
T −T
ε= . = h1 h 2
m h ch (Th1 − Tc1 ) Th1 − Tc1
.
. .
mc cc (Tc 2 − Tc1 )
T −T
ε= . = c 2 c1
mc cc (Th1 − Tc1 ) Th1 − Tc1
Figure 7.19 : Temperature profile for parallel exchanger
95
• For the counterflow exchanger :
.
. .
m h ch (Th1 − Th 2 )
T −T
ε= . = h1 h 2
m h ch (Th1 − Tc 2 ) Th1 − Tc 2
.
. .
mc cc (Tc1 − Tc 2 )
T −T
ε= . = c1 c 2
mc cc (Th1 − Tc 2 ) Th1 − Tc 2
Figure 7.20: Temperature profile for parallel exchanger
• In the general way the effectiveness is expressed as :
ΔT ( minimum fluid)
ε=
Maximum temperature difference in heat exchanger
• The effectiveness is usually written for parallel flow double pipe HE :

. . .
1 − exp[(−UA / mc cc )(1 + mc cc / m h ch )
ε= . .
1 + mc cc / m h ch
1 − exp[(−UA / Cmin )(1 + Cmin / Cmax )

ε=
1 + Cmin / Cmax
.
where C = m c = capacity rate
• The number of transfer units (NTU) is indicative of the size of the heat exchanger.
NTU = UA / Cmin
• Figure 7.21 to 7-26 presented effectiveness ratios for various heat exchanger
arrangements.
• Table 7.1 and Table 7.2 summarizes the effectiveness relations
96
Figure 7.21: Effectiveness for parallelflow Figure 7.22: Effectiveness for counterflow
exchanger performance exchanger performance
97
Figure 7.23 : Effectiveness for crossflow Figure 7.24: Effectiveness for crossflow
exchanger with one fluid exchanger with one fluid
mixed unmixed
Figure 7.25 : Effectiveness for 1-2 Figure 7.26 : Effectiveness for 2-4
parallel counterflow multipass counterflow
exchanger performance exchanger performance
98
Table 7.1: Heat exchanger effectiveness relations
99
Table 7.2: NTU relations for heat exchangers
Problems
1. Water flowing at a rate of 0.723 kg/s enters the inside of a countercurrent, double-
pipe heat exchanger at 300 K and is heated by an oil stream that enters at 385 K at
a rate of 3.2 kg/s. The heat capacity of the oil is 1.89 kJ/kg K, and the average
heat capacity of water over the temperature range of interest is 4.192 kJ/kg K. The
overall heat-transfer coefficient of the exchanger is 300 W/m2.K, and the area for
heat transfer is 15.4 m2. What is the total amount of heat transferred?
2. A shell and tube exchanger operates with two passes and four tube passes. The
shell fluid is ethylene glycol (c= 2.742 kJ/kg.0C), which enters at 1400C and
leaves at 80 0C with flowrate of 45000 kg/h. Water flows in tubes (4.174
kJ/kg.0C), entering at 35 0C and leaving at 85 0C. The overall heat transfer
coefficient for this arrangement is 850 W/m2.0C. Calculate the flow rate of water
required and the area of the heat exchanger.
3. A double pipe heat exchanger is used to heat an oil with c = 2.2 kJ/kg.0C from 50
0
C to 1000C. The other fluid having c = 4.2 kJ/kg.0C enters the exchanger at 160
0
C and leaves at 900C. The overall heat transfer coefficient is 300 W/m2.0C.
Calculate the area and effectiveness of the heat exchanger for a total heat transfer
rate of 600 kW.
100
REFERENCES
1. Felder, R.M., Rousseau, R.W., Elementary Principles of Chemical Processes,

Third Edition, John Wiley & Sons, 2000.
2. Frank P. Incropera and David P. De Witt, Introduction to Heat Transfer,
Fourth Edition, John Wiley & Sons, 2002.
3. Frank P. Incropera, David P. De Witt, Theodore, L. Bergman, Adrienne S.
Lavine, Fundamentals of Heat and Mass Transfer, Fourth Edition, John Wiley
& Sons, 2007.
4. J.P.Holman, Heat Transfer, Ninth Edition, Mc Graw Hill,2002.
5. Yunus A. Cengel, Heat Transfer: A Practical Approach, Second Edition, Mc

Graw Hill, 2003.
101
View publication stats

Overall PDF

Uploaded by

Copyright:

Available Formats

Overall PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Overall PDF

Uploaded by

Copyright:

Available Formats

See discussions, stats, and author profiles for this publication at: https://www.researchgate.

TEACHING MODULE: MATERIAL BALANCES, ENERGY BALANCES AND HEAT

Book · September 2010

Norazian Mohamed Noor Masitah Hasan

SEE PROFILE SEE PROFILE

Variation of particulate matter (PM10) in Malaysia View project

The user has requested enhancement of the downloaded file.

Biotechnology Lime and cement

Figure 1.1: Chemical process

1.2 Units and Dimensions

“What are units and dimensions and how do they differ?”

¾ Values could be added if UNITS are the same.

* Same dimension: length, different units -- can add

Table 1.2: American Engineering System Units

Physical Quantity Name of Unit Symbol

Table 1.3: Factors for unit conversions

Convert an acceleration of 1 cm/s2 to its equivalent in km/yr2.

1cm 1m 1km 3600 2 s 2 24 2 hr 2 365 2 day 2

Convert 400 in3/day to cm3/min.

1.4 Processes and Process Variables

A process is any operation or series of operations by which a particular objective is

1.3.1 Instructional Objectives

The objectives in studying this section are to be able to:

4. Convert a pressure expressed as a head of a fluid to the equivalent pressure

5. Convert a manometer reading into a pressure difference for an open end

6. Convert among temperatures expressed in K, °C, °F and °R.

A process is any operation or series of operations that cause a physical or chemical

Figure 1.2: Process streams with various process variables.

Table 1.4: Process variables

35 cm3 1.595 g = 55.825 g

and the volume of 9.3 lbm of CCl4 is

9.3 lbm 454 g cm3

ρH2O (4 °C) = 1 g/cm3

Volume = 215 kg lbm ft3 = 0.56 ft3

A liquid has a specific gravity of 0.5.

1.3.4: Flow rate

Figure 1.3: Flow rate

The volumetric flow rate of CCl4 in a 1.0-cm-diameter pipe is 100 cm3/min.

i. What is the mass flow rate?

ii. What is the molar flow rate?

= 159.5 g-CCl4/min x g mol-CCl4/153.838 g-CCl4

iii. What is the linear velocity of CCl4?

(a) What is the volumetric flow rate of the hexane?

(Answer: (a) 10 cm3/s; (b) 0.509 cm/s)

1.3.5: Chemical Composition

1.3.5(a): Moles and Molecular Weight

100 g CO2 1 mol CO2

= 2.273 mol CO2

ii. lb-moles CO2

2.273 mol CO2 1 lb-mol = 5.011 x 10-3 lb-mole CO2

2.273 mol CO2 1 mol C = 2.273 mol C

2.273 mol CO2 2 mol O = 4.546 mol O

2.273 mol CO2 1 mol O2 = 2.273 mol O2

4.546 mol O 16 g O = 72.7 g O

2.273 mol O2 32 g O2 = 72.7 g O2

viii. molecules of CO2

2.273 mol CO2 6.02 x 1023 molecules = 1.37 x 1024 molecules

(a) 50 g mol of CaCO3.

(Answer: (a) 0.11 lb mol; (b) 3.3 lb mol; (c) 1 lb mol)