Abstract

The present work shows a comprehensive methodology and design steps to recover energy from the automotive waste heat. A thermoelectric generator must be connected to a power converter in order to extract the maximum power from the generator and, also, satisfy different constrains to charge a battery. Starting from the electrical model of thermoelectric cells, it is evaluated their combination to realize a thermoelectric generator (TEG) comply with the automotive regulation, then considering input/output electric characteristics, it is evaluated the best converter topology to satisfy all constrains. Design steps and power dissipation estimation are deeply explained. TEG and power converter models are simulated in a model-based environment to allow the design of the control algorithms. The control system consists of nested control loops. Two maximum power point tracking (MPPT) algorithms are evaluated. The MPPT output is used as reference for a current control loop. The maximum power characteristic of a TEG has a quadratic behavior and working without the tracking of the maximum power point could drastically decrease the generated power from the TEG and the system efficiency. There are presented simulation results of the control algorithms and experimental data are shown in order to validate the design steps.

Keywords: automotive, thermoelectric generator, waste heat, storage, DC-DC converter, MPPT

## 1. Introduction

Nowadays the reuse of wasted energy is an important task for many fields, particularly for the automotive one. This aims at more efficient systems, reducing pollution and making the product more compatible with the environment, regulations, and consumers. Internal combustion engine vehicles do not totally convert the chemical energy into mechanical, the greater part of a typical fuel path is wasted in the coolant system and in the exhaust, and according to Chuang and Chau [1], only 25% of the energy flows into specific vehicle operations (Figure 1). The strategy is to convert the heat wasted energy into electricity by means of a thermoelectric generator (TEG), which is based on the Seebeck effect. Heat flows through

**108**

*Advanced Thermoelectric Materials for Energy Harvesting Applications*

[9] Nikkiso Co. Standard Specifications for UV LED, model VPS131, Vol. 1.2;

[10] Kneissl M et al. Development of UV-C LED Disifection. TECHNEAU;

[11] Akca H, Yasa Y, Ayaz R, Durusu A, Ajder A, Nakir I, Tanrioven M. Thermal management of power LED system. In: Proceedings of ICRERA; Milwaukee,

[12] Xi Y, Gessmann T, Shah JM, Kim JK, Schuberta EF. Junction and carrier temperature measurements in deepultraviolet light-emitting diodes using three different methods. Applied Physics Letters. 2005;**86**:031907

[13] Kneissl M, Rass J. A brief review of III-nitride UV emitter technologies and their applications. In: III-Nitride Ultraviolet Emitters. 1st ed. Vol. 227, Springer Series in Materials Science.

[14] Reed ML, Wraback M, Lunev A, Bilenko Y, Hu X, Sattu A, et al. Device self-heating effects in deep UV LEDs studied by systematic variation in pulsed current injection. Physica Status

[15] Riffat SB, Ma X. Thermoelectrics: A review of present and potential applications. Applied Thermal Engineering. 2003;**23**(8):913-935

Solidi. 2008;**5**(6):2053-2055

2016. pp. 1-25

WI, USA; 2014. pp. 760-764

2016

2010

[1] Khan A, Balakrishnan K, Katona T. Ultraviolet light-emitting diodes based on group three nitrides. Nature

[2] Sholtes KA, Lowe K, Walters GW, Sobsey MD, Linden KG, Casanova LM. Comparison of ultraviolet lightemitting diodes and low-pressure mercury-arc lamps for disinfection of water. Journal of Enviromental Technology. 2016;**37**:2183-2188

[3] Hirayama H, Maeda N, Fujikawa S, Toyoda S, Kamata N. Recent progress and future prospects of AlGaN-based high-efficiency deep-ultraviolet lightemitting diodes. Japanese Journal of Applied Physics. 2014;**53**:100209

[4] Inoue S, Tamari N, Taniguchi M. 150 mW deep-ultraviolet lightemitting diodes with large-area AlN nanophotonic light-extraction structure emitting at 265nm. Applied Physics

Letters. 2017;**110**:141106

[5] Kheyrandish A, Mohseni M,

Research. 2017;**122**:570-579

[6] Schur MS, Gaska R. Deep-

2010;**57**:12-25

Taghipour F. Development of a method for the characterization and operation of UV-LED for water treatment. Water

ultraviolet light-emitting diodes. IEEE Transactions on Electron Devices.

[7] Verszellesi G, Saguatti D, Meneghini M, Bertazzi F, Goano M, Meneghesso G, et al. Efficiency droop in InGaN/ GaN blue light-emitting diodes: Physical mechanisms and remedies. Journal of Applied Physics. 2013;**114**:071101

[8] Pernot C. Improving Deep UV LED Performance. Presented at the LED China Conference, Shanghai, China; 2016. http://www.semi.org.cn/semiconchina/ mail/160331/1/02-SEMICONCHINA-

DUV-LED-cyrilpernot.pdf

Photonics. 2008;**2**:77-84

#### Advanced Thermoelectric Materials for Energy Harvesting Applications

the TEG that, thanks to the difference of temperature (ΔT) between the hot and cold sides, generates electrical power that flows into the connected load. Figure 2 shows the quantities that interact with the TEG: heat, temperature, thermal resistance, interface material, and electrical quantities.

into the converter in order for each delta temperature to be the correct load. The converter output can be connected to the vehicle battery, so the generated power can be directly used for vehicle operations or stored for future needs. The described system allows to recover part of the wasted energy increasing the total efficiency of the vehicle. Figure 3 shows a conceptual diagram of the system from fuel combustion to wasted heat and the recovery system, TEG plus DC-DC, that generates

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

This chapter describes a comprehensive methodology for the design of an energy recovery system. Starting from physical characteristics of the generator, its electrical model is obtained. Since the generator under examination is a thermoelectric generator, all derived parameters will be expressed as a function of the delta temperature, the temperature difference between the hot and cold sides of the generator. Firstly, it is characterized and modeled a single thermoelectric cell, then it is designed the cells combination to make the TEG generator: the resulting equivalent electrical model of the generator is similar to a single cell with grater voltage and current parameters. The parameters obtained by the generator model are used as constraint for the DC-DC converter topology definition, and the design steps are described, highlighting the necessary trade-offs. After these stages it is possible to face the control algorithms. In order to design a control algorithm, there are two necessary information: what the target is and what the weapons are to reach the target. Finally, the single sections are put together in a single simulation environment where it is possible to test the global system changing temperature and load conditions.

Figure 4 shows the open-circuit voltage and the maximum power characteristics of the hypothetical thermoelectric cell. They are represented as a function of the delta temperature across the cell. The maximum allowed delta temperature is 180° C. The open-circuit voltage (VOC) grows linearly with the delta temperature (Eq. (1)), while the available maximum power (PMAX) has a parabolic behavior (Eq. (2)). The maximum open-circuit output voltage is almost 8 V, and the

