ТЕПЛООБМЕН И ТРЕНИЕ В КАНАЛАХ СУДОВЫХ ГАЗОТУРБИННЫХ ЭНЕРГЕТИЧЕСКИХ УСТАНОВОК ПРИ ИЗМЕНЕНИИ ДИНАМИКИ УВЕЛИЧЕНИЯ ТЕМПЕРАТУРЫ РАБОЧЕГО ТЕЛА
Аннотация и ключевые слова
Аннотация (русский):
Представлены результаты экспериментального исследования по изучению нестационарных эффектов и влияния нестационарности, вызванной резким увеличением температуры газового потока, на значения локальных коэффициентов трения и теплоотдачи в области начального участка цилиндрического канала. Причиной такого разрушения двигателя является локальный нагрев отдельных элементов, т. е. процесс теплоотдачи, происходящий за пределами камеры сгорания в жаровых трубах, где в нестационарных условиях на входе в турбину формируется вращающий её нагретый газовый поток. Кинематическая и тепловая структуры такого газового потока в рассматриваемый период времени испытывают на себе одновременное воздействие ряда дестабилизирующих факторов: неизотермичности, продольного градиента давления в условиях проявления эффектов динамической и тепловой нестационарности. Коэффициент трения Cf и коэффициент теплоотдачи St , являясь одними из важнейших параметров потока, характеризуют газодинамику течения потока и теплообменные процессы, протекающие в нём. Нестационарность существенно влияет на величины коэффициентов трения и теплоотдачи, которые в 2,0-2,5 раза увеличиваются или уменьшаются относительно своих квазистационарных аналогов. За пределами нестационарного процесса, когда температура рабочего тела становится постоянной, определяющую роль играет температурный напор, при наличии которого возникает явление ламинаризации теплового турбулентного пограничного слоя.

