СРАВНЕНИЕ РАЗЛИЧНЫХ ПОДХОДОВ К ОПТИМИЗАЦИИ ОДНОКАСКАДНЫХ ТЕРМОЭЛЕКТРИЧЕСКИХ МОДУЛЕЙ Драбкин И.А.1, Ершова Л.Б.2 1 Институт Химических Проблем Микроэлектроники, Москва, Россия 2 ЗАО «РМТ», Москва, Россия тел: +7-495-132-6817 факс: +7-495-132-5870 e-mail: [email protected] http://www.rmtltd.ru Долгое время все расчеты термоэлектрических охладителей проводились с использованием простейшей модели, предполагающей независимость термоэлектрических параметров от температуры (см., например, [1]). Такая модель рассматривает в качестве первичного элемента термоэлектрическую ветвь и дает простые алгебраические выражения для нахождения оптимальных режимов работы. Достоинством этой модели является простота и наглядность. Небольшие усложнения этой модели, связанные с введением эффективных термоэлектрических параметров [2,3] (далее будем ссылаться как на метод эффективных параметров, метод ЭП), позволяют точно учесть тепловые балансы на концах ветви, сохраняя простоту и наглядность в качестве основных преимуществ. Математически строгий учет температурных зависимостей был осуществлен с использованием метода максимума Понтрягина (метода МП) в работах [4,5]. Этот метод дает математически формализованные рецепты нахождения оптимальных решений, из-за чего связь между термоэлектрическими параметрами и оптимальными режимами выглядит весьма опосредованной. Сильной стороной этого метода является то, что в качестве первичного звена может выступать не ветвь, а термоэлемент. В данной работе проводится сравнение результатов расчетов термоэлектрических охладителей в рамках модели эффективных параметров и метода максимума и показывается, что они дают близкие результаты. Работа ветви термоэлемента обычно описывается уравнением теплопроводности с температурно-зависимыми коэффициентами в одномерном приближении. Для численного решения это уравнение второго порядка удобно привести к системе двух уравнений первого порядка, выбирая в качестве одного независимого переменного температуру T x , а в качестве второго – величину dT x T jT x , где T - теплопроводность, dx dT x T - термоэдс, j - плотность тока, и двойной знак перед членом с dx связан с выбором направления оси x . Такой выбор в качестве второго переменного q удобен тем, что на концах ветви его величина оказывается q , пропорциональную T 378 пропорциональной либо холодопроизводительности либо Q0 , теплопроизводительности Q в зависимости от того, теплопоглощающий или тепловыделяющий конец ветви рассматривается. Будем считать далее, что ветвь имеет единичное сечение и длину, горячему концу ветви соответствует начало координат, тогда холодному – единица. Плотность тока, независимо от его направления, будем считать j 0 , а термоэдс относительно материала контакта 0 , независимо от ее истинного знака. Тогда для ветви термоэлемента систему уравнений первого порядка, к которой приводится уравнение теплопроводности, можно записать в виде: dT j jT q , dx , (1) dq j j j 2 q T, dx dT jT x где - электропроводность, а q dx . Систему (1) надо решать при j граничных условиях q0 Q / sj , q1 Q0 / sj , где s - сечение ветви. В случае задачи Коши необходимо задать T и q на одном из концов ветви, и во всех языках программирования можно найти стандартные программы для решения такого типа задач (например, для Turbo Pascal он приведен в книге [6]). В термоэлектрических задачах обычно задаются температуры горячего Th и холодного Tc концов ветви. Для сведения такой задачи к задаче Коши обычно применяется так называемый метод «стрельбы», для которого также имеются стандартные программы [6]. Интерполяционные программы, позволяющие построить температурные зависимости термоэлектрических параметров по экспериментальным точкам также можно найти в любом языке программирования. Поэтому само по себе решение уравнения теплопроводности для термоэлектрической ветви не представляет значительных трудностей. Основная проблема термоэлектричества заключается в оптимизации полученных решений. Обычно необходимо оптимизировать ток, т.е. найти j opt , такой чтобы тепловой коэффициент : Q q0 , Q0 q1 (2) достигал своего минимума min . Наиболее математически последовательным методом оптимизации (1) является метод МП [7], в котором вводятся сопряженные к T и q переменные 1 и 2 , а затем строится Гамильтониан H : 379 dT dq , 2 dx dx из которого находятся уравнения движения для 1 и 2 : H 1 (3) d 1 H d 2 H , (4) dx T dx q Ищется решение, которое минимизирует (2), и, исходя из этого, строятся граничные условия для сопряженных переменных из условия трансверсальности. При необходимости в них учитываются дополнительные ограничения на работу ветви. В этом случае количество уравнений, необходимое для описания работы ветви удваивается, что, впрочем, не создает особых трудностей для их решения. Приравнивая Гамильтониан к нулю в любой () x можно найти оптимальное управление. При решении термоэлектрических задач, связанных с охлаждением, для заданных ветвей оптимизация сводится к определению одного параметра j opt , который дается для ветви выражением: 1 jopt где Z x 2T 0 q T d 1 2 1 2 0 1 ZT dT q T dx (5) T x 2 T x - термоэлектрическая эффективность. T x Важно отметить, что выражение несколько отлично от приводимого в работе [4]. Так как с помощью метода МП удобно оптимизировать не ветвь, а термоэлемент, то система (1) удваивается, при этом вводятся переменные T и q для каждого типа проводимости. Соответственно, в Гамильтониане (3) проводится суммирование по типам проводимости, количество сопряженных переменных 1 и 2 также удваивается в согласии с типом проводимости. Уравнения движения (4) пишутся отдельно для каждого типа, а в (5) суммирование по типам проводимости проводится и в числителе, и знаменателе. Перейдем теперь к методу ЭП. В случае температурно-независимых термоэлектрических переменных нет необходимости находить распределение температур вдоль ветви термоэлемента, достаточно рассмотреть простые алгебраические уравнения теплового баланса на концах ветви, из которых очень просто находится j opt . В случае температурно-зависимых термоэлектрических переменных было показано [2,3], что для уравнений теплового баланса на концах ветви можно получить аналогичные выражения: 380 1 2 j c eff T q1, 2 (6) 1 hTh j j 2 Rh eff T q0 , 2 где T Th Tc , а h , c , h , c для ветви единичной длины и сечения определяются следующим образом: 1 eff (7) 1 dx 0 ( T ) x cTc j c ( Tc ) h ( Th ) 1 Ty Tc 0 1 Ty Th 0 d ( T y ) dT dT dy d ( T y ) dT dT dy 1 1 dx ( Tx ) (8) y dx ( Tx ) (9) dy 0 dx 1 c 2 0 ( T y )dy y dy y ( Tx ) (10) dx (11) ( Tx ) Эти выражения отличаются от предложенных ранее Бурштейном [8] тем, что они являются точными, а не ограничиваются случаем постоянства теплового потока вдоль ветви. Из пяти величин (7)-(11) независимыми являются только три, т.к. c и h связаны соотношением: 1 y h 2 0 ( T y )dy 0 hTh cTc T T h T dT c T (12) а c и h связаны как: 1 h c 2 20 ( Tx )dx (13) Для вычисления этих параметров предварительно необходимо решить систему (1), а затем уже вычислить интегралы, входящие в (7) – (11). Следует подчеркнуть, что система (6), описывающая поведение ветви, является точной. Теплота Томсона проявляется в том, что значения эффективных параметров c и h отличаются друг от друга, а тот факт, что тепловые потоки от тепла Джоуля делятся не поровну между концами ветви, приводит к различию в величинах c и h . Оптимальный ток для ветви с заданными Th и Tc определяется как: j opt T M eff 1 где M eff равно: 381 (14) M eff 1 ZTeff , Z 2 (15) а эффективное значение температуры: Tc Teff c h Th h c 2 Tc c T c T h Th h 2 2 (16) Уравнения теплового баланса для термоэлемента аналогичны уравнениям для ветви: 1 cn cp Tc j j 2 cn cp effn effp T qn 1 q p 1, 2 (17) 1 hn hp Th j j 2 hn hp effn effp T qn 0 q p 0 , 2 где индекс n,p относится к типу проводимости ветви, а cn , cp , cn , cp , effn , effp описывается формулами аналогичными (7)-(11) с учетом типа проводимости. Соответственно , , являются суммами средних величин для n- и p- ветвей. Выражения (14)-(16) также справедливы для термоэлемента с вышеописанной заменой. В дальнейших вычислениях min , q0 , q1 находим непосредственно из решения уравнений (1) или (17), а не вычислялись их по каким-либо приближенным формулам, т.к. зная jopt , при решении уравнения теплопроводности можно получить их значения с желаемой степенью точности. Для численных расчетов температурные зависимости термоэлектрических параметров материалов твердых растворов на основе халькогенидов висмута-сурьмы аппроксимировались по экспериментальным точкам полиномами третьей степени по температуре. Для того, чтобы иметь возможность строить кривые для любой концентрации носителей заряда, коэффициенты разложения по температуре для различных концентрации также аппроксимировались полиномами третьей степени по термоэдс материала при комнатной температуре - 300 . Таким образом, для любой величины 300 в пределах 200 – 290 мкВ/К были построены температурные зависимости термоэлектрических параметров для температур 150 – 330 К. Все расчеты охладителей проводились с аппроксимироваными таким образом кривыми. Расчеты как методом МП, так и методом ЭП проводились методом последовательных приближений. Т.е. для выбранного из каких-либо 0 решалась система (1) или (17), соображений тока в нулевом приближении jopt вычислялось новое значение оптимального тока по (5) для метода Понтрягина и по (14) для метода эффективных значений. Далее с этим новым значением тока вновь решались уравнения и вновь вычислялись значения оптимального тока. 382 Вычисления проводились до тех пор, пока различие между последовательными приближениями отличались не более чем на 0,2%, так как это дает различие в холодильном коэффициенте, выходящее за рамки точности определения 0 термоэлектрических параметров. В качестве нулевого приближения для jopt удобно брать значение тока, вычисляемого в соответствии с (14), но среднее значение любого термоэлектрического параметра p вычисляется как: pTh pTc 2 а эффективное значение температуры в нулевом приближении: T T Teff h c 2 p (18) (19) Результаты расчетов приведены в таблице. Таблица Сравнение результатов расчета оптимальных режимов работы термоэлемента методом Понтрягина и методом эффективных параметров. Расчет по методу Расчет по методу 300 максимума эффективных 300 p , n , коэффициентов Tc , К Th , мкВ/К мкВ/К jopt , jopt , К 250 220 250 220 300 250 300 250 210 230 210 230 А/см2 31,126 21,88 26,986 18,751 210 230 230 250 4,185 3,012 4,156 3,061 А/см2 31,263 22,516 28,698 20,366 4,165 3,014 4,155 3,052 Из приведенной таблицы видно, что расчет и тем и другим способом дает для весьма близкие значения, отличающиеся друг от друга на доли процента. В то же время значения оптимальных токов заметно (почти на 10%) различны. Это связано с тем, что для рассматриваемого диапазона температур минимум в зависимости от тока очень пологий, и поэтому трудно требовать большей точности, учитывая многостадийный способ получения оптимальных токов. При практической разработке термоэлектрических модулей это не должно приводить к каким либо существенным отрицательным последствиям, т.к. отличие расчетных токов скажется лишь на числе ветвей модуля из-за различия в величинах холодопроизводительности термоэлемента. Вопрос точности вычислений тепловых коэффициентов весьма существенен при практических расчетах термоэлектрических модулей, но переоценивать точность расчетов не имеет смысла, т.к. сами по себе термоэлектрические ветви не идентичны по своим свойствам, поэтому для 383 каждой ветви (термоэлемента) имеется свой оптимальный ток. Благодаря пологой зависимости j это расхождение по величинам оптимальных токов не приводит к большому различию в величине , и поэтому достигаемой точности вычислений вполне достаточно. Так как в расчетах модулей обычно выступают температурные зависимости термоэлектрических параметров, аппроксимируемые полиномами, то можно такие же полиномиальные аппроксимации использовать и в отношении cn , cp , cn , cp , effn , effp в оптимальных режимах. Тогда расчеты модулей станут более наглядными и будет легче проследить связь между исходными свойствами материала и конечным результатом. ЛИТЕРАТУРА 1. Вайнер А.Л., ред. Термоэлектрические охладители, Радио и связь, Москва, 1983, стр. 127-132 2. Drabkin I., Dashevsky Z. Proc. of fifth European workshop on thermoelectrics, Pardubice, 1999, p. 154. 3. Драбкин И.А., Дашевский З.М., Термоэлектрики и их применение, С.Петербург,2000 г, стр. 292-297. 4. Семенюк В.А. Каскадный термоэлектрический охладитель как объект оптимального управления, Инженерно-физический журнал, 1984, т. 47, № 6, стр. 977-978. 5. Анатычук А.И., Семенюк В.А. Оптимальное управление свойствами термоэлектрических материалов и приборов, Черновцы, «Прут», 1992, стр. 159-177. 6. Немнюгин С.А. Turbo Pascal. Программирование на языке высокого уровня, «Питер Принт», 2005, стр. 544. 7. Понтрягин Л.С., Болтянский В.Г., Гамкрелидзе Р.В., Мищенко Е.Ф. Математическая теория оптимальных процессов. М, Наука, 1976, стр. 13-132. 8. Бурштейн А.И. Физические основы расчета полупроводниковых термоэлектрических устройств. М, Физматгиз, 1962, стр. 64-78. 384