**Discussion Panel**

[17] Shakouri A. Recent developments in semiconductor thermoelectric physics and materials. Annu. Rev. Mater. Res.. 2011;41:399–431. doi:10.1146/annurev-matsci-062910-100445 [18] Ouyang Z, Li D. Modelling of segmented high-performance thermoelectric generators with effects of thermal radiation, electrical and thermal contact resistances. Sci. Rep.

[19] Goupil C, Seifert W, Zabrocki K, Muller E, Snyder GJ. Thermodynamics of thermoelectric phenomena and applications. Entropy. 2011;3:1481–1517. doi:10.3390/e13081481

[20] Vargas A, Olivares MA, Camacho P. Thermoelectric system in different thermal and electrical configurations: its impact on the figure of merit. Entropy. 2013;15:2162–2180.

[21] Abdelkefi A, Alothman A, Hajj MR. Performance analysis and validation of thermoelectric energy harvesters. Smart Mater. Struct. 2013;22:095014. doi:10.1088/0964-1726/22/9/

[22] Nemir D, Beck J. On the significance of the thermoelectric figure of merit Z. J. Electron.

[23] Apertet Y, Ouerdane H, Glavatskaya O, Goupil C, Lecoeur P. Optimal working conditions for thermoelectric generators with realistic thermal coupling. EPL. 2012;97:28001.

[24] Vargas A, Olivares MA, Camacho P. Maximum power of thermally and electrically coupled thermoelectric generators. Entropy. 2014;16:2890–2903. doi:10.3390/e16052890

[25] Vargas A, Olivares MA, Mendez Lavielle F. Performance of composite thermoelectric generator with different arrangements of SiGe, BiTe and PbTe under different configura-

Mater. 2010;39:1897–1901. doi:10.1007/s11664-009-1060-4

tions. Entropy. 2014;16:2890. doi:10.3390/e17117387

2016;6:24123. doi:10.1038/srep24123. .

534 Thermoelectrics for Power Generation - A Look at Trends in the Technology

doi:10.3390/e15062162

doi:10.1209/0295-5075/97/28001

095014

#### **Non-Stationary Thermoelectric Generators Non-Stationary Thermoelectric Generators**

John G. Stockholm

Additional information is available at the end of the chapter Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/66421

#### **Abstract**

A review of theoretical publications on non-steady thermoelectrics is given. Review concerns different aspects of non-stationary and pulsed processes in thermoelectric materials and devices. Theoretical analysis of dynamic behaviour of thermoelectric devices, including analysis of small and large signals of thermoelectric generator, is given and details of concepts of quasi-equilibrium thermoelectricity are discussed as well. Special attention is paid to theoretical study of the non-routine regime of nonsteady thermoelectricity—fast-time dependence of thermoelectric properties when material or device is well out of equilibrium. Theoretical findings of fast-time dependence give reason to believe that it can increase the output electrical power of thermoelectric generator compared to stationary regime of operation. We also present experimental results obtained with first non-stationary thermoelectric generator prototype, which was designed for operation in fast-time dependence mode. Several research teams are presently making and testing devices to confirm that more electrical power can be obtained in AC mode (AC frequency about hundreds of kHz) than in DC mode. Descriptions with an analysis are given.

**Keywords:** thermoelectricity, electricity generation, quasi-stationary systems, out-ofequilibrium transient modes, ultra-fast conduction, microscopic processes

## **1. Introduction**

Operation of thermoelectric (TE) heat pumps used for cooling (TECs) and thermoelectric generators (TEGs) are closely related. Results of first theoretical studies of thermoelectric

© 2017 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2016 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

devices operating in non-steady regimes were published for thermoelectric coolers by Stilbans and Fedorovich in the year 1958 [1]. Then, Gray presented in the year 1960 [2] results of analytical consideration concerning small and large signals of thermoelectric generator. In 1961, Landecker [3] published results concerning fast transient behaviour of Peltier junctions during temperature measurements with thermocouples. Then, in the year 1971 Babin and Iordanishvili studied quasi-stationary thermoelectric generation and in the year 1983, Iordanishvili and Babin published a monograph [4] of non-stationary processes in thermoelectric and thermomagnetic energy conversion systems. It presents the fundamentals of pulsed current generation, pulsed thermoelectric generator, quasi-steady-state generator and single-pulse thermoelectric generator.

There is also non-steady, well-out-of-equilibrium, fast-time dependence. This subject has not been extensively developed and few papers were published on this subject. Such non-steady processes have the potential to transfer higher electrical powers, than in steady state. Currently ongoing experimental work is presented.

## **2. Theoretical aspects of pulsing**

#### **2.1. Two different types of non-stationary phenomena**

#### *2.1.1. Quasi-stationary thermoelectricity*

Quasi-stationary thermoelectricity, that is, quasi-equilibrium thermoelectricity, assumes a slight time dependence and is allowed for all the laws of the equilibrium, which means that this time dependence is only slow in comparison with the microscopic processes of heat transfer and electric conduction. This is considered in detail by Gray [2] and Iordanishvili and Babin [4]. Goupil [5, 6] studied non-stationary operation with Onsager approach and calculates for thermoelectric device the important parameters for transient regimes: time constant *τ* and capacitance *C*.

#### *2.1.2. Far from equilibrium fast-time dependence*

Thermoelectric device developed by Strachan, published by Aspen in the year 1992 [7, 8], was operated in MHz range; unfortunately, functional physics was never established and the unit operated in heat-pump mode and in electricity-generating mode, as well.

The first person who studied the fast-time dependence or well-out-of-equilibrium regime was Apostol [9]. His work can be considered as the beginning of new era in thermoelectric power generation. The physical mechanism is ultra-fast conduction, but this mechanism can be obtained only under certain circumstances, viz.:


## **2.2. Quasi-stationary thermoelectricity**

#### *2.2.1. Dynamic thermoelectric processes*

devices operating in non-steady regimes were published for thermoelectric coolers by Stilbans and Fedorovich in the year 1958 [1]. Then, Gray presented in the year 1960 [2] results of analytical consideration concerning small and large signals of thermoelectric generator. In 1961, Landecker [3] published results concerning fast transient behaviour of Peltier junctions during temperature measurements with thermocouples. Then, in the year 1971 Babin and Iordanishvili studied quasi-stationary thermoelectric generation and in the year 1983, Iordanishvili and Babin published a monograph [4] of non-stationary processes in thermoelectric and thermomagnetic energy conversion systems. It presents the fundamentals of pulsed current generation, pulsed thermoelectric generator, quasi-steady-state generator and

There is also non-steady, well-out-of-equilibrium, fast-time dependence. This subject has not been extensively developed and few papers were published on this subject. Such non-steady processes have the potential to transfer higher electrical powers, than in steady state. Currently

Quasi-stationary thermoelectricity, that is, quasi-equilibrium thermoelectricity, assumes a slight time dependence and is allowed for all the laws of the equilibrium, which means that this time dependence is only slow in comparison with the microscopic processes of heat transfer and electric conduction. This is considered in detail by Gray [2] and Iordanishvili and Babin [4]. Goupil [5, 6] studied non-stationary operation with Onsager approach and calculates for thermoelectric device the important parameters for transient regimes: time

Thermoelectric device developed by Strachan, published by Aspen in the year 1992 [7, 8], was operated in MHz range; unfortunately, functional physics was never established and the unit

The first person who studied the fast-time dependence or well-out-of-equilibrium regime was Apostol [9]. His work can be considered as the beginning of new era in thermoelectric power generation. The physical mechanism is ultra-fast conduction, but this mechanism can be


operated in heat-pump mode and in electricity-generating mode, as well.

single-pulse thermoelectric generator.

ongoing experimental work is presented.

**2. Theoretical aspects of pulsing**

*2.1.1. Quasi-stationary thermoelectricity*

constant *τ* and capacitance *C*.

*2.1.2. Far from equilibrium fast-time dependence*

obtained only under certain circumstances, viz.:

widen, during the time it takes to go around the circuit.

**2.1. Two different types of non-stationary phenomena**

538 Thermoelectrics for Power Generation - A Look at Trends in the Technology

The study of dynamic behaviour of thermoelectric devices was presented in 1960 by Gray [2]. The object is behaviour of TE energy conversion devices under dynamic conditions. TE generator and TE heat pump also are non-linear distributed-parameter systems. Therefore, it is difficult to exactly analyse their dynamic behaviour. In this study, the dynamic behaviour is investigated by means of linear small-signal-distributed-parameter models, for which useful analytical results are obtained. Analysis allows predicting the dynamic behaviour of thermoelectric device both in frequency-domain and in time-domain equalization algorithms. This can be used to assess quantitative effects of device parameters with regard to dynamic response of the device. It is possible to linearize system's equations in small-signal variations in order to obtain linear mathematical models, more amenable to analysis. This approach is of limited usefulness in computing large-signal dynamic behaviour, such as transient processes, when TE generator is turned on. The study of complete small-signal behaviour of TE generators is developed, the results are expressed in terms of frequency-domain transfer functions and procedures are developed, which allow direct calculation of time-domain transients for the case of input and disturbance time functions, which are Fourier transformable. Experimental observation of small-signal dynamic behaviour of simple heat pump is presented and it confirmed theoretical model. No measurements were made on TE generator, but the accuracy of description should be as good, because mathematical model and analytical techniques are identical.

Comparison of experimentally obtained large-signal static and small-signal dynamic behaviour of TE heat pump with corresponding theoretical results indicates a disagreement at large currents due to neglecting of temperature dependence of parameters.

Complete set of equations for the case of large-signal dynamic analysis for step changes in current is presented.

To conclude, all this work with many equations is the basis for all studies related to small and large signals and it is written in view for appropriate control of TE generators.

#### *2.2.2. Non-stationary processes in thermoelectricity*

The subject of non-stationary thermoelectricity is also presented in detail in 1983 by Iordanishvili and Babin [4]. Several very interesting examples of TE generators are presented, including the case of solid body fuel (SF), capable of producing uniform heat production throughout the volume.

Thus, SF must maintain solid-state properties after heat production process (combustion).

Using the model of semi-infinite rods, assume that at time *t* = 0 there is instantaneous heat production, after which SF temperature becomes equal to *Tω* (initial temperature of two bodies is *T*0). The system temperature as a function of position and time was calculated and the graph is given (**Figure 1**). Such a generator has high speed and stability of the output characteristics.

**Figure 1.** Heat flux *q* versus time *t*, is the average effective heat flow during time *t, q*0 is the average effective heat flow in stationary mode, /0 ≈ 3.

This graph *q*(*t*) shows that in this model, due to thermal conductivity, average, during time *t*, effective value of heat flow is three times higher than the average effective heat flow *q*0 in classic TEG, operating in stationary mode.

Consideration of non-stationary mode includes the following parameters: energy capacity, speed of response, response time improvement, temperature dependence of physical parameters and output characteristics, accounting for Peltier heat production, effectiveness of sidesurfaces thermal insulation and thermal stabilization based on the model of finite length. All this information can be very useful for researchers working in the field of transients.

#### **2.3. Thermoelectricity in systems far from equilibrium**

M. Apostol in the year 2000 has presented theory of ultrafast thermoelectric conduction [9]. The theory was developed and described [10–12]. It is about moving heat pulse along a rod of thermoelectric material of length 1 mm with on and off connection along a wire 1 cm long. The electrical power is extracted along this wire. A sketch of the pulses is below (**Figure 2**).

**Figure 2.** The pulse is deflated of heat at the cold end and another pulse is constructed at the hot end. *T* is the temperature, *δT* is the variation in temperature, *t*′ is time off, which is the time between 2 pulses.

The shape of the wave as a function of time *t* and distance *x* is *n*(*x*,*t*) [10]:

