Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Зысина-Моложен, Л. М. Теплообмен в турбомашинах

.pdf
Скачиваний:
10
Добавлен:
22.10.2023
Размер:
15.66 Mб
Скачать

Величину б* обычно называют толщиной вытеснения, а величину

6Т— толщиной

теплового

вытеснения.

 

При введении в (1.38) обозначений по (1.40) можно получить

6

 

 

 

 

 

J

(Pot/ — pu) dy = — Рои

8* + pot/6*.

(1.41)

о

 

 

 

 

 

Подставив соотношение (1.41) в уравнение (1.37), получим

зависимость

 

 

 

 

 

d

 

РОи

dU

 

(1.42)

dx [p o t/V * =

dx

 

или

 

 

 

 

 

Pot/2^

+

W * ^ +

2p 06**U

 

п dU

гр

0Х б’

= P°U dx-

( 0

<5

*

 

' 0

 

 

 

Обозначив

 

 

 

 

И? II

• f =

 

t / '

'

 

dx

 

 

— р( диду

( ди

){/=0

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

db**

+

б**

 

о к

-I-

U'

б* —

и

JJ* \

(1.43)

dx

 

р 1

и

и

т~0

т

Если, как это иногда делают,

ввести определение

t и to не

по (1.39),

 

а обозначить

 

 

 

 

 

 

 

 

 

 

 

t* = T ' - T w\

t'o =

r 0- T

w,

 

то уравнение (1.43) будет иметь более сложный вид:

 

db**

 

,

Ро

U’_

(26** +

б’ +

аоб***

 

_

Тxs)

dx

"т-

ро

U

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

U2

=

1_РоР_

 

 

U'и ) dy.

 

 

 

ао

* рТ'о’

 

 

 

Рассмотрим частный случай уравнения (1.43).

При обтекании поверхности потоком сжимаемого газа при практически одинаковых температурах стенки и потока, т. е. при

значении так называемого температурного фактора T J T

1,

31

физические свойства потока будут постоянными и р = р0; ро = 0. Тогда уравнение (1.43) приобретает вид:

^ + -^(26**+ 8*) = - ^ .

(1.44)

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

6**= оJ - H '- T - W «* = о/(>- т г ) “ У-

<145>

Как видно из сопоставления уравнений (1.43) и (1.44), в рас­ сматриваемом частном случае из уравнения интегрального соот­ ношения импульсов выпадает член, содержащий характеристики

теплового пограничного слоя 6Т, поэтому динамическая задача решается независимо от тепловой.

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

Интегральное соотношение энергии получается аналогичным

методом с помощью уравнений (IX) и (X).

 

 

 

Вводятся тождественные

преобразования:

 

 

 

д иТ*)

д vT*)

_

* Г д и)

,

д (ро) 1

 

дТ*

,

дТ* . , 4 g,

дх

ду

-

1 [

дх

+

ду

J

+ p U

дхдх

 

ду V м *

 

д(риТ'0)

а ( р О

_

дТ\

 

pv

 

 

 

 

дх

+

 

 

= р« дх

 

ду

+

 

 

 

 

т * Гд (ри) , д (ри) ]

 

 

(1-47)

 

 

 

То [ - g j - + —

 

J •

 

 

Соотношение (1.46) подставляется в уравнение (X), и получен­ ный результат вычитается из уравнения, получаемого при под­ становке (IX) в тождество (1.47).

Далее, используя все принятые ранее обозначения, а также вводя понятие толщины потери теплосодержания

6Т**

6 Т

(1.48)

о

32

и интегрируя в пределах теплового пограничного слоя, можно получить уравнение

Л

ср

■pot/o [б -б * + (1

(1.49)

Отсюда получаем интегральное соотношение энергии в развернутом виде

**

dbТ dx

\.

1

+| з +

 

?

1

т**

б* +

- 2- Гб —

,

и2

db***[

|_g***

LT_

 

ь /о

L dx

 

U

 

 

 

(

i

 

Я

(1.50)

 

P0cpU/o

 

 

 

'

Как видно, в уравнение (1.50) наряду с условными толщинами

теплового пограничного слоя 6Ти 8Т входят также условные тол­ щины динамического пограничного слоя б*, б**, б***, т. е. тепло­ вая и динамическая задача должны решаться также совместно.

Частные случаи уравнения (1.50):

U2

а) при малых скоростях ----— 1, тогда

2сЛ

 

Ро

 

Uo

тг+ «

 

 

Я

(1.51)

 

WpUt'o

'

б) при T J T q

1 и при М 0

имеем То = const, тогда

 

 

(1.52)

где

5. Представление интегральных соотношений пограничного слоя в переменных Дородницына

Из сравнения формул (1.44) и (1.52) с формулами соответ­ ственно (1.43) и (1.50) видно, насколько проще решается задача о расчете теплообмена и сопротивления при обтекании тел потоком несжимаемой жидкости. Кроме простоты формы уравнений очень

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

33

важно, что для несжимаемого потока тепловая и динамическая за­ дачи полностью разделяются, т. е. в уравнении (1.44) мы имеем дело только с характеристиками динамического пограничного слоя, а в уравнении (1.52) — с характеристиками теплового по­ граничного слоя.

Вследствие этого для определения сопротивления не нужно рассчитывать тепловой слой, а для расчета теплоотдачи не нужно рассматривать слой динамический, как это имело место в общем случае обтекания сжимаемым газом. Это обстоятельство весьма упрощает расчеты.

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

Наиболее широкое распространение получили преобразования Дородницына—Стюартсона.

А. А. Дородницын предложил в 1942 г. [42] преобразование координат, которое учитывает влияние сжимаемости и в то же время позволяет придать уравнениям для обтекания сжимаемым потоком форму, близкую к уравнениям пограничного слоя для обтекания несжимаемым потоком. Было предложено вместо обыч­ ных переменных х и у ввести в основную систему уравнений новые переменные:

6 = (1.53)

Здесь ро и ро — значения р и р соответственно в адиабатно и изоэнтропно заторможенном внешнем потоке.

Введем новые переменные в уравнения (VII), (IX), (X), учи­ тывая, что

д _ д£ д

. дг± _3___р_ ___,

_д_ .

дх

дх

дх дт]

р*

дх дт) ’

д _

_р_ _д_

(1.54)

 

~ д у ~

Ро *1 '

 

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

ft

p = p o (l-« S ) *- 1 .

(1-55)

Здесь

U

а 0

^2Г0 '

34

Соответственно из (1.55), используя (1.54),

можно получить

fe+i

( l _ a g ) - P ^ - ^ - .

(1.56)

- ^ ( l - a g ) * - 1 p0t/ - g - =

Поскольку внутри пограничного

слоя

0, что

следует

из уравнения (VIII), то, очевидно, соотношение (1.55) будет спра­ ведливо в любой точке пограничного слоя.

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

Р

Т о

Ро

Ро

Т

(1.57)

(l-ag)*/**-11

(1.58)

И м

О

Л

 

t * \

-\п

(1.59)

‘ ■ -“ ’ - К т т - ' М ' + т

 

Ро

 

 

Здесь

 

 

 

 

 

. а,- ~

=

 

(1.60)

 

V 2«о

У 2срТ'0

 

 

Используя полученные соотношения, можно переменные Дород­ ницына представить в виде:

| = J ( i _ ag)*-1 dx-

h = ( i

g)A-l l

___________ ^ -----------------

(1.61)

— ®о)

 

Если ввести

эти переменные в уравнение (VII) и

поделить

все члены уравнения на комплекс рр/р0, то можно получить урав­ нение движения в виде

3'

36

Если ввести

обозначение

 

 

 

 

____и____дп ,

 

 

и,

(1.63)

 

k

дх

-I"

 

 

(! -

‘А*-1

1 -— а 2 +

 

 

 

 

«о)

 

 

 

 

 

 

то можно представить уравнение (1.62) в виде

 

 

 

ди

, ~

ди

 

 

U

dl ^

 

Ж

 

1— об

 

 

+

v0 ЙГ)

 

 

II

п—1

ди

 

1-

“* + № - 1) ( 1 - 1

$

 

ждЧ-

(164)

Аналогичным образом вводя переменные Дородницына в урав­ нение (IX), можно получить следующий вид уравнения нераз­ рывности:

 

 

 

 

 

 

ди

.

d v

~

 

 

 

 

(1.65)

 

 

 

 

 

 

Ж

+

aTf

и'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

 

выражение

 

 

 

 

 

 

 

 

dUu

 

.

dUv

 

 

 

 

 

 

dU

~

ди

 

 

д\

 

1

5ri

 

= y(t + ^)+“ а Г + и ^ Г

 

( 1.66)

Используя уравнение

(1.65)

и учитывая, что

dU

= 0 , равенство

(1.66) можно

свести

 

к

выражению

 

 

 

 

 

 

 

 

 

 

dUu .

dUv

ди =

0.

 

 

 

(1.67)

 

 

 

 

 

ag '

Зт)

 

д1

 

 

 

 

 

Аналогично

можно

получить

 

 

 

 

 

 

ди2

 

дт

2и

.

~

ди

dv

 

ди ,

~

ди

(1.68)

ag +

"ЗтГ

 

 

 

U-3— =

и -==- + V

ай

 

 

 

 

 

ац

 

ag

1

 

Если ввести соотношение (1.68) в уравнение (1.64), то

(1.69)

Вычитая из уравнения (1.67) уравнение (1.69) и производя преобразования, аналогичные преобразованиям в п. 3, можно получить интегральное соотношение импульсов в переменных Дородницына:

 

u't 5** I 2

2

 

 

 

 

db** .

0

6*

То

г w

(1.70)

" 5 Г +

и

1— «8

1— d)

Т w

p0U2

 

36

Здесь

6*‘ = f l g r ( ' —

(1.71)

о

о

Введя переменные Дородницына в уравнение (X), разделив все члены уравнения на рр / р 0 и использовав выражение (1.63), получим уравнение энергии в виде

х U -

1—

дТ* . ~ дТ*

и —S.— v -з— =

д% 1 дт)

г,кЦк-\)

0 -« 8 )

1д

*ап X

Рг(1-«8)

[ > + ( Р г - | ) * ] [ . (1.72)

Используя уравнение сплошности (1.65), можно получить соотношения:

диТ0

WTо

*(*+£) +

~W~. +

ап

 

дТп

~ дТп

 

dTn

(1.73)

•Ь и- д% + и-an

 

d \ '

 

диТ* , dvT* ю

( ди ,

до

\

,

дТ* ,

г д г + ~ * г = т 1 а г + ^ г ; + и ж +

, ~

дТ*

дТ*

. ~ дТ*

(1.74)

+ " - * г = “ - а г + и ^ г -

 

Подставим соотношение (1.74) в уравнение (1.72):

диТ*

dvT*

 

 

 

ап

X

 

ап

 

 

 

 

 

Pr (1 — ай)*

1

 

(1 —

11

л

г

 

 

и2 -|

х U--------- /Г

N/----- Г + (Р [-

' )

г

(1.75)

1~ “ ’ + ( ^ - , ) ( ' - ^ ) ^ 1

 

 

"J

Далее, вычитая почленно из уравнения (1.73) уравнение (1.75) и интегрируя результат по т] в пределах теплового пограничного слоя, получим интегральное соотношение энергии для изоэнтропного течения в виде

аб

+ £

Ul_ I _ Ч _ \ _

то

Я

(1.76)

4

+ t ' J ~

tw

pcput;

37

Здесь

б

в

Сравнение уравнений (1.70) и (1.76) с уравнениями (1.44) и (1.52) показывает, что введение переменных Дородницына позво­ ляет разделить переменные теплового и динамического слоя, и сами уравнения становятся аналогичными уравнениям для об­ текания несжимаемым газом.

Использование координат Дородницына позволяет производить расчеты теплоотдачи и сопротивления при сверхзвуковых ско­ ростях без необходимости дополнительного учета влияния сжи­ маемости при числах М ^ 2,5 [118].

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

| | Г Л А В А

ЭЛЕМЕНТЫ ТЕОРИИ ПОДОБИЯ

6. Основные представления

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

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

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

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

39

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

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

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

7. Условия подобия

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

40