Parameters matching of the thermoelectric system parameters for cooling heat-loaded electronics elements

Cover Page

Cite item

Full Text

Abstract

A thermoelectric cooling and thermal control system for electronic devices is considered. Based on a mathematical model using the operating characteristics of a serial thermoelectric module as initial data, calculations of the energy characteristics of a thermoelectric cooling system were carried out, taking into account its thermal resistances. The calculation results are presented in the form of diagrams, which allow for a coordinated selection of the system’s thermal resistances, ensuring the specified values of the cooling capacity and temperature difference.

Full Text

1. Введение

Одним из важнейших факторов, необходимых для обеспечения ресурса и стабильности функционирования радиоэлектронной аппаратуры, является поддержание оптимального температурного режима всех её элементов. Для выполнения этой задачи используются различные методы и технические средства. Отвод теплоты от теплонагруженных элементов может осуществляться механизмами излучения и конвекции [1, 2], для более эффективного охлаждения за счет теплопроводности, используются современные материалы, имеющие коэффициент теплопроводности в несколько раз превышающий соответствующее значение для меди [3, 4]. В то же время эти методы зачастую не имеют достаточной эффективности охлаждения и требуемой точности поддержания заданного температурного режима. Термоэлектрические системы охлаждения и терморегулирования являются перспективными для обеспечения оптимальных тепловых режимов электроники [5–7]. Они дают возможность плавного и точного регулирования температуры, а также прецизионного термостатирования [8]. Термоэлектрические устройства имеют ряд преимуществ по сравнению с другими системами охлаждения и терморегулирования: высокую надежность, отсутствие движущихся частей, практически неограниченный ресурс работы, небольшие габариты и вес, малую инерционность, независимость от ориентации в пространстве.

Основной задачей термоэлектрической системы охлаждения и терморегулирования является поддержание заданной температуры электронного устройства, имеющего мощность тепловыделения Q. Выбор оптимальной конструкции системы охлаждения является достаточно сложной и неоднозначной задачей, поскольку зависит от целого набора исходных параметров, характеризующих отдельные элементы конструкции. Причем взаимное тепловое влияние этих элементов имеет нелинейный характер, поэтому наибольшая эффективность работы системы охлаждения достигается при согласованном наборе параметров, обеспечивающем минимизацию тепловых потерь. В работе [9] был проведен анализ влияния термических сопротивлений на характеристики процесса охлаждения, направленный на обоснование выбора оптимального термоэлектрического модуля (ТЭМ) для системы охлаждения. При этом было показано, что применение высокопроизводительного ТЭМ не дает автоматического повышения эффективности системы охлаждения, оно может быть обеспечено только при выполнении достаточно жестких требований к величине термических сопротивлений. Помимо выбора оптимального ТЭМ важной задачей при разработке термоэлектрической системы охлаждения и терморегулирования является согласованный выбор ее основных параметров. В настоящей работе представлена методика согласования, которая обосновывает выбор такого соотношения термических сопротивлений и рабочих характеристик ТЭМ, при которых обеспечиваются заданные значения холодопроизводительности и температурного перепада.

2. Конструкция термоэлектрической системы охлаждения и терморегулирования

Схема термоэлектрической системы охлаждения и терморегулирования приведена на рис. 1. Теплота от теплонагруженного элемента 1 с мощностью тепловыделения Q передается на холодную сторону ТЭМ через два тепловых контакта и теплораспределитель 2, которые имеют суммарное термическое сопротивление RТ. Теплораспределитель, обычно представляющий собой прямоугольную пластину из материала с высоким коэффициентом теплопроводности, необходим для выравнивания распределения тепловой мощности, поступающей от теплонагруженного элемента на поверхность термоэлектрического модуля при несоответствии их поперечных размеров. Термоэлектрический модуль 3 отводит тепловую мощность Q и обеспечивает при собственном энергопотреблении W заданный перепад температуры ΔT0 между основанием теплонагруженного элемента и внешней средой. Устройство отвода теплоты 4, имеющее термическое сопротивление RS, передает теплоту Q + W с горячей стороны ТЭМ во внешнюю среду.

 