Ключевые слова:
газовый поток, энергетическая установка, коэффициент трения, коэффициент теплоотдачи, нестационарность, нестационарные эффекты, неизотермичность
Текст
При эксплуатации газотурбинных двигателей выделяют четыре основных режима. Пусковому режиму отводится роль одного из наиболее важных и ответственных, поскольку за короткий промежуток времени в результате горения топлива происходит резкое повышение температуры рабочего тела и, как следствие, такие же резкие по величине изменения всех параметров двигателя, приводящие к повышению температуры обтекаемых поверхностей конструктивных элементов двигателя. В пусковом режиме двигателя возникают ситуации, когда происходит преждевременный отказ работы двигателя с последующим разрушением конструктивных элементов. Рассмотрим совместное влияние неизотермичности, эффектов нестационарности и продольного градиента давления, проявляющегося формированием турбулентного пограничного слоя на начальном участке цилиндрического канала, в условиях различной по величине динамики возрастания температуры газового потока. Экспериментальное исследование Эксперименты были выполнены на экспериментальном стенде, который представляет собой газодинамическую трубу разомкнутого типа (рис. 1). В её состав входят форкамера 2 с выравнивающими решётками 3 [1], на выходе которой установлено спрофилированное по кривой сопло Витошинского. К выходному фланцу сопла болтами крепится экспериментальный канал 4. Все перечисленные элементы газодинамической трубы изготовлены из нержавеющей стали 1Х18Н9Т. Выравнивающие решётки 3 совместно с соплом обеспечивают равномерное распределение скоростей и температур на входе в экспериментальный канал с коэффициентами искажения w1cp/w1max и T1cp/T1max не менее 0,97. Нагрев рабочего тела, притекающего в форкамеру, производился плазмотроном 1, изготовленным по однокамерной схеме, в котором для стабилизации дуги выполняется её вихревая крутка. Анод и катод плазмотрона изготовлены из меди и по резьбе вкручиваются с двух сторон в обойму из органического стекла. Обойма имеет тангенциальное отверстие, через него в рабочую камеру плазмотрона подаётся сжатый воздух. Через выходной фланец плазмотрон стыкуется с форкамерой. Экспериментальный канал - это имитационная модель жарового патрубка, представляющая собой цилиндрическую трубу диаметром D = 45 мм и толщиной стенок 0,08 мм. Труба собрана из семи взаимозаменяемых секций длиной D. В середине 1, 3, 5 и 7-й секций были организованы контрольные сечения, в которые смонтированы хромель-копелевые термопары 6 диаметром 40 мкм для измерения температуры стенок Тw. В тех же контрольных сечениях 3, 5 и 7-й секций были смонтированы отборные устройства для измерения распределения по длине канала динамического напора на его оси, статического давления Δрст и пристенных касательных напряжений трения τw [2, 3]. Перепады давлений в этих сечениях измерялись индуктивными датчиками давления. Температура измерялась микротермопарами, а перепады давлений - малоинерционными индуктивными датчиками ДМИ-0,1 в комплекте с вторичной аппаратурой УГ-УМ. Кроме этого, хромельалюмелевой термопарой 7 диаметром 40 мкм на входе в экспериментальный канал менялась температура ТО, которая изменялась в диапазоне 293-1 500 К с градиентом температуры до 12 000 К/с. Полученная информация аналогового ввода ADAM 4019 и RS-232/485 поступала на персональный компьютер (РС). Рис. 1. Схема экспериментального стенда Перераспределение подачи холодного воздуха в форкамеру, минуя плазмотрон, позволило получать необходимые по величине температурные режимы. Присутствие аргона в зоне дуги плазмотрона существенно снижает эрозию катода. Форкамера конструктино имеет радиальные отверстия, куда воздух поступает через регулировочный кран 10 для снижения рабочей температуры газа и обеспечения необходимого теплового диапазона исследуемых параметров. Для стабилизации дуги в плазмотрон из баллона через редуктор 11 и расходомер 8 подаётся аргон. Катод и анод плазмотрона изготовлены из меди М1 и установлены в корпус из органического стекла. Для этого с воздухом смешивался аргон - 1,5-2,0 % от величины расхода воздуха. Расходы воздуха и аргона измерялись расходомерами 8. Измерительные системы были метрологически обследованы [3-5], в результате чего установлено, что они во всём диапазоне измерений не имеют искажений амплитуды и фазы. С целью оценки точности выполненных измерений определены среднеквадратичные погрешности для температуры и скорости потока, а также коэффициентов трения и теплоотдачи, которые не превысили значений 0,75, 1,6, 9,2 и 9,5 % соответственно. Результаты экспериментальных исследований Переменность величины того или иного параметра во времени формирует нестационарность протекающего процесса. В нашем случае в пусковом режиме при включении плазмотрона резко повышается температура газового потока T0 на входе в экспериментальный канал. Результатом этого является возрастание температуры обтекаемой поверхности (т. е. стенки канала), которая вниз по направлению течения снижается из-за увеличения толщины пограничного слоя. Нестационарные эффекты. Первым и главным источником нестационарности является повышающаяся температура газового потока T0 (рис. 2) на входе в экспериментальный канал. В это время в газовом потоке с изменением температуры изменяются по величине физические свойства рабочего тела, т. е. интенсивно уменьшается плотность ρ0 и увеличивается вязкость μ0, в результате скорость потока w0 вне пограничного слоя увеличивается. В экспериментах температура T0 и скорость w0 потока увеличились более чем в 3 раза, и производные по времени составили величины температуры dT0 / dt = 12 000 К/с и скорости м/с2 (рис. 2). Одновременное присутствие временных производных температуры и скорости потока свидетельствует об одновременном проявлении эффектов тепловой и динамической нестационарности с похожим характером изменения величин производных во времени. В рассматриваемой ситуации число Re (рис. 2) уменьшается, что стимулирует нестационарные процессы [6-10]. Рис. 2. Изменения начальных условий во времени Рассматриваемый процесс интенсивного повышения температуры потока T0 можно разделить на два интервала. В первом, продолжительностью от 0 до 0,15 с (рис. 2), интенсивно повышается температура потока T0 и почти не изменяется по величине температура стенки канала Тw. Во втором, температура потока T0 принимает практически постоянную величину при нарастающем прогреве стенок канала. Переменность во времени величин температур T0 и Тw указывает на присутствие нестационарности в обоих временных интервалах. В первом рассматриваемом интервале одновременно с температурой T0 увеличивается скорость потока w0 с ускорением до 700 м/с2, т. е. кроме тепловой нестационарности присутствует и динамическая. В состав параметра динамической нестационарности z, кроме коэффициента трения Сf и толщины динамического пограничного слоя, входит комплекс определяющий характер его изменения и достигающий в этом интервале максимального значения. На начальном участке канала формируется пограничный слой с толщиной, увеличивающейся вниз по направлению течения. Поэтому рассматриваемый участок течения находится под небольшим воздействием продольного градиента давления, которое характеризуется производной скорости потока w0 по продольной координате, увеличивающейся во времени более чем в 4 раза. Параметр продольного градиента давления определяет комплекс и его величина во времени почти не изменяется. Поскольку производная изменяется по величине во времени, то это ничто иное, как ещё один вид нестационарности по производной . В то же время присутствует вторая производная (рис. 2). Температурный фактор в этом временном интервале уменьшается от 1 до 0,25, и, соответственно, с участием производной является источником ещё одного нестационарного эффекта. Примечательным является и тот факт, что все временные производные принимают своё максимальное значение на одной временной отметке. Параметр тепловой нестационарности формируется на основе производной Величины обоих параметров - тепловой и динамической нестационарности - имеют одинаковый характер изменения во времени с той лишь разницей, что параметр тепловой нестационарности достигает максимального значения быстрее, чем параметр динамической нестационарности (рис. 2). Анализ совместного развития комплексов параметров отрицательного продольного градиента давления и динамической нестационарности z показал, что первый из параметров увеличивается вниз по направлению течения, а второй в том же направлении уменьшается. Данный факт свидетельствует о том, что даже в цилиндрическом канале по мере нарастания пограничного слоя отрицательный продольный градиент давления стимулирует уменьшение величины параметра динамической нестационарности и ослабление эффектов нестационарности. Во втором временном интервале температура потока T0 достигает своей постоянной величины. Значения всех временных производных обретают нулевую величину, что говорит об отсутствии эффектов нестационарности. С другой стороны, здесь интенсивно прогревается обтекаемая поверхность с повышающейся температурой стенки канала Тw. Температурный фактор, достигнув минимального значения, меняет направление изменения своей величины в сторону единицы. Можно сделать предположение о возможности проявления нестационарности за счёт изменения температуры стенки канала Тw во времени. Коэффициент трения. Когда параметр динамической нестационарности принимает отрицательные значения, а температурный фактор ѱh - величину меньше 1 (как в наших экспериментах), то поток ускоряется [6-9]. Такая ситуация складывается в первом временном интервале. Здесь пристенные касательные напряжения трения τw увеличиваются и в момент времени t = 0,04 с принимают максимальное значение. Все производные по времени в этот же момент времени принимают тоже максимальные значения, и, следовательно, все эффекты и параметры нестационарности имеют также максимальную величину. В рассматриваемом интервале времени по направлению течения пристенные касательные напряжения трения уменьшаются, что показывает зависимость коэффициента трения Сf в функции характерного числа Re (рис. 3). Рис. 3. Зависимость коэффициента трения от числа Re Характер изменения величины коэффициента трения аналогичен характеру изменения пристенных касательных напряжений трения τw, т. е. коэффициент трения Сf во всех контрольных сечениях монотонно увеличивается до максимума, а затем снижается. На графике построена прямая, представляющая собой «стандартный» закон трения [11]: (1) Все опытные точки расположились выше «стандартного» закона трения. К концу рассматриваемого первого временного интервала, начиная с отметки времени t = 0,15 с, величина коэффициента трения не изменяется, т. е. нестационарный процесс закончился. К этому моменту времени сформировалась величина температурного напора около 700 К. Если учесть влияние неизотермичности, то за пределами временной отметки t > 0,15 с опытные точки группируются около «стандартной» зависимости. Обращает на себя внимание также и тот факт, что с увеличением характерного числа Re** коэффициент трения Сf становится менее чувствительным к воздействию дестабилизирующих факторов: нестационарности и продольному градиенту давления. Во временном интервале от 0 до 0,15 с активного увеличения коэффициента трения Сf экспериментальные результаты удалось аппроксимировать зависимостью с погрешностью 10 %: (2) Коэффициент теплоотдачи. На рис. 4 показана зависимость коэффициента теплоотдачи St в функции числа а также построена прямая 1, представляющая собой «стандартный» закон для турбулентных режимов течения. В первом временном интервале в условиях активного возрастания температуры потока T0 коэффициент теплоотдачи St больше своего квазистационарного аналога (прямая 1) только в интервале времени до 0,06 с. Рис. 4. Зависимость числа St от числа а - Х = 2,5; б - Х = 4,5; в - Х = 6,5; 1 - ; 2 - ; 3 - St максимально превышает величину аналога в 2,0-2,5 раза. В работах [8-9] показано, что динамическая и тепловая нестационарность взаимно противоположным образом воздействуют на процесс теплоотдачи. Так как первопричиной рассматриваемого нестационарного процесса является изменение температуры T0 и, соответственно, влияние сформировавшейся тепловой нестационарности, то в начальные мгновения времени коэффициент теплоотдачи St во всех контрольных сечениях канала, начиная с максимального значения, уменьшается. Причинами снижения числа St являются динамическая нестационарность и неизотермичность, приводящие к ускорению потока и формированию более устойчивой кинематической его структуры к внешним воздействиям. В этом интервале времени достигается максимальная величина температурного напора при постоянной скорости прогрева обтекаемой поверхности. Процесс снижения величины коэффициента теплоотдачи St происходит в каждом контрольном сечении канала по параллельным лучам 3 (рис. 4), которые представлены аппроксимационной зависимостью вида (3) За пределами первого временного интервала t > 0,15 с нестационарные эффекты отсутствуют. Процесс теплообмена происходит при температурном напоре около 900 К. Экспериментальные точки (рис. 4) сгруппировались около прямой линии 2, которая является «стандартным» законом теплоотдачи при ламинарных режимах течения. Число Re (по нему определяют режим течения), построенное по среднерасходной скорости, равно 18 800. При таком значении его режим течения является турбулентным. Таким образом, во втором временном интервале экспериментально получено явление ламинаризации теплового турбулентного пограничного слоя [12] при тепловом потоке, имеющем направление от газа к стенке. Выводы В пусковом режиме газотурбинной энергетической установки происходит резкое увеличение температуры газового потока и возникает ряд нестационарных эффектов, в результате действия которых формируются тепловая и динамическая нестационарности. В рассмотренном случае, как показали эксперименты, несмотря на турбулентный режим течения рабочего тела, коэффициенты теплоотдачи должны рассчитываться по зависимостям для ламинарных режимов течения. Динамический пограничный слой в рассмотренном пусковом режиме не испытывает на себе аномальных метаморфоз, и поэтому величину коэффициента трения Сf следует рассчитывать, как обычно, по зависимостям для турбулентных режимов движения рабочих сред.
Список литературы

