Advanced Environmental Engineering: Chapter 1. Mass and Energy Transfer
Advanced Environmental Engineering: Chapter 1. Mass and Energy Transfer
Major Theme
1. Introduction
1) American unit verse SI units
2) Law of conservation of mass
- Important tool for quantitatively tracking pollutants
- Pollutant have to go somewhere
(Transport pollutants from one medium to another)
3) Law of conservation of energy
- Global climate change
- Thermal pollution
- Dispersion of air pollution
2. Units of Measurement
- 1 -
Quantity Prefix Symbol
18
10 exa E
1015 peta P
1012 tera T
109 giga G
106 mega M
103 kilo k
102 hecto h
10 deca da
-1
10 deci d
-2
10 centi c
-3
10 milli m
-6
10 micro μ
10-9 nano n
10-12 pico p
- 2 -
× ×
×
×
or more simply
×
℃ and (5)
3. Material Balance
- 3 -
Figure 1. A material balance diagram.
Figure 2. A steady state conservative system. Pollutants enter and leave the
region at the same time rate.
- 4 -
mg/L and the tributary chloride concentration is 40.0 mg/L. Treating chlorides as
a conservative substance, and assuming complete mixing of the two streams, find
downstream chloride concentration.
<Solution>
Figure 3 Flow rate and chloride concentration for example stream and tributary.
which yields
- 5 -
ln ln ln
<Solution>
Input rate = QsCs + QwCw
= (5.0 m3/s×10.0 mg/L+0.5 m3/s×100.0 mg/L)×103 L/m3
= 1.0 × 105 mg/s
- 6 -
Output rate = QmCm = (Qs + Qw) × C, since Cm = C.
= (5.0 + 0.5)m3/s × C mg/L × 103 L/m3
= 5.5 × 103 mg/s
× × × ×
Decay rate = KCV =
×
= 23.1 × 103 C mg/s
So from Eq.(1.15),
1.0 × 105 = 5.5 × 103 C + 23.1 × 103 C = 28.6 × 103 C
×
×
Estimate the steady state concentration of HCHO in the air and how does the
results compare with the threshold for eye irritation of about 0.05 ppm?
Assuming complete mixing, 25℃, 1 atm.
<Solution>
- 7 -
Input rate = 50 smokers × 2 cigas/hr × 1.4 mg/ciga = 140.0 mg/hr
Output rate = 1000 m3/hr × C (mg/m3) = 1000 C mg/hr
Decay rate = KCV = (0.4/hr)×(C mg/m3)×(500 m3) = 200 C mg/hr
Input rate = Output + Decay rate
140 = 1000 C + 200 C = 1200 C
C = 0.117 mg/m3
Convert mg/m3 to ppm by using Eq.(1.7), MW of HCHO is 30,
× ×
This is more than enough to cause eye irritation.
- 8 -
At steady state, dC/dt = 0, which yields
∞ (17)
At unsteady state, from Eq.(1.16)
(18)
Let
(19)
Then
(20)
So, Eq.(1.18) becomes
(21)
(22)
or
(1.24)
- 9 -
Figure 7. Step function response for a complete-mix box model.
- 10 -
[Example 5] Cooling tower improvement
<Given Data>
1) Q5 = 0 (only air flow)
2) C4 = 0 (only pure water is evaporating)
3) Q1, Q̂ 6, C1, W5
<Question> Q4 and C2 = ?
<Solution>
▶ Overall Mass Balance : Q4 = Q̂ 1 – Q̂ 6 ①
▶ Component Balance : Accumulation = Input – Output
VdC/dt = (Q̂ 1C1 + W5) – (Q4C4 + Q6Ĉ 6) = 0 at steady state
and C4 = 0
∴ Q6C6 = Q̂ 1C1 + W5 → Ĉ 2 = Ĉ 6 = (Q̂ 1C1 + W5) / Q̂ 6 ②
- 11 -
<Given Data>
1) Q5 = 0 (only air flow)
2) C4 = 0 (only pure water is evaporating)
3) Q1, Q6, C1, W5
3) Recovery rate of filter = R = Qb/Qa → Qb = RQa ③
4) Rejection rate of filter = η = (Ca – Cb)/Ca = 1 – Cb/Ca
→ Cb = (1 – η)Ca = (1 – η)C2 ④
<Question> Q4 and C2 = ?
<Solution>
▶ Mass Balance :
Q1 = Q4 + Q6 → Q6 = Q1 – Q4 ⑤
Q6 = Qa – Qb ⑥
Eq③→Eq⑥ : Q6 = (1 – R)Qa → Qa = Q6 / (1 – R) ⑦
▶ Component Balance : Accumulation = Input – Output
V∙dC/dt = (Q1C1 + W5) – (Q4C4 + Q6C6) = 0 at steady state
- 12 -
and C4 = 0
Q6C6 = Q1C1 + W5 ⑧
Q6C6 = QaCa – QbCb = QaC2 – QbCb ⑨
Eq④→Eq⑨ : Q6C6 = QaC2 – (1 - η)QbC2 ⑩
Eq⑧→Eq⑩ : Q1C1 + W5 = QaC2 – (1 – η)QbC2 ⑪
Eq③→Eq⑪ : Q1C1 + W5 = QaC2 – (1 – η)RQaC2
= QaC2[1 – (1 – η)R] ⑫
Eq⑦→Eq⑫ : Q1C1 + W5 = Q6C2[1 – (1 – η)R] / (1 – R)
Let ζ = R / (1 – R) → Q1C1 + W5 = Q6C2(1 + ζη) ⑬
Rearrange Eq⑬ → C2 = (Q1C1 + W5) / Q6(1 + ζη) ⑭
≪ 1.0
C2 after improvement C2 1
= =
C2 before improvement Ĉ 2 1 + ζη
Q̂ 1(Ĉ 2 – C) 1 = W5 + Q4Ĉ 2
Q̂ 1 =
W5 + Q4Ĉ 2
⑮
Ĉ 2 –C 1
- 13 -
From Eq⑭,
–
W5 + Q1C1
C2 =
(1 + ζη)(Q1 Q4)
–
(1 + ζη)(Q1 Q4)C2 = W5 + Q1C1
Q1[(1 + ζη)C2 C1]– = W5 + (1 + ζη)Q4C2
Q1 =
W5 + (1 + ζη)Q4C2
(1 + ζη)C2 C1 – ⑯
Divide Eq ⑯ by ⑮, Q1 (Ĉ 2–C )[W
1 5 + (1 + ζη)Q4C2]
Q̂ 1
=
[(1 +ζη)C2 –C ](W1 5 + Q4Ĉ 2)
- 14 -
[Extra Example 6] RO system with feed-and-bleed mode operation
- 15 -
= (1 – η)(Q1C1 + C1Q1Q5 / Q6 – C3Q3Q5 / Q6)
= (1 – η)Q1C1 + (1 – η)C1Q1Q5 / Q6 – (1 – η)C3Q3Q5 / Q6
Multiply Q6 on both side and rearrange, then
Q2Q6C3 = (1 – η)Q1Q6C1 + (1 – η)Q1Q5C1 – (1 – η)Q3Q5C3
[Q2Q6 + (1 – η)Q3Q5]C3 = (1 – η)(Q5 + Q6)Q1C1 = (1 – η)Q1Q4C1
(1 – η)Q1Q4
C3 = × C1
Q2Q6 + (1 – η)Q3Q5
(1 – η)Q1(1 / R – 1)Q3
= × C1
(Q3/R)(Q1 – Q3) + (1 – η)Q3(Q3/R – Q1)
(1 – η)Q1(1 – R)Q3
= × C1
Q3(Q1 – Q3) + (1 – η)Q3(Q3 – RQ1)
(1 – η)(1 – R)Q1Q3
= × C1
Q3[(Q1 – Q3) + (1 – η)(Q3 – RQ1)]
(1 – η)(1 – R)Q1
= × C1
(Q1 – Q3) + (1 – η)(Q3 – RQ1)
C3 =
(1 – η)(1 – R)Q1
[1 – (1 – η)R]Q1 – ηQ3
× C1 ∥Ans ⑪
From Eq ②, C2 = C3 / (1 – η) ⑫
From Eq ⑪ and ⑫,
C2 =
(1 – R)Q1
[1 – (1 – η)R]Q1 – ηQ3
× C1 ∥Ans ⑬
From Eq ⑨ and ⑬,
- 16 -
(1 – R)Q1Q3
C6 = [Q1C1 – ] / Q6
[1 – (1 – η)R]Q1 – ηQ3
(1 – R)Q1Q3
= [Q1C1 – ] / (Q1 – Q3)
[1 – (1 – η)R]Q1 – ηQ3
=
Q1C1
Q1 – Q3
–
(1 – R)Q1Q3
(Q1 – Q3){[1 – (1 – η)R]Q1 – ηQ3}
∥Ans
- 17 -
4 Energy Fundamentals
- 18 -
atoms within molecules, and within atoms.
∴ E(Total) = U(Internal) + KE(Kinetic) + PE(Potential) (47)
▶ Specific Heat(c)
The amount of energy needed to raise the temperature of a unit mass of a
substance by 1 degree.
▶ British Thermal Unit(Btu)
the energy required to raise 1 lb of water by 1 ℉.
▶ Kilocalorie(kcal)
the energy required to raise 1 kg of water by 1 ℃
※ Conversion Factor : 1 kcal/kg℃ = 1 BTU/lb℉ = 4.184 kJ/kg℃
▶ Specific heat at constant volume(c ) & Specific heat at constant pressure(c )
v p
- 19 -
【Example 10】A Water Heater
∙Energy source : Water heater 5 kW
∙Boundary volume : 40 gallon
<Question> Time to raise the temperature 50 ℉ → 140 ℉ ?
<Assumption>
1) No energy loss during converting electrical energy to heat.
2) No energy losses from the tank.
3) No water withdrawn from tank.
<Solution>
∙ Rate of energy input = Power = 5 kW × Δt(hr) = 5∙Δt (kWhr)
∙ Energy output = 0
∙ Change in stored energy
=40gal×8.34lb/gal×1Btu/lb°F×(140-50)°F=30,060BTU=30×103BTU
※ specific gravity of water = 8.34 lb/gal
∙ By conservation, Input energy rate = Internal energy change rate,
5∙Δt (kWhr) × 3,412 ((BTU/hr)/kW) = 30 × 103 BTU
∴ Δt = 1.76 hr∥Answer
■ Phase Change
- 20 -
Energy released or absorbed in phase change = m L (52)
where m = mass, L = latent heat of fusion or vaporization.
※ Melting energy of Ice at 0 ℃ = 333 kJ/kg
※ Vaporization energy of water at 100 ℃ = 2,256 kJ/kg
※ 15℃ is a useful number that can be used to estimate the amount of energy
required to cause surface water on the Earth to evaporate. The value of 15℃
has been picked as the starting temperature since that is approximately the
current average surface temperature of the globe.
- 21 -
Figure 13 The latent heats of fusion and vaporization for water.
<Question>
1) Energy required to cause that much water to evaporate each year.
2) Compare this to world energy consumption at 2007(4.7×1017 kJ).
3) Compare it to the average rate at which sunlight strikes the surface of the
Earth(168 W/m2).
<Solution>
In Table 4, Heat of vaporization at 15 ℃ = 2,465 kJ/kg
The total energy required to vaporize all of that water is
1) Increase in stored energy
= 1 m/yr × 5.10 × 1014 m2 × 103 kg/m3 × 2,465 kJ/kg
= 1.25 × 1021 kJ/yr ∥ Answer
- 22 -
2) 4.7×1017 kJ ÷ 1.25 × 1021 kJ ≃ 2,700 Times∥ Answer
3) From Answer 1), the rate of sunlight striking on the Earth is
1.25×1021 kJ 1000 J yr day hr 1 W
= 78 W/m2
yr kJ 365 day 24 hr 3600 sec 1 J/s 5.10×1014 m2
- 23 -
■ Open System
Rate of change of stored energy = ṁ × c × ΔT
where ṁ =Q×d : the mass flow rate across the system boundary.
<Question>
1) ΔT of cooling water = 10 ℃, Q of cooling water?
2) Temperature of river water after receiving the heated cooling water?
<Solution>
1) Rate of change in internal energy = ṁ × c × ΔT
1,700 MWt = ṁ (kg/s) × 4,184(J/kg℃) × 10(℃) × 1(MWt/106 J/s)
ṁ = 40.6 × 103 kg/s = 40.6 m3/s ∥Answer
2) From Rate of change in internal energy = ṁ × c × ΔT
- 24 -
Figure 14 Cooling water energy balance for the 33.3 percent efficient, 1000
MWe power plant in Example 12.
- 25 -
2. The 2nd Law of Thermodynamics
The second law of thermodynamics says that there will always
be some waste heat; that is, it is impossible to devise a machine that can
convert heat to work with 100 percent efficiency
max
FIGURE 15 Definition of terms for a Carnot engine.
where T = absolute temperature,
K = ℃ + 273.15 (56)
R = ℉ + 459.67 (57)
- 26 -
FIGURE 16. A fuel-fired, steam-electric power plant.
- 27 -
Electricity
↑
Steam Turbine
Electricity
(Rankin Cycle)
↑ ↑
Combustion
Gas
→ Gas Turbine
(Brayton Cycle)
→ Steam
Generator
→ Off gas
to stack
- 28 -
∙ emission
Sulfur
=
130 g S
106 kJ input
10,800 kJ input
∙
kWe hr
= 1.4
g S
kWe hr ∙
∙ emission
SO 2
=
260 g SO2
106 kJ input
10,800 kJ input
∙
kWe hr
= 2.8
g SO2
∙
kWe hr
∙ Particulate
emission
=
13 g Part.
106 kJ input
10,800 kJ input
∙
kWe hr
= 0.14
g Part.
∙
kWe hr
∙ Coal
input
=
10,800 kJ input
∙
kWe hr
1 g coal
24 kJ
= 450
g coal
kWe hr ∙
∙ Carbon
emission
=
0.62 g carbon
g coal
450 g coal
kWe hr∙ = 280
g carbon
kWe hr∙
2) Removal efficiency to meet limitation
∙S % of coal = 2%
∙ Sefficiency
removal
= 1 -
1.4 g S by limitation
9.0 g S by raw emission
= 0.85 = 85%
∙ Fly ash
generated
= 10% × 70% ×
450 g coal
kWe hr∙ = 31.5
g fly ash
kWe hr ∙
∙ Part. removal
efficiency
= 1 -
0.14 g Part. by limitation
31.5 g S by raw emission
= 0.885 = 99.5%
- 29 -
FIGURE 17. Energy and mass balance for a coal-fired power plant generating 1
∙
kWe hr of electricity (see Example 13).
∙T = 5 ℃
c
<Question> ηmax?
5 + 273
<Solution> ηmax = 1 - = 0.08 = 8%
30 + 273
※ Real efficiency in real OTEC plant : 2∼3%
- 30 -
3. Conductive and Convective Heat Transfer
▶ Heat will be transferred from the hotter object to the colder one.
▶ Conductive heat transfer is usually associated with solids, as one molecule
vibrates the next in the lattice.
▶ Convective heat transfer occurs when a fluid at one temperature comes in
contact with a substance at another temperature.
▶ Radiation can take place even in the absence of any physical medium
between the objects.
- 31 -
℃
R = overall thermal resistance (m2- /W) or (hr-ft2- /BTU) ℉
※ The overall thermal resistance R is called the R-value.
【Example 15】Reducing Pollution by Adding Ceiling Insulation
∙Heating season = 8 month per year = 8 × 30 × 24 = 5,760 hr
∙T = 40 ℉
o
∙T = 70 ℉
i
<Question>
1) Annual money saving? and pay back period?
2) Annual reduction in particulate, and carbon emissions?
Basis : 1 million homes, coal plants in Example 13
<Solution>
1) Annual saving and pay back period
∙ Heat loss rate with the existing insulation
A(Ti - To) 1500 ft2 × (70-40) ℉ BTU
q =
R
=
℉
11 (hr-ft2- /BTU)
= 4090
hr
Energy
=
(4090 –1125)BTU kWe-hr 5,760 hr
= 5005
kWe-hr
saved hr 3412 BTU yr yr
∙ Annual saving
Annual
Saving
=
5005 kWe-hr
yr
8 ¢
kWe-hr
= 40,040
¢
yr
≃ 400
$
yr
- 32 -
∙ Saving energy for 1 million same home = 5×109 kWe-hr/yr
∙ Carbon
reduction
=
280 g C
∙
kWe hr
5×109 kWe-hr
yr
MT
106 g
= 1400×103
MT
yr
∙ SO2
reduction
=
2.8 g SO2
∙
kWe hr
5×109 kWe-hr
yr
MT
106 g
= 14×103
MT
yr
∙ Particulate
reduction
=
0.14 g Part.
∙
kWe hr
5×109 kWe-hr
yr
MT
106 g
= 0.7×103
MT
yr
- 33 -