$$m(\mathbf{x}, t) = \frac{V \times \delta\_{\mathbf{v}}}{\sqrt{2\pi \times \Lambda \times \nu \times \tau}} \times \exp\left(-\frac{\mathbf{x}^2}{2\Lambda \nu t}\right),\tag{1}$$

here *V* original volume of *δ* peak; *δn* is quasi-particle density in *δ* peak, Λ is mean free path of quasi-particles; *t* is time; *ν* is transport velocity.

The width of the pulse is expressed as:

$$I = \sqrt{\Lambda \times \nu \times \tau\_{\text{off}}} \,, \tag{2}$$

where *l* = width of pulse, *τ*off = time off.

**Figure 1.** Heat flux *q* versus time *t*, is the average effective heat flow during time *t, q*0 is the average effective heat

This graph *q*(*t*) shows that in this model, due to thermal conductivity, average, during time *t*, effective value of heat flow is three times higher than the average effective heat flow *q*0 in

Consideration of non-stationary mode includes the following parameters: energy capacity, speed of response, response time improvement, temperature dependence of physical parameters and output characteristics, accounting for Peltier heat production, effectiveness of sidesurfaces thermal insulation and thermal stabilization based on the model of finite length. All

M. Apostol in the year 2000 has presented theory of ultrafast thermoelectric conduction [9]. The theory was developed and described [10–12]. It is about moving heat pulse along a rod of thermoelectric material of length 1 mm with on and off connection along a wire 1 cm long. The electrical power is extracted along this wire. A sketch of the pulses is below (**Figure 2**).

**Figure 2.** The pulse is deflated of heat at the cold end and another pulse is constructed at the hot end. *T* is the tempera-

ture, *δT* is the variation in temperature, *t*′ is time off, which is the time between 2 pulses.

this information can be very useful for researchers working in the field of transients.

flow in stationary mode, /0 ≈ 3.

classic TEG, operating in stationary mode.

540 Thermoelectrics for Power Generation - A Look at Trends in the Technology

**2.3. Thermoelectricity in systems far from equilibrium**

He gives the relation of maximum power output ext max with pulses of ultrafast conduction compared to maximum power in DC mode dc max as:

$$P\_{\rm ext}^{\rm max} = 2\frac{l}{\sqrt{\Lambda I}} \times P\_{\rm dc}^{\rm max},\tag{3}$$

$$P\_{\text{ext}}^{\text{max}} = 60 \times P\_{\text{dc}}^{\text{max}}.\tag{4}$$

The ratio of 60 is theoretical value calculated with some assumptions. Even when the experimental ratio is much smaller than 60, then the output power ext max can be many times greater than the DC mode value. This can be the major breakthrough that we are searching for to 'beat' material limitation of *ZT* =2.

#### **3. Pulsed mode thermoelectric generators**

#### **3.1. Thermoelectric generator working on MHz frequency**

#### *3.1.1. Primary patents and papers*

Priority on thermoelectric energy conversion in high-frequency mode was filed by Aspen and Strachan in the years 1990–1991 [13]. In the year 1992, Aspen and Strachan published results of development and demonstration tests of MHz thermoelectric generator [7, 8]. This is the first example of thermoelectric device operating far from equilibrium.

#### *3.1.2. Description of thermoelectric device operating far from equilibrium*

Strachan develops a vibrator to break kidney stones. He discovered that his device could operate as a heat pump (to produce ice) or, with temperature difference between two sides, it could generate electrical power sufficient to operate a small fan. Experimental study was done in collaboration with Oxborrow. It was headed by John Stockholm of Marvel Thermoelectrics and was financed by automobile company PSA.

John Stockholm and M. Almeida from Supélec visited Scott Strachan in his office at the Technology Transfer Centre, University of Edinburgh, Scotland, in February 1996. The device built by Scott Strachan over a year ago was no longer operational, as it had degraded with time, but he showed a video and gave us a copy of it.

He was extremely open and gave a great deal of advice on how to make a unit.

A 'still' from the video is shown in **Figure 3**.

**Figure 3.** Scott Strachan operating the device.

Thermoelectric generator is located below his fingers, below the square white box, into which Scott Strachan pours hot water to generate Δ*T* across both sides of TE generator.

It is a small multi-layered device 30 × 2.5 mm (Δ*T* is across dimension 2.5 mm) consisting of stack: Al film—Ni film—piezoelectric film (PVDF)—ethyl cyanoacrylate adhesive—Fe film on BASF recording tape—Mylar, BASF recording tape—ethyl cyanoacrylate adhesive—Al film; repeat periodically (**Figure 4**). A device of several cubic centimetres in volume is vibrated by PZT at 500 kHz. It produces 0.14 W of electrical power from heat throughput of 3.7 W.

**Figure 4.** Strachan device-stack of different materials.

*3.1.2. Description of thermoelectric device operating far from equilibrium*

542 Thermoelectrics for Power Generation - A Look at Trends in the Technology

and was financed by automobile company PSA.

but he showed a video and gave us a copy of it.

A 'still' from the video is shown in **Figure 3**.

**Figure 3.** Scott Strachan operating the device.

Strachan develops a vibrator to break kidney stones. He discovered that his device could operate as a heat pump (to produce ice) or, with temperature difference between two sides, it could generate electrical power sufficient to operate a small fan. Experimental study was done in collaboration with Oxborrow. It was headed by John Stockholm of Marvel Thermoelectrics

John Stockholm and M. Almeida from Supélec visited Scott Strachan in his office at the Technology Transfer Centre, University of Edinburgh, Scotland, in February 1996. The device built by Scott Strachan over a year ago was no longer operational, as it had degraded with time,

Thermoelectric generator is located below his fingers, below the square white box, into which

It is a small multi-layered device 30 × 2.5 mm (Δ*T* is across dimension 2.5 mm) consisting of stack: Al film—Ni film—piezoelectric film (PVDF)—ethyl cyanoacrylate adhesive—Fe film on BASF recording tape—Mylar, BASF recording tape—ethyl cyanoacrylate adhesive—Al film; repeat periodically (**Figure 4**). A device of several cubic centimetres in volume is vibrated by PZT at 500 kHz. It produces 0.14 W of electrical power from heat throughput of 3.7 W.

