Discrepancy between Constant Properties Model and Temperature-Dependent Material Properties for Performance Estimation of Thermoelectric Generators

The efficiency of a thermoelectric (TE) generator for the conversion of thermal energy into electrical energy can be easily but roughly estimated using a constant properties model (CPM) developed by Ioffe. However, material properties are, in general, temperature (T)-dependent and the CPM yields meaningful estimates only if physically appropriate averages, i.e., spatial averages for thermal and electrical resistivities and the temperature average (TAv) for the Seebeck coefficient (α), are used. Even though the use of αTAv compensates for the absence of Thomson heat in the CPM in the overall heat balance, we find that the CPM still overestimates performance (e.g., by up to 6% for PbTe) for many materials. The deviation originates from an asymmetric distribution of internally released Joule heat to either side of the TE leg and the distribution of internally released Thomson heat between the hot and cold side. The Thomson heat distribution differs from a complete compensation of the corresponding Peltier heat balance in the CPM. Both effects are estimated quantitatively here, showing that both may reach the same order of magnitude, but which one dominates varies from case to case, depending on the specific temperature characteristics of the thermoelectric properties. The role of the Thomson heat distribution is illustrated by a discussion of the transport entropy flow based on the α(T) plot. The changes in the lateral distribution of the internal heat lead to a difference in the heat input, the optimum current and thus of the efficiency of the CPM compared to the real case, while the estimate of generated power at maximum efficiency remains less affected as it is bound to the deviation of the optimum current, which is mostly <1%. This deviation can be corrected to a large extent by estimating the lateral Thomson heat distribution and the asymmetry of the Joule heat distribution. A simple guiding rule for the former is found.


Introduction
Thermoelectric generator (TEG) materials convert a certain fraction of the heat passed through them into useful electrical power, as the charge carriers (holes/electrons) absorb the thermal energy and move from the hot side to the cold side, carrying entropy [1,2]. The transport entropy flux related to the convective heat transport is given by αj, with the Seebeck coefficient α(T) and current density j. Typically, a thermoelectric (TE) module consists of a series of pn leg pairs (thermocouples), electrically connected in series and thermally in parallel [3]. In steady-state conditions, the exact − jT dα dT ∂T ∂x = −ρ(T) j 2 (1) where the thermal conductivity κ, the electrical resistivity ρ and α are the three main temperature-dependent thermoelectric properties. Here, ∂ ∂x · κ(T) ∂T ∂x = κ(T) ∂ 2 T ∂x 2 + dκ dT ∂T ∂x 2 corresponds to the (negative) divergence of the Fourier heat flux, i.e., its local change; jT dα dT ∂T ∂x corresponds to the local Thomson heat absorption driven by the change of the convective entropy flux αj related to the temperature dependence of α(T), and ρ(T) j 2 corresponds to the local Joule heat dissipation. With a typical TE material, with κ(T) falling with T and the amount of α(T) rising with T, Thomson heat will be released and Fourier heat flow will grow from the hot to the cold side along a TE leg in TEG operation, where the current flow is driven by the thermo-voltage generated by the leg. Equation (1) is a second-order non-linear partial differential equation, which can be solved using numerical methods like finite element methods (FEMs) [5,6], finite volume methods (FVMs) [7][8][9] or finite difference methods (FDMs) [6]. However, these solution methods are costly and time-consuming.
On the other hand, when assuming constant properties of the TE properties, an approximate solution can be found analytically, as suggested by Ioffe [1]. This solution by the constant property model (CPM) involves a discrepancy from the exact results due to the underlying simplification. Moreover, the choice of the averaged constant properties to be obtained from the actual temperature-dependent data is not straightforward. As can be seen from Equation (1), the Thomson heat vanishes when the Seebeck coefficient (and with that the convective entropy flux αj) remains constant. Various models corrected the CPM to compensate for this "missing Thomson heat" [10][11][12][13][14][15][16] have been proposed. Meanwhile, Sandoz et al. [17] attempted to explain the use of the T-averaged Seebeck coefficient in predicting exact power in the CPM mathematically, but did not recognize the importance of the asymmetry in heat distribution for the prediction of efficiency.
In a previous study [18], on the physically appropriate choice of averages in the CPM, we highlighted that spatial averages (SpAv) for resistivities (electrical and thermal) and temperature averaging (TAv) for the Seebeck coefficient are essential for a meaningful CPM estimate. However, there is still a remaining deviation due to unconsidered local redistribution of internal heat release or absorption and of thermal conduction in the CPM, which is linked to a change in the T profile T(x) [1,12,18,19]. Here, we will analyze the individual heat contributions exemplarily for six representative thermoelectric materials that we considered previously [18] plus PbTe [20], as this is one of the best TE materials in practice and shows an especially large deviation between CPM and exact results.
Initially, the effect of the T dependence of each of the TE properties, α, ρ and κ, leading to locally shifted heat release and transport over the TE element, for performance estimation, is studied separately. Calculated maximum efficiency in the full temperature-dependent case, η max , is compared with tailored model materials, in order to separate and quantify the individual contributions. Model materials are defined by setting one or two of the three TE properties constant at its respective average while keeping the other properties T dependent. Next, we explain the physical origin of a relevant part of the discrepancy between CPM results and the real situation using a schematic plot of the convective entropy flux derived from an α(T) graph, alongside showing that the net Peltier/Thomson heat is correctly considered by the CPM when appropriate temperature averaging is used for α(T). Marked areas in the entropy flux diagram quantify the exchange of Peltier and Thomson heat, and with that, a correction for the related deviation in CPM efficiency estimation, dη max = η max −η max CPM η max CPM , is suggested and demonstrated.

