книги из ГПНТБ / Гидродинамика лопаточных машин и общая механика [сб. ст
.].pdfИтак, рассматривается прямоугольный капал с произвольным соотношением оторон, который характеризуется большим отяоте ■- нием длины к гидравлическому диаметру, а поток в нем - высоки ми числами Рейнольдса. Поэтому рассматривается случай стабили зированного турбулентного потока жидкости (га за). Физические свойства жидкости полагаются независящими от температуры.
Особенностью потока в прямоугольных каналах является на личие развитых 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» |
|
|
|
|
|
|
|
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 |
и |
Ш - |
; |
|
|
|
|
|
|
|
|
|
|
0Х |
|
Эу |
|
|
|
|
|
|
|
|
|
|
|
|
|
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
Полученные в результате обработки результатов расчета значе ния констант и показателей имеют следующие значения: