Прямолинейное равномерное течение газов с теплоотдачей в энергетических установках летательных аппаратов Текст научной статьи по специальности «Физика»

вычайных ситуаций : материалы науч.-прак. семинара. Железногорск, 2011. С. 107-113.

4. Система дистанционного мониторинга пожаров и земель агропромышленного комплекса / Н. А. Тес-тоедов, В. И. Лавров, В. В. Двирный и др. // Мониторинг, прогнозирование и моделирование опасных природных явлений и чрезвычайных ситуаций : материалы науч.-прак. семинара. Железногорск, 2011. С. 122-129.

5. Проектирование надежных систем связи / В. Г. Афанасьев, В. И. Верхотуров, В. А. Заславский и др. ; под ред. акад. М. Ф. Решетнева ; МГП «Раско». Томск, 1993.

6. Голенков В. К., Гаврилов А. И. Создание и использование информационных ресурсов на предприятии // НТИ. 1990. № 12. С. 22-23. Сер. 1.

7. Королева Н. А., Тютюнник В. М. Методика оценки уровня обеспечения информационной безопасности организации // НТИ. 2007. № 1. С. 15-17. Сер. 2.

8. Бессонов А. Б. Организационно-правовое обеспечение информационной безопасности организации // НТИ. 2008. № 9. С. 16-28. Сер. 1.

9. Совершенствование системы информационного обеспечения процесса создания космических аппаратов связи / Р. П. Туркенич, А. А. Носенков, В. В. Двирный, С. Е. Терпигин // Решетневские чтения : материалы 15 Междунар. науч. конф. Красноярск, 2011. Ч. 1. С. 101-102.

10. Носенков А. А. Техническая совместимость: практика, наука, проблемы : монография / Сиб. гос. аэрокосмич. ун-т. Красноярск, 2005.

N. A. Testoedov, V. V. Dvirnyj, A. A. Nosenkov, M. V. Elfimova ABOUT DOMESTIC SYSTEM OF SPACE MONITORING OF EMERGENCIES

The problem of development of domestic system of space monitoring of emergencies with a focus on its core issues: relevance, recourses requirements, realization potential.

Keywords: space monitoring, spacecraft, emergency, design, production.

© TecToegoB H. A., ,3bhphlih B. B., HoceHKOB A. A., En^HMOBa M. B., 2012

УДК 532.526.4

М. И. Толстопятов, А. А. Зуев, А. А. Кишкин, Д. А. Жуйков, В. П. Назаров

ПРЯМОЛИНЕЙНОЕ РАВНОМЕРНОЕ ТЕЧЕНИЕ ГАЗОВ С ТЕПЛООТДАЧЕЙ В ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ*

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

Ключевые слова: турбулентное течение, пограничный слой, коэффициент теплоотдачи, теплообмен.

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

установки, который для получения максимального КПД имеет довольно узкий диапазон величин.

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

*Исследование выполнено при поддержке Министерства образования и науки Российской Федерации, соглашение 14.132.21.1584.

Рис. 1. Распределение температурного и динамического пограничных слоев при турбулентном течении газового потока

Для достоверного определения распределения температурных полей в проточных и вспомогательных полостях течения рабочего тела в тепловом расчете двигателя необходимо использовать локальный коэффициент теплоотдачи. В отличие от среднеинтегрального коэффициента теплоотдачи, который дает усредненные данные по длине исследуемого участка, локальный коэффициент теплоотдачи позволяет выявить возможные местные зоны перегрева рабочего тела. Учет локальных тепловых нагрузок очень важен для таких полостей как полость между ротором и статором турбины турбонасосного агрегата (ТНА), подводящее и отводящее устройство газовой турбины ТНА.

При проведении аналитического исследования прямолинейного турбулентного течения газа или газовой смеси необходимо рассмотреть модель распределения температурного и динамического пограничных слоев при обтекании газовым потоком плоской, прямой, не проницаемой пластины (рис. 1). Толщина температурного пограничного слоя выше толщины динамического. Такое распределение толщин в некоторой степени характеризуется теплофизическими свойствами рабочего тела, а именно: у газов и газовых смесей динамическая вязкость не может быть больше температуропроводности, а их отношение всегда будет меньше единицы. В свою очередь, основное изменение скорости (и) от поверхности обтекаемого тела, где она равна нолю, до скорости в ядре потока (Ц), происходит в толщине динамического пограничного слоя (5). Аналогично, основное изменение температуры (Т) от температуры на поверхности теплообмена (Т0) до температуры в ядре потока (Т5), происходит в толщине температурного пограничного слоя (5,). Тогда выражение для определения толщины потери энергии примет следующий вид:

8, / ^ т \

ёу . (1)

** г и

8^ = I и

1 _ Т - Т0

Т8- Т у

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

1 _-

ёу +

ёу.

(2)

