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

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.


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 Energies 2020, 13, 4691 2 of 10 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 2 O 3 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.

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 energy transfer equation in the thermoelectric leg is [20]: 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 Ql,C of the TEC and the heat absorption expressions of the TEG can be obtained as: In addition, the heat source of the TEG is the sum of the heat loss of the TEC and the intrinsic heat source Qs; so, the total heat source of the TEG can be obtained as: Since both the TEG and TEC have the same voltage and internal current, the internal current of the integrated thermoelectric generation-cooling system is: 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 Qh,G, current I, and cold end temperature Tl,C of the system.
If the time term in Formula (1) is 0, the steady state heat absorption can be obtained as follows: The governing equation of the integrated thermoelectric generation-cooling system is [18,19]: where ρ, C p , κ, → 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]: 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 = α l,C IT l,C − κ l,C ∇T l,C 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: Since both the TEG and TEC have the same voltage and internal current, the internal current of the integrated thermoelectric generation-cooling system is: 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 Energies 2020, 13, 4691 4 of 10 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 l,C = α l,C IT l,C − 1 2 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:

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 3 , 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: The hot end of the TEC and the cold end of the TEG are in the same environment [22]: In addition, the thermoelectric material selected for the thermoelectric leg is Bi 2 Te 3 , and its performance parameters are shown in Table 1 [23].  [23].

Material Properties
The electrode material is copper, its thermal conductivity is k = 400 W/(m · K), the conductivity is σ = 5.88 × 10 8 S/m, and the thermal conductivity of the Al 2 O 3 ceramic plates is k = 35.3 W/(m · K).
Energies 2020, 13, 4691 5 of 10 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. Al O ceramic plates is 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 Th and Tc 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.

Steady State Result Analysis
Under the above-described dimensions and boundary conditions, the cold-end temperature Tl,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 Qs will change Tl,C, I, and Qh,G, as shown in Figure 3. As can be seen from Figure 3, Tl,C decreased with the increase of Qs, from the maximum value of 303 K to 282.97 K. In addition, the current I increased with the increase of Qs, from 1.1 A to 1.87 A. Qh,G increased with the increase of Qs, from 8.8 W to 22 W.
In this system, the actual heat absorption Qh,G of the TEG's hot end is the sum of Qs and Qh,C, as shown in formula (5). The increase in Qh,G is a matter of course. By comparing the difference between Qh,G and Qs in Figure 3, it can be found that the difference also increases with the increase of Qs, from 2.8 W to 3.8 W. This difference is greater than the TEC's cooling load Ql,C, which is 2.42 W. This

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. Al O ceramic plates is 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 Th and Tc 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.

Steady State Result Analysis
Under the above-described dimensions and boundary conditions, the cold-end temperature Tl,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 Qs will change Tl,C, I, and Qh,G, as shown in Figure 3. As can be seen from Figure 3, Tl,C decreased with the increase of Qs, from the maximum value of 303 K to 282.97 K. In addition, the current I increased with the increase of Qs, from 1.1 A to 1.87 A. Qh,G increased with the increase of Qs, from 8.8 W to 22 W.
In this system, the actual heat absorption Qh,G of the TEG's hot end is the sum of Qs and Qh,C, as shown in formula (5). The increase in Qh,G is a matter of course. By comparing the difference between Qh,G and Qs in Figure 3, it can be found that the difference also increases with the increase of Qs, from 2.8 W to 3.8 W. This difference is greater than the TEC's cooling load Ql,C, which is 2.42 W. This 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 Energies 2020, 13, 4691 6 of 10 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 .

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.
Energies 2020, 13, x FOR PEER REVIEW 6 of 10 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 Qh,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 Ql,C is a quadratic function of current I, and there is an optimal current to get the minimum value of Tl,C. When current I is less than the optimal value, Tl,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, Tl,C decreases with the increase of I, and Tl,C decreases with the increase of Qs.

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 (Qh,G) is about 20% larger than the intrinsic heat source (Qs).
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 Tl,C, I, and Th,G over time was obtained as shown in Figure 4.  It can be seen that in Figure 4, Tl,C, I, and Th,G all show a nonlinear change trend over time under the condition that the intrinsic heat source Qs is 9.68 W. Among them, the TEC's cold-end temperature Ql,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 Qh,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 Th,G at the hot end of the TEG increases gradually under the action of heat Qh,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 Th,G, so the current I and Th,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 Ql,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 Tl,C before 1.7 s corresponds to the situation where the current is small and the Peltier heat is less than Ql,C. When the current increases to the size of the Peltier heat greater than Ql,C, the cold-end temperature Tl,C of the TEC begins to drop until it reaches a stable state. 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 Energies 2020, 13, 4691 7 of 10 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 .

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.
Energies 2020, 13, x FOR PEER REVIEW 7 of 10 By observing the time when Tl,C, I, and Th,G reach the stable state in Figure 4, it can be found that current I reaches the stable state earlier than Tl,C and Th,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 Tl,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 Qh,C changes along with it, which leads to the changes of Qh,G and Th,G. That is, I reaches a stable state earlier than Th,G.