Scott Strachan pours hot water to generate Δ*T* across both sides of TE generator.

He was extremely open and gave a great deal of advice on how to make a unit.

Energy conversion efficiency calculated from the above values is equal to 4%. This is very high, considering that Δ*T* is below 75 K (boiling water and ambient air), compared to bismuth telluride, which for similar Δ*T* would give much less than 2%.

Two prototypes were made by Mark Oxborrow.

Prototypes are very difficult to make, in particular with 'super glue', if the sheet is badly placed, one has to start the assembly process all over again.

The one, he sent to Supélec, was destroyed during testing and John Stockholm terminated the study, as it was apparent that the device was not adaptable to powers required for automobiles.

Our objective was to test the device, then to make similar more simple devices to determine the key components, but we were not successful in even reproducing the device.

Many people showed an interest in Strachan's device. Details of TE generator were given, in presentation by Maynard and Mahan of Penn State at ICT2005 in Clemson SC [14]. Also, people from MIT showed interest in Strachan's device. But nobody was able to make the device following all advices from Scott Strachan. All those involved have doubts about the theoretical aspects described by Harold Aspen, who by profession was a patent attorney.

It is incredible that Strachan's work could not be reproduced and devices made to determine the key technical parameters. He was pioneer in the realization of far-from-equilibrium regime of TE device by pulsing in MHz range.

#### **3.2. Schroeder's pulsed TE generator**

In the year 1994, Jon Schroeder announced a ring-shaped TE generator pulsed at 60 Hz and he published a paper about it in the year 1999 [15].

A prototype was shown to visitors, but only measurement on one thermoelectric couple was made. This measurement was used to calculate the output electrical power of the ring, but thermoelectric generator did not work.

Only in 2004, TE generator prototype with bismuth telluride elements, pulsed at 200 kHz, was demonstrated. But only the output electrical power of prototype was measured. So, we did not know if pulsed-operated TE generator produced greater efficiency than standard DC-operated TE generator with bismuth telluride couples operated within temperature limits of bismuth telluride.

In 2008, Marin Nedelcu with a friend visited Jon Schroeder in Leander, Texas. A video was made showing Jon Schroeder operating the device.

The electrical conversion efficiency was calculated as exceeding two times the DC mode efficiency. A photograph taken from the video made by Marin Nedelcu is shown in **Figure 5**.

**Figure 5.** Schroeder's pulsed TE generator.

The unit weighs about 12 kg. Heat source is from propane in a bottle and cold source is from Blades facing downwards immersed in water.

Schroeder ring. It consists of a ring designed by Jon Schroeder, described in a patent [16]. The heat is produced in the centre by hot gases; gases can be combustion gases from natural gas. The heat is transferred by convection and radiation to the radial blades 'hot' (**Figure 6**).

**Figure 6.** Schematic of Schroeder's ring. The ring is 250 mm diameter. The cold blades 21 are oriented downwards and are cooled by forced convection. Voltage outlets 36, 37.

**Figure 7** shows cross section, with a central burner; hot blades marked 4 are connected to the hot side of the pieces of TE and cold blades are marked 3. The propane is burned in volume 14 and the hot combustion gases exit through the central hole 20.

**Figure 7.** Cross section of the unit.

Only in 2004, TE generator prototype with bismuth telluride elements, pulsed at 200 kHz, was demonstrated. But only the output electrical power of prototype was measured. So, we did not know if pulsed-operated TE generator produced greater efficiency than standard DC-operated TE generator with bismuth telluride couples operated within temperature limits of bismuth

In 2008, Marin Nedelcu with a friend visited Jon Schroeder in Leander, Texas. A video was

The electrical conversion efficiency was calculated as exceeding two times the DC mode efficiency. A photograph taken from the video made by Marin Nedelcu is shown in **Figure 5**.

The unit weighs about 12 kg. Heat source is from propane in a bottle and cold source is from

Schroeder ring. It consists of a ring designed by Jon Schroeder, described in a patent [16]. The heat is produced in the centre by hot gases; gases can be combustion gases from natural gas. The heat is transferred by convection and radiation to the radial blades 'hot' (**Figure 6**).

made showing Jon Schroeder operating the device.

544 Thermoelectrics for Power Generation - A Look at Trends in the Technology

**Figure 5.** Schroeder's pulsed TE generator.

Blades facing downwards immersed in water.

telluride.

This schematic shows that radial blades 4 are heated by hot gas; blades 3 facing downwards are air-cooled. Thermoelectric material: bismuth telluride (20 × 20 × 1 mm), alternatively n- and p-type, is placed between heated and cooled blades. With temperature difference between hot and cold blades, voltage is created at outlets 36 and 37 (**Figure 6**) that are connected to the primary winding of a transformer through MOSFETs, which operate in 100-kHz range.

Output connections of TE generator are connected to two circuits in parallel (not shown).

The secondary winding of a transformer increases the voltage. The switching between two circuits (primaries wound in different directions of a transformer) is done before opening the circuit, so one primary winding is always connected to TE generator. In this way, electrical current is pulsed at around 200 kHz, then rectified and can be converted into AC current, for example, 220 V.

## **4. Impedance of thermoelectric devices**

We examine values obtained from three different approaches of the impedance:




#### **4.1. Onsager equations**

Onsager equations are for stationary systems; equations can be extended to non-stationary conditions as presented by Goupil [5].

Assuming that the entropy per charge carrier expression is:

$$A = -e \times \alpha,\tag{5}$$

α is Seebeck coefficient, *e* is charge of the electron.

The others terms derive directly from the thermo-elastic coefficients of the electronic gas,

0 <sup>=</sup> 2 , = , = 2 , where *K* and *G* are, respectively, thermal and electrical

conductance of the system. We also have the specific heat:

$$C\_N = \frac{T}{N} \times \left(\frac{\partial \mathbf{S}}{\partial T}\right)\_N,\tag{6}$$