Рис. 1. Схема термоэлектрической системы охлаждения и терморегулирования: 1 — теплонагруженный элемент, 2 — теплораспределитель, 3 — термоэлектрический модуль, 4 — устройство отвода теплоты.

 

Термическое сопротивление теплораспределителя зависит от материала и размеров теплонагруженного элемента, для медного теплораспределителя 40 × 40 мм2 его минимальное термическое сопротивление составляет примерно 0.03 K/W при размере элемента 22.5 × 22.5 мм2 [10]. Значение термического сопротивления тепловых контактов зависит от теплопроводности наполнителя, толщины и площади зазора. При использовании стандартной термопасты КПТ-8 с коэффициентом теплопроводности λ = 0.85 Вт/мК величина термического сопротивления R = = δ/(λS) при толщине слоя δ = 0.1–0.2 мм составит 0.074–0.15 К/Вт для площади контакта S = 40 × 40 мм2 и 0.23–0.46 К/Вт для площади контакта S = 22.5 × 22.5 мм2. Значение R может быть существенно снижено при использовании паяных контактов, так низкотемпературный припой ПОСВ-50 имеет λ = 16 Вт/мК. Анализ возможных значений термического сопротивления устройства отвода теплоты RS приведена в [11, 12]. Характерный диапазон термического сопротивления устройств теплоотвода на основе воздушного теплообменника (кулера) составляет 0.1–0.5 К/Вт.

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

ΔT0 = RTQ - ΔTТЕМ(I,Q)+Rs (Q + U(I ,Q)I ),

где U, I – напряжение и сила тока электропитания ТЭМ. Первое и третье слагаемые в правой части формулы соответствуют перепаду температуры в теплораспределителе и устройстве теплоотвода, второе слагаемое – разность температур между горячей и холодной сторонами ТЭМ. Зависимости U(I, Q) и ΔTТЭМ(I, Q) определяются на основе рабочих характеристик ТЭМ с помощью интерполяционных полиномов, методика расчета этих зависимостей представлена в [11, 13]. Для серийно выпускаемых ТЭМ их рабочие характеристики обычно приводятся в технической документации производителя. Вычислительная модель, основанная на численном решении нелинейной алгебраической системы уравнений, позволяет рассчитать энергетические характеристики термоэлектрической системы охлаждения для заданных значений исходных параметров RS, RT, ΔT0 и силы тока I. Расчеты проведены на примере серийного ТЭМ «S-199–14–11» с максимальными значениями холодопроизводительности Q = 124.2 Вт и перепада температуры ΔT = 72.5 К, достигаемыми при силе тока Imax = 7.9 А и напряжении Umax = 25.3 В [14].

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

Величина холодопроизводительности численно равна тепловой мощности, которая будет отведена от теплонагруженного элемента при заданных параметрах системы охлаждения. Зависимость холодопроизводительности Q(I) для ΔT0 = –10 оС приведена на рис. 2 для следующих значений термических сопротивлений: 1RТ = 0.1 К/Вт, RS = 0.1 К/Вт, 2RТ = 0.3 К/Вт, RS = 0.1 К/Вт, 3RТ = 0.1 К/Вт, RS = 0.3 К/Вт, 4RТ = 0.3 К/Вт, RS = 0.3 К/Вт. Представленные зависимости Q(I) имеют максимумы, величина которых зависит от термических сопротивлений. Из сравнения графиков видно более существенное влияние изменения величины RS, это обусловлено тем, что через устройство теплоотвода во внешнюю среду помимо теплоты от теплонагруженного элемента дополнительно передается собственное тепловыделение ТЭМ. Значения силы тока, при которых достигаются максимумы Q, также в большей степени зависят от термического сопротивления RS. Значениям Q, которые ниже максимальных, соответствуют по два значения силы тока, более эффективный режим работы обеспечивается при меньшем значении силы тока, поскольку в этом случае собственное энергопотребление ТЭМ ниже.

 

