HEAT TRANSFER AND FRICTION IN CHANNELS OF MARINE GAS TURBINE POWER PLANTS UNDER CHANGING DYNAMICS OF INCREASING WORKING MEDIUM TEMPERATURE
Authors:
1.
Kazan State University of Architecture and Engineering
2.
Volga State University of Water Transport
(Candidate of Technical Sciences, Assistant Professor; Assistant Professor of the Department of Ship Power Plant Operation)
3.
Caspian Institute of Sea and River Transport, Branch of Volga Academy of Water Transport
Type:
Article
Pages:
from 50
to 57
Status:
Published
Received:
20.08.2018
Accepted:
25.08.2018
Published:
25.08.2018
Subject area:
UDK 536.24+532.52
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
Language:
Russian
Keywords:
gas flow, power plant, friction coefficient, coefficient of heat losses, non-stationarity, non-stationary effects, non-isothermality
Abstract (English):
He article presents the results of the experimental research of non-stationary effects and the influence of transiency caused by burst heating of the gas flow on values of local coefficients of friction and heat emission on the initial section of the cylindrical channel. The main reason of destruction of the main engine is local heating of some elements, i.e. the process of heat-transfer that occurs outside the combustion chamber in fire tubes, where at the entrance of turbine generates heated gas flow in non-isothermal environment. This heated gas flow rotates the turbine. Kinematic and heat structures of this gas flow at the studied period of time undergo the mutual effect of destabilizing factors, such as non-isothermal effect of the longitudinal pressure gradient caused by the uprising of the dynamic and heat transiency. Friction coefficient Cf and coefficient of heat losses St are the most important criteria of the gas flow, which characterize gas dynamics and heat exchange of the gas flow. Non-stationarity has an essential effect on the values of friction coefficient and heat loss coefficient, which increase or decrease in 2.0-2.5 times, compared to their quasisteady analogues. Beyond non-stationary process, when the temperature of working body becomes stationary, the key role is played by the thermal head which causes laminarization phenomenon of the of the heat turbulent boundary layer.
He article presents the results of the experimental research of non-stationary effects and the influence of transiency caused by burst heating of the gas flow on values of local coefficients of friction and heat emission on the initial section of the cylindrical channel. The main reason of destruction of the main engine is local heating of some elements, i.e. the process of heat-transfer that occurs outside the combustion chamber in fire tubes, where at the entrance of turbine generates heated gas flow in non-isothermal environment. This heated gas flow rotates the turbine. Kinematic and heat structures of this gas flow at the studied period of time undergo the mutual effect of destabilizing factors, such as non-isothermal effect of the longitudinal pressure gradient caused by the uprising of the dynamic and heat transiency. Friction coefficient Cf and coefficient of heat losses St are the most important criteria of the gas flow, which characterize gas dynamics and heat exchange of the gas flow. Non-stationarity has an essential effect on the values of friction coefficient and heat loss coefficient, which increase or decrease in 2.0-2.5 times, compared to their quasisteady analogues. Beyond non-stationary process, when the temperature of working body becomes stationary, the key role is played by the thermal head which causes laminarization phenomenon of the of the heat turbulent boundary layer.
Keywords:
gas flow, power plant, friction coefficient, coefficient of heat losses, non-stationarity, non-stationary effects, non-isothermality
gas flow, power plant, friction coefficient, coefficient of heat losses, non-stationarity, non-stationary effects, non-isothermality
При эксплуатации газотурбинных двигателей выделяют четыре основных режима. Пусковому режиму отводится роль одного из наиболее важных и ответственных, поскольку за короткий промежуток времени в результате горения топлива происходит резкое повышение температуры рабочего тела и, как следствие, такие же резкие по величине изменения всех параметров двигателя, приводящие к повышению температуры обтекаемых поверхностей конструктивных элементов двигателя. В пусковом режиме двигателя возникают ситуации, когда происходит преждевременный отказ работы двигателя с последующим разрушением конструктивных элементов. Рассмотрим совместное влияние неизотермичности, эффектов нестационарности и продольного градиента давления, проявляющегося формированием турбулентного пограничного слоя на начальном участке цилиндрического канала, в условиях различной по величине динамики возрастания температуры газового потока. Экспериментальное исследование Эксперименты были выполнены на экспериментальном стенде, который представляет собой газодинамическую трубу разомкнутого типа (рис. 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. Povh I. L. Aerodinamicheskiy eksperiment v mashinostroenii. L.: Mashinostroenie. Leningrad. otd-nie, 1974. 478 s.
2. Repik E. U., Kuzenkov V. K. Issledovanie novogo metoda opytnogo opredeleniya poverhnostnogo treniya v turbulentnom pogranichnom sloe // Inzhenerno-fizicheskiy zhurnal. 1980. T. 38. № 2. S. 197-200.
3. Volodin Yu. G., Marfina O. P., Bogdanov A. N., Cvetkovich M. S., Kuznecov A. B. Izmerenie kasatel'nyh napryazheniy treniya v nestacionarnom gazovom potoke // Datchiki i sistemy. 2009. № 2 (117). S. 34-36.
4. Nikiforov A. N., Fafurin A. V., Fesenko S. S., Husnutdinov Sh. N. Issledovanie dinamicheskih harakteristik pnevmometricheskih priemnikov // Truda metrologicheskih institutov SSSR. M.: Izd-vo standartov, 1977. Vyp. 182 (242). S. 84-88.
5. Volodin Yu. G., Gil'fanov K. H., Marfina O. P., Zakirov I. F., Kazakov A. A., Kuznecov A. B., Ryzhakova Zh. S. Eksperimental'noe issledovanie teplovoy inercionnosti mikrotermopar // Pribory. 2008. № 4. S. 52-55.
6. Volodin Yu. G., Marfina O. P. Raschet koefficientov treniya i teplootdachi pri nestacionarnom neizotermicheskom techenii neszhimaemogo gaza v osesimmetrichnyh kanalah // Izv. vyssh. ucheb. zavedeniy. Mashinostroenie. 2007. № 3. S. 21-26.
7. Volodin Yu. G., Marfina O. P. Matematicheskoe modelirovanie puskovyh rezhimov energeticheskih ustanovok. SPb.: Info-da, 2007. 128 s.
8. Volodin Yu. G., Fedorov K. S., Yakovlev M. V. Nestacionarnye effekty i trenie pri zapuske energeticheskih ustanovok // Izv. vyssh. ucheb. zavedeniy. Aviacionnaya tehnika. 2006. № 1. S. 34-36.
9. Volodin Yu. G., Fedorov K. S., Yakovlev M. V. Nestacionarnye effekty i teploobmen v puskovom rezhime energeticheskih ustanovok // Izv. vyssh. ucheb. zavedeniy. Aviacionnaya tehnika. 2006. № 4. S. 41-44.
10. Volodin Yu. G., Marfina O. P. Granicy primeneniya matematicheskoy modeli nestacionarnogo techeniya neszhimaemogo gaza v osesimmetrichnyh kanalah // Vestn. Tehnolog. un-ta. 2016. T. 19. № 6. S. 130-131.
11. Kutateladze S. S., Leont'ev A. I. Teplomassoobmen i trenie v turbulentnom pogranichnom sloe. M.: Energiya, 1972. 320 s.
12. Volodin Yu. G., Kirpichnikov A. V. Laminarizaciya turbulentnogo pogranichnogo sloya pri puske energeticheskoy ustanovki i parametry laminarizacii // Vestn. Tehnolog. un-ta. 2016. T. 19. № 17. S. 122-126.
