ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 132
Скачиваний: 0
ния тех или иных параметров течения на так называемые константы турбулентного обмена. Так, в работе [251 ] на основе интегрирова ния уравнения энергии для трехслойной схемы пограничного слоя получено выражение для фактора аналогии Рейнольдса
5 = |
Ч |
= Рг, + 5Рг, У Ц - {JL (1 - Рг,) [ £ + -§- <1 - РГ,)] + |
|
2 St |
|
|
|
+ ( |
w |
- |
1) + |
, n [ 1 + |
i ( |
^ |
- |
1) ] |
||
По данным работы [237] S |
Р х —0,57- по экспериментальным дан- |
|
||||||||||
ным S |
0,9 ч-1,1. |
В ра |
|
|
|
|
|
|
||||
боте [74] для этого пара |
|
|
|
|
|
|
||||||
метра |
получены значения, |
|
|
|
|
|
|
|||||
представленные на рис. 42, |
|
|
|
|
|
|
||||||
из |
которого видно, |
что |
S |
|
|
|
|
|
|
|||
зависит от М иф = |
TwITaw |
|
|
|
|
|
|
|||||
и значительно |
отличается |
|
|
|
|
|
|
|||||
от единицы при обтекании |
|
1,0 |
?,о |
з, |
|
|||||||
сжимаемым потоком. |
|
|
|
|||||||||
|
В работе |
[239] |
на ос |
|
Рис. 42. Зависимость S от ф |
|
|
|||||
нове обобщения |
результа |
|
|
|
||||||||
тов |
прецизионных |
опыт |
|
|
|
|
|
|
||||
ных данных предлагаются следующие формулы для расчета моди |
|
|||||||||||
фицированного фактора аналогии Рейнольдса а, представляющего |
|
|||||||||||
собой |
отношение |
коэффициента |
турбулентного |
обмена |
|
тепла |
|
|||||
к коэффициенту турбулентной вязкости а = К/А: |
|
|
|
|||||||||
|
1) |
при у+ < |
„+ |
|
|
|
|
|
|
|
|
|
|
г/? |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
а = Р г—0,5 |
|
|
|
|
|
|
2) |
при y t |
< |
у+ < y t |
|
п(у+ |
|
|
|
|
||
|
|
а- |
|
Рг°’5+ |
|
|
' У^т) |
|
|
|
||
|
|
|
1 -h (Рг0-5 — 1) sin |
] ‘ |
|
|
||||||
|
|
2 Рг0,5 |
|
|
|
|
|
|
||||
|
3) |
при у,+ > Ш |
|
|
|
|
|
|
|
|
а1.
При этом
У? = 4; УТ — Ю; y L yt + yt ytm — yt — yt
В работе [228] предлагается уточнение формулы Прандтля для длины пути перемешивания I путем введения демпфирующего множителя Ван-Дриста:
1) при у!б < 0,1
■f = 0,41 [ l - e x p ( - - g - ) ] - g - ;
9 |
131 |
Рис. 43. Зависимость 6л/6 от 1|) и М<ю |
Рис. 44. Зависимость |
т)0 от г|э при |
следующих |
значениях Мга: |
при R = 105 |
О — 6,8; • |
— 5,0; V — 3,0; |
X — 2,0; J |
— 1.5 |
|
2) |
|
при |
0,2 < |
у/8 < 0,6 |
|
- |
4 -= 0,41 1 — еХр ( — |^ ) ] |
— 1,53506 |
|||
|
|
|
+ 2,75625 ( - |---- 0,1 ) 3 |
1,88425 — |
|
3) |
при |
ylb > |
0,6 |
0,089. |
|
|
|
|
|
//8 = |
Отмеченное выше значительное влияние температурного фак тора на сопротивление и теплообмен пластины может быть объяс нено тем, что при изменении ф происходит существенная пере стройка пограничного слоя, которая характеризуется, в частности, изменением протяженности области вязкого подслоя. Из приве
денных на рис. 43 кривых [74] видно, |
что при ф = |
T J T aW — 2 и |
||||
Мсо *=» 5 толщина подслоя |
6Л может |
составлять |
50% |
толщины |
||
всего |
пограничного слоя |
8. Безразмерная |
толщина |
подслоя |
||
г|о = |
yvjv, являющаяся при обтекании пластины |
несжимаемым |
||||
потоком величиной постоянной и равной ri0 = |
10,65 (штриховая |
линия на рис. 44), с изменением температурного фактора в слу чае сжимаемого потока, как видно из рис. 44 [74], заметно изме няется. Линия на рисунке соответствует Мга г» 0.
Все эти данные говорят, что основные допущения, принимае мые при построении приближенных методов расчета турбулент ного пограничного слоя (Pr, = 1; т)0 = const; S = 2St/cf = 1; l = ху и т. д.), становятся практически неприемлемыми для рас чета теплообмена при значении ф, значительно отличающемся от единицы, и при больших значениях числа Мт .
21. Расчет турбулентного пограничного слоя при обтекании с продольным градиентом давления
Задача о градиентном обтекании поверхности потоком сжимае мого газа (обтекание крылового профиля и криволинейных поверх ностей, течение в межпрофильных каналах решеток турбинных и компрессорных лопаток и т. п.) значительно осложняется тем, что при сверхзвуковых течениях сильные продольные градиенты давления или изменение кривизны обтекаемой поверхности могут сопровождаться возникновением скачков уплотнения и вызывае мых ими отрывов пограничного слоя. Неоднородность основного потока, непостоянство давления поперек пограничного слоя и другие факторы, обусловленные взаимодействием скачка уплот нения с пограничным слоем, приводят к невозможности примене ния уравнений пограничного слоя в обычном виде. Описанные в литературе методы расчета как ламинарного, так и турбулент ного пограничного слоя у поверхностей, обтекаемых сжимаемым газом с продольным градиентом давления, основаны на допущении
133
об идеализированном (бесскачковом и безотрывном) обтекании. Это нужно иметь в виду, используя эти методы для решения прак тических задач.
Опытный материал по турбулентному обтеканию поверхно стей с продольным градиентом давления еще более ограничен, чем для случая обтекания плоской пластины.
В настоящее время появился ряд опытных данных, указываю щих на то, что при некоторых условиях (очень больших ускоре ниях потока, очень большом повышении температуры потока и т. п.) возможно возникновение явления обратного перехода от турбулентного режима течения в пограничном слое к ламинарному. Имеются опытные данные, наглядно иллюстрирующие обратную перестройку профилей скорости в пограничном слое и соответ ствующие изменения коэффициентов сопротивления. Не вдаваясь в подробный физический анализ этого явления, который сейчас еще недостаточно ясен, можно дать следующую качественную оценку этого явления.
Напишем интегральное соотношение импульсов в виде
т |
Дб** + 6** (2 + Я )U^' + |
d p ' |
РV* |
d x |
d x _ |
Пусть профиль скорости в пограничном слое соответствует закону одной седьмой. Тогда 8*/8** 1,28. Соответственно допустим, что для т можно использовать выражение [118]
—0,25
т = 0,0128pt/2 ( - ^ —)
Введя эти значения в интегральное соотношение импульсов после преобразований получим
d R * *
d x
у |
_ |
0,0128_______ у_ |
2,29UR** |
|
2,29 (R**)1,25 |
<Ш
(IV.98)
d x '
Если принять, что обратный переход происходит при тех же значениях R**, что и обычный переход от ламинарного к турбу
лентному течению, т. |
е. при R** |
360, |
то |
|
|
d R** |
у |
|
|
v d U |
(IV.99) |
d x |
825U — 3,55-10~6— |
Тр Чх ' |
|||
Поскольку при обратном переходе |
|
величина |
dR** |
||
|
d должна |
уменьшаться, т. е. dR** < 0, то условием обратного перехода
будет
-jjs- > 3,5 • 10~®.
Существующие в настоящее время полуэмпирические методы
расчета теплового и динамического пограничных слоев при гра
134
диентном обтекании можно разбить на две группы. К первой группе относятся методы, базирующиеся на решении дифференциальных уравнений в частных производных [уравнения (VII)—(X) в гл. I] с привлечением эмпирических соотношений для зависимостей так называемых эффективных коэффициентов переноса (вязкость, теплопроводность, число Прандтля и т. д.). Ко второй группе относятся методы, часто называемые эмпирическими, основанные на решении обыкновенных дифференциальных уравнений (обычно интегральные соотношения импульсов и энергии).
В настоящей монографии рассмотрены только те методы рас чета, которые могут быть применены или уже применяются для расчета теплового и динамического пограничного слоя в элемен тах проточной части турбомашин.
Из работ первой группы в этом плане особого внимания заслу живает метод Патанкора и Сполдинга [226]. Эта работа является дальнейшим развитием идей метода Сполдинга и Чи [241 ], описан ного в предыдущем параграфе. Широкое применение ЭВМ позво лило авторам произвести подробный численный эксперимент и сопоставить результаты расчетов с экспериментальными данными с целью отбора гипотез, определяющих коэффициенты турбулент ного обмена, использование которых дает наилучшее совпадение расчетных результатов с непосредственными измерениями для конкретных задач. На основе этих сопоставлений проведено дальнейшее развитие метода и усовершенствование основных физических гипотез, а также усовершенствование математического аппарата с точки зрения повышения его точности и экономичности. Метод численного решения основной системы дифференциальных уравнений пограничного слоя является конечно-разностным, он подразумевает использование ломаных профилей из прямолиней ных отрезков. Независимой переменной в поперечном сечении пограничного слоя служит безразмерная функция тока. Непре рывное изменение эффективной вязкости поперек погранич ного слоя характеризуется модифицированной гипотезой ВанДриста [252]:
\Ьи = И + Р*У 1 — ехр - У |
ди I |
(IV. 100) |
|
д у г |
|||
|
|
где величина
ди
P*V ду
называется турбулентным аналогом вязкости.
Метод решения задачи запрограммирован на ЭВМ, программа записана на алгоритмическом языке «Фортран-IV», подробное ее
изложение приведено в книге |
[226], недавно вышедшей на рус |
ском языке. |
. . . . . . . |
Эмпирические методы решения (методы второй группы) в основ ном строятся на решении интегральных соотношений импульсов
и энергии. Наибольшее распространение из них получили так называемые однопараметрические методы расчета в силу своей простоты и удовлетворительной точности при решении сравни тельно простых практических задач, возникающих в авиации, турбостроении и других отраслях техники. Сущность этих мето дов заключается, как уже говорилось выше, в предположении, что все многообразие профилей скорости в пограничном слое обтекаемого тела можно представить в виде однозначной функции некоторого форм-параметра f, изменяющегося вдоль обтекаемой поверхности. Здесь также сначала были созданы методы расчета пограничного слоя при обтекании поверхности потоком несжимае мого газа, которые впоследствие тем или иным путем распростра нялись на сжимаемый поток газа.
Так, известный метод Солодкина для расчета обтекания сжи маемым потоком представляет собой обобщение известного одно параметрического метода расчета турбулентного пограничного слоя в несжимаемом газе Лойцянского [115]. Распространение идей этого метода расчета на тепловой пограничный слой шло двумя путями [60, 150].
В работе [60] разработан приближенный полуэмпирический метод, позволяющий рассчитывать единообразно все три области теплового пограничного слоя (ламинарного, переходного и тур булентного режимов течения), возникающие при обтекании по верхности сжимаемым газом. Рассматривается обтекание поверх ности плоским потоком вязкой сжимаемой жидкости при наличии теплообмена в предположении постоянства температуры поверх
ности. Интегральное соотношение энергии |
|
|||
|
■+ <£ |
и |
То Nil* |
(IV.101) |
d\ |
_L |
|||
|
и |
Тц) Р г R* |
|
решается в переменных Дородницына. Здесь переменные | и т)
определяются по (1.53), 6Т — по (1.77), t и to — по (1.39); Nu* = ах/Х\ Pr = v*la*.
Использовав соотношения (1.54), (1.56), (1.61), можно уравне
ние (IV. 101) привести к |
виду |
ft |
|
|
|
I и ' я** |
£I ш 0 |
Nu, |
|
||
2\ f t - I |
(IV.101a) |
||||
+ 7Г6т |
|
■оЗ) |
Pr Rjc |
Вводятся параметры:
сV
/т ~Ц~ От
.
(IV. 102)
и предполагается, что эти параметры, изменяющиеся вдоль обте каемой поверхности, однозначно определяют все характеристики теплового пограничного слоя. Принимается допущение, что в пер
136