Subscriber access provided by University of Florida | Smathers Libraries
Article
Performance Evaluation of Two 25 kW Residential Wood Pellet Boiler Heating Systems Kui Wang, Mauro Masiol, Devraj Thimmaiah, Yuanyuan Zhang, and Philip K. Hopke Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b01868 • Publication Date (Web): 16 Oct 2017 Downloaded from http://pubs.acs.org on October 17, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
1
Performance Evaluation of Two 25 kW Residential
2
Wood Pellet Boiler Heating Systems
3
Kui Wang,1.2 Mauro Masiol,1 Devraj Thimmaiah, 1 Yuanyuan Zhang, 1 Philip K. Hopke1,2*
4
1. Center for Air Resources Engineering and Science
5
2. Department of Chemical and Biomolecular Engineering
6
Clarkson University
7
Potsdam, NY 13699 USA
8
ABSTRACT
9
A significant increase in the use of wood pellets for residential space heating has occurred over
10
the past decade. The performance of two modern residential wood pellet boilers (designated PB
11
and WPB) were evaluated including boiler thermal efficiency, thermal energy storage (TES) tank
12
discharge efficiency, and system efficiency. A correlation applicable to both systems between the
13
boiler thermal efficiency (ηth, in %) and the boiler output load (χ, in %) was found in the form of
14
ηth = 52.69 ln χ – 137.7 with R2 = 0.79 (for 25 < χ < 75). This equation provides an easy,
15
accurate estimation of the boiler thermal efficiency in field operations. The boiler thermal
16
efficiency decreased with time and this decline was determined using a Mann-Kendall trend
17
analysis with Sen’s slope. This decrease was primarily the result of fouling in the heat
18
exchanger and thus, this analysis identifies the need for manual cleaning of the heat exchanger
19
tubes to restore maximal system performance. The evaluation of the TES tank performance
20
found that the TES tank discharge efficiency was correlated with a dimensionless function of
21
tank inlet Reynolds number (Red) and temperature differences in the tank and inlet and outlet
22
pipes. Overall system efficiency showed a seasonal average of 62.8%, 62.0%, and 75.8% for
23
three heating seasons of the PB system. These results provide a comprehensive performance
24
evaluation of these wood pellet boiler heating systems in the field over an extended period of
25
operation.
26 27
Keywords: wood pellet boiler; thermal efficiency; Sen’s slope; TES tank performance
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
28 29 30
Page 2 of 23
Nomenclature Btotal: energy to the tank from boiler, kJ
∆Td:
cp,w: water specific heat, kJ/(kg·K) Echa: energy to building during tank charge, kJ Edis:
energy to building discharge, kJ
during
tank
HHV: pellet fuel higher heating value, MJ/kg
TES tank inlet/outlet temperature difference, °C
∆Td,ave: average TES tank inlet/outlet temperature difference during discharge, °C ∆Ts: TES tank top and bottom temperature difference, °C
Mw: pellet fuel moisture content (w.b.)
∆Ts,ave: average TES tank top and bottom temperature difference during discharge, °C
N: boiler nominal load, kW
ε: boiler stoker feed fraction
Q: boiler loop flow rate, m3/h
ηth: boiler thermal efficiency, %
Q*: building loop flow rate, m3/h
ηsys: system efficiency, %
Red: inlet pipe Reynolds number
κ100: boiler stoker feed rate at 100%, kg/h
tc: time of a complete burn cycle, h
ξ: TES tank flow factor
td: boiler off time or tank discharge time, h
ρw: water density, kg/m3
T1: boiler output water temperature, °C
φ: TES tank overall efficiency (cycle basis), %
T1*: TES tank output water temperature, °C T2: boiler return water temperature, °C T2*: TES tank return water temperature, °C
φdis: TES tank discharge efficiency, % χ: boiler output load, %
∆t: boiler flow meter recording interval, s
ACS Paragon Plus Environment
Page 3 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
31 32
Energy & Fuels
1. Introduction Because of the increasing interest in the use of renewable energy, biomass has been
33
increasingly utilized over the past decade. According to IEA, bioenergy provides 10% of world
34
primary energy supply currently1 and biofuels can provide up to 27% of world transportation fuel
35
by 2050.2 Generally, biomass refers to any material that is derived directly or indirectly from
36
plants.3 Biomass, especially wood and wood pellets, is becoming more widely used as an
37
alternative energy source for domestic space heating due to its relatively high energy density and
38
uniform quality compared to cordwood or chips.4 In the European Union, wood pellet use is
39
growing rapidly due to government policies to meet renewable energy targets.3,4,5
40
Wood pellets offer both economic and climate benefits. Wood can be harvested sustainably
41
and provide a substantial local energy supply as well as offering extra employment to collect and
42
process the wood. Chau et al.6 provided a detailed techno-economic evaluation of a small scale
43
wood pellet heating system. Wood pellet boiler technology has also improved over the past
44
decades particularly driven by European interest in utilizing a local, renewable fuel. Modern
45
wood pellet boilers have high efficiency and low pollutant emissions. Fiedler7 reviewed the state
46
of the art of small-scale pellet combustion units for domestic use and the individual EU country
47
regulations where more stringent regulations are likely to appear in the future. Öhman et al.8
48
studied the slagging tendencies in residential wood pellet boilers related to ash sintering
49
temperatures. They found that wood pellets have a low sintering tendency (sintering
50
temperatures higher than 1100 °C) mainly due to the low alkali and Si content in the fuel.9 Hence,
51
modern wood pellet boilers utilize moving grates or sliding grate ash removal systems.7 Better
52
combustion chamber and flow geometry can be designed and optimized by applying
53
computational fluid dynamics (CFD) simulations.10,11 In addition, studies have shown that when
54
firing with non-woody biomass pellets, better boiler control such as fuel feeding, air intake, and
55
ash management is required.12
56
A key index for assessing the performance of wood pellet boiler is its thermal efficiency.
57
There are many factors affecting the thermal efficiency of wood pellet boilers. Chandrasekaran et
58
al. characterized the efficiency of two wood pellet boilers13 and wood stoves burning multiple
59
types of pellets14. Efficiencies ranged from 70% to 91% for different boilers running with wood
60
pellets and from 75% to 85% when burning with non-woody pellets. Verma et al.15,16 found that
61
the performance of a domestic pellet boiler is a function of the boiler operational loads. Carlon
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
62
et al.17 reported a non-linear relationship between thermal efficiency and a modified load factor.
63
Several factors lead to the efficiency-load relationships. The most important factor is the excess
64
air ratio that is usually higher under low load conditions. In addition, low load conditions cannot
65
generate sufficient turbulence for complete combustion, which leads to low combustion
66
efficiency and temperatures.16 Heat transfer is also reduced during low load conditions. Thus, it
67
is better for the boiler to operate at high or full load condition so that maximum thermal
68
efficiency can be achieved.
69
Depending on the building load and the ambient temperature, the building may not require all
70
of the heat produced by the boiler when it operates at full load. The use of thermal energy
71
storage (TES) enables the wood pellet boiler to operate for longer periods at higher loads and,
72
consequently, reduces the number of boiler ignitions and avoids short-cycling during warmer
73
weather. TES in the form of large insulated water tanks was initially utilized in solar water
74
heating systems because the stratification inside the tank can help to extract more useful heat
75
from the tank.18 Many experimental and numerical studies have been performed to study the
76
performance of TES tank in terms of temperature stratification, efficiency, and operational
77
behavior.19-25
78
The present work systematically evaluated the performance of two modern wood pellet
79
heating systems in terms of: boiler thermal efficiency at various loads, TES tank discharge
80
efficiency and overall system efficiency based on long-term monitoring from 2015 to 2017.
81
Currently, there are no such systematic evaluations of in-use wood pellet heating systems for
82
such an extended period.
83 84
2. Material and Methods
85
2.1 The Two Systems
86
In this project, two nominally identical 25 kW wood pellet boilers were installed in two
87
homes in northern New York State with different configurations (as shown in Figure 1). In each
88
boiler, wood pellets are fed by a stoker screw auger into the combustion chamber from the top of
89
the burner pot. A hot air gun ignites the pellets. Primary air enters the combustion chamber from
90
the bottom of the burner pot. Hot combustion gas passes through a 3-stage fire tube heat
91
exchanger to transfer heat into jacket water and exhausts through the chimney. The induced draft
92
fan at the end of the flue gas pathway controls the primary air flow through a Lambda sensor (O2
ACS Paragon Plus Environment
Page 4 of 23
Page 5 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
93
control) (calibrated with an external combustion analyzer) and a negative pressure sensor (draft
94
control). Within each fire tube, there is a twisted turbulator to improve heat transfer as well as to
95
clean the inside surfaces of the fire tubes by moving the turbulators up and down. The residual
96
ash is removed from the bottom ash collector. The boiler is controlled by software that records
97
its operational data.
98 99
2.2 Pellet Fuel Characterization
100
The pellets used in this study are PFI (Pellet Fuels Institute26) certified premium pellets.
101
Samples were collected during bulk delivery to the pellet bins. The pellets were characterized
102
for heating value, moisture, composition, and ash content. Detailed analytical methods and
103
results are listed in Table 1. The pellets have a mean calorific content of 19.12 MJ/kg (or 8221
104
Btu/lb) higher heating value (HHV), 3.88% moisture, 0.58% ash, and 46.64% carbon.
105 106
Table 1. Summary of analytical methods and results Parameter
Method
Result
HHV
ASTM E 711-87 (2004)
19.12 MJ/kg (a.r.)
Moisture
ASTM E 871-82 (2006)
3.88% (w.b.)
Ash
ASTM D 1102-84 (2007)
0.58% (d.b.)
Carbon
ASTM E 777-87 (2004)
46.64% (d.b.)
Hydrogen
ASTM E 777-87 (2004)
6.4% (d.b.)
Nitrogen
ASTM E 778-87 (1996)
0.14% (d.b.)
Oxygen
By difference
46.23% (d.b.)
Chlorine
ASTME 776-87 (2009)
0.00388% (d.b.)
Sulfur
ASTM E 775-87 (2004)
0.00689% (d.b.)
107 108
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
109 110
2.3 Experimental Methods Figure 1 shows the two system configurations as operated during the 2016-17 heating season.
111
These systems had each undergone modification over time and detailed histories of each system
112
are presented in the Supporting Information. The PB (Parishville boiler) system was installed in
113
a modified shipping container near a family house in Parishville NY, USA (44.63 N, 74.83 W,
114
271 m a.s.l.) replacing their existing outdoor wood boiler (OWB). Hot water from the boiler was
115
stored in two identical TES tanks (A.O. Smith Co., volume of each = 0.45 m3) and then delivered
116
to the house for space heating through an underground circulation system. The tanks provided
117
36.0 liter per kW (2.8 gallons per kBtu/h) storage, above manufacturer recommendations of 25.7
118
liter per kW (2 gallons per kBtu/h). Return water from the building flowed to the bottom of the
119
tanks. This house has a total heating area of 406.5 m2 with an estimated peak load of 27 kW
120
(92000 Btu/h). The tanks store hot water from the boiler at temperature ranging from 55°C
121
(131°F) to 75°C (167°F) that is used for a radiant floor heating system with return water
122
temperatures ranging from 35°C (95°F) to 45°C (113°F), depending on outdoor temperature and
123
the thermostat settings.
124
The other nominal 25 kW boiler was installed in the basement of a family house located in
125
West Potsdam, NY, USA (44.69 N, 75.09 W, 122 m a.s.l.). An existing oil boiler was integrated
126
into the system as an auxiliary boiler. This system is denoted as WPB (West Potsdam boiler)
127
system. A single 210-gallon TES tank (0.79 m3) was installed in the system. In this system, the
128
TES tank was installed in a 2-pipe configuration such that both the boiler and the heat demand
129
are on the same side of the tank. In addition, WPB system had a 0.1 m3 domestic hot water
130
(DHW) tank that provided domestic hot water. The additional oil boiler and the DHW tank
131
made both the system control and the boiler behavior different from PB system.
132
For both systems, flow and temperature were monitored in two loops: boiler to tank loop and
133
tank to building loop to calculate heat output from the boiler and the TES tank, respectively. A
134
certified ultrasonic flow meter (EF-10 series, Spire Metering Technology) and two temperature
135
probes were installed in each loop. The boiler was controlled by two temperature sensors in the
136
TES tank (bottom ON and top OFF).
137
To determine the energy input rate, the boiler stoker feed auger was calibrated at several feed
138
rates. The pellets were collected in a weighed bucket and after a fixed time, the bucket was
139
weighed. From these results, the maximum feed rates (κ100) were determined to be 12.98 kg/h
ACS Paragon Plus Environment
Page 6 of 23
Page 7 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
140
and 12.17 kg/h for PB and WPB, respectively. The constant speed augers are modulated by
141
operating them only a part of the time (feed fraction, ε). The actual feed fractions are recorded at
142
1 second intervals by each boiler. Using the fuel characterization data, the heat input into the
143
boiler can be calculated as the product of the mass flow rate and the calorific content of the
144
pellets. It was found that the augers were quite reproducible over time such that multiple
145
calibrations produced essentially identical relationships between time of auger operation and
146
mass flow rates of fuel delivered into the boiler.
147
The operation of these boiler-thermal storage systems was characterized in terms of “system
148
cycles” that began with the start-up of the boiler, steady-state combustion, boiler shut down, and
149
finally depletion of heat from the thermal storage system to the building until the boiler restarted
150
because the control temperature at the top of the thermal storage tank had been reduced to the
151
setpoint value. Within these system cycles, there are “boiler cycles.” A complete boiler cycle
152
includes three phases: start-up, steady state combustion, and burnout/shutdown.
153
For PB system, measurements were made over three heating seasons denoted as PB-1, PB-2,
154
and PB-3. PB-1 included data from February 1, 2015 to April 11, 2015 representing 493 cycles;
155
PB-2 included data from December 1, 2015 to April 1, 2016 with a total of 392 cycles; and PB-3
156
included data from January 5, 2017 to April 5, 2017 with a total of 372 cycles. During PB-1 and
157
PB-2, the thermal storage tanks were continuously connected to the circulation loop that
158
conducted heat to the building (See Supplemental Figure SI-3). Water was continuously
159
circulated through this loop to avoid freezing in the underground line. Before PB-3, a 3-way
160
valve was installed (Figure 1, top) such that the circulation loop was separated from the thermal
161
storage tanks when there was no demand for heat from the house. Local weather data obtained
162
from Weather Underground27 were used in this study.
163
For WPB system, major changes had been made over the prior heating seasons to address
164
operational problems that had been identified (see supplemental information). Thus, only WPB-
165
3 data were used for the performance analysis. WPB-3 ran from December 28, 2016 to April 12,
166
2017 and included a total of 644 cycles.
167 168 169 170
2.4 Efficiency Calculations Because of the slower start-up and shut down associated with solid fuel combustion, the boiler thermal efficiency was calculated based on each boiler cycle. Overall efficiency was
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 23
171
calculated based on system cycles. For each boiler cycle, the boiler thermal efficiency (ηth) was
172
calculated as the ratio of heat output to heat input over each cycle time (tc) =
∑ , ( – ) ∙ Δ ∑(HHV ∙ ∙ )
× 100%
(1)
173
where ∆t is ultrasonic flowmeter data sampling interval (10 s). The boiler stoker feed fraction ε
174
was recorded every 1 s.
175
System efficiency is defined based on system cycles. It is the total energy delivered into the
176
building during a system cycle divided by the total energy input into the boiler for each boiler
177
cycle, as #$# =
178
∑&' ∗ , ( ∗ – ∗ ) ∙ Δ ∑(HHV ∙ ∙ )
× 100%
where td is the time when boiler is off until next boiler cycle.
ACS Paragon Plus Environment
(2)
Page 9 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
179
180 181
Figure 1. Schematic of the PB (top) and WPB (bottom) systems with data acquisition devices (S-
182
1: flue gas O2 control signal, S-2: boiler pump control signal; S-3: Boiler ON signal; S-4: Boiler
183
OFF signal)
184
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
185
3. Results and Discussion
186
3.1 Boiler Thermal Efficiency
187
Page 10 of 23
During a boiler cycle, the instantaneous efficiency of the boiler changes. Low efficiency
188
occurs during start-up due to inefficient combustion during the ignition and stabilization stages
189
when there are also high CO and particle emissions.28 As the combustion stabilizes, the boiler
190
thermal efficiency reaches its highest value during its steady state operation. Low efficiencies
191
again occur as the boiler shuts down and the residual fuel burns out. Thus, boiler thermal
192
efficiency over each cycle is an appropriate metric for boiler performance.
193
Figure 2 shows the regression between boiler thermal efficiency and boiler output load
194
percentage for both boiler systems. For the PB system, due to boiler software updates and
195
frequent setting change of TES tank ON/OFF temperatures during PB-2 and PB-3, a relatively
196
low r2 = 0.68 was determined. However, when these parameters remained fixed in PB-1, the r2
197
rose to 0.86. For the WPB system, the operation of DHW significantly affected the boiler cycling
198
characteristics, as explained in the supporting information. Because the DHW tank represented a
199
low thermal mass, very little energy was dissipated into the domestic hot water. When the boiler
200
was ON because of demand from the thermal storage tank and the DHW was also demanding
201
heat, the DHW tank acts as an independent TES tank. However, the low heat dissipation in the
202
DHW tank resulted in very high boiler return water temperatures that approached the safety limit
203
set temperature of the boiler. Therefore, the boiler reduced its output quickly sometimes causing
204
boiler short cycling. The problem was reduced during WPB-3 by defining wider temperature set-
205
points for DHW tank but controls alone could not eliminate the problem.
206 207
A general regression equation based on four datasets for the two systems was determined by merging the data and fitting an equation as follows: ηth = 52.69 ln χ – 137.7
(for 25< χ i), n is the number of total data points, and sgn(xj - xi ) is the
272
sign function defined as
273
+1, sgn234 − 3 6 = : 0, −1,
(7)
Mann30 and Kendall31 have shown that when n ≥ 18, S is approximately a normal distribution
274 275
34 − 3 > 0 34 − 3 = 0 34 − 3 < 0
with expectation of 0 and variance as >
(-)
?(? − 1)(2? + 5) − ∑B 8 ( − 1)(2 + 5) = 18
ACS Paragon Plus Environment
(8)
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 23
276
where m is the number of tied groups and ti denotes the number of ties of extent i. A tied group is
277
defined as a set of samples with the same value. The standard normal test (or Z-test) statistic ZS
278
is computed as
279
-−1 , if - > 0 I GJ> (-) DE = 0, if - = 0 H -+1 G , if - < 0 FJ> (-)
280
usually used to estimate the slope of the trend (also called Sen’s slope analysis). Figure 4 shows
281
the analysis results for PB and WPB systems.
282
(9)
ZS > 0 means increasing trends while ZS < 0 means decreasing trends. A Theil-Sen’s line is
A decreasing trend in boiler thermal efficiency can be observed in Figure 4(a) when there
283
was no manual cleaning performed. High boiler thermal efficiency of around 85% was initially
284
observed after the boiler installation. An average rate of -0.022 % per boiler cycle was
285
determined for PB-1. Therefore, if the boiler must be cleaned manually when the thermal
286
efficiency has dropped by 10%, PB has to be cleaned after 455 boiler cycles. Given an average
287
of 7 boiler cycles per day, the boiler needs to be cleaned after ~65 days of continuous operation.
288
Figure 4(b) shows a relatively constant trend of boiler thermal efficiency for PB-2, which
289
suggests that the boiler can only maintain around 70% thermal efficiency based on its internal
290
cleaning mechanism. After two manual cleanings were performed during PB-3, the boiler
291
thermal efficiency significantly improved compared to PB-2. For WPB-3, an average boiler
292
thermal efficiency rate of -0.038% per boiler cycle was determined, which means WPB must be
293
manually cleaned after 263 boiler cycles or ~38 days of continuous operation given the same
294
evaluating criterion as PB. WPB-3 had the higher rate of decline in the boiler thermal efficiency.
295
One factor for this boiler could be the high moisture content in the basement. This old basement
296
has little ventilation and standing water was often observed on the basement floor on rainy days.
297
When the boiler was manually cleaned at the beginning of WPB-3, the cleaning turbulators had
298
to be replaced due to substantial fouling and corrosion. In a high moisture environment, it is
299
much easier to have condensation in the boiler, which facilitates boiler fouling and corrosion
300
ACS Paragon Plus Environment
Page 15 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
(a). PB-1
(b). PB-2
(c). PB-3
(d). WPB-3
301 302
Figure 4. Sen’s slope of boiler thermal efficiency for PB and WPB systems (c.i. – confidence
303
interval)
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
304 305
process. Thus, operation of a dehumidifier in these circumstances may improve the boiler performance.
306 307
3.3 Thermal Storage Tank Performance
308
The TES tank top and bottom temperatures are measured by the ON and OFF temperature
309
sensors inserted into sensor wells in the tank as shown in Figure 1. The TES tank top and bottom
310
temperature difference (∆Ts) can be used to assess tank stratification since a well stratified tank
311
has a larger ∆Ts. Figure 5 shows tank behavior during a complete system cycle in PB-1.
312 313
Figure 5. TES tank temperature distribution during one complete boiler burn cycle in PB-3
314 315
The vertical dashed lines divide the system cycle into charge and discharge periods. During
316
the TES tank charging period, the boiler heats the tank to the set temperature. Then the boiler
317
shuts down, and the tank discharges heat to the building. The amount of energy discharged is
318
determined by the stratification inside the tank since a well stratified tank has higher exergy than
ACS Paragon Plus Environment
Page 16 of 23
Page 17 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
319
a completely mixed tank.32 The TES tank can continuously supply energy to the building without
320
interruption depending on the building demand.
321 322
A well stratified tank can discharge more energy than a completely mixed tank, as shown in earlier studies.19-24 The overall tank efficiency (φ) is defined as φ = (Echa + Edis)/Btotal × 100%
323
(10)
Simultaneously, tank discharge efficiency (φdis) is defined as φdis = Edis/(Btotal – Echa) × 100%
(11)
324
where Echa and Edis are energy delivered to the building during charge and discharge periods,
325
respectively; Btotal is the total energy delivered to the tank during a complete boiler cycle.
326
Tank discharge efficiency, φdis, is a good indicator of the stratification inside the tank when
327
the boiler is off. Studies have shown that two major factors affect the stratification: tank
328
geometry (inlet/outlet positions,20 tank aspect ratio,19 baffle plates,23 etc.) and operational
329
conditions (temperature difference,33 flow rates,33 operation modes,21,22 etc.). The present
330
evaluation focuses on these operational conditions because the tanks had a fixed geometry during
331
the entire monitoring period and the TES tank was operated with continuous flow on the building
332
side for PB system during the first two heating seasons.
333
Figure 6(a) shows TES tank top and bottom average temperature difference (∆Ts,ave) over
334
each discharge period as a function of corresponding discharge time td (or boiler off time) for
335
PB-1 data. The correlation equation is shown below (R2 = 0.83) ∆Ts,ave = 66.75(td)-1 – 150.0(td)-2 + 119.3(td)-3 – 5.965
ACS Paragon Plus Environment
(12)
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 23
336 337
Figure 6. (a): TES tank top and bottom average temperature difference as a function of discharge
338
time for PB-1; (b): Tank discharge efficiency correlated with dimensionless flow factor for PB-1
339 340
Longer discharge times indicate less heat demand from the building (i.e. warmer weather)
341
that result in smaller temperature differences in the tank since the TES tank supply and return
342
water temperatures will be closer in temperature. Then, no clear thermocline may develop in the
343
tank.34 Consequently, when the building is demanding heat, a larger ∆Ts,ave will occur due to the
344
larger temperature difference between the TES tank supply and return water temperatures. Figure
345
6(a) shows that stratification of a TES tank in the heating system is highly dependent on the
346
building heat load. This result suggests that maintaining stratification inside the tank at low
347
building heat loads can be attained by decreasing the flow rate to create a larger return water
348
temperature difference from the tank top temperature. Figure 6(b) shows TES tank discharge efficiency φdis is correlated with a dimensionless
349 350
number ξ by the following equation (r2 = 0.73) ln φdis = 8.3 – 175.3ξ
(13)
351
where φdis is in %, and ξ is expressed in pipe Reynolds number (Red) and temperature differences
352
as ξ = (∆Ts,ave /∆Td,ave)0.35(Red)-0.41
(14)
353
where ∆Ts,ave is the average TES tank top and bottom temperature difference during discharge
354
period and ∆Td,ave is the average TES tank supply and return temperature difference during the
355
same discharge period. Figure 6(b) indicates that an increase in Red and decrease in ∆Ts,ave
ACS Paragon Plus Environment
Page 19 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
356
/∆Td,ave will lead to an exponential increase in tank discharge efficiency. For a specific TES tank,
357
increasing the Red means an increased flow rate, which corresponds with previous studies that
358
the TES tank discharge efficiency (or extraction efficiency) increases with increasing flow
359
rate.19,24,33 Thus, inlet flow diffusers and baffles23,24,35 are used to minimize flow mixing and
360
maintain a high level of stratification. The ratio of ∆Ts,ave /∆Td,ave is an indicator of the amount of
361
heat delivered to the building. The maximum value of ∆Ts,ave /∆Td,ave is 1 when there is no heat
362
delivered to the building or the return water temperature is equal to the temperature at the bottom
363
of TES tank. Small ratios suggest large TES tank supply and return temperature differences,
364
which further indicate large amounts of heat consumed by the building.
365
For PB-2 and PB-3, correlations between ∆Ts,ave and td, φdis and ξ have much lower R2 values
366
compared to PB-1 data as shown in Figure 6. This result is caused by the constant TES tank
367
ON/OFF setting maintained over the whole heating season for PB-1, but a variety of TES tank
368
ON/OFF settings were explored during PB-2 and PB-3 that in turn reduced the correlation. For
369
WPB-3, the TES tank stratification (∆Ts) was constantly disturbed by the DHW tank operation.
370
Therefore, the correlation was worse than for PB-2 and PB-3. In addition, the 2-pipe
371
configuration prohibits the determination of the exact amount of energy delivered into the TES
372
tank since part of the energy produced while the boiler is on can be directed to the DHW tank.
373
Therefore, no correlation for TES tank discharge efficiency can be established. An additional
374
flowmeter on the DHW tank loop would have been required to permit this calculation.
375
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
376 377
3.4 System Efficiency System efficiency evaluates the efficiency of the whole system based on the heat delivered to
378
the building and the heat input into the boiler. The average seasonal system efficiency values
379
were 62.8%, 62.0%, and 75.8% for PB-1, PB-2, and PB-3, respectively. Low system efficiencies
380
during PB-1 and PB-2 was mostly due to the low thermal efficiencies and low TES tank
381
discharge efficiencies when there was no 3-way valve installed to prevent the constant
382
circulation of the building circulator as shown in Figure 1 (as discussed in supplemental
383
information). After effective boiler cleaning and system modification, the system efficiency in
384
PB-3 significantly increased. System efficiencies for WPB could not calculated because the
385
current measurement system did not permit the determination of the amounts of heat input
386
transferred into both the TES tank and the building based on the 2-pipe configuration system.
387 388
4. Conclusions and Suggestions
389
4.1 Conclusions
390
The objective of this study was to evaluate the performance of two 25 kW residential wood
391
pellet boilers in terms of boiler thermal efficiency, thermal energy storage (TES) tank discharge
392
efficiency, and system efficiency based on long term field site monitoring results. A general
393
correlation between boiler thermal efficiency (ηth, in %) and boiler output load percentage (χ,
394
in %) was determined to be ηth = 52.69 ln χ – 137.7 with r2 = 0.79 (for 25 < χ < 75) for both
395
systems. Evaluation of thermal efficiency bins showed that most of time, the PB boiler was
396
operating above 70% efficiency with seasonal average thermal efficiencies of 73.6%, 70.8%, and
397
79.0% for PB-1, PB-2, and PB-3, respectively. However, the WPB boiler was operating below
398
70% thermal efficiency most of the time caused by the DHW induced load modulation and
399
severe boiler fouling. A Sen’s slope analysis on the boiler thermal efficiency can determine the
400
boiler thermal efficiency decreasing rate and provide information on when additional cleaning is
401
needed. The results showed that the PB boiler had a thermal efficiency drop of 0.022% per boiler
402
cycle and could maintain around 70% boiler thermal efficiency without manual cleaning. The
403
WPB boiler had a thermal efficiency drop of 0.038% per boiler cycle and must be cleaned every
404
38 days (given an average of 7 boiler cycles per day). TES tank performance showed a
405
correlation of the tank discharge efficiency with a dimensionless number that is a function of
406
inlet pipe Reynolds number and temperature differences in both the tank and pipe inlet/out. The
ACS Paragon Plus Environment
Page 20 of 23
Page 21 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
407
system efficiency evaluation revealed that effective boiler cleaning and system modification can
408
enhance overall system performance. Overall system efficiency showed seasonal averages of
409
62.8%, 62.0%, and 75.8% for PB-1, PB-2, and PB-3, respectively. Future work will be focused
410
on continuing monitoring and performing dynamic process simulation of the system for further
411
optimization. The boiler emissions will be discussed separately in a future manuscript.
412 413 414
4.2 Suggestions for Improvement (1). High moisture content in the combustion air can accelerate wood pellet boiler fouling, as
415
shown in WPB installation. Therefore, it is important to ensure a low moisture content
416
environment for the wood pellet boiler such as adding a dehumidifier to high moisture content
417
environment.
418
(2). The 2-pipe configuration does not necessarily perform better than the traditional 4-pipe
419
configuration. It is important to recognize the effect of the thermal mass of the different heat
420
demanding units connected to the system such that the boiler short-cycling and load modulation
421
shown in WPB system can be avoided. An alternative 3-pipe configuration might solve the
422
problem.
423
(3). Frequent manual cleaning of the heat exchanger is essential to maintaining high overall
424
system performance for these wood pellet heating systems. The decline in the thermal efficiency
425
can serve as an indicator of the need for additional cleaning.
426 427
Acknowledgements
428
This work was supported by the New York State Energy Research and Development Authority
429
(NYSERDA) under agreements #32977 and 63007. We want to express our thanks to the
430
homeowners of the boiler systems for their willingness to work with us to assess the
431
effectiveness of this technology. John Siegenthaler of Appropriate Designs and Khaled Yousef
432
of Pyramid Energy Engineering Services, PLLC are acknowledged for their advice on system
433
design and modifications. We also want to acknowledge the cooperation of the boiler
434
manufacturer in helping to address issues as they were identified.
435 436 437
ACS Paragon Plus Environment
Energy & Fuels
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
438
Page 22 of 23
References (1) International Energy Agency (IEA), available at https://www.iea.org/topics/renewables/subtopics/bioenergy/. (2) International Energy Agency (IEA), available at https://www.iea.org/newsroomandevents/pressreleases/2011/april/biofuels-can-provide-up-to27of-world-transportation-fuel-by-2050-iea-report-.html. (3) Demirbas, A. Biomass resource facilities and biomass conversion processing for fuels and chemicals. Energy Convers Manage 2001, 42, 1357-78. (4) Goh, C.S.; Junginger, M.; Cocchi M.; Marchal, D.; Thrän, D.; Hennig, C.; et al. Wood pellet market and trade: a global perspective. Biofuels Bioprod Biorefin 2013, 7, 24-42. (5) Parikka, M. Global biomass fuel resources. Biomass Bioenergy 2004, 27, 613-20. (6) Chau, J.; Sowlati, T.; Sokhansanj, S.; Preto, F.; Melin, S.; Bi, X. Techno-economic analysis of wood biomass boilers for the greenhouse industry. Appl Energy 2009, 86, 364-71. (7) Fiedler, F. The state of the art of small-scale pellet-based heating systems and relevant regulations in Sweden, Austria and Germany. Renew Sustain Energy Rev 2004, 8, 201-21. (8) Öhman, M.; Boman, C.; Hedman, H.; Nordin, A.; Boström, D. Slagging tendencies of wood pellet ash during combustion in residential pellet burners. Biomass Bioenergy 2004, 27, 585-96. (9) Vega-Nieva, D.J.; Torres, L.; Tabares, J.; Morán J. Measuring and Predicting the Slagging of Woody and Herbaceous Mediterranean Biomass Fuels on a Domestic Pellet Boiler. Energy Fuels 2016, 30, 1085-95. (10) Porteiro, J.; Collazo, J.; Patiño, D.; Granada, E.; Gonzalez, J.C.M.; Míguez, J.L. Numerical modeling of a biomass pellet domestic boiler. Energy Fuels 2009, 23, 1067-75. (11) Chaney, J.; Liu, H.; Li, J. An overview of CFD modelling of small-scale fixed-bed biomass pellet boilers with preliminary results from a simplified approach. Energy Convers Manage 2012, 63, 149-56. (12) Carvalho, L.; Wopienka, E.; Pointner, C.; Lundgren, J.; Verma, V.K.; Haslinger, W.; et al. Performance of a pellet boiler fired with agricultural fuels. Appl Energy 2013, 104, 286-96. (13) Chandrasekaran, S.R.; Laing, J.R.; Holsen, T.M.; Raja, S.; Hopke, P.K. Emission characterization and efficiency measurements of high-efficiency wood boilers. Energy Fuels 2011, 25, 5015-21. (14) Chandrasekaran, S.R.; Hopke, P.K.; Newtown, M.; Hurlbut, A. Residential-scale biomass boiler emissions and efficiency characterization for several fuels. Energy Fuels 2013, 27, 4840-9. (15) Verma, V.K.; Bram, S.; Gauthier, G.; Ruyck, J.D. Performance of a domestic pellet boiler as a function of operational loads: Part-2. Biomass Bioenergy 2011, 35, 272-9. (16) Verma, V.K.; Bram, S.; Delattin, F.; Ruyck, J.D. Real life performance of domestic pellet boiler technologies as a function of operational loads: A case study of Belgium. Appl Energy 2013, 101, 357-62. (17) Carlon, E.; Schwarz, M.; Golicza, L.; Verma, V.K.; Prada, A.; Baratieri, M.; et al. Efficiency and operational behaviour of small-scale pellet boilers installed in residential buildings. Appl Energy 2015, 155, 854-65. (18) Hollands, K.G.T.; Lightstone, M.F. A review of low-flow, stratified-tank solar water heating systems. Sol Energy 1989, 43, 97-105. (19) Lavan, Z.; Thompson, J. Experimental study of thermally stratified hot water storage tanks. Sol Energy 1977, 19, 519-24.
ACS Paragon Plus Environment
Page 23 of 23
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
(20) Shah, L.J.; Furbo, S. Entrance effects in solar storage tanks. Sol Energy 2003, 75, 337-48. (21) Fernandez-Seara, J.; Uhia, F.J.; Sieres, J. Experimental analysis of a domestic electric hot water storage tank. Part I: Static mode of operation. Appl Therm Eng 2007, 27,129-36. (22) Fernandez-Seara, J.; Uhia, F.J.; Sieres, J. Experimental analysis of a domestic electric hot water storage tank. Part II: dynamic mode of operation. Appl Therm Eng 2007, 27,137-44. (23) Levers, S.; Lin, W. Numerical simulation of three-dimensional flow dynamics in a hot water storage tank. Appl Energy 2009, 86, 2604-14. (24) Han, Y.M.; Wang, R.Z.; Dai, Y.J. Thermal stratification within the water tank. Renew Sustain Energy Rev 2009, 13, 1014-26. (25) Rosen, M.A.; Dincer, I. Exergy methods for assessing and comparing thermal storage systems. Int J Energy Res 2003, 27, 415-30. (26) Pellet Fuels Institute, available at http://www.pelletheat.org/joining-the-standards-program. (27) Weather History & Data Archive, available at https://www.wunderground.com. (28) Win, K.; Persson, T. Emissions from residential wood pellet boilers and stove characterized into start-up, steady operation, and stop emissions. Energy Fuels 2014, 28, 2496-505. (29) Sippula, O.; Hytönen, K.; Tissari, J.; Raunemaa, T; Jokiniemi, J. Effect of wood fuel on the emissions from a top-feed pellet stove. Energy Fuels 2007, 21, 1151-60. (30) Mann, H.B. Nonparametric tests against trend. Econometrica, 1945, 13, 245–259. (31) Kendall, M.G. Rank Correlation Methods. 1975, Griffin, London, UK. (32) Rosen, M. A. The exergy of stratified thermal energy storages. Sol Energy, 2001, 71, 173185. (33) Hariharan, K.; Murthy, M.V.K. Temperature stratification in hot-water storage tanks. Energy 1991, 16, 977-82. (34) Abdoly, M.A.; Rapp, D. Theoretical and experimental studies of stratified thermocline storage of hot water. Energy Convers Manage 1982, 22, 275-85. (35) Altuntop, N.; Arslan, M.; Ozceyhan, V.; Kanoglu, M. Effect of obstacles on thermal stratification in hot water storage tanks. Appl Therm Eng 2005, 25, 2285-98.
ACS Paragon Plus Environment