Исследование процесса теплоотдачи в системе охлаждения супер электронно-вычислительной машины с тепловым насосом

Целью оптимизации систем охлаждения супер ЭВМ является интенсификация процесса теплоотдачи в окружающую среду при повышении температуры окружающей среды и снижении требований к кондиционированию и вентиляции ЦОД (центров обработки данных), для достижения энергосбережения ЦОД в целом. При этом необходимым является сохранение качества термического менеджмента во избежание повышения температур процессоров и срабатыванию системы тепловой защиты процессоров. В[1], [2]рассматриваются явления тепло- и массобмена в радиоэлектронной аппаратуре и методы проектирования систем охлаждения электронных устройств с применением естественной и вынужденной конвекции в воздушной среде. В данной статье поставлена задача оптимизации системы охлаждения процессоров супер ЭВМ с тепловым насосом на обратном цикле Стирлинга путем численного моделирования процессов теплоотдачи в исследуемом объекте. Системой охлаждения с тепловым насосом на обратном цикле Стирлинга называется система, в которой охлаждаемые электронные компоненты соединены через магистраль передачи тепла с тепловым насосом, а тот, в свою очередь, соединен через магистраль передачи тепла с радиатором, рассеивающим тепловую энергию в окружающую среду. Тепловым насосом на обратном цикле Стирлинга называется в данном случае компрессионный тепловой насос, представляющий собой двигатель Стирлинга, работающий в режиме холодильной машины, с электрическим приводом, с воздухом или гелием в качестве рабочего тела[3,4,5].

В процессе работы тепловой насос потребляет электроэнергию. Соотношение вырабатываемой тепловой энергии и потребляемой электрической называется коэффициентом трансформации, или коэффициентом преобразования теплоты COP - (coefficientofperformance) и служит показателем эффективности теплового насоса [6]. Для вычисления COP используется следующая формула:

, (1)

Где COP - коэффициент преобразования;

А - работа, совершенная насосом, Дж;

Qin-теплота, перенесенная тепловым насосом из источника тепла, на перенос которой затрачена работа A, Дж.

Для тепловых насосов коэффициент трансформации СОР>1. При этом тепловой насос создает лифт температур ДФ.

Рассмотрим тепловой насос на основе двигателя Стирлинга, работающего в режиме холодильной машины, б-типа. Такой насос содержит два рабочих цилиндра, "холодный" и "горячий", соединенных рекуператором. К "холодному" цилиндру термически присоединена упаковка процессора, а "горячий" цилиндр омывается охлаждающей жидкостью и магистралью передачи тепла соединен с радиатором, рассеивающим тепловую энергию в окружающую среду.При этом термические сопротивления в модели термического сопротивления для системы охлаждения на участках упаковка процессора - холодный цилиндр и горячий цилиндр - охлаждающая жидкость уменьшают лифт температур теплового насоса и приводят к повышению термического сопротивления системы охлаждения в целом. Для оценки повышения термического сопротивления системы охлаждения произведем численное моделирование процессов теплопередачи на участках упаковка процессора - холодный цилиндр и горячий цилиндр - охлаждающая жидкость для диапазона коэффициентов теплоотдачи и коэффициента трансформацииCOP.

Экспериментальный тепловой насос на основе двигателя Стирлинга, работающего в режиме холодильной машины, б-типа, содержит два рабочих цилиндра, "холодный" и "горячий", представлен на рис. 1. На рисунке цифрами обозначены: холодный цилиндр 1, упаковка процессора 2, рабочая камера с охлажденным газом 3, поршень холодного цилиндра 4, горячий цилиндр 5, рабочая камера с нагретым газом 6, поршень холодного цилиндра 7. Поверхность горячего цилиндра омывается охлаждающей жидкостью.

Рис. 1. -Холодный (А) и горячий (Б) цилиндры теплового насоса.

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

В работах [7-10] был предложен метод трехмерного расчета тепловых потоков в системах охлаждения процессоров. Процесс теплопередачи в случае холодного цилиндра идет по пути полупроводниковый кристалл процессора-упаковка процессора- холодный цилиндр--камера с охлажденным газом. Процесс теплопередачи в случае горячего цилиндра идет по пути камера с нагретым газом - горячий цилиндр - охлаждающая жидкость. Процесс теплопередачи в зонах контактов будет описан через граничные условия второго и третьего рода соответственно.