TEG cell characteristics as a function of delta temperature; on the left the open-circuit voltage and on the right

power useful for vehicle operations.

DOI: http://dx.doi.org/10.5772/intechopen.86232

2.1 TEG and its electrical model

Figure 4.

111

the available maximum power.

2. Design methodology of an energy recovery system

The increase of delta temperature also increases the available power, but for each delta temperature, there is only one maximum power that can be generated connecting a precise load. A simple approach is to connect a DC-DC converter after the TEG. The DC-DC converter can set the current flow; thus, from the TEG point of view, it acts as a variable load. Specific control algorithms can be implemented

Figure 1.

Quantity of energy used for the vehicle operation.

#### Figure 2.

TEG quantities. Qin is the heat quantity to flow into TEG, Qout is the heat quantity to flow out of TEG,Tq is the temperature of heat source,Ta is the temperature of ambient (cold source),Th is the temperature of TEG hot-side surface,Tc is the temperature of TEG cold-side surface,TIM is the thermal interface material, Rh is the heat resistance of TIM between TEG and heat source, Rc is the heat resistance of TIM between TEG and ambient, P is the electrical power, V is the voltage, I is the current, and RL is the load resistance.

Figure 3. Diagram of the wasted heat energy recovery system.

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

into the converter in order for each delta temperature to be the correct load. The converter output can be connected to the vehicle battery, so the generated power can be directly used for vehicle operations or stored for future needs. The described system allows to recover part of the wasted energy increasing the total efficiency of the vehicle. Figure 3 shows a conceptual diagram of the system from fuel combustion to wasted heat and the recovery system, TEG plus DC-DC, that generates power useful for vehicle operations.

#### 2. Design methodology of an energy recovery system

This chapter describes a comprehensive methodology for the design of an energy recovery system. Starting from physical characteristics of the generator, its electrical model is obtained. Since the generator under examination is a thermoelectric generator, all derived parameters will be expressed as a function of the delta temperature, the temperature difference between the hot and cold sides of the generator. Firstly, it is characterized and modeled a single thermoelectric cell, then it is designed the cells combination to make the TEG generator: the resulting equivalent electrical model of the generator is similar to a single cell with grater voltage and current parameters. The parameters obtained by the generator model are used as constraint for the DC-DC converter topology definition, and the design steps are described, highlighting the necessary trade-offs. After these stages it is possible to face the control algorithms. In order to design a control algorithm, there are two necessary information: what the target is and what the weapons are to reach the target. Finally, the single sections are put together in a single simulation environment where it is possible to test the global system changing temperature and load conditions.

#### 2.1 TEG and its electrical model

Figure 4 shows the open-circuit voltage and the maximum power characteristics of the hypothetical thermoelectric cell. They are represented as a function of the delta temperature across the cell. The maximum allowed delta temperature is 180° C. The open-circuit voltage (VOC) grows linearly with the delta temperature (Eq. (1)), while the available maximum power (PMAX) has a parabolic behavior (Eq. (2)). The maximum open-circuit output voltage is almost 8 V, and the

#### Figure 4.

TEG cell characteristics as a function of delta temperature; on the left the open-circuit voltage and on the right the available maximum power.

the TEG that, thanks to the difference of temperature (ΔT) between the hot and cold sides, generates electrical power that flows into the connected load. Figure 2 shows the quantities that interact with the TEG: heat, temperature, thermal resis-

The increase of delta temperature also increases the available power, but for each delta temperature, there is only one maximum power that can be generated connecting a precise load. A simple approach is to connect a DC-DC converter after the TEG. The DC-DC converter can set the current flow; thus, from the TEG point of view, it acts as a variable load. Specific control algorithms can be implemented

TEG quantities. Qin is the heat quantity to flow into TEG, Qout is the heat quantity to flow out of TEG,Tq is the temperature of heat source,Ta is the temperature of ambient (cold source),Th is the temperature of TEG hot-side surface,Tc is the temperature of TEG cold-side surface,TIM is the thermal interface material, Rh is the heat resistance of TIM between TEG and heat source, Rc is the heat resistance of TIM between TEG and ambient, P is

the electrical power, V is the voltage, I is the current, and RL is the load resistance.

tance, interface material, and electrical quantities.

Advanced Thermoelectric Materials for Energy Harvesting Applications

Figure 1.

Figure 2.

Figure 3.

110

Diagram of the wasted heat energy recovery system.

Quantity of energy used for the vehicle operation.

maximum power is a little more than 10 W. The coefficients mv and mp are specifics of the selected TEG cell:

$$V\_{OC} = m\_v \cdot \Delta T \tag{1}$$

$$P\_{MAX} = m\_p \cdot \Delta T^{\prime} \tag{2}$$

Figure 5 shows the voltage-current and the power-voltage characteristics at a single delta temperature. When the TEG output is closed in short circuit, the voltage is zero, and the output current reach is local maximum (Isc). It is important, with the purpose of generating the available maximum power, to know, for each delta temperature, the correct voltage-current couple. In Eq. (3) the relation between the voltage to generate the maximum power (Vmpp) and the open-circuit voltage is pointed out, namely, the delta temperature. In Eq. (4) the current to generate the maximum power (Impp) is derived as a function of the already calculated parameters:

$$V\_{mpp} = \frac{V\_{OC}}{2} = \frac{m\_v}{2} \cdot \Delta T \tag{3}$$

$$I\_{mpp} = \frac{P\_{MAX}}{V\_{mpp}} = \frac{2 \cdot m\_p}{m\_v} \cdot \Delta T \tag{4}$$

From the voltage-current characteristic, it is also possible to extract the value of the TEG equivalent series resistance (Eq. (5)). Figure 6 shows the voltage-current characteristic at different delta temperature: different slopes mean different equivalent resistance:

$$R\_{eq} = \left| \frac{I\_{sc}}{-V\_{oc}} \right| \tag{5}$$

generator with a variable series resistance, both proportional of the delta

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

V0

R0

P0

eq ¼ Req �

max ¼ V<sup>0</sup> oc <sup>2</sup>

Cs Cp

=R<sup>0</sup>

More TEG cells can be combined in order to increase the total generated power. Combination can be done by putting the cells in series, increasing output voltage, or, in parallel, increasing the output current. Due to regulatory constraint, it is better to keep the maximum voltage under 60 V. Furthermore, the TEG pack design is strictly related to available space, point of contact of the exhaust pipe, and heat cooling passing through the TEG (it will cause different available powers between near cells) [2, 3]. The TEG pack for this application is constituted by six cells in series and four in parallel. Figure 8 shows the equivalent electric model of the TEG pack constituted by 24 cells. Under the hypothesis of a good pack design, Eqs. (7)– (12) give a general relation between electrical parameters and the number of cells in

