Математическая модель жидкостного погружного охлаждения вычислительных устройств

ISSN 2079-3316 ПРОГРАММНЫЕ СИСТЕМЫ: ТЕОРИЯ И ПРИЛОЖЕНИЯ №1(28), 2016, с. 187-199 УДК 004.3'124

А. А. Ахременков, А. М. Цирлин

Математическая модель жидкостного погружного охлаждения вычислительных устройств

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

Ключевые слова и фразы: жидкостное охлаждение, кипящий хладагент, естественная циркуляция.

Введение

Повышение мощности и компактности современных вычислительных устройств значительно повысило требования к системе охлаждения. Применение диэлектрических хладагентов позволяет погрузить вычислительный узел в жидкость [1—3]. Эффективность погружного жидкостного охлаждения значительно превышает эффективность воздушного, поскольку жидкость обладает большой теплоемкостью и коэффициентом теплоотдачи. Применение низкокипящих диэлектрических жидкостей позволяет еще более повысить эффективность систем охлаждения, т.к. при фазовом переходе теплоемкость и коэффициент теплоотдачи резко возрастают. В данной работе рассмотрены две математические модели погружных систем охлаждения: с использованием фазового перехода и без него.

1. Жидкостная система охлаждения погружного типа

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

© © ©

А. А. Ахременков, А. М. Цирлин, 2016

Институт программных систем имени А. К. Айламазяна РАН, 2016 Программные системы: теория и приложения, 2016

Рис. 1. Структура системы погружного жидкостного охлаждения конвективного типа

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

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

1.1. Математическая модель

Будем предполагать, что емкость теплоизолирована. Подвод теплоты происходит только от охлаждаемого вычислительного устройства, высота тепловыделяющего устройства равна Ь# < Ь, где Ь [м] — общая высота установки, (х) [°С], ад(х,т) [ВтС-1 • м-1]—функции температуры поверхности и коэффициента теплоотдачи для вычислительного устройства, зависящие от текущей координаты по высоте х [м] и от скорости потока т [м-е-1].

Введем эффективный коэффициент теплоотдачи г Lh

(1) / ан(x,w)dx = ан(w), ан(x,w) = 0 Уж ф [0,Lh].

Jo

Для контакта жидкость-металлическая поверхность зависимость ан (w) [Вт C-1] имеет вид [4]

(2) a(w) = S (350 + 2000Vw),

где S [м2] — площадь контакта.

Плотность теплового потока q+(x) [Вт • м-1], отводимого в сечении х, равна

(3) Ч+(х) = ан(x,w)(tH(х) - t+ (х)),

где t+(x) [°C] — температура охлаждающей жидкости в сечении х.

Суммарная тепловая нагрузка q [Вт] (мощность охлаждаемого вычислительного узла) задана и равна

/¿я

(4) / q+(x)dx = q.

Jo

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

. f to + ^ f0X q+m^ х Ф [0, LH]

(5) f+(x)=\ to + JL, x Ф [LH,L] ,

где to [°C] — температура охлаждающей жидкости, подаваемой в секцию с вычислительным узлом, д [кг • с-1] — расход жидкости, С [Дж • кг-1 C-1] — удельная весовая теплоемкость жидкости.

Продифференцируем по х выражение (5) на отрезке [0, Lh] с учетом (3), предполагая, что ан(x,w) = ан(w) = S(350 + 2000^/w)/Lh, tn(x) = tn = const. Получим

(6) + t+(x) - tH = 0. ан (x) ax

Решение уравнения (6) с граничным условием i+(0) = to имеет вид

С* JJ ( w) X

(7) t+(x)= tH + (t0 - tH)e ^ .

Подставив полученное решение в равенство (4), найдем значение to:

(8) to = tn--7-Яан (w)Lh-

w 1 — e

-у

В секции охлаждения находится теплообменник, охлаждающий жидкость. Он характеризуется функциями температуры tc (х) [°C] и коэффициента теплоотдачи ac(x,w) [Вт C-1 • м-1]. Плотность теплового потока q-(x) [Вт • м-1], отбираемого в сечении х, равна

(9) q-(x) = ac(x,w)(t-(х) — tc(х)),

где t-(x) [°C] — температура охлаждающей жидкости в сечении х. Количество теплоты, отводимой в секции ,L

(10) / q-(x)dx = q.

JL-Lc

Температура охлаждающей жидкости в сечении х, в секции охлаждения равна

