Расширенные методы Z-метрии

Расширенные методы Z-метрии
В.Володин, Г.Громов, Л.Ершова, Д.Кондратьев, А. Огрызко, А. Рогов
РМТ Лтд. Россия 119991, Москва, Ленинский проспект 53.
тел: 095-132-6817 факс: 095-132-5870
e-mail: [email protected] http://www.rmtltd.ru
Аннотация
Представляется устройство DX3065 для измерения параметров термоэлектрических
охладителей (ТЭО), а также сопутствующее программное обеспечение.
Прибор дает
возможность решить три проблемы, существующие в Z-метрии: необходимость малости
величины электрического тока для каждого измеряемого ТЭО, достаточность времени
эксперимента для выхода ТЭО в стационарный режим и учет неадиабатичности измерений с
помощью поправок и биполярного усреднения. Прибор позволяет измерить сопротивление,
константу времени и термоэлектрическую добротность ТЭО. Устройство DX3065 может
использоваться как в качестве оперативного средства для экспресс-сертификации ТЭО и
контроля качества ТЭО в массовом производстве.
Введение
Идея Z-метрии [1,2,3,4] подразумевает выполнение комплекса условий. К тепловым
условиям относится незначительность Джоулевых потерь и взаимодействия с окружающей
средой. Временное условие – измерение ТЭО в стационарном режиме, ибо математическая
база метода есть балансные стационарные уравнения [2]. В приборе DX3065 реализуются все
три условия благодаря возможности сканирования по току, расчетным поправкам и
считыванию телеметрических данных во времени.
Прибор DX3065 позволяет получить:
•
Сопротивление на переменном токе (R)
•
Временную константу ТЭО, характеризующую выход на режим при данных условиях
измерения
•
Термоэлектрическую добротность (Z)
1
Прибор позволяет испытывать различные типы однокаскадных ТЭО.
Кроме того, с помощью устройства DX3065 возможно оценить значение Z двухкаскадного
ТЭО, а также временную константу и сопротивление многокаскадных модулей.
Работа построена следующим образом.
В Разделе 1 обсуждается метод, задействованный в приборе DX3065 при измерении
термоэлектрической добротности Z. Приведено математическое обоснование по усреднению
результата измерений
для различных полярностей тока [4], что иллюстрируется
экспериментальными данными. Экспериментально оценивается корректность вводимых
поправок.
В разделе 2 описывается применение DX3065 для измерения константы времени ТЭО.
Обсуждаются экспериментальные результаты.
Раздел 3 охватывает существенные технические аспекты измерения сопротивления модуля и
напряжения термо-ЭДС. Показан внешний вид прибора и программной оболочки.
Завершает работу заключение.
1. Измерение термоэлектрической добротности:
биполярное усреднение и сканирование по току
При измерении термоэлектрической добротности Z ТЭО, модуль исследуется в незамкнутой
системе. Имеет место теплообмен с окружающей средой, тепло Джоуля не равно нулю.
Следовательно, необходимо учесть некоторые поправки [3] для того, чтобы применить
выражения Хармана [2] для получения числа Z Z=α2/kR=α2/κρ [5] (α - коэффициент Зеебека;
κ и ρ - теплопроводимость и электрическое удельное сопротивление модуля соответственно)
и максимально достижимую ТЭО разность температур ∆Tmax.
В работе предлагается единая теоретико-экспериментальная основа для выполнения условий
измерения числа Z ТЭО в рамках метода Хармана.
2
Обычно уравнения теплового баланса в однокаскадном случае записываются следующим
образом:
1 2

αIT0 − 2 I R − k' ∆T = a 1 ( Ta − T0 ) / N

1
αIT1 + I 2 R − k' ∆T = a 2 ( T1 − Ta ) / N
2