Req ¼ mr � ΔT þ qr (6)

oc ¼ VOC � Cs (7)

eq (9)

(8)

temperature:

Figure 7.

113

Figure 6.

Equivalent resistance as a function of temperature difference.

DOI: http://dx.doi.org/10.5772/intechopen.86232

TEG equivalent resistance as a function of delta temperature.

series (Cs) and in parallel (Cp):

Figure 7 shows the equivalent series resistance calculated at different delta temperature. Eq. (6) shows the equivalent series resistance as a function of delta temperature. The equivalent electric model of a TEG cell is a variable voltage

Figure 5. Current-voltage and power-voltage characteristics.

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

Figure 6. Equivalent resistance as a function of temperature difference.

Figure 7. TEG equivalent resistance as a function of delta temperature.

generator with a variable series resistance, both proportional of the delta temperature:

$$R\_{eq} = m\_r \cdot \Delta T + q\_r \tag{6}$$

More TEG cells can be combined in order to increase the total generated power. Combination can be done by putting the cells in series, increasing output voltage, or, in parallel, increasing the output current. Due to regulatory constraint, it is better to keep the maximum voltage under 60 V. Furthermore, the TEG pack design is strictly related to available space, point of contact of the exhaust pipe, and heat cooling passing through the TEG (it will cause different available powers between near cells) [2, 3]. The TEG pack for this application is constituted by six cells in series and four in parallel. Figure 8 shows the equivalent electric model of the TEG pack constituted by 24 cells. Under the hypothesis of a good pack design, Eqs. (7)– (12) give a general relation between electrical parameters and the number of cells in series (Cs) and in parallel (Cp):

$$V\_{\rm ac}^{\prime} = V\_{\rm OC} \cdot \mathcal{C}\_{\rm s} \tag{7}$$

$$R\_{eq}^{\prime} = R\_{eq} \cdot \frac{\mathbf{C\_r}}{\mathbf{C\_p}} \tag{8}$$

$$P'\_{\text{max}} = \left(V'\_{oc}\right)^2 / R'\_{eq} \tag{9}$$

maximum power is a little more than 10 W. The coefficients mv and mp are specifics

Advanced Thermoelectric Materials for Energy Harvesting Applications

Figure 5 shows the voltage-current and the power-voltage characteristics at a single delta temperature. When the TEG output is closed in short circuit, the voltage is zero, and the output current reach is local maximum (Isc). It is important, with the purpose of generating the available maximum power, to know, for each delta temperature, the correct voltage-current couple. In Eq. (3) the relation between the voltage to generate the maximum power (Vmpp) and the open-circuit voltage is pointed out, namely, the delta temperature. In Eq. (4) the current to generate the

maximum power (Impp) is derived as a function of the already calculated

Vmpp <sup>¼</sup> VOC

Impp <sup>¼</sup> PMAX Vmpp

<sup>2</sup> <sup>¼</sup> mv

From the voltage-current characteristic, it is also possible to extract the value of the TEG equivalent series resistance (Eq. (5)). Figure 6 shows the voltage-current characteristic at different delta temperature: different slopes mean different equiv-

Req <sup>¼</sup> Isc

Figure 7 shows the equivalent series resistance calculated at different delta temperature. Eq. (6) shows the equivalent series resistance as a function of delta temperature. The equivalent electric model of a TEG cell is a variable voltage

 

�Voc

 

<sup>¼</sup> <sup>2</sup> � mp mv

VOC ¼ mv � ΔT (1) PMAX <sup>¼</sup> mp � <sup>Δ</sup>T<sup>2</sup> (2)

<sup>2</sup> � <sup>Δ</sup><sup>T</sup> (3)

� ΔT (4)

(5)

of the selected TEG cell:

parameters:

alent resistance:

Figure 5.

112

Current-voltage and power-voltage characteristics.

Advanced Thermoelectric Materials for Energy Harvesting Applications

$$I\_{sc}' = I\_{sc} \cdot \mathbf{C}\_p \tag{10}$$

is possible to combine these two topologies. From the combination of a buck and

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

In the boost and buck cascaded topology, there are two stages: the first one is a boost stage, and the second is a buck stage. Both stages have their own inductor (Figure 11). If the input voltage is lower than the output voltage, only the first stage works, while the high-side switch belonging to the second stage is always closed; on the other hand, if the input voltage is higher than the output voltage, only the second stage works, while the low-side switch belonging to the first stage is always open. As a third working option, if the input voltage is near the output voltage, both stages could work: the boost stage slightly increases the input voltage, and the buck

In the non-inverting buck-boost topology, the two stages are compressed, and

there is only one inductor that is shared by them (Figure 12). The first leg is

boost converters, there are two topologies:

DOI: http://dx.doi.org/10.5772/intechopen.86232

• Non-inverting buck-boost (NIBB)

stage regulates and complies with the desired output voltage.

• Boost and buck cascaded

Figure 10.

Figure 11.

Figure 12.

115

Boost and buck cascaded topology.

Non-inverting buck-boost (NIBB) topology.

Buck converter topology.

$$V\_{mpp}' = V\_{mpp} \cdot \mathcal{C}\_s \tag{11}$$

$$I\_{mpp}' = I\_{mpp} \cdot \mathbf{C}\_p \tag{12}$$

#### 2.2 DC-DC topology evaluation

Usually the aim of a DC-DC converter is to regulate an unregulated input voltage for a specific load; in this application the final goal is to take the maximum power from the generator at the larger possible temperature excursion, namely, manage a wide input voltage and set the current flows. A battery is connected to the output, so the output voltage must be slightly higher than the battery voltage, in order to allow the recharge. The battery voltage is set at 12 V. The input voltage is supposed to be in the range from 0 to 60 V.

Every DC-DC topology has its own relation between input and output voltages. In this test case, commonly used topologies, like buck or boost converter, are excluded. A buck converter (Figure 9) works only with input voltage higher than output; it means only with high delta temperature. A boost converter (Figure 10) works only with input voltage lower than output; it means only with low delta temperature. In order to manage input voltage lower and higher than the output, it

Figure 8. TEG cell configuration electrical equivalent.

Figure 9. Boost converter topology.

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

is possible to combine these two topologies. From the combination of a buck and boost converters, there are two topologies:

• Boost and buck cascaded

I 0

Usually the aim of a DC-DC converter is to regulate an unregulated input voltage for a specific load; in this application the final goal is to take the maximum power from the generator at the larger possible temperature excursion, namely, manage a wide input voltage and set the current flows. A battery is connected to the output, so the output voltage must be slightly higher than the battery voltage, in order to allow the recharge. The battery voltage is set at 12 V. The input voltage is supposed

Every DC-DC topology has its own relation between input and output voltages.

