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

книги из ГПНТБ / Гидродинамика лопаточных машин и общая механика [сб. ст

.].pdf
Скачиваний:
12
Добавлен:
23.10.2023
Размер:
5.18 Mб
Скачать

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

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

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

дых

поверхностях

(р и с .і). Выделенные

зоны потока

(турбулентное

ядро

и ламинарный

слой) сопрягаются

по некоторой

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

рую

можно назвать

уравнением границы

ламинарного

слоя, я энв -

ние которой является решающим для речения тепловой части зада­ чи.

 

В

турбулентных слоях распределение скорости хорошо описы­

вается

степенной

формулой Прандтля

 

 

 

 

 

 

 

 

 

 

 

 

W

-

-

W

(I)

 

где

W

-

текущее

энппение

скорости;

 

А -толщина турбулент­

ного

слоя;

у

-

линейная

координата,

 

нормальная

к омываемой

поверхности;

\ѴД -

значение

скорости

на

внешней границе

турбу­

лентного

слоя.

 

 

 

 

 

 

 

 

 

 

 

Поскольку

рассматривается стабилизированный

потек,

толщи­

ны турбулентных слоев для прямоугольного канала равны: по оси

Х -1

 

,

по оси

у - h . .

Направление координатных

осей

по­

казано

 

па

рно. 1,

причем t

- полуширина капала/

К

-полу-

f-О* ^

 

 

 

отношения

 

 

wwuia

капала. Замет им ,• что

 

 

рассматриваются в дальнейшем как безразмерные координаты. Та­ ким образом, пс аналогии с уравнением (I) может быть записано следующее выражение для распределения скорости в турбулентном ядре потока:

 

 

 

w

= w * x'l y'1,

 

 

 

(2 )

 

 

ЩѲ

VV.f

- скорость

в центре

канала

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

 

(х=Р,

у = h ) .

 

 

 

 

 

 

 

 

 

 

 

 

Формула (?)

справедлива,

естественно,

для левой нижней чеч

верти канала (рисЛ ) ,

то

еоть

при

OstX&l

и

О é-ÿ < |

,

рассмотрением которой в силу симметрии структуры потока можно

 

ограничиться.

 

 

 

принимает равным 1 /7 ,

 

 

 

 

Показатель

п

обычно

хотя в

лите

-

ратуре имеются указания,

что

он иожет

варьироваться в пределах

 

r t= ÿ - r

в

зависимости

от

числа

Re

[й] .

 

 

 

 

На рис 2а и 2в приведены для сравнения картины изотах дей­

ствительного турбулентного течения в прямоугольном

канале

[ 5 ]

по данным

Иикурадзе и

рассчитанная

по формуле (2)

при

 

 

Я-1/7 . Как следует из данных рисунков, формула (2) правильно от­ ражает физику явления, хотя действительный профиль скоростей от личается большей заполненностью, что является следствием наличия

вторичных

течений,

направленных к углам канала,

Однако,

учитывая

приближенный характер

решения задачи,в

дальнейшем

иопо чъэуется

формула

(2 ).

 

 

 

 

 

 

 

 

 

Формула распределения скорости удобно выразить через средин

расходную скорость

жидкости

 

 

 

 

 

 

 

 

 

 

 

VVcj> =

G

 

 

(8 )

 

 

 

 

 

4P

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь '

G

-

весовой

расход жидкости

через канал,

у

-удельный

вес жидкости.

 

 

 

 

 

 

 

 

 

Цели

пренебречь

площадью канала,

занятой ламинарный слоем

ввиду малости,

то выражение

для среднерасходной

скорости может

быть получено

путем

интегрирования (?)

яо всей

площади

канала

 

 

 

 

 

I

г

 

 

 

 

 

 

 

 

W e? = J j

 

W d x d(j

 

(ч)

 

 

 

 

 

(п + і)г

 

 

 

 

 

 

 

 

 

 

 

 

о о

Р и с I

W - Û.9GA

W -0,913

W =0,B55

W =0,811

W = 0,748

W = 0,655

g

-

і)0

 

 