is the isothermal compressibility, *S* is the entropy:

Non-Stationary Thermoelectric Generators http://dx.doi.org/10.5772/66421 547

$$\mathcal{Z}\_T = \frac{1}{N} \times \left(\frac{\partial N}{\partial \mu}\right)\_S,\tag{7}$$

with *N* the number of carriers and *μ* the electrochemical potential. After calculation of the eigenvalues and eigenvectors, we get the complete expressions of the transient response of the system. In most of the cases, the thermal and electrical time constants, respectively, thermal and elec, are mixed together in a complicated way, but two extreme cases can be considered:

– Under very weak coupling, < < /, we find classical results, elec = 0 and:

$$
\tau\_{\text{thermal}} = \tau\_0 \frac{l}{L};
\tag{8}
$$

– Under strong coupling condition, > > /, thermal time constant diverges and electrical time constant becomes:

$$
\sigma\_{\rm elec} = \frac{C\_N}{\chi\_T e^2 \alpha^2 T} \sigma\_0 = \frac{C\_N}{\alpha^2 G T}. \tag{9}
$$

Since *G* is the resistive part of the impedance, the reactive part gives the expression of the socalled 'thermoelectric capacitance':

$$C\_{\rm TE} = \frac{C\_N}{a^2 T}.\tag{10}$$

Thermoelectric time constant of the system is then finally:

Output connections of TE generator are connected to two circuits in parallel (not shown).

We examine values obtained from three different approaches of the impedance:

thermoelectric modules was proposed by García-Cañadas [17].


Assuming that the entropy per charge carrier expression is:

2

conductance of the system. We also have the specific heat:

is the isothermal compressibility, *S* is the entropy:

α is Seebeck coefficient, *e* is charge of the electron.

, =

example, 220 V.

flow.

0 <sup>=</sup> 2 , =

**4.1. Onsager equations**

conditions as presented by Goupil [5].

 

**4. Impedance of thermoelectric devices**

546 Thermoelectrics for Power Generation - A Look at Trends in the Technology

The secondary winding of a transformer increases the voltage. The switching between two circuits (primaries wound in different directions of a transformer) is done before opening the circuit, so one primary winding is always connected to TE generator. In this way, electrical current is pulsed at around 200 kHz, then rectified and can be converted into AC current, for



Onsager equations are for stationary systems; equations can be extended to non-stationary

a

The others terms derive directly from the thermo-elastic coefficients of the electronic gas,

, æ ö ¶ = ´ ç ÷ è ø ¶ *<sup>N</sup>*

*T S <sup>C</sup>*

*N*

, (5)

, where *K* and *G* are, respectively, thermal and electrical

*N T* (6)

*A e* = ´-

$$
\pi\_{\rm TE} = R\_{\rm re} C\_{\rm TF} \,\tag{11}
$$

Reference [5] gives all the intermediate steps and the result is thermoelectric capacitance:

$$C\_{\rm TE} = \frac{\Im N \mathbf{k}}{a^2 T}.\tag{12}$$

Assuming standard sample with parameters close to room temperature, (*V′* = 1 cm3 , *n* = 5 × 1019 cm-3, *T* = 400 K, *α* = 160 µV K−1), we get:

$$C\_{\rm TE} = \frac{3Nk}{a^2T} \approx 202\text{ F.}\tag{13}$$

Complete internal resistance of the system is now replaced by internal impedance:

$$\mathbf{R}\_{\rm in} = \mathbf{R}\_{\rm TE} \;/\ \mathbf{C}\_{\rm TE} + \mathbf{R}.\tag{14}$$

This result shows that the internal resistance is larger at zero frequency (stationary conditions), // means that RTE and CTE are electrically in parallel, (f=0) = RTE +R , than it is for high frequencies, in( = 0) = TE + , then it is for high frequencies, since *R*TE term becomes shortcircuited by *C*TEG giving in( > > 0) = .

TE <sup>≈</sup> ≈ <sup>10</sup>−4 Ohm; TE ≈ TETE ≈ 2 × 10−1 s, which gives low-pass filtering at really low frequencies. Depending on precise values of internal isothermal resistance, low-pass filtering effect may be found in the range from 1 kHz to 1 MHz.

#### **4.2. Spectroscopic impedance**

Spectroscopy analysis of thermoelectric devices has been performed by García-Cañadas [17] and Gao Min [18]. They have studied thermodynamics of thermoelectric phenomena by frequency-resolved methods, which show that TE devices have a capacitive component that can be exploited to improve performance.

#### **4.3. Basics of pulsed-operating thermoelements**

Apostol and Cune [11] show how electrical circuit of thermoelectric device can be represented with electrical resistance in series and capacitor in parallel.

## **5. Experimental devices**

#### **5.1. Extraction of electrical power**

Electrical power from DC-operated TE systems is most commonly extracted on electrical load resistance. Recently, pulse-width modulation (PWM) with matched load has been used, especially, when powers and voltages are small. In the case of pulsed systems, we have already AC system, so PWM is the logical way to extract electrical power.

When we have DC-operated thermoelectric generator, then power is extracted on resistive load. The maximum power is extracted, when external load resistance *RL* is equal to the resistance of power source *RS*. It was shown in 2001 by M. Nedelcu and J. Stockholm that when electrical current was pulsed at 50 kHz, then the electrical output power was constant at < . This observation showed that small electrical resistance, such as impedance, might be of interest to extract electrical power efficiently.

Equivalent electrical circuit contains electrical resistance of thermoelectric material, capacitance in parallel, to which one can add an inductance composed of the primary winding of a transformer to extract electrical power at voltages of usable values. Today, there are no mathematical models to dimension such circuits.

Spectroscopic impedance measurements made on commercial thermoelectric modules (e.g., 40 × 40 mm with 127 couples) have shown very high values of capacitance in Farad range. We are dealing with super-capacitors, which can be operated at more than 100 kHz. Several teams are designing and making different units to check ultra-fast conduction theory.

## **5.2. Devices**

