Проблемы машиностроения и надежности машин, 2023, № 2, стр. 103-112

Динамические тепловые процессы, происходящие в устройствах, обогреваемых на холоде

А. А. Шульженко 1*, М. Б. Модестов 1

1 Институт машиноведения им. А.А. Благонравова РАН
Москва, Россия

* E-mail: aa_shulzhenko.01@mail.ru

Поступила в редакцию 01.09.2022
После доработки 29.11.2022
Принята к публикации 20.12.2022

Полный текст (PDF)

Аннотация

В ходе работы были проведены математические эксперименты по изучению динамических тепловых процессов, происходящих при обогреве технических устройств, выносимых на холод. Новизной является получение объемных представлений о тепловых процессах, протекающих в этих устройствах во времени: обогреве их на уровне нижней границы рабочего диапазона температур, использовании ранее накопленной тепловой энергии в теплом помещении и при повышении ее за счет дополнительного нагрева ткаными электронагревателями в период нахождения устройств в теплом помещении. Настоящая статья может интересовать разработчиков нагревательных систем, в которых применяются автономные источники электроэнергии.

Ключевые слова: тканый электронагреватель, обогреваемые устройства, обогрев аккумулятора, автономные источники электроэнергии, остаточное тепло, математическое моделирование

Многие технические устройства, работающие на открытом воздухе, испытывают воздействие холода. Например, кино-, фотоаппаратура, измерительные устройства, которые для производства работ выносят на холод из теплых помещений, аккумуляторы, находящиеся в двигательных отсеках транспортных средств и т.д. Работоспособность этих устройств, зависит от обеспечения температурной стабильности, т.е. от соблюдения температурного рабочего диапазона. Эксплуатация устройств на холоде без дополнительных мер по теплозащите приводит к росту погрешностей при измерениях, к затруднениям при эксплуатации, или к полной потере работоспособности [15].

Только применяя тепловые способы защиты можно добиться нормальной эксплуатации устройств. Обычно применяются, как пассивные (дополнительные слои тепловой изоляции), так и активные (дополнительные нагревательные устройства) средства тепловой защиты объекта.

Что касается использования только пассивных средств тепловой защиты, то их применение приведет к некоторому продлению поддержания условий, обеспечивающих нормальную работу технических средств, однако, их возможности в плане тепловой защиты, особенно в регионах Севера, ограничены [3, 69].

Целью настоящей статьи является рассмотрение динамических тепловых процессов, происходящих в устройстве при его нагреве на холоде, с применением не только активных, но и активных в сочетании с пассивными средствами тепловой защиты, при использовании ранее накопленного тепла. Изучение динамических тепловых процессов, происходящих при таких условиях, позволит не только получить информацию об их протекании в устройствах, но и дать рекомендации по эксплуатации обогреваемых технических устройств в условиях холода.

В качестве активных средств тепловой защиты будем применять тканые электронагреватели (далее нагреватели), а в качестве источников электроэнергии для них – автономные источники. Как правило, автономные источники обладают ограниченными запасами электроэнергии, поэтому очень важно рационально расходовать электроэнергию для обогрева устройств. Для снижения энергетических затрат будем рассматривать и обогрев, обеспечивающий работу обогреваемым устройствам в области нижних границ температур их рабочего диапазона. Такой режим работы является особенно актуальным в регионах Севера, в районах со стабильно низкими зимними температурами, значительно ниже нижней границы температур рабочего диапазона.

Для обеспечения возможности сопоставления, будем считать, что при всех экспериментах используется одно и то же устройство. Также, во всех случаях нагрев будет проводиться одним и тем же нагревателем, обладающим неизменными мощностными параметрами. Так как толщина нагревателя крайне мала, то будем считать, что нагреватель является внутренним источником тепла, расположенным по поверхности обогреваемого устройства.

Пусть, с момента включения нагреватель работает непрерывно.

Будем изменять моменты включения нагревателя относительно времени выноса устройства из тепла на холод и уровни температуры на устройстве перед выносом.

Пусть число БИО, характеризующее соотношение внутреннего и внешнего теплового сопротивления устройства Bi < 0.1. Сначала рассмотрим уже известные данные о динамике тепловых процессов, происходящих при нагреве и охлаждении таких объектов. Изменение температуры будем отслеживать в одной точке, например, на поверхности устройства.

Аналитическое выражение, описывающее общий характер изменения температуры предварительно нагретого тела при его охлаждении, можно записать как

