§ 20] ВЛИЯНИЕ ОБДУВА ДИСКА ПРИ ТУРБУЛЕНТНОМ РЕЖИМЕ 89
Используя уравнение неразрывности (1.2), исключим v2.
Почленное интегрирование второго интеграла в соотно шении (4.54) тогда дает
hh.li rfv^dz=~u J №+vr')dz+j
о о о
Подставив это выражение в (4.54), получим после пере группирования членов искомое интегральное соотношение импульсов для меридионального направления:
Г |
h |
h |
h |
|
г J(U—vr)vrdz |
+ r^=f (.u~vr)dz-}-f &dz = r^. |
|
L |
о |
о |
0 |
|
|
|
(4.55) |
Заметим, что верхний предел интегрирования может быть |
|||
взят |
любой от 8 до |
со, так как |
подинтегральные выраже |
ния |
равны нулю при |
h ^>8. |
|
Аналогично из второго уравнения (4.53) выводится инте гральное соотношение импульсов для окружного направле ния; при этом получается уравнение (4.2). Заметим, что в слу чае отсутствия обдува диска (£7 = 0) уравнение (4.55) пре вращается в уравнение (4.1). Теперь выберем соответствующие профили скоростей в пограничном слое.
Профиль окружных скоростей в пограничном слое будет выражаться той же формулой (4.3), что и для случая без
обдува:
1
В то же время, для того чтобы радиальная составляющая на внешней границе пограничного слоя удовлетворяла усло
вию (2.29):
р\ |
= ± |
3=5 |
Га 1г_3 а ’ |
90 |
ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА |
[гл. IV |
профилю радиальных скоростей следует придать |
вид: |
|
|
|
'7 = [“(1-|)+^т](тГ- |
O-S6) |
||
При a — Q |
получаем формулу (4.3) для F. |
|
|||
Как и ранее |
[см. (4.13)] вблизи диска (г -> 0) |
|
|||
|
|
|
vr |
Ъг |
|
|
|
|
----------- = — |
— = а, |
|
|
|
|
— |
|
|
и имеют место формулы (4.14) |
для составляющих сил тре |
||||
ния. Вычисление |
интегралов, входящих в соотношения им |
||||
пульсов |
(4.55), (4.2) дает |
|
|
||
|
6 |
|
|
|
|
J* |
v^dz — (0,207ю2а2 -|~ 0,266асва + 0,304а2) г28, |
|
|||
■о |
|
|
|
|
|
5 |
|
|
|
|
|
J* v/vlfdz — 0,0681 (га))2а§ + 0,0292а<о г2о, |
|
||||
о |
г |
|
|
|
|
|
= 0,0278 (гы)2 5. |
|
|
||
J v*dz |
|
|
|||
о |
|
|
|
|
|
Подставляем эти |
значения в соотношения (4.55), (4.2), |
пред- |
|||
8 |
|
виде |
|
|
|
ставив — в |
|
|
|
||
Тогда получится после почленного деления одного уравнения
на другое, что величина р исключается, и будем |
иметь ква |
дратное уравнение для определения а: |
|
а2 + 0,0294 — а = 0,0262-4-1,0557 (-Y. |
|
Таким образом, |
|
а = ]/”0,0262+ 1,0559^ — 0,0147 |
(4.57) |
Величина р определится из второго интегрального соотно шения (4,2), которое принимает вид:
А |
,* |
4),1342-^ + 0,3131а = 0,0225(1 +а2)8 р |
§ 20] |
ВЛИЯНИЕ ОБДУВА ДИСКА ПРИ ТУРБУЛЕНТНОМ РЕЖИМЕ |
91 |
|||||||||||
Отсюда получаем: |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
р = 0,0481 |
|
(1 + °2)% |
|
|
|
(4.58) |
|||||
|
|
0,1342— 4-0,3131а |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
(О |
‘ |
|
|
|
|
|
|
а для составляющих напряжения трения в соответствии с |
фор |
||||||||||||
мулами |
(4.4) и (4.6) будем |
иметь: |
|
|
|
|
|
|
|
||||
/а>г2\0,2 |
т |
|
|
|
л о/ |
|
а |
|
|
|
\0,а |
||
— I |
р (<ог)2 |
= 0,048 la (1-Н а2) |
\ |
0,1342-4-0,3131а |
) |
, |
|||||||
\ ч / |
|
' |
1 |
' |
|
со |
1 |
|
|
’ |
|||
/№\0,2 |
|
|
|
|
0„/ |
а |
|
|
|
\0,2 |
|||
— -^5 = — 0,0481 (14-а2)0’8 |
0,1342 — 4~0,3131а) . |
||||||||||||
\ v / р (cor)2 |
|
|
' |
1 ' |
\ |
со |
1 |
’ |
|
/ |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
(4.59) |
|
Интегрируя т?, |
получим формулу для с№. Для |
случая |
|||||||||||
а = 0 |
получаем, естественно, результат Кармана. |
Теперь |
|||||||||||
можно |
вычислить |
коэффициент |
увеличения |
тг{. |
в |
результате |
|||||||
обдува: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
('сф)а =о |
|
|
|
|
|
|
||
представленный на |
рис. |
7 |
(стр. |
31). |
Как видно, |
при неиз |
|||||||
менной угловой скорости крутящий момент быстро увеличи вается с увеличением скорости набегающего потока, быстрее, чем для случая ламинарного режима.
Заметим при этом, что для случая малых ы, особенно вблизи предельного случая отсутствия вращения (<о == 0), вы веденные выше соотношения неудобны для непосредственного использования. Поэтому их необходимо преобразовать, видо изменив параметры, а именно: если обозначить
(4.60)
то получим следующие формулы:
а* = |/~0,0147^4- 1,0559 —0,0147,
(4.57а)
=* |
0,0481 |
0,13424- *0,3131a . |
(4.58а) |
|
|
|
92 |
ТУРБУЛЕНТНОЕ |
ТЕЧЕНИЕ |
ВОКРУГ ДИСКА |
[ГЛ. IV |
|||||
Вычисления показывают, |
что |
для |
очень |
малых |
значений |
по |
|||
лучается асимптотическое |
решение |
|
|
|
|
||||
|
|
_ 0,226 |
|
|
|
(4-61) |
|||
|
|
CiI~ |
|
|
|
|
|
||
В цитированной работе Труккенбродта ') решения при |
|||||||||
водятся |
в виде, |
удобном |
для получения |
результатов и |
для |
||||
других |
степенных профилей скоростей |
в |
пограничном |
слое |
|||||
вращающегося |
диска. Это |
необходимо |
|
для расчетов |
при |
||||
больших значениях чисел Рейнольдса, где уже не сохраняется степень 7? Для профиля скоростей.
!) См. сноску 2) на стр. 87.
ГЛАВА V
ТЕПЛООТДАЧА ВРАЩАЮЩЕГОСЯ В СВОБОДНОМ ПРОСТРАНСТВЕ ДИСКА
§ 21. Нагревание вязкой жидкости вращающимся диском. Точное решение уравнения баланса энергии
Используя значения скоростей из точного решения урав нений Навье-Стокса для вращающегося диска (§ 4), можно определить поле температур, решая уравнение баланса энер
гии (1.7) при заданных граничных условиях. |
|
|
Впервые такое решение получено |
И. А. Кибелем 1), |
ко |
торый заметил, что уравнение (1.7) |
можно удовлетворить, |
|
если распределение температур искать в виде |
|
|
T(z, r)=Tl(z)-\-rzT2(z), |
(5.1) |
|
где 7\ и Т2 — функции одного г. Действительно, вставляя (5.1) в уравнение (1.7) и сделав замену переменных (2.6) для ско ростей, получим, собирая члены при одинаковых степенях г, следующие два уравнения:
т" — Рр< + 2/М = - F'2 — G'2, ] |
|
|
, |
> |
(5.2) |
т" — РН^ = — 4т2—4Д2 — 2АГ, |
J |
|
относительно переменных |
|
|
’.<о=^лю. >.(0=^;-,» |
(5.3) |
|
гдеР = —?-------число Прандтля, а штрихи обозначают про-
А.
изводные по переменной С.
Ч Киб ел ь И. А., Нагревание вязкой жидкости вращающимся диском, Прикл. матем. и мех., XI, 1947, № 6.
94 |
ТЕПЛООТДАЧА ВРАЩАЮЩЕГОСЯ ДИСКА |
[ГЛ. V |
Следует заметить, что распределение (5.1) соответствует |
||
лишь |
определенному виду граничных условий — квадратич |
|
ному распределению температур по радиусу диска, так как,
согласно |
(5.1): |
|
|
|
7’(О) = 7’1(О) + Т2(О)А |
(5.4) |
|
В случае, когда число |
Прандтля равно единице, |
задача |
|
решается |
в квадратурах. |
При этом функция T2(zj |
опреде |
ляется в |
конечном виде через составляющие скорости, па |
||
раллельные плоскости диска. Чтобы в этом убедиться, запи шем первое уравнение (5.2) (при Р=1) в виде
+ 2Лта = т" + F'2 + G'2,
умножим первое уравнение (2.7) на F: FF'H-4- FF2 — FG2 = F"F,
а второе уравнение (2.7) на G-\-k, где k — произвольная постоянная:
(G 4- k) G'H-\- 2FG (G 4- k) = (G + k) G",
и сложим полученные три уравнения. Тогда будем иметь:
яh +1 2 + G2) + ko] + 2F [t2 +1 (G2 + F2) + kG] =
=[т2 + |(О2 + ^2) + /гОр.
Решением этого уравнения является выражение
^ + 4(G2 + ^2) + AG = 0. |
(5.5) |
Постоянная k определяется по заданному граничному условию.
Если диск теплоизолирован, то
-^■ = 0 при С= 0,
Лг
и из (5.5) следует
G(0) G' (0) + &G' (0) = 0,
т. е. & = —1, и получаем решение
t2 = G —1(G2 + F2). |
(5.6) |
§ |
21 j наГрёвайиё |
вйзКой жидкости вращающимся диском |
93 |
||||||
|
В случае, если |
температура |
стенки постоянна, т, |
е. |
|
||||
|
|
|
т2 = 0 |
при |
С = 0, |
|
|
||
то |
k = — у и |
распределение безразмерных температур |
т2 |
||||||
имеет вид: |
|
^ = ^0 —i(G2+F2). |
(5.7) |
||||||
|
|
|
|||||||
|
В общем случае квадратичного распределения темпера |
||||||||
тур (5.4) по радиусу |
диска, когда |
|
|
|
|||||
|
|
|
|
* (0) = Тю4- г2х2й, |
(5.8) |
||||
распределение температур |
т2(£) |
имеет |
вид: |
|
|
||||
|
|
|
т2(С) = ^)-Н20О(С). |
(5.9) |
|||||
Следует заметить, что решение |
(5.5) |
предполагает, |
что |
от |
|||||
счет температур |
выбран таким, |
чтобы |
т2(оо) = 0. |
из вто |
|||||
|
После того |
как т2 определено, можно найти |
|||||||
рого уравнения |
(5.2): |
|
|
|
|
|
|
||
|
|
|
с |
\ |
|
|
|
|
|
|
|
(У Н Л I |
X |
|
|
|
|
||
о/
|
|
|
|
с |
|
/ |
с |
\ |
|
|
|
D — 4 J* (т2 —|— 3F2) exp I |
— *J |
Н Л] Л |
л. |
||||
|
|
|
|
о |
\ |
о |
/ |
|
|
|
|
|
|
|
|
||||
где постоянные С и D определяются из граничных условий. |
|||||||||
При вращении теплоизолированного |
диска |
|
|
||||||
|
|
|
УЭ- = О |
при |
С = 0, |
|
|
||
поэтому £) = 0, |
а |
условие |
т1(со) = 0 дает |
|
|
||||
ОО |
|
с |
\ |
|
|
|
|
|
|
C=4jО |
(оу н л1 |
I х |
/ |
|
|
|
|
||
|
- |
с |
|
|
с |
\ |
Л d(^7,6. |
||
|
X |
J” (т2 + ЗЛ2) exp I |
— у Н I |
||||||
|
-о |
|
|
\ |
о |
/ |
|
|
|