Traditional fossil fuel power generation process typically has low efficiency. Large amount of the energy loss in Rankine cycle steam turbines (ST) is due to the temperature difference between the combustion flame temperature ∼2250 K (adiabatic) and the high pressure steam temperature up to 900 K. However, some of this energy can be harvested using solid-state thermoelectric (TE) power generators which are placed into the gap between the flame temperature and the steam temperature that produce additional electrical power. This study investigates the potential placement of TE on water tube wall inside a boiler at a coal-fired power plant. Three-dimensional (3D) numerical model of a simplified TE module is developed, and hot gas temperature and steam temperature from the boiler are used as boundary conditions at the hot side and cold side of the TE. The numerical results are compared with analytical calculations. The 3D effects of the thermal spreading in the TE module are investigated. Parameters such as TE leg cross section area and TE fill factor are examined in order to maximize the electrical power production of the TE without sacrificing the boiler efficiency (i.e., reducing the steam temperature). The study also looks into the various locations inside the boiler that have good potential for TE installation.

Introduction

TE generators utilize TE effect to directly convert thermal energy into electric. TE generators produce voltage when there is a temperature difference between the hot side and the cold side of TE generators.

TE generators have been receiving relatively a lot of attention for the waste heat recovery applications. Masahide et al. [1], Yodovard [2], Hsu et al. [3], and Kumar et al. [4] applied the TE technology to recover heat from vehicle exhaust and engine and convert it into electrical power to provide auxiliary power for the vehicle. TE generators have also been used to generate electricity from solar power [58]. Chen et al. [9] presented an analysis of system efficiency and evaluation of impact of TE power cycle application to thermal energy systems.

Yazawa et al. [10] proposed adding a TE power generator into the gap between the flame temperature and the steam temperature. Yazawa et al. developed a generic model and analyzed energy economy for a combined TE generator on top of a ST cycle. The temperature or this application will be much higher than those of the waste heat recovery applications which have been heavily investigated. The TE materials can be tuned for the targeted high temperatures. Yazawa and Shakouri [11] reported that the TE module can be designed for optimum by changing the element thickness to match the external thermal resistances.

The primary objective of this study is to investigate the potential placement of TE on water tube wall inside a boiler at a coal-fired power plant. The present study looks into the 3D effects of the thermal spreading in the TE module. Parameters such as TE leg cross section area and TE fill factor are examined in order to maximize the electrical power production of the TE without sacrificing the boiler efficiency (i.e., reducing the steam temperature).

Model Description

Figure 1 presents the schematic of the simplified TE module studied. The model only consists of a single TE leg. The TE module consists of outer shells, electrodes, and TE leg. The materials for the outer shells and the electrodes are molybdenum and spattered molybdenum, respectively. The TE module hot side is exposed to the hot gas, the cold side is attached to the boiler tube surface. The dimension of the TE module and its material properties are listed in Table 1. The fill factor (fractional area coverage of TE element relative to cross section are of heat flow) of the TE leg is 40%.

Fig. 1
(a) Schematic of the simplified TE module and (b) the computational geometry used in numerical model
Fig. 1
(a) Schematic of the simplified TE module and (b) the computational geometry used in numerical model
Close modal
Fig. 2
Thermal conductivity of materials
Fig. 2
Thermal conductivity of materials
Close modal
Table 1

Material properties and baseline configuration

