ОПТИМИЗАЦИЯ КОНСТРУКЦИИ СИСТЕМЫ ОХЛАЖДЕНИЯ ЭЛЕКТРОННОГО МОДУЛЯ
Оптимизация конструкции системы охлаждения электронного модуля
61
УДК 62-712
Г. В. БИРЮЛИН, В. И. ЕГОРОВ, С. Ю. МУРОВ, И. В. НЕВСКИЙ
ОПТИМИЗАЦИЯ КОНСТРУКЦИИ СИСТЕМЫ ОХЛАЖДЕНИЯ ЭЛЕКТРОННОГО МОДУЛЯ
Рассматривается конструкция биметаллического (медь—алюминий) паяного радиатора, эффективно охлаждающего элементы электронной аппаратуры с высокой неравномерностью тепловых потоков. Представлена методика оптимизационного расчета конструктивных параметров радиатора.
Ключевые слова: принудительное воздушное охлаждение, радиатор, конечноэлементный метод, оптимизационный расчет.
Функционирование мощных электронных изделий малых размеров требует использования охлаждающих устройств, минимизация массы которых является важной и актуальной задачей. Для принудительного воздушного охлаждения применяются биметаллические ребристые радиаторы [1], что позволяет в областях с низкой плотностью теплового потока использовать металлы с малым удельным весом.
В современных электронных устройствах применяются изделия электронной техники с контактной площадью порядка 4·10–4 м2 и мощностью тепловыделений до 300 Вт. Площадь основания ребристых радиаторов, используемых для обеспечения требуемой температуры корпуса силового элемента, часто значительно превосходит контактную площадью элемента. В этом случае при большом значении коэффициента теплоотдачи оребренной поверхности локальным источником теплоты создается неравномерность температурного поля основания радиатора. В данных условиях проводить расчеты конструкции радиатора, основываясь на известных из литературы [2] значениях тепловых сопротивлений, недопустимо.
Плотности теплового потока в зоне расположения локального источника теплоты и в областях, удаленных от источника, существенно различаются. Поэтому в области, близкой к источнику теплоты, используют высокотеплопроводные материалы. На краях основания радиатора можно применить материал (сплав алюминия), имеющий меньшую плотность и теплопроводность по сравнению с медными сплавами.
Такое техническое решение можно реализовать, используя технологию создания паяных радиаторов, разработанную в БНТ „Прибой“ (Санкт-Петербург) [1]. Такие радиаторы состоят из спаянных между собой ребер и планок. Основание радиатора образуют планки и торцы ребер. Каждое ребро запаяно в два параллельных основания. Толщина основания есть толщина планки. Как планки, так и ребра могут быть медными и алюминиевыми в различных сочетаниях, например, по краям — из сплава алюминия, а в центре — из меди.
Рассмотрим вариант оптимизации конструкции системы охлаждения электронного модуля в целях минимизации его веса. Для обеспечения требуемого температурного режима электронных компонентов используется принудительное воздушное охлаждение. Несущими элементами конструкции (см. рисунок) являются пластины, образующие каналы прямоугольного сечения. На входе и выходе каналов (длиной L = 0,25 м) закреплены вентиляторы. На двух противоположных пластинах (основаниях радиатора) размещены тепловыделяющие элементы. Ширина оснований фиксирована. На основании 1 расположены два транзистора мощностью по 175 Вт, отстоящие на расстояние 10 мм. Удельная мощность тепловыделений каждого из них составляет 8·105 Вт/м2. Кроме того, на основании 1 размещены также три транзистора мощностью по 15 Вт. Размеры контактных площадок транзисторов 24×9 мм. На основании 2 расположены три преобразователя мощностью по 45 Вт, общее тепловыделение которых составляет 140 Вт. Размеры контактных площадок преобразователей 60×70 мм. При
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
62 Г. В. Бирюлин, В. И. Егоров, С. Ю. Муров, И. В. Невский температуре воздуха на входе в каналы tвх = 50 ºС необходимо обеспечить предельную температуру контактных площадок преобразователя (95 ºС) и транзистора (80 ºС).
Тепловыделяющие элементы
L=250 мм
H1 Ребра Планки
H=105 мм
Воздушный H2 канал
db 105 мм
Площадь оснований радиатора значительно превосходит площадь контакта локальных
источников теплоты. Самые мощные источники теплоты находятся в центре основания, по-
этому при расчетах появляется возможность применения алюминиевых ребер в крайних об-
ластях. В рассматриваемой конструкции радиатора толщина всех ребер одинакова, также
одинакова и толщина всех планок.
Искомыми параметрами являются: толщины двух оснований, ребер и планок между ни-
ми, количество медных ребер и планок, количество алюминиевых ребер и планок.
Как показали предварительные оценки, плотность теплового потока на основании 2 в 50 раз
меньше, чем на основании 1. Поэтому на начальном этапе расчетов используем сплав алюми-
ния для всех планок основания 2 и сплав меди — для планок основании 1. Все ребра медные.
Учитывая высокую плотность теплового потока, выбираем толщину ребра d = 1 мм. Началь-
ная толщина обоих оснований 8 мм.
Для вычисления коэффициента конвективного теплообмена ребер необходимо опреде-
лить скорость воздуха в каналах между ними. В установившемся режиме потери давления в ка-
налах радиатора (∆Pк) компенсированы давлением, создаваемым вентиляторами (∆Pв). Величи-
на ∆Pк соответствует точке пересечения двух кривых: напорной характеристики вентилятора
∆Pв=f1(G) и зависимости перепада давления в каналах от расхода воздуха G — ∆Pк=f2(G):
∆Pк
=
0,
25ρV
2
⎛⎜⎝1+
48νL b2V
⎞ ⎠⎟
,
где b — ширина канала (зазор между ребрами), L — длина воздушного канала, V — скорость
воздуха в канале, ρ — плотность воздуха, ν — кинематическая вязкость воздуха, G=VA —
расход воздуха во всех каналах, A — площадь проходного сечения.
При варьировании ширины канала b зависимость ∆Pк изменяется, вследствие чего из-
меняются и значения V и G, результаты расчета которых для трех значений b при d = 1 мм
приведены в табл. 1 (где п — количество ребер).
Таблица 1
b, мм
n V, м/с G, м3/ч
1
48 1,5
23
1,5 38 2,2 42
2
32 2,8
60
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
Оптимизация конструкции системы охлаждения электронного модуля
63
Для определения коэффициента теплообмена на поверхности ребра (αр) при L>40b используем формулы для „длинного канала“ [3]:
αр=4,12λв/b, где λв— теплопроводность воздуха.
Тогда тепловая проводимость ребра вычисляется как [4]
σp = L а проводимость всех n ребер — как
⎛ 2α р λмd ⋅ th ⎝⎜⎜ H
2αр λмd
⎞ ⎟⎟⎠
,
⎛ σрn = nLth ⎝⎜⎜ H
2αр λрd
⎞ ⎟⎟⎠
2αрλрd ,
где λм — теплопроводность медного ребра; Н — высота ребра; λp принимает значения теп-
лопроводности либо меди, либо алюминия.
Тепловая проводимость радиатора определяется выражением
σрад=W(1–exp(–σрп/W)), где W=Gνc — водяной эквивалент, здесь с — удельная теплоемкость воздуха.
Перегрев оснований радиатора (∆tосн) относительно температуры воздуха на входе канала (tвх) определяется соотношением
∆tосн =tосн–tвх=Ф/σрад,
где tосн — температура оснований, Φ — мощность теплового потока.
Результаты расчетов величин σр и σрп, σрад и ∆tосн ( при Φ = 395 Вт) приведены в табл. 2—4.
b, мм
1,5 2
Таблица 2
σр, Вт/К Медь Алюминий
1,84 1,43 1,56 1,23
Таблица 3
b, мм
n
σрп, Вт/К Медь Алюминий
1,5 38 69,9
54,3
2 32 49,9
39,4
b, мм
1,5 2
1,5 2
σрп, Вт/К
σрад, Вт/К
Медь
69,9 12,4
49,9 16,7
Алюминий
54,3 12,2
39,4 15,8
Таблица 4 ∆tосн, К
32,2 24,0
32,8 25,3
Наименьшее значение ∆tосн относительно tвх составляет 24 К для радиатора с медными ребрами и шириной канала 2 мм. Такое значение перегрева получено при равномерно распределенной мощности локальных источников теплоты. Учет реального расположения локальных источников теплоты производился с помощью программ расчета трехмерных температурных полей.
Для оптимизации конструктивных параметров радиатора был использован метод конечных элементов. Создана модель радиатора с локальными источниками теплоты на основаниях и возможностью варьирования толщины оснований, количества и физических свойств ребер. Модель состоит из 50 000 конечных элементов. На поверхности ребра задано постоянное значение коэффициента теплообмена αр, тогда переменная температура воздуха в каналах рассчитывается по формуле
l
tв (l, h) = tвх + ∫ Φ(l, h)dl (cρG) ,
0
где l и h — координаты соответственно вдоль воздушного потока в канале и перпендикулярно ему.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
64 Г. В. Бирюлин, В. И. Егоров, С. Ю. Муров, И. В. Невский
Предварительный анализ показал, что температурное поле, создаваемое локальными источниками теплоты на двух основаниях, удовлетворяет принципу суперпозиции, поэтому расчеты для двух оснований были проведены раздельно. Сначала была осуществлена оптимизация толщины основания. Затем была выбрана оптимальная ширина канала, после чего исследованы возможности замены меди на алюминиевый сплав в крайних ребрах и планках. По результатам расчетов установлено, что температурное поле на основании 1 неравномерно: в центре основания максимальный перегрев составляет 28 К, а на краях — 10 К. Температурное поле ребра в зоне действия мощного локального источника также неравномерно: у основания 1 перегрев составляет 18 К, а вблизи основания 2 — 5 К. Таким образом, учет расположения локальных источников теплоты с использованием трехмерной модели позволяет более точно найти конструктивные параметры радиатора и определить области, в которых теплопроводность материалов может быть уменьшена.
В результате получены следующие оптимальные параметры конструкции: толщина основания 1 (Н1 — см. рисунок) 10 мм; толщина основания 2 (Н2) 5 мм; толщина ребер 1 мм, ширина канала 1,5 мм; количество медных ребер — 26, алюминиевых — 12; количество медных планок — 27, алюминиевых — 51. Температура воздуха в зоне мощных транзисторов не превышает 79 ºС, а на входе в канал — 50 ºС. При замене в крайних областях до 30 % ребер и планок перегрев на контактных площадках транзисторов увеличился на 2 К, что позволило уменьшить массу радиатора на 25 % по сравнению с радиатором с медными ребрами.
Рассмотренный вариант оптимизации конструкции системы охлаждения на примере паяного биметаллического радиатора позволяет получить значительный выигрыш в массе для схем с интенсивными локальными источниками теплоты, однако для выбора оптимальной конструкции требуется численное моделирование температурных полей.
СПИСОК ЛИТЕРАТУРЫ
1. Муров Ю. М. Роль базовых несущих конструкций в конкурентоспособности российского приборостроения // Сб. материалов VI отраслевой науч.-техн. конф. „Стратегия развития радиоэлектронного комплекса России“, 12—14 сент. 2007 г. Нижн. Новгород, 2007.
2. Дульнев Г. Н. Тепло- и массообмен в радиоэлектронной аппаратуре. М.: Высш. школа, 1984.
3. Петухов Б. С. Теплообмен и сопротивление при ламинарном течении жидкости в трубах. М.: Энергия, 1967. 412 с.
4. Ройзен Л. И., Дулькин И. Н. Тепловой расчет оребренных поверхностей. М.: Энергия, 1977. 256 с.
Гавриил Владимирович Бирюлин
Владимир Иванович Егоров
Святослав Юрьевич Муров Игорь Владимирович Невский
Сведения об авторах — аспирант; Санкт-Петербургский национальный исследовательский
университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: gavrila@bk.ru — канд. техн. наук, доцент; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: KTF@grv.ifmo.ru — ОАО Бюро новой техники „Прибой“, Санкт-Петербург; инженер — студент; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: nevskigor@mail.ru
Рекомендована кафедрой компьютерной теплофизики и энергофизического мониторинга
Поступила в редакцию 06.07.11 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
61
УДК 62-712
Г. В. БИРЮЛИН, В. И. ЕГОРОВ, С. Ю. МУРОВ, И. В. НЕВСКИЙ
ОПТИМИЗАЦИЯ КОНСТРУКЦИИ СИСТЕМЫ ОХЛАЖДЕНИЯ ЭЛЕКТРОННОГО МОДУЛЯ
Рассматривается конструкция биметаллического (медь—алюминий) паяного радиатора, эффективно охлаждающего элементы электронной аппаратуры с высокой неравномерностью тепловых потоков. Представлена методика оптимизационного расчета конструктивных параметров радиатора.
Ключевые слова: принудительное воздушное охлаждение, радиатор, конечноэлементный метод, оптимизационный расчет.
Функционирование мощных электронных изделий малых размеров требует использования охлаждающих устройств, минимизация массы которых является важной и актуальной задачей. Для принудительного воздушного охлаждения применяются биметаллические ребристые радиаторы [1], что позволяет в областях с низкой плотностью теплового потока использовать металлы с малым удельным весом.
В современных электронных устройствах применяются изделия электронной техники с контактной площадью порядка 4·10–4 м2 и мощностью тепловыделений до 300 Вт. Площадь основания ребристых радиаторов, используемых для обеспечения требуемой температуры корпуса силового элемента, часто значительно превосходит контактную площадью элемента. В этом случае при большом значении коэффициента теплоотдачи оребренной поверхности локальным источником теплоты создается неравномерность температурного поля основания радиатора. В данных условиях проводить расчеты конструкции радиатора, основываясь на известных из литературы [2] значениях тепловых сопротивлений, недопустимо.
Плотности теплового потока в зоне расположения локального источника теплоты и в областях, удаленных от источника, существенно различаются. Поэтому в области, близкой к источнику теплоты, используют высокотеплопроводные материалы. На краях основания радиатора можно применить материал (сплав алюминия), имеющий меньшую плотность и теплопроводность по сравнению с медными сплавами.
Такое техническое решение можно реализовать, используя технологию создания паяных радиаторов, разработанную в БНТ „Прибой“ (Санкт-Петербург) [1]. Такие радиаторы состоят из спаянных между собой ребер и планок. Основание радиатора образуют планки и торцы ребер. Каждое ребро запаяно в два параллельных основания. Толщина основания есть толщина планки. Как планки, так и ребра могут быть медными и алюминиевыми в различных сочетаниях, например, по краям — из сплава алюминия, а в центре — из меди.
Рассмотрим вариант оптимизации конструкции системы охлаждения электронного модуля в целях минимизации его веса. Для обеспечения требуемого температурного режима электронных компонентов используется принудительное воздушное охлаждение. Несущими элементами конструкции (см. рисунок) являются пластины, образующие каналы прямоугольного сечения. На входе и выходе каналов (длиной L = 0,25 м) закреплены вентиляторы. На двух противоположных пластинах (основаниях радиатора) размещены тепловыделяющие элементы. Ширина оснований фиксирована. На основании 1 расположены два транзистора мощностью по 175 Вт, отстоящие на расстояние 10 мм. Удельная мощность тепловыделений каждого из них составляет 8·105 Вт/м2. Кроме того, на основании 1 размещены также три транзистора мощностью по 15 Вт. Размеры контактных площадок транзисторов 24×9 мм. На основании 2 расположены три преобразователя мощностью по 45 Вт, общее тепловыделение которых составляет 140 Вт. Размеры контактных площадок преобразователей 60×70 мм. При
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
62 Г. В. Бирюлин, В. И. Егоров, С. Ю. Муров, И. В. Невский температуре воздуха на входе в каналы tвх = 50 ºС необходимо обеспечить предельную температуру контактных площадок преобразователя (95 ºС) и транзистора (80 ºС).
Тепловыделяющие элементы
L=250 мм
H1 Ребра Планки
H=105 мм
Воздушный H2 канал
db 105 мм
Площадь оснований радиатора значительно превосходит площадь контакта локальных
источников теплоты. Самые мощные источники теплоты находятся в центре основания, по-
этому при расчетах появляется возможность применения алюминиевых ребер в крайних об-
ластях. В рассматриваемой конструкции радиатора толщина всех ребер одинакова, также
одинакова и толщина всех планок.
Искомыми параметрами являются: толщины двух оснований, ребер и планок между ни-
ми, количество медных ребер и планок, количество алюминиевых ребер и планок.
Как показали предварительные оценки, плотность теплового потока на основании 2 в 50 раз
меньше, чем на основании 1. Поэтому на начальном этапе расчетов используем сплав алюми-
ния для всех планок основания 2 и сплав меди — для планок основании 1. Все ребра медные.
Учитывая высокую плотность теплового потока, выбираем толщину ребра d = 1 мм. Началь-
ная толщина обоих оснований 8 мм.
Для вычисления коэффициента конвективного теплообмена ребер необходимо опреде-
лить скорость воздуха в каналах между ними. В установившемся режиме потери давления в ка-
налах радиатора (∆Pк) компенсированы давлением, создаваемым вентиляторами (∆Pв). Величи-
на ∆Pк соответствует точке пересечения двух кривых: напорной характеристики вентилятора
∆Pв=f1(G) и зависимости перепада давления в каналах от расхода воздуха G — ∆Pк=f2(G):
∆Pк
=
0,
25ρV
2
⎛⎜⎝1+
48νL b2V
⎞ ⎠⎟
,
где b — ширина канала (зазор между ребрами), L — длина воздушного канала, V — скорость
воздуха в канале, ρ — плотность воздуха, ν — кинематическая вязкость воздуха, G=VA —
расход воздуха во всех каналах, A — площадь проходного сечения.
При варьировании ширины канала b зависимость ∆Pк изменяется, вследствие чего из-
меняются и значения V и G, результаты расчета которых для трех значений b при d = 1 мм
приведены в табл. 1 (где п — количество ребер).
Таблица 1
b, мм
n V, м/с G, м3/ч
1
48 1,5
23
1,5 38 2,2 42
2
32 2,8
60
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
Оптимизация конструкции системы охлаждения электронного модуля
63
Для определения коэффициента теплообмена на поверхности ребра (αр) при L>40b используем формулы для „длинного канала“ [3]:
αр=4,12λв/b, где λв— теплопроводность воздуха.
Тогда тепловая проводимость ребра вычисляется как [4]
σp = L а проводимость всех n ребер — как
⎛ 2α р λмd ⋅ th ⎝⎜⎜ H
2αр λмd
⎞ ⎟⎟⎠
,
⎛ σрn = nLth ⎝⎜⎜ H
2αр λрd
⎞ ⎟⎟⎠
2αрλрd ,
где λм — теплопроводность медного ребра; Н — высота ребра; λp принимает значения теп-
лопроводности либо меди, либо алюминия.
Тепловая проводимость радиатора определяется выражением
σрад=W(1–exp(–σрп/W)), где W=Gνc — водяной эквивалент, здесь с — удельная теплоемкость воздуха.
Перегрев оснований радиатора (∆tосн) относительно температуры воздуха на входе канала (tвх) определяется соотношением
∆tосн =tосн–tвх=Ф/σрад,
где tосн — температура оснований, Φ — мощность теплового потока.
Результаты расчетов величин σр и σрп, σрад и ∆tосн ( при Φ = 395 Вт) приведены в табл. 2—4.
b, мм
1,5 2
Таблица 2
σр, Вт/К Медь Алюминий
1,84 1,43 1,56 1,23
Таблица 3
b, мм
n
σрп, Вт/К Медь Алюминий
1,5 38 69,9
54,3
2 32 49,9
39,4
b, мм
1,5 2
1,5 2
σрп, Вт/К
σрад, Вт/К
Медь
69,9 12,4
49,9 16,7
Алюминий
54,3 12,2
39,4 15,8
Таблица 4 ∆tосн, К
32,2 24,0
32,8 25,3
Наименьшее значение ∆tосн относительно tвх составляет 24 К для радиатора с медными ребрами и шириной канала 2 мм. Такое значение перегрева получено при равномерно распределенной мощности локальных источников теплоты. Учет реального расположения локальных источников теплоты производился с помощью программ расчета трехмерных температурных полей.
Для оптимизации конструктивных параметров радиатора был использован метод конечных элементов. Создана модель радиатора с локальными источниками теплоты на основаниях и возможностью варьирования толщины оснований, количества и физических свойств ребер. Модель состоит из 50 000 конечных элементов. На поверхности ребра задано постоянное значение коэффициента теплообмена αр, тогда переменная температура воздуха в каналах рассчитывается по формуле
l
tв (l, h) = tвх + ∫ Φ(l, h)dl (cρG) ,
0
где l и h — координаты соответственно вдоль воздушного потока в канале и перпендикулярно ему.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3
64 Г. В. Бирюлин, В. И. Егоров, С. Ю. Муров, И. В. Невский
Предварительный анализ показал, что температурное поле, создаваемое локальными источниками теплоты на двух основаниях, удовлетворяет принципу суперпозиции, поэтому расчеты для двух оснований были проведены раздельно. Сначала была осуществлена оптимизация толщины основания. Затем была выбрана оптимальная ширина канала, после чего исследованы возможности замены меди на алюминиевый сплав в крайних ребрах и планках. По результатам расчетов установлено, что температурное поле на основании 1 неравномерно: в центре основания максимальный перегрев составляет 28 К, а на краях — 10 К. Температурное поле ребра в зоне действия мощного локального источника также неравномерно: у основания 1 перегрев составляет 18 К, а вблизи основания 2 — 5 К. Таким образом, учет расположения локальных источников теплоты с использованием трехмерной модели позволяет более точно найти конструктивные параметры радиатора и определить области, в которых теплопроводность материалов может быть уменьшена.
В результате получены следующие оптимальные параметры конструкции: толщина основания 1 (Н1 — см. рисунок) 10 мм; толщина основания 2 (Н2) 5 мм; толщина ребер 1 мм, ширина канала 1,5 мм; количество медных ребер — 26, алюминиевых — 12; количество медных планок — 27, алюминиевых — 51. Температура воздуха в зоне мощных транзисторов не превышает 79 ºС, а на входе в канал — 50 ºС. При замене в крайних областях до 30 % ребер и планок перегрев на контактных площадках транзисторов увеличился на 2 К, что позволило уменьшить массу радиатора на 25 % по сравнению с радиатором с медными ребрами.
Рассмотренный вариант оптимизации конструкции системы охлаждения на примере паяного биметаллического радиатора позволяет получить значительный выигрыш в массе для схем с интенсивными локальными источниками теплоты, однако для выбора оптимальной конструкции требуется численное моделирование температурных полей.
СПИСОК ЛИТЕРАТУРЫ
1. Муров Ю. М. Роль базовых несущих конструкций в конкурентоспособности российского приборостроения // Сб. материалов VI отраслевой науч.-техн. конф. „Стратегия развития радиоэлектронного комплекса России“, 12—14 сент. 2007 г. Нижн. Новгород, 2007.
2. Дульнев Г. Н. Тепло- и массообмен в радиоэлектронной аппаратуре. М.: Высш. школа, 1984.
3. Петухов Б. С. Теплообмен и сопротивление при ламинарном течении жидкости в трубах. М.: Энергия, 1967. 412 с.
4. Ройзен Л. И., Дулькин И. Н. Тепловой расчет оребренных поверхностей. М.: Энергия, 1977. 256 с.
Гавриил Владимирович Бирюлин
Владимир Иванович Егоров
Святослав Юрьевич Муров Игорь Владимирович Невский
Сведения об авторах — аспирант; Санкт-Петербургский национальный исследовательский
университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: gavrila@bk.ru — канд. техн. наук, доцент; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: KTF@grv.ifmo.ru — ОАО Бюро новой техники „Прибой“, Санкт-Петербург; инженер — студент; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра компьютерной теплофизики и энергофизического мониторинга; E-mail: nevskigor@mail.ru
Рекомендована кафедрой компьютерной теплофизики и энергофизического мониторинга
Поступила в редакцию 06.07.11 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 3