Energy balance (Notes - Prof Harding)
Energy balance (Notes - Prof Harding)
0)
1 Energy balances
1.1 Introduction Felder (7.0)
1.1.1 Background
As for mass balances, the conservation laws extend in an analogous way for energy. This
chapter will cover basic introductory concepts and then energy balances over non-reacting
systems as well as over reacting systems.
Before being able to calculate an energy balance it is important to understand the terms needed
for such [Move to chapter 1?]. These include terms which may already be familiar to you. For
a full list of terms as well as shorter definitions, see the Glossary.
https://youtu.be/oclirLhr4k0?list=PLDb9sXA9R_c7LeqVujyRRfR80HvWp9rPv
Video 1.2: Steady state (YouTube)
https://youtu.be/HRLCOW_wWWM?list=PLDb9sXA9R_c7LeqVujyRRfR80HvWp9rPv
Video 1.3: Systems: Open, closed and isolated systems (YouTube)
https://youtu.be/oAdSZ0MhLbY?list=PLDb9sXA9R_c7LeqVujyRRfR80HvWp9rPv
Video 1.4: State vs path functions (YouTube)
Iso-
https://youtu.be/vVs5H6cvT1c?list=PLDb9sXA9R_c7LeqVujyRRfR80HvWp9rPv
Video 1.5: Iso- (isothermal, isochoric, isobaric, …) (YouTube)
1.1.2 Heat
Heat (Q, in units of kJ, kJ/s, kW or kJ/kg) is the disorganized transfer of energy that results due
to a temperature gradient. Heat is not included in the energy of a system but rather as the
transfer of energy across the system boundary. Heat (transfer) can be through conduction (no
movement of fluids), convection (movement of fluids) or radiation (no temperature gradient
needed). More information on this is given in a separate chapter.
Relating heat to energy balances and thermodynamics, for a closedm system heat, energy and
work are related by:
DE = Q + W 1.1
Where:
E = Energy
Q = Heat, and
W = Work
And DE = Ef – E0, is the change in energy from the initial to final states, inside the system.
1.1.3 Work
Work is typically the transfer of energy through movement due to a force. It is typically broken
down into the following subcategories:
- Expansion work (W∆" ) [mechanical work], also called PV Work (WPV), is the work needed to
change a systems volume.
- Flow work (WFlow) is the work performed when a stream exits or enters a system.
- Shaft work (WS) is the work due to mechanical work changing the pressure, e.g. from pumps
compressors, turbines etc.
- Electrochemical work (Wel) is the work performed due to such things as batteries or fuel cells.
- Other forms of work (Wother) might include electromagnetic work, sound work, etc.
Where:
Wn is all the non-flow work
Solution:
1.1.4 Energy
Energy is often referred to as the “ability to perform work”. However, the definition needs to
further extend to include the concept of being able to release heat. Therefore, the definition of
energy can be said to be either the ability to perform work or to release heat.
- Kinetic energy (EK) is the energy due to movement of the centre of mass of an object.
- Rotational energy (ER) is the energy when the centre of mass is at rest but there is rotation.
- Potential energy (EP) is the energy due to an elevated height relative to a reference point.
- Electric energy (Eel) is energy stored in such things as a capacitor.
- Surface energy (ES) is the intermolecular energy creating a surface.
- Internal energy (U) is the energy of molecules. It can be seen as the sum of thermal energy,
chemical bonding energy and latent energy on a molecular level.
For the sake of this course, we will only be looking at kinetic energy (EK), potential energy
(EP) and internal energy (U).
For a system with no potential, kinetic, rotational etc. energy, i.e. no movement or energy due
to height or other external factors, this can be simplified to:
E=U 1.6
Therefore:
DU = Q + W 1.7
where the internal energy term is an ordinary differential, while Q and W are partial
differentials.
For a constant volume (closed) system, where all W terms are also zero, we can simplify this
as:
dU = δQ 1.9
Or
DU = Q 1.10
Example 1.2:
Solution:
1.1.6 Enthalpy
NEED TO INCLUDE CP DEFINITION IN HIMMELBLAU (5.xxx)
Enthalpy is defined as:
H = U + PV 1.11
where:
H: Enthalpy (J)
U: Internal Energy (J)
P: System pressure (N/m2)
V: System volume (m3)
Enthalpy is measured according to some known value. This is known as the basis. Typically,
the basis is chosen as elements in their standard state (at 1 atm and 298.15K) à H = 0.
Given a stream of CO(g), flowing at 100 mol/min at a temperature of 298.15 K and pressure of 1 bar, calculate
the enthalpy.
Solution:
Basis: Elements in standard state at 298.15K and 1 bar
Enthalpy flow:
𝑚𝑜𝑙 𝑘𝐽
𝐻̇ = 100 × (−110.53)
𝑚𝑖𝑛 𝑚𝑜𝑙
𝑘𝐽
𝐻̇ = −11 053
𝑚𝑖𝑛
Converting:
𝑘𝐽 1000 𝐽 1𝑚𝑖𝑛
𝐻̇ = −11. 053 × 10# × × = −1.84 × 10$ 𝐽/𝑠
𝑚𝑖𝑛 𝑘𝐽 60𝑠
𝐻̇ = −1.84 𝑥 10$ 𝑊 = −184 𝑘𝑊
The negative value indicates that energy is being released, i.e. exothermic
In addition to 25°Cand 1atm basis entails elements in standard state having H0f = 0 kJ/mol (this
accounts for formation and breaking of bonds in a reactor).
It is possible to measure these properties under other conditions for use as a basis. However,
25°C and 1 atm is the most convenient for measurement.
Definition:
"#
C.! = / 0 "$ %&,(
1.12
Where:
𝑥2 = composition
P = Pressure
i.e., constant composition and pressure
Where:
T1 = reference temperature (K)
T2 = new temperature (K)
$% $% "#
∫$& C.! dt = ∫$& / "$ 0 dt = H|$% − H|$& 1.14
%&,(
i.e., if we don’t have Cp data, we can use the enthalpy data (if available)
NOTE:
If Cp is constant, i.e., NOT dependent on T, then:
$%
$%
5 C.! dt = C.! T;$ = C.! (T+ − T, ) = C.! ∆T
$& &
If Cp is a function of temperature the integral becomes more complex and integral tables may
be needed.
Kevin Harding ([email protected]) 6
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – Introduction Felder (7.0)
Given a stream of CO(g), flowing at 100 mol/min at a temperature of 350 K and pressure of 1 bar, calculate
the enthalpy.
Solution:
∴ 𝐻̇ = −11.053 × 10! 𝑘𝐽. 𝑚𝑖𝑛"# + 𝑛̇ 𝐶4$ ∆𝑇
𝐻̇ = −11.053 × 10! 𝑘𝐽. 𝑚𝑖𝑛"# + [100𝑚𝑜𝑙. 𝑚𝑖𝑛"# × 29.14𝐽. 𝐾 "# 𝑚𝑜𝑙"# × (350 − 298.15)𝐾]
𝐻̇ = −11.053 × 10! 𝑘𝐽. 𝑚𝑖𝑛"# + 151 090.9𝐽. 𝑚𝑖𝑛"#
𝐻̇ = −11.053 × 10! 𝑘𝐽. 𝑚𝑖𝑛"# + 151.0909𝑘𝐽. 𝑚𝑖𝑛"#
𝐻̇ = −10.902 × 10! 𝑘𝐽. 𝑚𝑖𝑛"#
NOTE:
From previous examples, we see that H = -11.053 x 103 kJ/min when T = 298.15K, and when
T = 350 K, H = -10.902 x 103 kJ/min. Therefore, energy is taken in when the reaction proceeds
at a higher temperature of 350K
(More negative to less negative à energy taken in)
Logically, to move from -10°C to -5°C, energy is taken in as we are heating the substances.
Therefore, a change from a negative larger magnitude value to a negative smaller magnitude
value, energy is absorbed or taken in.
Note:
- We cannot measure specific heat directly; we measure enthalpy relative to some basis.
- In practice, enthalpies are measured, and a polynomial model is fitted to the data in order to
obtain specific heat. The model has different constants depending on which units of
temperature, say, were used to make the enthalpy measurements
Cp vs Cv
1.1.11 Thermodynamics
- Zeroth Law
- First Law of Thermodynamics
- Second Law of Thermodynamics
Kevin Harding ([email protected]) 7
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – Steam Tables
https://youtu.be/oZBKAF7_Z7E?list=PLDb9sXA9R_c7LeqVujyRRfR80HvWp9rPv
Video 1.7: Thermodynamic laws (YouTube)
The data is often compiled into three main parts for water (steam):
- Saturated Water and Steam;
- Superheated and Supercritical Steam; and
- Further Properties of Water and Steam
While colloquially called steam tables, most compilations also include a section on the
properties of other materials such as Mercury, Ammonia and CFCs.
Felder 7.5
Solution:
(a) 909 kJ/kg
(b) 2743 kJ/kg
(c) 0.06665 m3/kg
(d) 0.0010030 m3/kg
(e) 3052 kJ/kg
(f) 2563 kJ/kg
(g) 572 kJ/kg
We have saturated steam at 10 bar with 𝑉= = 0.02m3/kg. How much liquid and vapor do we have?
Solution:
For saturated steam:
P = 10 bar; Ts = 179.9°C; 𝑉=' = 0.1944 m3/kg
But the specific volume of the system is only 0.02 m3/kg. Therefore, the system is not only steam; there must
be some liquid present in equilibrium with the steam.
By interpolation:
(0.1128 × 10)* − 0.1114 × 10)* )𝑚# E𝑉=( − 0.01114 × 10)* F𝑚#
𝑘𝑔 𝑘𝑔
=
(10.03 − 7.920)𝑏𝑎𝑟 (10 − 7.920)𝑏𝑎𝑟
(0.1128 × 10)* − 0.1114 × 10)* )𝑚# /𝑘𝑔 × (10 − 7.920)
𝑉=( = + 0.01114 × 10)* 𝑚# /𝑘𝑔
(10.03 − 7.920)
𝑉=( = 1.125 × 10)* 𝑚# /𝑘𝑔
Assuming 1 kg of steam:
𝑀+ 𝑉I+ = 𝑀( 𝑉I( + (𝑀+ − 𝑀( ) 𝑉I'
𝑀+ 𝑉I+ = 𝑀( 𝑉I( + 𝑀+ 𝑉I' − 𝑀( 𝑉I'
𝑀( E𝑉I( − 𝑉I' F = 𝑀+ E𝑉I+ − 𝑉I' F
Steam quality derivation
(0.02 − 0.1944)𝑚# /𝑘𝑔
𝑀( = 1𝑘𝑔
(1.125 × 10)# − 0.1944)𝑚# /𝑘𝑔
𝑀( = 0.9024𝑘𝑔
Therefore:
Mg = (1 – 0.9024) kg AND Mg = 0.0976 kg
Find the enthalpy of steam at T = 175°C and P = 5 bar. Interpolate between saturation point and 200°C.
Solution:
At saturation: T = 151.8°C; P = 5 bar; hg = 2749 kJ/kg
And at T = 200°C and P = 5 bar: hg = 2857 kJ/kg
Therefore:
2857 − 2749
ℎ' = K L . (175 − 151.8) + 2749
200 − 151.8
hg = 2801 kJ/kg
Steam quality
The quality of steam is defined as the mass fraction (or percentage) of vapour in the mixture.
-
x = -, 1.15
-
And, by definition:
V = mVA 1.17
Divide by mT:
1 --
0 1 . -.
0 1 , -,
0
--
= --
+ --
1.19
1 1
A = 0. - . + 0, - ,
V 1.20
- - - -
We obtain:
-- - -
--
− x = -- − -, 1.21
- -
-- 3 -,
1−x= --
1.22
-
1 − x = -. 1.23
-
And re-arranging:
A=V
V A. − x. V
A. + x. V
A/ 11
A−V
V A. = x. DVA/ − VA. E 12
01 301.
x= 1 , 30
0 1.
13
[4 x diagrams here]
Example 1.8:
Solution:
In Out
When considering energy this can be in the form of energy due to material flowing through the
system, work or added heat.
Min Mou
t
Q
Kevin Harding ([email protected]) 12
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – General Energy Balance Equation
where:
M = Mass (kg or kg/s)
W = Work (J or J/s)
Q = Heat (J or J/s)
Rearranging:
0 = Energy from flow45 − Energy from flow678 + Q̇ + Ẇ 1.27
But this is only true for steady state – if not at steady state:
(Change in Energy) = (Energy from flow45 − Energy from flow678 ) + Q̇ + Ẇ 1.28
Energy is going to be defined as the internal energy (U) of the system/flow terms:
9:
98
= (U45 − U678 ) + Q̇ + Ẇ 1.29
9:
98
= ((H − PV)45 − (H − PV)678 ) + Q̇ + Ẇ 1.31
Rearranging:
9:
98
= (H45 − H678 ) + Q̇ + Ẇ + Ẇ(0 1.32
9:
98
= ∑54;, H4 + Q̇ + Ẇ + Ẇ(0 1.33
What about kinetic and potential energy? These are forms of energy due to movement (flow
terms) or change with time (‘/dt’ terms).
Or:
1%
9A:= BC = DEF
% H%
98
= ∑54;, Ṁ G /H
A+
+
+ φ0 + Q̇ + Ẇ + Ẇ(0 1.35
4
Typically, the kinetic and potential energy terms are ignored as they do not influence
temperature and work terms are written as a single term, so:
9:
= ∑54;,DMH A E + Q̇ + Ẇ + Ẇ(0 1.36
98 4
Or as a molar basis:
9:
= ∑54;, /NH0 + Q̇ + Ẇ + Ẇ(0 1.37
98 − 4
Adiabatic process
9:
= ∑54;,DMHA E + Ẇ + Ẇ(0 1.41
98 4
Sign Conventions
[Q]: What does positive or negative Q represent?
Heat Transfer
[A]: Positive Q is heat transferred from hot to cold, i.e., out of the system.
Therefore, definition of heat transfer is opposite in Energy Balance calculations to Heat Transfer
calculations.
BEWARE: When looking at both types of approaches in a single problem ensure the correct sign for Q.
Example 1.9:
Using Simpson’s Rule, determine the enthalpy from 20°C to 200°C given the data below.
Solution:
Note:
h = constant value of 30
7 points à 6 intervals (even number)
Therefore, we can use Simpson’s
Now:
*""
𝐻= $ 𝐶𝑝(𝑥)
*"
5 2
Simpson’s Rule: ∫6 𝑓(𝑥) ≈ [𝑓(𝑥3 ) + 4𝑓(𝑥* ) + 2𝑓(𝑥# ) + ⋯ + 2𝑓(𝑥4)* ) + 4. 𝑓(𝑥4)3 ) + 𝑓(𝑥4 )]
#
Therefore:
*""
ℎ
𝐻= $ 𝐶𝑝(𝑥) ≈ [𝐶𝑝3 + 4𝐶𝑝* + 2𝐶𝑝# + 4𝐶𝑝7 + 2𝐶𝑝$ + 4𝐶𝑝8 + 𝐶𝑝9 ]
*" 3
𝐻 ≈ 5305.6𝐽/𝑚𝑜𝑙
NOTE:
If we were asked to find the integral up to 230°C, this would have been an odd number of intervals and
Simpson’s would not work. Therefore, would need to use Simpson’s as above AND then add the last portion
by another technique, e.g., Trapezoidal Rule.
Example 1.10:
Solution:
Example 1.11:
Solution:
Example 1.12:
Solution:
Add FH approximation
Example 1.13:
Solution:
Example 1.14:
Solution:
The Joback method (often named Joback/Reid method) predicts eleven important and
commonly used pure component thermodynamic properties from molecular structure only.
[expand on this]
Example 1.15:
Solution:
Example 1.16:
Solution:
3,
Rln(PJ ) 1
TI = Y + [
∆HHK! TJ
where:
TB = the normal boiling point, K
R = the ideal gas constant, 8.314 J/K.mol
P0 = the vapor pressure at a given temperature, atm
DHvap = the heat of vaporization of the liquid, J/mol
T0 = the given temperature, K
https://youtu.be/vRS4jdgYjPs
Video 1.13: Evaporation vs boiling (YouTube)
Example 1.17:
Solution:
Example 1.18:
Solution:
Example 1.19:
Solution:
O
HHV (BTU/lb) = 14.544C + 62.028. dH − e + 4050S
8
Where C, H, O and S are the weight fractions of Carbon, Hydrogen, Oxygen and Sulphur
respectively. These are typically obtained from flue gas analyses.
Example 1.20:
Solution:
1
Lowry, HH (ed), 1945. Chemistry of Coal Utilization, Ch. 4, Wiley: New York
Kevin Harding ([email protected]) 19
School of Chemical and Metallurgical Engineering, University of the Witwatersrand, Johannesburg (South Africa)
ORCID - research profile ResearchGate - publications & lecture notes YouTube - lecture videos
Energy balances – Non-Steady State Energy Balances
Example 1.21:
Solution:
Example 1.22:
Solution:
1.11 Psychometry
LB
Lh MN gMh
P = PJ d1 − e ≈ PJ . exp d− e
TJ RTJ
where:
P = Atmospheric pressure (Pa)
P0 = Sea level standard atmospheric pressure (101 325 Pa)
L = Temperature lapse rate (0.0065 K/m)
h = height above sea level (m)
T0 = sea level standard temperature (298.15K)
g = Standard gravity (9.80665 m/s2)
M = molar mass of dry air (0.0289644 kg/mol)
R = universal gas constant (8.31447 J/mol.K)
Therefore:
(P.RJSST).(J.J+RPSVV)
0.0052h (R.W,VVX).(J.JJST) (9.80665). (0.0289644). h
P = 101325 d1 − e ≈ PJ . exp Y− [
298.15 (8.31447). (298.15)
Additional Reading
Balmer, RT, 2011. Modern Engineering Thermodynamics, Academic Press.
Elliott, JR, Lira, CT, 2012. Introductory Chemical Engineering Thermodynamics, 2nd ed.,
Pearson International Edition.
Felder RM, Rousseau RW, 2000. Elementary Principles of Chemical Processes, 3rd ed., Wiley.
Himmelblau, DM, 2012. Basic Principles and Calculations in Chemical Engineering, 8th ed.,
Prentice-Hall International.
Koretsky, MD, 2013. Engineering and Chemical thermodynamics, 2nd ed., Wiley.
Reklaitis, GV, 1983. Introduction to Material and Energy Balances, Wiley.
Rogers, G, Mayhew, Y, 1992. Engineering Thermodynamics – Work and Heat Transfer, 4th
ed., Pearson – Prentice Hall.
Rogers, G, Mayhew, Y, 1995. Thermodynamic and Transport Properties of Fluids (SI Units),
5th ed., Blackwell Publishing.
Sandler, SI, 2006. Chemical, biochemical, and Engineering Thermodynamics, 4th ed., Wiley.
Skogestad, S, 2009. Chemical and Energy Process Engineering, CRC Press, Taylor & Francis
Group.
Problems
$
Example 1.23: Understanding energy - No work done on or by the system b) 𝐸'! ()* = 𝐸'+ ,)(()- = 𝑣,)(()- #
$
A piston (diameter 5 cm; height 10 cm) containing a certain volume of air is placed 2 m Example 1.25: Water from a waterfall Nitrogen is stored in a 10 L flask at an initial pressure of 5.921 atm. Calculate the heat
above the ground. that is transferred when the nitrogen is cooled from 90 to 30°C.
a) Define a suitable reference point. Water falls from a 210 m waterfall. Given:
b) What is the potential energy of the gas in the piston? a) What is the potential energy of the water at the top of the waterfall? Cp [kJ/mol. °C] = 0.0290 + 0.2199 × 10-5T + 0.5723 × 10-8T2 – 2.871 × 10-12T3
c) The piston is moved and placed on the floor. What is the potential energy of the b) Assuming that all the energy is converted to kinetic energy, what is the water’s
gas at its new position? speed at the bottom of the waterfall? Solution:
d) While still on the ground, what is the potential energy of the gas if the piston is c) The terminal velocity of water is given as roughly 10 m/s (depending on the size From an energy balance that assumes no kinetic or potential energy:
0
compressed to halve its original height?* of the droplets). How does this compare to the answer for (b)? 𝑑𝑈
= I𝐻 +𝑄 +𝑊
e) Discuss the importance of defining a reference point. d) Assuming that the water droplet is only travelling at 10 m/s when it reaches the 𝑑𝑡
/1$
bottom of the waterfall, what happened to the remaining energy?
f) How does a system differ from and a reference point?
e) If the temperature of the water at the top of the waterfall is 25°C, calculate the For a closed system with no work added, this becomes:
temperature of the water at the bottom of the waterfall. Assumption: No water 𝑑𝑈
Solution: evaporates. =𝑄
a) The reference point is taken as the floor. 𝑑𝑡
f) If water evaporates, would the temperature in (e) be higher or lower? Explain.
Rearranging and rewriting for the definition of U:
b) Assuming:
- Mass of piston can be neglected Solution: 𝑄 = 𝑛 S 𝐶3 𝑑𝑇
- No heat transfer between system and surroundings a) 𝐸'! ()* = 𝑔ℎ()* = 9.8 𝑚𝑠 "# × 270 𝑚 = 2058 𝑚# 𝑠 "# 𝑜𝑟 𝑁𝑚/𝑘𝑔 𝑜𝑟 𝐽/𝑘𝑔
A stream of water (15 kg/hr) is heated from 25°C to T2 at a constant pressure of 2 bar Temperature
such that the final enthalpy is 335630.9 kJ/hr. Calculate the specific volume of the water
at the new temperature.
Solution:
a) Liquid water
Given:
b) Superheated steam