Role of the T Dependence of Material Properties in Performance Estimation
Since a generalized temperature dependence study for all types of T dependence is quite elaborate, a comparative study based on seven well-known and representative TE materials [20][21][22][23][24] was conducted. To understand the role of the T dependence of each of α, ρ and κ in performance estimation, the calculated maximum efficiencies when all properties are considered as T dependent (referred to as "real case" or "exact" from now on) were compared with the calculated efficiencies of model materials. These model materials have the same T dependence as the real materials for one or two of the three thermoelectric transport properties, while the remaining properties are kept constant; these materials are denoted as two temperature-dependent property (2TD) materials and 1TD materials, respectively. The constants used to define the model materials were obtained using the spatial averages (SpAv; for electrical and thermal resistivity) at a current density corresponding to the maximum efficiency of the real material and the temperature average (TAv; for the Seebeck coefficient). The SPAv and TAv of a T-dependent quantity p for a hot side temperature T h and a cold side temperature T c are given by [1,12,18,25] where ∆T = T h − T c and L are the length of the TE leg. The exact efficiency using T-dependent properties was obtained using the 1D solution algorithm developed in [18] by calculating Here, P is the output power, V is the net output voltage which is given by the Seebeck voltage where A is the area of the TE leg and ρ SpAv = 1 L L 0 ρ(T(x))dx. I = jA is the current passing through the TE material due to the generated voltage. The efficiency (η) is given by the ratio of output power to the input heat flow (Q in ) as in Equation (5), where Q in is given by Q in consists of the Fourier heat flow −κ h ·A· dT dx h (including the fraction of Joule and Thomson heat contributions released in the leg which is flowing to the hot side) plus the Peltier heat (I·α h · T h ) absorbed at the hot side. The suffix h indicates the hot side values, i.e., κ h = κ(T h ) and α h = α(T h ). As the spatial averages depend on T(x), which in turn varies with current, they were formed pre-assuming the optimum current of the real materials. For brevity, the efficiency was also calculated at the optimum current of the real material. The optimum current in the numerical calculation was obtained by finding the current where dη dI becomes zero. The relative deviation (RD) of the calculated maximum efficiency between the 2TD model materials and the real materials, δη model max = , is shown in Figure 1a. Here, and in the following, for brevity, we will use δ and d to denote a relative and absolute deviation, respectively. The comparison shows how strongly each of the contributing T dependences alone would shift efficiency. Obviously, the T dependence of ρ will affect the calculated efficiency to a lower extent than Entropy 2020, 22, 1128 4 of 18 α(T) and κ(T) will do for some materials (middle section of Figure 1a); the asymmetry of Joule heat generation mostly plays a minor role. However, this does not hold for all materials and it does not mean that the RD between the CPM and a real material due to asymmetric distribution of Joule heat, Q h J , would be insignificant, as all of the three identified effects will act simultaneously when comparing the CPM and the real case. Although the effects of the T dependence of α(T) and κ(T) are much larger for some materials, they often partly cancel each other. A comparison of the real Joule heat partial T profiles in Figure 1b shows a considerable asymmetry, in correlation to the deviations in the ρ(T) = const. case for SnSe and PbTe (Figure 1a, mid); however, the RD contribution related to the profiles in Figure 1b is larger as they contain an asymmetry due to the asymmetry of axial heat conduction linked to κ(T), in addition to the asymmetry of Joule heat generation which alone is represented by Figure 1a. Calculation of the partial T profiles is explained in Appendix A.2. It should be noted that unlike for α(T), where the absence of the T dependence means an absence of Thomson heat, the absence of the T dependence of ρ just means that there is no local asymmetry in Joule heat generation, whereas the amount of Joule heat that appears remains unchanged. Both symmetrically or asymmetrically released Joule heat will contribute, together with Thomson heat, to the effect of a T dependence of κ(T) that consists in shifting the distribution of the inner reversible and irreversible heat towards the hot and cold sides. Accordingly, the magnitude of the effect of a T dependence of κ(T) will scale with the total amount of inner heat.
heat, max J = maẋ J ḣJ h , would be insignificant, as all of the three identified effects will act simultaneously when comparing the CPM and the real case. Although the effects of the T dependence of ( ) and ( ) are much larger for some materials, they often partly cancel each other. A comparison of the real Joule heat partial profiles in Figure 1b shows a considerable asymmetry, in correlation to the deviations in the ( ) = . case for SnSe and PbTe (Figure 1a, mid); however, the RD contribution related to the profiles in Figure 1b is larger as they contain an asymmetry due to the asymmetry of axial heat conduction linked to ( ), in addition to the asymmetry of Joule heat generation which alone is represented by Figure 1a. Calculation of the partial T profiles is explained in Appendix A.2. It should be noted that unlike for ( ) , where the absence of the T dependence means an absence of Thomson heat, the absence of the T dependence of just means that there is no local asymmetry in Joule heat generation, whereas the amount of Joule heat that appears remains unchanged. Both symmetrically or asymmetrically released Joule heat will contribute, together with Thomson heat, to the effect of a T dependence of ( ) that consists in shifting the distribution of the inner reversible and irreversible heat towards the hot and cold sides. Accordingly, the magnitude of the effect of a T dependence of ( ) will scale with the total amount of inner heat.
When or is kept constant, there can be large discrepancies, as seen from the scatter in the left and right section of Figure 1a. Switching off Thomson heat results in a change from non-constant to constant convective entropy flux linked to a different partition of reversible (Peltier + Thomsonbound) heat to both sides of the leg. When setting ( ) = ., net Fourier heat transmitted does not change as the thermal resistance of the leg is fixed by the definition of the SpAv. Rather, the observed differences are merely due to a changed lateral distribution of Thomson and Joule heat. Comparing this to Figure 2a reveals that a large RD for ( ) = . correlates to strongly nonlinear T profiles linked to ( ) (see profiles for = 0); see also Appendix A.2, Figure A1a, where SnSe, Bi2Te3 and PbTe have significantly different h and c and Figure A2a, showing that the weight of Joule and Thomson heat to in is comparably large for these materials. The dominating effect of the T dependence of and on the estimated performance is also seen by comparing the T profiles of the model cases with the real temperature profile of n-type Mg2(Si,Sn) (referred to as n-Mg2X), Figure 2b. All profiles are calculated for the optimum current for maximum efficiency of the real material. Here, in addition to the 2TD materials, 1TD materials were also involved.
( ) and ( ) play a dominating role in the shaping of the temperature profile, which is reflected by the closeness of the ( ) ≠ ., ( ) ≠ . case to the real material. When α or κ is kept constant, there can be large discrepancies, as seen from the scatter in the left and right section of Figure 1a. Switching off Thomson heat results in a change from non-constant to constant convective entropy flux linked to a different partition of reversible (Peltier + Thomson-bound) heat to both sides of the leg. When setting κ(T) = const., net Fourier heat transmitted does not change as the thermal resistance of the leg is fixed by the definition of the SpAv. Rather, the observed differences are merely due to a changed lateral distribution of Thomson and Joule heat. Comparing this to Figure 2a reveals that a large RD for κ(T) = const. correlates to strongly non-linear T profiles linked to κ(T) (see T profiles for j = 0); see also Appendix A.2, Figure A1a, where SnSe, Bi 2 Te 3 and PbTe have significantly different κ h and κ c and Figure A2a, showing that the weight of Joule and Thomson heat to Q in is comparably large for these materials.
The dominating effect of the T dependence of κ and α on the estimated performance is also seen by comparing the T profiles of the model cases with the real temperature profile of n-type Mg 2 (Si,Sn) (referred to as n-Mg 2 X), Figure 2b. All profiles are calculated for the optimum current for maximum efficiency of the real material. Here, in addition to the 2TD materials, 1TD materials were also involved. α(T) and κ(T) play a dominating role in the shaping of the temperature profile, which is reflected by the closeness of the α(T) const., κ(T) const. case to the real material.
(red dotted line) deforms the T profile considerably but hardly increases the heat input (Equation (6)) compared to the real material, as the SpAv of ( ) maintains an unchanged thermal resistance of the TE leg. We can conclude that replacing the dependence of ( ) and ( ) by adequate constants will, although significantly changing the profile, influence the inflowing heat and thus efficiency to a much lower extent due to compensating effects. The RD of CPM efficiency in effect arises mainly from a redistribution of internal Joule and Thomson heat due to considerable deformation of the T profile by neglecting the T dependence of ( ) and ( ) and local redistribution of reversible heat generation as a consequence of neglect of the T dependence of the convective entropy flux.
The procedure to calculate the partial profiles is described in Appendix A.3. From the close coincidence of combinations of the real partial T profiles to the T profiles of the 1TD and 2TD model materials, as evident from Figure 2b, we can conclude that the contributions from each of the effects (Thomson heat, Joule heat, T dependence of ) to the total ( ) behave in good approximation and are independent and additive (a small note on this is given in the Appendix A.1.). The reason for the overall weak cross-coupling between the contributing effects is the small amplitude of the partial profiles ∆ Joule , ∆ Thomson , ∆ ( ) compared to the overall ∆ but also the fact that ∆ Thomson and ∆ ( ) often partially compensate. Therefore, the profiles of a real material and the CPM may also be quite close to each other for some materials. It is evident that the shape of ( ) and ( ) affects the temperature profile much more than that of ( ) but this does The effects of the 2TD cases on the overall inflowing Fourier heat and thus on the efficiency of n-Mg 2 X from Figure 1a (red dots) can be discussed in terms of the hot side slopes of the corresponding temperature profiles (red lines) in Figure 2b when comparing between cases with the same κ(T). The downward dT dx h for the 2TD material with α(T) = const.(red solid line) indicates an increase in the inflowing Fourier due to missing Thomson heat, compared to the actual case (dark green line).
Simultaneously, but only partly compensated in the Q in balance by missing Thomson heat, less Peltier heat is absorbed at the hot side and therefore the efficiency is overestimated (Figure 1a left side, red dot). The 2TD κ(T) = const.(red dotted line) deforms the T profile considerably but hardly increases the heat input (Equation (6)) compared to the real material, as the SpAv of κ(T) maintains an unchanged thermal resistance of the TE leg. We can conclude that replacing the T dependence of α(T) and κ(T) by adequate constants will, although significantly changing the T profile, influence the inflowing heat and thus efficiency to a much lower extent due to compensating effects. The RD of CPM efficiency in effect arises mainly from a redistribution of internal Joule and Thomson heat due to considerable deformation of the T profile by neglecting the T dependence of κ(T) and α(T) and local redistribution of reversible heat generation as a consequence of neglect of the T dependence of the convective entropy flux.
When comparing the 1TD and 2TD model materials, additionally a shift of the SpAv values of ρ and κ as a consequence of different T profiles, as well as coupling effects among the individual contributions, play a role, but only to a very minor extent, as proven by the close coincidence of their profiles to combinations of the individual partial T profiles of the real material, see Figure 2b (pink and cyan lines). The latter represent the physical contributions to the real temperature profile, ∆T Joule , ∆T Thomson and ∆T κ(T) and are plotted by symbols and lines in Figure 2b. They sum up, The procedure to calculate the partial profiles is described in Appendix A.3. From the close coincidence of combinations of the real partial T profiles to the T profiles of the 1TD and 2TD model materials, as evident from Figure 2b, we can conclude that the contributions from each of the effects (Thomson heat, Joule heat, T dependence of κ) to the total T(x) behave in good approximation and are independent and additive (a small note on this is given in the Appendix A.1.). The reason for the overall weak cross-coupling between the contributing effects is the small amplitude of the partial T profiles ∆T Joule , ∆T Thomson , ∆T κ(T) compared to the overall ∆T but also the fact that ∆T Thomson and ∆T κ(T) often partially compensate. Therefore, the T profiles of a real material and the Entropy 2020, 22, 1128 6 of 18 CPM may also be quite close to each other for some materials. It is evident that the shape of α(T) and κ(T) affects the temperature profile much more than that of ρ(T) but this does not mean that the asymmetry of Joule heat distribution between the hot and cold side would contribute insignificantly to the difference of the inflowing heat between the CPM case and a real material. The redistribution of Joule heat affects the maximum efficiency to a relevant extent along with the redistribution of Thomson heat. Thus, we can split the RD of the maximum efficiency according to the physical origin-redistribution of Peltier-Thomson heat and Joule heat-as δη max = Depending on the slope ratio of κ(T) and α(T), the efficiency discrepancy due to Joule heat asymmetry, δη maxJ , will vary considerably between different materials and may change sign from case to case, as observed in [18]. Now, let us proceed to understand in more detail how the absence of Thomson heat in the CPM will affect the efficiency calculation. We will see that it is partially and usually not entirely compensated by the difference in Peltier heat between a real material and its CPM approximation.