,
(1.1)
где N число ветвей модуля, а1 , a2 – коэффициенты теплоотдачи с холодной и горячей сторон
ТЭО ( a1 ≠ a 2 ), k ′ эффективная теплопроводность ветви с учетом газовой среды (воздуха) и
теплового излучения между ветвями.
Коэффициент k ' описывается так:
k' = k ( 1 + bth )
,
(1.2)
где
bth = Bcond + Brad
Величины Bcond и
Brad
,
(1.3)
- поправки, учитывающие дополнительную теплопроводность по
воздуху и излучение между ветвями, соответственно:
Bcond =
κ air 1
( − 1) ,
κ β
(1.4)
Здесь коэффициент заполнения ветвями есть:
β=
Ns
S
,
(1.5)
где S – площадь холодной стороны ТЭО;
Brad = γ
S
σTa3 ( 1 − β ) ,
Nk
(1.6)
где σ - константа Больцмана, γ - испускательная способность материала керамической
подложки. Выражение (1.6) может рассматриваться, как некая приближенная оценка.
Конвекционный вклад в теплопроводность между ветвями отсутствует, как следует из
критериев Грасгофа и Прандтля для характерных размеров зазоров между ветвями [6].
В работе [8] уравнения (1.1) решались в рамках следующих приближений:
3
a2
a1
<< k ′
<< k ′ ,
N
N
и I << k ′ ,
(1.7)
α
И было найдено следующее [7,8]:

 I 2 RN 
a1 a 2
a + a2 
αI (a 2 − a1 ) ,
 +
T = Ta  1 +
1 + 1
 + Ta
N (a1 + a 2 )k ′  a1 + a 2 
4 nk ′ 
2k ′ (a1 + a 2 )

∆T =
αI 
a1 − a 2 I 2 R
I 2 RN 
 Ta +
 +
k′ 
a1 + a 2  k ′(a1 + a 2 ) 2
(1.8)
(1.9)
Требованиям (1.7) должно быть уделено серьезное внимание. Они описывают ограничение
на взаимодействие измеряемого модуля с окружающей средой и на величину электрического
тока, при котором проводятся измерения. Требование малости (1.7) означает, что количество
тепла Пельтье, должно быть малым по сравнению с теплоотдачей в среду. Кроме того, мы
видим, что (1.9) содержит величину
I 2 RN
- характеристику теплоты Джоуля.
a1 + a 2
Существенно, что ограничение для тока по теплоте Джоуля
I <<
( a1 + a2 )Ta
I
~ ( a1 + a 2 )Ta max
NR
U max
(1.10)
Играет более важную роль, чем в неравенстве для тока (1.7). Это проиллюстрировано на рис.
(1.1a) и (1.1б).
4
Критерии Пельтье и Джоуля для ТЭО с разной высотой ветви
2,5
Высота ветви, мм
2,0
I (a/k), A
I(Джоуль), А
1,5
1,0
0,5
0,0
1
2
3
4
5
6
Рис.1.1а. Сравнительная диаграмма требований Пельтье и Джоуля, накладываемых на значения
измерительного электрического тока для ТЭО, отличающихся высотой ветвей. Испытываемые
ТЭО - 1MC04-004-(05 – 15).
Критерии Пельтье и Джоуля
для ТЭО с разным сечением ветви
3
2,5
2
ширина ветви, мм
I(a/k), A
I(Джоуль), A
1,5
1
0,5
0
1
2
3
4
Рис.1.1б. Сравнительная диаграмма требований Пельтье и Джоуля, накладываемых на значения
измерительного электрического тока для ТЭО, отличающихся сечением ветвей. Испытываемые
ТЭО - 1MT(03 – 07)-008-13.
Как видно из рис. (1.1a,б) требованием малости тока должен служить критерий Джоуля
(1.10).
Аппаратная нижняя граница измерений не учитывается ни одним из вышеприведенных
математических рассуждений и будет рассмотрена далее.
С учетом (1.9) мы получаем следующие выражения для термо-ЭДС и отношения величин
напряжений:
5
α 2 IN 
I 2 RN
 Ta +