Transient Analysis of Different Qs
The results of the steady-state analysis show that the size of Qs has a great influence on the coldend temperature Tl,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 Tl,C. Therefore, the variation of Tl,C and I was analyzed under different Qs.
Transient analysis was carried out of the integrated thermoelectric generation-cooling system with Qs of 6.05 W, 9.68 W, and 14.52 W, and the variation trend of Tl,C and I over time under different intrinsic heat sources Qs was obtained, as shown in Figure 5. It can be seen that in Figure 5, under a different intrinsic heat source Qs, 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 Tl,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 Ql,C, while the Peltier heat under a large current is larger than the cooling load Ql,C.
Although the change of Qs does not change the trend of Tl,C and I, it does change the maximum value and corresponding time of Tl,C. As shown in Figure 5, when Qs was 6.05 W, 9.68 W, and 14.52 W, the maximum value of Tl,C was 312.48 K, 309.02 K, and 306.51 K, respectively. In addition, the time corresponding to the maximum of Tl,C was 2.6 s, 1.7 s, and 1.1 s, respectively. In other words, with the increase of Qs, Tl,C can reach its maximum value in a shorter time. In order to further analyze this result, transient results under different Qs conditions were processed. The curve of the maximum value of the Tl,C and its corresponding time change with Qs were obtained as shown in Figure 6.
In Figure 6, Tlmax represents the maximum value of Tl,C and tm represents its corresponding time. 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. As can be seen from Figure 6, the maximum value of the cold-end temperature Tlmax and the corresponding time tm decrease with the increase of Qs.
According to formula (3), Tl,C is determined by the magnitude of Peltier heat and Fourier heat. In the case of constant Ql,C, the current I directly determines the change of the TEC's temperature. Figure 7 shows the variation curves of Tl,C and I within 0-4 s under different Qs.
It can be seen from Figure 7 that, although the maximum value and corresponding time of Tl,C are different under different Qs, the corresponding current is between 0.79 A and 0.94 A. Combined with the previous analysis of the Tl,C variation trend, the current I directly determines the change of Tl,C. Therefore, the increasing rate of current I determines the magnitude of tm. 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 Qs 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 Tl,C reach its maximum value.

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. 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 . As can be seen from Figure 6, the maximum value of the cold-end temperature Tlmax and the corresponding time tm decrease with the increase of Qs.
According to formula (3), Tl,C is determined by the magnitude of Peltier heat and Fourier heat. In the case of constant Ql,C, the current I directly determines the change of the TEC's temperature. Figure 7 shows the variation curves of Tl,C and I within 0-4 s under different Qs.
It can be seen from Figure 7 that, although the maximum value and corresponding time of Tl,C are different under different Qs, the corresponding current is between 0.79 A and 0.94 A. Combined with the previous analysis of the Tl,C variation trend, the current I directly determines the change of Tl,C. Therefore, the increasing rate of current I determines the magnitude of tm. 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 Qs 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 Tl,C reach its maximum value.

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. 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.

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.