Modern gas turbines operate at turbine inlet temperatures well above the melting point of their combustor’s materials. This is possible due to very sophisticated cooling systems and ceramic shielding of the parts that experience the highest thermal loads. Cooling is realized by extracting part of the air leaving the compressor and directing it to the hot gas path components, while the remaining air enters the combustion chamber. The latter air stream is in turn split in two, one taking part directly in the combustion reaction and the other used to cool the walls of the chamber. To compute the power generated from thermoelectric generators installed in a gas turbine combustion chamber, the temperature distribution in the walls of the chamber must be first calculated. This can be achieved only if both the chamber geometry and the air and fuel mass flow rates are known.
The multi-level modeling approach adopted in the current work is depicted in
Figure 1. On a first level, a full gas turbine thermodynamic model, with the respective secondary air system, is built and validated. This model computes the flow rate, pressure and temperature of the working fluid in all major stations of the gas turbine process. This way, the air flow rate used for cooling and sealing in the turbine expander is estimated along with its temperature and pressure. The same model delivers values of the air and fuel flow rates entering the combustion chamber, along with their inlet temperature and the outlet temperature of the exhaust gases. On a second level, a separate model of the combustion chamber is developed to compute the mass flow rates of the air streams in its interior and the temperature distribution in its walls. The wall temperature distribution is first computed in a conventional combustion chamber with ceramic heat shields as a benchmark case. Subsequently, the heat shields are replaced by high temperature thermoelectric generators (TEGs) and the wall temperature distribution is computed again. This result is fed in a final model that computes the power generated from the TEGs for the available temperature difference across them. For each TEG configuration (length, fill factor, etc.), the detailed combustion chamber model is used anew to compute the new wall temperature distribution. The modeling approach is depicted graphically in
Figure 1, whereas the details of the models are given in the following sections.
2.2. Combustor Model
The aim of the current work is to study the power generation potential of thermoelectric generators when they are integrated in the liner wall of a gas turbine combustor. The power generated from TEG depends on the temperature difference across them. This temperature distribution will result from a study of the heat transfer conditions in the combustion chamber.
To simplify the analysis and make the findings easier to evaluate, a tubo-annular combustion chamber is modeled as a co-flow heat exchanger with a complex wall structure (
Figure 5). The inlet conditions of the combustion chamber are taken from the gas turbine model, described in the previous section. Combustion takes place in the inner side of the chamber and cooling air flows in the outer annulus.
The temperature distribution in the combustor is computed based on the recommendations of Lefebvre and Ballal [
21]. The applied model computes the convective and film cooling heat transfer coefficients along with radiative heat transfer from the hot gas side and assumes one-dimensional heat transfer through the inner combustor wall. The outer wall of the chamber is taken as perfectly insulated.
The studied chamber has seven film cooling slots at a distance of 5 cm from each other and a cooling hole at the end of each film cooling region. The distance of the slots and holes is chosen so that the film of air is not perturbed. Ideal mixing is assumed between the hot combustion gases and the film cooling air. Furthermore, it is assumed that the hot gases’ compositions remain the same along the length of the combustion chamber. The cooling air that does not flow through the film cooling slots and holes is mixed with the main stream at the end of the chamber (not shown in
Figure 5). A sketch of the chamber geometry along with the geometrical assumptions [
20,
22] made for the film cooling are presented in
Figure 5 and in
Table 5.
The depth of film cooling
s determines the liner wall diameter variation. After every film cooling slot, the liner wall inner diameter increases as shown in the following equation:
In contrast, the outer casing diameter remains constant along the axial length, resulting in a reduction of the annulus flow area after every film cooling slot.
2.3. Thermoelectric Generator Model
As part of the combustor liner wall, the thermoelectric generator replaces the ceramic heat shield layer, as depicted in
Figure 6. It is placed between the hot combustion products, which become its high temperature heat source, and the metal liner, which is its heat sink. Two thin ceramic insulation layers, on which the semiconductor materials are fitted, are still present to protect the TEG from the aggressive hot combustion gases and insulate its cold junction from the liner metal wall.
The introduction of TEG in the combustion chamber does not affect its basic geometry. The diameter of the metal liner wall remains the same as in the case with simple ceramic heat shield. The number of junctions in the TEG and its length “
L” (see
Figure 6) depend on its materials and the techno-economic optimization of the heat transfer system.
For the application of energy conservation at the hot and cold junctions, four energy transport mechanisms are taken into account, which result in Equations (
2) and (
3) for the hot and the cold junction, respectively [
23,
24].
In these equations,
and
are heat input from the heat source to the hot junction and the heat output from the cold junction to the heat sink, respectively.
represents internal resistive heating, while
is the electrical energy output from thermoelectric effect
. Finally,
corresponds to the heat conducted from the hot to the cold junction. It may be assumed that the Joule heating is equal for both junctions, since the temperature profiles for the studied materials are nearly linear [
23,
25]. The generated electrical power and the efficiency of the TEG can be thus given from Equations (
4) and (
5).
For the computation of the temperature distribution in the TEG and the liner, the same convective and radiative heat transfer models are applied as in the reference case. The heat balance statement for the two TEG junctions can be developed to a closed system of five independent equations with the five wall temperatures as unknowns (see
Figure 6). The system is solved under the assumption that the total wall thickness is small compared to its length (i.e., the same heat exchange surface is assumed both at the cold and the hot end).
2.4. Thermoelectric Module Optimization Model
Up to date, TEGs have demonstrated efficiencies (see Equation (
5)) in the region of 2–3% [
4]. It is thus expected that the generated power from the proposed system will be relatively low, compared to that generated by the host gas turbine. Hence, a simple power output optimization of the TEG would not make much sense if it was not supported by a techno-economic optimization. In other words, materials and geometrical characteristics have to be optimized on a techno-economic basis. On the other hand, the economic evaluation of high temperature TEGs is challenging, because they have not reached commercial maturity yet. Under this premise, the TEG optimization in this work relies on the work of LeBlanc et al. [
3] and the material prices reported in it. For the same reason, the following analysis is more an attempt to identify the most important parameters that affect the financial viability of TEGs, rather than providing absolute numbers for their economic performance.
The figure of merit of a thermoelectric material relates its properties to the device efficiency but does not take into account costs. Much more appropriate for the intended optimization is the use of the specific power generation costs (Equation (
11)).
The installation and operation of TEGs in the configuration in question will incur only investment costs, while no additional operating and maintenance costs will be present. The investment costs include the volumetric (
) and areal (
) TEG module costs and the heat exchanger costs for the energy conversion system [
3], and can be expressed by Equation (
12).
where
F represents the fill factor of the TEG module,
is the overall heat transfer coefficient of the heat exchanger and
is the specific costs of heat exchanger (
$/Wk). The fill factor is defined as the ratio of the area covered by the active thermoelectric material to the plate area of the TEG module (cf. [
3]). In the current work, the heat exchanger is the combustion chamber and the respective costs have been neglected. In addition to the thermoelectric system cost, balance of plant and installation costs are estimated each as a 10% of the total system costs [
26]. The resulting total cost of the TEG is thus 20% higher than that resulting from Equation (
12).The optimization objective is to minimize
G while the TEG length, its fill factor
F and its dimensionless figure m are being varied. The latter is defined as the ratio between the electrical load resistance R
and the internal resistance of the module R. According to LeBlanc et al. [
3] (Equation (
13)), the value of m that minimizes
G is the load matching condition
m = 1. Hence, the two final design parameters to be optimized are
L and
F.
In general, as the TEG leg length or its fill factor decrease, less material is used and the costs are reduced. However, reducing the leg length also reduces the temperature drop between the TEG junctions, resulting in less power output, which increases
G. At the same time, manufacturing feasibility must be taken into account as well, since the thermoelectric device cannot be arbitrarily small. As stated in [
3], there is not a minimum value for
G, but there is a point of diminishing returns at finite
L and
F. The task of the presented optimization is to find this point. This is done by fixing the TEG fill factor
F at values between 0 and 1 and computing the respective optimal leg length based on Equation (
13) from [
3].
For each combination of the two varying parameters, L and F, the power output, the cost and G are determined and compared. From the thermoelectric materials presented in the work of LeBlanc et al. [
3], only nanostructured bulk silicon [
27] and Al-doped ZnO [
28] qualify for the expected wall temperatures in a gas turbine combustion chamber, which are presented in
Table 6. The costs for each material are computed according to LeBlanc et al. [
3].
Once these parameters, the design of the TEG, its costs and the power it generates are defined, and the net present value and the internal rate of return for the respective investment are computed. The former is computed with Equation
14, where
N is the lifetime of the investment in years,
CF is the net cash-flow of each year and “
i” is the discount rate used to project each cash-flow to the starting year of the investment. The internal rate of return (IRR) is the discount rate value that results to a net present value equal to zero.
The annual cash flows are the difference between the generated income and the incurred costs. In the system in question, income is generated by selling electricity, whereas only the initial installation cost is accounted for (see Equation (
12)).