Файл: Зысина-Моложен, Л. М. Теплообмен в турбомашинах.pdf

где звездочка соответствует параметрам торможения, позволило продвинуть задачу о расчете турбулентного обтекания в сжимае­ мом потоке. На базе использования переменных Дородницына был разработан ряд полуэмпирических методов расчета теплооб­ мена и сопротивления в турбулентном потоке [84, 150, 60 и др. ], позволивших приближенно решать большой круг практических задач.

Об ограниченности этого пути решения говорилось в пп. 5 и 12. Однако для задач турбостроения, в которых даже для пер­ спективных машин М 2,5 ч-3,0, это ограничение не является существенным и использование переменных Дородницына пло­ дотворно.

3. Третий путь заключается в попытках решения основной полной системы дифференциальных уравнений при условии за­ мыкания ее эмпирическими соотношениями, построенными на той или иной формулировке гипотез Прандтля или Кармана с исполь­ зованием так называемых констант турбулентного обмена. Эти константы получаются из экспериментов для обтекания несжимае­ мой жидкостью без теплообмена или при значении температур­ ного фактора ф = T J T aw 1. Подавляющее большинство иссле­ дований по влиянию сжимаемости и температурного фактора на турбулентное обтекание проводится применительно к безградиентному обтеканию наиболее простых объектов — пластин и конусов. Различные полуэмпирические решения этой задачи строятся на основе использования гех или иных гипотез для опре­ деления констант турбулентного обмена и их связи друг с другом или для определения турбулентного трения и его связи с усреднен­ ными скоростями.

Как уже говорилось выше, кроме проблемы внутренних за­ кономерностей в турбулентной части пограничного слоя возни­

в нем и влияния на них числа М и температурного фактора ф. Решению этих двух вопросов для турбулентного течения и посвящено большинство полуэмпирических и эмпирических иссле­

дований в этой области.

В известной работе Ван-Дриста для обтекания сжимаемым

124

где

du

~dxp

В формуле (IV.86) величина cpAIK = Рг, иногда называется турбулентным числом Прандтля.

На основе этих допущений получается известное соотношение между коэффициентами теплообмена и сопротивления трения

В результате получается так называемыйсинус-логарифми- ческий профиль скорости и соответствующие формулы для рас­ чета теплообмена и сопротивления.

Аналогичные допущения принимаются и в ряде других работ. В работе Сполдинга предлагается принять допущение о существо­ вании особой формы гипотезы Прандтля

где

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

Дейслер в ряде своих работ пытался найти подход к решению вопроса о структуре турбулентного пограничного слоя пластины. Он рассматривает двух- и трехслойную схемы пограничного слоя

ипринимает допущение, что в пристенной области при т] < 26

имеет место соотношение

125

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

Cflcf ис при R —>оо, при этом авторы полагают, что для сжимаемого потока сохраняется соотношение

(IV.94)

На рис. 36 приведено сопоставление шести различных мето­ дов расчета по влиянию сжимаемости на Cf для теплоизолирован­ ной пластины. Как видно, имеет место практическое совпадение данных. Кривая Кутателадзе и Леонтьева лежит ниже других. Все эти теоретические решения в целом хорошо соответствуют опытным значениям коэффициентов сопротивления. Это позво­ ляет считать , что для расчета теплоизолированной пластины не имеет существенного значения, какая структура турбулентного пограничного слоя кладется в основу теоретического метода рас­ чета.

Совершенно другая картина получается при интенсивном теп­ лообмене поверхности, т. е. при значениях температурного фак­ тора TwITaw, значительно отличающихся от единицы. На рис. 37 приведено сопоставление расчетов по влиянию сжимаемости на сопротивление трения по методам [240, 241, 199]. Как видно,

126

расхождение между результатами расчетов с уменьшением темпе­ ратурного фактора увеличивается. Наилучшее соответствие с экспериментальными данными (нанесены точками), полученными

для 0,6 ф ^ 1,0; 2 < М с Ю ; 7- 10е==£ R*, 2 • 109 [244],

имеет метод Хеббора и Паранжпи [199], в котором влияние сжи­ маемости на Cf выражается простой формулой

даже два метода Ван-Дриста, один из которых основывается на гипотезе Прандтля для пути перемешивания, а другой — на ги­ потезе Кармана, дают значительно отличающиеся результаты. Это различие еще больше, когда теплообмен происходит в средах

счислами Рг, значительно отличающимися от единицы.

Внастоящее время наиболее популярными методами расчета профилей скорости и сопротивления при безградиентном обтека­

нии пластины потоком сжимаемого газа являются методы Спол­ динга и Чи [241] и метод Коулза [185].

В работе [241] на основании анализа результатов исследова­ ний двадцати различных авторов вводится модифицированная

127

форма гипотезы Прандтля в виде (IV.91) и производится числен­ ное интегрирование системы уравнений движения методом ко­ нечных разностей с помощью ЭВМ. На основе этих вычислений построены номограммы, позволяющие быстро выполнять расчет

(рис. 39, 40).

Рис. 38. Зависимость St/StHC от температурного фактора при 2,5 3,5 и 0,8-106 ^ ; R e g sg 1,4 10»:

1 — метод Эккерта; ■2, 3 — [251 ]; 4 — [241]; 5 — ме-

тод Винклера и Чи

С помощью рис. 39 по заданным ф = T J T и М определяются значения условных функций Fc и Fr :

1

находится локальное или среднее значение коэффициента сопро­ тивления трения.

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

128

Рис. 39. Номограмма для определения условных функций Fc и FR

5)

Рис. 40. Расчетные значения су по методу Сполдинга и Чи: а — локальные; б — средние

9 Л. М. Зысина-Моложен и др.

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

В работе Коулза [185] логарифмический профиль скорости в пограничном слое, полученный при обтекании пластины несжи-

Рис. 41. Номограмма для определения коэффициента сопротивления трения: а — R** = Ю3; б — R**' = 104;

-------------- адиабатное течение;--------------- неадиабатное течение

маемым потоком газа, уточняется для сжимаемого потока путем введения дополнительного члена

1

где П — параметр, зависящий от градиента давления; 6С— зна­ чение у, при котором

Зависимость (IV.96) получила название закона следа.

Формула для расчета коэффициента сопротивления получается в работе [185], как обычно для логарифмического профиля ско­ рости, в неявном виде. Однако недавно Майлз и Ким [220] сравнили экспериментальные данные 18 работ с расчетами по ме­ тоду [185] в диапазоне изменения параметров R** = 103-т-104; TwITm = 0 ,5 -т-8,0 и установили, что отклонения не превышают 7,3% для адиабатного течения и 15,3% для неадиабатного течения. На рис. 41 приведена номограмма, полученная в работе [220], позволяющая быстро определять изменения Cf с изменением чисел

130