k′ 
a1 + a 2
Uα =
 (a1 − a 2 )I 2 RNα
 +
2 k ′(a1 + a 2 )

(1.11)
Uα α 2 
I 2 RN  (a1 − a 2 )Iα ,
 Ta +
+
=
U R k ′R 
a1 + a 2  2 k ′(a1 + a 2 )
(1.12)
или

Uα
I 2 RN
= Z'  Ta +
UR
a1 + a 2

 (a1 − a 2 )Iα ,
 +
 2 k ′(a1 + a 2 )
(1.13)
где Z ′ = α .
2
k ′R
Таким образом, Z (=α2/kR) может быть найдена как
Z=
U α

1
( 1 + bth )( 1 + br ) + b A  ,

Ta ( 1 + bT ) U R

(1.14)
где
2
1) bT = 1 I RN - поправка к температуре окружающей среды
Ta a1 + a 2
2) bA =
(a2 − a1 )Iα
(a2 + a1 )2k
- поправка, учитывающая асимметрию холодная-горячая сторона модуля;
3) bth = Bcond + Brad - поправка, учитывающая дополнительную теплопроводность по воздуху и
излучение между ветвями (в соответствие с (1.4)-(1.6));
4) br =
r
- поправка на сопротивление r подводящих проводов
NR
Благодаря поправкам, приведенным выше, выражение (1.14) дает возможность учесть
влияние воздействия измерений в DX3065.
Сравнительный анализ вклада поправок в результирующее значение Z был ранее
осуществлен в работах [4,8]. В работе [8] были также рассмотрены две проблемы: измерение
Z для системы ТЭО+радиатор и двухкаскадных ТЭО.
6
Безусловно, допущения, сделанные в теории, накладывают многие ограничения на
применения метода измерения Z в смонтированных сборках ТЭО. Однако можно выделить
два возможных выхода:
1. Учет нежелательных эффектов с помощью поправок (1.14). Нужно, однако иметь в
виду, что всем поправкам присуща некая доля неопределенности, так как известны лишь
приближенные значения величин, входящих в расчет.
2. В работе [4] было рекомендовано так называемое биполярное усреднение. Покажем,
что этот метод позволяет исключить поправку bA (1.14). Рассмотрим выражения (1.12) (1.14).
Как видим, bA представляет собой линейную функцию тока. Обозначая полярности (+) и
(-), имеем:
 Uα

UR

I 2 RN  (a1 − a 2 )Iα
α2 
 =
 Ta +
+
a1 + a 2  2 k ′(a1 + a 2 )
 + k ′R 
(1.15a)
Uα

U R

α2 
I 2 RN
 =
 Ta +
a1 + a 2
 − k ′R 
 (a 2 − a 1 )Iα
 +
 2k ′(a 1 + a 2 )
(1.15b)
Суммируя выражения (1.15a) и (1.15b), получаем:
 Uα

U R

U 

I 2 RN
 +  α  = 2 Z ′  Ta +
a1 + a 2
+ U R −




(1.16)
Таким образом, величина Z может быть найдена как:



1
U
α
Z=
1 + bth (1 + br )


T a ( 1 + bT )  U R  averaged


(
)
(1.17)
На рис. 1.2a,б можно видеть результаты биполярного тестирования и усреднения
(Z
av
=
1
Ta
Uα 
) в зависимости от тока с учетом и без учета поправок (1.14).


 U R  averaged