Выразив из уравнения (О

W*

через

ѴѴср , получим окончат

тельную формулу распределения скорости в турбулентном ядре

W=(n +l)2 w<LpX,lÿr -

(5)

Уравнение границы ламинарного

слоя ищется исходя из следующих

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

лентное

происходит при одинаковом отношении

коэффициентов вяэкооти

(турбулентного

и молекулярного) по всему

профилю границы.

 

На

рис.З

схематично показаны структура

потока в рассматривае­

мой канале, направления

координатных осей

и касательные напряжения

/f

,

действующие ва

элемевт жидкости.

С

учетом обозначений

рис.З для составляющих скорости справедливы оледуюшие поотноиения:

ѴѴХ= Wy = D ;

 

 

 

Wz-VV(X,y)

и определяется no (5),

 

ÔVV*

ôWx

d\Vx

= o

Ce)

Зх

" äy

э х

 

a w a _ aw«

awa

= o

 

3x

~ 3y

 

32

 

 

^ ï = o .

Составляющие вязкого напряжения, действующие на элемент,опреде ляютоя по известным формулам работы (5) , которые с учетом уравне - ний (б) выглядят следующим образом:

Тхж=<*

дѵѴк

ЗѴѴг \ -

Ыiî'zT

öZ

ах

)

^ dX

 

Т у г -

3Wu

ЭѴѴ2 \ =

оѴѴ

М - --

 

ѳ г

ö'y

'

э а

 

 

Полнив

вязкое напряжение, действувдее на элемент

где

- коэффициент динаиичесиой вязкоотн жидкости.

D

соответствии о теорией турбулентности составляющие тур<3у -

лентного напряжения, дейотвующкгв ne влемепт, определяются соотно-

■ѳвиякя ( 5 ]

!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'і'* г

= - p w ^ v v i

, .

 

 

 

 

 

 

'Гуг

= - p W j j W i

,

 

 

 

 

где р

- плотность жидкости,

w it, ѵѴу,

Wz

-пулюаниошшз оостявдя-

юцие скоростей

ѴѴ*

,

WjAVj.

 

 

 

 

 

 

Пульсационнье составляющие

Wj

и

VV'ÿ

пропорциональны произ­

водный

ÜW

и

Ш -

;

 

 

 

 

 

 

 

 

 

 

 

Эу

 

 

 

 

 

 

 

 

 

 

 

 

 

wi‘=î

Э\Ѵ

О,MX aw

 

(9)

 

 

 

 

 

'а х

 

 

ѳх" ’

 

 

 

 

w!,4 sayaw

o.-a Щ- a w-

 

 

 

При написании уравнений (9) в соответствии о теорией турбулент­

ности принято, что длина пути смешения

L

и

в направлениях

соответствующих коорш чатных осой пропорциональна координатам от

твердых

поверхностей

X

Я

ij

, причеи коэффициент пропорциояадь-

лости равен 0,4 ( по ькспвриѵізнтальным данный Кариона).

 

Пульсаиионная составляюшвя

Wz

имеет

порядок максимальной

из величин

ѴѴІ

или

 

Wÿ

и может быть определена как

их оуі«а

 

 

 

W

 

 

 

 

 

 

 

 

с ™ ,

Используя уравнения (ѳ), (9) и (ІО) полное турбулентное напря­ жение, действующее на элемент, определяется следующим обрааом:

(H )

По аналогии о уравнением (7) дли турбулентного напряжении может быть записано соотношение:

Г г -

л,т ( Ш +

)

(I?)

L