Peltier-Thomson Heat Balance and the Resulting Uncertainty in CPM Efficiency
Consider a TE material with constant κ and a linearly increasing α(T) curve (which is typical for a TE material below the peak zT temperature), as schematically shown in Figure 3. In a TE material under current flow, the convective entropy flux is given by  Depending on the slope ratio of ( ) and ( ), the efficiency discrepancy due to Joule heat asymmetry, max J , will vary considerably between different materials and may change sign from case to case, as observed in [18]. Now, let us proceed to understand in more detail how the absence of Thomson heat in the CPM will affect the efficiency calculation. We will see that it is partially and usually not entirely compensated by the difference in Peltier heat between a real material and its CPM approximation.

Peltier-Thomson Heat Balance and the Resulting Uncertainty in CPM Efficiency
Consider a TE material with constant and a linearly increasing ( ) curve (which is typical for a TE material below the peak zT temperature), as schematically shown in Figure 3. In a TE material under current flow, the convective entropy flux is given by ̇( ) = ( ). Hence, in a TE leg with a current flow I, the convective entropy flow ̇( ) = ( ) is directly linked to the temperature dependence of the Seebeck coefficient.

S(T) curve (indicated by slant lines) is
where τ = T dα dT is the Thomson coefficient. This area represents the net Thomson heat generated in the TE leg, . Q τ , which is directly linked to the variation of the convective entropy flow over the leg. The reversible heat balance shows that the loss of Peltier heat in the sample equals released Thomson heat plus produced gross electrical power.
. Q τ and P 0 are counted here as positive when going out of the system. Part of the Thomson heat will flow back, as a contribution to the overall Fourier heat flow, to the hot side. For simplification we assume that Thomson heat that is released at any point in the leg will flow out to the closer side. This is physically not strict but sufficient to qualitatively illustrate the relevant effect of undercompensation of the difference in Peltier heat exchanged at the hot side in a real material compared to the CPM by Thomson heat flowing back to the hot side, i.e., compensation of d Tdα. The relevant question on the Seebeck value α τ,ex , from which the integration gives the correct amount of . Q τ,h (and its corresponding temperature T τ,ex with α τ,ex = α(T τ,ex )), will be touched on below.
In the CPM, the Peltier heat at the hot side is given by IαT h , while at the cold side it is IαT c , where α = α TAv is the temperature average of α(T) (see Equation (2)). Therefore, the following equation holds: i.e., Peltier heat is completely balanced by electrical production. From Equations (9) and (10), it is obvious that globally the explicit absence of Thomson heat in the CPM is taken care of correctly by the use of temperature averaged α in the CPM, i.e., With this choice of α as the CPM value, the gross power generated is exactly the same in the CPM as in the real material, at the same current. On the other hand, it implies that, typically, considerably less Peltier heat is absorbed at the hot side in the CPM case than in reality, whereas back-flowing Thomson heat partly compensates the actually higher Peltier heat intake. Figure 3 visualizes with the green triangle that this compensation is incomplete, i.e., d Accordingly, more Thomson heat is leaving at the cold side. It is evident that this holds not only for a linear but also for a left-or right-hand bowed Seebeck curve.
In a less typical case with strongly asymmetric heat conduction, i.e., κ(T) strongly increasing with T, or if α(T) forms a significant maximum, this typical tendency could reverse, but mostly it leads to underestimation of the inflowing heat in the CPM case Q in CPM and hence to overestimation of the efficiency by the CPM. With p-Mg 2 X, a particular example is given in Appendix A.3.2 ( Figure A2c) where, with α(T) weakly changing between T c and T h but peaking inside, this compensation can also be almost perfect, or, as for SnSe (Figures 4, 5b and 6), overcompensation may even occur. Figure 4a shows that for most materials, opt, changes for about 1% or less and, consequently, also the deviation of the output power, remain small. However, for PbTe, opt, reaches 10%. Then the deviation of output power, max , may grow in absolute amount to be as large as ̇i n , doubling its effect. Whereas the contribution to ̇i n , due to opt, usually remains insignificant, it becomes relevant for PbTe where it compensates half of ̇i n related to the distribution of inner heat at an unchanged current, dQ̇π τ h (I opt,η ) − dQ̇J h (I opt,η ), see Equation (12) and black stars in Figure 4a.  is plotted and shows that both effects reach the same order of magnitude. Typically, both contributions partly compensate. Furthermore, no general behavior can be observed in their mutual relation over the materials, as in some cases clearly one effect dominates, in others the other. As seen from Figure 1b, usually, more Joule heat is released to the hot side than to the cold side in a real material, whereas there are symmetric amounts in the CPM case. This contributes to an underestimation of the efficiency in the CPM case, max J > 0. On the other hand, as explained, the Peltier-Thomson balance tends to an overestimation, max < 0, thus, both effects counteract and partially compensate. From Figure 4a, it can be seen that the CPM overestimates the efficiency compared to the real case for all selected materials except Bi2Te3, which has an exceptionally higher be neglected as it cancels out completely. Only the intervals outside, ( c ; τ,ex ) or ( τ,ex ; h ) , have to be considered. Among both intervals, the side has to be chosen where ( ) is a monotonous function in the relevant temperature interval, where it is closer to linearity, and possibly where τ,ex is closer to h or c .
Applying Equation (11) accordingly to the chosen interval, the integration for ̇ can be substituted by one for 0 , e.g., for the cold side: This facilitates practical execution as ( ) is mostly known as a low-order polynomial, thus integration could be done analytically.
If the Thomson profile is not known, half of the leg length, 2 , can be taken as a first guess of the position for the calculation of ̇, h . The corresponding temperature is marked in the diagrams. This can be a quite good estimate when the Thomson profile is close to symmetric, as for PbTe (see Figure A2b), but may fail greatly when Thomson heat is strongly asymmetric, as for SnSe. On the contrary, an entropy consideration of Thomson heat in the TE leg (see Appendix A.4.) leads to a rule of thumb for τ,ex that is Indeed, it applies well for all example materials involved here. With this rule, approximation of ̇, h is facilitated considerably, as just a crossing point of (T) with its TAv has to be found.   Equation (12)) and a first guess by the position.
Only occasionally, e.g., when ( ) is close to linear, the 2 position works well for correction but fails for most materials as it does not take into account the asymmetry of heat sources and heat conduction. Similarly, models suggesting half of the Thomson heat on either side for correcting the CPM results [14-16,26,27] will mostly not work sufficiently. The correction employing the Q in , will accordingly have a lower slope in reality than for the CPM, equivalent to a lower maximum position I opt,η . Thus, usually, the CPM will overestimate the optimum current, δI opt,η = I opt,η −I CPM opt,η I CPM opt,η < 0, and hence will overestimate output power at maximum efficiency (δP η max < 0), which adds to the overestimate of maximum efficiency: δη max = δP η max − δ . Q in , amplifying the effect of δ . Q in (see Figure 4a). Hence, for a quantitative analysis, we have to consider three contributions to the (absolute) deviation of where, similar to the outflowing Thomson heat, outflowing Joule heat is also counted as positive and d . Q h J is due to the Joule heat asymmetry at the hot side. Asymmetry of Joule heat distribution and heat conduction will, with falling κ(T), as for PbTe and SnSe, favor heat release to the cold side. This will likewise contribute to a higher . Q in and steeper . Q in (I), amplifying the same trend as for reversible heat, or will counteract it with rising κ(T). Thus, asymmetry of Joule heat distribution will add to the mispoint in I opt,η CPM . Figure 4a shows that for most materials, I opt,η changes for about 1% or less and, consequently, also the deviation of the output power, remain small. However, for PbTe, δI opt,η reaches 10%. Then the deviation of output power, δP η max , may grow in absolute amount to be as large as δ Q h πτ is plotted and shows that both effects reach the same order of magnitude. Typically, both contributions partly compensate. Furthermore, no general behavior can be observed in their mutual relation over the materials, as in some cases clearly one effect dominates, in others the other.
As seen from Figure 1b, usually, more Joule heat is released to the hot side than to the cold side in a real material, whereas there are symmetric amounts in the CPM case. This contributes to an underestimation of the efficiency in the CPM case, δη maxJ > 0. On the other hand, as explained, the Peltier-Thomson balance tends to an overestimation, δη maxπτ < 0, thus, both effects counteract and partially compensate. From Figure 4a, it can be seen that the CPM overestimates the efficiency compared to the real case for all selected materials except Bi 2 Te 3 , which has an exceptionally higher κ h compared to the cold side ( Figure A1a in Appendix A.2) together with high Joule release (Figure 4b) and almost compensation of the Peltier-Thomson balance. Thus, the Joule contribution dominates, leading to an underestimation of the efficiency. Additionally, SnSe behaves somewhat differently from the general trend, with a falling α(T) curve (Appendix A.2 Figure A1b) and the over-resistivity at the cold side (Appendix A.2 Figure A1c). Moreover, κ h is much lower than κ c . As an effect, Joule heat is preferentially led to the cold side; consequently, hot side Joule heat is greatly overestimated in the CPM (Figure 4b), but as the relative contribution of Joule heat to Q in is small (Figure A2a), the resulting trend towards the overestimation of performance in the CPM remains moderate. On the other hand, as seen from Appendix A.3.2 Figure A1b, Thomson heat is absorbed in the leg as α(T) for SnSe is a falling curve and is mainly bound to the hot side. As seen from Figure 4b, for SnSe, the hot side Peltier-Thomson heat will, unlike for most of the other materials, be overestimated by the CPM. However, the resulting underestimation of efficiency in the CPM will be overcompensated by the counteracting Joule heat distribution.
The first four materials in our list (see Figure 4a) show a minor discrepancy of the CPM with reality. Although Joule heat asymmetry is contributing comparably, from case to case, the dominating source of discrepancy is mostly the uncompensated Peltier heat according to Equation (11)). It is particularly relevant in the cases of n-Mg 2 X, Mg 2 Si and PbTe, which have larger Thomson contributions ( Figure A2a), leading to larger discrepancies of the CPM efficiency estimate.