7
Z*1000, 1/K
Биполярное тестирование без поправок
2,8
2,7
2,6
2,5
2,4
2,3
2,2
2,1
2
1,9
Z(-)
Z(+)
Zav
0
20
40
60
80
100
I, мA
Рис.1.2a.Результирующая термоэлектрическая добротность при сканировании по току без
поправок для однокаскадного ТЭО (1MC06-032-15).
Z*1000, 1/K
Биполярное тестирование с учетом поправок
3
2,9
2,8
2,7
2,6
2,5
2,4
2,3
2,2
2,1
2
Z(-)+corr
Z(+)+corr
Zav+corr
0
20
40
60
80
100
I, мA
Рис.1.2б. Результирующая термоэлектрическая добротность при сканировании по току с учетом
поправок (1MC06-032-15).
Вооруженные возможностью сканировать ток, мы можем оценить правильность поправки на
асимметрию с помощью следующей формулы:
 Uα

UR
  Uα 
(a − a1 )Iα
 − 
 = 2
 +  U R  − k ′(a1 + a 2 )
8
(1. 18)
На рис. 1.3 сопоставлены поправка на асимметрию, полученная из измеренных отношений
напряжений (левая часть (1.18)) и теоретическое выражение для данной поправки (правая
часть (1.18)).
Сравнение значений поправки на асимметрию,
полученных теоретически и экспериментально
0,02
0,00
-0,02
-0,04
-0,06
-0,08
-0,10
-0,12
-0,14
D-ratioU(+-)
bA*2
0
20
40
60
80
100
I, mA
Рис. 1.3. Пример сравнения теоретически и экспериментально полученной зависимости поправки на
асимметрию от тока (ТЭО 1MC06-032-15). Конвекционный коэффициент теплоотдачи 4Вт/м2K,
площади холодной и горячей сторон 64мм2 и 80мм2 соответственно, удельная теплопроводность
материала ветвей 1.35Вт/мK.
Видно, что расчетная поправка на асимметрию не оправдывает себя при малых токах
(примерно меньше 10мА). Исключение этой поправки с помощью биполярного усреднения,
действительно, разумно.
По-видимому, отклонение на малых токах может быть объяснено чисто аппаратными
явлениями. К этому выводу позволяет прийти тот факт, что для разнообразия тестированных
однокаскадных ТЭО рассматриваемое отклонение затухает на токах выше 10мA. Как
следствие мы предлагаем выбирать токи для измерения величины Z не ниже 5мA.
Можно показать, что вышеизложенная теория справедлива для системы ТЭО+радиатор в
случае a0 / aheat sin k << 1 , что, фактически, всегда имеет место.
Простейшая система уравнения для двухкаскадного ТЭО [8]:
9
1 2

αIT0 − 2 I R − k' 1 ( T1 − T0 ) =

1
αIT2 + I 2 R − k' 2 ( T2 − T1 ) =
2

a1
( Ta − T0 )
N1
,
a2
( T2 − Ta )
N2
(1.19)
где T0 , T1 и T2 - температуры, соответственно, холодной, промежуточной и горячей
подложки, a1 и a2 полные коэффициенты теплоотдачи в среду с холодной и горячей сторон.
Измерение Z для двухкаскадного ТЭО допускается лишь при выполнении требований [8]:
a1 a2
N s
=
= A = const , β i = i = const , k '1 = k ' 2 = k '
Si
N1 N 2
(1.20)
Поэтому биполярное тестирование также имеет смысл для двухкаскадного ТЭО с точки
зрения уточнения полученных данных.
К сожалению, математическая обработка двухкаскадного случая, подобная
выше
приведенной для однокаскадного случая (см. уравнения (1.1) – (1.14)), дающая
непосредственную корреляцию между средней температурой модуля и температурой
окружающей среды не столь проста и до сих пор кажется невозможной. В двухкаскадном
случае с учетом (1.20) мы получаем [8]:
ZT = ( 1 + bth )( 1 + br )
Здесь T =
T2 + T0
2
Uα
.
UR
(1.21)
есть средняя температура, которую полагают равной температуре
окружающей среды. Поправка br удовлетворяет (1.14).
bth = Bcond + Aconv + ( Brad + Arad )
(1.22)
Параметры Bcond , Bconv , Brad описаны в (1.4-1.6). Конвекция и излучение с внешних
поверхностей есть [8]:
A conv =
al
(1.23a)
κβ
10
A rad =
γ
σTa3 l
κβ
(1.23б)
На рис. 1.4a,б можно видеть результаты биполярного тестирования и усреднения
(Z
av
=
1
Ta
Uα 
) для двухкаскадного ТЭО в зависимости от тока с учетом и без учета


 U R  averaged
