Bidirectional Modeling of Thermoelectric Module Using MATLAB/Simulink for Circuit Simulations

This paper presents a complete model for thermoelectric modules (TEMs) for use in circuit simulations using MATLAB/Simulink software. The purpose of this study is to develop a single block similar to other circuit elements in the Simulink library by combining both the thermoelectric generator (TEG) and thermoelectric cooler (TEC) module parameters and their operations. Since TEG and TEC are both thermoelectric devices, mostly whose main material, such as Bi₂Te₃, is the same for similar temperature ranges, a TEG can be operated as TEC and vice versa a TEC can be operated as TEG. In some systems, such as PV-thermoelectric (TE), fuel cell-TE, and battery temperature management of electric vehicles, modules can be used in both modes. In addition, in some cases of these systems, the TEM may need to be operated as a generator and sometimes as a cooler. The bidirectional operation of a TEM is a cost-saving approach because the two types of TEMs do not need to be installed together. This study proposes a novel model, whereas the studies to date in the literature have not been fully compatible with the bidirectional usage of modules. The proposed model deals with the polarities of both the thermal and electrical ports of the TEM to determine the mode of operation, that is, TEG or TEC. The maximum and performance parameters for both the TEG and TEC operations are included in the module block to observe and control the device under the desired operating conditions. It can be seen that the simulations with temperature-dependent parameters gave very close results to the graphics of the datasheet. However, simulations in which Thomson heat is not ignored have more deviations compared to the datasheet. Moreover, the simulations with constant material properties provided better results than those with Thomson heat in the cases of higher <inline-formula> <tex-math notation="LaTeX">$\Delta \text{T}$ </tex-math></inline-formula> levels. When <inline-formula> <tex-math notation="LaTeX">$\Delta \text{T}$ </tex-math></inline-formula> is 10&#x00B0;K, <inline-formula> <tex-math notation="LaTeX">$COP_{max}$ </tex-math></inline-formula> without Thomson heat was found to be 3.399, which is the same as the datasheet. Meanwhile, <inline-formula> <tex-math notation="LaTeX">$COP_{max}$ </tex-math></inline-formula> with Thomson heat was found to be 3.7147 and <inline-formula> <tex-math notation="LaTeX">$COP_{max}$ </tex-math></inline-formula> with constant properties was 3.8057. However, when <inline-formula> <tex-math notation="LaTeX">$\Delta \text{T}$ </tex-math></inline-formula> is 50&#x00B0;K, <inline-formula> <tex-math notation="LaTeX">$COP_{max}$ </tex-math></inline-formula> is 0.247 without Thomson heat, which is the same as the datasheet, but <inline-formula> <tex-math notation="LaTeX">$COP_{max}$ </tex-math></inline-formula> is 0.341 with Thomson heat, and <inline-formula> <tex-math notation="LaTeX">$COP_{max}$ </tex-math></inline-formula> is 0.268 with constant material properties. In addition, it is shown that under variable polarities of thermal and electrical ports, bidirectional operation of the TEM is successfully achieved, and the TEM is operated as both TEG and TEC in the same simulation.

Abstract licence: CC BY-NC-ND

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.