(1)
$U = {{U}_{1}}{{K}_{1}}\exp ( - {{k}_{1}}t),$
где U = T – Tср; T – текущее значение температуры; Tср – температура окружающей среды; U1 = T1Tср – стационарный температурный напор в теплом помещении; T1 – температура теплого помещения, до которой нагревается объект; t – относительное время протекания теплового процесса; ${{k}_{1}}$ – коэффициент; K1 – коэффициент, отражающий особенности структуры устройства.

Характер изменения температуры устройства при его нагреве можно представить в виде

(2)
$U = {{U}_{2}}{{K}_{2}}[1 - \exp ( - {{k}_{2}}t)],$
где U2 = T2Tср – стационарный температурный напор нагревателя; T2 – максимальная температура нагрева нагревателя; ${{k}_{2}}$ – коэффициент; K2 – коэффициент, отражающий особенности структуры нагревателя.

Предположим, что устройство, равномерно разогретое до температуры помещения, в котором оно находилось, вынесли на холод и, одновременно, включили нагреватель. Тепловые процессы во времени, используя (1), (2), можно представить как

(3)
$U = {{U}_{1}}{{K}_{1}}\exp ( - {{k}_{1}}t) + {{U}_{2}}{{K}_{2}}[1 - \exp ( - {{k}_{2}}t)].$

Выражения (1)–(3) были выведены и многократно использовались в широком круге работ [69].

Рассмотрим еще один случай, который можно описать с помощью выражений (1)–(3), когда нагреватель включался в период нахождения устройства в помещении. Устройство, нагретое теплом помещения и дополнительно нагревателем, в момент времени t1 выносится на холод. Этот процесс можно представить в виде

(4)
$\left\{ {\begin{array}{*{20}{l}} {U = {{U}_{1}} + {{U}_{2}}{{K}_{2}}[1 - \exp ( - {{k}_{2}}t)]\quad при\quad 0 \leqslant t \leqslant {{t}_{1}},} \\ \begin{gathered} U = \left\{ {{{U}_{1}} + {{U}_{2}}{{K}_{2}}[1 - \exp ( - {{k}_{2}}{{t}_{1}})]} \right\}\exp ( - {{k}_{1}}t) + \hfill \\ + \;{{U}_{2}}{{K}_{2}}[1 - \exp ( - {{k}_{2}}t)]\quad при\quad {{t}_{1}} \leqslant t \leqslant {{t}_{3}}, \hfill \\ \end{gathered} \end{array}} \right.$
где t1 – относительное время выноса объекта из помещения на холод; t3 = 1 – полное относительное время протекания теплового процесса.

На основании выбранных тепловых параметров устройства, нагревателя и данных внешней среды построим графики температурных распределений с использованием выражений (1)–(4), при разном времени включения нагревателя и разных температурах устройства в момент их выноса на холод. Проведем их анализ.

На рис. 1 показано изменение температуры на поверхности устройства во времени.

Рис. 1.

График распределения температур во времени на устройстве при различных условиях: J1–J9 – линии графика 19 соответственно.

На рис. 1, линия 1, происходит охлаждение устройства с момента его выноса из теплого помещения, t1 – условие J1. Снижение от температуры равной температуре помещения U1, до уровня температуры внешней среды происходит по экспоненте (1).

Условие J2 отражает нагрев от Tср до T2 в области поверхности тела, при включении нагревателя в момент выноса устройства из теплого помещения (рис. 1, линия 2). U2 – стационарный температурный напор нагревателя (2).

Условие J3, описанное выражением (3), отображает результирующее изменение температуры при одновременном воздействии тепловых условий J1 и J2 [610]. В нашем случае наблюдается сначала экспоненциальное снижение температуры до уровня, лежащего ниже уровня нижней границы температур рабочего диапазона, а затем повышение суммарной температуры (рис. 1, линия 3) до уровня стационарного температурного напора нагревателя.

В следующем эксперименте J4 – уровень температуры нагретого устройства в момент начала охлаждения t2 соответствует нижней границе рабочего диапазона устройства U3 (рис. 1, линия 4). Далее оно охлаждается до температуры среды аналогично (3).

Параметры нагревателя остаются такими же, как и в первом случае. Однако, включение нагревателя, расположенного на поверхности устройства, охлажденного до температуры внешней среды, происходит в момент t2, когда температура устройства соответствует температуре нижней границы температур рабочего диапазона устройства – условие J5 (рис. 1, линия 5).

Суммарное изменение температуры устройства, условие J6, при одновременном действии охлаждения по условию J4 и нагрева нагревателем, условие J5, продемонстрировано на рис. 1, линия 6.

Сопоставление результатов, приведенных на рис. 1 линии 3 и 6, показывает, что при включении нагревателя сразу после выноса устройства из теплого помещения (условие U3) наблюдается снижение температуры, но оно незначительно на небольшое время выходит за пределы температур рабочего диапазона устройства. Этому способствует использование остаточного тепла, накопленного устройством. Обеспечение более длительного нахождения устройства в диапазоне рабочих температур, в первом случае, обеспечивает более правильную работу устройства.

Оценим значение остаточного тепла. Расчеты показывают, что накопленная остаточная тепловая энергия в ряде случаев позволяет сэкономить значительную энергию используемого для обогрева аккумулятора, до ~12% от его расчетной емкости [1]. Так как процесс охлаждения устройства при выносе его из тепла во внешнюю среду достаточно быстрый, то рационально включать нагреватель сразу, после выноса его из теплого помещения, чтобы максимально использовать накопленную в помещении (или отсеке) тепловую энергию. Использование устройства через длительный срок после выноса его на холод приводит к его остыванию, потере накопленной энергии, и при дальнейшем его нагреве необходимо затратить дополнительно энергию автономного источника электропитания нагревателя для достижения уровня температур на устройстве, соответствующего уровню температур, обеспечивавшегося остаточным теплом.

Теперь рассмотрим случай, когда нагреватель включается, когда устройство находится еще в теплом помещении (условие J7). Работа нагревателя показана на рис. 1, линия 7.

Изменение температуры согласно условию J8: сначала – поддержание температуры на постоянном уровне в теплом помещении, а затем – экспоненциальное охлаждение (рис. 1, линия 8).

Суммарное изменение температуры при одновременном воздействии тепловых условий (рис. 1, линии 7 и 8) – условие J9, описанное в формуле (4), показано на рис. 1, линия 9). Условие J9 демонстрирует, что если нагреватель включается, еще до выноса устройства из тепла, то это обеспечивает устройству дополнительный нагрев относительно температуры помещения. Максимальное дополнительное тепло будет достигаться при выходе нагревателя на уровень стационарного температурного напора в период нахождения в обогреваемом помещении.

