Основная задача расчета термоэлектрогенератора состоит в определении его выходных электрических параметров и в нахождении температурных полей во всех его элементах для стационарного режима работы. При этом основными исходными данными являются термоэлектрические свойства применяемого полупроводникового материала и температурные ограничения, относящиеся к полупроводниковому термоэлектрическому материалу или к отдельным конструктивным деталям термоэлектрогенератора. Иногда эти ограничения накладываются и на элементы источника тепла.
В зависимости от того, готова ли (спроектирована) конструкция термоэлектрогенератора или предстоит ее только спроектировать, основная задача расчета распадается на две взаимопротивоположные задачи: либо по заданной геометрии конструкции и подводимому к термоэлектрогенератору потоку тепла- определить его температурный режим и выходные электрические параметры, либо по максимально допустимому температурному режиму определить его геометрию, подводимое тепло и выходные электрические параметры.
Несмотря на простоту формулировок, решить эти задачи трудно. Более того, окончательные математические зависимости, получаемые при решении каждой из этих задач, формально могут значительно различаться между собой. Это объясняется тем, что решение первой задачи связано большей частью с проведением проверочного расчета, термоэлектрогенератора, а решение второй задачи — с оптимизацией и выбором его основных параметров, исходя из тех или иных ограничений. Естественно, что эти различия носят лишь формальный характер, и то в окончательных результатах. Исходные основные уравнения для первой и второй задач имеют совершенно одинаковый вид. Это должно следовать из физической сущности самого процесса преобразования энергии.
Основное при расчете термоэлектрогенераторов — изучение распределения температур по термоэлектрическому материалу — наиболее важного термодинамического параметра в любых процессах преобразования тепловой энергии. В самом общем' случае распределение температур в преобразователе описывается уравнением теплопроводности [1]
где у — плотность; t — время; — теплоемкость вещества; Т — температура; — вектор переноса теплопроводностью. В этом уравнении первое слагаемое представляет собой изменение энтальпии вещества в единичном объеме в единицу времени; второе слагаемое — результирующий перенос тепла путем теплопроводности в единичном объеме в единицу времени. Величина обозначает удельную производительность источников тепла в единичном объеме.
В декартовой системе уравнение теплопроводности записывается в вцде
где А — коэффициент теплопроводности термоэлектрического материала.
В дальнейшем будем рассматривать только стационарные режимы термоэлектрогенератора, когда температурное поле не меняется во времени и d(ycpT)/dt=0. Нестационарные режимы связаны главным образом с процессом нагрева и расхолаживания термоэлектрогенератора, а также с переходом с одного режима работы на другой и т. д. Стационарный режим, при котором термоэлектрогенератор вырабатывает электрическую энергию, описывается уравнением Пуассона
Для решения этого уравнения необходимы соответствующие граничные условия, определяющие тепловое состояние на геометрических границах тела. Они могут быть разнообразными. Обычно различают три рода граничных условий:
1.на поверхности тела поддерживается заданное распределение температур;
2.на границах имеется заданное распределение удельных тепловых потоков;
3.поверхность тела окружает среда с заданной температурой й коэффициентом теплоотдачи.
Кроме этих трех классических случаев возможны следующие два важных рода граничных условий:
4/изучаемое тело может подвергаться тепловому воздействию- контактной теплопроводности другого тела;
5.при тепловом воздействии на поверхность тела происходит одно-, временное превращение его массы (изменение состояния), фронт которого перемещается во времени, например в процессах плавления, затвердевания и т. д.
Все граничные условия с математической точки зрения могут быть записаны в виде функциональной зависимости названных выше величин или в виде дифференциальных уравнений. В последнем случае при расчете термоэлектрогенераторов приходится решать уже систему дифференциальных уравнений.
Пр ц решении практических вопросов, и в частности расчета термо- генераторов, достаточно пользоваться нестрогими и общими решениями.
Причем последние могут быть настолько громоздкими и неудобными, что для практического их применения невозможно получить всю необходимую информацию. Поэтому сделаем допущения, которые будут влиять на точность количественных решений в приемлемых границах, оговаривая получаемые погрешности.
Во-первых, для описания распределения температур по отдельному термостолбику воспользуемся одномерным приближением. Это означает, что передачами тепла между боковыми поверхностями термостолбиков можно пренебречь. Действительно, если градиенты температур по термостолбикам в продольном направлении составляют до нескольких сот и тысяч град/см, то в поперечном направлении они не превышают нескольких град/см.
Во-вторых, все параллельные ветви находятся в одинаковых температурных условиях, т. е. будем пренебрегать паразитными перетечками тока.
Рис. 4.1. Последовательная цепь термоэлементов в термоэлектрогенераторе.
Таким образом, можно считать, что изменение температурных условий работы термоэлементов происходит только по длине последовательной цепи, например по оси х (рис. 4.1). Направление теплового потока в отдельных термоэлементах параллельно оси у.
В качестве внутренних источников тепла будем учитывать только эффекты Томсона и Джоуля, т. е.
где / — плотность тока; р — удельное сопротивление; т — коэффициент Томсона.
Тогда дифференциальное уравнение, описывающее распределение температур по отдельному термостолбику в стационарном режиме, примет вид
Здесь знак «плюс» перед вторым членом соответствует выделению тепла, а знак «минус» — поглощению. В дальнейшем изложении знак минус опускаем. Решение этого уравнения позволяет определить распределение температур и величину теплового потоца в любом сечении по всей длине термоэлемента. Возможное решение приведено в главе 3.
В случае многоэлементного термоэлектрогенератора, состоящего из т последовательно соединенных термоэлемёнтов, формально можно записать систему из т уравнений (4.2.1) и решать ее при граничных условиях на горячих и холодных спаях, заданных в конечных разностях. Однако наличие большого количества последовательно соединенных термоэлементов и постоянство плотности тока по любым сечениям любого термоэлемента позволяют сделать предельный переход от системы множества уравнений в конечных разностях к системе нескольких дифференциальных уравнений.
В этой системе одним уравнением является дифференциальное уравнение (4.2.1). В общем случае минимальное число остальных уравнений должно быть не менее двух, описывающих распределение температур в зависимости от координаты х по горячим и холодным спаям термоэлектрогенератора или, как их называют, по теплоконтактным поверхностям:
где а, b, с — параметры.
Конкретный вид последних уравнений определяется источником тепла и условиями охлаждения.
Для решения задачи по определению температурного поля в термоэлектрогенераторе проинтегрируем уравнение (4.2.1) дважды. После первого интегрирования в пределах от Тх до Т с учетом произвольной температурной зависимости параметров , т, р имеем
где
Учитывая, что
где — удельный тепловой поток в сечении холодного спая (холодной теплоконтактной поверхности), после второго интегрирования получим
Уравнеия (4.2.2) — (4.2.4) позволяют определить удельные» тепловые потоки в любом сечении термоэлементов по оси у. В частности, распределение удельного теплового потока по холодному спаю определяется интегральным выражением
где l — длина (высота) термоэлементов (по оси у), а из уравнения
где Тг и Тх — функции координаты х.
Для определения интеграла воспользуемся следующим соотношением:
где и — коэффициенты термо-э. д. с. при температуре Т и Тх соответственно.
Удельные тепловые потоки, определяемые выражениями (4.2.5) и
, не учитывают выделение и поглощение теплот Пельтье на спаях. Так как в термоэлектрогенераторах теплота Пельтье на горячем спае всегда поглощается, а на холодном спае всегда выделяется, то подводимые и отводимые от него удельные тепловые потоки будут соответственно равны
В случае наличия тепловыделения в месте коммутации термо элементов
где и — сопротивления коммутации соответственно горячего и холодного спаев, отнесенные к единице площади теплоконтактной поверхности.
Уравнения (4.2.5) — (4.2.9) совместно с граничными условиями, описывающими распределение температур или удельных тепловых потоков по оси х, позволяют определить тепловой режим термоэлектрогенератора. В более общем случае для этого требуется совместное решение уравнений (4.2,5) — (4.2.9) с уравнениями, описывающими распределение температур по источнику тепла и холодильнику, при соответствующих граничных условиях. Наиболее простые граничные условия получаются, когда заданы распределения температур горячего и холодного спаев термоэлементов (граничное условие первого рода):
Очень часто встречаются граничные условия второго рода, особенно для энергетических установок, где применяются ядерные реакторы. В этих установках задается распределение тепловыделения по теплоконтактной поверхности, т. е. распределение удельного потока по источнику тепла вдоль оси х:
К этому условию обычно добавляется требование, чтобы максимум распределения температуры по горячему спаю термоэлектрогенератора не превосходил допустимой величины. Это требование, естественно, должно соблюдаться во всех случаях, при любых граничных условиях.
Граничные условия могут быть и нелинейными, например, для космических установок с теплоотводом за счет лучистого теплообмена, который подчиняется закону Стефана—Больцмана:
где — коэффициент черноты; — постоянная Стефана—Больцмана.
Часто нагрев и охлаждение термоэлектрогенератора осуществляется с помощью применения контуров теплоносителей. В этом случае процесс нагрева и охлаждения описывается дифференциальными уравнениями:
где GT, ст, Т г.т и Gx> cXf Т хл— соответственно расход, удельная теплоемкость, температура горячего и холодного теплоносителей; Н — ширина термоэлектрогенератора по тракту теплоносителя.
Практически всегда с горячими и холодными спаями термоэлементов в непосредственном контакте находятся конструктивные элементы. В этом случае имеем граничные условия третьего рода.
В общем случае граничные условия столь разнообразны, что совершенно невозможно дать точный математический анализ работы термо- электрогенератора в общем виде. Все это усложняется тем обстоятельством, что даже уравнение (4.2.1) нельзя решить ни в аналитическом виде, ни в квадратурах.