Modeling of an Integrated Thermoelectric Generation–Cooling System for Thermoelectric Cooler Waste Heat Recovery

Abstract

This paper focuses on the problem of thermoelectric cooler waste heat recovery and utilization, and proposes taking the waste heat together with the original heat source as the input heat source of the integrated thermoelectric generation–cooling system. By establishing an analytic model of this integrated thermoelectric generation–cooling system, the steady-state and transient thermal effects of this system are analyzed. The steady-state analysis results show that the thermoelectric generator’s actual heat source is about 20% larger than the intrinsic heat source. The transient analysis results prove that the current of thermoelectric power generation and the cold end temperature of the system show a nonlinear change rate with time. The cold end temperature of the system has a maximum value. Under different intrinsic heat sources, this maximum value can be reached between 1 s and 2.5 s.

1. Introduction

In recent years, thermoelectric modules (TEMs) have attracted much attention because of their stability, direct energy conversion, and having no moving parts in terms of energy management, temperature control, etc. [1,2,3]. Thermoelectric devices are divided into thermoelectric generators (TEGs), thermoelectric coolers (TECs), etc. TEGs use the Seebeck effect to convert thermal energy into electric energy, and use the Peltier effect to realize the function of a heat pump [4,5,6].

To date, many scholars have done theoretical and applied research on thermoelectric modules. Lam et al. [7] obtained the analytical solutions of transient temperature distributions within semiconductor elements by solving the governing heat conduction equation using the superposition technique with the aid of a special transformation function. The performance of a thermoelectric cooler (TEC) under the influence of the Thomson effect, Joule heating, and combined radiation and convection cooling is investigated in this study. In addition, Rahman et al. [8] presented a comprehensive analysis of a dual, micro-gap thermionic–thermoelectric hybrid energy converter by using a self-consistent iterative algorithm considering the energy balance condition. The maximum conversion efficiency of the TEG reached 9.52% by optimizing the system. Shen et al. [9] presented a three-dimensional finite element analysis that couples the high-temperature proton exchange membrane fuel cell (HT-PEMFC) and thermoelectric generator into a single model. Through the analysis and optimization of this model, it was found that a TEG with only two pairs of thermocouples can improve its conversion efficiency by 1.2%. Gong et al. [10] optimized the TEC’s input current, leg geometry, and contact layer by taking into account the material’s temperature characteristics, the Thomson effect, and contact layer resistance to achieve maximum cooling performance and high operational reliability. Cuce et al. [11] investigated the impacts of nanofluid use in TECs on cooling power and the main performance parameters through a comprehensive experimental methodology. Their results indicated that the best enhancement in indoor temperature difference compared to water cooling is obtained for 1% Al₂O₃ at 30 °C ambient temperature with 26% for without load conditions.

In general, the output power of a TEG is enough to drive a TEC in refrigeration, and more and more research is being conducted on thermoelectric generation–cooling systems [12,13,14]. Kwan et al. [15] significantly improved the TEG–TEC control technique by including components that ensure the TEG’s output power can truly be returned to the original power source. In addition, steady-state analysis of the TEG–TEC model with input waste heat parameter scanning shows that the TEG mode can only function in a certain range. Lin et al. [16] proposed a new design of combined TEC–TEG systems, where two single-stage TEGs are employed to separately power the hot stage and cold stage of the TEC. The analysis results of the system show that the cooling capacity of the new design is enhanced by 75% and the maximum temperature drop is elevated by 76.8% compared with the original design.

Analysis of the thermoelectric modules has found that the TEC’s cooling process will release much more heat than its cooling load. Scholars have put forward some solutions for the recovery and utilization of this part of the heat. Teffah et al. [17] proposed a new TEG–TEC system, in which the TEG plays the role of a partial heatsink for the TEC by transferring this waste heat to the total system heatsink and converting an amount of this heat into electricity by a phenomenon called the Seebeck effect, for this thermoelectric module. The analysis results of this system show that the TEG can effectively convert this portion of heat into electrical energy. However, in this case, the output power of the TEG is very small, and since the TEG and TEC have the same hot-end temperature, the overall temperature of the TEC increases and the cooling capacity decreases.

In view of this situation, this paper proposes to collect the thermoelectric cooler’s waste heat and to use it together with the inherent heat source of the integrated thermoelectric generation–cooling system as its input heat source for the purpose of energy recycling. By establishing an analytic model of this integrated thermoelectric generation–cooling system, the steady-state and transient thermal effects of the system are analyzed. As the heat released by the thermoelectric cooler presents a nonlinear change, the heat absorption of the thermoelectric generator will present a nonlinear change, which will have a great impact on the transient performance of the system. Therefore, it is necessary to conduct transient response analysis on the integrated thermoelectric generation–cooling system.

2. Theory and Simulation

2.1. Theoretical Model

A schematic diagram of the integrated thermoelectric generation–cooling system is shown in Figure 1. The output electric energy of the TEG is taken as the input electric energy of the TEC. At the same time, the heat loss of the TEC is collected and used together with the intrinsic heat source Q_s as the heat source of the TEG. In this way, not only can the TEC waste heat be recovered and utilized, but also the output power of TEG is improved. The hot end of the TEG, Q_s, and the “heat accumulation” are direct thermal connections.

The governing equation of the integrated thermoelectric generation–cooling system is [18,19]:

$$\rho C_P\frac{\partial T}{\partial t}=\nabla (\kappa \nabla T)+\frac{{\overrightarrow{J}}^2}{\sigma }-\beta \overrightarrow{J}\cdot \nabla T$$

(1)

where ρ, C_p, κ, $\overrightarrow{J}$, σ, and β are the density, heat capacity, thermal conductivity, current density, electrical conductivity, and Thomson coefficient, respectively. It should be noted that the thermal conductivity, electrical conductivity, and Thomson coefficient of thermoelectric materials are functions of the temperature.

The energy transfer equation in the thermoelectric leg is [20]:

$$Q=\alpha IT\pm \kappa A\cdot \nabla T$$

(2)

where α and A are, respectively, the Seebeck coefficient and the cross-sectional area. The energy transmission direction is shown in Figure 1.

According to Equation (2), the heat absorption Q_l,C of the TEC and the heat absorption expressions of the TEG can be obtained as:

$$Q_{l,C}={\alpha }_{l,C}I{T}_{l,C}-{\kappa }_{l,C}\nabla T_{l,C}$$

(3)

$$Q_{h,G}={\alpha }_{h,G}I{T}_{h,G}+{\kappa }_{h,G}\nabla T_{h,G}.$$

(4)

In addition, the heat source of the TEG is the sum of the heat loss of the TEC and the intrinsic heat source Q_s; so, the total heat source of the TEG can be obtained as:

$$Q_{h,G}=Q_s+Q_{h,C}=Q_s+{\alpha }_{h,C}I{T}_{h,C}-{\kappa }_{h,C}A\cdot \nabla T_{h,C}.$$

(5)

Since both the TEG and TEC have the same voltage and internal current, the internal current of the integrated thermoelectric generation–cooling system is:

$$I=\frac{({\alpha }_{h,G}{T}_{h,G}-{\alpha }_{l,G}{T}_{l,G})-({\alpha }_{h,C}{T}_{h,C}-{\alpha }_{l,C}{T}_{l,C})}{R_G+R_C}.$$

(6)

It can be seen from Equations (3)–(6) that the TEG and TEC have a strong coupling relationship in terms of heat and current. The temperature distribution of the TEC is determined by the current and heat absorption of the TEC. In addition, the temperature distribution of the TEC can change the current in the system by changing the heat absorption of the TEG. The complex coupling relationship between thermoelectric power generation and thermoelectric refrigeration results in the complex nonlinear variation trend of the heat absorption Q_h,G, current I, and cold end temperature T_l,C of the system.

If the time term in Formula (1) is 0, the steady state heat absorption can be obtained as follows:

$$Q_{h,G}={\alpha }_{h,G}I{T}_{h,G}-\frac{1}{2}{I}^2{R}_G+K_G(T_{h,G}-T_{l,G})$$

(7)

$$Q_{l,C}={\alpha }_{l,C}I{T}_{l,C}-\frac{1}{2}{I}^2{R}_C-K_C(T_{h,C}-T_{l,C})$$

(8)

where $R_G$ and $R_C$ are the resistance of the TEG and TEC, respectively; and $K_G$ and $K_C$ are the thermal conductivity of the TEG and TEC, respectively.

The relationship between the TEC’s hot-end temperature T_h,C, the TEG’s cold-end temperature T_l,G, the ambient temperature T_∞, and the equivalent convective heat transfer coefficient h_eff is:

$$\begin{array}{l}{Q}_{h,C}=h_{eff}(T_{h,C}-T_{\infty })\\ Q_{l,G}=h_{eff}(T_{l,G}-T_{\infty })\end{array}.$$

(9)

2.2. Numerical Simulation Model

In the numerical simulation, the logarithm of the p–n junction for the TEG and TEC is 18 pairs, the thermoelectric leg size is 1 × 1 × 1 mm³, and the leg spacing is 1 mm.

In the numerical analysis of the integrated thermoelectric generation–cooling system, the following assumptions are made:

(1)

one can assume that all of the other surfaces in the model have thermal insulation except for the cold and hot ends [21];

(2)

one can ignore the contact resistance and thermal resistance at the interface;

(3)

one can assume that there is no heat loss in the heat transfer process.

The boundary conditions are as follows:

$$\begin{array}{l}{Q}_s=9.68\mathrm{W}\\ Q_{l,C}=2.42\mathrm{W}\end{array}.$$

(10)

The hot end of the TEC and the cold end of the TEG are in the same environment [22]:

$$\begin{array}{l}{T}_{\infty }=300\mathrm{K}\\ h_{eff}=7712{\mathrm{W}/(\mathrm{m}}^2\cdot \mathrm{K})\end{array}.$$

(11)

In addition, the thermoelectric material selected for the thermoelectric leg is Bi₂Te₃, and its performance parameters are shown in

Table 1. Material properties of Bi₂Te₃ [23].

	Material Properties
P-type	$\kappa =1.472\times [1-1.29(T-300)\times {10}^{-3}+1.35{(T-300)}^2\times {10}^{-5}]\mathrm{W}/(\mathrm{m}\cdot \mathrm{K})$
	$\alpha =2.207\times [1+1.55(T-300)\times {10}^{-3}-3.15{(T-300)}^2\times {10}^{-6}]\times {10}^{-4}\mathrm{V}/\mathrm{K}$
	$1/\sigma =8.826\times [1+5.88(T-300)\times {10}^{-3}+8.93{(T-300)}^2\times {10}^{-6}]\times {10}^{-6}\mathrm{\Omega }\cdot \mathrm{m}$
N-type	$\kappa =1.643\times [1-9.8(T-300)\times {10}^{-4}+1.56{(T-300)}^2\times {10}^{-5}]\mathrm{W}/(\mathrm{m}\cdot \mathrm{K})$
	$\alpha =-2.23\times [1+5.62(T-300)\times {10}^{-4}-4.56{(T-300)}^2\times {10}^{-6}]\times {10}^{-4}\mathrm{V}/\mathrm{K}$
	$1/\sigma =8.239\times [1+4.7(T-300)\times {10}^{-3}+2.67{(T-300)}^2\times {10}^{-6}]\times {10}^{-6}\mathrm{\Omega }\cdot \mathrm{m}$

[23].

The electrode material is copper, its thermal conductivity is $k=400\mathrm{W}/(\mathrm{m}\cdot \mathrm{K})$, the conductivity is $\sigma =5.88\times {10}^8\mathrm{S}/\mathrm{m}$, and the thermal conductivity of the ${\mathrm{Al}}_2{\mathrm{O}}_3$ ceramic plates is $k=35.3\mathrm{W}/(\mathrm{m}\cdot \mathrm{K})$.

All numerical simulations were performed in COMSOL Multiphysics 5.5.

Figure 2 shows the difference in temperature variation between the numerical predictions here and the experimental results reported in a previous study [24]. It can be seen from Figure 2 that the numerical results obtained by the proposed numerical model closely agree with the experimental data. The maximum temperature difference of T_h and T_c between the numerical predictions and the measurement data is less than 2 K at t = 75 s and 1.5 K at t = 50 s, respectively. Hence, the validity of the proposed numerical model can be ensured.

3. Results and Analysis

3.1. Steady State Result Analysis

Under the above-described dimensions and boundary conditions, the cold-end temperature T_l,C is 294.58 K and the current I is 1.4 A, indicating that the integrated thermoelectric generation–cooling system can realize the refrigeration function. Based on the proposed model, a change in the intrinsic heat source Q_s will change T_l,C, I, and Q_h,G, as shown in Figure 3.

