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

книги из ГПНТБ / Борисенко А.И. Аэродинамика и теплопередача в электрических машинах

.pdf
Скачиваний:
65
Добавлен:
25.10.2023
Размер:
17.94 Mб
Скачать

получаем уравнение энергии в виде, характеризующем процессы внутри рассматриваемого элемента,

 

dT

dw

д*Т

(1-68)

 

РсV dt

Р Ж

ді2

 

d9 I Ddw

Так как в

силу уравнения неразрывности

 

dw

V dS

p dv ,

d t ' v d l ~

, I

= 0' T0 член

 

~ТШ’==Ж 7Г (где ^=1/р-удель-

кый объем) характеризует работу деформации элемента (здесь это работа только сил давления, так как силы вязкости не проявляются). Таким образом, тепло, подве­ денное к газу (Я д2Т/дР), идет на увеличение внутренней энергии рcydT/dt и на работу расширения р dw/dl объ­ ема элемента.

Если нет теплообмена (адиабатическое движение), нет сил трения и жидкость не совершает механической работы, то интеграл уравнения энергии для струйки (1-66) носит название уравнения Бернулли и имеет вид:

сѵТ +

+ —- + 0 2 = const.

(1-69)

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

1 Г + ^ + у = С0П8+'-

(1'70)

В этой записи все члены имеют размерность длины. Смысл этого уравнения: для всех точек идеальной жид­ кой струйки сумма высот — скоростной (wzl2g), пьезо­ метрической (p/pg) и геодезической (у), отмериваемых от одной горизонтальной плоскости — есть величина по­ стоянная.

Давление р, входящее в уравнение (1-69), в гидрав­ лике называют гидростатическим, а величину р + рш2/2

гидродинамическим давлением.

Во многих задачах (в частности, при рассмотрении вентиляции электрических машин) влияние изменения величины у ничтожно; замечая еще, что p = pRT и R =

сѵСѵ, можно написать:

срТ + - ^ - = c o n st= < ? p7'0.

(1-71)

40

Температура Т0 = Т -\-~^г носит название темпера­

туры заторможенного потока.

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

2

2

 

 

0W\

0Wo

 

 

Р . + - 2- = Р г + - 2^ +

ЛРі,2-

0 "7 2)

Подчеркнем, что уравнение Бернулли справедливо для

одной и той же струйки

(линии тока); константа в нем

имеет различное значение для каждой линии тока.

В общем случае трехмерного

движения

уравнение

энергии в прямоугольной

системе

координат,

связанной

с частицей, имеет вид:

 

 

 

(1-73)

В векторной форме

pcH 4 r s p * v ( 4 r + wv:r ) + />(Ѵ«0 = ЯѵіТ + pD, (1-74)

где

(xD = p,

 

d w x

 

dwn

dwx

d w t

(div w)a -j-

dy

+

d x УЧ

ö z +

dx +

 

dw

 

dwv

 

(1-75)

+( dy

+

dz }

 

 

 

— так называемая функция рассеивания (диссипации); это работа сил трения, переходящая в тепло, которое рас­ сеивается в массе жидкости.

Физический смысл этого уравнения: тепло, подведен­ ное к жидкости ХѴ2Т и выделяющееся вследствие работы сил вязкости pD, идет на увеличение внутренней энергии рс dTjdt и на работу расширения p(Vw).

41

Динамический пограничный слой. Отсутствие сколь­ жения на твердых границах («прилипание») приводит к большим градиентам скорости вблизи этих границ и зна­ чительных касательных напряжений т = ц dw/dn внутри тонкого пограничного слоя. В этом слое влияние сил вязкости соизмеримо с влиянием сил инерции, тогда как вне его силы вязкости пренебрежимо малы по сравнению с другими силами (инерции, давления, массовыми) и те­ чение близко к потенциальному.

Расстояние от стенки до места, где скорость течения «почти равна» скорости внешнего течения, называется толщиной пограничного слоя б. Эта величина является несколько неопределенной.

Рассмотрим

двухмерный установившийся поток

2= 0, d/dz=0;

dfdt = 0 ) вязкой несжимаемой жидкости

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

 

 

 

dwx

 

I

dwu

 

 

 

 

 

 

 

 

 

 

дх

 

*

ду

 

 

 

 

 

 

 

Wx

dwx

ШУ

dwx

__

 

1

др

1

г (

d2Wx

1

dswx

(1-76)

 

дх 1

ду

 

 

р

дх

 

у дх2

Т

ду2

wx

dwy

dwy

~~

 

1

др

1

V ( д*щ

1

d2Wy

 

 

дх

 

ду

 

Р

ду

 

1дх2 П

 

ду2

 

 

Пусть

ось

X направлена

вдоль

границы, а

ось у

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

x=lU У = цЬ,

где £ и г| — безразмерные величины.

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

S)=- j^

i f^f'-Л|.

0

О

42

Если скорость wx изменяется от значения нуль у стен­ ки до некоторого значения wox на внешней границе слоя, которое сравнимо с некоторой характерной для данного потока скоростью w (например, со скоростью невозму­ щенного потока), то производная dwx/dl мало зависит от толщины слоя и будет порядка w. Следовательно, попе­ речная скорость wy имеет внутри слоя величину порядка

W Ö / L

Тогда в левой части второго уравнения (1-76) полу­ чим: порядок величины первого члена wxdwxjdx, равный w2!t\ порядок величины второго члена wydwx/dy, также равный w2jl.

В правой части

этого

уравнения

получим,

принимая

в качестве

масштаба давления

скоростной напор

р&у2/2 ,

порядок

величины первого члена

~ ~ —j57-.

равный

w2[l;

 

 

 

 

 

 

 

 

Р

дх

 

 

д

dw-

порядок

величины

второго

члена ѵ

d2w,

=

12

 

 

дЬ

равный vw/F; порядок величины третьего члена

d2w~

ѵ - ~ =

д2

д

dwx , равный ѵда/82.

 

 

 

 

 

 

 

ârj

 

dr;

 

 

 

 

 

 

 

 

 

 

Рассматривая

порядок

величин всех членов

уравне­

ния по отношению

к порядку

величины ш2//, получаем,

что член

V•

d2wx

s

 

 

 

ѵш>

W2

 

V

 

1

- дх2

будет порядка

 

 

Т

 

wl

 

Re

а порядок величины члена ѵ-

d2w~

будет

 

 

W2

ду2

 

Т

___lj_ _ 1_

 

 

 

 

 

 

 

 

 

 

 

— ‘ ö 2

Re

 

 

 

 

 

 

 

 

 

 

 

 

Сравнивая между собой члены, характеризующие вли-

d2wx

яние вязкости, можно сказать, что ѵ- ду2 будет в /2/82

d2w

раз больше ѵ—-^а—. Опыт показывает, что при больших

Re отношение

Ь/1 мало. Поэтому при больших числах Ре

,,

 

d2wx

• •

можно пренебречь членом ѵ

 

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

Аналогичной оценкой порядка величины членов третье­ го уравнения получим:

dwy

w2S

wtl

dwy

w2S

dp_

~T

дх

~~г~’

~~ду~

~ ~ ~ ;

P dy

 

d2wy

 

vwS

d2Wy

\w

 

 

дх2

 

 

öy2

 

 

Так как все члены второго уравнения очень малы по сравнению с членом др/ду, то можно принять при боль­ ших Re

дрІду = О,

(1-78)

т. е. давление внутри пограничного слоя одинаково вдоль нормали к контуру тела и равно давлению на внешней границе слоя на рассматриваемой нормали.

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

с а н ы в ф о р м е

dwx

 

 

 

 

 

I

ÖWy — п-

 

)1

dwx I1

дх

1

ду

 

1

dwx

 

1 dp 1

d2wx

[ ( 1 - 7 9 )

а* Н

äy -

 

P дх

д у •

 

1

 

dp

- 0 .

 

1

 

ду

 

)

Результаты численного решения этих уравнений для плоской пластинки в предположении, что профили ско­ ростей в пограничном слое подобны во всех сечениях x = const, т. е. wx/w0=f(y/Ö) , и сопоставление их с опыт­ ными данными приведены на

рис. 1-15.

Хотя найденное реше­ ние асимптотическое, но, как видно из графика, уже при

т V

т- е- ПРИу =

= 4 х /У Я ех, скорость'практически равна скорости невоз-

Рис. 1-15. Распределение скоро­ стей в пограничном слое на пла­ стинке при ламинарном режиме течения (сплошная линия — теория).

44

мущенного потока іѵ0. Эту величину можно принять за толщину пограничного слоя для сечения х, в котором сосредоточено все влияние вязкости при ламинарном об­ текании пластинки.

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

Напряжение трения на стенке при этом течении

dw — ? *=о = Г* ----PW’xW'y,

где второй член характеризует молекулярную вязкость, связанную с молярным объемом импульсов.

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

области порядок величины скорости

равен ш * = ]/ то/р

и соответствующее число Re= y0w jv

должно быть

по­

рядка единицы; отсюда y0=v/w *.

 

 

Таким образом (рис. 1-16), в непосредственной

бли­

зости от стенки

к ней прилегает тонкая прослойка жидко­

сти, в которой

усредненная скорость мала и

меняется

по линейному закону от нуля до значения ю* =

0/Р-

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

Рис. 1-16. Развитие пограничного слоя на пластине.

45

ное, но с линейным законом усредненных скоростей. Вне области вязкого подслоя (резкой границы между вязким подслоем и остальным потоком жидкости провести, ко­ нечно, нельзя) основную роль играют турбулентные на­ пряжения pw'xw'y.

Из величин, определяющих движение жидкости с плот­ ностью р, создающее напряжение то, можно составить единственную комбинацию, имеющую размерность dwjdy,

а именно У Xolply, и поэтому должно иметь место равен­ ство

d w ___L

i / " _I°_,

d y ~ x y

V

p

где x — так называемая универсальная постоянная, кото­ рая определяется только опытом (по результатам испы­ таний в трубах).

Интегрирование дает логарифмический закон распре­

деления скоростей:

 

да = ^ -1 п ^ * -+ С ,

(1-80)

где г/ш*/ѵ — так называемое безразмерное расстояние от стенки.

Из опыта для труб: х = 1/2,4 и С=5,84.

Для плоской пластинки толщина турбулентного по­ граничного слоя

8 = 0,37//R£T.

Переход от ламинарной формы движения к турбу­ лентной зависит от числа Rex=T2j0x/v. С увеличением ско­ рости при всех прочих одинаковых условиях область пе­ рехода приближается к переднему краю хорошо обте­ каемого тела. Структура пограничного слоя конечной толщины дана на рис. 1-16.

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

46

Рис. 1-17. Интегральные соотношения получаются при­ менением законов сохранения к элементу пограничного

ABCD.

ные выше и подтвержденные опытом, показывают, что при ламинарном движении, например, уже на безразмер­

ном расстоянии порядка 1 /|/R e разница между скоро­ стью потенциального обтекания и скоростью в погранич­ ном слое практически исчезает. Эта идея пограничного слоя оказалась весьма плодотворной. Она привела к за­ мене дифференциальных уравнений для каждой точки потока в пограничном слое интегральными соотношения­ ми, которые получаются в результате применения зако­ нов сохранения к элементу длины слоя. Метод интеграль­ ных соотношений по своему замыслу, оперирующий с не­ которым образом усредненными параметрами, оказывает­ ся сравнительно малочувствительным к точности приня­ того в расчет закона распределения продольной скорости по нормали. Расчеты показывают, что во многих случаях достаточно, чтобы принятое распределение было только более или менее правдоподобным.

Изменение количества движения в объеме ABCD (рис. 1-17) должно быть равно секундному импульсу сил,

складывающихся

из силы трения r0dx

на

грани

АС и

результирующей силы давления, равной

р

dS

,

ах, как

это видно из рис.

1-17.

 

 

 

Это дает:

 

 

 

 

гг

d Г

,

d Г 2 ,

dd

 

Р З Г

J 9 * ° J y = - % - р -3F

о

 

9

'

47

или, замечая, что при р = const получаем p-f-рДОд /2 =const

и, следовательно, dpfdx — — pwa dwjdx,

5

 

5

 

 

- ^ - J p wx{w0 — wx) d y + - ^ ~ [ р(ш0 — ts»,)dj/ =

v

(1-81)

о

 

b

 

 

В дальнейших вычислениях для несжимаемой

жидко­

сти удобно первый

интеграл записать в виде

хю~0 §**, а

второй — в виде ш05*,

где

 

 

6*

=

 

 

 

 

5,

оо

 

(1-82)

 

 

 

dy.

О

Знак 6, оо означает, что с некоторого значения б прак­ тически безразлично, до какого расстояния от стенки про­ изводят интегрирование.

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

При этом вместо (1-82) получаем:

 

 

dè*

w\

(23**-1-8*) =

to

(1-83)

dx

w0

Po^o

 

 

 

 

Это основное интегральное соотношение теории по­ граничного слоя было выведено Карманом и называется уравнением импульсов.

Уравнение энергии Л. С. Лейбензон получил анало­ гичным путем, применяя закон сохранения энергии к эле­ менту ламинарного пограничного слоя. В ней фигуриру­ ет толщина потери энергии

5, оо

 

 

wx \

(1-84)

l - - f

К

wo

)

 

О

 

 

Условные толщины пограничного слоя. Толщины вы­ теснения б*, потерь импульса б**, энергии б*** значи-

48

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

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

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

При наиболее грубом приближении — линейном рас­ пределении скорости (wx= w0yl6) в пограничном слое не­ сжимаемой жидкости имеем:

6* = 6/2; б** = 6/6; 6*** = 6/4.

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

скоростей (рис.

1-18) в пограничном слое

следует, что

в месте отрыва

производная от скорости

по нормали

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

Так как для ламинарного слоя S//~ 1 / у Re; W yj W x ^ 1 / У Re

и при изменении Re изменяются только поперечные мас­ штабы длины и скорости, а движение сохраняется подоб-

4—233

49

Соседние файлы в папке книги из ГПНТБ