(11) i-(x) = { С+ Стд ¡0 9-(0dl, х J[L — -—-] .

Предполагая, что

ac(х, w) = ac(w) = 5(350 + 2000\f(w))/—c, tc(x) = tc = const,

продифференцируем по x выражение

Cq dt- . .

(12) -f--^ — t-(x)+ tc =0

ac ( x) x

с граничным условием -( — — —c) = 0 имеет вид

ac

(13) t-(x)= tc + (t о — tc) e c— .

Подставив полученное решение в равенство (10), найдем значение tc:

(14) fc = 10--/ ac (!,)Lc-

we °я — 1

Разность температур холодной и горячей частей жидкости приводит к различию плотностей холодной р-(х) [кг • м-3] и горячей р+ (х) [кг • м-3] частей жидкости, что порождает движущую силу. Перепад давления Ар [Па] соответственно равен

(15) Ар = даЬ [ (р-(х) - р+(х))йх,

Jo

где да [м • с-"2] — ускорение свободного падения.

Плотность жидкости линейно зависит от температуры

(16) р(^=ро + Р1(г + 273),

здесь р\ [К-1] — коэффициент термического расширения.

Подставим выражения (16), (13), (7) в формулу (15) и проинтегрируем полученное выражение:

[Ь , ч , ,т т ч -т &0 - ^с)дС ( »с_ \

г-(х)Лх = го(ь -Ьс)+ ^Ьс + (——Ст~ е ас - Ч ,

7о ас (х) V /

I t + (х)1х = НЬН + (Ь - Ьн)(^о + -

(tо -1н)дС / - °д

ан (т

М V

е -1

)

откуда

Ар = д аЬр1

/ Ь-(х)(1х - Ь +(х)(1х

о - о

Создаваемое давление совершает работу, связанную с созданием скоростного потока Арск [Па] и преодолением гидравлического сопротивления А тр [Па]:

(18) АР = АР ск + АР тр.

Согласно [5] потери давления, связанные с созданием скоростного потока, равны

(19) Арск = т2--/2,

где т — скорость потока в секции охлаждения, р- [кг • м-3] — средняя плотность в секции охлаждения. Потеря давления, связанная с трением, равна

Ь т2--

Артр = X

1е 2

где X — коэффициент трения

х = ,

Пе

Ъ

где Ке - критерий Рейнольдса:

(22) Ке = т с1,е р~/р.

Здесь р [Па-с] — динамический коэффициент вязкости. Подставим выражения (21) и (22) в (20), получим

57 Ьт

(23) Дтр =

Расход жидкости, протекающей со скоростью т в трубе диаметра ¿е, равен

(24) д = 0.25т— ¿2е р-.

Выразим из этого выражения т(д) и подставим в выражения (18)-(20). Найдем Ар(д):

(25) А ( ) ^д2 + 114 Ь д

(25) Др(9) = —щ- + р

Подставим полученное выражение в уравнение (17) и приведем его к итерационной форме

(26) ^ ^ Г ' 8д2 114Ьд

9г+1 =1 ^ 9 °Ьр\ J 1-{х)(кс — J 1+(х)<кс

-р +

+ 91,

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

Ч

90 = 5С'

что соответствует расходу, обеспечивающему перепад температуры в 5°0.

Решение полученного уравнения итерационным методом позволяет найти значение расхода д, после чего по формулам (8), (14) могут быть найдены температуры Ь о, Ьс.

1.2. Пример

Решим задачу (26) для следующих начальных данных. Характеристики охлаждающей жидкости:

р(г) = 929.325 — 0.89625£[кг • м-3], С = 1605[Дж/кг К],

р = 0.004555[Па • с].

Рис. 2. Зависимость между отбираемой мощностью д [Вт] и высотой столба хладагента Ь [м] при фиксированных значениях Ьс = 30°С, Ьн = 57°С

Характеристики теплоотводящей секции:

5 = 0.05[м2], ас = 120[Вт С-1], ¿е = 0.05[м].

Характеристики тепловыделяющей секции:

ан = 120[Вт С-1], гя = 57[°С].

Коэффициент сжатия для решения итерационного уравнения (26) брался равным 7 = 0.01. На рис. 2 представлена зависимость отбираемой мощности от высоты емкости при фиксированном значении температуры 1с = 30° С.

2. Жидкостная система охлаждения погружного типа с использованием кипящего хладагента

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

ООО ООО ООО ООО ООО

:

Рис. 3. Структура системы погружного жидкостного охлаждения с кипящем хладагентом

и конденсатора, в котором получившийся пар конденсируется и стекает обратно в емкость.

2.1. Математическая модель

Будем предполагать, что в стационарном режиме температура жидкости равна температуре кипения Ь в [°С] и вся теплота, выделяемая вычислительным устройством, уходит на испарение. Расход жидкости при испарении существенно меньше, чем при жидкостном охлаждении (с одинаковой тепловой нагрузкой Ч), следовательно, высота столба хладагента Ь ни как ни влияет на работу системы охлаждения. Поскольку температура вычислительного устройства Ь н, температура охлаждающей жидкости Ь в, удельный коэффициент теплоотдачи ав [Вт С-1 • м-2] и площадь контакта Яв [м2] заданы, то предельная отводимая мощность Ч [Вт] равна

Ч = Явав (Ък — Ьв).

В этом случае массовый поток пара дп [кг • с 1] равен

п = Ч/ ,

где г [Дж • кг 1] — удельная теплота парообразования охлаждающей жидкости. Объемный расход пара ду [м3 • с-1] равен

ду = дпРп,

где рп [кг • м 3] — плотность пара.

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

Удельный коэффициент теплоотдачи ав при пузырьковом кипении равен [5]

2/3"

0,075 + 0, 75 ' ' '

(30) ав

Ад2

г/а(273 + Ьв)'

где рж [кг • м-3] — плотность жидкости, А [Вт • м-1 • К-1] — коэффициент теплопроводности, V = р/рж [м2 • с-1] — кинематическая вязкость, <г[Н • м-1] — коэффициент поверхностного натяжения.

В конденсатор пар поступает при температуре кипения и, конденсируясь, отдает теплоту охлаждающей жидкости. Будем предполагать, что вся жидкость после конденсации попадает в емкость при температуре кипения. Количество теплоты дк [Вт], отводимое при конденсации, равно

(31) дк = кк(гв - гс),

где Ьс [°С] — температура охлаждающей жидкости, Кк [Вт С-1] — коэффициент теплопередачи,

(32) Кк = ^ + ^ + ,

ак + ас + аж

и ак [Вт • м-2 С-1 ] — удельный коэффициент теплоотдачи для конденсирующего пара, ас [Вт • м-1 С-1] — теплопроводность стенки теплообменника, 6 [м] — толщина стенки, аж[Вт • м-2 С-1] — удельный коэффициент теплоотдачи нагреваемой жидкости, Як [м2] — площадь контакта.

Удельный коэффициент теплоотдачи для конденсирующего пара ак в конденсаторе из пучка вертикальных труб высотой Н [м] равен [5]

(33) ак = 1,15 4 АЗрж9'

м(^в - tc)Н

где р [Па-с] — динамическая вязкость, да — ускорение свободного падения. Площадь контакта в конденсаторе, содержащем пучок из п труб диаметра с! [м] и высоты Н [м], равна

(34) Як = пж(Ш.

Из условия теплового баланса Ч = Чк можно определить температуру охлаждающей жидкости

(35) ±с = г в — -¡Ч-.

Кк

3. Сравнение двух типов охлаждения

Для сравнения двух типов охлаждения рассмотрим систему жидкостного охлаждения со следующими параметрами: Яв=0.05 м2, Ьс = 25°С, Ь о = 35.8°С, соответствующее значение тепловой нагрузки равно ч = 432 Вт, высота ванны Ь = 0.6 м, максимальная температура охлаждающей жидкости на выходе из холодильника равна 46.9°С.

Для следующих исходных данных рассчитаем характеристики системы охлаждения:

Яв = 0.05[м2], Л = 0.124[Вт • м-1 • К-1], рж = 911[кг • м-3],

V = — =2.8 * 10-6[м2 • с-1], а = 20 * 10-3[Н • м-1], Рж

рп = 1.2[кг • м-3], гн = 60°С, гв = 50°С, г = 90[кДж • кг-1].

Найдем удельный коэффициент теплоотдачи а в по формуле (30):

ав = 1470 [Вт С-1 • м-2],

соответствующее значение тепловой нагрузки (27) равно 735 [Вт], коэффициент теплопередачи при конденсации (33) Кк = 93.6 [ВтС-1], максимально отбираемая теплота = 2340 [Вт].

Для сравнения двух систем зафиксируем температуры охлаждающей жидкости Ьс = 25°С, вычислительного устройства £ н = 60°С, тепловыделяющей площади в вычислительном устройстве и холодильнике/конденсаторе. Коэффициент теплоотдачи при кипении (30) в 4.2 раза больше, чем при нагреве (2), разность температур при кипении равна 10°С, а при нагреве — 25° С. Таким образом, теплосъем в системе с кипящем хладагентом в 1.7 раз больше. При конденсации паров коэффициент теплоотдачи (33) в 5.2 раза больше, чем коэффициент теплообмена в холодильнике, разность температур при конденсации

равна 25°С, а при охлаждении — 15°С. Следовательно, эффективность охлаждения при конденсации в 25/15 * 5.2 = 8.7 раза больше, чем при охлаждении в системе с некипящим хладагентом.

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

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

[1] В. С. Лытарь. «Реализация суперкомпьютеров с массивно-параллельной архитектурой и прямым жидкостным охлаждением на примере стенда RSC Petastream», Наука и образование в XXI веке, Сборник научных трудов по материалам международной научно-практической конференции. Т. 15 (31 октября 2014 г.), ООО «Консалтинговая компания Юком», Тамбов, 2014, с. 88-90. t 187

[2] Ю. Хомуткий. «Жидкостное охлаждение в ЦОД: опыт, выводы, перспективы», Журнал сетевых решений LAN, 2013, №4, с. 42-47. t 187

[3] В. П. Киц, А. А. Новиков, С. А. Пахомова. «Эффективность жидкостного охлаждения центрального процессора персонального компьютера», Инженерный вестник, 2015, №1, с. 1-2. t 187

[4] Х. Кухлинг. Справочник по физике, Мир, М., 1985. t 189

[5] А. Г. Касаткин. Основные процессы и аппараты химической технологии, Изд. 9-е, Химия, М., 1973, 750 с. t 191,196

Рекомендовал к публикации Программный комитет

Четвёртого национального суперкомпьютерного форума НСКФ-2015

Об авторах:

Андрей Александрович Ахременков ИПС им. А.К. Айламазяна РАН, с.н.с. ИЦСА. Области исследований: термодинамика при конечном времени, промышленное моделирование

e-mail: andrei@svp.botik.ru

Анатолий Михайлович Цирлин ИПС им. А.К. Айламазяна РАН, г.н.с. ИЦСА. Области исследований: теория оптимального управления, условия оптимальности усредненных задач, термодинамика при конечном времени, термодинамический подход к анализу экономических систем

e-mail: tsirlin@sarc.botik.ru

Пример ссылки на эту публикацию:

А. А. Ахременков, А. М. Цирлин. «Математическая модель жидкостного погружного охлаждения вычислительных устройств», Программные системы: теория и приложения, 2016, 7:1(28), с. 187-199. URL: http://psta.psiras.ru/read/psta2016_1_187-199.pdf

Andrei Akhremenkov, Anatoly Tsirlin. Mathematical model of liquid immersion cooling system for supercomputer.

Abstract. In this paper, a model of the cooling system for supercomputer with their direct contact with the refrigerant. The first part considers the capabilities of the system with natural circulation of the liquid refrigerant. The second for the case of boiling refrigerant in the heater and condensed in the condenser. In all cases there are no pumps to pump the refrigerant. In the third part a comparison of these two systems. (In Russian).

Key words and phrases: immersion cooling system, boiling refrigerant, natural circulation.

References

[1] V. S. Lytar'. "Implementation of supercomputers with massively parallel architecture and direct liquid cooling considering RSC Petastream stand", Nauka i obrazovaniye v XXI veke, Sbornik nauchnykh trudov po materialam mezhdunarodnoy nauchno-prakticheskoy konferentsii. V. 15 (31 oktyabrya 2014 g.), OOO "Konsaltingovaya kompaniya Yukom", Tambov, 2014, pp. 88—90 (in Russian).

[2] Yu. Khomutkiy. "Liquid cooling in the data center: experience, conclusions, perspectives", Zhurnal setevykh resheniy LAN, 2013, no.4, pp. 42—47 (in Russian).

[3] V. P. Kits, A. A. Novikov, S. A. Pakhomova. "The efficiency of PC CPU liquid cooling", Inzhenernyy vestnik, 2015, no.1, pp. 1—2 (in Russian).

[4] Kh. Kukhling. Handbook of physics, Mir, M., 1985 (in Russian).

[5] A. G. Kasatkin. Basic processes and apparatuses of chemical technology, Izd. 9-ye, Khimiya, M., 1973 (in Russian), 750 p.

Sample citation of this publication:

Andrei Akhremenkov, Anatoly Tsirlin. "Mathematical model of liquid immersion cooling system for supercomputer", Program systems: theory and applications, 2016, 7:1(28), pp. 187-199. (In Russian).

URL: http://psta.psiras.ru/read/psta2016_1_187- 199.pdf

© A. A. Akhremenkov, A. M. Tsirlin, 2016 © Ailamazyan Program System Institute of RAS, 2016 © Program systems: Theory and Applications, 2016