Рис. 2. Зависимости Q(I) для ΔT0 = –10 оС: 1RТ = = 0.1 К/Вт, RS = 0.1 К/Вт, 2RТ = 0.3 К/Вт, RS = 0.1 К/Вт, 3 – RТ = 0.1 К/Вт, RS = 0.3 К/Вт, 4 – RТ = 0.3 К/Вт, RS = 0.3 К/Вт.

 

Значительное влияние на холодопроизводительность термоэлектрической системы оказывает степень охлаждения теплонагруженного элемента. Для сравнения на рис. 3 приведена зависимость Q(I) для ΔT0 = –20 оС и тех же значений термических сопротивлений: 1RТ = 0.1 К/Вт, RS = 0.1 К/Вт, 2RТ = 0.3 К/Вт, RS = 0.1 К/Вт, 3RТ = 0.1 К/Вт, RS = 0.3 К/Вт, 4RТ = 0.3 К/Вт, RS = 0.3 К/Вт. В данном случае характер зависимостей не изменился, но значения максимумов кривых от 27% (RS = 0.1 К/Вт) до 47% (RS = 0.3 К/Вт) ниже, чем на рис. 2. Следует отметить, что производители термоэлектрических систем в сопроводительной документации обычно предлагают в качестве оптимального параметра электропитания ТЭМ уровень напряжения или тока, составляющий 75% от максимальной величины [15, 16]. Рассчитанные зависимости Q(I) дают возможность определить с учетом термических сопротивлений оптимальные значения силы тока, которые могут существенно отличаться от рекомендации производителя. Так при RS = 0.3 К/Вт зависимости Q(I), приведенные на рис. 2 и 3 (синие линии), имеют максимумы в диапазоне силы тока 4.1–4.3 А, что существенно ниже рекомендованного значения I = 0.75 Imax ≈ 6 А.

 

Рис. 3. Зависимости Q(I) для ΔT0 = –20 оС: 1RТ = = 0.1 К/Вт, RS = 0.1 К/Вт, 2RТ = 0.3 К/Вт, RS = 0.1 К/Вт, 3 – RТ = 0.1 К/Вт, RS = 0.3 К/Вт, 4 – RТ = 0.3 К/Вт, RS = 0.3 К/Вт.

 

Энергетическая эффективность системы охлаждения определяется холодильным коэффициентом ε, который равен отношению холодопроизводительности Q к величине потребляемой электрической мощности W. Влияние термических сопротивлений на холодильный коэффициент исследовано в [17], где было установлено, что зависимости ε(I) имеют максимумы, величина и расположение которых существенно зависят от значений RS и RТ. Максимальные значения характеристик Q(I) и ε(I) характеризуют потенциал применяемого ТЭМ при заданных ΔT0, RТ и RS, поскольку позволяют определить предельно достижимый предел энергетической эффективности термоэлектрической системы охлаждения. Поэтому дальнейший анализ энергетических характеристик проведен применительно к максимальным значениям Qmax и εmax.

Результаты расчетов холодопроизводительности для значений ΔT = –10 оС и –20 оС приведены на плоскости RSRТ в виде изолиний, на которых указаны соответствующие им значения Qmax в Ваттах (рис. 4). Данная диаграмма дает общую картину зависимости максимальной холодопроизводительности от термических сопротивлений термоэлектрической системы охлаждения и позволяет оценивать величины RS и RТ, при которых обеспечивается заданный температурный режим теплонагруженного элемента с известной мощностью тепловыделения. Так, например, для ΔT0 = –20 оС значение Q = 60 Вт может быть достигнуто при RТ = 0.1 К/Вт и RS ≈ 0.05 К/Вт, эти значения отмечены на рисунке штриховыми линиями, а место пересечения на изолинии выделено кружком. При этом каждой точке на изолинии соответствует своя согласованная пара значений RS и RТ. В то же время холодопроизводительность 60 Вт не может быть получена при RS > 0.08 К/Вт для ΔT0 = –20 оС и RS > 0.13 К/Вт для ΔT0 = –10 оС. Термическое сопротивление теплораспределителя в меньшей степени ограничивает величину холодопроизводительности, величина Q = 60 Вт не может быть получена при RT > 0.29 К/Вт для ΔT0 = –20 оС и RT > 0.46 К/Вт для ΔT0 = –10 оС. Более сильное ограничивающее влияние величины RS на холодопроизводительность объясняется увеличенным перепадом температуры в устройстве теплоотвода при отводе суммарной тепловой мощности Q + W. Из проведенного анализа следует, что представленная диаграмма позволяет определить диапазоны значений RS и RT, при которых достигается заданная холодопроизводительность термоэлектрической системы охлаждения. Определение силы тока, обеспечивающей заданную холодопроизводительность, производится с помощью зависимостей Q(I) (рис. 2).

 