Таким образом, на обогреваемом устройстве при выбранных условиях может обеспечиваться тепловой процесс по сравнению с условиями J3 (представлены на рис. 1, линия 3) с изменениями температуры, не выходящими за пределы температур рабочего диапазона в течение всего времени его протекания (рис. 1, линия 9).

Как отмечалось, процессы, показанные на рис. 1, соответствуют нагреву и охлаждению устройства, но только в одной точке. В нашем случае выбрана точка на его поверхности. Реально происходит более сложная картина. Для ее демонстрации необходимо иметь более полное представление о тепловых процессах, происходящих в устройстве. Использовать аналитические методы для описания тепловых процессов даже для упрощенного вида устройства и внешних условий становится затруднительным [8, 9].

Поэтому, для получения более реальных данных о процессах, происходящих в устройстве при его нагреве и охлаждении, в соответствии с поставленной целью, было проведено математическое моделирование методом конечных разностей объемных тепловых процессов, протекавших во времени.

При формулировании математической задачи учитывались следующие требования и допущения. Во всех случаях рассматривалось одно и то же устройство в виде параллелепипеда. Тепловые свойства устройства по его объему при различных условиях воздействия предполагались постоянными. Исследовались тепловые процессы, происходившие при различных условиях выноса устройства из тепла, при различных моментах включения электронагревателей, и, для большей полноты представления, продемонстрированы варианты с использованием и при отсутствии пассивных средств тепловой защиты.

В качестве электронагревателя применялся тканый электронагреватель со свойственными ему теплофизическими параметрами и очень малой толщиной. Тепловые свойства материала нагревателя близки к тепловым свойствам обогреваемого устройства. Нагреватель был размещен по всей поверхности объекта обогрева. Для описания объемных тепловых процессов, происходящих в устройстве и одновременно отражающих работу электронагревателя, использовалось трехмерное уравнение Фурье с внутренними источниками.

При моделировании температурные режимы на границах обогреваемого объекта определялись с помощью граничных условий третьего рода.