Стационарное уравнение теплопроводности - уравнение Лапласав декартовых координатах:

, (2)

где: - температура, - координаты;

зададим граничные условия для холодного цилиндра:

для (граничные условие второго рода):

,

,

где - тепловой поток с поверхности полупроводникового кристалла процессора, - мощность теплового потока, выделяемая с площади полупроводникового кристалла процессора;

для (граничное условие третьего рода на основании закона Ньютона-Рихмана):

,

, (3)

,

где - тепловой поток с поверхности камеры в охлажденный газ, - средний коэффициент теплоотдачи, - коэффициент теплопроводности, - температура охлажденного газа, - мощность теплового потока, выделяемая с площади камеры;

условие сохранения энергии:.

Система дифференциальных уравнений решалась методом конечных элементов.

Для численного решения зададим граничные условия, а именно:

мощность теплового потока ;

Внутренний диаметр камеры 38 мм и высота камеры 12 мм, внешний диаметр цилиндра 46 мм;

Температура газа в камере260К.

Коэффициент теплопроводности упаковки процессора постоянный и равен коэффициенту теплопроводности меди, коэффициент теплопроводности холодного цилиндра постоянный и равен коэффициенту теплопроводности алюминия.

зададим граничные условия для горячего цилиндра:

для (граничное условие третьего рода):

,

, (4)

,

где - тепловой поток на поверхность рабочей камеры из нагретого газа,

- средний коэффициент теплоотдачи в рабочей камере, - коэффициент теплопроводности цилиндра, - температура нагретого газа, - мощность теплового потока, выделяемая с учетом COP из камеры на ее площадь ;

условие сохранения энергии: .

для (граничное условие третьегорода):

,

, (5)

,

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

Внутренний диаметр камеры 38 мм и высота камеры 12 мм, внешний диаметр цилиндра 46 мм, материал - алюминий;

Температура охлаждающей жидкости300К;

.

Коэффициент теплопроводности горячего цилиндра постоянный и равен коэффициенту теплопроводности алюминия.

Полученное в результате численного моделирования распределение температурного поля по сечению по оси симметрии цилиндров теплового насоса в зависимости от величины коэффициента теплоотдачи и коэффициента трансформации представлены на рис 2.

Рис. 2. -Распределение температурного поля по сечению по оси симметрии цилиндров теплового насоса; для холодного цилиндра в зависимости от величины коэффициента теплоотдачи: а) 500, б) 1000, в) 2000; для горячего цилиндра в зависимости от коэффициента трансформации COPг) 6,5, д) 3,25, е) 2,2.

Результаты численного моделирования зависимости потерь температурного лифта теплового насоса системы охлаждения на участке упаковка процессора - холодный цилиндр от коэффициента теплоотдачи в камере холодного цилиндра приведены на рис. 3.

Рис. 3. -График зависимости потерь температурного лифта теплового насоса системы охлаждения на участке упаковка процессора - холодный цилиндр от коэффициента теплоотдачи в камере холодного цилиндра.

Результаты численного моделирования зависимости потерь температурного лифта теплового насоса системы охлаждения на участке горячий цилиндр - охлаждающая жидкость от коэффициента трансформацииCOP приведены на рис. 4.

Рис. 4. -- График зависимости потерь температурного лифта теплового насоса системы охлаждения на участке горячий цилиндр - охлаждающая жидкость от коэффициента трансформации COP.

Методом численного моделирования процесса теплоотдачи получены зависимости потерь температурного лифта теплового насоса системы охлаждения на участках упаковка процессора - холодный цилиндр и горячий цилиндр - охлаждающая жидкость для диапазона коэффициентов теплоотдачи и коэффициента трансформации COP.

Из результатов численного эксперимента можно сделать вывод, что наибольший вклад в потери температурного лифта теплового насоса системы охлаждения на участках упаковка процессора - холодный цилиндр. Таким образом, система охлаждения с использованием теплового насоса имеет максимальный вклад в общее термическое сопротивление в месте контакта с упаковкой процессора, аналогично системам воздушного и водяного охлаждения. Таким образом, термический контакт упаковка процессора - холодный цилиндр требует применения методов оптимизации термического сопротивления, аналогичных описанным в работах [8-10].

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

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

Литература