поправок (1.21) – (1.23).
Биполярное тестирование без поправок
2,7
Zx1000, 1/K
2,5
2,3
2,1
Z(-)
Z(+)
Zav
1,9
1,7
1,5
0
10
20
30
40
50
60
I, мA
Рис.1.4a. Результирующая термоэлектрическая добротность при сканировании по току без
поправок для двухкаскадного ТЭО (2MC06-010-10).
Биполярное тестирование с поправками
3,2
Zx1000, 1/K
3,0
Z(-)+corr
Z(+)+corr
Zav+corr
2,8
2,6
2,4
2,2
2,0
0
10
20
30
40
50
60
I, мA
Рис.1.4б. Результирующая термоэлектрическая добротность при сканировании по току с
поправками (1.21) для двухкаскадного ТЭО (2MC06-010-10).
Как же ответить на вопрос: какой должен быть ток измерений при определении величины Z?
Для измерения однокаскадных и двухкаскадных ТЭО ответ таков. Следует взять значение
11
усредненной Z, соответствующее токам, при которых измерения на той и другой полярности
не испытывают нефизического возрастания.
Для однокаскадного случая 1МС06-032-15, рассмотренного выше, мы исключаем
нефизические значения Z при значениях тока менее 5мА и получаем искомое значение Z.
Для ТЭО 1MC06-032-15 (см. рис. 1.2б) Z=2.6·10-3K-1.
Для приводившегося примера двухкаскадного ТЭО аналогичная процедура позволяет найти
Z=2.4·10-3K-1 (см. рис.1.4б).
2. Измерение константы времени
Как показано в разделе 1, вся математика измерений Z основана на балансных уравнениях,
что справедливо лишь для времен, много больших чем характеристическое время системы. С
учетом этого должны выбираться временные ограничения прибора.
Нестационарное охлаждение для единичной ветви и однокаскадного ТЭО изучалось в ряде
работ [9,10].
Было показано [11], что решение соответствующих дифференциальных
уравнений может быть записано в виде суперпозиции бесконечного числа экспонент:
∞
∆T (t , x ) = ∑ ( AnU n ( x ) )e − m n t + ∆Tst ( x )
(2.1)
i =1
Здесь ∆T(t,x)=T-Ta, T - температура точки ветви ТЭО с обобщенной координатой х в момент
времени t, Un и mn – собственные функции и собственные значения, An – тепловые
амплитуды, ∆Tst(x) – решение в стационарном случае.
Анализ решения (2.1) показывает, что процесс охлаждения можно условно разделить на два
этапа: нерегулярный и регулярный[11]. Первый этап детерминирован начальными условиями
и описывается суммой множества экспонент. Эта фаза довольно быстро затухает, в случае
если теплопроводность ветвей ТЭО достаточно велика. По прошествии первого этапа
динамикa описывается одной экспонентой (регулярный режим). Это явление и нашло свое
применение в DX3065.
12
C учетом вышесказанного для оценок имеет смысл пользоваться следующим качественным
соотношением для временной константы ветви τ = 1 / mmin [3,12]:
τ0 ~
l 2 dc
κ
(2.2)
,
где l, d, κ, c высота ветви, плотность, теплопроводимость и удельная теплоемкость материала
ветви соответственно.
Из качественных соображений для однокаскадного ТЭО с учетом теплоемкости С объекта на
холодной стороне можно записать:
τ~
lC
,
κsN
(2.3)
Выражение (2.3) есть не что иное, как отношение теплоемкости охлаждаемого объекта к
полной теплопроводности всех ветвей.
В работах [10,12] было также показано, что характеристическое время обнаруживает
зависимость от электрического тока. DX3065 дает возможность проверить это.
Процедура обработки временных данных в приборе такова.
Процесс выхода в стационарный режим однокаскадного ТЭО регистрируем посредством
измерения термо-ЭДС U α :
U α = α∆T
(2.4)
Измерения проводятся для обеих полярностей. Время сбора данных и временной шаг можно
варьировать. Полученные экспериментальные данные затем интерполируются с помощью
функции:
U α (t ) = Ustα (1 − e −t / τ )
(2.5)
Экспоненциальная регрессия основана на методе наименьших квадратов. Результат
обработки - получение значения временной константы τ и стационарного напряжения
термо-ЭДС Ustα , необходимого в измерениях по методу Хармана.
На рис. 2.1 представлены динамические кривые для однокаскадного ТЭО 1MC06-032-15 при
различных значениях тока.
13
45
2мА
12мА
22мА
32мА
40
35
Ua, мВ
30
25
20
15
10
5
0
0,00
2,00
4,00
6,00
8,00
10,00
t, с
Рис.2.1. Выход в стационарное состояние для однокаскадного ТЭО в зависимости от значения тока
(1MC06-032-15).
На рис. 2.2 мы видим зависимость констант времени от тока для того же ТЭО (рис. 2.1) и для
ТЭО с в три раза меньшей высотой ветви.
Константа времени, с
3,0
2,5
2,0
Высота ветви=1.5мм
1,5
Высота ветви=0.5мм
1,0
0,5
0,0
0
20
40
60
80
100
I, мA
Рис.2.2. Значения констант времени двух однокаскадных ТЭО, полученная при разных токах
(1MC06-032-15 и 1MC06-032-05).
На следующем рисунке представлено нестационарное поведение двух ТЭО отличающихся
числом ветвей для одного и того же значения тока (10мА). Отличие временных констант, повидимому, связано с небольшим варьированием отношения (2.2).
14
14
N=64, τ=2.50с
12
1MC06-008-15
1MC06-032-15
Ua, мВ
10
8
6
N=16, τ=2.26с
4
2
0
0
2
4
6
8
10
t, c
Рис.2.3. Временное поведение двух однокаскадных ТЭО с различным числом ветвей (1MC06-03215 и 1MC06-032-15).
Для многокаскадных модулей отношение (2.4) не применимо. Однако константа времени
может быть оценена по (2.4), и способ получения стационарного напряжения Зеебека тот же,
что и для однокаскадного модуля. На рис. 2.4 показан пример временного поведения
двухкаскадного ТЭО.
20
Ua, мВ
15
8мA
16мA
24мA
36мA
44мA
52mA
10
5
0
0,0
10,0
20,0
30,0
40,0
50,0
t, с
Рис.2.4. Выход на режим двухкаскадного ТЭО (2MC06-010-10) при разных значениях тока.
15
3. Технические аспекты
Внешний вид DX3065 представлен на рис. 3.1.
Fig.3.1. Внешний вид DX3065.
Корпус DX3065 выполнен из алюминиевого сплава. При измерении ТЭО он выполняет
функцию пассивного термостата. Температура корпуса измеряется с помощью цифрового
термометра с точностью не хуже 0.1°.
Измеряемый ТЭО помещается в этот корпус.
Подсоединение модуля происходит через специальные разъемы.
Упрощенная функциональная схема DX3065 показана на рис. 3.2.
TEC
R S -2 3 2
d r iv e r
R
S - 2 3
Рис.3.2. Функциональная схема DX3065
Рассмотрим основы измерений двух ключевых параметров (ACR и Uα).
3.1 Измерение сопротивления на переменном токе (ACR)
На ТЭО подается переменный ток малой амплитуды.
16
T
T
+E /2
m
t
0
Up
1
t
1
Un
2
t
1
Up
t
3
t
4
Un
2
2
Un
2 n -2
n -1
t
Up
t
2 n -1
2 n
Un
n
n
Рис. 3.3 Временная диаграмма напряжений
В отсутствии сигнала выходное напряжение усилителя равно Em/2 (рис. 3.3). Во время
измерения сопротивления выходное напряжение измеряется с помощью 12-битового ADC
прежде, чем ток Im реверсируется. Точки измерения помечены на рисунке
ti. Падения
напряжения Upi и Uni соответсвуют положительной и отрицательной полярности.
Сопротивление рассчитывается по следующей формуле
n
R=
∑ ( Upi − Uni )
(3.1)
i =1
2 ⋅ I m ⋅ AV ⋅ n
Здесь AV – коэффициент усиления усилителя; n – полное число выборок за измерение.
3.2 Измерение термо-ЭДС
При измерении полного напряжения и напряжения Uα на модуль периодически подается
малый ток (скважность равна двум).
Необходимо два последовательных сеанса измерений для получения полного напряжения и
термо-ЭДС при токе той и другой полярности.
Измерительный ток
0
Полное напряжение и термо-ЭДС ТЭО
Uα
i
U
i
Uα
i+ 1
Uα
U
i+ 1
n + 1
U
n + 1
0
Рис. 3.4 Временная диаграмма измерения напряжений
17
3.3 Программное обеспечение
Основное окно программы DX3065 дано на рис. 3.5.
Рис.3.5. Основное окно программы DX3065
4. Заключение
Представлена усовершенствованная версия Z-метра, позволяющая измерять
термоэлектрическую добротность ТЭО, константу времени, характеризующую выход на
стационарный режим, и сопротивление на переменном токе, что играет существенную роль
при сертификации, контроле качества, обнаружении неисправностей модулей и может быть
использовано как производителями, так и потребителями ТЭО.
18
Литература
1. Putley E.H. Thermoelectric and Galvanometric effects in lead selenide and telluride. Proc. Phys. Soc. B., 68, 1955,
35.
2. Harman T.C., Special techniques for measurement of thermoelectric properties. J.Appl. Phys., 29, 1959, 1373.
3. Babin V.P., Gorodetskiy S.M. Thermoelectric modules quality testing by a manufacturer. Proc. IV Int. Conf. on
Thermoelectrics, 1995, 338.
4. Buist R.J. Methodology for Testing Thermoelectric Materials and Devices.CRC Handbook of Thermoelectrics. CRC
Press LLC, 1995, 701.
5. Лукишкер Э.M, Вайнер A.Л., Сомкин M.H., Водолагин В.Ю. Термоэлектрические охладители. Москва,
“Радио и связь”, 1983, 175.
6. Дульнев G.N. Теплообмен в радиоэлектронных приборах. Ленинград, Госэнергоиздат, 1963, 288.
7. И.А.Драбкин. Характеристики термоэлектрических модулей.Санкт-Петербург, 2002,146.
8. Gromov G.,Kondratiev D., Rogov A., Yershova L. Z-meter: Easy-to-use Application and Theory. Proc. Of VI Eur.
Workshop on Thermoelectrics, Freiburg, 2001.
9. Gray P.E. The Dynamical Behaviou of Thermoelectric Devices. N.Y., London, 1960.
10. Астахова Е.И., Бабин В.П., Равич Ю.И. Расчет и измерение постоянной времени охладающего термоэлемнта
в регулярном режиме. Инженерно-физический журнал, 62, 1992, 284.
11. Кондратьев Г.M. Регулярный тепловой режим. Москва, ГИТТЛ, 1954, 408.
12. Анатычук Л.И. Термоэлементы и термоэлектрические устройства. Киев, Наукова думка, 1979,766.
19