As can be seen from Figure 3, T_l,C decreased with the increase of Q_s, from the maximum value of 303 K to 282.97 K. In addition, the current I increased with the increase of Q_s, from 1.1 A to 1.87 A. Q_h,G increased with the increase of Q_s, from 8.8 W to 22 W.

In this system, the actual heat absorption Q_h,G of the TEG’s hot end is the sum of Q_s and Q_h,C, as shown in formula (5). The increase in Q_h,G is a matter of course. By comparing the difference between Q_h,G and Q_s in Figure 3, it can be found that the difference also increases with the increase of Q_s, from 2.8 W to 3.8 W. This difference is greater than the TEC’s cooling load Q_l,C, which is 2.42 W. This suggests that the TEC produces a lot of waste heat, which can be collected for the TEG’s power generation to boost its output.

As Q_h,G increases, it can be seen from formula (4) that the current I in the system also increases. It can be seen from formula (8) that Q_l,C is a quadratic function of current I, and there is an optimal current to get the minimum value of T_l,C. When current I is less than the optimal value, T_l,C decreases with the increase of I [25]. In Figure 3, the maximum value of current I is 1.87 A, which is far from the optimal current value. Therefore, T_l,C decreases with the increase of I, and T_l,C decreases with the increase of Q_s.

3.2. Transient Analysis

It can be seen from the steady-state analysis results that the integrated thermoelectric generation–cooling system has a good refrigeration capacity. Since the TEC’s dissipated heat is also the TEG’s heat source, the actual heat source of the TEG (Q_h,G) is about 20% larger than the intrinsic heat source (Q_s).

The transient performance of the system was analyzed to see whether extremely high or extremely low temperatures would occur during the process of the system reaching a stable state. Transient analysis was conducted on the integrated thermoelectric generation–cooling system, and the variation trend of T_l,C, I, and T_h,G over time was obtained as shown in Figure 4.

It can be seen that in Figure 4, T_l,C, I, and T_h,G all show a nonlinear change trend over time under the condition that the intrinsic heat source Q_s is 9.68 W. Among them, the TEC’s cold-end temperature Q_l,C has a maximum value of 309.02 K, which is obtained at 1.7 s, 14.54 K higher than the steady-state value of 294.48 K. The current I of the integrated thermoelectric generation–cooling system and the hot-end temperature Q_h,G of the TEG showed a trend of gradual increase, reaching a stable state at 8 s and 14 s, respectively.

It should be noted that the boundary condition of the TEG in this system is the heat flow boundary condition; the temperature T_h,G at the hot end of the TEG increases gradually under the action of heat Q_h,G, as shown in Figure 4. It can be seen from formula (6) that the current I of the system is determined by the temperature difference between the TEG and TEC. At the beginning of the system’s operation, the temperature difference between hot and cold ends of TEC is very small, and the current I is mainly determined by the temperature difference between hot and cold ends of TEG. The temperature difference of the TEG increases with the increase of T_h,G, so the current I and T_h,G have the same change trends, both of which are gradually increasing. In the TEC, the Peltier effect is the main effect of refrigeration. With the Seebeck coefficient of thermoelectric materials unchanged, the Peltier heat increases with the increase of current I. The first term on the right side of formula (3) is Peltier heat. When the Peltier heat is less than Q_l,C, only a negative value of the second term to the right of Equation (3) can make the equation true. In this case, the temperature at the cold end of the TEG is higher than that at the hot end. In Figure 4, the increase of T_l,C before 1.7 s corresponds to the situation where the current is small and the Peltier heat is less than Q_l,C. When the current increases to the size of the Peltier heat greater than Q_l,C, the cold-end temperature T_l,C of the TEC begins to drop until it reaches a stable state.

By observing the time when T_l,C, I, and T_h,G reach the stable state in Figure 4, it can be found that current I reaches the stable state earlier than T_l,C and T_h,G. According to Formula (6), the current I is jointly determined by the temperature of the TEG and TEC. The current I reaches a stable state at about 8 s. At this time, the equilibrium state is the dynamic equilibrium state caused by the synchronous increase of the temperature difference between hot and cold ends of TEG and TEC. That is, the increase of the TEG open-circuit voltage caused by the increase of the TEG temperature difference is the same as the increase of the TEC reverse potential difference caused by the increase of the TEC temperature difference. Thus, the molecular value in formula (6) remains unchanged and the current I reaches a stable value. It should be noted that even after current I reaches a stable value, the temperature change of the TEC caused by current I will continue for a certain time. This is because the thermal change in Joules caused by the current change cannot be completed immediately on the cold and hot end face of the TEC. As a result, I reaches a stable state earlier than T_l,C. Combined with formula (2) and formula (5), it can be seen that in the process of TEC temperature change, the heat dissipated by its hot end Q_h,C changes along with it, which leads to the changes of Q_h,G and T_h,G. That is, I reaches a stable state earlier than T_h,G.

