Вышедшие номера
Математическое моделирование нестационарных режимов работы термоэлектрических источников тока
Дударев Ю.И.1, Максимов М.З.1
1Сухумский физико-технический институт АН Республики Абхазия, Сухуми, Абхазия
Поступила в редакцию: 14 апреля 1997 г.
Выставление онлайн: 20 мая 1998 г.

Термоэлектрический способ преобразования энергии, основанный на использовании эффектов Зеебека, Томсона и Пельтье, широко используется уже несколько десятилетий в различных областях науки, техники и малой энергетики. Традиционно как в холодильной технике, так и в термоэлектрических генераторах использовались стационарные режимы работы. Нестационарные режимы охлаждения и генерации тока стали изучаться значительно позднее (Л.С. Стильбансом и Н.А. Федорович [1], Е.К. Иорданишвили и его сотрудниками [2] и др.) и сразу же были выявлены некоторые преимущества таких режимов, особенно на начальной стадии функционирования, когда существенным является различие в инерционности тепловых и электрических процессов. Термоэлектрические источники питания представляют собой многослойные структуры, температурные поля в которых определяются решением системы нелинейных дифференциальных уравнений теплопроводности с соответствующими граничными условиями - условиями подвода и отвода тепла. В общем виде в одномерном приближении такая система выглядит следующим образом: (d)/(d x)gl[varkappai(T) (d Ti)/(d x)gr]+ ji(t)(dalphai(T))/(d Ti) (d TiTi)/(d x) +j2i(t)rhoi(T) =Ci(T)(d Ti)/(d t), (d T1)/(d x)- (Q1(t))/(varkappa0(0))+b1l(T-T1(t)r)= f1(t,T), x=0, (d TL)/(d x)- (Q2(t))/(varkappa0(L))+b2l(T-T2(t)r)= f2(t,T), x=L, T(x,0)=TH(x), (1) где j - плотность тока; varkappa, alpha - коэффициенты теплопроводности и термоэдс; rho - удельное электрическое сопротивление; C - объемная теплоемкость, i - номер слоя. Обычно теплофизические характеристики материалов в каждом слое можно представить в виде sigma(x,T)=sigma0(x)+Deltasigma(x,T), (Deltasigma)/(sigma0)<1. (2) Система (1) должна быть дополнена уравнениями для определения плотности тока j, которые зависят от электрической схемы соединения термоэлементов. В соответствии с методикой, изложенной в работе [3], решение этой задачи удобно представить в виде T(x,t)=T0(x,t)+U(x,t), (3) где T0(x,t) - температура системы без учета зависимости теплофизических свойств и тепловых источников от температуры и токовых эффектов. Переходя от многослойной системы к эффективной однородной на основе метода ВКБ [4], т. е. заменяя x-> Pi(x)=x0(dx)/(sqrt(a0)(x))sqrt и используя аппарат функций Грина [4], можно получить [b] T0(x,t)&=Pi(L)0TH(Pi')G(Pi,t;Pi',0)dPi' &  - t+0glgl( (Q1sqrt(a0)/()(0))sqrt Для определения U(x,t) будем иметь интегральное уравнение [b] U(x,t)&=-t+0 glPi(L)0Psi(Pi',tau) G(Pi,t;Pi',tau)dPi' [2mm]&  + sqrt(a0(0))sqrt где Psi=&-(d)/(d x) gl[Deltavarkappa(x,T)(d T)/(d x)gr]+ Delta C(x,T)(d T)/(d t) [2mm]&- j(t) (dalpha)/(d T)(d T)/(d x)T- j2(t)rho(T); [2mm]fk=&- (Deltavarkappak)/(varkappa0k) (d T)/(d x)+omegak(t,T), k=1,2, Qk, omegak - известные функции; верхний предел t+ означает, что интегрирование по tau проводится до t+varepsilon с последующим переходом к пределу при varepsilon->0; G(Pi(x),t;Pi(x'),tau) - функция Грина системы (1), для ее приближенного определения применялся метод ВКБ [b] G(Pi,t;Pi',tau)&=(1)/(2pi i) delta+iбесконечностьdelta-iбесконечность N(Pi,Pi',P) [2mm]&  x exp[p(t-tau)]dp, (7) N&~= ()/()[(p)sqrt [b] R&=sqrt(p)sqrt Pi>Pi'. При Pi<Pi' необходимо поменять местами Pi и Pi'. Уравнение (6) легко решается методом последовательных приближений, за нулевое принимают T0(x,t). [!tb] Зависимость параметра W/WmaxT от времени, 1 - расчет по предложенной методике, 2 - эксперимент, 3 - точное решение. Следует отметить, что применение асимптотических методов [5], в частности приведенного метода сращивания (ПСМ) [6], позволяет существенно упростить расчетные соотношения. При малых временах (t)/(Pi20(L))<(1)/(4) и небольшой плотности тока j<0.1(varkappa)/(alpha L) влиянием токовых эффектов на температуру можно пренебречь, что значительно сокращает вычислительную процедуру. Эффективные характеристики для многослойных систем могут быть получены с помощью правила сумм для собственных значений исходных задач [7]. На рисунке показано изменение во времени относительной электрической мощности термогенератора W/WmaxT, полученное путем расчета по предлагаемой методике (кривая 1), точным решением (кривая 3) и экспериментально (кривая 2). Экспериментальные данные получены А.А. Соколовым (СФТИ). Сравнение этих результатов показывает их удовлетворительное согласие, что позволяет рекомендовать разработанный метод для расчета и анализа характеристик конкретных термоэлектрических устройств. Математическое моделирование нестационарных режимов работы термоэлектрических генераторов показало, что для ускорения вывода термогенератора на заданный энергетический режим необходимо снижать толщину переходных слоев со стороны горячих спаев и применять в них материалы с высокими коэффициентами температуропроводности. Для поддержания длительного функционирования системы следует увеличивать ее теплоемкость, температуру или теплотворную способность источника тепла. От этих же факторов существенно зависит и амплитуда выходного электрического сигнала.
  1. Стильбанс Л.С., Федорович Н.А. // ЖТФ. 1958. Т. 28. N 3. С. 489--492
  2. Иорданишвили Е.К., Бабин В.П. Нестационарные процессы в термоэлектрических и термомагнитных системах преобразования энергии. М.: Наука, 1983. 216 с
  3. Дударев Ю.И., Максимов М.З. // ТВТ. 1988. Т. 26. N 4. С. 824--827
  4. Морс Ф.М., Фешбах Г. Методы теоретической физики. М.: ИЛ, 1958. Т. 1, 2
  5. Дударев Ю.И., Кашин А.П., Лозбин В.И., Марченко О.В. // ИФЖ. 1982. Т. 42. N 3. С. 492--493
  6. Кашин А.П., Кварацхелия Т.М., Максимов М.З., Чиковани З.Е. // ТМФ. 1989. Т. 78. N 3. С. 392--399
  7. Дударев Ю.И., Кашин А.П., Максимов М.З. // ИФЖ. 1985. Т. 48. N 2. С. 333--334

Подсчитывается количество просмотров абстрактов ("html" на диаграммах) и полных версий статей ("pdf"). Просмотры с одинаковых IP-адресов засчитываются, если происходят с интервалом не менее 2-х часов.

Дата начала обработки статистических данных - 27 января 2016 г.