Под температурой понималась разность U между текущими значениями температуры T и температурой окружающей среды Tср = –50°С, U = T – Tср. Будем считать, что

${{U}_{3}} = {{T}_{{нгр}}} - {{T}_{{ср}}},$
где Tнгр – температура нижней границы рабочего диапазона устройства.

Будем считать, что во всех экспериментах Tнгр= –10°С. Тогда разность температур на уровне нижней границы среды будет равна U3 = 40°С.

Для удобства рассмотрения настроим максимальные значения на шкалах температур, на рис. 2а–5а, на уровень U = 100°С, а на всех остальных, на рис. 2б, в, г–рис. 5б, в, г, на уровень нижней границы рабочего диапазона температур устройства, U = 40°С.

Рис. 2.

Температурное распределение по объему обогреваемого устройства в разное время при его начальном нагреве до температуры помещения: (а) – в момент выноса устройства из теплого помещения; в моменты времени, соответствующие значениям координаты t рис. 1: (б) – 0.4; (в) – 0.5; (г) – 0.6.

Сформулированная математическая задача имела вид

(5)
$\left\{ {\begin{array}{*{20}{l}} {{{U}_{{x1}}} = {{T}_{{x1}}} - {{T}_{{ср}}} = 0\quad при\;t = 0,\;x{{1}_{0}} \leqslant x1 \leqslant x{{1}_{1}},} \\ {{{U}_{{x2}}} = {{T}_{{x2}}} - {{T}_{{ср}}} = 0\quad при\;t = 0,\;x{{2}_{0}} \leqslant x2 \leqslant x{{2}_{1}},} \\ {{{U}_{{x3}}} = {{T}_{{x3}}} - {{T}_{{ср}}} = 0\quad при\;t = 0,\;x{{3}_{0}} \leqslant x3 \leqslant x{{3}_{1}},} \\ { - \lambda \frac{{\partial {{U}_{{x1}}}}}{{\partial x1}} = - \alpha ({{T}_{{x1}}} - {{T}_{{sr}}})\quad при\;t > 0,\;x1 = x{{1}_{{0\;}}},\;x1 = x{{1}_{1}},} \\ { - \lambda \frac{{\partial {{U}_{{x2}}}}}{{\partial x2}} = - \alpha ({{T}_{{x2}}} - {{T}_{{sr}}})\quad при\;t > 0,\;x2 = x{{2}_{0}},\;x2 = x{{2}_{1}},} \\ { - \lambda \frac{{\partial {{U}_{{x3}}}}}{{\partial x3}} = - \alpha ({{T}_{{x3}}} - {{T}_{{sr}}})\quad при\;t > 0,\;x3 = x{{3}_{0}},\;x3 = x{{3}_{1}},} \\ \begin{gathered} \rho с\frac{{\partial U}}{{\partial t}} = \lambda \left( {\frac{{{{\partial }^{2}}{{U}_{{x1}}}}}{{\partial {{{(x1)}}^{2}}}} + \frac{{{{\partial }^{2}}{{U}_{{x2}}}}}{{\partial {{{(x2)}}^{2}}}} + \frac{{{{\partial }^{2}}{{U}_{{x3}}}}}{{\partial {{{(x3)}}^{2}}}}} \right) + {{q}_{{v}}} \hfill \\ при\;t > 0,\;x{{1}_{0}} < x1 < x{{1}_{1}},\;x{{2}_{0}} < x2 < x{{2}_{1}},x{{3}_{0}} < x3 < x{{3}_{1}}. \hfill \\ \end{gathered} \end{array}} \right.$

На основании поставленной математической задачи (5) была разработана программа [10], с помощью которой получены результаты моделирования тепловых процессов, происходящих в устройстве при его нагреве на холоде во времени и в объеме.

Следует отметить, что при численных экспериментах использовались фиксированные временные точки для обеспечения возможности сопоставления температурных изменений в устройстве при различных условиях теплового воздействия.

На рис. 2 показан тепловой процесс, происходящий при нагреве устройства нагревателем, расположенным по всей его поверхности, при одновременном остывании этого устройства под воздействием внешней среды с уровня температуры, соответствующего оставшемуся теплу от нахождения его в теплом помещении. Нагреватель включали одновременно с выносом объекта из теплого помещения. Приведенные характерные точки, отражающие тепловое состояние устройства, примерно соответствовали точкам рис. 1, со значениями временной координаты t: 0.1, 0.4, 0.5, 0.6. Эти же точки использовались при всех экспериментах с объемными тепловыми процессами и соответствовали временным значениям, показанным на рис. 2–5, графиках: (а), (б), (в) и (г) соответственно.

Рис. 3.

Температурное распределение по объему обогреваемого устройства в разное время при его начальной температуре, равной температуре нижней границы рабочего диапазона, Uнгр= 40°С: (а) – в момент времени начала охлаждения и одновременно включения нагревателя; в моменты времени, соответствующие значениям координаты t рис. 1: (б) – 0.4; (в) – 0.5; (г) – 0.6.

Рис. 4.

Температурное распределение по объему обогреваемого устройства в различное время при включении нагревателя еще в теплом помещении: (а) – в момент выноса устройства из теплого помещения; в моменты времени, соответствующие значениям координаты t рис. 1: (б) – 0.4; (в) – 0.5; (г) – 0.6.

Рис. 5.

Температурное распределение по объему устройства в разное время при использовании теплоизоляционного слоя на поверхности нагревателя: (а) – в момент выноса устройства из теплого помещения; в моменты времени, соответствующие значениям координаты t рис. 1: (б) – 0.4; (в) – 0.5; (г) – 0.6.

На рис. 2а показано температурное распределение по устройству в момент его выноса из теплого помещения. В этот момент температура равномерно распределена по всему объему и равна температуре устройства в помещении, U1 = 80°C. Она существенно превышает температуру нижней границы рабочей зоны. То есть, устройство имеет накопленное тепло. Через некоторое время, в ходе охлаждения температура большей части объема устройства выходит практически полностью за пределы рабочей зоны (рис. 2б). В этот период нагреватель еще не вышел на максимальный режим нагрева. Однако, уже на следующем рис. 2в, зафиксированном спустя еще некоторое время, видно, что устройство опять почти полностью находится в пределах зоны рабочих температур. На рис. 2г видно, что распределение температур по объему объекта еще больше выравнивается.

На рис. 3 показан тепловой процесс, также происходящий при нагреве устройства нагревателем. В момент выноса устройства во внешнюю среду оно имело равномерно распределенную температуру, соответствующую температуре нижней границы рабочей зоны, U3 = 40°С (рис. 3а). В этот момент включался нагреватель, находящийся на его поверхности.

На рис. 3б, в можно наблюдать, что температурное распределение по устройству значительно выходит за границы рабочей зоны. И только на рис. 3г температурное распределение выравнивается и почти полностью соответствует рабочему диапазону температур. То есть, данный результат демонстрирует ухудшение в распределении температур по объему устройства и во времени.

Температурные изменения, показанные на рис. 2, 3, происходящие в устройстве, аналогичны характеру температурных изменений, представленных на рис. 1, линии 3 и 6. Меньший выход температур за границы рабочего диапазона зафиксирован на рис. 2. Эти изменения температур аналогичны изменениям, продемонстрированным на рис. 1, линия 3. Основными причинами, приводящими к снижению уровня и длительности выхода температур за границу рабочего диапазона, является включение нагревателя в момент выноса устройства из теплого помещения, что позволяет использовать тепловую энергию, запасенную в теплом помещении.

Теперь рассмотрим, что произойдет, если нагреватель включить при нахождении устройства еще в помещении и дополнительно, равномерно его нагреть.

На рис. 4 отчетливо видно, что степень выхода температур за границы рабочего диапазона, по сравнению с тепловыми процессами, показанными на рис. 2, стала еще меньше. Аналогичная картина наблюдается и при сравнении рис. 1, линия 9 и линия 3. Отличие в распределении температур по сравнению с рис. 2 также объясняется ростом температуры за счет дополнительного к температуре помещения нагреву нагревателем устройства во время пребывания в помещении, что обеспечивает дополнительный подъем температуры в период нахождения устройства вне помещения (рис. 4б, в, г). Температура устройства в момент его выноса равна U = 95°С.

Следует отметить, что процессы, показанные на рис. 4, имеют распределение температур по объему устройства, незначительно выходящее за пределы границ температур рабочего диапазона (рис. 4б), хотя на рис. 1, линия 9, вся кривая, отслеживающая тепловое состояние в одной точке на поверхности устройства, полностью находится в пределах границ температур рабочего диапазона. То есть, более полное представление о степени нахождения или ненахождения устройства в области рабочих температур при протекании переходных процессов дает объемное тепловое поле.

Рассмотрим распределение температур по объему устройства во времени при одновременном использовании активных и пассивных средств тепловой защиты при нахождении его на холоде. Начальная температура устройства в момент его выноса из помещения равнялась, как и в предыдущем случае, U = 95°C.

Сопоставляя тепловые процессы, показанные на рис. 4, 5, хорошо видно, что на рис. 5, наблюдается улучшение результатов обогрева устройства за счет использования теплоизоляционного слоя. В рассматриваемом случае, активный обогрев в сочетании с применением теплоизоляционного слоя обеспечивает практически в течение всего времени работы в условиях холода распределение температур по объему устройства, не выходящее за пределы нижних температур рабочего диапазона (рис. 5б).

Даже небольшой дополнительный теплоизоляционный слой позволяет добиться существенного снижения степени охлаждения, приводящего за одно и то же время, к более равномерному распределению тепла во времени и более высокому уровню нагрева устройства при тех же энергетических затратах со стороны нагревательной системы.

Выводы. Используя режим поддержания температуры устройства в области нижней границы рабочего диапазона температур можно добиться снижения энергетических затрат автономного источника электропитания на обогрев устройства в условиях сильного холода. Однако необходимо учитывать, что при таком нагреве присутствует возможность выхода устройства за границу температур рабочего диапазона.

Снижения периода выхода за границы температур рабочего диапазона во время протекания динамических тепловых процессов при обогреве устройства на холоде можно добиться при использовании остаточного тепла при выносе устройства из теплого помещения.

Этот эффект можно усилить, если включать нагреватели еще в момент нахождения устройства в теплом помещении, т.к., в этом случае, дополнительно к теплу, накопленному в помещении, добавляется энергия от нагревателей.

Полученные в ходе моделирования данные о тепловых процессах в устройстве показывают, что в реальных условиях для настройки тепловой системы: нагреватель–объект–внешняя среда, нужно использовать методы измерения, позволяющие учесть температурные изменения, происходящие в нем и по объему, и во времени.

Применение слоев пассивной теплозащиты повышает эффективность обогрева устройств, то есть позволяет обеспечить им при меньших энергетических затратах более устойчивое соблюдение температурного режима в пределах температур рабочего диапазона устройства и по объему, и во времени.

Список литературы

  1. Болотовский В.И., Вайсгант З.И. Эксплуатация, обслуживание и ремонт свинцовых аккумуляторов. Л.: Энергоатомиздат, 1988. 208 с.

  2. Сайт Компании “Выбор”: https://wybor-battery.com/stati/ehkspluataciya-svincovo-kislotnyh-akkumulyatornyh-batarej-pri-otricatelnyh-temperaturah/19.12.20 г.

  3. Александров В.И., Большаков Ю.Н., Макарихин И.В., Салихов И.И. Особенности эксплуатации аккумуляторных батарей в северных условиях // Актуальные исследования. 2021. № 39 (66). С. 11.

  4. Сайт Компании “ПромПрибор” г. Москва. URL:https//www.промприбор.москва (дата обращения 07.03.2022).

  5. Сайт ООО СП “Унибелус” URL:https://www.unibelus.by/itm/K88341 (дата обращения 07.03.2022).

  6. Бухмиров В.В. Тепломассообмен. Учеб. пособие. Иваново: Ивановский государственный энергетический университет им. В.И. Ленина, 2014. 360 с.

  7. Лобасова М.С., Финников К.А., Миловидова Т.А., Дектерев А.А., Серебренников Д.С., Минаков А.В., Кузоватов И.А., Васильев В.В. Тепломассообмен [Электронный ресурс]: Курс лекций. Красноярск: ИПК СФУ, 2009.

  8. Лыков А.В. Теория теплопроводности. М.: Высшая школа, 1967. 600 с.

  9. Жуковский В.С. Основы теории теплопередачи. Л.: Энергия, 1969. 224 с.

  10. Дмитриев А.С. Введение в нанотеплофизику. М.: БИНОМ. Лаборатория знаний, 2015. 790 с.

  11. Шульженко А.А., Модестов М.Б. Программа моделирования пространственных тепловых процессов в аккумуляторе на холоде. Свидетельство о государственной регистрации программ для ЭВМ № 2021612869, Москва, 2021.

Дополнительные материалы отсутствуют.