3.3. Transient Analysis of Different Q_s

The results of the steady-state analysis show that the size of Q_s has a great influence on the cold-end temperature T_l,C and current I of the integrated thermoelectric generation–cooling system. In the transient analysis and theoretical analysis of the system, current I is the direct factor affecting the T_l,C. Therefore, the variation of T_l,C and I was analyzed under different Q_s.

Transient analysis was carried out of the integrated thermoelectric generation–cooling system with Q_s of 6.05 W, 9.68 W, and 14.52 W, and the variation trend of T_l,C and I over time under different intrinsic heat sources Q_s was obtained, as shown in Figure 5.

It can be seen that in Figure 5, under a different intrinsic heat source Q_s, the variation trend of each parameter is the same. The current I tends to gradually increase to a stable state, which is caused by the gradual increase of the TEG temperature difference. The cold-end temperature T_l,C of the TEC increased first and then decreased before reaching a stable state. This is caused by the fact that the Peltier heat under a small current is smaller than the cooling load Q_l,C, while the Peltier heat under a large current is larger than the cooling load Q_l,C.

Although the change of Q_s does not change the trend of T_l,C and I, it does change the maximum value and corresponding time of T_l,C. As shown in Figure 5, when Q_s was 6.05 W, 9.68 W, and 14.52 W, the maximum value of T_l,C was 312.48 K, 309.02 K, and 306.51 K, respectively. In addition, the time corresponding to the maximum of T_l,C was 2.6 s, 1.7 s, and 1.1 s, respectively. In other words, with the increase of Q_s, T_l,C can reach its maximum value in a shorter time. In order to further analyze this result, transient results under different Q_s conditions were processed. The curve of the maximum value of the T_l,C and its corresponding time change with Q_s were obtained as shown in Figure 6.

In Figure 6, T_lmax represents the maximum value of T_l,C and t_m represents its corresponding time.

As can be seen from Figure 6, the maximum value of the cold-end temperature T_lmax and the corresponding time t_m decrease with the increase of Q_s.

According to formula (3), T_l,C is determined by the magnitude of Peltier heat and Fourier heat. In the case of constant Q_l,C, the current I directly determines the change of the TEC’s temperature. Figure 7 shows the variation curves of T_l,C and I within 0–4 s under different Q_s.

It can be seen from Figure 7 that, although the maximum value and corresponding time of T_l,C are different under different Q_s, the corresponding current is between 0.79 A and 0.94 A. Combined with the previous analysis of the T_l,C variation trend, the current I directly determines the change of T_l,C. Therefore, the increasing rate of current I determines the magnitude of t_m. It can be seen from formula (6) that the current I and the temperature difference between TEGs have the same change trend. Obviously, the rate of increase of the temperature difference between TEGs increases as Q_s increases. Correspondingly, I with the same variation trend as the temperature difference between TEGs can reach a specific value in a shorter time, which makes T_l,C reach its maximum value.

4. Conclusions

In this paper, the waste heat of a TEC is collected as part of the input heat source of an integrated thermoelectric generation–cooling system for the purpose of waste heat recovery and utilization. In addition, steady-state and transient analyses are carried out for the system. The content of this paper can be summarized as follows:

(1)

A theoretical model of the integrated thermoelectric generation–cooling system is established. In this model, the heat source of a TEG presents a nonlinear changing trend, resulting in the nonlinear changing trend of the current and cold-end temperature.

(2)

Under the steady-state condition, the actual heat source of the thermoelectric generator is about 20% larger than the intrinsic heat source. At the maximum value of the intrinsic heat source, the cold-end temperature of the thermoelectric cooler reaches the minimum value of 282.97 K.

(3)

The transient analysis of integrated thermoelectric generation–cooling system shows that the cold-end temperature in the system has a maximum value. Under different intrinsic heat sources, this maximum value can be reached between 1 s and 2.5 s.