Refining the CPM Efficiency Estimate
Having identified the effects causing a systematic uncertainty in the CPM efficiency estimation, they can be accordingly corrected.
We want to analyze how this can be done practically for the Thomson contribution, δ . Q h πτ , by calculating the uncompensated Peltier heat at the hot side. Therefore, we discuss the approach for example materials with dissimilar α(T) characteristics.
The values of . Q π,h and . Q CPM π,h are known from T h , α h and α, for a given current, where, as a first approximation, I CPM opt,η is used. We have seen that the Thomson heat flowing to the hot side is strictly calculated from the partial T profile ∆T Thomson h . We apply this route to form a reference for an approximate estimation to be developed and, because of this, we omit a numerical calculation of exact T profiles. As derived from Equation (10), we obtain the uncompensated Peltier-Thomson heat from d CPM . Neglecting any deviation of current, this can be illustrated in the α(T) diagram based on our interpretation of areas by amounts of reversible heat, see Figure 3. Thus, we aim for a good approximation of the green marked area in Figure 3 by an appropriateand simple approximation. The problem splits into two aspects: finding the temperature T τ,ex above which the inner Thomson heat is conducted to the hot side and finding a close approximation of the integral. As α(T) may be quite different (see Figure A1b), we meet various situations, represented by different ∆T Thomson (x) temperature profiles ( Figure A2b), among them typical ones with a single maximum according to Thomson heat flowing out to both sides, but also less typical ones with a single minimum (Thomson heat flowing in from both sides) or even two extrema (for Bi 2 Te 3 ) where Thomson heat is released to the cold side but absorbed from the hot side. A rule to treat all of the cases likewise is needed. Figure 5a,b and Figure A2c,d accordingly show scenarios where α(T) contains almost linear intervals along with strongly bowed ones, where α(T) is monotonous or contains a maximum, where α h and α are far from each other or close together or where α(T) crosses the α horizontal once or twice. The position of the extrema (maxima or minima) of ∆T Thomson (x) is marked in each diagram by a brown line. Accordingly, the area corresponding to the uncompensated heat might be more complex than is shown in Figure 3, e.g., see Figure 5a. The area to the left of the α(T) curve to the α-axis from this point up to the hot side α h (marked by a red border) represents . Q τ,h . The fact that the respective area also contains negatively counted parts when α(T) goes through a maximum is also taken into account. Accordingly, the upper slim boat-shaped area in Figure 5a counts as negative; symbolically, it is mirrored in the green area.
However, in such a case, the integration can be simplified, switching from the hot to the cold side, as Q τ,c and with Equation (10), Q π,h − . Q π,c − P 0 . Note that if there are two extrema of ∆T Thomson (x), then we have two T τ,ex values where the Thomson heat between both can be neglected as it cancels out completely. Only the intervals outside, (T c ; T τ,ex ) or (T τ,ex ; T h ), have to be considered. Among both intervals, the side has to be chosen where α(T) is a monotonous function in the relevant temperature interval, where it is closer to linearity, and possibly where T τ,ex is closer to T h or T c .
Applying Equation (11) accordingly to the chosen interval, the integration for . Q τ can be substituted by one for P 0 , e.g., for the cold side: This facilitates practical execution as α(T) is mostly known as a low-order polynomial, thus integration could be done analytically.
If the Thomson T profile is not known, half of the leg length, L 2 , can be taken as a first guess of the position for the calculation of . Q τ,h . The corresponding temperature is marked in the diagrams. This can be a quite good estimate when the Thomson T profile is close to symmetric, as for PbTe (see Figure A2b), but may fail greatly when Thomson heat is strongly asymmetric, as for SnSe. On the contrary, an entropy consideration of Thomson heat in the TE leg (see Appendix A.4.) leads to a rule of thumb for T τ,ex that is Indeed, it applies well for all example materials involved here. With this rule, approximation of . Q τ,h is facilitated considerably, as just a crossing point of α(T) with its TAv has to be found. Figure 6 shows the remaining efficiency deviation, δη corr max , corrected by the uncompensated Peltier-Thomson heat calculated from the α(T) graph using the L 2 position, using T τ,ex according to the extremum (maximum) position of ∆T Thomson (x) but neglecting the current deviation δI opt,η , as well as corrected by the exact deviation d CPM . The efficiency estimate by the CPM is greatly improved when the ∆T Thomson (x) extremum position is used(red dots). Only occasionally, e.g., when α(T) is close to linear, the L 2 position works well for correction but fails for most materials as it does not take into account the asymmetry of heat sources and heat conduction. Similarly, models suggesting half of the Thomson heat on either side for correcting the CPM results [14][15][16]26,27] will mostly not work sufficiently. The correction employing the ∆T Thomson (x) peak position is close to the exact numerical correction for most materials as this position considers the asymmetry exactly. The difference between both cases is merely due to the change of the optimum current which is as yet unconsidered by the graphical correction. The remaining discrepancy is due to Joule heat asymmetry.
Whereas we have used exact numerical calculations to demonstrate the principle of the Thomson correction method and to show that the rule α(T τ,ex ) ≈ α holds well, the suggested practical procedure for the correction of d . Q h, I=const πτ described here, which is based on an analysis of the physical effects behind the deviation of CPM performance estimates, is limited to basic algebraic operations which can be instantaneously calculated by any table calculation software.