in TE TE *ZR C R* = + // . (14)

This result shows that the internal resistance is larger at zero frequency (stationary conditions), // means that RTE and CTE are electrically in parallel, (f=0) = RTE +R , than it is for high frequencies, in( = 0) = TE + , then it is for high frequencies, since *R*TE term becomes short-

TE <sup>≈</sup> ≈ <sup>10</sup>−4 Ohm; TE ≈ TETE ≈ 2 × 10−1 s, which gives low-pass filtering at really low frequencies. Depending on precise values of internal isothermal resistance, low-pass filtering

Spectroscopy analysis of thermoelectric devices has been performed by García-Cañadas [17] and Gao Min [18]. They have studied thermodynamics of thermoelectric phenomena by frequency-resolved methods, which show that TE devices have a capacitive component that

Apostol and Cune [11] show how electrical circuit of thermoelectric device can be represented

Electrical power from DC-operated TE systems is most commonly extracted on electrical load resistance. Recently, pulse-width modulation (PWM) with matched load has been used, especially, when powers and voltages are small. In the case of pulsed systems, we have already

When we have DC-operated thermoelectric generator, then power is extracted on resistive load. The maximum power is extracted, when external load resistance *RL* is equal to the resistance of power source *RS*. It was shown in 2001 by M. Nedelcu and J. Stockholm that when electrical current was pulsed at 50 kHz, then the electrical output power was constant at < . This observation showed that small electrical resistance, such as impedance, might

Equivalent electrical circuit contains electrical resistance of thermoelectric material, capacitance in parallel, to which one can add an inductance composed of the primary winding of a transformer to extract electrical power at voltages of usable values. Today, there are no

circuited by *C*TEG giving in( > > 0) = .

can be exploited to improve performance.

**4.3. Basics of pulsed-operating thermoelements**

with electrical resistance in series and capacitor in parallel.

AC system, so PWM is the logical way to extract electrical power.

be of interest to extract electrical power efficiently.

mathematical models to dimension such circuits.

**4.2. Spectroscopic impedance**

**5. Experimental devices**

**5.1. Extraction of electrical power**

effect may be found in the range from 1 kHz to 1 MHz.

548 Thermoelectrics for Power Generation - A Look at Trends in the Technology

Jon Schroeder built the first high-frequency-operated generator, but, unfortunately, he has not been able to make a second unit. There are explanations. The ring design is the ideal design, but it has major flaws: when a generator is operated and stopped, the ring expands and contracts creating shear stresses on a thermoelectric material; these stresses degrade bismuth telluride material and interfaces with copper. The second flaw is that all thermoelectric discs must be soldered in one operation; chances of success are very small. So, prototypes being built now have column structure, where TE elements are in line and compressed between heat sources, which are placed laterally or in line, when liquid is used.

#### *5.2.1. Column structure*

Several teams are working on devices based on the column structure.

#### *5.2.1.1. Column design by Nedelcu*

A team headed by Marin Nedelcu in Bucharest started by building ring units with individual water cooling for each bismuth telluride disc, but following problems of water leaks abandoned the ring and resorted to column structure. M. Nedelcu used in column design bismuth telluride discs: *D* = 24 mm and 1 mm thick; he built several columns, the latest is shown in **Figure 8**.

**Figure 8.** Thermoelectric column built by Nedelcu.

Nedelcu has built eight columns the same as this one. Each column has seven TE couples plus one at the end, so that both ends are cold. Approximate dimensions of cold blocks are 30 × 35 × 8 mm with two holes with a diameter of 4 mm. Hot plate is attached to the top flat surfaces; water flows through the holes through the cold side, which are connected by rubber u-bends, not shown. The water is distilled and has high electrical resistance. The assembly is on-going and testing is planned.

## *5.2.1.1.1. Commercial module design*

Recently, Nedelcu built TE generator using four commercial electricity generating modules connected electrically in parallel, with total electrical resistance around 30 mOhm. Modules are assembled between water-cooled plate at *T* ~38°C and heated plate at *T* ~188°C, heating power with 170 V and 5.9 A = 1003 W.

Electrical current output was pulsed using MOSFETs at around 200 kHz. The output current from the transformer is rectified. The load is 100 W filament electrical light bulb. Measured voltage is 210 V and current is 0.4 A, so electrical power output is 84 W. So, overall efficiency (including heat losses) is 8.4%; this is about two times higher conversion efficiency, when operated in DC mode.

These measurements seem to confirm that pulsing can improve efficiency due to ultra-fast conduction as predicted by Apostol. When dealing with such breakthrough results, one must be extremely cautious and obtain results from another laboratory.

This TE device made by Nedelcu is very different from Schroeder ring (120 slices of bismuth telluride with an area of 4 cm2 , thickness of 1 mm, all in series electrically); here we have about 71 couples of bismuth telluride (unknown dimensions, can be estimated: area 10 mm2 ).

#### *5.2.1.2. Pulsed TE generators by Krzysztof Wojciechowski team*

A team headed by Krzysztof Wojciechowski of AGH Cracow and Institute of Electron Technology, Cracow, designed and built several small prototypes to detect high-conversion efficiencies by pulsing thermoelectric circuit. First prototype had water-cooled plate at the top and the bottom, see **Figure 9**.

**Figure 9.** Unit with 18 bismuth telluride discs (n- and p-type).

**Figure 10** shows a smaller unit with six TE discs with a diameter of 24 mm and a thickness of 1 mm, aluminium water-cooled plates above and below, electrical heaters at the front and back. Detail of column electrical heater 'red' at top and bottom without the water cooled plate, that would be located on the right side and on the left side (**Figure 11**).

**Figure 10.** Unit with six bismuth telluride discs—electrical heaters front and back, compact water coolers, see **Figure 12**.

**Figure 11.** Unit with three TE discs of an outside diameter 24 mm and thickness 2 mm. The heaters are electrical strips on the left and right sides. The cooling is with compact cooler as shown in **Figure 12**.