ParameterValueUnit
Geometry
TE leg (length × width × thickness)6.3 × 6.3 × 1.1mm
Outer shell (width × length × height)10 × 10 × 0.25mm
Electrode thickness0.1mm
Tube3mm
Materials
Outer shells thermal conductivity (molybdenum)(see Fig. 2)W/m K
Electrodes thermal conductivity (spattered molybdenum)(see Fig. 2)W/m K
TE leg thermal conductivity (β)1.5W/m K
TE leg electric conductivity (σ)250001/Ω m
TE density (ρ)8200kg/m3
TE leg Seebeck coefficient (S)2 × 10−4V/K
Tube thermal conductivity (carbon steel)(see Fig. 2)W/m K
Boundary conditions
Hot side temperature (Thot)1680K
Hot side heat transfer coefficient (hhot)246W/m2 K
Cold side temperature (Tcold)647K
Cold side heat transfer coefficient (hcold)722W/m2 K
ParameterValueUnit
Geometry
TE leg (length × width × thickness)6.3 × 6.3 × 1.1mm
Outer shell (width × length × height)10 × 10 × 0.25mm
Electrode thickness0.1mm
Tube3mm
Materials
Outer shells thermal conductivity (molybdenum)(see Fig. 2)W/m K
Electrodes thermal conductivity (spattered molybdenum)(see Fig. 2)W/m K
TE leg thermal conductivity (β)1.5W/m K
TE leg electric conductivity (σ)250001/Ω m
TE density (ρ)8200kg/m3
TE leg Seebeck coefficient (S)2 × 10−4V/K
Tube thermal conductivity (carbon steel)(see Fig. 2)W/m K
Boundary conditions
Hot side temperature (Thot)1680K
Hot side heat transfer coefficient (hhot)246W/m2 K
Cold side temperature (Tcold)647K
Cold side heat transfer coefficient (hcold)722W/m2 K

Solution Methodology

As shown in Fig. 1(b), the computational geometry of the TE module and boiler tube assembly is discretized into structured-hexagonal cells. The finite element analysis solver ansys is employed to solve for the temperature in each individual cell, thus providing temperature distributions throughout the TE module and boiler tube assembly. The governing equation of the model is
(1)
The performance of the TE leg is then analyzed based on the temperature distributions obtained from the numerical model. The generic power output, wTE, for a given hot side TE leg temperature, Th, and cold side TE leg temperature, Tc, is calculated as
(2)
where
(3)

and ZT, the dimensionless figure-of-merit of the TE material and Z is defined as Z = σS2/β. T is the mean temperature across the TE leg.

The efficiency of the TE is defined as
(4)

where q is the heat flux through the TE leg.

Boundary Conditions.

Figure 3 presents the boundary conditions applied to the assembly. The hot side temperature (Thot) is the temperature of the hot gas inside the boiler chamber; and the heat transfer coefficient is hhot. The cold side temperature (Tcold) is the temperature of the steam inside the steam tube, and the heat transfer coefficient is hcold. The values of the heat transfer coefficients and the temperatures used in this study are listed in Table 1. These values are obtained from results of a computational fluid dynamics simulation of an industrial boiler of a 520 MW coal-fired power plant unit by Lou et al. [12], which is shown in Fig. 4. Please note that the temperature range in this plot has been adjusted to show more even color distributions. The maximum temperature inside the boiler is 2399 K. Since the single TE leg module being studied is only a section of the actual TE module which consists of many TE legs, symmetry boundary condition is imposed on the side surfaces of the TE module and boiler tube. The symmetry boundary conditions sets gradient of the temperature at the wall to be zero (∂T/∂n = 0). The outer surfaces of the electrodes and the TE leg are assumed to be perfectly insulated (q = 0).

Fig. 3
Boundary conditions
Fig. 3
Boundary conditions
Close modal
Fig. 4
Temperature at the center plane of boiler
Fig. 4
Temperature at the center plane of boiler
Close modal

Material properties and baseline configuration are listed in Table 1.

Results and Discussions

Figure 5 presents the temperature distribution on the TE module for the baseline configuration (40% fill factor and 0.1 mm leg thickness). As expected, the temperature at the center of the top shell, where the TE leg is located, is lower than the temperature near the edges. And the opposite occurs on the shell's bottom surface. The average hot side temperature of TE leg is 1437.6 K, and the average cold side temperature is 1016.0 K. Its figure-of-merit, ZT, is 0.81; and, the generic power output per TE leg unit area, wTE,leg, is 4.70 × 104 W/m2. The generic power output per TE module unit area, wTE, is 1.87 × 104 W/m2. The efficiency of the module, ηTE, is 6.9%.