1. Повх И. Л. Аэродинамический эксперимент в машиностроении. Л.: Машиностроение. Ленинград. отд-ние, 1974. 478 с.

2. Репик Е. У., Кузенков В. К. Исследование нового метода опытного определения поверхностного трения в турбулентном пограничном слое // Инженерно-физический журнал. 1980. Т. 38. № 2. С. 197-200.

3. Володин Ю. Г., Марфина О. П., Богданов А. Н., Цветкович М. С., Кузнецов А. Б. Измерение касательных напряжений трения в нестационарном газовом потоке // Датчики и системы. 2009. № 2 (117). С. 34-36.

4. Никифоров А. Н., Фафурин А. В., Фесенко С. С., Хуснутдинов Ш. Н. Исследование динамических характеристик пневмометрических приемников // Труда метрологических институтов СССР. М.: Изд-во стандартов, 1977. Вып. 182 (242). С. 84-88.

5. Володин Ю. Г., Гильфанов К. Х., Марфина О. П., Закиров И. Ф., Казаков А. А., Кузнецов А. Б., Рыжакова Ж. С. Экспериментальное исследование тепловой инерционности микротермопар // Приборы. 2008. № 4. С. 52-55.

6. Володин Ю. Г., Марфина О. П. Расчёт коэффициентов трения и теплоотдачи при нестационарном неизотермическом течении несжимаемого газа в осесимметричных каналах // Изв. высш. учеб. заведений. Машиностроение. 2007. № 3. С. 21-26.