Thermal and electrical measurements showed that produced voltage was insufficient to operate MOSFET switching. Once again, this confirms the difficulty to have sufficient heat flux through bismuth telluride material, as bismuth telluride has thermal conductivity of about 1.5 W/(m K), si in bismuth telluride disc with a cross section of 4.5 cm2 and a thickness of 1 mm, heat flux power at Δ*T* = 100 K is equal to 67.5 W (**Figure 12**).

#### *5.2.1.3. Unit with column structure by Glick et al*

ber u-bends, not shown. The water is distilled and has high electrical resistance. The assem-

Recently, Nedelcu built TE generator using four commercial electricity generating modules connected electrically in parallel, with total electrical resistance around 30 mOhm. Modules are assembled between water-cooled plate at *T* ~38°C and heated plate at *T* ~188°C, heating

Electrical current output was pulsed using MOSFETs at around 200 kHz. The output current from the transformer is rectified. The load is 100 W filament electrical light bulb. Measured voltage is 210 V and current is 0.4 A, so electrical power output is 84 W. So, overall efficiency (including heat losses) is 8.4%; this is about two times higher conversion efficiency, when

These measurements seem to confirm that pulsing can improve efficiency due to ultra-fast conduction as predicted by Apostol. When dealing with such breakthrough results, one must

This TE device made by Nedelcu is very different from Schroeder ring (120 slices of bismuth

A team headed by Krzysztof Wojciechowski of AGH Cracow and Institute of Electron Technology, Cracow, designed and built several small prototypes to detect high-conversion efficiencies by pulsing thermoelectric circuit. First prototype had water-cooled plate at the top

**Figure 10** shows a smaller unit with six TE discs with a diameter of 24 mm and a thickness of 1 mm, aluminium water-cooled plates above and below, electrical heaters at the front and back. Detail of column electrical heater 'red' at top and bottom without the water cooled plate, that

71 couples of bismuth telluride (unknown dimensions, can be estimated: area 10 mm2

, thickness of 1 mm, all in series electrically); here we have about

).

be extremely cautious and obtain results from another laboratory.

*5.2.1.2. Pulsed TE generators by Krzysztof Wojciechowski team*

**Figure 9.** Unit with 18 bismuth telluride discs (n- and p-type).

would be located on the right side and on the left side (**Figure 11**).

bly is on-going and testing is planned.

550 Thermoelectrics for Power Generation - A Look at Trends in the Technology

power with 170 V and 5.9 A = 1003 W.

operated in DC mode.

telluride with an area of 4 cm2

and the bottom, see **Figure 9**.

*5.2.1.1.1. Commercial module design*

A team with John Stockholm, James Glick and Brad Vier are also building a TE generator unit, which has a column structure. First prototype P1 was built with five couples of bismuth telluride discs with a diameter of 24 mm and a thickness of 1 mm. It was operated with propane combustion gases and cooled by ambient air. The objective was to measure temperatures at the interface of bismuth telluride.

A schematic of pulsed-width modulation electrical circuit is shown in **Figure 13**.

**Figure 12.** Detail of compact water cooler with an outside diameter of 25 mm placed on the axis of the column.

**Figure 13.** Schematic of PWM controller electrical circuit.

A photograph in **Figure 14** shows a lateral view of the inside of TE generator.

**Figure 14.** Lateral view of inside of P1 prototype TE generator, left side at the bottom is the infrared burner, combustion gases go upwards; in the centre at top, column structure-heated blades; to the left, cooled blades by ambient air blown upwards (fan is not shown).

**Figure 15** shows the top view of TE generator: on the left are heated blades and on the right are blades cooled by air, the white material is to insulate thermally thermoelectric column from hot combustion air. Bismuth telluride disc with copper caps is shown in **Figure 16**.

**Figure 15.** Top view of inside of P1 prototype.

**Figure 12.** Detail of compact water cooler with an outside diameter of 25 mm placed on the axis of the column.

A photograph in **Figure 14** shows a lateral view of the inside of TE generator.

**Figure 14.** Lateral view of inside of P1 prototype TE generator, left side at the bottom is the infrared burner, combustion gases go upwards; in the centre at top, column structure-heated blades; to the left, cooled blades by ambient air

**Figure 13.** Schematic of PWM controller electrical circuit.

552 Thermoelectrics for Power Generation - A Look at Trends in the Technology

blown upwards (fan is not shown).

**Figure 16.** Bismuth telluride disc with copper caps.

When the unit was operated, the propane combustion gases heated and forced ambient air cooled. However, measured temperature difference at the interface with bismuth telluride was insufficient to operate TE generation, but valuable information was obtained to improve the design.

Stockholm presented in 2014 results concerning different aspects of pulsed TE generator operation, such as high-efficiency pulse-width modulator thermoelectric generator with an active load [19] and transient thermoelectric generator: an active load story [20].

Development is on-going. A prototype designed for measurements is being made.

**Figure 17** shows the 10 kW transformer to be installed. **Figure 18** shows four Hi-Z-14 modules.

Results will be published as soon as they have been repeated to ensure validity.

**Figure 17.** PWM 10 kW transformer dimensions 250 × 190 × 140 mm high.

**Figure 18.** Four Hi-Z-14 modules on a water-cooled plate.

## **6. Potential of non-stationary far-from-equilibrium TE generators**

Very few people have studied the microscopic phenomena of non-stationary, far-fromequilibrium thermoelectrics. Marian Apostol is the main contributor to the subject with many papers. He has shown great potential of this phenomenon. Why is it that 15 years after his first paper nobody has an industrial unit? There are many reasons. People don't believe new phenomena!

Technology of these systems is difficult, because of constraints: high heat fluxes and very low electrical resistances. Electrical currents can be enormous, more than 1000 A and voltages can reach kV. All this is confirmed by the numerous problems encountered by all experimentalists. Strachan had great voltage problems. High-frequency pulsing is not stable. Marian Apostol compares these systems to a sphere on top of a mountain! Prototypes are being built and tested to confirm these high-energy conversion efficiencies.