^ (й х

9у />

 

где (М - эффективный коэффициент турбулентной вязкости,

іИэ сопоставления ( I I ) и (1?) легко найти выражения для коэф­

фициента турбулентпой

вязкости (ѵт

или для отношения

:

Z .

- - '

IГ (

«

' )

(13)

о1

= ' 0,16(М "

эуѵ

9W

 

 

 

 

ЭХ

Эу

 

Выражение (13) преобразуется о учетом принятого закона распре­ деления скорости (5 ), так что окончательно оно выглядит:

 

 

 

 

 

 

 

х п У"

(Ій)

 

£ = - D . * n ( H aRe( l t j ) - = ^

 

 

 

 

 

 

 

 

 

h

 

 

Здесь число

Re определено по средперасходной

скорости

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

диаметру

канала

d r

.как это

обычно принято

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

тепло гидравлических

процессов в

каналах, то есть

 

 

 

 

,PWq>dr

 

 

 

 

 

Re = — JJ---- I

 

 

 

 

 

d

-

M

 

 

 

 

 

 

 

d r ' ~ l + b

 

 

 

 

 

Полученное выражение (.Ій) является в конечном счете искомым

уравнением границы ламинарного олоя,

воли положить его левую часть

равной

некоторой

постоянной величине,

 

причем эта величина должна

иметь

порядок 10,

так как переход ламинарного течения в турбулент­

ное происходит при (ит >> ()'

 

.

 

 

 

 

Запишем уравнение (ій ) для границы ламинарного слоя на уров­

не середины пысоты канала,

то

есть

при

ÿ = 1 :

 

 

- не

Здесь индекс „Л” оанпчает, что параметр берется на трави­ не ламинарного слои, а индекс “ К. " - на уровне середины высоты канала ( или ÿ = I ), Таким образом, ХлЬ =Хль/£,

абсцисса границы ламинарного слоя или относительная толщина лами­

нарного слоя, па уровне середины

высоты канала. При написании вто­

рого приблизительного равенства (Т5) учтено, что

Хм, «s: K/h

 

Используя допутепие о постоянстве значения

(лту^и '

вдоль

всей

границы ламинарного слоя,

приравниваем

(і'О

и (І5)

и полу­

чаем окончательное выражение для уравнения

границы ламп парного

слоя:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(іб )

 

Остановимся на определении величины X-.lt

.Ее легко опре­

делить,

если предположить, что

в районе

середины

высоты

канала

точение одномерно и изменение

скорости

имеет место только по

оси

X

. Это допущение оправдано для

высоких

каналовЛроме

того, принятый степенной закон распределения скорости обегпечи-

вает

довольно заполненный профиль

в оредней части.

 

 

 

 

Тогда напряжения трения на твердой поверхности

Lw

в

рассматриваемом сечении ( y=h

),

с

одной стороны,

определи-

ртся

линейным распределением скорости

в

ламинарном слое

 

 

 

 

 

 

 

 

 

UV)

 

где

 

ѴѴм,

- скорость

жидкости на

границе ламинарного

слон при

и,

о другой

стороны,

законом Блазиуса

 

 

 

 

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

Tw = 0,03225 ( пт i)3'5f W 'P'Àw

(ІѲ)

Ae

 

Вакон трения Блазиуса проверен экспериментально в отроком дна-

паэоне чисел

Re

(вплоть

до Re = l07 )

и справедлив для

гидрав­

лически

гладки*

поверхностей.

 

 

 

 

 

Приравняв

(І7) я (І8)

и реинв полученное

равенство

относи­

тельно

ѴѴдн

.найдем

 

 

 

 

 

 

 

Wflh = 0,03225(0 +if '5 | т ^

О* f r

/■»лЬ

( 19)

 

 

 

 

 

г\е

 

 

 

С другой

стороны,

скорость

\ѴлЬ

мощно найти

из (5 ), полонив

Х=Хд|,

И

i( = h

:

 

 

 

 

 

 

 

 

 

WAh = ( n + l ) 2VVcP

('

 

 

(20)

Приравняв (19) и (20) и решив полученное равенство относи - тельно ХдН , найдем окончательное выражение для толщины лаг мипарного слоя на уровне середины высоты капала

Хд|\ г

0,00806 ( n - tl /’5 Re" 5.(|

(21)

Для показателя

1

 

П =- 7

 

 

\\h

221

 

 

Re675R F

 

 

 

 

Возвращаясь

к уравнению (іб ), следует отметить,

что оно не

допускает явного решения относительно одной из координат грани­ цы ламинарного слоя в зависимости от другой. Поэтому была сде­ лана попытка аппроксимировать решение (16) приближенными зави­ симостями. С этой целью граничная кривая была разбита на Четыре участка (ри с.І).

I

Тогда применительно к участку I, например, уравнение (іб ) удобно переписать в виде

Ггг)

Зависимость (22)

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

граничной

кривей

X* с

линейной

координатой

ÿ

вдоль высоты

канала

при

ij 4 I

 

Учитывая,

что Хл^У

 

, особенно при

і]

стремящихся

к I,в т о ­

рым слагаемым в скобках (22) ыокно пренебречь и рассматриѵіать поду­

ченное уравнение

как решение

(22) относительно

первого приближения

 

 

 

 

V

_ п

- -2~

 

 

 

 

 

 

 

 

 

 

Л*

 

ntl

 

 

 

 

 

 

 

 

 

.

Я д ,," 3

 

Хл/х ль

 

 

 

 

 

 

Дей отвита дьное значение

 

можно рассматривать

как

ре­

шение

первого приближения,

 

умноженное

на некоторый попра*очный

коэффициент

£

,

то есть

 

 

 

 

 

 

 

 

 

 

 

 

 

ТГ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 22)

 

Предварительный анализ

позволяет установить общий вид поправ­

ки

б,

Из уравнений (21) я (22)

ясно, что

£,

зависит

от

R e ,

^

 

и

§

> KP0“ 6 того,

предполагается,

что эта зави­

симость

носит стейе.ной характер. При

Cj ^

I поправка

f

,

так как в атом случае рвение первого прйближенкя мало отличает-

оя от действительного

зьачѳния

Хд/X h .

-Таким ибіэтэом

 

е, = И

с ,he

(1 /ь )

g

,

(2 0

где

константа С, и

показатели степени

L, , И.,

N, устанавли­

ваются по результатам обработки точных решений (22).

 

Применяя аналогичные рассуждения к участку П,

уравнение

(іб )

переписывается в

виде:

 

1

 

 

 

Іі

.1

 

 

 

 

 

За

1

h+ l

 

(25)

 

XAh

' X

+ h

 

 

 

который дает зависимость абсциссы граничной кривой

ОТ линейной координаты вдоль ширины канала

X

.Аналогично участ­

ку I

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

 

 

 

 

 

 

 

 

 

 

JL - Ц . ' Г

"

5 ’ ""' с

 

(2 6 )

 

 

 

Хль

" U

)

 

х

Ьг

 

где

 

 

 

 

и

/ В \м г

_ Hz

 

 

 

 

 

 

 

 

 

 

 

 

£2=і + CaRe

 

 

X

 

(27)

 

Учитывая симметрию граничной кривой на участках III и ІУ участ­

кам I и П

соответственно, легко

найти,

что

на

учаетке ID( U ÿ s 2 )

 

 

 

 

 

 

 

 

-

 

 

причем

 

е! = 1+ C,Re (* )"

 

 

 

 

 

 

 

 

 

и на участке

ІУ

 

 

 

 

 

 

 

У

 

 

h

 

\ n+î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ХлП = г - І І I

 

я

 

 

 

 

Значения

констант

С,

и

с 2

 

 

 

Мі,

N2

были получены в

результате

обработки точных решений (22)

и (25). Решение этих уравнений проводилось численным методом нй~ БЭСМ-4. При этом геометрия канала и числа Re варьировались в широком диапазоне

0 , 2 é - £ - £ 2

5 - ю Ч р .е * И 0 6

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

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