Fig. 5
Temperature distribution on TE module in hot gas at 1640 K
Fig. 5
Temperature distribution on TE module in hot gas at 1640 K
Close modal
Due to the heat flow resistance increase caused by the presence of the TE module, the power output per unit area of the steam cycle, wsteam, reduces from 2.26 × 105 W/m2 to 1.26 × 105 W/m2. The total work output, wtotal(wTE + wsteam), is 1.44 × 105 W/m2. The percentage of gain or loss the work output per unit area after the TE module installation, Δw, is calculated as follows:
(5)

where wo is the work output per unit area of the steam cycle without the installation of the TE module. For this baseline configuration, Δw is −36.2%.

Effects of TE Leg Thickness.

The TE module size was changed to study the effects of the TE leg thickness (d). Only the leg thickness was changed; the fill factor (fractional area coverage of TE element relative to cross section are of heat flow) was not changed, i.e., the length and the width of the TE leg were kept the same. The results are summarized in Table 2.

Table 2

Temperature and power output of TE modules with different leg thickness


Leg thickness, d (mm)
0.30.550.8251.1 (baseline)1.652.2
Th (K)1299.01356.31402.61437.61486.81519.5
Tc (K)1117.01063.41016.0977.8920.7880.1
ZT0.810.810.810.810.800.80
wTE,leg (W/m2)2.70 × 1043.81 × 1044.43 × 1044.70 × 1044.75 × 1044.55 × 104
wTE (W/m2)1.07 × 1041.51 × 1041.76 × 1041.87 × 1041.89 × 1041.80 × 104
ηTE (%)2.94.55.96.98.49.3
wsteam (W/m2)1.82 × 1051.60 × 1051.41 × 1051.26 × 1051.03 × 1058.76 × 104
wo (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 1052.26 × 1052.26 × 105
wtotal (W/m2)1.93 × 1051.75 × 1051.59 × 1051.44 × 1051.22 × 1051.06 × 105
Δw (%)−14.8−22.6−29.9−36.2−46.0−53.3

Leg thickness, d (mm)
0.30.550.8251.1 (baseline)1.652.2
Th (K)1299.01356.31402.61437.61486.81519.5
Tc (K)1117.01063.41016.0977.8920.7880.1
ZT0.810.810.810.810.800.80
wTE,leg (W/m2)2.70 × 1043.81 × 1044.43 × 1044.70 × 1044.75 × 1044.55 × 104
wTE (W/m2)1.07 × 1041.51 × 1041.76 × 1041.87 × 1041.89 × 1041.80 × 104
ηTE (%)2.94.55.96.98.49.3
wsteam (W/m2)1.82 × 1051.60 × 1051.41 × 1051.26 × 1051.03 × 1058.76 × 104
wo (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 1052.26 × 1052.26 × 105
wtotal (W/m2)1.93 × 1051.75 × 1051.59 × 1051.44 × 1051.22 × 1051.06 × 105
Δw (%)−14.8−22.6−29.9−36.2−46.0−53.3

TE efficiency increases as the TE leg thickness increases due to the higher temperature difference between the top and the bottom of the TE leg (Th and Tc). However, increasing the TE leg thickness increases the heat flow resistance and, thus, reduces the steam cycle work output. The steam cycle's work output loss is greater than the work output gained from the TE leg, so the net work output is a loss. The smallest loss is obtained from the TE with leg thickness of 0.3 mm, which has 14.8% total work output loss.

Effects of Adding Fin on TE Module.

In order to increase heat flow through the TE module, a fin is added to its outer shell as shown in Fig. 6. The fin width is 3 mm and the fin height is varied from 2 mm to 10 mm. Temperature distribution on the TE module with fin with the height of 6 mm is shown in Fig. 7. From the results listed in Table 3, it is seen that the increasing the fin length increases the work output of both the TE and the steam cycle as more heat is able to pass through the module due to the presence of the fin. For fin height of 4 mm and higher, there is a total work output gain.

Fig. 6
TE module with fin
Fig. 6
TE module with fin
Close modal
Fig. 7
Temperature distribution on TE module with fin in hot gas at 1640 K
Fig. 7
Temperature distribution on TE module with fin in hot gas at 1640 K
Close modal
Table 3

Temperature and power output of TE module with fin (40% fill factor and 0.3 mm leg thickness)


Fin height (mm)
246810
Th (K)1385.81434.41460.81482.71494.6
Tc (K)1176.31209.71229.31242.61250.9
ZT0.850.880.900.910.92
wTE,leg (W/m2)3.57 × 1044.10 × 1044.35 × 1044.68 × 1044.82 × 104
wTE (W/m2)1.42 × 1041.63 × 1041.73 × 1041.86 × 1041.91 × 104
ηTE (%)3.33.53.63.73.8
wsteam (W/m2)2.09 × 1052.24 × 1052.33 × 1052.39 × 1052.43 × 105
wo (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 1052.26 × 105
wtotal (W/m2)2.10 × 1052.26 × 1052.35 × 1052.41 × 1052.45 × 105
Δw (%)−1.46.210.714.016.0

Fin height (mm)
246810
Th (K)1385.81434.41460.81482.71494.6
Tc (K)1176.31209.71229.31242.61250.9
ZT0.850.880.900.910.92
wTE,leg (W/m2)3.57 × 1044.10 × 1044.35 × 1044.68 × 1044.82 × 104
wTE (W/m2)1.42 × 1041.63 × 1041.73 × 1041.86 × 1041.91 × 104
ηTE (%)3.33.53.63.73.8
wsteam (W/m2)2.09 × 1052.24 × 1052.33 × 1052.39 × 1052.43 × 105
wo (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 1052.26 × 105
wtotal (W/m2)2.10 × 1052.26 × 1052.35 × 1052.41 × 1052.45 × 105
Δw (%)−1.46.210.714.016.0

One of the design criteria is to minimize the effects of the TE on the steam cycle work, i.e., minimize the steam cycle work loss. However, depending on the size of the fin, the heat flow increase due to the fin can result in a heat flux to the steam side to be higher than that of the bare tube (wo), as seen in the case of TE with fin height 6 mm, 8 mm, and 10 mm in Table 3. In order to have a fair evaluation of the potential work output increase by adding the TE module, only TE design which yields steam work output (wsteam) equal or less than that of the bare tube (wo) are considered. In other words, there will be no additional steam work output gained. The steam work output (wsteam) of the TE with fin height of 4 mm is 2.24 × 105 W/m2 very close to the work output of the bare tube (wo = 2.24 × 105 W/m2), and its work output gain is 6.2%. The fin height is further optimized such that wsteam matches wo. The optimized fin height is 4.3 mm and the corresponding percentage of work output gain in 7.1%.

Table 4 also shows the results for different TE fill factors. Similar to the 40% fill factor case, the fin height is optimized such that there is no work output loss for the steam side. The work output gain (which ranges from 6.6% to 7.5%) is slightly higher for TE with higher fill factor. The TE with higher fill factor design requires thicker leg. The higher fill factor and thicker leg means that the TE will require more TE leg material to produce; thus, will cost more to produce. The impact of fill factor is moderate, so the smaller fill factor is preferred in order to reduce the mass of the TE material used significantly. Fill factor of 40% and leg thickness of 0.43 mm are chosen.

Table 4

Temperature and power output of TE module (fin height 4.3 mm, 3 mm height)


Fill factor
30%40%50%60%
Thickness, d (mm)0.210.30.390.47
Th (K)1434.61439.81441.91444.2
Tc (K)1220.41213.41209.41205.8
ZT0.890.880.880.88
wTE,leg (W/m2)5.33 × 1044.17 × 1043.38 × 1042.95 × 104
wTE (W/m2)1.55 × 1041.65 × 1041.70 × 1041.75 × 104
ηTE (%)3.33.53.63.7
wsteam (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 105
wo (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 105
wtotal (W/m2)2.41 × 1052.42 × 1052.43 × 1052.43 × 105
Δw (%)6.67.17.57.5

Fill factor
30%40%50%60%
Thickness, d (mm)0.210.30.390.47
Th (K)1434.61439.81441.91444.2
Tc (K)1220.41213.41209.41205.8
ZT0.890.880.880.88
wTE,leg (W/m2)5.33 × 1044.17 × 1043.38 × 1042.95 × 104
wTE (W/m2)1.55 × 1041.65 × 1041.70 × 1041.75 × 104
ηTE (%)3.33.53.63.7
wsteam (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 105
wo (W/m2)2.26 × 1052.26 × 1052.26 × 1052.26 × 105
wtotal (W/m2)2.41 × 1052.42 × 1052.43 × 1052.43 × 105
Δw (%)6.67.17.57.5

Effects of TE Module Locations.

The hot side temperature of the TE module was varied to represent different locations of TE module installation inside the boiler. Figure 8 shows the average gas temperature near the boiler wall at various locations. The temperatures chosen are 1500 K, 1300 K, and 1150 K, to represent the locations as the gas move downstream from combustion zone. The cold side temperature was not changed.

Fig. 8
Average gas temperature near the wall at various boiler locations
Fig. 8
Average gas temperature near the wall at various boiler locations
Close modal

TE configuration with fill factor of 40% and fin width of 3 mm is used. The TE leg thickness used are 0.3 mm, 0.4 mm, 0.5 mm, and 0.55 mm for the 1680 K, 1500 K, 1300 K, and 1150 K gas temperature, respectively. The TE fin height was optimized for each gas temperature and the results are summarized in Table 5. The percentage of work output gain ranges from 7.1% in the hotter section to 4.4% in the cooler section.

Table 5

Temperature and power output of TE module (40% fill factor, 3 mm fin width)


Gas temperature (K)
1680150013001150
Leg thickness (mm)0.30.40.50.55
Fin height (mm)4.34.43.83.4
Th (K)1439.81284.41104.1977.1
Tc (K)1213.41067.6927.2842.8
ZT0.880.780.730.61
wTE,leg (W/m2)4.17 × 1042.88 × 1041.54 × 1048.08 × 103
wTE (W/m2)1.65 × 1041.14 × 1046.10 × 1033.21 × 103
ηTE (%)3.53.42.82.1
wsteam (W/m2)2.26 × 1051.62 × 1051.07 × 1057.36 × 104
wo (W/m2)2.26 × 1051.62 × 1051.07 × 1057.36 × 104
wtotal (W/m2)2.42 × 1051.74 × 1051.13 × 1057.68 × 104
Δw (%)7.17.05.64.4

Gas temperature (K)
1680150013001150
Leg thickness (mm)0.30.40.50.55
Fin height (mm)4.34.43.83.4
Th (K)1439.81284.41104.1977.1
Tc (K)1213.41067.6927.2842.8
ZT0.880.780.730.61
wTE,leg (W/m2)4.17 × 1042.88 × 1041.54 × 1048.08 × 103
wTE (W/m2)1.65 × 1041.14 × 1046.10 × 1033.21 × 103
ηTE (%)3.53.42.82.1
wsteam (W/m2)2.26 × 1051.62 × 1051.07 × 1057.36 × 104
wo (W/m2)2.26 × 1051.62 × 1051.07 × 1057.36 × 104
wtotal (W/m2)2.42 × 1051.74 × 1051.13 × 1057.68 × 104
Δw (%)7.17.05.64.4

Conclusions

A 3D numerical model was used to investigate the potential placement of TE modules inside a boiler at a coal-fired power plant. Adding TE by itself directly on the tube outer surface increases the heat flow resistance and, thus, reduces the steam cycle work output. A fin can be added to the TE module to increase the heat flow and to avoid the steam cycle work output loss. Different TE configurations (fill factor and leg thickness, and fin height) were examined. The TE module size was optimized for different locations with different flue gas temperatures inside the boiler. Our analysis shows that TE with ZT ∼ 0.61–0.88 at 843–1440 K range can increase power generated in different areas of the boiler by 4.4–7.1%. The overall output power of the plant can increase by ∼6.5%.

Nomenclature

d =

thickness (m)

m =

electrical resistance ratio

S =

Seebeck coefficient (V/K)

w =

power per unit area (W/m2)

β =

thermal conductivity (W/m K)

η =

efficiency

ρ =

density (kg/m3)

σ =

electrical conductivity (1/Ω K)

References

1.
Masahide
,
M.
,
Michio
,
M.
, and
Masaru
,
O.
,
2001
, “
Thermoelectric Generator Utilizing Automobile Engine Exhaust Gas
,”
Therm. Sci. Eng.
,
9
, pp.
17
18
.
2.
Yodovard
,
P.
,
2001
, “
The Potential of Waste Heat Thermoelectric Power Generation From Diesel Cycle and Gas Turbine Cogeneration Plants
,”
Energy Source
,
23
(
3
), pp.
213
224
.
3.
Hsu
,
C. T.
,
Huang
,
G. Y.
, and
Chu
,
H. S.
,
2011
, “
Experiments and Simulations on Low Temperature Waste Heat Harvesting System by Thermoelectric Power Generators
,”
Appl. Energy
,
88
(
4
), pp.
1291
1297
.
4.
Kumar
,
S.
,
Heister
,
S. D.
,
Xu
,
X.
,
Salvador
,
J. R.
, and
Meisner
,
G. P.
,
2013
, “
Thermoelectric Generators for Automotive Waste Heat Recovery Systems. Part I: Numerical Modeling and Baseline Model Analysis
,”
J. Electron. Mater.
,
42
(
4
), pp.
665
674
.
5.
Amatya
,
R.
, and
Ram
,
R. J.
,
2010
, “
Solar Thermoelectric Generator for Micropower Applications
,”
J. Electron. Mater.
,
39
(
9
), pp.
1735
1740
.
6.
Jarman
,
J. T.
,
Khalil
,
E. E.
, and
Khalaf
,
E.
,
2013
, “
Energy Analyses of Thermoelectric Renewable Energy Sources
,”
Open J. Energy Effic.
,
2
(
4
), pp.
143
145
.
7.
Eswaramoorthy
,
M.
,
Shanmugam
,
S.
, and
Veerapan
,
A. R.
,
2013
, “
Experimental Study on Solar Parabolic Dish Thermoelectric Generator
,”
Int. J. Energy Eng.
,
3
(3), pp.
62
66
.
8.
Maslamani
,
T. M.
,
Omer
,
A. I.
, and
Majid
,
M. A.
,
2014
, “
Development of Solar Thermoelectric Generator
,”
Eur. Sci. J.
,
10
(
9
), pp.
123
134
.
9.
Chen
,
M.
,
Lund
,
H.
, and
Rosendahl
,
L. A.
,
2010
, “
Energy Efficiency Analysis and Impact Evaluation of the Application of Thermoelectric Power Cycle to Today's CHP Systems
,”
Appl. Energy
,
87
(
4
), pp.
1231
1238
.
10.
Yazawa
,
K.
,
Koh
,
Y. R.
, and
Shakouri
,
A.
,
2013
, “
Optimization of Thermal Topping Combine Steam Turbine Cycles for Energy Economy
,”
Appl. Energy
,
109
, pp.
1
9
.
11.
Yazawa
,
K.
, and
Shakouri
,
A.
,
2012
, “
Optimization of Power and Efficiency of Thermoelectric Devices With Asymmetric Thermal Contact
,”
J. Appl. Phys.
,
111
(
2
), p.
024509
.
12.
Lou
,
X.
,
Zhang
,
Q.
,
Wu
,
B.
,
Heim
,
J.
, and
Zhou
,
C. Q.
,
2012
, “
Development of a Virtual Power Plant Boiler for Training
,”
ASME
Paper No. HT2012-58426, pp.
1085
1091
.