In 2014, measurements on the unit built by Schroeder were made in his workshop, on a unit not designed for any measurements inside the unit. The only possible measurements were the amount of propane burnt and electrical power output. Measured values of conversion efficiencies exceeded by a factor of more than 2, the values, of equivalent DC-operated unit with similar temperature differences.

Recently, Nedelcu has built a TE device with four commercial electricity-generating modules Hi-Z14 (56 × 56 mm) with 31 bismuth telluride couples. They are electrically in parallel so that the electrical resistance of the load is minimum about 40 mOhm placed between heated plate at 188°C and cold plate at 38°C. He claims that the conversion efficiency exceeds two times the DC conversion efficiency from thermal to electrical power. All these claims need to be repeated in a laboratory and fully documented in particular measurements of electrical pulses and power on primary circuit must be made.

There are also potential sources of degradation at interfaces. Anatol Casian at the University of Moldova is concerned about possible charge carriers' depletion.

Today, several research groups, on very small budgets, are working on devices to confirm experimentally that very high-performing generators are a reality; as soon as there is confirmation, research money will be available. Many theoretical researches are necessary before any transfer to industry.

## **7. Conclusions and outlook**

insufficient to operate TE generation, but valuable information was obtained to improve the

Stockholm presented in 2014 results concerning different aspects of pulsed TE generator operation, such as high-efficiency pulse-width modulator thermoelectric generator with an

**Figure 17** shows the 10 kW transformer to be installed. **Figure 18** shows four Hi-Z-14 modules.

active load [19] and transient thermoelectric generator: an active load story [20].

554 Thermoelectrics for Power Generation - A Look at Trends in the Technology

Results will be published as soon as they have been repeated to ensure validity.

**Figure 17.** PWM 10 kW transformer dimensions 250 × 190 × 140 mm high.

**Figure 18.** Four Hi-Z-14 modules on a water-cooled plate.

Development is on-going. A prototype designed for measurements is being made.

design.

The operation in non-stationary modes near equilibrium has been presented. Two books, one by Gray (1960) [2] and one by Iordanishvili [4], cover extensively the subject. Recently, Nedelcu et al. [21], Min [18] and García-Cañadas [17] have done experimental spectroscopic analysis of thermoelectric systems, which confirm that TE device is equivalent electrically with a capacitance in parallel with TE generator, which has internal electrical resistance. This is confirmed by theoretical papers by Apostol [9–11] and Goupil [5].

One important advantage of operating in a non-stationary mode is to obtain thermoelectric generators with better performances, than with DC operation. The present-day DC-operated TE generators are limited by *ZT* of TE material, which is less than 2.

Measurements need to be confirmed by other laboratories. The main interest is energy conversion efficiency. Prototypes have mainly used bismuth telluride, because it is the most suitable for generators operating with the cold side around ambient, but any thermoelectric material could be used.

Prototypes are being built and tested to confirm these high-energy conversion efficiencies, which are predicted by ultra-fast conduction theory of M. Apostol [9].

## **Acknowledgements**

The author wishes to thank all the people who have contributed to his interest in non-stationary TE systems, especially Nikolai Kukhar, Valeriy Zhmurko, Lyudmyla Vikhor, Oleg Luste, Scott Strachan, Jon Schroeder, Marian Apsotol, Marin Nedelcu and Christophe Goupil.

## **Author details**

John G. Stockholm

Address all correspondence to: johngstockholm@gmail.com

Marvel Thermoelectrics, Joachim du Bellay, Vernouillet, France

## **References**


tance in parallel with TE generator, which has internal electrical resistance. This is confirmed

One important advantage of operating in a non-stationary mode is to obtain thermoelectric generators with better performances, than with DC operation. The present-day DC-operated

Measurements need to be confirmed by other laboratories. The main interest is energy conversion efficiency. Prototypes have mainly used bismuth telluride, because it is the most suitable for generators operating with the cold side around ambient, but any thermoelectric

Prototypes are being built and tested to confirm these high-energy conversion efficiencies,

The author wishes to thank all the people who have contributed to his interest in non-stationary TE systems, especially Nikolai Kukhar, Valeriy Zhmurko, Lyudmyla Vikhor, Oleg Luste, Scott

[1] Stilbans LS, Fedorov NA. On the work of cooling thermoelements in the transient

[2] Gray P. The dynamic behavior of thermoelectric devices. The Technology Press, MIT &

[3] Landecker K, Findlay AW. Study of the fast transient behaviour of Peltier junctions.

regime. Journal Tekhn. Physics. 1958;XXVIII:489–492 (in Russian).

John Wiley; 1960 (Library of Congress No. 60-12981).

Solid-State Electronics. 1961;3:239–260.

Strachan, Jon Schroeder, Marian Apsotol, Marin Nedelcu and Christophe Goupil.

by theoretical papers by Apostol [9–11] and Goupil [5].

556 Thermoelectrics for Power Generation - A Look at Trends in the Technology

material could be used.

**Acknowledgements**

**Author details**

John G. Stockholm

**References**

TE generators are limited by *ZT* of TE material, which is less than 2.

which are predicted by ultra-fast conduction theory of M. Apostol [9].

Address all correspondence to: johngstockholm@gmail.com

Marvel Thermoelectrics, Joachim du Bellay, Vernouillet, France


[19] Stockholm J, Wojciechowski K, Leszcznski J, Witek K, Guzdec P, Kulawik J. Thermoelectric generator with an active load. In: Materials Today Proceedings 12th European

[20] Stockholm J, Goupil C, Maussion P, Ouerdan H. Transient thermoelectric generator: an active load story. Journal of Electronic Materials. Proceedings of International Confer-

[21] Nedelcu M, Apostol J, Stockholm JG. Pulsed thermoelectric machine. Materialy

Conference on Thermoelectrics, Madrid; 2014.

558 Thermoelectrics for Power Generation - A Look at Trends in the Technology

ence on Thermoelectrics; Nashville, TE 2014.

Ceramiczne/Ceramic Materials. 2014;66:170–174.
