CC BY
6Электронный журнал «Труды МАИ». Выпуск № 49
www.mai.ru/science/trudy/
УДК 536.33.614.841
Теплообмен в рабочем участке стенда с газоразрядными
источниками излучения
В.Н. Елисеев, В.А. Товстоног В.А., Мьо Тан Аннотация
Приведена математическая модель расчёта радиационного теплообмена в рабочем участке стенда с газоразрядными источниками излучения, позволяющая определять потоки излучений на поверхностях, образующих внутренний контур рабочего участка.
Ключевые слова
теплообмен; излучение; модель; расчёт; материал; конструкция; спектр; испытания; эксперимент.
Особенности радиационного теплообмена в рабочем участке стенда и математическая модель его расчёта
Многие вопросы, связанные с определением теплофизических свойств современных материалов и работоспособностью элементов конструкции летательных аппаратов в условиях интенсивного теплообмена, успешно решаются на установках и стендах радиационного нагрева [1-5]. На таких стендах удается нагревать образцы материалов и конструкций из них достаточно больших размеров, что исключает влияние краевых эффектов на тепловой режим объекта испытания (ОИ) (расслоение слоистого материала на его концах, выход газовых фракций через торцы и др.).
Достижимая плотность потока падающего излучения на стендах с галогенными
ттг «-» 2 „
лампами накаливания (ГЛН) ограничена величиной 100 -120 Вт/см2 [3]. Более высокий уровень плотности потока излучения обеспечивается на стендах с водоохлаждаемыми
тт 2
газоразрядными источниками излучения (ГИИ) (до 300 Вт/см2, [5]). Указанные источники обладают меньшей инерционностью, а величина нагреваемой поверхности зависит лишь от располагаемой мощности стенда.
Практическое использование ГИИ в целях тепловых испытаний связано с некоторыми особенностями, которые отличают их от ГЛН. По ряду причин, например, из-за
возможности возникновения пробоя между электродами, расположение ГИИ над поверхностью объекта испытаний отличается заметной дискретностью, что приводит к неравномерному распределению потока падающего излучения на нагреваемой поверхности. Существующие математические модели сопряженных задач радиационно-кондуктивного теплообмена (РКТ) применительно к установкам с ГИИ и алгоритмы их решения, как правило, не учитывают особенностей влияния спектральных характеристик излучения на температурное состояние нагреваемых объектов. Отсутствуют обоснованные критерии идентичности температурных полей в объектах испытаний, определяемые с использованием спектральных или интегральных характеристик. Отсутствует также какое-либо обоснование принципов выравнивания неравномерности температурного поля испытуемых объектов, связанное с дискретностью источников излучения, путем целенаправленного изменения оптических свойств нагреваемой поверхности.
Формулировка математической модели радиационного теплообмена, учитывающая реальные характеристики источников излучения и нагреваемого объекта, - необходимый этап подготовки тепловых испытаний с использованием источников излучения. Рассмотрим схему тепловых испытаний плоских объектов (рис.1), учитывая, что в приведенных схемах форма рефлектора может быть разная, а число источников излучения равно единице или кратно трем. Полагаем, что температура поверхности рефлектора 2 поддерживается постоянной, а температура поверхности объекта 4 есть функция времени и координат х, у. На внешней поверхности трубчатых ГИИ задана спектральная плотность потока излучения.
Для расчёта радиационного теплообмена воспользуемся зональным методом, в соответствии с которым все поверхности в системах на рис.1 представляются совокупностью элементарных площадок, между которыми происходит радиационный теплообмен через непоглощающую (диатермичную) среду. При этом учитывается наличие поверхностей ГИИ, затеняющих излучение между отдельными участками поверхностей системы.
Одна из основных задач расчёта теплообмена излучением в приведённых системах -определение результирующего потока излучения, обеспечивающего требуемое температурное состояние нагреваемого объекта. При этом расчёты радиационного и радиационно-кондуктивного теплообмена существенно упрощаются, если система, в которой рассматривается теплообмен, - замкнутая. В действительности зона нагрева может быть и незамкнутой в силу конструктивных или технологических причин (отверстия, зазоры между отдельными конструктивными элементами или открытые участки между нагревательным блоком и поверхностью объекта испытания).
Рис.1. Схема тепловых испытаний образцов материалов: а -система с тремя ГИИ; б - система с шестью ГИИ; 1- водоохлаждаемые ГИИ трубчатой формы; 2 - водоохлаждаемый рефлектор; 3 - условные поверхности, замыкающие системы а,б; 4 - испытуемый образец; 5 - термопары; 6- нагреваемая поверхность образца;
7 - тыльная поверхность образца.
Несмотря на это, каждую, даже незамкнутую систему, можно считать замкнутой (условно замкнутой, вводя фиктивные замыкающие поверхности), придавая всем поверхностям (в том числе и фиктивным) определенные оптические свойства и температуру (или тепловой поток). Тогда для каждой из этих поверхностей может быть рассчитано количество подводимой и отводимой энергии излучения. При таком подходе открытым участкам в условно замкнутых системах могут быть приписаны свойства абсолютно черной поверхности с нулевой температурой.
В задачах радиационного нагрева модель, позволяющая на каждом временном шаге определить потоки излучения на поверхностях нагреваемого объекта, представляет собой систему уравнений для определения спектральных потоков излучения, падающих на все выделенные элементы замкнутой геометрической системы (рис.1)
N N N
йпад (л) = е (л)ж + (л)ж = е (л) 6г (л)ж +
] =1 к=N. +1 ] =1
} 1 } (1)
N. N N У 7
+ (л)е;ад(л)ж, + еди(л)ж + еКк(Л)ОГ(лж ,
]=1 к=N+1 к+1
и спектрального потока результирующего излучения /-1 площадки
(л) = ^р- (л)¥г = А (л)е;ад (л) _ ^ (л) о^д) (2)
где
С F
о соб0(л) = а соб'° (л) f =-^-
v 7 v ) г л5(ес2/лТ) _ 1)
(3)
- спектральный поток собственного излучения /-го элемента площадью нагреваемого объекта, определяемый функцией Планка В. (л,Т ) = аС°б'°(Л), С = 3,742 • 10 _16Вт • м2, С = 0,01439м • К;
ои (Л) = аи (Л) ^ (4)
- спектральный поток излучения участка поверхности Гк , расположенного на источнике излучения; аИ (л) - спектральная плотность потока излучения с внешней поверхности источника; _диффузные угловые коэффициенты, определяемые с учётом возможного затенения одних площадок другими; ^(л) и ¿>(л) _ спектральная
отражательная и излучательная способности у'-го и к-го участков поверхности; N - число участков поверхности в замкнутой системе, на которых заданы их температуры; N2 - число участков, на которых заданы потоки их собственного излучения.
В отличие от системы уравнений, предназначенных для расчёта радиационного теплообмена в рабочем участке с ГЛН, система уравнений (1) не содержит слагаемого, учитывающего собственное излучение их водоохлаждаемых оболочек из кварцевого стекла.
Площадки, расположенные на водоохлаждаемом рефлекторе и ограждающих поверхностях, хотя и входят в число площадок Ы\, имеют, как правило, сравнительно низкую температуру и вносят малый вклад в общий баланс потока излучения.
При использовании для решения системы уравнений (1) - (4) приближения спектральных полос, найденные низ них значения 0гпад(Лт) и ОР^(Лт) позволяют определить
интегральные потоки падающего и результирующего излучений.
м
опад = у опад (л )дл ; (5)
т=1
м
орез = у оРез (л )лл , (6)
^г V т / т 5 ^ '
т=1
где М - число полос в спектре излучения: ЛЛт - ширина спектральной полосы, приведенная к длине волны Лт .
Собрав в выражении (1) слагаемые, содержащие неизвестные потоки падающего излучения в левой части, а известные величины в правой, запишем его в виде
N. N
дг(л) - 2я} (лм,дпад(л) - 2як (лк0г(л) =
1=1
к=N +1
(7)
=Е ^ (л) дгчлр, + е ди (Л)9к,г
к=Ы1+1
1=1
Обозначив для упрощения записи д™д (л) как д , известные величины в системе уравнений (7) через сг и коэффициенты перед д , через а , систему (7) запишем в
компактном виде
N
а
^а д = с.,
(8)
1
где
а, 1 = 81 - я (л)р = 8 1 - [1 - а (лр для 1 < i < ы15
а
-я (л)р , =-[1 - а (л) - ^ (лр, для (N +1)
+1) < I < N;
(9)
_ [1 при 1 = 1 1 10 при I Ф1
-8 - символ Кронекера;
N¡
2^ (л)д;06 0 (л)р.. для 1 </<n
1=1
N
2 е; (л) Рк,, для (n1 +1)< i n
с = 4
(10)
к=N +1
Метод решения системы линейных алгебраических уравнений (8)-(10) хорошо известен и описан в литературе, в частности, в работе [6].
Определив спектральные потоки падающего на г-е площадки излучения, которые в приближении спектральных полос соответствуют потокам на средней для диапазона лл
длине волны Лт, т.е. д = д™д (Ли), находят спектральные потоки результирующего
излучения дрез (Ли ) для любой площадки системы по формуле (2).
В завершение расчёта радиационного теплообмена в замкнутой системе, на части которой (1 < г < N) заданы температуры, а на остальных (N + 1) < г < М) потоки излучения источников, по формулам (5) и (6) могут быть вычислены интегральные
епад /ОРез
и о .
После определения потоков результирующего излучения на этом же временном шаге решается задача нагревания объекта испытания. С этой целью формируется краевая задача теплопроводности, в которой в качестве одного из граничных условий используется найденное выше значение результирующего потока излучения на фронтальной поверхности объекта.
Процесс вычисления имеет итерационный характер и осуществляется с помощью существующих программ, основанных на использовании метода конечных элементов. В ряде случаев недостаточно определить только потоки падающего и результирующего излучений и требуется выделить отдельно вклад в общий баланс этих потоков, поток, поступающий на ОИ, например, только от рефлектора с целью выбора его геометрических параметров (формы, размера), обеспечивающих наилучшее, в смысле заданного критерия, решение задачи моделирования теплообмена на поверхности ОИ. Один из подходов к решению такой задачи и алгоритму ее решения приведен в работах [7,8].
Физическая модель газоразрядного источника излучения
Решение приведённой выше системы уравнений радиационного теплообмена в рабочем участке стенда с ГИИ (уравнения 8-10) предполагают, что известны ответы на два следующих вопроса. Первый - каким образом задать в качестве известного параметра поток излучения ГИИ на его внешней поверхности? Второй - какова прозрачность источника для его же излучения, идущего от рефлектора к ОИ? Учитывая достаточно большие размеры поперечного сечения водоохлаждаемых ГИИ (20-22 мм) при расстоянии между ними 40-50 мм, важность ответа на второй вопрос очевидна, поскольку они могут создавать заметное затенение излучения в рассматриваемой геометрической системе, а это неизбежно отразится на алгоритме вычисления угловых коэффициентов р , входящих в систему уравнений (810).
Характеристики водоохлаждаемых ГИИ, полученные на основе экспериментов, включая и зависимость плотности потока их излучения, приведены в работе [5].
Результаты измерений плотности потока излучения на расстоянии 21 мм от оси источника в его центральном сечении в непрерывном режиме работы ГИИ представлены в виде следующих аппроксимационных формул
Тип ГИИ ДТП 10/200 ДТП 10/220 ДК с ТВ 6000
Плотность потока излучения, (и),
= 288 - 4,77и + 0,0243и2 = 459 - 6,09и + 0,0235и2 = 213 - 2,94и + 0,011и2
кВт
м
Среднеквадратичная погрешность аппроксимации не превышает 5%. Общий вид и схема устройства ГИИ показаны на рис.2.
Рис.2. Общий вид и схема водоохлаждаемого трубчатого источника излучения
Излучение от разрядного канала плазмы проходит слой воды и стекла, частично поглощаясь в них. Излучение, вышедшее из источника, может вновь падать на его поверхность после отражения от рефлектора и других поверхностей системы, образующих внутренний контур рабочего участка. В связи с этим для оценки возможного затенения излучения источником необходимо принять допущение о степени прозрачности ГИИ для излучения, падающего на его внешнюю поверхность. С этой целью выполнялись расчёты радиационного теплообмена в рабочих участках с шестью расположенными в ряд ГИИ и рефлекторами различной формы [10]. Расчёты проводились для трёх моделей источников излучения: непрозрачные ГИИ, частично прозрачные и полностью прозрачные ГИИ. В модели частично прозрачного ГИИ варьировали диаметр его непрозрачной части. Результаты расчёта сопоставлялись с экспериментом.
На рис.3 представлены, в качестве примера, результаты расчёта и эксперимента для рефлектора в форме трапеции (рис.4).
2
Анализ серийных результатов расчёта с ГИИ, диаметр внешней оболочки которых составлял 20 мм, а мощность изменялась в диапазоне от 40 до 130 КВт, показал, что наилучшим соответствием эксперименту обладает модель частично прозрачного ГИИ с диаметром разрядного канала 5,5 ...6,5 мм. Именно в этом интервале выбирался в дальнейшем диаметр затеняющего цилиндра при расчёте угловых коэффициентов.
Рис.3. Распределение плотности суммарного падающего излучения по ширине (х2) испытуемого образца для трёх моделей ГИИ: 1 - полностью прозрачные ГИИ; 2 - частично прозрачные ГИИ; 3 - непрозрачные ГИИ; ■ - эксперимент.
Рис.4. Схема рабочего участка с трапецевидным рефлектором.
Другой подход к определению потока излучения с поверхности одиночного ГИИ и оценке его теплового режима основан на использовании понятия внутреннего коэффициента полезного действия (КПД)
Яо (и) F
77Ш) = ^ } внут, (11)
/V у р(ц) , ( )
где (и) - плотность потока излучения на внутренней поверхности горячей оболочки
(горелки) ГИИ - ,Рвнут; Р(и) - электрическая мощность источника излучения, зависящая от величины питающего напряжения и.
Исследования излучательных характеристик газоразрядных источников излучения, выполненные В.М. Градовым, в частности, в работе [9] позволили определить интегральную плотность потока излучения и потока , передаваемого через внутреннюю поверхность
кварцевой оболочки путём теплопроводности, их сумму, внутренний КПД, а также температуру плазменного разряда в его центре и спектр источника. Для источника излучения (ГИИ) ДТП 10/200 его основные характеристики приведены в табл.1
Таблица 1
Р, Вт о д0, Вт/см2 о дт, Вт/см2 о дЕ, Вт/см2 0 Р Тшах, К
128 2674 226,5 2900,5 0,92 8885
59 1210 133,2 1343,2 0,9 8220
Зависимость внутреннего КПД источника ДТП 10/200 с внутренним радиусом разрядного канала 3,5 мм от электрической мощности приведена на рис. 5.
КПД 0,92
0,91 0,90 0,89 0,88 0,87
40 50 60 70 80 90 100 110 120
Р(кВт)
Рис.5. Зависимость внутреннего КПД источника ДТП 10/200 от его мощности
Полагая, что вся срабатываемая в ГИИ электрическая мощность в конечном итоге передается его внутренней оболочке путём излучения и теплопроводности от плазмообразующего газа Р = (а. + а ) Ь ,
V 10 1т у внут7
можно для заданной мощности Р определить плотность потоков энергии
Рл( Р)
ао = Р^(Р); (12)
ь
внутр
[1 - л(Р)]Р
ат =1 /( ^ . (13)
ь
внутр
Величины плотностей потоков излучения ао и теплопроводности ат на внутренней поверхности оболочки ГИИ используются в граничных условиях при расчёте её температурного состояния и служат исходной информацией при выборе режима охлаждения ГИИ [10].
Расчёты радиационного теплообмена с использованием интегральных характеристик излучения экономнее по затратам времени, но требуют соблюдения определенных условий, обеспечивающих идентичность температурных полей, найденных тем или другим методом. Достаточным условием идентичности температурных полей в нагреваемом образце при прочих одинаковых условиях теплообмена может служить равенство потоков результирующего излучения на нагреваемой поверхности, рассчитанных с использованием
спектральных (О^ ) и интегральных (Ор^ ) характеристик излучения
(*)=о:(*)> *=
или равенства их плотностей, если Ь = Ь = •••Ь
а::, (?) = (*),
где п - число участков на нагреваемой поверхности ОИ.
Рассмотренная модель теплообмена в рабочем участке стенда с водоохлаждаемыми газоразрядными источниками излучения позволяет на этапе подготовки физического эксперимента:
- определить необходимые для расчёта радиационного теплообмена исходные данные в виде зависимостей плотности потока излучения на внешней поверхности ГИИ и его внутреннего КПД от питающего напряжения;
- выбрать величину «затеняющего» диаметра источника излучения, необходимого для расчёта угловых коэффициентов;
- выполнить расчёт падающего и результирующего излучений на всех поверхностях, образующих внутренний контур рабочего участка, с учётом их спектральных характеристик.
Библиографический список
1.Гофин М.Я. Жаростойкие и теплозащитные конструкции многоразовых аэрокосмических аппаратов. - М. ЗАО «ТФ Мир», 2003.-671 с.
2. Баранов А.Н. Теплопрочностные испытания натурных конструкций летательных аппаратов. Экспериментальные исследования, Т.3.- М. Изд-во МГТУ им. Н.Э. Баумана, 2002, - С.156-198.
3. Синярев Г.Б. Общие принципы и организация стендовых тепловых испытаний // Вопросы теплообмена и тепловых испытаний конструкций. - 1982. - Труды МВТУ . - № 392. - С.3-16.
4. Елисеев В.Н., Товстоног В.А., Чирин К.В., Селезенев В.А. Установка лучистого нагрева широкого применения // Вестник БССР. Сер. Физико-энергетических наук.- 1990. -№ 2. - С.93-97.
5. Елисеев В.Н., Товстоног В.А. Характеристики источников излучения и излучательных систем высокоинтенсивного нагрева // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2001. - № 4. - С.3-32.
6. Зигель Р., Хауэлл Дж. Теплообмен излучением. - Мир, 1975. - 934 с.
7.Мьо Тан, Баслык К.П., Товстоног В.А., Елисеев В.Н. Алгоритм расчёта радиационно-кондуктивного теплообмена в установках с газоразрядными источниками излучения // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2007. - № 4. - С.33-46.
8. Мьо Тан. О выборе рациональной формы рефлектора для нагревателя с газоразрядными источниками излучения // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2. -С.117-120.
9. Градов В.М., Щербаков А.А. Расчёт излучательных характеристик дуговых и ксеноновых разрядов // Оптика и спектроскопия. - 1979. - т.47, № 4. - С.655-642.
10. Мьо Тан. Разработка методического и алгоритмического обеспечения тепловых испытаний материалов и элементов конструкции в стендах с газоразрядными источниками излучения. Дисс.. .канд.техн. наук.- М., 2008. - 181 с.
Сведения об авторах
Елисеев Виктор Николаевич, профессор Московского государственного технического университета имени Н.Э. Баумана, д.т.н. 2-ая Бауманская ул., дом 5, Москва, 105005; ешаП: у.п.еНвееу @ gmail.com
Товстоног Валерий Алексеевич, профессор Московского государственного технического университета имени Н.Э. Баумана, д.т.н. 2-ая Бауманская ул., дом 5, Москва, 105005; ешаП:1;оу81;опу@шаП.ги
Мьо Тан, выпускник Московского государственного технического университета имени Н.Э. Баумана, к.т.н.
2-ая Бауманская ул., дом 5, Москва, 105005; ешаП:1;оу81;опу@шаП.ги
