11.0 Mathematical Modeling of Thermoelectric Cooling Modules
11.1 INTRODUCTION: The operation of thermoelectric cooling devices may be described mathematically and device performance can readily be modeled on a personal computer. Since the semiconductor material used in module fabrication has several temperaturedependent properties, temperature effects on module operation must be considered if a realistic model is to be developed. We have not attempted to provide a highly detailed description of computer modeling techniques, but rather to present the basic algebraic expressions needed to simulate thermoelectric module performance.
Ferrotec America has performed a comprehensive analysis of many thermoelectric cooling modules over a wide temperature range. This study has resulted in the development of mathematical models that may be used to reliably predict module performance under typical operating conditions. Data presented herein is based on module operation in a normal air atmosphere with thermally conductive grease (heat sink compound) used at both hot and cold module interfaces. These conditions are applicable to the majority of thermoelectric cooling applications. It should be noted that for modules having metabolized external surfaces, slight performance improvement may be observed if the modules are mounted with solder as opposed to thermal grease. In addition, when modules are operated in a vacuum, a small to moderate performance increase may be seen, especially in the case of multistage devices.
11.2 TEMPERATUREDEPENDENT MATERIAL PROPERTIES: There are a number of parameters associated with thermoelectric materials and modules that normally would have to be considered in a mathematical model. However, since actual module test data was used to derive several important coefficients, certain factors may be neglected thereby simplifying the modeling process. Elements that must be incorporated into the model include the module’s effectiveSeebeck coefficient (SM), Electrical Resistance (RM), and Thermal Conductance (KM).
The values of SM, RM, and KM can be expressed mathematically by polynomial equations. The specified equation coefficients, applicable over a range of 100°C to +150°C, were derived from an industrystandard 71couple, 6ampere module. Other module configurations easily can be modeled by applying a simple correction factor as described in paragraph 11.2.4. Note that when using the various equations, temperature values must be stated in degrees Kelvin.
An alternative method for estimating temperaturedependent module properties, which may be useful under certain circumstances, involves the use of tabulated module data. Values representing average SM, RM, and KM characteristics for selected modules over a wide temperature range will be found in Appendix A at the end of this manual. Although somewhat less accurate than using calculated values, this method provides a relatively simple approach to predicting module performance.
11.2.1 SEEBECK COEFFICIENT: When a temperature differential is maintained across a thermoelectric device, a voltage can be detected at the input terminals. The magnitude of the resultant voltage, called the Seebeck emf, is proportional to the magnitude of the temperature difference. The Seebeck coefficient, as a function of temperature, can be expressed as a third order polynomial:
S_{M} = s_{1} + s_{2}T + s_{3}T^{2} + s_{4}T^{3}
T is the average module temperature in °K
Coefficients for a 71cpl, 6amp module
s_{1} = 1.33450 × 10^{2}
s_{2} = 5.37574 × 10^{5}
s_{3} = 7.42731 × 10^{7}
s_{4} = 1.27141 × 10^{9}
The above polynomial expression represent the Seebeck coefficient when the temperature difference across the module is zero (DT = Th – Tc = 0). When DT>0, the Seebeck coefficient must be evaluated at both temperatures Th and Tc using the expressions:
S_{MTh} or S_{MTc} = s_{1}T +  s_{2}T^{2}  +  s_{3}T^{3}  +  s_{4}T^{4} 
2  3  4 
S_{M} = (S_{MTh} – S_{MTc}) / DT
S_{MTc} is the module’s Seebeck coefficient at the cold side temperature Tc
11.2.2 MODULE RESISTANCE: The electrical resistance of a thermoelectric module, as a function of temperature, can be expressed as third order polynomials for the two conditions (a) and (b):
(a) when DT = 0:  R_{M} = r_{1} + r_{2}T + r_{3}T^{2} + r_{4}T^{3}  
(b) when DT > 0: 



R_{MTh} is the module’s resistance at the hot side temperature Th
R_{MTc} is the module’s resistance at the cold side temperature Tc
T is the average module temperature in °K
Coefficients for a 71cpl, 6amp module
r_{1} = 2.08317
r_{2} = 1.98763 × 10^{2}
r_{3} = 8.53832 × 10^{5}
r_{4} = 9.03143 × 10^{8}
11.2.3 MODULE THERMAL CONDUCTANCE: The thermal conductance of a thermoelectric module, as a function of temperature, can be expressed as third order polynomials for the two conditions (a) and (b):
(a) when DT = 0:  K_{M} = k_{1} + k_{2}T + k_{3}T^{2} + k_{4}T^{3}  
(b) when DT > 0: 



K_{MTh} is the thermal conductance at the hot side temperature Th
K_{MTc} is the thermal conductance at the cold side temperature Tc
T is the average module temperature in °K
coefficients for a 71cpl, 6amp module
k_{1} = 4.76218 × 10^{1}
k_{2} = 3.89821 × 10^{6}
k_{3} = 8.64864 × 10^{6}
k_{4} = 2.20869 × 10^{8}
11.2.4 PARAMETER CONVERSIONS FOR OTHER MODULE CONFIGURATIONS: The SM, RM, and KM parameters shown are calculated for a 71couple, 6ampere thermoelectric module. If a new or different module configuration is to be modeled, it is necessary to apply a conversion factor to each of these parameters, as follows:





R_{new} is the electrical resistance of the new module
K_{new} is the thermal conductance of the new module
N_{new} is the number of couples in the new module
I_{new} is the optimum or maximum current of the new module
11.3 CALCULATION OF THERMOELECTRIC MODULE PERFORMANCE: There are five variable parameters applicable to a thermoelectric module that affect its operation. These parameters include:
I –  the input current to the module expressed in amperes 
Vin –  the input voltage to the module expressed in volts 
Th –  the hot side temperature of the module expressed in °K 
Tc –  the cold side temperature of the module expressed in °K 
Qc –  the heat input to (or heat pumped by) the module expressed in watts 
In order to calculate module performance it is necessary to set at least three of these variables to specific values. Two common calculation schemes involve either (a) fixing the values of Th, I, and Qc or, (b) fixing the values of Th, I and Tc. For the computeroriented individual, a relatively straightforward calculation routine can be developed to incrementally step through a series of fixed values to produce an output of module performance over a range of operating conditions.
11.4 SINGLESTAGE MODULE CALCULATIONS: These equations mathematically describe the performance of a singlestage thermoelectric module as illustrated in Figure (11l). When entering numerical data, do not forget that temperature values must be expressed in degrees Kelvin (°K). Calculations of the various parameters should be performed in the order shown.
Figure (111)
a) The temperature difference (DT) across the module in °K or °C is:
DT = T_{h} – T_{c}
b) Heat pumped (Qc) by the module in watts is:
Q_{c} = (S_{M} × T_{c} × I) – (0.5 × I^{2} × R_{M}) – (K_{M} × DT)
c) The input voltage (Vin) to the module in volts is:
V_{in }= (S_{M} × DT) + (I × R_{M})
d) The electrical input power (Pin) to the module in watts is:
P_{in} = V_{in} × I
e) The heat rejected by the module (Qh) in watts is:
Q_{h} = P_{in} + Q_{c}
f) The coefficient of performance (COP) as a refrigerator is:
COP = Q_{c} / P_{in}
11.5 HEATING MODE OPERATION: Thermoelectric modules may be operated in the heating mode by reversing the polarity of the applied DC power. When used in this manner, the TE module functions as a “heat pump” and heating efficiencies in excess of 100 percent may be realized under certain conditions. A rapid increase in temperature occurs when heating a smallmass object, and care must be taken to avoid overheating either the module or object. In the heating mode, illustrated in Figure (112), the heat sink and object effectively are in opposite positions whereby the object is now at temperature (Th) and the heat sink is at temperature (Tc).
Figure (112) a) Heat flow to the object (Qh) is given by the expression:
Q_{h} = (S_{M} × T_{h} × I) + (0.5 × I^{2} × R_{M}) – (K_{M} × (T_{h} – T_{c}))
b) The coefficient of performance as a heater (COPH) is:
COPH = Q_{h} / P_{in}
11.5.1 Heating mode performance of a standard 71Couple, 6Ampere module is presented graphically in Figures (113) and (114). These graphs illustrate module performance at a heat sink temperature of 25°C.
Figure (113)
Figure (114)
11.6 OTHER THERMOELECTRIC DEVICE ATTRIBUTES: There are many other properties of thermoelectric devices that can be described mathematically. Several characteristics that might be of interest for specific situations follow. Remember that temperature values must be expressed in °K.
a) The maximum heat pumping capacity (Qmax) in watts of a thermoelectric module is given by the following expression. Note that DT =0 at the maximum Qc condition and, therefore, Tc = Th.
Q_{max} =  S_{M}^{2} × T_{C}^{2} 
2 × R_{M} 
b) The maximum temperature differential (DTmax) in °K may be expressed as shown below. To obtain an accurate DTmax value, however, it will be necessary to perform an iterative series of calculations comparing Tc to DTmax at a fixed value of Th.
DT_{max} =  S_{M}^{2} × T_{C}^{2} 
2 × R_{M} 
c) The FigureofMerit (Z) is a measure of the overall performance of a thermoelectric device or material. Z always is higher for raw thermoelectric semiconductor material than for an actual module functioning within a thermal system. Since an operating module is affected by interface, conductive, convective, and other losses, the effective FigureofMerit is less than that of the raw material. The FigureofMerit may be expressed:





p is the electrical resistivity of the material in ohmcm
k is the thermal conductivity of the material in w/cm°K
d) The optimum current (Iopt) in amperes required to produce the maximum heat removal rate (Qmax) is:


l is the length (height) of an individual thermoelectric element in centimeters
R is the resistance of an individual thermoelectric element in ohms.
11.7 MODELING OF COMPLETE THERMOELECTRIC COOLING SYSTEMS: Information presented in the foregoing paragraphs describes the mathematical modeling of thermoelectric modules as opposed to complete thermal systems. By incorporating the module calculations into a more sophisticated system model, it is possible to accurately simulate the overall thermal performance. Two heat leak sources that must not be overlooked in a complete thermal model include (a) heat conduction between the cooled object and heat sink, and, (b) heat conduction through clamping screws, if any, that physically connect the heat sink and cooled object.
Heat conduction between the heat sink and object generally involves the transfer of heat through the air gap surrounding the module mounting area. The actual heat leak value can be calculated using the equation in paragraph 8.4.1 where area (A) is the “open” surface area not covered by the thermoelectric modules, distance (x) is the width of the air gap, and thermal conductivity (K) is the value for air.
Heat conduction through the clamping screws also can be calculated by means of the same equation. In this case, area (A) is the crosssectional area of all mounting screws calculated from the screw’s pitch diameter and (K) is the thermal conductivity of the screw material.
Table of Contents
Do You Need More Information?
Contact your local Ferrotec thermal solutions representative.