In this test case, commonly used topologies, like buck or boost converter, are excluded. A buck converter (Figure 9) works only with input voltage higher than output; it means only with high delta temperature. A boost converter (Figure 10) works only with input voltage lower than output; it means only with low delta temperature. In order to manage input voltage lower and higher than the output, it

V0

Advanced Thermoelectric Materials for Energy Harvesting Applications

I 0

2.2 DC-DC topology evaluation

to be in the range from 0 to 60 V.

Figure 8.

Figure 9.

114

Boost converter topology.

TEG cell configuration electrical equivalent.

sc ¼ Isc � Cp (10)

mpp ¼ Vmpp � Cs (11)

mpp ¼ Impp � Cp (12)

• Non-inverting buck-boost (NIBB)

In the boost and buck cascaded topology, there are two stages: the first one is a boost stage, and the second is a buck stage. Both stages have their own inductor (Figure 11). If the input voltage is lower than the output voltage, only the first stage works, while the high-side switch belonging to the second stage is always closed; on the other hand, if the input voltage is higher than the output voltage, only the second stage works, while the low-side switch belonging to the first stage is always open. As a third working option, if the input voltage is near the output voltage, both stages could work: the boost stage slightly increases the input voltage, and the buck stage regulates and complies with the desired output voltage.

In the non-inverting buck-boost topology, the two stages are compressed, and there is only one inductor that is shared by them (Figure 12). The first leg is

Figure 10. Buck converter topology.

Figure 11.

Boost and buck cascaded topology.

Figure 12. Non-inverting buck-boost (NIBB) topology.

dedicated for buck operation, while the second leg for the boost operation. As in the previous topology, it is possible to manage input voltage that is lower or higher than the output voltage. The three working regions are:

the boost leg (low-side switch); the command of the high-side switch is

ering in which region the converter is working (Eq. (14)):

8 >>>>>><

DOI: http://dx.doi.org/10.5772/intechopen.86232

>>>>>>:

Vo <sup>¼</sup> Vi 1 � DB

Vo <sup>¼</sup> Vi � DA 1 � DB

Vo <sup>¼</sup> Vi � DA 1 � DB

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

It is known from Section 2.1 that to extract the maximum power from the TEG, its output voltage must be half of the open-circuit voltage; in the following Vmpp is considered as the effective input voltage of the DC-DC converter. Figure 14 shows the relation between Vmpp and delta temperature, highlighting the working regions. The general relation between input and output voltages could be simplified consid-

Vo ¼ Vi � DA, DB ¼ 0 ð Þ Buck region

In order to define switching frequency and inductor value, the two operation modes of a DC-DC converter are introduced: continuous current mode (CCM) and discontinuous current mode (DCM). The DCM operation is characterized by a current of the inductor that is zero during the switching period, while in the CCM, the inductor current never reduces to zero (Figure 15). The DCM operation is particularly a disadvantage for the buck region. In order to avoid the DCM operation for the buck working region, Eq. (15) must be verified. The value of Lcrit is calculated in Eq. (16), where ΔIL is the current ripple on the inductor and can be considered a percentage of the output current, typically 20–30% of the output

Lcrit <sup>¼</sup> DA � Vi � Vo

Converter input voltage as a function of the temperature difference, highlighting working regions.

ΔIL � f SW

, DA ¼ 1 ð Þ Boost region

ð Þ Buck � Boost region

L>Lcrit (15)

(13)

(14)

(16)