Conclusions
From the study of 2TD and 1TD model materials with one or two selected properties among α, ρ and κ set as constant, which results in both redistribution of heat between the hot and cold side of the element and the change of spatial averages, we see that in some examples, large deviations in efficiency δη model max arise as a consequence of considerable modification of the T profile. In comparison to the efficiency deviation between the CPM and real materials δη max which conserve the spatial property averages and are mostly below 2%, this shows that a change of spatial averages due to an arbitrary modification of the T profile may contribute a strong shift to the efficiency estimate. Thus, conservation of the leg's thermal and electrical resistance is essential for a valid efficiency estimate. However, the shift mainly remains low if only ρ(T) is switched to constant. Nevertheless, it cannot be concluded from this that the temperature dependence of the electrical resistivity plays a minor role in the efficiency estimation by the CPM. The 2TD and 1TD model materials lead to quite good approximations of the partial T profiles ∆T Joule (x), ∆T Thomson (x) and ∆T κ(T) (x).
It is shown that the deviation of a CPM-based efficiency estimate, δη max , is not just due to the absence of Thomson heat in the CPM, as the choice of the temperature average of α(T) as a CPM parameter mainly compensates for the absence of Thomson heat. Rather, the discrepancy in efficiency determination in the CPM is shown to be, to a major extent, due to the excess unaccounted heat at the hot side in the CPM δ . Q in , which usually leads to overestimation of performance, and, to a minor extent, due to a shift of the optimum current δI opt,η and, consequently, of the produced electrical power at maximum efficiency, δP η max . In most cases, the change of the optimum current is small. In materials with rising α(T), less of the released Thomson heat flows back to the hot side than would compensate for the reduced hot side Peltier heat absorption assumed by the CPM. This systematic undercompensation tends towards a higher actual heat intake at the hot side compared to the CPM, thus overestimating efficiency when the CPM is used. Asymmetry of Joule heat usually has an opposite influence but is overcompensated in most cases.
In order to correct for the Peltier-Thomson heat-related deviation δ . Q h πτ , a graphical illustration in terms of convective entropy flow based on the α(T) curve is given. It confirms that the rule for the splitting of Thomson heat to the sides α(T τ,ex ) ≈ α, which results from an entropy consideration, holds well. This enables a valid approximation of δ . Q h πτ with a simple algebraic procedure, omitting the exact numerical calculation of the temperature profile. Although a considerable deformation of the T profile caused by the T dependence of κ(T) is observed, it will affect the deviation between the real situation and its CPM approximation simply via a local shift of the thermal and electrical resistivity but will not explicitly contribute to the inflowing heat balance . Q in . In summary, the performance of a TE material does not only depend on its averaged material parameters but also on local asymmetry of Thomson and Joule heat, driven by the T dependence of the TE properties. In particular, Thomson heat can show highly asymmetric distribution. Thus, TE device efficiency can be varied beyond the averaged properties, represented by a figure of merit.  Since SnSe has much higher resistivity, the scale for it is given on the right y-axis. All the raw experimental data taken from the literature [20][21][22][23][24] were fitted with appropriate polynomials (usually 3rd or 4th order). For SnSe, a 9th order polynomial fit was used owing to the complex T dependence and hence shows an unusually high max. However, this does not affect the physics discussed and hence these fitted data were used throughout the manuscript. Table A1. Temperature range of analysis for all materials of Figure A1.  (1) to each of the other corresponding terms, assuming isothermal boundary conditions and fixing all coefficients in the equation according to the total profile ( ). Thus, solving for the respective partial profile reduces to a double integration, where the first step provides the total amount of each partial heat contribution to the thermal balance.