7. Володин Ю. Г., Марфина О. П. Математическое моделирование пусковых режимов энергетических установок. СПб.: Инфо-да, 2007. 128 с.

8. Володин Ю. Г., Федоров К. С., Яковлев М. В. Нестационарные эффекты и трение при запуске энергетических установок // Изв. высш. учеб. заведений. Авиационная техника. 2006. № 1. С. 34-36.

9. Володин Ю. Г., Федоров К. С., Яковлев М. В. Нестационарные эффекты и теплообмен в пусковом режиме энергетических установок // Изв. высш. учеб. заведений. Авиационная техника. 2006. № 4. С. 41-44.

10. Володин Ю. Г., Марфина О. П. Границы применения математической модели нестационарного течения несжимаемого газа в осесимметричных каналах // Вестн. Технолог. ун-та. 2016. Т. 19. № 6. С. 130-131.

11. Кутателадзе С. С., Леонтьев А. И. Тепломассообмен и трение в турбулентном пограничном слое. М.: Энергия, 1972. 320 с.

12. Володин Ю. Г., Кирпичников А. В. Ламинаризация турбулентного пограничного слоя при пуске энергетической установки и параметры ламинаризации // Вестн. Технолог. ун-та. 2016. Т. 19. № 17. С. 122-126.


Войти или Создать
* Забыли пароль?