Рис. 4. Диаграмма холодопроизводительности, на кривых указаны значения Qmax в Ваттах.

 

Как было показано выше, заданная холодопроизводительность может быть достигнута в определенных диапазонах термических сопротивлений теплораспределителя и устройства теплоотвода, но при этом возможны различные сочетания пар значений RS и RТ. Для выбора значений термических сопротивлений, соответствующих наибольшей эффективности системы охлаждения, на рис. 5 приведена диаграмма, на изолиниях которой указаны значения εmax. Значение холодильного коэффициента определяется по выбранной с помощью рис. 4 паре значений RS и RТ для заданного значения Q, так для выделенной на рис. 4 точки соответствующее ей значение εmax = 1.24. Изолинии εmax на большей части плоскости RSRТ имеют меньший угол наклона к оси RS по сравнению с изолиниями Qmax, поэтому наиболее высокое значение холодильного коэффициента будет получено при RТ = 0, так в рассматриваемом случае изолинии Qmax = 60 Вт соответствует значение εmax = 1.33. Следует отметить, что на практике нулевые значения термического сопротивления не могут быть получены, поэтому при анализе следует исходить из значений RS и RТ для реальных конструкций теплораспределителя и устройства теплоотвода с массогабаритными параметрами, соответствующими техническим требованиям к термоэлектрической системе охлаждения. Дополнительно диаграммы холодопроизводительности и холодильного коэффициента позволяют рассчитать собственное энергопотребление ТЭМ, определяемое из соотношения W = Q/ε.

 

Рис. 5. Диаграмма холодильного коэффициента, на кривых указаны значения εmax.

 

Заключение

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

Конфликт интересов

Автор заявляет об отсутствии конфликта интересов.

×

About the authors

Е. N. Vasil’ev

Institute of Computational Modelling, Siberian Branch, Russian Academy of Sciences

Author for correspondence.
Email: ven@icm.krasn.ru
Russian Federation, Krasnoyarsk