Material
As the partial profiles can have opposite signs in amplitude and partially compensate for many of the common TE materials (however, not always), the profiles of a real material and the CPM may be quite close, as in the example of n-type Mg2X, Figure 2b. Since SnSe has much higher resistivity, the scale for it is given on the right y-axis. All the raw experimental data taken from the literature [20][21][22][23][24] were fitted with appropriate polynomials (usually 3rd or 4th order). For SnSe, a 9th order polynomial fit was used owing to the complex T dependence and hence shows an unusually high zT max . However, this does not affect the physics discussed and hence these fitted data were used throughout the manuscript. Table A1. Temperature range of analysis for all materials of Figure A1. The partial T profiles are each found by equating κ(T) ∂ 2 T ∂x 2 in Equation (1) to each of the other corresponding terms, assuming isothermal boundary conditions and fixing all coefficients in the equation according to the total T profile T(x). Thus, solving for the respective partial T profile reduces to a double integration, where the first step provides the total amount of each partial heat contribution to the thermal balance.

Material Temperature Range of Analysis
As the partial T profiles can have opposite signs in amplitude and partially compensate for many of the common TE materials (however, not always), the T profiles of a real material and the CPM may be quite close, as in the example of n-type Mg 2 X, Figure 2b.
As both Joule and Thomson heat, after appearing inside the leg, will flow out, physically, as Fourier heat, we have to consider in this discussion the pure Fourier heat Q F,h = K∆T (with K = κ −1 −1 A/L), which is merely related to the thermal resistance of the leg and is constant along the leg, separately from the Joule-and Thomson-related contributions. Accordingly, Q in is composed of The real Joule-and Thomson-related contributions, − . Q τ,h and . −Q J,h , to the inflowing hot side heat are calculated by splitting the overall temperature profile T(x) into additive partial T profiles, each related to one of the individual physical contributions. Partial Thomson T profiles of example materials are plotted in Figure A2b. Q J,h , respectively. Figure A2a shows the relative contribution of each heat to . Q in : .
This comparison reveals that Joule and Thomson heat contribute about 1-5% to . Q in , usually flowing out, with their contributions being roughly of the same order. Figure  In order to illustrate example situations of the distribution of Peltier and Thomson heat along the leg, α(T) graphs for p-Mg 2 X and Bi 2 Te 3 are given in Figure A2c,d, respectively. Due to the bowed shape of the α(T) graph and relatively close values of α h to α c for p-Mg 2 X, the difference between . Q π,h and . Q CPM π,h is almost negligible, but . Q τ,h amounts to more than twice the amount of . Q π,h − . Q CPM π,h . Nevertheless, this did not affect the efficiency deviation δη max too much, as . Q τ,h is quite small in absolute terms. In the case of Bi 2 Te 3 , . Q CPM π,h is even higher than . Q π,h again due to the curved shape of α(T), affecting the position of α TAv . However, . Q τ,h almost completely compensates for this Peltier heat difference, keeping the influence on the efficiency deviation negligible.