1 � DB (DBÞ:

current:

Figure 14.

117


Practically the inductor is shared only in the buck-boost working region; on the other only one leg is active, and the other has a fixed state. In order to reduce weight, number of components, and loss sources, the NIBB topology is the perfect candidate for this application. In the next section, design steps and loss estimation of the NIBB topology are presented.

#### 2.3 NIBB design

In the non-inverting buck-boost topology, a buck and a boost are fused together and share a single inductor. A technique for improving the converter's efficiency is the synchronous rectification. The two versions of the NIBB converter are also called two-switch and four-switch non-inverting buck-boost. The two-switch NIBB (2SW-NIBB) has diodes, while in the four-switch NIBB (4SW-NIBB), the diodes are replaced by power MOSFETs (Figure 13). This topology will be used for the DC-DC converter, and in the following, the design steps are reported.

Below the starting specifications for the design are reported:


The steady-state relation between input and output voltages of a 4SW-NIBB is outlined in Eq. (13). DA and DB are the commands of the switches, namely, duty cycles; their range is from 0 to 1. DA is the duty cycle of the buck leg high-side switch; the command of the low-side switch is 1 � DA (DA Þ. DB is the duty cycle of

Figure 13. NIBB with synchronous rectification (4SW-NIBB) circuit.

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

the boost leg (low-side switch); the command of the high-side switch is 1 � DB (DBÞ:

$$V\_o = \frac{V\_i \cdot D\_A}{1 - D\_B} \tag{13}$$

It is known from Section 2.1 that to extract the maximum power from the TEG, its output voltage must be half of the open-circuit voltage; in the following Vmpp is considered as the effective input voltage of the DC-DC converter. Figure 14 shows the relation between Vmpp and delta temperature, highlighting the working regions. The general relation between input and output voltages could be simplified considering in which region the converter is working (Eq. (14)):

$$\begin{cases} V\_o = \frac{V\_i}{1 - D\_B}, D\_A = 1 \text{ (Boost region)}\\\\ V\_o = \frac{V\_i \cdot D\_A}{1 - D\_B} \text{ (Buck} - \text{Boost region)}\\\\ V\_o = V\_i \cdot D\_A, D\_B = 0 \text{ (Buck region)} \end{cases} \tag{14}$$

In order to define switching frequency and inductor value, the two operation modes of a DC-DC converter are introduced: continuous current mode (CCM) and discontinuous current mode (DCM). The DCM operation is characterized by a current of the inductor that is zero during the switching period, while in the CCM, the inductor current never reduces to zero (Figure 15). The DCM operation is particularly a disadvantage for the buck region. In order to avoid the DCM operation for the buck working region, Eq. (15) must be verified. The value of Lcrit is calculated in Eq. (16), where ΔIL is the current ripple on the inductor and can be considered a percentage of the output current, typically 20–30% of the output current:

$$L \succ L\_{crit} \tag{15}$$

$$L\_{crit} = D\_A \cdot \frac{V\_i - V\_o}{\Delta I\_L \cdot f\_{SW}} \tag{16}$$

Figure 14. Converter input voltage as a function of the temperature difference, highlighting working regions.

dedicated for buck operation, while the second leg for the boost operation. As in the previous topology, it is possible to manage input voltage that is lower or higher than

Practically the inductor is shared only in the buck-boost working region; on the

In the non-inverting buck-boost topology, a buck and a boost are fused together and share a single inductor. A technique for improving the converter's efficiency is the synchronous rectification. The two versions of the NIBB converter are also called two-switch and four-switch non-inverting buck-boost. The two-switch NIBB (2SW-NIBB) has diodes, while in the four-switch NIBB (4SW-NIBB), the diodes are replaced by power MOSFETs (Figure 13). This topology will be used for the

The steady-state relation between input and output voltages of a 4SW-NIBB is outlined in Eq. (13). DA and DB are the commands of the switches, namely, duty cycles; their range is from 0 to 1. DA is the duty cycle of the buck leg high-side switch; the command of the low-side switch is 1 � DA (DA Þ. DB is the duty cycle of

other only one leg is active, and the other has a fixed state. In order to reduce weight, number of components, and loss sources, the NIBB topology is the perfect candidate for this application. In the next section, design steps and loss estimation

DC-DC converter, and in the following, the design steps are reported. Below the starting specifications for the design are reported:

the output voltage. The three working regions are:

Advanced Thermoelectric Materials for Energy Harvesting Applications

• Buck (Vi > Vo)

• Boost (Vi < Vo)

2.3 NIBB design

• Vo ≈ 12 V.

Figure 13.

116

• Buck-boost (Vi ≈ Vo)

of the NIBB topology are presented.

• Input voltage 0 < Vi < 60 V.

• Maximum output power 300 W.

NIBB with synchronous rectification (4SW-NIBB) circuit.

• Output power is proportional to input voltage.

#### Advanced Thermoelectric Materials for Energy Harvesting Applications

Figure 16 shows the relation between Lcrit and the temperature difference considering different switching frequency. The use of high switching frequency allows the use of smaller inductor but means higher switching losses. The switching frequency and the inductor value are a trade-off between reduced component dimension and losses. The last important parameter for the inductor selection is the saturation current that is calculated in Eq. (17), where μ is the converter efficiency that as a safety factor can be set to 80%:

$$I\_{L, \text{sat}} = \frac{V\_o \cdot I\_o}{\mu \cdot V\_{in}} + \frac{\Delta I\_L}{2} \tag{17}$$

Input capacitor selection closely follows the method used for the output capacitor. The maximum RMS current flowing in CIN occurs in buck mode that is calculated in Eq. (20). As for the output capacitance, the input voltage ripple has two contributions: the first one related to input capacitor ESR and the second due to

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

Io � Vo Vin � Cin � f SW

2

2

2

2

At first approximation inductor and capacitors are due to the component's equivalent series resistance (ESR). In Eqs. (22) and (23), the input capacitor losses and output capacitor losses are calculated, respectively; inductor losses are calcu-

PCin ð Þ¼ ΔT ESRCin � I

PCo ð Þ¼ ΔT ESRCo � I

can be used in order to reduce switching losses (Figures 17 and 18):

PLð Þ¼ ΔT ESRL � I

Pcondð Þ¼ ΔT 2 � Rdson � I

Loss estimation of the passive components. Continuous line input capacitor; dashed line is the output capacitor;

The MOSFET losses have two components: conduction and switching losses. Conduction losses are calculated in Eq. (25). Rdson is the equivalent resistance of the MOSFET when it is closed. For the switching losses, it is necessary to separate the two legs: Eq. (26) for the buck leg and Eq. (27) for the boost leg. tR is the rise time and tF is the fall time, parameters of the selected switch. While output current is a constraint of the application, the switching frequency is a degree of freedom that

1 12 ΔI 2 L <sup>s</sup> � � (20)

<sup>1</sup> � Vo

Vin � � (21)

Cin,RMS ð Þ ΔT (22)

Co,RMS ð Þ ΔT (23)

<sup>o</sup>ð Þ ΔT (24)

<sup>o</sup>ð Þ ΔT (25)

capacitance (Eq. (21)):

2.3.1 Loss estimation

lated in Eq. (24):

Figure 17.

119

dash-dotted line is the inductor. Parameters of Table 1.

ICin,RMS ¼

DOI: http://dx.doi.org/10.5772/intechopen.86232

ΔVCin ¼ Io � ESR þ

Vo VIN I 2 <sup>o</sup> <sup>1</sup> � Vo Vin � � <sup>þ</sup>

The capacitor selection is based on maximum RMS current that is flown into it and the voltage ripple. The output capacitor selection is based on the maximum RMS current in the capacitor during the boost working region and the output voltage ripple specification. Maximum RMS current in the output capacitor is calculated in Eq. (18). Input voltage is the maximum input voltage for the boost working region. The output voltage ripple has two contributions: the first one due to capacitor equivalent series resistance (ESR) and the other related to capacitance (Eq. (19)):

$$I\_{\rm Co,RMS} = I\_o \sqrt{\frac{V\_o}{V\_{IN}}}\tag{18}$$

$$
\Delta V\_{Co} = \frac{I\_o V\_o}{V\_{in}} \cdot ESR + \frac{I\_o}{C\_o \cdot f\_{SW}} \left(\frac{V\_o - V\_{in}}{V\_o}\right) \tag{19}
$$

Figure 15. Inductor current waveform: (left) CCM and (right) DCM.

Figure 16. Minimum CCM inductance as a function of temperature difference.

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

Input capacitor selection closely follows the method used for the output capacitor. The maximum RMS current flowing in CIN occurs in buck mode that is calculated in Eq. (20). As for the output capacitance, the input voltage ripple has two contributions: the first one related to input capacitor ESR and the second due to capacitance (Eq. (21)):

$$\mathbf{I}\_{\text{C}\_{\text{in}},\text{RMS}} = \sqrt{\frac{V\_o}{V\_{IN}} \left[ I\_o^2 \left( 1 - \frac{V\_o}{V\_{in}} \right) + \frac{1}{12} \Delta I\_L^2 \right]} \tag{20}$$

$$
\Delta V\_{C\_{in}} = I\_o \cdot ESR + \frac{I\_o \cdot V\_o}{V\_{in} \cdot C\_{in} \cdot f\_{SW}} \left(1 - \frac{V\_o}{V\_{in}}\right) \tag{21}
$$

#### 2.3.1 Loss estimation

Figure 16 shows the relation between Lcrit and the temperature difference considering different switching frequency. The use of high switching frequency allows the use of smaller inductor but means higher switching losses. The switching frequency and the inductor value are a trade-off between reduced component dimension and losses. The last important parameter for the inductor selection is the saturation current that is calculated in Eq. (17), where μ is the converter efficiency

Advanced Thermoelectric Materials for Energy Harvesting Applications

IL,sat <sup>¼</sup> Vo � Io μ � Vin

ICO,RMS ¼ Io

� ESR þ

<sup>Δ</sup>VCO <sup>¼</sup> IoVo

Inductor current waveform: (left) CCM and (right) DCM.

Minimum CCM inductance as a function of temperature difference.

Vin

þ ΔIL

> ffiffiffiffiffiffiffiffi Vo VIN

> > Vo � Vin Vo � �

r

Io Co � f SW

The capacitor selection is based on maximum RMS current that is flown into it and the voltage ripple. The output capacitor selection is based on the maximum RMS current in the capacitor during the boost working region and the output voltage ripple specification. Maximum RMS current in the output capacitor is calculated in Eq. (18). Input voltage is the maximum input voltage for the boost working region. The output voltage ripple has two contributions: the first one due to capacitor equivalent series resistance (ESR) and the other related to capacitance

<sup>2</sup> (17)

(18)

(19)

that as a safety factor can be set to 80%:

(Eq. (19)):

Figure 15.

Figure 16.

118

At first approximation inductor and capacitors are due to the component's equivalent series resistance (ESR). In Eqs. (22) and (23), the input capacitor losses and output capacitor losses are calculated, respectively; inductor losses are calculated in Eq. (24):

$$\mathbf{P\_{C\_{in}}}(\Delta T) = E \mathbf{S} \mathbf{R\_{C\_{in}}} \cdot \mathbf{I\_{C\_{in,RMS}}^2}(\Delta T) \tag{22}$$

$$\mathbf{P\_{C\_o}}(\Delta T) = E \mathbf{S} \mathbf{R\_{C\_o}} \cdot \mathbf{I\_{C\_{o\_{\text{RMS}}}}}(\Delta T) \tag{23}$$

$$\mathbf{P}\_{L}(\Delta T) = E \mathbf{S} \mathbf{R}\_{L} \cdot \mathbf{I}\_{\mathbf{o}}^{2}(\Delta T) \tag{24}$$

The MOSFET losses have two components: conduction and switching losses. Conduction losses are calculated in Eq. (25). Rdson is the equivalent resistance of the MOSFET when it is closed. For the switching losses, it is necessary to separate the two legs: Eq. (26) for the buck leg and Eq. (27) for the boost leg. tR is the rise time and tF is the fall time, parameters of the selected switch. While output current is a constraint of the application, the switching frequency is a degree of freedom that can be used in order to reduce switching losses (Figures 17 and 18):

$$P\_{cond}(\Delta T) = \mathcal{Z} \cdot R\_{dom} \cdot \mathbf{l}\_o^2(\Delta T) \tag{25}$$

#### Figure 17.

Loss estimation of the passive components. Continuous line input capacitor; dashed line is the output capacitor; dash-dotted line is the inductor. Parameters of Table 1.

current are used for the maximum power point tracking; output voltage and current

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

There are many MPPT algorithms that are suitable for use with a TEG, most of which have been adapted from those used by photovoltaic systems [4]. The most common are perturb and observe (P&O) and incremental conductance (IC) algorithms [5–8]. P&O algorithm is an iterative algorithm, its working principle is to introduce a perturbation in the operating voltage and current, and then the operating power is observed and compared with the power of the previous step. If the power difference is positive (negative), subsequently the MPP will be reached if the perturbation is stepped in the same (opposite) direction. The flowchart for P&O is shown in Figure 19. A P&O algorithm cannot determine if the MPP is reached: when it does reach it, it will pass the operating voltage beyond which will result in a decrease of power; in response the algorithm will reverse the tracking direction. Therefore, in employing a P&O algorithm, the operating point oscillates around the MPP. In order to reduce the oscillation, it is possible to make the perturbation size smaller. However, reducing perturbation size means reducing the algorithm dynamics and increasing the time to reach the MPP. A trade-off between

steady-state performance and dynamic response is the use of variable perturbation

In order to improve the steady-state error, the incremental conductance (IC) algorithm was implemented. It tracks by stepping the voltage like the P&O, and it operates by incrementally comparing the ratio of the derivative of conductance (dI=dV) with the instantaneous conductance (I=V). The advantage of using IC method with respect to the P&O is that it can determine the MPP; in fact, at MPP, the derivative of power with respect to voltage (dP=dV) is 0. The relation between conductance and power-voltage ratio is expressed in Eq. (28). The basic rules for IC are written in Eq. (29). The IC algorithm flowchart is reported in Figure 20:

are used for the other control loops.

DOI: http://dx.doi.org/10.5772/intechopen.86232

2.4.1 Maximum power point tracking

size [9, 10].

Figure 19.

121

P&O algorithm flowchart.

#### Figure 18.

Loss estimation of the switches. Continuous line is the conduction losses; dashed lines are the switching losses. Parameters of Table 1.


#### Table 1.

Summary of the converter design.

$$P\_{sw}(\Delta T) = f\_{sw} \cdot (\mathbf{t}\_R + \mathbf{t}\_F) \cdot V\_{mpp}(\Delta T) \cdot I\_o(\Delta T) \tag{26}$$

$$P\_{sw}(\Delta T) = f\_{sw} \cdot (t\_R + t\_F) \cdot V\_o \cdot I\_o(\Delta T) \tag{27}$$

The final design is a trade-off between inductor and capacitor values, namely, dimension and weight, and the losses. Table 1 shows the summary of the design. Device's parameters are comparable with real components available on the market. These parameters are used for the control algorithms and simulation sections.

#### 2.4 Control algorithms

The control system consists of nested control loops. In order to extract the maximum power from the TEG, two different maximum power point tracking (MPPT) algorithms are evaluated. The MPPT output is used as reference for voltage and current control loops. Analog variables are input and output voltage and current of the converter; commands are the duty cycles of the legs. Input voltage and

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

current are used for the maximum power point tracking; output voltage and current are used for the other control loops.

#### 2.4.1 Maximum power point tracking

There are many MPPT algorithms that are suitable for use with a TEG, most of which have been adapted from those used by photovoltaic systems [4]. The most common are perturb and observe (P&O) and incremental conductance (IC) algorithms [5–8]. P&O algorithm is an iterative algorithm, its working principle is to introduce a perturbation in the operating voltage and current, and then the operating power is observed and compared with the power of the previous step. If the power difference is positive (negative), subsequently the MPP will be reached if the perturbation is stepped in the same (opposite) direction. The flowchart for P&O is shown in Figure 19. A P&O algorithm cannot determine if the MPP is reached: when it does reach it, it will pass the operating voltage beyond which will result in a decrease of power; in response the algorithm will reverse the tracking direction. Therefore, in employing a P&O algorithm, the operating point oscillates around the MPP. In order to reduce the oscillation, it is possible to make the perturbation size smaller. However, reducing perturbation size means reducing the algorithm dynamics and increasing the time to reach the MPP. A trade-off between steady-state performance and dynamic response is the use of variable perturbation size [9, 10].

In order to improve the steady-state error, the incremental conductance (IC) algorithm was implemented. It tracks by stepping the voltage like the P&O, and it operates by incrementally comparing the ratio of the derivative of conductance (dI=dV) with the instantaneous conductance (I=V). The advantage of using IC method with respect to the P&O is that it can determine the MPP; in fact, at MPP, the derivative of power with respect to voltage (dP=dV) is 0. The relation between conductance and power-voltage ratio is expressed in Eq. (28). The basic rules for IC are written in Eq. (29). The IC algorithm flowchart is reported in Figure 20:

Figure 19. P&O algorithm flowchart.

Pswð Þ¼ ΔT fsw � ð Þ� tR þ tF Vmppð Þ� ΔT Ioð Þ ΔT (26) Pswð Þ¼ ΔT fsw � ð Þ� tR þ tF Vo � Ioð Þ ΔT (27)

The final design is a trade-off between inductor and capacitor values, namely, dimension and weight, and the losses. Table 1 shows the summary of the design. Device's parameters are comparable with real components available on the market. These parameters are used for the control algorithms and simulation sections.

Loss estimation of the switches. Continuous line is the conduction losses; dashed lines are the switching losses.

Advanced Thermoelectric Materials for Energy Harvesting Applications

Parameter/component Value

fsw 30 kHz L 30 μH L ESR 20 mΩ Cin 660 μF ESR Cin 20 mΩ Cout 330 μF ESR Cout 40 mΩ Rds,ON 5 mΩ tR 350 ns tF 200 ns

The control system consists of nested control loops. In order to extract the maximum power from the TEG, two different maximum power point tracking (MPPT) algorithms are evaluated. The MPPT output is used as reference for voltage and current control loops. Analog variables are input and output voltage and current of the converter; commands are the duty cycles of the legs. Input voltage and

2.4 Control algorithms

Summary of the converter design.

Figure 18.

Table 1.

120

Parameters of Table 1.

Figure 20. IC algorithm flowchart.

$$\frac{dP}{dV} = \frac{d(VI)}{dV} = I + V \cdot \frac{dI}{dV} \tag{28}$$

$$\begin{cases} \qquad \frac{dI}{dV} = -\frac{I}{V}, \text{at } MPP \\\\ \qquad \frac{dI}{dV} > -\frac{I}{V}, \text{left of } MPP \\\\ \qquad \frac{dI}{dV} < -\frac{I}{V}, \text{right of } MPP \end{cases} \tag{29}$$

control loop. The average current of the inductor is related to the output current, and the stationary relation between input and output current is reported in Eq. (31):

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

Ii <sup>¼</sup> Io � DA 1 � DB

The system-state variables of a NIBB converter are reported in Eq. (32); the

<sup>x</sup> <sup>¼</sup> iL vC

� x þ

Vo

DB

The system is characterized by two control inputs (dA, dB), but as shown in Section 2.3, it is possible to consider them one at a time if the system is in buck or boost region. In order to have only one control variable also in the buck-boost region, it is possible to fix one leg, i.e., boost leg, and command only the other. Another way is to use the same command: the two duty cycles (dA ¼ dB). The buck transfer function between duty cycle and current is reported in Eq. (35), and the boost transfer function between duty cycle and current is reported in Eq. (36). The

Gdi <sup>¼</sup> <sup>R</sup>

tional integrative (PI) controller, the current reference comes from the MPPT algorithm, and the output current is the feedback. Figure 21 shows the system

ð Þ 1 þ sRCout

� sRC <sup>þ</sup> <sup>2</sup> Rð Þ 1 � DB

The proposed controller for the average current control loop is a simple propor-

The complete system is simulated in a model-based environment; thermoelectric generator is modeled using the mathematical model described in Section 2.1, while for the DC-DC converter and the load, it is used as a package to use circuital and power components inside the environment; values of Table 1 are used. Sense blocks are placed before and after the converter, and measures are used in the control block in order to generate the correct commands. The complete system functional blocks are reported in Figure 22. For the purpose of validating the single control algorithm, the blocks are tested separately. Figure 23 shows a detail of the P&O MPPT algorithm simulation: the algorithm's output is the equivalent resistance

sL � dBðÞ�<sup>s</sup>

sC � iLðÞ�<sup>s</sup>

DB sL � voð Þ<sup>s</sup>

1 sRC � voð Þ<sup>s</sup>

<sup>2</sup> <sup>þ</sup> sL <sup>þ</sup> <sup>s</sup><sup>2</sup>RLCo

The Laplace form of linearized model of the converter is reported in Eq. (34):

DA=<sup>L</sup> 0 � �

� � (32)

� Vi (33)

model of the system is expressed in matrixial form in Eq. (33):

<sup>x</sup>\_ <sup>¼</sup> <sup>0</sup> �DB=<sup>L</sup> DB <sup>C</sup> �1 =R�<sup>C</sup>

Vin

IL

effect of the output capacitor series resistance is neglected:

Gdi <sup>¼</sup> Vin

transfer function and the controller behavior.

2.5 Simulations and design validation

123

ð Þ 1 � DB

iLðÞ¼ s

8 >>><

DOI: http://dx.doi.org/10.5772/intechopen.86232

>>>:

voðÞ¼� s

" #

sL � dAðÞþ<sup>s</sup>

sC � dBðÞþ<sup>s</sup>

(31)

(34)

(35)

(36)

In practical application the equality condition rarely exists. Moreover, the conductance instantaneous value is related to perturbation size; because of that a marginal error (ε) is introduced. The IC algorithm considers the MPP reached when the operating point is within a certain error margin (Eq. (30)):

$$\left| I + V \cdot \frac{dI}{dV} \right| < \varepsilon \tag{30}$$

Using nested control system, the output perturbation of the MPPT algorithm is used as a reference for the internal control loops that are described in the next subsection.

#### 2.4.2 Current loop

A change of the TEG operating point and a change of converter output current are alike. In order to control the output current, an average current control loop is used. The output of the MPPT algorithm is used as a reference for the current

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

control loop. The average current of the inductor is related to the output current, and the stationary relation between input and output current is reported in Eq. (31):

$$I\_i = \frac{I\_o \cdot D\_A}{1 - D\_B} \tag{31}$$

The system-state variables of a NIBB converter are reported in Eq. (32); the model of the system is expressed in matrixial form in Eq. (33):

$$\mathbf{x} = \begin{bmatrix} \dot{\mathbf{z}}\_L \\ \mathbf{z}\_C \end{bmatrix} \tag{32}$$

$$
\dot{\boldsymbol{\omega}} = \begin{bmatrix} \mathbf{0} & -\overline{\boldsymbol{D}\_{\mathbb{R}} \!/ \!} \\ \frac{\mathbf{0}}{\mathbf{C}} & -\mathbf{1} \!/ \_{\mathbb{R}\cdot\!\!/} \end{bmatrix} \cdot \boldsymbol{\infty} + \begin{bmatrix} \boldsymbol{D}\_{\mathbb{A}} \!/ \_{\mathbb{L}} \\ \mathbf{0} \end{bmatrix} \cdot \mathbf{V}\_{i} \tag{33}
$$

The Laplace form of linearized model of the converter is reported in Eq. (34):

$$\begin{cases} \begin{aligned} \dot{u}\_L(\boldsymbol{\varsigma}) &= \frac{V\_{in}}{sL} \cdot d\_A(\boldsymbol{\varsigma}) + \frac{V\_o}{sL} \cdot d\_B(\boldsymbol{\varsigma}) - \frac{\overline{D\_B}}{sL} \cdot \boldsymbol{\nu}\_o(\boldsymbol{\varsigma}) \\\\ \boldsymbol{\upsilon}\_o(\boldsymbol{\varsigma}) &= -\frac{I\_L}{sC} \cdot d\_B(\boldsymbol{\varsigma}) + \frac{\overline{D\_B}}{sC} \cdot \dot{\boldsymbol{\iota}}\_L(\boldsymbol{\varsigma}) - \frac{1}{sRC} \cdot \boldsymbol{\nu}\_o(\boldsymbol{\varsigma}) \end{aligned} \end{cases} \tag{34}$$

The system is characterized by two control inputs (dA, dB), but as shown in Section 2.3, it is possible to consider them one at a time if the system is in buck or boost region. In order to have only one control variable also in the buck-boost region, it is possible to fix one leg, i.e., boost leg, and command only the other. Another way is to use the same command: the two duty cycles (dA ¼ dB). The buck transfer function between duty cycle and current is reported in Eq. (35), and the boost transfer function between duty cycle and current is reported in Eq. (36). The effect of the output capacitor series resistance is neglected:

$$G\_{di} = \frac{R}{(1 + sRC\_{out})} \tag{35}$$

$$G\_{di} = \frac{V\_{in}}{\left(1 - D\_B\right)} \cdot \frac{sRC + 2}{R\left(1 - D\_B\right)^2 + sL + s^2 RLC\_o} \tag{36}$$

The proposed controller for the average current control loop is a simple proportional integrative (PI) controller, the current reference comes from the MPPT algorithm, and the output current is the feedback. Figure 21 shows the system transfer function and the controller behavior.

#### 2.5 Simulations and design validation

The complete system is simulated in a model-based environment; thermoelectric generator is modeled using the mathematical model described in Section 2.1, while for the DC-DC converter and the load, it is used as a package to use circuital and power components inside the environment; values of Table 1 are used. Sense blocks are placed before and after the converter, and measures are used in the control block in order to generate the correct commands. The complete system functional blocks are reported in Figure 22. For the purpose of validating the single control algorithm, the blocks are tested separately. Figure 23 shows a detail of the P&O MPPT algorithm simulation: the algorithm's output is the equivalent resistance

dP dV <sup>¼</sup> d VI ð Þ

Advanced Thermoelectric Materials for Energy Harvesting Applications

8 >>>>>>><

>>>>>>>:

dI dV ¼ � <sup>I</sup>

dI dV <sup>&</sup>gt; � <sup>I</sup>

dI dV <sup>&</sup>lt; � <sup>I</sup>

the operating point is within a certain error margin (Eq. (30)):

subsection.

122

Figure 20.

IC algorithm flowchart.

2.4.2 Current loop

� � � � dV <sup>¼</sup> <sup>I</sup> <sup>þ</sup> <sup>V</sup> � dI

<sup>V</sup> , at MPP

<sup>V</sup> , left of MPP

<sup>V</sup> , right of MPP

� � �

In practical application the equality condition rarely exists. Moreover, the con-

ductance instantaneous value is related to perturbation size; because of that a marginal error (ε) is introduced. The IC algorithm considers the MPP reached when

> <sup>I</sup> <sup>þ</sup> <sup>V</sup> � dI dV

used as a reference for the internal control loops that are described in the next

Using nested control system, the output perturbation of the MPPT algorithm is

A change of the TEG operating point and a change of converter output current are alike. In order to control the output current, an average current control loop is used. The output of the MPPT algorithm is used as a reference for the current

dV (28)

� <sup>&</sup>lt; <sup>ε</sup> (30)

(29)

Figure 21. System transfer functions and controller transfer function.

Figure 24 shows the average current loop behavior with variable reference. The reference changes with steps like the behavior of the MPPT algorithm output. The reported current is the inductor current. In order to realize a nested control loop, it is important that the inner current loop band is larger than the outer MPPT loop. A NIBB prototype is realized using the designed components. In order to vali-

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power…

date the functionality and loss estimation, some power tests are carried out.

selected components have slowly better parameters.

3. Conclusions

125

Figure 24.

Figure 25.

Efficiency as a function of input voltage.

Average current control.

DOI: http://dx.doi.org/10.5772/intechopen.86232

Figure 25 shows the overlapping of the theoretical efficiency and the real efficiency obtained with open-loop commands. The difference from real data and estimation is due to the discrepancy of parasitic values: estimation plot is referred to Table 1; real

A comprehensive methodology for the design of a wasted heat energy recovery

system is given. The general approach is to model separately the single stages starting from physical characteristic and derivate the resulting design constraints. Starting from the characteristic of a thermoelectrical generator cell, the mathematical model and the equivalent electrical model are obtained; after this step it was possible to design the TEG pack. From this point the generator is considered a general input for the power converter. Different topologies are analyzed, and all the fundamental design steps are reported for a non-inverting buck-boost converter: as part of the design, a loss estimation analysis was carried out. The control system is described as nested control loops. Two maximum power point tracking algorithms

Figure 22. Wasted heat energy recovery system block diagram.

Figure 23. Maximum point tracking.

placed after the TEG. The operating point is fixed; the value of the resistance seen by the TEG is initialized at a wrong value, after some iteration; and according to the maximum power transfer theorem, the maximum power point is reached. During the steady-state operation, the resistance value continuously changes its value as seen in Section 2.4.1.

Harnessing the Automotive Waste Heat with Thermoelectric Modules Using Maximum Power… DOI: http://dx.doi.org/10.5772/intechopen.86232

Figure 24. Average current control.

Figure 25. Efficiency as a function of input voltage.

Figure 24 shows the average current loop behavior with variable reference. The reference changes with steps like the behavior of the MPPT algorithm output. The reported current is the inductor current. In order to realize a nested control loop, it is important that the inner current loop band is larger than the outer MPPT loop.

A NIBB prototype is realized using the designed components. In order to validate the functionality and loss estimation, some power tests are carried out. Figure 25 shows the overlapping of the theoretical efficiency and the real efficiency obtained with open-loop commands. The difference from real data and estimation is due to the discrepancy of parasitic values: estimation plot is referred to Table 1; real selected components have slowly better parameters.