References

  1. Martyushev S.G., Sheremet M.A. Two factors affecting the cooling rate of fuel elements in sealed units // Mikroelektronika. 2014. V. 43. No 5. P. 390–398. doi: 10.7868/S0544126914050056
  2. Boutina L., Bessaih R. Numerical simulation of mixed convection air-cooling of electronic components mounted in an inclined channel // Applied Thermal Engineering. 2011. V. 31. P. 2052–2062. https://doi.org/ 10.1016/j.applthermaleng.2011.03.021
  3. Glinskii I.A., Zenchenko N.V. Computer simulation of the heat distribution element for high-power microwave transistors // Russian Microelectronics. 2015. V. 44, No 4. P. 236–240. doi: 10.1134/S1063739715040058
  4. Zuev S.M., Prokhorov D.A., Maleev R.A., Debelov V.V., Lavrikov A.A. Application of Graphene in The Cooling System of a Personal Electronic Computer // Russian Microelectronics. 2021. V. 50. No 6. P. 404–411. https://doi.org/10.1134/S1063739721050097
  5. Chang Y.W., Chang C.C., Ke M.T., Chen S.L. Thermoelectric air-cooling module for Chang electronic devices // Applied Thermal Engineering. 2009. V. 29. № 13. P. 2731–2737. https://doi.org/10.1016/j.applthermaleng.2009.01.004
  6. Shtern M.Yu., Shtern Yu.I., Sherchenkov A.A. Thermoelectric systems for providing thermal modes of computer technology // Izvestija Vysshykh Uchebnykh Zavedenii. Elektronika. 2011. No 4. P. 30–38.
  7. Vasil’ev E.N. Thermoelectric cooling of heat-loaded electronics // Russian Microelectronics. 2020. V. 49. No 2. P. 123–131. https://doi.org/10.1134/S1063739720020092
  8. Zagorodnov A.P., Yakunin A.N. Precision thermostatting of a resonator on volume acoustic waves. Modeling and synthesis of a control system // Zhurnal Radioelektroniki. 2013. No 10. P. 1–14, http://jre.cplire.ru/iso/oct13/12/text.pdf
  9. Vasil’ev E.N. On the importance of thermal resistances of the cooling system when choosing a thermoelectric module // Technical Physics. 2023. V. 68. No 5. P. 574–579, https://doi.org/10.21883/tp.2023.05.56062.13-23
  10. Vasil’ev E.N. Calculation of the thermal resistance of a heat distributer in the cooling system of a heat-loaded element // Technical Physics. 2018. V. 63. No 4. P. 471–475. https://doi.org/10.21883/JTF.2018.04.45714.2312
  11. Vasil’ev E.N. Determination of thermoelectric cooling modes of heat-loaded electronics // Russian Microelectronics. 2020. V. 49. No 4. P. 278–284, https://doi.org/10.31857/S0544126920030072
  12. Vasil’ev E.N. Optimization of thermoelectric cooling regimes for heat-loaded elements taking into account the thermal resistance of the heat-spreading system // Technical Physics. 2017. V. 62. No 9. P. 1300–1306, https://doi.org/10.1134/S1063784217090286
  13. Vasil’ev E.N. Calculation and optimization of thermoelectric cooling modes of thermally loaded elements // Technical Physics. 2017. V. 62. No 1. P. 90–96. https://doi.org/10.21883/JTF.2017.01.44022.1725
  14. https://crystalltherm.com/upload/iblock/5af/1knj372x6ho3v82cwzufypmaktsqm4ph/TM_S_199_14_11_L2_SPEC.pdf. Accessed September 13, 2024.
  15. http://kryotherm.ru/ru/assembly-instructions.html. Accessed September 13, 2024.
  16. http://ecogenthermoelectric.com/ru/texnicheskaya-podderzhka.html. Accessed September 13, 2024.
  17. Vasil’ev E.N. The Effect of thermal resistances on the coefficient of performance of a thermoelectric cooling system // Technical Physics. 2021. V. 66. No 6. P. 815–819. https://doi.org/10.1134/S1063784221050248

Supplementary files

Supplementary Files
Action
1. JATS XML
2. Fig. 1. Schematic diagram of a thermoelectric cooling and thermal control system: 1 — heat-loaded element, 2 — heat distributor, 3 — thermoelectric module, 4 — heat removal device.

Download (102KB)
3. Fig. 2. Q(I) dependencies for ΔT0 = –10 oC: 1 – RT = 0.1 K/W, RS = 0.1 K/W, 2 – RT = 0.3 K/W, RS = 0.1 K/W, 3 – RT = 0.1 K/W, RS = 0.3 K/W, 4 – RT = 0.3 K/W, RS = 0.3 K/W.

Download (87KB)
4. Fig. 3. Q(I) dependencies for ΔT0 = –20 оС: 1 – RT = 0.1 K/W, RS = 0.1 K/W, 2 – RT = 0.3 K/W, RS = 0.1 K/W, 3 – RT = 0.1 K/W, RS = 0.3 K/W, 4 – RT = 0.3 K/W, RS = 0.3 K/W.

Download (84KB)
5. Fig. 4. Cooling capacity diagram, the curves show Qmax values ​​in Watts.

Download (118KB)
6. Fig. 5. Coefficient of performance diagram, with εmax values ​​indicated on the curves.

Download (130KB)

Copyright (c) 2024 Russian Academy of Sciences

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

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

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».