Appendix A.4. Thomson Heat Distribution and Entropy
With the TEG leg, we discuss the entropy flow in a reversible system of Peltier heat transport and Thomson heat exchange which is running on a non-equilibrium temperature background mainly fixed by the continuous flow of Fourier heat. As released Thomson heat will be transported as Fourier heat but is small in relation to the Fourier heat background (see Figure A2a), which is driven by the temperature difference and the thermal resistance of the TE leg, we will treat the variation of the temperature profile by the conducted Thomson heat as insignificant for the following consideration.
In the steady state, the entropy of the system remains constant; there is a continuous entropy production by the dissipative heat transport from hot to cold and the balancing continuous entropy export by transmitted Fourier heat (plus a negligible fraction arising from outflowing Joule heat). Assuming ideal outer current leads with α = 0, there is no other entropy exchange at the hot and cold sides.
difference between ̇, h and ̇, h CPM is almost negligible, but ̇, h amounts to more than twice the amount of ̇, h −̇, h CPM . Nevertheless, this did not affect the efficiency deviation δη max too much, as ̇, h is quite small in absolute terms. In the case of Bi2Te3, ̇, h CPM is even higher than ̇, h again due to the curved shape of ( ) , affecting the position of TAv . However, ̇, h almost completely compensates for this Peltier heat difference, keeping the influence on the efficiency deviation negligible.  In the CPM, we have a constant convective entropy flow αI throughout the element, equal to the absorbed and released entropy rate αI by absorption and release of Peltier heat at the terminals. In a real material, the absorbed entropy rate α h I equals the convective entropy flow at the hot side, and, likewise, the amount of α c I at the cold side. The variation of α along the leg drives local Thomson heat production d . Q τ = T dα dT IdT = TIdα, contributing an entropy flow increment d . S = Idα. Thomson heat flows to the hot and cold sides and the related total entropy exchange is (α h − α c )I = ∆αI. It distributes by the fraction x h to the hot and cold sides: Driven by the gradient of the partial Thomson temperature profile, all Thomson heat released at one side of a maximum (or minimum) of this profile will be exchanged to this side of the leg. With the temperature T τ,ex of this position and its Seebeck coefficient α τ,ex = α(T τ,ex ), the shares of the entropy exchange which are bound to each of the sides are Q π − IV 0 , we can conclude that α τ,ex ∆T = V 0 = α∆T, thus α τ,ex = α (A5) This gives us a rule for the temperature intervals over which the Thomson heat is flowing to either side of the leg. Consequently, Thomson heat has to be integrated from the crossing point of the curve of the Seebeck coefficient α(T) with its temperature average α. As a reversible approximation, this result is approximate and not strict as we have neglected here that dissipative processes are involved when Thomson heat is conducted to the leg sides. Below we will analyze these changes and find that these are small, and thus the rule stated here on the position of α τ,ex , although not strict, is a good guide for estimates of the distribution of Thomson heat. Indeed, as observed by comparison to exact numerical calculations, this rule is almost perfectly fulfilled for all the example materials.
Within this reversible approximation, the Thomson heat flowing to the hot side is obtained as . Q τ,h = T h (α h − α)I. This would be equivalent to a complete compensation of the Peltier heat difference between reality and the CPM, i.e., the vanishing axial redistribution of reversible heat which is consistent with the simplifying assumption that the Thomson heat flowing to the outside is transmitted free of dissipation, i.e., equivalent to reversible heat. Here, the (additional) T gradient related to the flow of Thomson heat is neglected, whereas an underlying T profile related to an independent heat flow (here, the background of Fourier heat transfer) does, in effect, not contribute to its dissipation.
We will see below that this happens as Thomson heat flowing to different sides will contribute almost compensating shares to the entropy balance. What is neglected here is that the Thomson heat itself when flowing to the ends of the leg will dissipate, according to the slight shift of the inner T profile it is causing. Above, this T offset was separated and called the partial T profile due to Thomson heat, ∆T Thomson (x). Additionally, this omission will contribute to a weak deviation from the position rule α τ,ex = α. The dissipative part of the entropy transport to the sides of the leg is related to the T drop or step-up between the location where an increment of Thomson heat d . Q τ is released and the side temperature, T h or T c . The entropy increment is released over a segment of the leg with the T increment