книги из ГПНТБ / Дорфман Л.А. Гидродинамическое сопротивление и теплоотдача вращающихся тел
.pdfГЛАВА I
ОСНОВНЫЕ УРАВНЕНИЯ
§ 1. Дифференциальные уравнения движения
Внастоящей главе будут рассмотрены уравнения для течений несжимаемой вязкой жидкости при малых темпера турных напорах, когда можно пренебречь влиянием темпе ратуры на плотность р, вязкость р и теплопроводность К жидкости.
Вдальнейшем по необходимости будут введены соответ
ствующие усложнения в уравнениях, связанные с условиями сжимаемости и интенсивного подогрева.
Для осесимметричных движений жидкости наиболее есте ственно представить уравнения движения в цилиндрической
системе координат г, ср, |
|
г. |
|
При этом |
ввиду осевой |
симме |
|||||||||
трии производные |
встречающихся в |
уравнениях величин по <р |
|||||||||||||
будут |
|
равны |
нулю. |
|
|
|
|
|
|
|
|
|
|||
|
Поэтому в случае несжимаемой вязкой жидкости при |
||||||||||||||
постоянстве |
физических |
|
констант |
р |
р, |
к уравнения |
Навье- |
||||||||
Стокса |
осесимметричного |
|
движения |
без учета массовых сил |
|||||||||||
примут вид1): |
|
|
|
|
|
|
|
|
|
|
|||||
dvr |
! |
|
|
dvr |
r |
dvr |
|
|
_ |
1 |
др |
( |
|
||
dt |
' |
|
V>' |
dr |
""Г”dz |
|
r |
|
|
ftir |
1 |
|
|||
|
|
|
|
|
|
. (d2vr |
. |
1 dvr |
|
vr . d2vr |
|
||||
|
|
|
|
|
|
|
V |
|
|
' |
7 ~dr |
T2'^^ |
(71) |
||
dvv |
|
|
dva , |
dva |
|
v„,vr |
|
|
|
||||||
|
|
|
|
|
|
|
|||||||||
dt |
+ Vr |
----- 1 |
~5--- - — |
|
|
|
|
|
|||||||
|
1 |
r |
dr |
1 |
2 |
dz |
1 |
|
r |
|
|
|
|
|
|
|
!) |
Лойцянский |
Л. |
|
Г., |
Аэродинамика пограничного слоя, |
Гостехиздат, 1941, стр. 52.
10 |
|
ОСНОВНЫЕ |
УРАВНЕНИЯ |
[ГЛ. I |
||||
дщ |
, |
dvs |
_ |
1 |
др |
|
|
|
dt |
r dr |
дг |
~ |
Р |
дг |
|
(1-1) |
|
|
|
|
|
|
1 |
дг |
\ |
|
|
|
|
|
|
г |
д& Г |
|
|
|
В этих уравнениях vr, |
|
и |
vz |
представляют собой соот |
ветственно радиальную, тангенциальную (окружную) и осевую
составляющие |
вектора скорости, р — давление, |
v — кинема |
тическую вязкость. |
|
|
К системе |
уравнений Навье-Стокса (1.1) |
присоединим |
еще уравнение неразрывности, которое для несжимаемой жидкости примет вид:
dvr |
|
vr . |
dvz _ |
(1-2) |
|
dr |
' |
г -т" дг |
|||
|
|||||
Решение систем (1.1), |
(1.2) |
требует задания начальных |
и граничных условий. Начальные условия определяются зна
чениями |
искомых величин в начальный момент |
времени |
t~t0, |
они имеют смысл для неустановившихся движений |
|
жидкости. |
|
|
Граничные условия определяют поле скоростей |
на гра |
|
ницах, |
а также давление на свободной поверхности жидкости |
и на поверхностях раздела. В частности, важное значение имеет такое граничное условие, как условие прилипания
жидкости к твердой стенке, т. е. отсутствие относительной скорости скольжения жидкости по поверхности твердого тела.
Для компонентов касательных напряжений в случае осе симметричного течения будем иметь:
§ 2] |
УРАВНЕНИЯ РЕЙНОЛЬДСА |
11 |
|
§ 2. Уравнения Рейнольдса |
|
для осредненного турбулентного движения |
||
Известно, |
что с увеличением скорости |
поток жидкости |
теряет устойчивый, упорядоченный характер, наступает не упорядоченное движение, в котором начинают играть суще ственную роль нерегулярные пульсации. Упорядоченное дви жение, наблюдающееся при малых скоростях, называется ламинарным, а неупорядоченное, наблюдающееся при боль
ших скоростях — турбулентным.
Явления, связанные с потерей устойчивости ламинарного движения, будут рассмотрены в дальнейшем на конкретных примерах потоков вокруг вращающихся осесимметричных тел. Здесь мы лишь заметим, что в случае развитого турбу лентного течения уравнения (1.1), (1.2) остаются в силе только для мгновенных действительных скоростей. Однако ввиду нерегулярности мгновенных скоростей и случайности характера пульсаций уравнения (1.1), (1.2) невозможно про интегрировать. Поэтому Рейнольдс предложил их преобра зовать в уравнения осредненного турбулентного движения. Для этого действительная скорость представляется в виде суммы двух слагаемых: осредненной скорости с составляю
щими vr, |
vr |
vz |
и |
скорости |
турбулентных пульсаций |
||||
с составляющими v', v', |
vr: |
|
|
|
|
||||
|
|
|
г |
|
z |
|
|
|
|
vr = vr-j-v'r, |
v,f = v.fy-uf, |
|
(1.4) |
||||||
Сделаем |
подстановку |
этих соотношений в уравнения (1.1), |
|||||||
(1.2) и произведем |
их |
осреднение во времени (обозначается |
|||||||
чертой сверху), полагая, |
что |
операция |
осреднения |
обладает |
|||||
следующими свойствами: |
|
|
|
|
|
||||
|
|
т'^0’ |
|
|
|
4г=4Ь ■ |
(1-5) |
||
где $ — одно |
из |
переменных |
г, |
ср, г, |
t; ср, ф— любые две |
||||
функции этих |
переменных. |
|
|
|
|
||||
Тогда получим |
следующую |
систему |
уравнений осреднен |
ного осесимметричного движения вязкой несжимаемой жидко сти— уравнения Рейнольдса — в цилиндрической системе
12 |
|
|
|
|
|
ОСНОВНЫЕ |
УРАВНЕНИЯ |
[ГЛ. I |
|||||
координат: |
|
|
|
|
|
|
|
|
|
|
|
||
д^г I |
- <4 , |
- д”г |
|
|
|
|
|
|
1 |
||||
dt 1 |
г |
dr |
1 |
г |
|
dz |
|
г |
|
|
|
|
|
__ |
1 др , |
(d2vr |
, |
1 dvr |
|
vr |
, |
d2vr\ |
|||||
|
p |
dr |
|
|
|
dr2 |
' |
r |
dr |
|
r2 |
|
dz2 ) |
dvv |
- |
dvv |
|
- |
|
dv„ |
|
vrv„ |
|
|
|
|
|
~dT^!~Vr~dr |
1“Vz ~дГ + ~~Г~ |
|
|
|
|
||||||||
|
|
|
d‘-v,f |
( |
1 |
dVy |
v,f |
' |
d2v^ |
|
|||
|
|
|
dr2 |
' |
r |
dr |
r2 |
|
dz2 / |
(1-6) |
|||
|
|
|
|
|
|
0 |
----- |
|
d |
____ |
v v |
||
|
|
|
|
— dr \ r v/ |
|
dz \ ? г/ — 2 rr <p , |
|||||||
|
|
1 |
dp |
\ |
|
\ dr2 |
I |
1 |
dr |
I |
d2ve \ |
||
|
|
p |
dz |
' |
|
' |
r |
|
dz2 |
/ |
|||
dvr ^Vr'dvz—r. |
|
|
|
|
|
|
|
||||||
dr |
r |
|
dz |
~ |
|
|
|
|
|
|
|
||
Сопоставляя систему уравнений (1.6) с уравнениями (1.1), |
|||||||||||||
замечаем, |
что |
к |
вязким |
членам |
уравнений |
(1.1) добавляются |
|||||||
члены, |
соответствующие |
турбулентным |
напряжениям трения, |
возникающим благодаря турбулентным пульсациям скорости.
§ 3. Уравнение баланса энергии
Не останавливаясь на выводе1), приведем уравнение баланса энергии для осесимметричного потока несжимаемой жидкости:
+Ч2(»’+ 2Ш+2(»3+
, |
(д^ |
v^2 , |
(dv^2 ' |
(dvs |
|
|
\ dr |
г ) |
dz ) |
\ dr + dz ] J ’ |
(1 -7> |
1)Лойцянский Л. Г., Аэродинамика пограничного |
слоя, |
||||
Гостехиздат, |
1941, |
стр. 53. |
|
|
|
§ 3] |
УРАВНЕНИЕ БАЛАНСА ЭНЕРГИИ |
13 |
где |
Т — температура, ср — теплоемкость при |
постоянном |
давлении. Член в фигурных скобках соответствует диссипации (рассеиванию) энергии в результате трения.
Вслучае турбулентного режима можно из уравнения (1.7) получить соответствующее осредненное уравнение баланса энергии.
Без учета диссипации оно примет вид:
Врассматриваемом случае постоянства физических кон стант р, ср, р, X уравнение баланса энергии можно решить,
определив сначала поле скоростей по уравнениям движения. При решении уравнения баланса энергии следует удовлетво рить заданным начальным и граничным условиям для распре деления температуры. В частности, аналогично условию при липания одним из граничных условий для температур является условие отсутствия скачка температур между обтекаемой поверхностью и прилипающими частицами жидкости.
ГЛАВА II
ЛАМИНАРНОЕ ТЕЧЕНИЕ ВОКРУГ ВРАЩАЮЩЕГОСЯ
ВСВОБОДНОМ ПРОСТРАНСТВЕ ДИСКА
§4. Ламинарное движение, создаваемое
вращающимся диском.
Точное решение уравнений Навье-Стокса
|
Пусть |
бесконечная |
плоская |
пластинка |
z = 0 вращается |
|||||
в |
вязкой жидкости вокруг |
оси |
г = 0 с постоянной угловой |
|||||||
скоростью |
(О. |
|
|
|
|
|
|
|
|
|
|
Рассматриваем движение жидкости в полупространстве |
|||||||||
2^>0. Граничные условия задачи будут иметь вид: |
|
|||||||||
|
■пг = 0, |
— |
|
и. = 0 при |
2 = 0, |
} |
(2-1) |
|||
|
|
|
|
|
|
|
|
|
||
|
Vr — Q, |
Uf = 0 |
|
При |
2 = 00. |
|
||||
|
Скорость vz |
при |
2 = оо |
не |
может быть равна нулю, так |
|||||
как диск |
действует |
в |
качестве центробежного вентилятора |
|||||||
и |
вызывает подсос, |
|
создавая |
отрицательное |
значение уг |
|||||
в |
бесконечности |
и |
радиальное |
движение |
от центра, |
в осо |
||||
бенности вблизи |
диска. |
|
|
|
|
|
||||
|
Оценим |
вначале |
толщину S |
слоя жидкости, |
увлекаемого |
диском вследствие трения1). Пусть направление, в котором вдоль диска скользит поток и которое параллельно каса тельному напряжению на стенке тот, образует с окружным направлением угол &. Тогда радиальная составляющая каса тельного напряжения равна тот sin &. Она уравновешивается центробежной силой отбрасываемого пластинкой потока, т. е. пропорциональна ргш28. С другой стороны, окружная соста вляющая касательного напряжения totcos& пропорциональна
Ч Прандтль Л., Гидроаэромеханика, перев. с нем. ИЛ, 1949, стр. 436.
§ 4] |
ТОЧНОЕ РЕШЕНИЕ |
УРАВНЕНИЙ |
НАВЬЕ-СТОКСА |
15 |
|
градиенту |
окружной составляющей |
скорости около стенки, |
|||
т. е. |
(— — знак пропорциональности). Исключая каса |
||||
тельное напряжение тст из |
соотношений |
|
|||
|
тст sin & — pro)23, |
j |
(2.2) |
||
|
|
п |
У |
j |
|
|
Jot COS it — U -g- |
|
|||
и принимая, что угол О |
не зависит от радиуса, что |
под |
|||
тверждается наблюдениями, |
мы получим: |
|
|||
|
8 |
|
|
|
(2-3) |
следовательно, толщина увлекаемого диском слоя одинакова
на всех радиусах.
Отсюда следует, что касательное напряжение
тот — ргш28 — ргш j/vio,
а момент сопротивления диска радиуса R, который пропор
ционален |
произведению |
касательного напряжения, |
площади |
и плеча, |
будет |
~ р/?4ш "j/w. |
(2.4) |
|
Л1 ~ |
Из этих рассуждений следует, что для интегрирования уравнений (1.1), (1.2) целесообразно ввести вместо z без
размерное расстояние от стенки С.— |
т. е. принять |
v |
(2-5) |
Очевидно, что составляющие vr и vv пропорциональны гш, где множители пропорциональности есть функции С, а соста
вляющая |
вследствие уравнения неразрывности должна быть |
|||
пропорциональна 8w ='|/Ча>. |
|
|
||
Таким образом, целесообразно сделать следующую замену |
||||
переменных *): |
|
|
|
|
•ur = r^F^\ |
^ = rt»G(Q, |
vs = y^H^), | |
(2 6) |
|
|
|
р = — pvuP (С). |
/ |
|
!) К а г m a n Т h., |
Laminate und turbulente Reibung, |
ZAMM, |
16 |
ЛАМИНАРНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА |
[ГЛ. II |
Подстановка (2.6) приводит уравнения (1.1), (1.2) к виду (штрихами обозначены производные по С):
F2 —G2+F'H=FZ'
2FG + G'H=G",
(2.7)
НН' = Р'-^-Н", 2F + H' = Q,
а граничные условия будут:
F = 0, |
0=1, |
Н = 0 при |
F = 0, |
0 = 0 |
(2-8) |
при |
Заметим, что при такой замене вязкие члены в правой части уравнений (1.1) сводятся только к последним слагаемым, представляющим производные по направлению, нормальному
к плоскости |
диска, так что уравнения движения тождественны, |
||||
уравнениям пограничного слоя. |
|
|
|||
Первые два уравнения (2.7) совместно с последним дают |
|||||
возможность |
определить F, |
О, |
И, |
а третье — значения Р. |
|
Если И— с |
при С-> оо, |
то можно формально построить!) |
|||
разложения по |
степеням |
для |
F, |
О и Н, удовлетворяю |
щие дифференциальным уравнениям и условиям на бесконеч ности. Нетрудно определить, что первыми членами разло жений будут
|
|
Ai + B2 |
е-2Л |
Л(^+В2) |
_3ес _ |
|
|
2с2 |
|
4с< |
|
|
|
Д 2 _|_ Д2 |
|
|
|
|
|
--Л4^(17Л2+В2)е-4О:+ •••’ |
|||
О = Ве-^ |
В(42-г№) с_3л, , |
|
|||
12В |
|
|
(2.9). |
||
|
|
|
|
||
|
|
|
|
|
|
|
|
|
' |
18c6 |
-!-•••- |
Н = — с 4------ е-с- |
2c3 |
|
|
||
|
1 |
с |
|
|
|
т |
6c5 |
° |
288c7 |
(17Л2 +fi2) e |
4e’ + • • ■ |
Остается определить |
только |
значения постоянных А, В и с. |
9 Cochran W. Q., The flow due to rotating disc, Proc. Cambr. Philos. Soc., t. 30, № 3, 1934.
§ 41 |
|
ТОЧНОЕ |
РЕШЕНИЕ УРАВНЕНИЙ |
НАВЬЕ-СТОКСА |
|
17 |
||||||||||||
|
С другой стороны, можно получить |
формальное разло |
||||||||||||||||
жение |
вблизи |
С = 0, |
удовлетворяющее |
(2.7) |
и |
граничным |
||||||||||||
условиям |
при £ = 0: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Г, |
у |
|
1 |
r2 |
|
1 |
, |
уЗ |
1 ,'iyi |
|
1 |
у5 |
, |
|
|
1 |
|
|
Т7 |
■ а.(у, |
|
С |
|
-g |
|
|
bo-, |
-gQ- п0, |
|
|
|
|
|
||||
|
|
(1 |
|
arfin |
\ |
/ |
bn |
I |
апЬп |
\ |
Г7 |
I |
|
|
|
|||
|
|
|
________re |
I |
______ ®_ |
0 0 |
|
’ |
|
|
||||||||
|
|
|
360 |
|
90 |
/ |
’К 315 |
|
1260 / |
|
|
|
|
|||||
о=1+*ос+4 |
|
|
|
(«<л - щ4 - 4 v— |
■ |
(2.10) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
f ао |
| |
^0 |
ув |
| |
( |
ао |
60 |
|
й0^0 \ р | |
|
|
|
|
||||
|
\ 90 |
' |
45 / ’ |
' |
\ 315 |
315 |
|
252 / |
|
|
‘ |
’ |
|
|||||
Н = -а^. |
|
|
|
|
|
'О’* + W |
|
|
|
|
|
|
|
|
||||
где |
а0, |
^ — неизвестные постоянные. |
находятся из условия |
|||||||||||||||
|
Значения постоянных А, В, с, а0, Ьо |
|||||||||||||||||
непрерывного |
сращивания значений |
функций F, О и Н и |
||||||||||||||||
их |
производных, |
получаемых |
из |
разложений (2.9) |
и |
(2.10). |
||||||||||||
Если задаться |
приближенными значениями |
а0 |
и |
Ьо, |
то |
с по |
||||||||||||
мощью |
(2.10) |
находятся величины |
F, |
G и Н и их |
произ |
|||||||||||||
водные для малых значений С Затем методом Адамса можно |
||||||||||||||||||
продолжить |
решения |
дифференциальных |
уравнений |
(2.7) |
||||||||||||||
к большим значениям С. Сравнением полученных значений |
||||||||||||||||||
функций и их производных в |
точках С= 1,9 и |
2,5 со зна |
||||||||||||||||
чениями, которые даются разложениями (2.9), Кокрэн1) чис |
||||||||||||||||||
ленным путем уточнил величины постоянных а0 и Ьо и нашел |
||||||||||||||||||
значения с, |
А и В. Его результаты |
дают: |
|
|
|
|
|
|||||||||||
|
|
|
ао = О,51О; Ьо = — 0,616; |
|
|
|
|
|
(2.И) |
|||||||||
|
|
|
|
с = 0,886; |
А — 0,934; |
В = 1,208. |
|
|
||||||||||
|
Значения функций F, О, И, Н', G' и Р, |
вычисленные |
||||||||||||||||
Кокрэном, приводятся в таблице 1. Первые три функции |
||||||||||||||||||
графически представлены на рис. 1. |
Из графиков видно, что |
|||||||||||||||||
расстояние от стенки, на котором окружная скорость течения |
||||||||||||||||||
понижается до половины окружной скорости стенки, равно |
||||||||||||||||||
§0,5 ~ j/v/d). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
4) См. |
сноску на |
стр. |
16. |
ГОС. |
ПУБЛИЧНА |
|
I’-X |
||||||||||
|
|
|
|
|
|
|
|
|
НАУЧН-ТЕХНЙЧЕСИАЯгУ«/| |
2 Зак. 944. Л. А. Дорфман |
БИБЛИОТЕКА СССР | |
18 |
ЛАМИНАРНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА |
[ГЛ. II |
Таблица 1
Значения функций, определяющих распределение скоростей и давления вблизи диска,
|
|
|
вращающегося в неподвижной жидкости |
|
||||
|
|
|
|
(по В. |
Кокрэну) |
|
|
|
_ |
1/ ш |
F |
О |
-н |
—р |
F' |
-G' |
|
C=z у |
— |
|||||||
|
0 |
|
0 |
1,000 |
0 |
0 |
0,510 |
0,616 |
|
0,1 |
|
0,046 |
0,939 |
0,005 |
0,092 |
0,416 |
0,611 |
|
0,2 |
|
0,084 |
0,878 |
0,018 |
0,167 |
0,334 |
0,599 |
|
0,3 |
|
0,114 |
0,819 |
0,038 |
0,228 |
0,262 |
0,580 |
|
0,4 |
|
0,136 |
0,762 |
0,063 |
0,275 |
0,200 |
0,558 |
|
0,5 |
|
0,154 |
0,708 |
0,092 |
0,312 |
0,147 |
0,532 |
|
0,6 |
|
0,166 |
0,656 |
0.124 |
0,340 |
0,102 |
0,505 |
|
0,7 |
|
0,174 |
0,607 |
0,158 |
0,361 |
0,063 |
0,476 |
|
0,8 |
|
0,179 |
0,561 |
0,193 |
0,377 |
0,032 |
0,448 |
|
0,9 |
|
0,181 |
0,517 |
0,230 |
0,388 |
0,006 |
0,419 |
|
1,0 |
|
0,180 |
0,468 |
0,266 |
0,395 |
—0,016 |
0,391 |
|
1,1 |
|
0,177 |
0,439 |
0,301 |
0,400 |
—0,033 |
0,364 |
|
1,2 |
|
0,173 |
0,404 |
0,336 |
0,403 |
—0,046 |
0,338 |
|
1,3 |
|
0,168 |
0,371 |
0,371 |
0,405 |
—0,057 |
0,313 |
|
1,4 |
|
0,162 |
0,341 |
0,404 |
0,406 |
—0,064 |
0,290 |
|
1,5 |
|
0,156 |
0,313 |
0,435 |
0,406 |
—0,070 |
0,268 |
|
1,6 |
|
0,148 |
0,288 |
0,466 |
0,405 |
—0,073 |
0,247 |
|
1,7 |
|
0,141 |
0,264 |
0,495 |
0,404 |
—0,075 |
0,228 |
|
1,8 |
|
0,133 |
0,242 |
0,522 |
0,403 |
—0,076 |
0,210 |
|
1,9 |
|
0,126 |
0,222 |
0,548 |
0,402 |
—0,075 |
0,193 |
|
2,0 |
|
0,118 |
0,203 |
0,572 |
0,401 |
—0,074 |
0,177 |
|
2,1 |
|
0,111 |
0,186 |
0,596 |
0,399 |
—0,072 |
0,163 |
|
2,2 |
|
0,104 |
0,171 |
0,617 |
0,398 |
—0,070 |
0,150 |
|
2,3 |
|
0,097 |
0,156 |
0,637 |
0,397 |
—0,067 |
0,137 |
|
2,4 |
|
0,091 |
0,143 |
0,656 |
0,396 |
-0,065 |
0,126 |
|
2,5 |
|
0,084 |
0,131 |
0,674 |
0,395 |
-0,061 |
0,116 |
|
2,6 |
- |
0,078 |
0,120 |
0,690 |
0,395 |
—0,058 |
0,106 |
|
2,8 |
0,068 |
0,101 |
0,721 |
0395 |
—0,052 |
0,089 |
|
|
3,0 |
|
0,058 |
0,083 |
0,746 |
0,395 |
—0,046 |
0,075 |
|
3,2 |
|
0,050 |
0,071 |
0,768 |
0,395 |
—0,040 |
0,063 |
|
3,4 |
|
0,042 |
0,059 |
0,786 |
0,394 |
—0,035 |
0,053 |
|
3,6 |
|
0,036 |
0,050 |
0,802 |
0,394 |
—0,030 |
0,044 |
|
3,8 |
|
0,031 |
0,042 |
0,815 |
0,393 |
—0,025 |
0,037 |
|
4,0 |
|
0,026 |
0,035 |
0,826 |
0,393 |
—0,022 |
0,031 |
|
4,2 |
|
0,022 |
0,029 |
0,836 |
0,393 |
—0,019 |
0,026 |
|
4,4 |
|
.0,018 |
0,024 |
0,844 |
0,393 |
—0,016 |
0,022 |
|
оо |
|
0 |
0 |
0,886 |
0,393 |
0 |
0 |