The Onsager-de Groot-Callen transport theory, implemented as a network model, is used to simulate the transient Harman method, which is widely used experimentally to determine all thermoelectric transport coefficients in a single measurement setup. It is shown that this method systematically overestimates the Seebeck coefficient for samples composed of two different materials. As a consequence, the figure of merit is also overestimated, if the thermal coupling of the measurement setup to the environment is weak. For a mixture of metal and semiconductor particles near metal percolation the figure of merit obtained by the Harman method is more than 100% too large. For a correct interpretation of the experimental data, information on composition and microstructure of the sample are indispensable.
IntroductionThermoelectric materials are important for energy harvesting, especially from waste heat [1,2]. In order to optimize the conversion into electricity, it is desirable to predict device properties by efficient computer simulations. This is one goal of the present paper. More fundamentally, we are going to point out that memory effects render the global response of composite materials, which are common among modern nanostructured thermoelectrics [3], highly complex.The theoretical description of transport processes goes back to Onsager [4,5]. Applied to thermoelectrics, this became known as the Onsager-de Groot-Callen theory [6]. A special case is the so-called constant property model (CPM), where the electric and heat conductivity, σ and κ, and the Seebeck coefficient α are assumed to be constant. This model has been studied in detail analytically in one-dimension [7,8] and will serve as a reference system for validation in this paper.Since analytic calculations are restricted to simple compounds, numerical models have been developed in order to describe inhomogeneous materials. Although these models are also based on the Onsager-de Groot-Callen theory, most of them do not fully describe the thermoelectric effects, since they do not include Joule heat and/or Peltier heat [9,10]. They were used to calculate either the heat and electrical conductances or the Seebeck coefficient [11][12][13]. More complex models [14] exist in the framework of drift-diffusion models, and they are applied e.g. for the simulation of generators in complex geometries [15].In this paper a simple and coherent way to discretize the Onsager transport theory for thermoelectric materials will be used, which includes all relevant effects and time-dependencies. It can be seen as a version of the finite difference method reviewed in [16]. Originally it was designed for the investigation of transport processes in nano-particle configurations [13,17,18], which were mapped to a network model. In this paper it will be applied to disordered bulk systems rather than particle agglomerates.The model, derived in section 2, enables us to study thermoelectric effects in geometries, which can not be solved analytically. It is validated by comparing simulation...