Первый член выражения (2) отмечает изменение температуры и скорости в пределах от поверхности тела до толщины динамического пограничного слоя, второй член также, фиксирует изменение температуры и скорости в пределах от толщины динамического до толщины температурного пограничных слоев. Очевидно, что во втором члене выражения (2) изменения скорости нет, т. е.:

и = и или — = 1.

и

(3)

Тогда выражение для определения толщины потери энергии (2), в соответствие с выражением (3), запишется в виде

8

,;

=1

0

8,

Л

1—

ёу +

1|1

ёу.

Обозначим второй член выражения (4) через 5/ и проинтегрируем его в границах температурного пограничного слоя:

=і і'

А

Т - Т

Т8- То у

ёу.

(5)

где 8, - расстояние, на котором температурное ядро течения оттесняется наружу вследствие падения температуры в температурном пограничном слое.

Для аппроксимирования изменения профиля скорости [1], применяем «закон 1/7» степени. По аналогии с [2], используем отношение динамического и температурного пограничных слоев и обозначаем его

через

А:

А_^.

5

(6)

Тогда выражение для определения толщины потери энергии преобразуется к виду

8/ \1/7 ( / \1/7 Д

5*:_м

5/ (

1 -

1 -

У

А- 5

У

Д-5

1/7 Л

ёу +

(7)

ёу.

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

7

**

5/ф _"

5- і 9-А + -

1/7

-9

72

(8)

Выразим закон теплоотдачи чрез локальный коэффициент в виде критерия Стантона:

Б/ _

%о

у -

у=0

р-Ср- и- (Т5-То ) р- Ср- и - (Тр-То )

(9)

X

р - Ср - и

ду

Т - То V Т5-То

У=о

д_

дУ

[ Т - То ^ _ и [V -V'!

V Т5 - То у у_о а2 -V и3/4 15/'и )

(1о)

3/2 Л1/4 с-1/4 3/4

а, -А -5 -V

где аг - коэффициент ламинарного подслоя, полученный в работе [3]; здесь

а1 = 12,559. (11)

Выразив толщину динамического пограничного слоя из выражения (8) через толщину потери энергии и подставив в выражение (10), получим:

_д_

ду

[ Т - То Л

V Т5 - То У у=о

о,343 - и

3/4

(

\ 1/4

(12)

А1/4 -а,3/2 Vм

/ф

9-А +-

,1/7

-9

С учетом выражения (12), перепишем закон теплоотдачи в виде критерия Стантона:

Б/ _-

р- Ср- и

3/4

о, 343 -и

(

\1/4

А1/4 -а,3/2 Vм

/ф

9-А+~трт-9 V А1/7 .

(13)

Подставив выражение (13) в интегральное соотношение уравнения энергии, записанное в естественной криволинейной системе координат [4], получим интегральное соотношение уравнения энергии пространственного пограничного слоя для случая течения газа или газовой смеси, когда критерий Прандтля меньше единицы:

и

Н ф

д_

дф

дН,

дф

**

-5/ф +'

д(5*; )

д;

1

дН,

Решив выражение (9), видим, что необходимо взять производную на поверхности теплообмена от функции распределения температуры в температурном пограничном слое. Для того что бы определить производную, необходимо воспользоваться другой моделью распределения параметров в пограничном слое. По аналогии с [4] используем двухслойную модель с ламинарным подслоем, для которой производная для выражения (9) запишется в виде выражения

о,343 - X- V 9-А+дту-9

7: Л1/4 гг1/4 3/2 3/4 я**1/4

Ср - р- А - и - а, - V - 5*

д;

1/4

ф Я** — _-5/; _

(14)

/ф

фо

-(■+-2)

р-Ср - (Т5-То)

Принимая относительную характерную толщину:

**

Г _ 1-5;

**

є 5/ф

(15)

где е - тангенс угла скоса донных линий тока по от**

ношению к потенциальному течению, 8?ф - толщина

**

потери энергии в продольном направлении, 8^ -

толщина потери энергии в поперечном направлении. Переписав выражение (14), получим

X

1 д ,, з д(*;-е)

- 5/ф+-

Нф дф Н; д;

;

1 дН

; 0**

;- 5/ф +-

Нф -Н; дф ^ Нф -Н; д,

ф 0**

- 5/ф - е_

ф ;

о,343- X- і 9-А+-Г-Т-9

л 1//

1/4

(16)

п . 1/4 гл/4 3/2 3/4 я**1/4

Ср - р- А - и - аі - V - 5/(

/ф

фо

-(1)

Б/ _

о,956- Х4/3 - 9-А + Дрг-9

^ 4/5 т т-1/5 4/5 675 7Т75 3/5 1/5.

Ср -и -р - а, - А - V -ф

(2о)

Чтобы выражение привести к более простому виду, умножим числитель и знаменатель выражения (20)

на динамическую вязкость (ц4/5), а также воспользуемся безразмерным критерием подобия Рейнольдса:

„ р • и • ф) („ ц Ср

Яе = ------ | и критерием Прандтля: I Рг =------

Б/ _ о,4о662 -

[ 9А877 - 9А1/7 + 7 Л ч Рг4-Яе - Д8/7а,6 у

1/5

(21)

р-Ср• (Т8-Т0 у

При прямолинейном течении, поперечное течение д 0

отсутствует, т. е. — = 0 , а так же для случая прямо-

д—

линейного течения коэффициенты Ламе определяют-

дН— дН ф

ся Н— = 1, Нф = 1, тогда —— = 0, —ф = 0. Для

дф д—

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

1/4

д Ч 9-Д+Тш-9 '

—8** =------------------------------- (17)

дф ‘ф Ср • р- Д1/4 • и1/4 • аг3/2 • у3/4 • 8**1/4 " ^ ^

В интегральном соотношении уравнения энергии для прямолинейного равномерного течения (17) разделим переменные и проинтегрируем от нуля до текущего значения переменных, и определим толщину потери энергии:

8*ф 0,343 • Ь (9 • - + -^-9! ф

0 ^ё= Ср • р^ Д1/4 • и1/4 -аг3/2 • У3/4 0ёф. (18)

После интегрирования получим:

1,195• Х4/5 • (9• Д + -1у-9'] • ф4/5

8** =-------------------Д------------ (19)

*ф п 4/5 4/5 гг1/5 л 1/5 5/6 3/5 ' к ’

Ср • р • и • Д • аг • V

Окончательное выражение для определения толщины потери энергии необходимо подставить в выражение для вычисления локального коэффициента теплоотдачи (13):

1/5

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

А _ ^7=. (22)

^Рг

Выражение (21) с учетом формулы (22) примет

вид

Б/ _ о,4о662

[ 7-Рго,258 - 9-Рго,333 + 9 Л

Рг4-Яе - а,6

1/5

(23)

Для подтверждения методики получения локального коэффициента теплоотдачи в виде критерия Стантона воспользуемся аналогией Рейнольдса о процессах переноса тепла и импульсов, применяя выражение [5]:

Б/ _1 С/~. 2

(24)

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

На графике видно, что только при Рг = 1 сохраняется равенство 81 =1 С/ . В диапазоне Рг = 0,7...1, что 2

характерно для реальных газов и газовых смесей, среднее расхождение с теориями других авторов составляет не более 10 %. Стоит отметить, что сравнивать теории в области малых чисел Прандтля невозможно, так как зависимости, полученные различными авторами, предназначены для рабочих тел с числом Прандтля не намного отличным от единицы.

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

Библиографические ссылки

1. Шлихтинг Г. Теория пограничного слоя : пер. с нем. М. : Наука, 1969.

2. Кейс В. М. Конвективный тепло- и массообмен : пер. с англ. М. : Энергия, 1972.

тогда

+

3. Зуев А. А., Вахитов Р. Р., Мелкозеров М. Г. Уравнение теплоотдачи для турбулентного температурного пограничного слоя на основе двухслойной модели // Проблемы и перспективы развития авиации, наземного транспорта и энергетики «АНТЭ-о7» : материалы Междунар. науч.-технич. конф. (11-12 дек. 2оо7, г. Казань). Казань : Изд-во Казан. гос. тех. ун-та. 2оо7. Т. 2.

4. Интегральное соотношение уравнения энергии температурного пространственного пограничного слоя / А. А. Зуев [и др.] // Вестник Рыбинской гос. авиационной технол. академии имени П. А. Соловьева. 2о1о. № 2 (17). С. 37-42.

5. Романенко П. Н. Теплообмен и трение при градиентном течении жидкостей. М. : Энергия, 1964.

1.2 2St/Cf 10

".

N х

' *ч 1 к ч ч >1 1^ TeOpID

О " ", г „ ^1 h,:,— " " -і *- і к "■ 1 -ь У - —' 'і

ОД

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1 Pi-

—«— Полученная теоретическая зависимость —

■ - М.Рубезин —'

....Т.Карман, Г.Рейхард

■ - А.Кольбурн

- Л.Прандтль-Л.Тейлор - М.Широнов

Рис. 2. График сравнительного анализа аналогии между трением и теплообменом в турбулентных потоках при м » о, Яе=ю7

M. I. Tolstopyatov, A. A. Zuev, A. A. Kishkin, D. A. Zhuikov,V. P. Nazarov

RECTILINEAR UNIFORM GAS FLOW WITH THERMAL EFFICIENCY IN PROPULSION

EQUIPMENT OF AIRCRAFTS

The article considers a model of distribution of thermal and dynamic boundary layers distribution in the process of turbulent flow of the gas stream. Expression for determination of energy thickness is presented and expression for the determination of local thermal efficiency coefficient in the form of Stanton criterion. Comparative analysis of the theories of different authors is presented as well.

Keywords: turbulent flow, boundary layer, thermal efficiency coefficient, heat exchange.

© Толстопятов М. И., Зуев А. А., Кишкин А. А., Жуйков Д. А., Назаров В. П., 2012

13В