Professional Documents
Culture Documents
Review PDF
Review PDF
MECH3429
AIR CONDITIONING AND REFRIGERATION:
Yuguo Li
OVERVIEW
Try to understand this one:
OF 56 5`
A lot of practical tips:
Building services industry in Hong Kong
http://ibse.hk/about-bse.htm Dr Sam Hui
Government: ArchSD, BD, DSD, EMSD, EPD, FSD, Housing Dept, WSD
Utilities: CLP Power, HK Electric, Towngas, MTRC
Professional Societies: HKIE (BS Div), CIBSE HK Branch, ASHRAE HK Chapter
Hong Kong codes: CAP123 Building Ordinance &Regulations, COP from ASD, BSD, EPD,
FSD, EMSD and power companies, standards from BS, ASHRAE etc
Partners: Clients, developers, building owners, architects (Authorised Persons), civil
engineers, structural engineers, quantity surveyors, urban designers/planners, main
contractors (builders), sub-contractors/suppliers
Jobs: Consulting/design engineers, contracting engineers, BS coordinator (main/civil
contracting firms), project manager/engineer (for developers), government engineer, sale
engineer, facility management or maintenance engineer
Consulting firms: Arup, Wong& Ouyang, WSP (formally WSP/Parsons Brinckerhoff), J.
Roger Preston Limited, Meinhardt etc.
The world is big – go to Mainland, South East Asia, elsewhere.
Engineers can also do research – Doing master in Europe and USA, why not PhD here in
HKU and elsewhere (you can get scholarship, e.g. HK$17k or HK$25k in HK). Innovation is
future.
What are the options to keep a space cool?
1. Place ice in the room
2. Open windows if outside is cool
3. Have a fan? Desk fan or even a simple fan?
4. Buy an air-conditioner? But how does an air conditioner work?
𝑘
Table 2. Properties of water at atmospheric pressure 𝛼= 𝑃𝑟 =
𝜌𝑐𝑝 𝛼
T (kg/m3) cp (kJ/kg.K) (kg/s.m) k (W/m.K) ℎ𝑓𝑔
(oC) (cm2/s) (cm2/s) (kJ/kg)
Air loop
Secondary
components
Primary
components
Chapter 1 - intro
33 oC, 28 oC Control Loop
33 oC 27 oC
12 oC 10 oC
25 oC
12 oC 49 oC 35 oC
35 oC
7 oC 38 oC 29 oC
7 oC 3 oC 29 oC
13 oC 13 oC
𝑘𝐴 𝑇0 −𝑇𝐿 𝐿
We write 𝑄 = 𝑇0 − 𝑇𝐿 = , where 𝑅𝑡ℎ = 𝑘𝐴 is the equivalent thermal
𝐿 𝑅𝑡ℎ
resistance, or 𝑇0 − 𝑇𝐿 = 𝑄𝑅𝑡ℎ
Calculate the rate of heat flow through the wall, and the temperature of the
interface between the two layers making the wall.
Calculate the rate of heat flow through the wall, and the temperature of the
interface between the two layers making the wall.
Calculate the rate of heat flow into the room, and the temperature of the outer
surface of the wall.
Solution:
Consider a unit area of the wall.
1 1 1 1 𝐿
𝑅𝑐𝑖 = 𝐴ℎ = 1×10 = 0.1 K/W; 𝑅𝑐𝑜 = 𝐴ℎ = 1×25 = 0.04 K/W; 𝑅𝑤 = 𝐴𝑘 =
0.08
1×1.6
= 0.05 K/W
Substitute 𝑄𝑜 and 𝑄𝑖 into the energy balance equation, gives 𝑇𝑠𝑜 = 37.37 oC.
Then we can find that 𝑄𝑜 = −184.2 W/m2, suggesting that some of radiation
absorbed by the wall is lost to the ambient air by convection.
Example 6. A concrete wall of thickness 8 cm and thermal conductivity 1.6
W/mK, absorbs solar radiation at a steady rate of 300 W/m2. The heat transfer
coefficient between the outer surface of the wall and the ambient at 30 oC is 25
W/m2K. The heat transfer coefficient between the room air at 20 oC and the
inner surface of the wall is 10 W/m2K.
Calculate the rate of heat flow into the room, and the temperature of the outer
surface of the wall.
Solution:
Consider a unit area of the wall.
1 1 1 1 𝐿
𝑅𝑐𝑖 = 𝐴ℎ = 1×10 = 0.1 K/W; 𝑅𝑐𝑜 = 𝐴ℎ = 1×25 = 0.04 K/W; 𝑅𝑤 = 𝐴𝑘 =
0.08
1×1.6
= 0.05 K/W
Substitute 𝑄𝑜 and 𝑄𝑖 into the energy balance equation, gives 𝑇𝑠𝑜 = 37.37 oC.
Then we can find that 𝑄𝑜 = −184.2 W/m2, suggesting that some of radiation
absorbed by the wall is lost to the ambient air by convection.
Water vapor’s concentration in low
Dry air is the air that is free of any moisture.
The water vapor is generally present in moist air
in a superheated state.
The pressure of water vapor in typical ambient
air is relatively low. For example, at 30oC the
water vapor pressure in air that is fully
saturated with water vapor is about 4.24 kPa,
which is the saturated pressure of water vapor
at 30oC from the steam table.
Tibet
Chapter 2 Figure 1. The saturation vapor pressure as a function of its temperature.
T (oC) Saturation pressure Enthalpy of saturated Enthalpy of vaporization Enthalpy of
Pg (kPa) liquid (kJ/kg) (kJ/kg) saturated vapor, ℎ𝑔
(kJ/kg)
0.010 0.6113 0.0 2501.3 2501.3
2 0.7056 8.4 2496.6 2505.0
5 0.8721 21.0 2489.5 2510.5
10 1.228 42.0 2477.7 2519.7
15 1.705 63.0 2465.9 2528.9
20 2.338 83.9 2454.2 2538.1
25 3.169 104.9 2442.3 2547.2
30 4.246 125.8 2430.4 2556.2
35 5.628 146.7 2418.6 2565.3
40 7.383 167.5 2406.8 2574.3
45 9.593 188.4 2394.8 2583.2
50 12.35 209.3 2382.8 2592.1
55 15.76 230.2 2370.7 2600.9
60 19.94 251.1 2358.5 2609.6
65 25.03 272.0 2346.2 2618.2
70 31.19 293.0 2333.8 2626.8
75 38.58 313.9 2321.4 2635.3
80 47.39 334.9 2308.8 2643.7
85 57.83 355.9 2296.0 2651.9
90 70.13 376.9 2283.2 2660.1
95 84.55 397.9 2270.2 2668.1
𝑃𝑡𝑜𝑡𝑎𝑙 = 𝑃𝑖
𝑖=1
The humidity ratio, 𝜔 is defined as the mass of water vapor in a given volume of
mixture to the mass of dry air in the same volume.
∆𝑚 ∆𝑚 /∆𝑉 𝜌
𝜔 = ∆𝑚 𝑣 = ∆𝑚 𝑣//∆𝑉 = 𝜌 𝑣 ;
𝑎 𝑎 𝑎
where ∆𝑚𝑣 and ∆𝑚𝑎 the mases of vapor and dry air respectively in a volume of ∆𝑉.
𝜌𝑣 and 𝜌𝑎 the densities of vapor and dry air respectively.
Both 𝑐𝑎𝑚 and ℎ𝑔𝑚 are evaluated at mean temperature and mean humidity ratio of the
two states.
𝜔 2 +𝜔 1 𝑡 2 +𝑡 1
𝑐𝑎𝑚 = 𝑐𝑝𝑎 + 𝜔𝑐𝑝𝑤 = 𝑐𝑝𝑎 + 𝑐𝑝𝑤 ; and ℎ𝑔𝑚 = ℎ𝑔𝑜 + 𝑡𝑐𝑝𝑤 = ℎ𝑔𝑜 + 𝑐𝑝𝑤
2 2
Specific volume of moist air 𝑣 (m3/kg dry air) can be obtained by using the ideal gas
equation for dry air
𝑅𝑎 𝑇
𝑣= , we can also obtain
𝑃𝑎
𝑃𝑣 𝑅𝑎 0.622 𝑃 − 𝑃𝑎
𝜔= = ;
𝑃 − 𝑃𝑣 𝑅𝑣 𝑃𝑎
𝑅 𝑇 𝜔 𝑅 𝑇 𝜔
𝑣 = 𝑃𝑎 1 + 0.622 , so 𝑣 = 𝑃𝑎 1 + 0.622 .
http://www.flycarpet.net/en/PsyOnline
Solution
From the steam table, the saturation vapor pressure at 25oC is 3.166 kPa.
Solution
From the steam table, the saturation vapor pressure at 25oC is 3.169 kPa.
Solution
Assume dew point is 𝑡𝑑𝑝 =6oC, from the steam table, the saturation
pressure at this temperature is 0.9346 kPa.
Solution
Assume dew point is 𝑡𝑑𝑝 =6oC, from the steam table, the saturation pressure at
this temperature is 0.9433 kPa.
The saturation vapor pressure at the air dry-bulb temperature of 22oC is 2.670 kPa.
𝑃𝑣 0.9433
At 𝑡𝑑𝑏 = 22 oC, 𝑡𝑑𝑝 =6oC, = = = 35.33%, which is the max
𝑃𝑠 2.670
allowance RH.
Psychrometric processes for heating and air conditioning
Mixing of two moist air streams
The two streams 1 and 2 are mixed to become a joint stream 3. The change in
enthalpy is much greater than those in kinetic energy and potential energy, which
are ignored here. We also assume adiabatic process at steady state conditions and
constant pressure.
Mass balance of air: 𝑚𝑎1 + 𝑚𝑎2 = 𝑚𝑎3 , All there symbols
represent respectively dry air mass flow rates of the three streams.
Mass balance of water: 𝑚𝑎1 𝜔1 + 𝑚𝑎2 𝜔2 = 𝑚𝑎3 𝜔3
Energy balance of moist air: 𝑚𝑎1 ℎ1 + 𝑚𝑎2 ℎ2 = 𝑚𝑎3 ℎ3
𝑚𝑎1 ℎ 3 −ℎ 2 𝜔 3 −𝜔 2
Thus, = = , suggesting that 1-2-3 is a line.
𝑚𝑎2 ℎ 1 −ℎ 3 𝜔 1 −𝜔 3
If we assume that the specific heat capacities of the moist air are equal in the three
streams, then we can obtain he following appropriate relation for the mixed
temperature
𝑚𝑎1 𝑚𝑎2
𝑡3 = 𝑡1 + 𝑡 ;
𝑚𝑎3 𝑚𝑎3 2
Sensible heating and cooling
Heating may be obtained by an electrical resistance heater or a tubular coil
carrying a hot fluid. The air can be cooled by passing a cold fluid like chilled
water or a refrigerant through a similar coil.
During sensible heating and cooling, the moisture content of the air remains
constant. This means that we need to ensure that the coil surface temperature is
above the dew-point of the entering air.
Hence 𝑄𝑖𝑛 = 𝑚𝑎 ℎ2 − ℎ1
Assume that the specific heat capacity of moist air is constant, 𝑄𝑖𝑛 =
𝑚𝑎 𝑐𝑝𝑚 𝑡2 − 𝑡1
The porous structure distributes water in the form of a thin film, and
increases the contact area between air and water to facilitate evaporation.
Under ideal conditions, the state of air during adiabatic evaporative cooling follows a constant
wet-bulb temperature in the psychrometric chart. However, in practice, the moist air leaving
the cooler at 2 is not cooled to the wet-bulb temperature, which is the lowest possible
temperature to which the air could be cooled.
Humidification of air
Assume that all the moisture sprayed is retained in the air stream. Consider the control
volume as shown.
Consider the change of enthalpy of the two points 1 and 2 in case the dry-bulb
temperature is unchanged.
ℎ2 𝑡 − ℎ1 𝑡 = 𝜔2 − 𝜔1 ℎ𝑔 𝑡 ;
ℎ −ℎ
hence if 𝑞 = 𝜔 2 −𝜔1 = ℎ𝑤 = ℎ𝑔 𝑡 , i.e. the line 1 − 2 is the constant dry-bulb
2 1
temperature line passing through 1.
The tube surface temperature of the actually cooling coil is not uniform over
different tube rows. This results in the non-linear variation of the mean air
temperature, i.e. curve 1-2.
𝑚 𝑎𝑏
A bypass factor is defined 𝑏 = .
𝑚𝑎
For the assumed adiabatic mixing process, we obtain
𝑚 𝑚
𝑡2 = 𝑚𝑎𝑏 𝑡1 + 𝑚𝑎𝑐 𝑡𝑑 ; (note: 𝑡𝑑 equals to 𝑡4 in the figure).
𝑎 𝑎
𝑡 −𝑡
Using the above three equations, we have 𝑏 = 𝑡2 −𝑡𝑑 .
1 𝑑
Where SHR is the sensible heat ratio, which may be obtained directly from the
protractor in the psychrometric chart for the process 1-4 as shown for S2.
Space condition line
A room is a special space where heat/mass transfer processes occur
It is customary to list all energy and moisture flows into the space which the air
conditioning system needs to remove in two groups, i.e. sensible cooling load
and the latent cooling load.
The sensible cooling load (total sensible heat gain) 𝑄𝑠𝑡 : heat flow through
walls, roofs, and windows; appliances generating heat, lighting and
occupants within the space; unintended air infiltration.
The latent cooling load 𝑄𝑙𝑡 : the moisture released by appliances and
occupants within the space, and moist air infiltration. The total rate of
moisture gain by the air is 𝑚𝑤𝑡 .
Consider a control volume with a single inlet port 1 and a single outlet port 2,
surrounding the space.
ℎ −ℎ 𝑄𝑠𝑡 +𝑄𝑙𝑡
𝑞 = 𝜔 2 −𝜔1 = 𝑚 𝑤𝑡
;
2 1
Thus, for fixed values of the sensible and latent cooling loads, and the rate of
moisture gain, the enthalpy-moisture ratio, 𝑞 is constant.
Moreover, the state point of supply air at 1 on the psychrometric chart must lie on
∆ℎ
a straight line drawn through 2, parallel to the direction of 𝑞 in the ∆𝜔 -protractor.
Using the earlier results, the sensible cooling load 𝑄𝑠𝑡 and the latent cooling load 𝑄𝑙𝑡 may
be expressed by
𝑄𝑠𝑡 = 𝑚𝑎 𝑐𝑎𝑚 𝑡2 − 𝑡1 and 𝑄𝑙𝑡 = 𝑚𝑎 ℎ𝑔𝑚 𝜔2 − 𝜔1
𝑄𝑡 = 𝑄𝑠𝑡 + 𝑄𝑙𝑡 .
Again, for the temperature range from 0 oC to 60oC of a typical psychrometric chart, the
following mean values may be used. 𝑐𝑎𝑚 = 1.02 kJ/kgK and ℎ𝑔𝑚 = 2555 kJ/kg.
𝑄𝑠𝑡
The sensible heat ratio (SHR) is SHR= .
𝑄𝑡
Application of the psychrometric processes
Single-zone air conditioning system: summer air conditioning systems
A filter, a cooling and dehumidifying coil, and a fan are arranged in series
and are connected using ducts through which air flows.
For ventilation, a portion of the air withdrawn from the space at 2 is
discharged to the ambient at 3, and replaced with an equal mass of
outdoor air at 4. Again, we assume steady state condition.
Point 2 and condition line of the space– this is the space condition in term
of the temperature and humidity ratio or other moist air properties. If the
sensible and latent cooling loads of the space are estimated, then SHR
can be calculated. The condition line of the space is obtained by drawing
a straight line through Pont 2 in the direction of the line on the protractor
pointing towards the SHR-value for the space.
Point 4 – the state of outdoor air is easily located by knowing two
thermodynamic properties of the ambient air. The dry air mass fraction of
the outdoor air to compensate for the return air discharged at 3 is usually
specified.
Point 5 – the mixture of the outdoor air of 4 and the return air of 2. This
uses the adiabatic mixing analysis early, and 5 can be obtained by
dividing the straight line 2-4 in the inverse ratio of the dry air mass flow
rates of the two air streams.
The curved line 5-1, or the coil condition line gives the state of air as it flows
through the cooling and dehumidifying coil. This can be obtained using the
coil manufacturer’s data or the computer-based heat and mass transfer
analyses of the cooling coil.
The air undergoes a slight increase in temperature as it flows through the fan, which
for all practical purposes, could be neglected. However, we could add the energy
input to the fan as a sensible cooling load to the space. In the basic air conditioning
system, the cooling and dehumidifying coil is the sole air processing unit, and
therefore only one property of the air could be controlled. In most systems, this
property would be the dry-bulb temperature.
Example 6. In a summer air conditioning system, 25% of the return air from a space at 30 oC dry-bulb temperature, and
22 oC web-bulb temperature is exhausted and an equal quantity of outdoor air at 34 oC dry-bulb temperature, and 28 oC
web-bulb temperature is mixed with the remaining return air.
The mixture passes over a cooling coil whose coil surface temperature (apparatus dew-point) is 8 oC, and the bypass
factor is 0.25.
The mass flow rate of dry air to the space is 0.8 kg/s; and the pressure is constant at 101.3 kPa.
Calculate (a) dry-bulb temperature and the relative humidity of the air leaving the cooling coil; (b) enthalpy-humidity
ratio of the cooling process, and (c) the refrigeration capacity of the cooling coil.
Example 6. In a summer air conditioning system, 25% of the return air from a space at 30 oC dry-bulb temperature, and
22 oC web-bulb temperature is exhausted and an equal quantity of outdoor air at 34 oC dry-bulb temperature, and 28 oC
web-bulb temperature is mixed with the remaining return air.
The mixture passes over a cooling coil whose coil surface temperature (apparatus dew-point) is 8 oC, and the bypass
factor is 0.25.
The mass flow rate of dry air to the space is 0.8 kg/s; and the pressure is constant at 101.3 kPa.
Calculate (a) dry-bulb temperature and the relative humidity of the air leaving the cooling coil; (b) enthalpy-humidity
ratio of the cooling process, and (c) the refrigeration capacity of the cooling coil.
Solution
Identify state 1 (space air) and state 2 (outdoor air), and divide line 1-2 such that (length 1-3)/(length 2-3) = 1/3. We
have 𝑡𝑑𝑏3 = 31℃ and 𝑡𝑤𝑏3 = 23.5℃.
For the cooling coil process, the bypass factor is 0.25, meaning only 75% of the air passes through the ideal straight
line (3-4), to finally achieve the apparatus dew-point of 8oC, i.e. state 4. The cooled air at state 4 mixes with the
bypass air at state 3, to produce the supply air at state 5. To locate 5, we divide line 3-4 so that (length 4-5)/(length
3-5) =1/3.
Solution
A schematic diagram of the system is shown in Figure 10a. We locate state 2 of the space air (24oC db-temp and 50%RH).
(a) Read state 1 properties on the psychrometric chart. wb-temp = 12.8oC, and RH = 78%.`
23
The mass flow rate of dry air is obtained by 𝑄𝑠 = 𝑚𝑎 𝑐𝑎𝑚 𝑡2 − 𝑡1 , hence 𝑚𝑎 = = 2.505 kg/s.
1.02 24−15
(a) The refrigeration capacity of the cooling coil is given by
𝑄𝑟𝑒𝑓 = 𝑚𝑎 ℎ2 − ℎ1 = 2.505 48 − 36.3 = 29.3 kW, where the enthalpies are obtained from the chart. The total
heat load on the space is 23+6.13= 29.13 kW. This value should equal to 𝑄𝑟𝑒𝑓 by energy conservation.
(a) The apparatus dew-point or the ideal coil surface temperature is the temperature at the point of intersection 3, of the line
1-2 and the saturation line. We read 𝑡𝑑 = 𝑡3 = 10.5oC.
The bypass factor of the coil is
𝑡1 −𝑡𝑑 15−10.5
𝑏= = = 0.33.
𝑡2 −𝑡𝑑 24−10.5
Steady heat transfer through multi-layered
structures
In the isothermal plane method, we assume that the heat flow paths
are parallel only through the insulation and the studs.
Substituting the density and specific heat capacity of air at standard conditions
of 1 bar and 15oC,
𝑄𝑠𝑠[kW]
𝑄 [kW]==1.23𝑉
1.23𝑖𝑛𝑉𝑖𝑛
𝑡𝑜 −𝑡𝑜𝑡𝑖−; 𝑄
𝑡𝑖𝑙 [kW]
; 𝑄𝑙 [kW] = 3020
= 3010𝑉𝑖𝑛 𝜔𝑜 𝑉
−𝑖𝑛𝜔𝑖𝜔𝑜 − 𝜔𝑖
50
Example 1. An exterior wall of a building has a 10cm thick layer of face brick on the outside, followed by a layer of 20
cm thick concrete. A 15 cm thick layer of mineral wool insulation is sandwiched between the concrete and a layer of
plywood of thickness 10 mm. The wall is 10 m long and 3 m high. The outside and inside heat transfer coefficients are 30
W/m2K and 9 W/m2K respectively. The outside and inside air temperature are 22 oC and -15oC respectively. Calculate
(a) The total thermal resistance;
(b) The overall heat transfer coefficient, and Layer, i 𝑘𝑖 ,W/mK 𝐿𝑖 , mm 𝑅𝑖 , m2K/W
(c) The total heat transfer rate through the wall. Brick 0.81 100 0.123
Concrete 1.8 200 0.11
Mineral wool 0.035 150 4.29
Plywood 0.095 10 0.105
Solution: Consider unit area of the wall where the constituent layers are in series.
1
The thermal resistance of a layer is 𝑅𝑖 = 𝐴 𝐾 . Here we obtain the data as follows
𝑖 𝑖
(note: we obtain the thermal conductivity values from ASHRA 2013
Fundamentals).
1 1 1 1
Also 𝑅𝑖𝑛 = = = 0.11 m2K/W; 𝑅𝑜𝑢 = = = 0.033 m2K/W
ℎ 𝑖𝑛 9 ℎ 𝑜𝑢 30
Total thermal resistance 𝑅𝑡𝑜𝑡 = 0.033 + 0.123 + 0.11 + 4.29 + 0.105 + 0.11 =
4.77 m2K/W;
1 1
The overall heat transfer coefficient is 𝑈𝑜 = = = 0.21W/ m2K.
𝑅𝑡𝑜𝑡 4.77
The total heat transfer rate through the wall is 𝑄𝑡𝑜𝑡 = 𝐴𝑈𝑜 𝑇𝑖 − 𝑇𝑜 = 10 × 3 ×
51
0.21 × 22 + 15 = 233W.
Example 2. A building has an exterior wall with 3.8 cm × 14 cm framing that makes 20% of its
area. The wall consists of the following layers of materials:
9.5 cm thick gypsum board on the inside of the framing, fiber glass in the spaces between the
framing, 12.7 mm thick sheathing next to the insulation and framing, 38 cm thick layer of
expanded polystyrene, 100 cm thick layer of brick on the outside.
The air temperature inside is 20oC and the outside ambient temperature is -10oC. The wind speed
is 12 km/h. The inside heat transfer coefficient is 8 W/m2K. Calculate (a) the total thermal
resistance; (b) the average heat transfer coefficient, and (c) the average heat transfer rate through
the wall using (i) the parallel path method, and (ii) the isothermal plane method.
52
Solution
𝐿
It is convenient to first compute the unit thermal resistance 𝑅𝑖 = 𝑘𝑖 , considering the unit
𝑖
cross sectional area of each heat transfer path.
For wind speed 13 km/h or 3.6 m/s, we can find from reference that ℎ𝑜 = 25 W/m2K.
The unit thermal resistance of the indoor and outdoor air films, including contributions
due to thermal radiation are
1 1 1 1
𝑅𝑖 = ℎ = 8 = 0.125 m2K/W; 𝑅𝑜 = ℎ = 25 = 0.04 m2K/W
𝑖 𝑜
The space loads are time varying. The detailed energy and moisture flows in a
building are shown in Figure 3. Clearly, a detailed analysis of all the transient
process is complex and also tedious. In real design, load estimation is done by
computer software. Our aim here is to describe the physical principles that form the
basis of these computerized procedures of load estimation.
Note that ASHRAE simplified the wind heating load estimation, by ignoring the
solar radiation gain.
55
Outdoor design conditions
The design is normally done to meet the first the extreme weather conditions expected at a
location. The weather records for a large number of locations around the world have been
compiled and processed to develop a series of weather related outdoor design conditions. They
are tabulated in ASHRAE Fundamentals.
Table 1. Design weather data for winter heating systems
Location Lat deg. Long. deg DB (oC) DB (oC) 99% Wind speed Wind speed
99.6% (m/s) 1% (m/s) 5%
24oC DB
50% RH
57
Human body thermal balance
The human body metabolism generates heat. The metabolic heat generation rate, usually in
met unit in thermal comfort (1 met = 58.2 W/m2 of body area), depending on the activity
level, age and health conditions.
Temperature regulation of the body is achieved through the control of blood flow rate to the
skin. As the environmental temperature goes up, the blood flow rate to the skin increases to
raise the skin temperature, which in turn, increases the heat transfer rate to the environment.
This process is effective until the skin temperature reaches the core body temperature of 37
oC. At this stage, sweating is initiated, which transfer the metabolic heat to the surroundings
by evaporation.
58
Please rate your overall thermal sensation:
Please rate your overall thermal comfort:
59
Indoor design conditions
The ASHARE thermal comfort chart based
on the effective temperature concept is
shown in Figure.
In the middle of a thermal comfort zone, a person wearing the prescribed clothing would have a
neutral thermal sensation. The middle of the winter and summer comfort zones are 22oC /50%
RH, and 25oC/50%RH respectively. 60
PMV = exp(Met) * L
61
Courtesy of Innova
Courtesy of Innova
62
Indoor air quality (IAQ)
Airborne pollutants are generated within the conditioned space, and can also be transported
from the ambient or neighboring rooms. One major pollutant sources is due to the human body.
Some pollutants can be removed by air cleaning, e.g. filtration for removing fine particles.
In general, ventilation, i.e. supply of outdoor “clean” air can be used to dilute the indoor air
pollutants. This also means that the ventilation air intakes should not be close to high pollution
areas, such as roads with heavy traffic.
Ventilation air needs to be cooled and dehumidified (e.g. in Hong Kong), the amount of
ventilation air should be as small as possible while maintaining a reasonable IAQ.
The minimum ventilation rates are given. The values depend on the number of occupants, and
their activities, as well as the indoor materials.
Table 3. Minimum ASHRAE recommended ventilation rates.
Application Function Design Ventilation air flow rate
occupancy/100m2 per person (L/s)
Office Offices 7 10
Conference room 50 10
Restaurant Lounge 100 15
Dining room 70 10
Kitchen 20 7.5
Retail store Shops, malls 20 1 (L/sm2)
Sport area Ballrooms 100 13
63
Gymnasiums 30 10
Bath effect
1cu ft =0.028 m3
Economy
effect
64
Yaglou, 1937
Before 1900: After 1900: After 1990:
Miasma Influenza pandemic 1918-1919 Sick building syndrome
Lavoisier 1774 O2 Carrier 1919 “manufactured US Indoor Air Act of
Lavoisier 1775 CO2 weather” 1991
Polluted cities Paris and Fanger thermal comfort Eurovent 2002
London, late 18th and Energy crisis 1973 SARS epidemics 2003
19th century WorldVent 2006
EuroVent
2002
WorldVent
20 2007
ASHVE
1914
15 Bilings
1895 Comfort ASHRAE
Yaglou
10 1936 1973 ASHRAE
Flugge
1905 1989/2004
5
ASA ASHRAE
1946 1981
0
1875 1900 1925 1950 1975 2000 2025
Year
65
Internal heat sources in buildings (people, lights, and equipment)
Table 4. Heat and moisture released by people of different activities
Activity level Total heat (W) Sensible Latent Radiant/sensible heat (%)
Adult male Adjusted for heat (W) heat (W) Low air speed High air
M/F/children 0.5 m/s speed 2 m/s
Seated at theater 115 95 65 30 60 27
Moderate office 140 130 75 55 58 38
work
Walking, standing 160 145 75 70 58 38
Light bench work, 235 220 80 140 49 35
factory
Heavy work, 440 425 170 255 54 19
factory
66
Cooling load in an
apartment building
in HK
LI, Z., Chen, W., Deng, S. & Lin, Z. 2006. The characteristics of space cooling load and indoor humidity control for
67
residences in the subtropics. Building and Environment, 41, 1137-1147.
Transient effect
We demonstrate difficulty of thermal storage. Assume the room air temperature, 𝑇𝑎 is constant.
The interior or exterior surface absorbs a fraction 𝛼 of the incident solar radiation or lighting
radiation, 𝐼 𝑡 , and there is a combined convection and thermal radiation (ℎ𝑠 ). The wall is thin with
a high conductivity, and assume its temperature is uniform.
𝑑𝑇
𝜌𝑐𝑙𝐴 𝑠 = 𝐴𝛼𝐼 𝑡 − 𝐴ℎ𝑠 𝑇𝑠 − 𝑇𝑎 ; (1)
𝑑𝑡
where A is the wall surface area, l is the thickness, c is the specific heat capacity, and 𝜌 is the
density of wall material.
The heat flux per unit area entering the room is 𝑞𝑖 = ℎ𝑠 𝑇𝑠 − 𝑇𝑎 (2)
𝑑𝑞𝑖 ℎ𝑠 𝛼ℎ𝑠
Thus + 𝑞𝑖 = 𝛽𝐼 𝑡 , (3) where = and 𝛽 = = 𝛼 .
𝑑𝑡 𝜌𝑐𝑙 𝜌𝑐𝑙
We divide the period 𝑡 into 𝑡 = 𝑁∆𝜏. Assume each interval ∆𝜏 is one hour, and assume 𝐼 𝜏 is constant
during each time interval.
Solution (4) gives 𝑞𝑖 𝑡 = 𝑞𝑖0 𝑒 −𝑁∆𝜏 + 𝑛=0 𝛼𝐼𝑛 𝑒 −𝑛∆𝜏 − 𝑒 − 𝑛+1 ∆𝜏
𝑁−1
(6)
Where 𝐼𝑛 is the short wave radiation at time internal 𝑛 = 0, 𝑁 − 1.
In equation (6), we let 𝑒 −𝑛∆𝜏 − 𝑒 − 𝑛+1 ∆𝜏 = 𝐴𝑛 , i.e. the transfer function coefficients.
This equation is applicable to a 24-hour period, where 𝑁 = 24, ∆𝜏 = 1 hour. We define radiant time
factors
1−𝑒 −∆𝜏 1−𝑒 −∆𝜏
𝑒 −𝑛∆𝜏 ; and we can find out 𝑛=0 𝑟𝑛 = 𝑒 −𝑛∆𝜏 = 1
23 23
𝑟𝑛 =
1−𝑒 − 𝑁∆𝜏 1−𝑒 − 𝑁∆𝜏 𝑛=0
23
Hence 𝑞𝑖 𝑡 = 𝑛=0 𝛼𝐼𝑛 𝑟𝑛 .
This is an amazing formula, showing that the instantaneous heat flux (cooling load) at time 𝑡 can be
calculated using the radiant time factors, which depends only on the geometrical and thermal properties.
69
We write 𝛼𝐼𝑛 as 𝐼𝑠,𝑡−𝑛 , the radiant heat gain rate during the 𝑛-th hour before the current hour
𝑡, and 𝑟𝑛 is the radiant time factor for hour 𝑛. Hence
23
𝑞𝑖,𝑡 = 𝑛=0 𝑟𝑛 𝐼𝑠,𝑡−𝑛 .
Note that this approaches works for a thin wall, which can be assumed to have a uniform
temperature distribution.
In case, a thick wall does not have a uniform temperature distribution, the basic 1D heat
conduction equation can also be analyzed, though the procedure becomes complex. The
results are similar as follows.
Where 𝑐𝑛 is the so-called conduction time factor for hour 𝒕 or interval 𝒏. 𝑄𝑒,𝑡−𝑛 is the
heat imput at the exterior surface of the wall during the 𝑛-th hour before the current hour.
𝑄𝑒,𝑡−𝑛 = 𝑈𝐴 𝑡𝑠𝑎,𝑡−𝑛 − 𝑡𝑎 ; where 𝑈 is the overall heat transfer coefficient for the wall, 𝐴
surface area, 𝑡𝑎 indoor air temperature (constant). 𝑡𝑠𝑎 is the so-called sol-air temperature.
70
The concept of sol-air temperature for an exterior surface is used to simplify the surface heat gain
calculation due to the combined effect of solar radiation, and convection/thermal radiation.
𝑄𝑒 = 𝑞𝑐𝑜𝑛𝑑 = 𝑞𝑠𝑜𝑙 + ℎ𝑐 𝑇𝑎𝑎 − 𝑇𝑠 + ℎ𝑟 𝑇𝑠𝑘𝑦 − 𝑇𝑠 ; or
𝑄𝑒 = 𝑞𝑠𝑜𝑙 + ℎ𝑐 𝑇𝑎𝑎 − 𝑇𝑠 + ℎ𝑟 𝑇𝑎𝑎 − 𝑇𝑠 − ℎ𝑟 𝑇𝑎𝑎 − 𝑇𝑠𝑘𝑦 = 𝑞𝑠𝑜𝑙 + ℎ𝑜 𝑇𝑎𝑎 − 𝑇𝑠 − 𝜀∆𝑅;
where 𝜀∆𝑅 = ℎ𝑟 𝑇𝑎𝑎 − 𝑇𝑠𝑘𝑦 and ℎ𝑜 = ℎ𝑐 + ℎ𝑟 . We let 𝑄𝑒 = ℎ𝑜 𝑇𝑠𝑎 − 𝑇𝑠 by defining a new
variable 𝑇𝑠𝑎 .
Thus, ℎ𝑜 𝑇𝑠𝑎 − 𝑇𝑠 = 𝑞𝑠𝑜𝑙 + ℎ𝑜 𝑇𝑎𝑎 − 𝑇𝑠 − 𝜀∆𝑅, where 𝜀∆𝑅 = ℎ𝑟 𝑇𝑎𝑎 − 𝑇𝑠𝑘𝑦
𝑞𝑠𝑜𝑙 −𝜀∆𝑅
We find the new variable, i.e. sol-air temperature 𝑇𝑠𝑎 = 𝑇𝑎𝑎 + ℎ𝑜
The small term 𝜀∆𝑅 = ℎ𝑟 𝑇𝑎𝑎 − 𝑇𝑠𝑘𝑦 is generally considered to be a small correction term for
the difference between the sky and air temperature.
Thus, the sol-air temperature is an effective driving temperature that incorporates the
contributions of solar radiation, long wave radiation, and convection to the energy interactions at
an external surface.
71
The cooling load calculation by the RTS and CTS methods
The above derived radiant time factor series (RTS) and conduction time factor series (CTS) can
be used to convert heat gains on surfaces to cooling loads by accounting for the time delays
caused by the thermal storage of walls, roofs and surfaces within a zone.
Values for RTS and CTS for different walls and roofs are available in ASHRAE 2013
Fundamentals.
Table 6. Conduction time series (CTS) for two walls and one roof
Curtain wall: mass = 20.9 kg/m2; thermal capacity = 20.4 kJ/m2K; U=0.429 W/m2K
Hour 0 1 2 3 4 5 6 7 8 9 10 11
𝑐𝑛 (%) 25 57 15 3 0 0 0 0 0 0 0 0
Hour 12 13 14 15 16 17 18 19 20 21 22 23
𝑐𝑛 (%) 0 0 0 0 0 0 0 0 0 0 0 0
Brickwall: mass = 304 kg/m2; thermal capacity = 253.5 kJ/m2K; U=0.348 W/m2K
Hour 0 1 2 3 4 5 6 7 8 9 10 11
𝑐𝑛 (%) 2 2 2 3 5 6 7 7 7 7 6 6
Hour 12 13 14 15 16 17 18 19 20 21 22 23
𝑐𝑛 (%) 5 5 5 4 4 3 3 3 3 2 2 1
Metal deck roof: mass = 57.6 kg/m2; thermal capacity = 57.2 kJ/m2K; U=0.297 W/m2K
Hour 0 1 2 3 4 5 6 7 8 9 10 11
𝑐𝑛 (%) 0 10 22 20 14 10 7 5 4 3 2 1
Hour 12 13 14 15 16 17 18 19 20 21 22 23
72
𝑐 (%) 1 1 0 0 0 0 0 0 0 0 0 0
Building energy estimation and modelling: Degree-day method with the following assumptions
The heating load depends only on the indoor and outdoor dry-bulb temperatures
The efficiency of the heating system is constant
The internal heat gain due to solar radiation, people, lights and equipment are constant
Air infiltration affects only the sensible heating load
Adding (1) and (2), we obtain the total energy loss due to envelop heat transfer and infiltration,
𝑄𝑙𝑜𝑠𝑠 = 𝑈𝑒 𝐴𝑒 + 𝑚𝑖 𝑐𝑝 𝑇𝑖𝑛 − 𝑇𝑎𝑚𝑏 = 𝑈𝑜 𝐴𝑜 𝑇𝑖𝑛 − 𝑇𝑎𝑚𝑏 ; where 𝑈𝑜 𝐴𝑜 is the effective heat transfer conductance.
The rate of heat gain due to solar energy, people, lights and equipment 𝑄𝑖𝑛𝑡 = 𝑄𝑠𝑜𝑙 + 𝑄𝑝𝑒𝑜 + 𝑄𝑙𝑖𝑔ℎ𝑡𝑠 + 𝑄𝑒𝑞𝑢𝑖𝑝
At steady state the heating load, 𝐿ℎ = 𝑄𝑙𝑜𝑠𝑠 − 𝑄𝑖𝑛𝑡 = 𝑈𝑜 𝐴𝑜 𝑇𝑖𝑛 − 𝑇𝑎𝑚𝑏 − 𝑄𝑖𝑛𝑡
𝑄𝑖𝑛𝑡
Define a balance temperature. 𝐿ℎ = 𝑈𝑜 𝐴𝑜 𝑇𝑏𝑎𝑙 − 𝑇𝑎𝑚𝑏 , and 𝑇𝑏𝑎𝑙 = 𝑇𝑖𝑛 − . The heating load is zero when 𝑇𝑏𝑎𝑙 = 𝑇𝑎𝑚𝑏 .
𝑈𝑜 𝐴𝑜
The heating degree-day method is based on the observation that the heating system needs to supply heat only when 𝑇𝑎𝑚𝑏 <
𝑇𝑏𝑎𝑙 . Thus the total amount of heat to be supplied by the heating system during the heating season is
𝑡𝑠 +
𝐸ℎ𝑒 = 𝑈𝑜 𝐴𝑜 0
𝑇𝑏𝑎𝑙 − 𝑇𝑎𝑚𝑏 𝑑𝑡; (note: the superscript “+” indicates only positive values included.)
For a period, such a month,𝐸ℎ𝑒,𝑚 = 𝑈𝑜 𝐴𝑜 𝐷𝐷 𝑚, where 𝐷𝐷 𝑚 is the total heating degree-days for the month.
𝑈𝑜 𝐴𝑜 𝐷𝐷 𝑚
If the furnace efficiency is 𝑓 , the amount of fuel energy required 𝐸𝑓𝑢,𝑚 = .
𝑓
A similar cooling degree days can be defined. Thus we have heating and cooling degree days (HDD and CDD).
73
Every degree counts
74
Example 5. A thin vertical metal wall of a building at a location with a northern latitude of 40
faces 35 east of south. The measured direct beam and diffuse solar radiation intensities at the
location on August 15 at 10 am solar time are 580 W/m2, and 148 W/m2 respectively. The
reflectivity of the ground surrounding the wall is 0.3. With such information, we can estimate the
direct radiation incident on unit area of the wall surface is 331.8 W/m2; the diffuse radiation
incident on unit area is 133.57 W/m2; and the ground-reflected radiation falling on unit area of
surface is 89.4 W/m2. Thus the total solar radiation incident on unit area of the wall is 554.8
W/m2. (Note the solar radiation and solar shading calculations are beyond the scope of the
lecture).
The average emissivity of the wall surface is 0.85. The ambient temperature and the inside air
temperature are 28oC and 23oC respectively. The overall external and internal heat transfer
coefficients are 35 W/m2K and 8.5 W/m2K respectively. Assume that the heat capacity and the
thermal resistance of the wall are negligible,
Calculate
(a) Sol-air temperature; and
(b) The temperature of the wall.
75
Example 5. A thin vertical metal wall of a building at a location with a northern latitude of 40 faces 35 east of south. The
measured direct beam and diffuse solar radiation intensities at the location on August 15 at 10 am solar time are 580 W/m 2,
and 148 W/m2 respectively. The reflectivity of the ground surrounding the wall is 0.3. With such information, we can estimate
the direct radiation incident on unit area of the wall surface is 331.8 W/m2; the diffuse radiation incident on unit area is 133.57
W/m2; and the ground-reflected radiation falling on unit area of surface is 89.4 W/m2. Thus the total solar radiation incident on
unit area of the wall is 554.8 W/m2. (Note the solar radiation and solar shading calculations are beyond the scope of the
lecture).
The average emissivity of the wall surface is 0.85. The ambient temperature and the inside air temperature are 28 oC and 23oC
respectively. The overall external and internal heat transfer coefficients are 35 W/m2K and 8.5 W/m2K respectively. Assume
that the heat capacity and the thermal resistance of the wall are negligible,
Calculate
(a) Sol-air temperature; and
(b) The temperature of the wall.
Since the thermal capacity and the thermal resistance of the wall are negligible, the net heat
flow rate
76
Example 8. The wall of a building is 8 m long and 3 m high. It has two equal sized
windows of length 2.5 m and height 2 m located in it. The overall heat transfer
coefficient for the wall is 0.246 W/m2K. The overall heat transfer coefficient for the
window is 3.02 W/m2K. The outdoor and indoor design temperatures are -10oC and
20oC, respectively. Calculate the total heat load due to heat loss through the wall.
77
Example 9. For a two-story building, the infiltration rate is 0.068
m3/s. The indoor design conditions are 22oC and 40% relative
humidity. The outdoor ambient air is saturated at -10 oC. The
pressure is 101 kPa. Calculate the sensible and latent heat loads
due to infiltration of ambient air.
78
Example 9_new. For a three-story small house (70m2 each story, 8.5 m tall), the
infiltration rate is 1 ACH. The indoor design conditions are 24oC and 50% relative
humidity. The design outdoor air condition is 33.2oC DB and 26.3oC WB. Calculate
the sensible and latent heat loads due to infiltration of ambient air.
79
A lot of temperatures
1. Sensible air temperature
2. Dew point temperature
3. Dry bulb temperature
4. Wet bulb temperature
5. Sol-air temperature
6. Saturation temperature
7. Mean radiant temperature (not discussed)
8. Operative temperature (not discussed)
80
Chapter 4
Couvreur, V., Ledder, G., Manzoni, S., Way, D. A.,
Muller, E. B., & Russo, S. E. (2018). Water transport
through tall trees: A vertically explicit, analytical
model of xylem hydraulic conductance in
stems. Plant, cell & environment, 41(8), 1821-1839.
The total pressure entering the return air duct decreases progressively as it passes through the
grill, return duct and filter. There is a sharp increase across the fan due to the work input by
the fan.
The heating and cooling coils also result in a significant loss. The air flow divides into three
streams. Each air stream loses total pressure due to friction in the duct, fittings and diffuser.
Figure 1. (a) A typical duct network of an air conditioning system; (b) the total pressure (total head)
distribution.
Total pressure distribution
At any section of a duct system, the total mechanical energy per unit mass of fluid, 𝐸𝑚 (J/kg) is
𝑝 𝑣2
𝐸𝑚 = + + 𝑔𝑧; where the three RHS terms are potential energy due to static pressure (N/m2),
𝜌 2
kinetic energy due to fluid velocity and gravitational potential energy due to vertical height.
As fluid flows through a duct work, a fraction of its mechanical energy is converted to internal
thermal energy due to fluid friction, resulting in a small rise of the fluid temperature. We usually
call the change in the mechanical energy between any two sections as an energy loss, 𝐸𝐿 (J/kg), or
expressed by the head loss ∆𝐻 and total pressure loss ∆𝑃 (Pa)
ASHRAE duct fitting database (2012) contains an extensive data set of loss coefficients for nearly all air duct fittings.
ASHRAE Handbooks also contain a limited list. Here we cited some representative values for the typical fittings.
For fittings like contractions with unequal inlet and outlet areas, the loss coefficient may be changed from section o (outlet)
to i (inlet) by applying the mass balance equation.
1 𝐴𝑖 2
∆𝑝𝑡,𝑖 = 𝑐𝑖 𝑝𝑣𝑖 = 𝑐𝑖 𝜌𝑉𝑖2 , where 𝑐𝑖 = 𝑐𝑜 .
2 𝐴𝑜
For diverging and converging flow junctions, the total pressure loss (mechanical energy) through the straight (main) section
is
1
∆𝑝𝑡,𝑠 = 𝑐𝑠 𝑝𝑣𝑠 = 𝑐𝑠 𝜌𝑉𝑠2 ;
2
where 𝑐𝑠 – local loss coefficient referenced to section s [-]; 𝑝𝑣𝑠 - velocity pressure at the selected referenced section s [Pa];
𝑉𝑠 -velocity at the selected referenced section s [m/s].
In converging sections, mechanical energy is exchanged between two fluid streams moving at different velocities due to
turbulent mixing. Consequently, under some flow conditions, the mechanical energy per unit mass of the slower moving
stream can increase due to mixing. This results in a negative loss coefficient for the lower velocity stream.
Total pressure loss in duct sections
The total pressure loss in a duct section consisting of straight sections and fittings can
be obtained
𝐿𝑛 1 2 1 2
∆𝑃𝑡𝑜𝑡 = 𝑛 𝑓𝑛 𝐷 2 𝜌𝑉𝑛 + 𝑛 𝑐𝑛 2 𝜌𝑉𝑛 ;
𝑛
The loss coefficient 𝑐𝑛 and the friction coefficient 𝑓𝑛 have to be referenced to the
appropriate sections. It is possible to have a computer code doing all the calculations.
Total pressure – the rise in total pressure cross the fan. ∆𝑝𝑡 = 𝑝𝑡,𝑜 − 𝑝𝑡,𝑖 , where i and o denote
conditions at the inlet and outlet respectively.
Static pressure - the rise in the static pressure across the fan. ∆𝑝𝑠 = 𝑝𝑠,𝑜 − 𝑝𝑠,𝑖 .
Ideal power input – required to operate the fan under the given conditions. Assuming the
𝑚 𝑝𝑡,𝑜 −𝑝𝑡,𝑖
temperature rise of air to be negligible, and the density to be constant, 𝑊𝑖𝑑 = ,
𝜌
where 𝑚 is the mass flow rate of air.
Fan efficiency – the ratio of the ideal power required to the actual shaft power input of the
𝑊
fan, 𝑓 = 𝑊 𝑖𝑑 . Usually fans are driven by electric motors through belt drives. If the efficiency
𝑠ℎ
𝑊
of the electric motor is 𝑚 , then the required electric power input 𝑊𝑒𝑙 = 𝑖𝑑 .
𝑚 𝑓
(a) Centrifugal
Axial
(a)
In an axial fan, the air is moved axially by the propeller blades, mounted on a concentric
hub, rotated by an electric motor. The kinetic energy imparted to the air is converted to a
pressure rise by the stationary blades located downstream of the propeller. The stationary
blades also help improve the efficiency of the fan by reducing the swirl of the air. As axial
flow fans are located axially in the flow duct, they do not require a change in flow direction
of the air. Axial flow fans are used mainly in high volume flow applications.
(b)
(a)
(b)
A centrifugal fan consists of an impeller with a series of blades attached to a hub, driven by
an electric motor. The power is transferred from motor to the fan through a belt drive with
pulleys attached to the motor and the impeller. Air enters at the center of the impeller in a
direction normal to the plane of the figure as shown above. The centrifugal force created by
the rotating blades accelerates the air in the radial direction. The kinetic energy gained by
the air is converted to a high static pressure in the diffusing section of the fan.
Fan-duct network interaction
For the simple duct network in Figure 5, we have duct sections and fittings before and after the
fan. For duct sections 1-2 and 3-4:
𝐿 1 1
∆𝑃1−2 = 𝑝0 − 𝑝2 = 𝑛 𝑓𝑛 𝐷𝑛 2 𝜌𝑉𝑛2 + 𝑛 𝑐𝑛 2 𝜌𝑉𝑛2 ;
𝑛 1−2
𝐿𝑛 1 2 1 2
∆𝑃3−4 = 𝑝3 − 𝑝0 = 𝑛 𝑓𝑛 𝐷 2 𝜌𝑉𝑛 + 𝑛 𝑐𝑛 2 𝜌𝑉𝑛 . (5)
𝑛 3−4
where 𝑝0 is the uniform pressure of the space to which the inlet and outlet of the duct system are
connected at locations 1 and 4, respectively.
𝑄
The velocity 𝑉𝑛 = 𝐴 at a section with an area of 𝐴𝑛 . Substituting for the velocities in term of flow
𝑛
rate in (5), we have
∆𝑃1−2 = 𝑝0 − 𝑝2 = 𝑐12 𝑄2 and ∆𝑃3−4 = 𝑝3 − 𝑝0 = 𝑐34 𝑄2 ; where 𝑐12 and 𝑐34 are terms involving
various friction factors, duct dimensions and the loss coefficients of the fittings, and the two terms
can be assumed constant for a particular system.
We have ∆𝑝𝑑𝑠 = 𝑝3 − 𝑝2 = 𝑐12 + 𝑐34 𝑄2 , and the variation of the total pressure loss in the duct
system ∆𝑝𝑑𝑠 with the volume flow rate 𝑄 is the system curve. Ideally, the system curve is parabola.
The point of intersection of the fan characteristic and the system curve
gives the fan pressure and flow rate under steady operating conditions,
i.e. the operating point of the system.
Many AHUs operate under part-load conditions much of the year except during peak
load conditions. Using variable speed drive can reduce energy use. In a VAV system,
both damper control and variable speed drive are used, hence VAV systems can reduce
significantly fan power under part-load conditions.
Design methods for duct systems
The main requirements of a duct system are:
It supplies the specified rates of air flow to the different zones;
It is economical in initial cost, operating cost and also the cost of the building space
The noise level is acceptable.
100
Equal friction method
Using the system shown in Figure 7 as an example, here is how the method works.
(i) The flow rates to the zones 𝑄3 ,𝑄4 and 𝑄5 are estimated from the thermal loads of the zones.
If the air density is assumed constant, then the flow rate from the fan is 𝑄 = 𝑄3 + 𝑄4 + 𝑄5 .
(ii) The layout of the duct network is based on the location of the zones, and AHU. For each
duct section, the length, and the number and type of fittings are specified.
(iii)The pressure for each zone is also specified. Usually, all the zones are at the same pressure.
(iv)An initial unit frictional pressure loss 𝑝0 is chosen. This value is applied to all duct sections
of the duct network.
(v) Obtain the diameter and the air velocity from the duct friction chart based on the flow rate 𝑄
and pressure loss.
(vi)As the value of 𝑝0 is usually chosen based on past design experience, the duct sizes obtained
in (v) may not always deliver the specified flow rates to the three zones. This has to be
checked by calculating the total pressure loss from 1 to exits at 6, 8 and 10, ∆𝑝𝑡,1−6,
∆𝑝𝑡,1−8,and ∆𝑝𝑡,1−10 . The three total pressure losses must be equal to deliver the specified
flow rates.
(vii)If the pressure losses are not equal, then dampers may have to be installed at the exit of
some ducts to increase the pressure drop artificially to achieve the desired flow rates to the
zones. An alternative approach is to change the diameter of some branches of the network to
obtain the desired pressure losses.
Figure 2. Friction chart for
round ducts with curves for
constant diameter and
constant velocity.
(Temperature = 20oC, density
1.2041 kg/m3, and roughness
= 0.09 mm.
Calculate
(a) The static pressure and velocity pressure at section 1
(b) The total frictional pressure loss from 1 to 10;
(c) The pressure loss in the fittings from 1 to 10.
(d) The total static and velocity pressure at section 10
(e) The static and velocity pressure heads at section 1
(f) The total mechanical energy loss per unit mass of air from 1 to 10; and
(g) The rate of mechanical energy loss from 1 to 10.
Example 1. The supply duct of an air conditioning system, shown in Figure 8 delivers 0.8 m3/s of air at 20oC dry-bulb temperature. The diameter of the duct is 0.4 m, and the
total length of the straight duct sections from 1 to 10 is 80 m. The fan generates a total pressure of 125 Pa at section 1.
The loss coefficient for each right angle bend is 0.1. The friction coefficient is 0.0186. The density of air is 1.2 kg/m 3.
Calculate
(a) The static pressure and velocity pressure at section 1
(b) The total frictional pressure loss from 1 to 10;
(c) The pressure loss in the fittings from 1 to 10.
(d) The total static and velocity pressure at section 10
(e) The static and velocity pressure heads at section 1
Solution
(f) The total mechanical energy loss per unit mass of air from 1 to 10; and
(g) The rate of mechanical energy loss from 1 to 10. 𝑄 𝑄 4×0.8
The air velocity through the duct is 𝑉 = = 1 = = 6.37 m/s
𝐴 𝜋𝐷2 3.14×0.42
4
1 1
(a) Velocity pressure at section 1 is 𝑝𝑣1 = 𝜌𝑉 2 = × 1.2 × 6.372 = 24.34 Pa
2 2
Hence the static pressure 𝑝𝑠1 = 𝑝𝑡1 − 𝑝𝑣1 = 125 − 24.34 = 100.9 Pa
𝐿 1 1
(b) and (c) The total pressure loss ∆𝑝𝑡𝑜𝑡 = 𝑓 𝜌𝑉 2 + 𝑛 𝑐𝑛 2 𝜌𝑉𝑛
2
𝐷 2
80 1 2 1
∆𝑝𝑡𝑜𝑡,1−10 = 0.0186 × × 1.2 × 6.37 + 4 × 0.1 × 1.2 × 6.372 = 90.56 + 9.73 =
𝐷0.4 2 2
100.3 Pa;
(d) The total pressure at section 10, 𝑝𝑡,10 = 𝑝𝑡1 − ∆𝑝𝑡𝑜𝑡,1−10 = 125 − 100.3 = 24.7Pa
𝑝𝑠10 = 𝑝𝑡10 − 𝑝𝑣1 = 24.7 − 24.34 = 0.36 Pa
𝑝𝑠1 100.9 𝑝𝑣1
(e) Static pressure head at 1 𝐻𝑠1 = = = 8.57 m; and the velocity heat at 1 is 𝐻𝑣1 = =
𝜌𝑔 1.2×9.81 𝜌𝑔
24.34
= 2.06 m;
1.2×9.81
(f) The mechanical energy loss per kg of air from 1 to 10 is
∆𝑝𝑡𝑜𝑡,1−10 100.3
∆𝐸𝑚,1−10 = = = 83.58 J/kg;
𝜌 1.2
(g) The rate of mechanical energy loss from 1 to 10
∆𝑊𝑙𝑜𝑠𝑠 = 𝑄∆𝑝𝑡𝑜𝑡,1−10 = 0.8 × 100.3 = 80.24 W. Thus the mechanical energy loss or conversion to
heat results in a rise in temperature of air of about 0.08 oC.
Example 2. The dimension of a duct system, supplying air at 20oC to two zones 5 and 6.
The flow rates to the two zones are 1.4 m3/s and 1.0 m3/s respectively. The fan generates a
total pressure of 285 kPa at the entrance section 1.
Calculate
(a) (i) the velocity pressures in the duct sections; (ii) the total pressures at all sections from
1 to 6.
(b) If a damper is installed just upstream of section 5 to make the total pressures at 5 and 6
equal, calculate the required pressure loss through the damper.
Example 2. The dimension of a duct system, supplying air at 20oC to two zones 5 and 6. The flow rates to the two zones are
1.4 m3/s and 1.0 m3/s respectively. The fan generates a total pressure of 285 kPa at the entrance section 1. Calculate (a) (i)
the velocity pressures in the duct sections; (ii) the total pressures at all sections from 1 to 6. (b) If a damper is installed just
upstream of section 5 to make the total pressures at 5 and 6 equal, calculate the required pressure loss through the damper.
Solution: The air velocities in the duct sections 1-2, 3-5 and 4-6 are given by
𝑄12 4×2.4 𝑄35 4×1.4 𝑄46 4×1.0
𝑉12 = = = 8.49 m/s; 𝑉35 = = = 7.13 m/s; 𝑉46 = = = 7.95 m/s
𝐴12 3.14×0.62 𝐴35 3.14×0.52 𝐴46 3.14×0.42
The velocity pressure in the duct sections are
21 2 1
𝑝𝑣12 = 𝜌𝑉12 = 0.5 × 1.2 × 8.492 = 43.2Pa; 𝑝𝑣35 = 𝜌𝑉35 = 0.5 × 1.2 × 7.132 = 30.5Pa;
2 2
1 2
𝑝𝑣46 = = 0.5 × 1.2 × 7.952 = 37.9Pa.
𝜌𝑉46
2
The frictional pressure losses in the straight duct sections (find the friction f values)
∆𝑃𝑓12 = 53.5 Pa; ∆𝑃𝑓35 = 42.7 Pa; ∆𝑃𝑓46 = 120.8 Pa
The loss for the 90 bend 4-6, diverging tee-junction 2-3-4
2 1 1
∆𝑝𝑏𝑒𝑛𝑑−46 = 𝑐𝑜 𝜌𝑉46 = × 0.11 × 1.2 × 7.952 = 4.17 Pa;
2 2
1 2 1
∆𝑝𝑠𝑡𝑟𝑎𝑖𝑔ℎ𝑡−23 = 𝑐𝑠𝑡𝑟𝑎𝑖𝑔ℎ𝑡 𝜌𝑉12 = × 0.14 × 1.2 × 8.492 = 6.05 Pa;
2 2
1 2 1
∆𝑝𝑏𝑒𝑛𝑑−24 = 𝑐𝑏𝑒𝑛𝑑 𝜌𝑉12 = × 1.2 × 1.2 × 8.492 = 51.89 Pa;
2 2
The total pressure at different sections are
𝑝𝑡1 = 285 Pa; 𝑝𝑡2 = 𝑝𝑡1 − ∆𝑝𝑓12 = 285 − 53.5 = 231.5Pa;
𝑝𝑡3 = 𝑝𝑡2 − ∆𝑝𝑠𝑡𝑟𝑎𝑖𝑔ℎ𝑡−23 = 231.4 − 6.05 = 225.35Pa; 𝑝𝑡4 = 𝑝𝑡2 − ∆𝑝𝑏𝑒𝑛𝑑−24 = 231.4 − 51.89 =
179.5a;
𝑝𝑡5 = 𝑝𝑡3 − ∆𝑃𝑓35 = 225.35 − 42.7 = 182.65Pa; 𝑝𝑡6 = 𝑝𝑡4 − ∆𝑃𝑓46 = 179.5 − 120.8 = 58.7Pa;
Solution:
𝑄 1.6
The outlet velocity of the air is given by 𝑉𝑜 = 𝐴 = 0.62×0.32 = 8.06 m/s.
𝑜
1 1
The velocity pressure at the outlet of the fan is 𝑝𝑣𝑜 = 2 𝜌𝑉𝑜2 = 2 × 1.2 × 8.062 = 39.0 Pa
Solution: The ideal power input is 𝑊1 = 𝑄1 𝑝𝑡1 = 4.5 × 11.5 = 51.75 kW.
If the fan speed increases from 180 rpm to 200 rpm (density is constant)
𝜔 200
𝑄2 = 2 𝑄1 = × 4.5 = 5.0 m3/s
𝜔1 180
𝜔 2 200 2
𝑝𝑡2 = 2 𝑝𝑡1 = × 11.5 = 14.2 kPa
𝜔1 180
𝜔2 3 200 3
𝑊2 = 𝑊1 = × 51.75 = 71 kW
𝜔1 180
If the speed remains constant but the ambient air temperature increases to 40 oC, we assume the air is an
ideal gas, 𝑝 = 𝜌𝑅𝑇, hence
𝜌2 𝑇1 𝑇1 273+20
= ; , 𝜌2 = 𝜌1 = 1.2 × = 1.123 kg/m3.
𝜌1 𝑇2 𝑇2 273+40
The flow rate is independent of the density, so remains constant. For the same fan seed
𝜌2 1.123
𝑝𝑡2 = 𝑝𝑡1 = × 11.5 = 10.76 kPa;
𝜌1 1.2
𝜌2 1.123
𝑊= 𝑊1 = × 51.75 = 48.4kW.
𝜌1 1.2
Example 5. Figure 5.1 shows a supply air duct system and all ducts are circular. The flow rates
and duct lengths are indicated in the figure. The loss coefficient for the duct exists is 0.6.
(a) Use the equal friction method to size all ducts;
(b) Determine where dampers should be located to achieve the desired air flow rates to the
spaces;
(c) Can the pressure losses in the duct runs be balanced by changing the duct diameters?
We choose air velocity in section 1-2 is 9 m/s
Table 4.1 The chosen duct diameters and other parameters for each section
(b) Based on the flow rates and duct diameters, we obtain the loss coefficients of the two diverging tee-
junctions, and the 90 bend from tabulated data set.
𝑐2,3 = 0.13; 𝑐2,11 = 2.4; 𝑐4,5 = 0.14; 𝑐4,9 = 1.5; 𝑐6,7 = 0.11.
The total pressure losses in different duct runs are
𝑃1 − 𝑃8 = 𝑃12 + 𝑃23 + 𝑃34 + 𝑃45 + 𝑃56 + 𝑃67 + 𝑃78 + 𝑃𝑒𝑥𝑖𝑡 = 13.95 + 0.13 × 48.3 + 11.2 + 0.14 ×
42.7 + 18.6 + 0.11 × 28.42 + 7.4 + 0.6 × 28.42 = 83.6Pa;
𝑃1 − 𝑃12 = 𝑃12 + 𝑃2,11 + 𝑃11,12 + 𝑃𝑒𝑥𝑖𝑡 = 13.95 + 2.4 × 48.3 + 6.5 + 0.6 × 23.3 = 150 Pa;
𝑃1 − 𝑃10 = 𝑃12 + 𝑃23 + 𝑃34 + 𝑃4,9 + 𝑃9,10 + 𝑃𝑒𝑥𝑖𝑡 = 13.95 + 0.13 × 48.3 + 11.2 + 1.5 × 42.7 +
5.58 + 0.6 × 32.64 = 120.6 Pa.
Solution continued
The fan has to deliver the largest pressure loss of 150 Pa for section1-12. To balance the pressure losses in the
three duct runs, dampers have to be installed along section 7-8 and 9-10 to artificially introduce additional
pressure losses, i.e. with magnitudes of ∆𝑃7−8 = 150 − 83.6 = 66.4Pa, and ∆𝑃9−10 = 150 − 120.6 =
29.4Pa.
(c) We can reduce the duct 7-8 and duct 9-10 diameters to increase their pressure losses. After some
calculations, the solutions are for section 7-8, the diameter is reduced to 0.37 m so that velocity = 13.4 m/s
and pressure in 1-8 is 120 Pa. For section 9-10, the diameter is reduced to 0.42 m so that velocity = 14 m/s
and pressure in 1-10 is 145 Pa. Note that these changes bring the unit pressure loses in 7-8 and 9-10 to the
boundary of the recommended region in the friction chart.
𝑄12 = 1.45; 𝑄34 = 1.45; 𝑄56 = 1.9; 𝑄78 = 2.35; 𝑄9,10 = 0.45; 𝑄11,12 = 0.45
To apply the equal friction method, we shall use a unit pressure drop of 0.90 Pa/m, which is closer to the lower bound of the
recommended values.
Section 𝐿, m 𝑄, m3/s 𝐷, m 𝑉, m/s 𝑃𝑣𝑒𝑙 , Pa 𝑃𝑙𝑜𝑠𝑠 , Pa
1-2 9.2 1.45 0.524 6.7 27.2 8.28
3-4 4.6 1.45 0.524 6.7 27.2 4.14
5-6 4.6 1.9 0.581 7.2 31.0 4.14
7-8 9.5 2.35 0.629 7.55 34.4 8.55
9-10 9.2 0.45 0.337 5.04 15.3 8.28
11-12 9.2 0.45 0.337 5.04 15.5 8.28
Table 6.1 The chosen duct diameters and other parameters for each section
(b) Based on the flow rates and duct diameters, we obtain the loss coefficients of the two diverging tee-juctions, and the 90
bend from tabulated data set.
𝑐2,3 = 0.11; 𝑐4,5 = 0.36; 𝑐10,5 = −0.28; 𝑐6,7 = 0.25; 𝑐12,7 = 0.21.
The total pressure losses in different duct runs are
𝑃1 − 𝑃8 = 𝑃12 + 𝑃23 + 𝑃34 + 𝑃45 + 𝑃56 + 𝑃67 + 𝑃78 = 8.28 + 0.11 × 27.2 + 4.14 + 0.36 × 31 + 4.14 + 0.25 ×
34.4 + 8.55 = 48Pa;
𝑃9 − 𝑃8 = 𝑃9,10 + 𝑃10,5 + 𝑃56 + 𝑃67 + 𝑃78 = 8.28 − 0.28 × 31 + 4.14 + 0.25 × 34.4 + 8.55 = 21 Pa;
𝑃11 − 𝑃8 = 𝑃11,12 + 𝑃12,7 + 𝑃78 = 8.28 + 0.21 × 34.4 + 8.55 = 24 Pa.
The fan has to deliver the largest pressure loss of 48 Pa for section1-8. To balance the pressure losses in the three duct runs,
dampers have to be installed along section 9-10 and 11-12 to artificially introduce additional pressure losses, i.e. with
magnitudes of ∆𝑃9−10 = 48 − 21 = 27Pa, and ∆𝑃11−12 = 48 − 24 = 24Pa.
Example 7. A circular supply air duct system is shown schematically in Figure 7.1. The flow rate through
section 1-2 is 40 m3/s, and the air velocity is 18 m/s. The section 3-4 with a 90
Bend is 12 m long. The flow rate through section 3-4 is 24 m3/s. Determine the diameter of section 3-4 using
the static regain method.
1 1
4𝑄1 2 4×40 2
Solution: Section 1-2 diameter 𝐷1 = = = 1.68 m
𝜋𝑉2 3.14×18
Applying the modified Bernoulli’s equation between 2 and 4
1 1 𝐿34 1 1
𝑝𝑠2 + 𝜌𝑉22 = 𝑘𝑠 𝜌𝑉22 + 𝑓34 𝜌𝑉32 + 𝑘𝑏 𝜌𝑉32 + 𝑝𝑠4 + 𝜌𝑉32 ; (1)
2 2 𝐷3 2 2
Where 𝑘𝑠 and 𝑘𝑏 are the loss coefficients of the straight section of the tee-junction 2-3-5 and the 90 bend of
section 3-4.
1 1
For static regain method, 𝑝𝑠2 = 𝑝𝑠4 . As 𝑄1 = 𝜋𝐷12 𝑉2 ; 𝑄2 = 𝜋𝐷22 𝑉3
4 4
Equation (1) becomes
2
1−𝑘𝑠 𝑄1 1+𝑘𝑏 𝐿34 1−𝑘𝑠 40 2 1+𝑘𝑏 12𝑓34
= + 𝑓34 ; with numerical values = + (2)
𝐷14 𝑄2 𝐷24 𝐷25 1.684 24 𝐷24 𝐷25
Equation (2) can be numerically solved for the diameter 𝐷2 . However, the loss coefficient 𝑘𝑠 for the tee-
𝑄2 𝐴2 𝐷2 2
junction is a function of and = . The loss coefficient for the 90 bend is a function of the diameter
𝑄1 𝐴1 𝐷1
𝐷2 .
The main steps in a numerical solution are as follows: Guess an initial value of 𝐷2 . Obtain 𝑘𝑠 and 𝑘𝑏 . Obtain
the friction factor 𝑓34 . Iterate using simple methods such as bifurcation, then we can obtain 𝐷2 = 142.6 cm.
Design of room air distribution systems
Non-dimensional parameters –
Reynolds number, Re, is defined as the ratio between the time scales due
t UL
to overall convection and due to viscous diffusion. Re d
tc
Rayleigh number is defined as the ratio of the time scales due to thermal
g TL3
diffusion and due to thermal convection Ra
Mixing ventilation
Displacement ventilation In the mixed convection literature, Ar is written as Gr/Re2, where Grashof
number is defined as:
2 g TL3
Gr ; Ra = Gr Pr, where Pr is the Prantl number, and Pr =
2
.
Estimate the Reynolds number and the Archimedes number for the two ventilation systems
respectively. Assume the following air properties: = 1.84×10-5 Pa∙s, = 1.189 kg m-3.
Parameters Ranges
Occupied zone air temperatures 16oC-26 oC
Occupied zone air temperature gradient <3 K m-1
Surface temperature of radiators <80 oC
Supply air temperature 9 oC-50 oC
Supply air velocity
- mixing ventilation 2.0-8.0 m s-1
- displacement ventilation 0.1-0.4 m s-1
Archimedes number
- Buoyancy-controlled mixing ventilation <0.01
- Buoyancy-controlled displacement ventilation >0.5
Water distribution systems
Also called hydronic systems, have several practical advantages over air distribution:
(a) the size of hear source is smaller;
(b) water pipes requires less space than air ducts; and
(c) for heating applications, higher temperature are more practical with water than with air.
(d) Air distribution is commonly used in residential and small commercial applications where
the distances from the heat source to the conditioned spaces are usually short. For large
central air conditioning systems, it is more common to use hot water and chilled water as
distribution media.
Where 𝑊𝑖𝑛 [J/kg] is the work input by the pump per unit mass of fluid. The mechanical energy
loss per unit mass, 𝑊𝑙𝑜𝑠𝑠 includes frictional losses I the pipe and dynamic losses in the fittings.
Note for the open loop system, one part is in direct contact with air, as in the cooling tower
water loop of the air conditioning system as shown in Figure 15.
note that change in elevation for the close loop system is generally negligible. The chilled water
loop is an example of a close loop system.
Head losses
∆𝑃𝑓
Frictional loss in straight pipes ∆ℎ𝑓 = =
𝜌𝑔
𝐿 𝑉2
𝑓 𝐷 2𝑔
𝑉2
Dynamic loss in fittings ∆ℎ12 = 𝑘 2𝑔
𝐿 𝑉2
Total head loss ∆ℎ𝑙𝑜𝑠𝑠 = 𝑓 +
𝐷 2𝑔
𝑉2 𝐿 𝑉2
𝑘 = 𝑓 +𝑘
2𝑔 𝐷 2𝑔
𝑉2 𝐿𝑒𝑞 𝑉2
Equivalent length 𝑘 2𝑔
=𝑓 𝐷 2𝑔
, we
𝑘𝐷 𝑘
have 𝐿𝑒𝑞 = 𝑓
=
𝑓
𝐷.
Typical profiles of total pressure head, work input and efficiency with flow rate are shown in Figure 18. As the flow rate
increases, the head decreases progressively until at the maximum flow rate the head becomes zero. The power input at
zero flow rate is used to stir the fluid in the pump and is all converted into thermal internal energy. The power input
increases near linearly with flow rate.
The ideal work input to the fluid per unit mass of fluid by the pump
∆𝑃𝑡
𝑊𝑖𝑑𝑒𝑎𝑙 = = 𝑔∆ℎ𝑡 , where 𝜌 is the density of fluid, and ∆ℎ𝑡 is the total head.
𝜌
𝑚∆𝑃𝑡
The ideal power input to the fluid by the pump 𝑊𝑖𝑑𝑒𝑎𝑙 = = 𝑄𝜌𝑔∆ℎ𝑡 ;
𝜌
The actual work input per unit mass, 𝑊𝑎𝑐𝑡 is larger than the ideal work.
𝑊𝑖𝑑𝑒𝑎𝑙 𝑔∆ℎ𝑡 𝑄𝜌𝑔∆ℎ𝑡
The pump efficiency 𝑡 = 𝑊𝑎𝑐𝑡
= 𝑊𝑎𝑐𝑡
= 𝑊𝑎𝑐𝑡
Cavitation
𝑣𝑠2
The available NPSH: NPSH𝑎𝑣𝑎 = ℎ𝑠 + 2𝑔 − ℎ𝑣 ;
where ℎ𝑠 and 𝑣𝑠 are the head and the velocity at the
impeller entrance and ℎ𝑣 is the vapor pressure of the
liquid at the prevailing temperature.
Design of water distribution systems
A typical system consists of straight pipes, fittings like bends and tee-joints, valves, heat
exchangers, an expansion tank, and pumps. For close systems, two arrangements are
possible, i.e. the two-pipe direct-return systems, and two-pipe reverse-return systems.
Direct return pipe system - the terminal units in the zones A, B and C are supplied with
hot water or chilled water. Its drawback is that the available pressure to the three terminal
units are different due to the different lengths of the supply and return pipes to and from
the units. More water would be supplied to A, e.g. if the A pipes are much shorter. If the available pressures
differ significantly, then it would be difficult to ensure the desired flow rates to the three units, even with the help of
control valves.
Reserve-return system - the total lengths of the piping from the boiler or the chiller to the
three terminal units, including the supply and return pipe lengths, are approximately
equal. Therefore, the available pressures and the flow rates to the three units will not
differ by much. Its disadvantage is the additional pipe required.
(a) direct-return system (b) reverse-return system.
Each load is one or more cooling coils in AHU
http://www.automatedbuildings.com/news/apr02/art/htmn/htmn.htm
Design of water distribution systems
An open-loop system has at least one section where the water is in contact with
atmosphere air. Therefore, the pressure at any other point in the pipe network can be
determined by knowing pressure difference between the point and the air-water contact
surface.
A closed loop system is normally operated with the entire system above atmospheric
pressure to prevent air leakage into the system. Moreover, the total volume of closed loop
systems change with temperature due to the different coefficients of thermal expansion of
water and the pipe materials. An expansion tank is installed in close-loop systems to
accommodate such volume changes, and to provide a reference pressure to prevent
cavitation in the pump. The design procedures for expansion tanks is given elsewhere.
(a) direct-return system (b) reverse-return system.
Pipe network design
Guidelines are available in ASHRAE Handbook 2013 Fundamentals. The factors that
determine the maximum and minimum velocities in water pipes are the noise level, erosion
level, installation cost and operating cost.
If the pipes are too small, the resulting high velocities lead to unfavorable noise levels, erosion
levels and pumping costs. Too large pipes lead to excessive installation cost. For pipes of
nominal diameter 50 mm or less, the recommended velocity limit is 1.2 m/s. For pipe
diameters >50mm, a unit pressure loss limit of 400 Pa/m is recommended.
The main steps in the design of pipe networks are listed below:
1. Determine the layout of the pipe network with all the necessary fittings;
2. Size the pipes in all sections of the system based on the anticipated maximum flow rate;
3. Determine the head loss through the different loops of the pipe network;
4. Obtain the head requirement of the pump based on the largest head loss from step 3 above;
5. Adjust the pipe diameters of the branches so that all the loops in step 3 have the same head
loss as that used to selected the pump in step 4.
Example 12. Water at 70oC flows at the rate of 3.2 L/s through a horizontal steel pipe of
nominal diameter of 50 mm and length 45 m. The absolute roughness of commercial steel is
0.000046. The viscosity and density of water at 70oC are 0.406 mPas and 977.7 kg/m3.
Calculate (a) the total head loss through the pipe; (b) the velocity head, and (c) the change in
static head.
Solution: The actually inner diameter of the pipe with a nominal diameter 50 mm is 52.5 mm.
For steel pipe and the Reynolds number, 𝑓 = 0.0206.
4𝑄 4×3.2×10−3
Velocity 𝑉 = = 3.14× 52.5×10−3 2 = 1.479 m/s
𝜋𝐷2
𝑓𝐿 𝑉2 0.0206×45×1.4792
The friction head loss ∆𝐻𝑓𝑟 = 𝐷 2𝑔 = 52.5×10−3×2×9.81 = 1.97m;
𝑣2 1.4792
The constant velocity head in the pipe is 𝐻𝑣 = 2𝑔 = 2×9.81 = 0.111m;
For a straight pipe, the change in the total head is equal to the friction head loss,
𝑉12 𝑉22
As = , hence the change in static head is ∆𝐻𝑠𝑡 = ∆𝐻𝑓𝑟 = 1.97m.
2𝑔 2𝑔
Dyamic visocity 0.0000182 1.82E-05 1.82E-05
12 35 46
Length 45 40 70
Diameter 0.6 0.5 0.4
V 8.49 7.13 7.95
1/2*L/D*V^2 3243.6045 2440.171 6636.263
Re 3.36E+05 2.35E+05 2.10E+05
Epsilon 0.0003 0.0003 0.0003
epsilon/D 0.0005 0.0006 0.00075
f 0.0165 0.0175 0.0182
Delta P 53.519474 42.703 120.78
Re = 1.91x105
Example 13. Water at 20oC is distributed through the pipe network shown in Figure 13.1. The
lengths, diameters, and flow rates of the different pipe sections are indicated in the figure.
(a) Calculate the total pressure losses from 1-6, 1-8 and 1-0.
(b) If the pressures at the pump suction 0 and the exit sections 8, 6 and 10 are atmospheric,
estimate the required delivery pressure of the pump at 1;
(c) If the efficiency of the pump is 75%, calculate the required power input to the pump.
4𝑄𝑖
Solution: The water velocity in a pipe section: 𝑉𝑖 =
𝜋𝐷𝑖2
𝑓𝐿 𝑉2
The friction loss in the pipe section: ∆𝐻𝑓𝑟 =
𝐷 2𝑔
From Moody diagram – not required in this lecture
Table 2.1 Friction head losses in pipe sections
The loss coefficients for the fittings are obtained. Hence for the two tee-junctions
𝑘2,3 = 0.9; 𝑘2,7 = 1.23; 𝑘4,5 = 0.9; 𝑘4,9 = 1.4;
For the 90 elbows, in section 7-8, 𝐷 = 35mm, 𝑘 = 1.26; in section 9-10, 𝐷 = 35mm, 𝑘 = 1.26
(A) The corresponding velocity heads are the same as in Table 2.1. Hence
𝐻𝑡1 − 𝐻𝑡6 = ∆𝐻12 + ∆𝐻23 + ∆𝐻34 + ∆𝐻45 + ∆𝐻56 = 1.244 + 0.9 × 0.236 + 3.517 + 0.9 × 0.558 + 6.336 = 11.8m;
𝐻𝑡1 − 𝐻𝑡8 = ∆𝐻12 + ∆𝐻27 + ∆𝐻𝑒𝑙𝑏𝑜𝑤 + ∆𝐻78 = 1.244 + 1.23 × 0.236 + 1.26 × 0.495 + 9.809 = 11.97m;
𝐻𝑡1 − 𝐻𝑡8 = ∆𝐻12 + ∆𝐻23 + ∆𝐻34 + ∆𝐻49 + ∆𝐻𝑒𝑙𝑏𝑜𝑤 + ∆𝐻9,10 = 1.244 + 0.9 × 0.236 + 3.517 + 1.4 × 0.558 +
1.26 × 0.345 + 5.775 = 11.96m.
(b) The total head loss through the three sections are nearly equal. We can select the pump to deliver 9.5 L/s of water at a
total head of about 12 m. The ideal work input under these operating conditions is 𝑊𝑖𝑑𝑒𝑎𝑙 = ∆𝐻𝜌𝑔 𝑄 = 12 × 998 × 9.81 ×
9.5 × 10−3 = 1116W, and
𝑊 1116
(C) 𝑊𝑎𝑐𝑡 = 𝑖𝑑𝑒𝑎𝑙 = = 1488W.
𝑡 .75
Example 15. A water-cooled condenser of a water chiller plant rejects 145 kW of heat to the
atmosphere through a cooling tower as shown in Figure 15.1. The rise in temperature of the
cooling water in the condenser is 6oC. The total length of piping from the pump exit 1 to spray
header 2 of the cooling tower is 14 m. The length of piping from the cooling tower sump 3 to
the pump inlet 4 is 12 m.
The vertical distance between the spray header and the water surface in the sump is 1.2 m. The
pipe inlet in the sump has a square inlet. The head loss through the condenser is 2.2 m. The
head loss through the spray nozzles in the cooling tower is negligible. The pipe network uses
schedule-40 steel pipes with threaded fittings.
(a) Size the pipes for this system;
(b) Calculate the rise in total head required at the pump. Header
集水池
Pump exit 1
Pump inlet 4
Solution: First we calculate the water flow rate through the condenser.
𝑄 145
𝑄 = 𝑉𝜌𝑐𝑤 ∆𝑡; 𝑉 = 𝜌𝑐 ∆𝑡 = 1000×4.2×6 = 5.75 L/s
𝑤
Using the guideline from ASHRAE, using the friction chart, we choose the following pipe diameter
velocity and unit pressure loss for the pipe. The diameter satisfies the limit on velocity and unit
pressure loss.
Flow rate 5.75 L/s; diameter 70 mm, velocity 1.495 m/s, unit pressure loss 333 Pa/m
We then calculate the total pressure loss across each of the pipe sections. The relevant loss
coefficients for the fittings:
For 𝐷 = 70mm pipe, k-globe valve = 6.3, k-90 elbow = 0.83, k-square inlet = 0.5
The friction coefficient is 0.02089. Gravity acceleration 9.81 m/s2
The pressure in the cooling tower between 2 and 3 is atmospheric and 𝐻2 − 𝐻3 = 1.2m. We also
have
𝐻1 − 𝐻4 − 𝐻2 − 𝐻3 = 3.772 + 1.466 = 5.238m;
100m
Atm head in water
10m
DOI:10.1016/j.applthermaleng.2010.08.008