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

книги из ГПНТБ / Дорфман Л.А. Гидродинамическое сопротивление и теплоотдача вращающихся тел

.pdf
Скачиваний:
28
Добавлен:
29.10.2023
Размер:
10.79 Mб
Скачать

§ 18]

СРАВНЕНИЕ

С ЭКСПЕРИМЕНТАЛЬНЫМИ данными

79

Интересно то

обстоятельство, что показатель

степени

в этой формуле равен соответствующему показателю в фор­ муле сопротивления гладкой плоской пластинки1). Кривая, соответствующая формуле (4.35), хорошо согласуется с имеющимися опытными данными (рис. 3).

Нетрудно убедиться, что для толщины пограничного слоя общую формулу можно представить в виде

у = 0,05+ 5 (lg R)-2’7,

(4.36)

а для коэффициента а получится

 

а

0,33+3 (lg R)~2’*5.

(4.37)

Уменьшение величины а с увеличением числа R подтверж­

дается опытами по

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

относи­

тельного движения жидкости в пограничном слое вращаю­ щегося диска2).

§ 18. Сравнение с экспериментальными данными

Сравнение теоретических значений коэффициента момента сопротивления с опытными приводилось выше; оно изобра­ жено на рис. 3.

Сопоставим теперь профили скоростей для турбулентного режима, полученные замером, с теоретическими. Характерные результаты замеров изображены на рис. 18.

Вычислим вначале характерные толщины *о и **о для окружной составляющей скорости. В случае степенного про­

филя Кармана,

согласно (4.3) и (4.8):

 

 

G =

v,„

4-

/ v хТ

/

г \

—=1—£7,

S = 0,525г (-з- )

•(?=_),

 

Г»

 

\ Г2ш /

\

6 )

поэтому

 

 

 

 

«

 

1

 

 

1

V = J -^dz = Zj G(^)rf?=g S = 0,0656r^)

.(4.38)

оо

J)См., например, формулу (100) на стр. 664 книги Л. Г. Лойцянского «Механика жидкости и газа», Гостехиздат, 1957.

2) Buss man F., Die Orenzschichtbewegung auf rotierender Scheiben, Forschung auf d. Oebiete d. Ing. Wes., t. 21, № 9, 1931.

80

ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА

[гл. IV

Аналогично вычисляется

®:*о

 

 

«

 

 

-

=/**8

0(1 — O)dz = -^ 3 = 0,051r (tv)5 ■ <4-39)

 

о

 

 

 

Для величины

получаем:

 

 

 

Н*

о

 

 

 

= у= 1,285,

 

в то время

как

замеренное1) значение № при

R = -г~ =

= 323 000 оказалось равным 1,4.

 

Для логарифмического профиля скоростей, согласно (4.17):

G =-----In $,

поэтому будем иметь:

о =у, о =_^1-т},

1—^1—yJ, (4,40)

Согласно формуле Гольдштейна (4.25) для того же числа

R = 323 000

будем

иметь

Н* =1,24. В то же время

расчет

по формулам

§

17

для

уточненного

решения дает

более

близкое к опытному значение числа ,*/7

Н = 1,295.

 

Вычислим теперь для того же числа Рейнольдса R = 323 000

соответствующие

профили скоростей

и нанесем на

график

(рис. 23), беря в качестве безразмерного расстояния от диска величину $ = .г**/8 Как видно, окружная составляющая, опре­ деляемая степенным и логарифмическими профилями, близка к экспериментальным значениям!). Для радиальной соста­ вляющей при степенном профиле скоростей получается боль­ шое отличие от замеренного распределения радиальных ско­ ростей. Наилучшее согласование получается для случая уточ­ ненного решения с логарифмическим профилем, за исключе­ нием области, непосредственно примыкающей к стенке до максимума F.

Логарифмический профиль дает также хорошее согласова­ ть

ние с опытными данными по углу Ф8 = arctg = , составляемому г > S °

1) См. сноску на стр. 19.

§ 18]

Сравнение с экспериментальными данными

81

Рис. 23. Распределение радиальных и окружных составляющих скорости вблизи вращающегося диска при турбулентном режиме; дг=2100 об/мин,

/• = 146 мм (R—323 000.)

6 Зак. 944. Л. А. Дорфман

82

 

ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ

ВОКРУГ ДИСКА

[ГЛ. IV

вектором

касательной

скорости на внешней границе погра­

ничного

слоя с окружным направлением.

 

 

Для степенного профиля Кармана

 

 

 

 

у

 

 

 

Ф8 = arctg~ — arctgа —= arctg 7а — 48°36/,

 

 

z->8 и

А

 

 

для

 

 

1-67

 

 

решения Гольдштейна

 

 

 

 

Ф8 = arctg

= arctg а = arctg 4г = 18°22',

 

 

 

z->&

и

6

 

в то время как опытные значения Ф8 изменяются в пределах от 18 до 23°. Уточненное решение для логарифмического профиля скоростей дает такие же значения Ф8. Например, при R = 323 000 получаем а = 0,373, т. е. Ф6 = 20°25', что удовлетворительно согласуется с замерами (см. рис. 18).

Опыты Коба и Саундерса1) также показали (см. рис. 35 на стр. 117), что имеется согласование теоретических данных по Карману с замерами полной касательной скорости при турбулентном режиме.

§ 19. Влияние шероховатости поверхности диска

Шероховатость поверхности вращающегося диска проя­ вляется в заметном увеличении трения. Однако, как и в слу­ чае течения в трубках и вдоль плоской стенки, если высота бугорка шероховатости значительно меньше толщины лами­ нарного подслоя /г<^ол, то трение не увеличивается, так как бугорки обтекаются в этом случае без отрывов и вихреобразований. Это режим без проявления шероховатости.

Если же высота бугорков k имеет тот же порядок вели­ чины, что и 8Л, то начинает проявляться влияние шерохова­

тости— переходный или промежуточный режим-

Наконец, существенно проявляется шероховатость в слу­ чае, когда бугорки шероховатости выходят заметно за пре­ делы ламинарного подслоя, &^>олОтрывное обтекание бугорков приводит к тому, что структура потока и трение становятся независимыми от числа Рейнольдса и зависят лишь от относительной шероховатости k/R. Это режим развитой

шероховатости.

*) Cobb Е. С. and Saunders О. A., Heat transfer from a rotating disk, Proc. Roy. Soc, ser. A, t. 236, № 1206, 1956.

§ 19]

ВЛИЯНИЕ

ШЕРОХОВАТОСТИ ПОВЕРХНОСТИ ДИСКА

83

Для

расчета

сопротивления вращающегося шероховатого

дискаJ) на режиме развитой шероховатости обратимся к ло­

гарифмическому

профилю скорости, который наблюдался

в шероховатых трубах:

 

 

— = Д1п^- + В, Л = 2,5, 5 = 8,5,

(4.41)

где ks

обозначает высоту бугорка равнозернистой и

макси­

мально

плотной,

так называемой песочной шероховатости.

на кривых соответствуют индексам функций.

На внешней границе пограничного слоя скорость дости­ гает значения

ч |г=8 = U = М In

+ 5г»„

 

так что величина Y, определяемая

формулой (4.20),

здесь

принимает вид:

 

<4-42)

г=^==1пг+4-

Kg

zi

 

Для интегралов, входящих в систему уравнений (4.1),

(4.2), справедливы формулы (4.22).

До подстановки их

в уравнения (4.1), (4.2) аппроксимируем их степенными фор­ мулами зависимости от /г8/В с помощью формулы (4.42). Как

х) Дорфман Л. А., Журнал технической физики, т. 28, № 2,

6*

84

ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА

[гЛ. IV

видно из рис. 24, где кружками обозначены точные значения

величин, а

прямые показывают результат аппроксимации,

в интервале

In = = 2 = 5 получаем следующие формулы:

(4.43)

Подставив эти выражения в интегральные соотношения импульсов (4.1) и (4.2), после сокращения на множитель <о2р получим:

0,064 А. (г3 8&з’19 8-°д9а2) -0,114 г28°-72й®’28 =

 

 

= — 0,0091аг3£°,288 ~0,28 /1 Ч-a2,

0,067 A (A 8fe°’19 8~O,1S) =

(4.44)

 

dr v

ь

'

 

 

 

= 0,0091 Л®’28 8~®’28 У1

4- а2.

Решение ищем в

виде

 

 

 

a — aoksr^, 8 = Pofes/-1-'1,

(4.45)

причем показатели р и м находим из условия, чтобы в урав­ нениях (4.44) сократились множители, являющиеся степе­ нями ks. Нетрудно определить, что искомыми значениями являются

51 = ^ = 0,043.

В результате после дифференцирования и сокращений по­

лучаем из (4.44) уравнения

 

О,236а2ро’81 — 0,11 4 о’72 = — 0,0091аор(Г0,28/1 + а2,

I „

а0 = 0,0287Ро-1’09/1 + «2-

J

§ 19]

ВЛИЯНИЕ

ШЕРОХОВАТОСТИ ПОВЕРХНОСТИ ДИСКА

85

Отсюда имеем:

 

 

 

 

 

 

а0 —0,512 (1 + а2)-0’022,

'

(4-47)

 

 

 

1

;

 

 

ро = 0,0712(1+а2)а^.

 

 

Величина (1 —а2)

0,022

пренебрежимо мало отличается

от 1,

поэтому

с достаточным

приближением можно считать

 

 

 

 

а0 = 0,512.

 

 

Например, даже при таком большом значении r/ks, как 100, разность (1-(-а2)-0,022— 1 равна ~0,01. Таким образом, в формулах (4.47) с большой точностью можно полагать

(4.48)

Момент сопротивления двух сторон вращающегося диска будет тогда равен

- з

22И = 4к/?2р J vrv4 dz

= 4тг№р [2a/W . о,0335/г0,19 Г0,19]г=л

что дает для коэффициента момента

/Ь ч 0,268

Г

/ь ч 0,0861 °>365

Cjf = 0,1015(^.)

[1 +0,262

. (4.49)

Нетрудно заметить (рис. 25), что эту формулу можно упро­ стить, представив в виде

 

 

 

,

/й.\0’273

.

 

 

(4.50)

 

 

 

у Cjf=0,108(^-)

 

 

Так

как

интервал In-г- от 2 до

5

 

соответствует

интер-

валу

 

 

D

то указанные формулы справедливы

125 + -£-< 3000,

в этом

 

«g

Зависимость

толщины пограничного

интервале Rfks.

слоя

в

этом,

исходя из

(4.45) — (4.48),

определяется

фор­

мулой

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

J

г

В

Г

Ik \°’086 Пз.оэ / k \°>№3

-

(4-51)

(

1^0,07121 1 +0,262(^J

J

(+)

86

ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА

[гл. IV

или упрощенно

(4.52)

Некоторые экспериментальные данные о трении шерохо­ ватого диска, вращающегося в большом объеме жидкости, приведены в работе А. А. Канаева *). Стальной диск радиуса /? = 0,135 м и толщиной /> = 0,01 м был покрыт искус­ ственной песочной шероховатостью с высотой бугорка

Рис. 25. Зависимость коэффициента момента сопротивления шероховатого диска от величины /?/frg.

ks = 0,45 мм. Диск прокручивался в воде, заполнявшей со­ суд больших размеров. При этом, как видно из рис. 4 ци­ тируемой работы, максимальное достигнутое значение коэф­ фициента момента трения диска и обода равно с^ = 0,027. Если исключить трение обода, исходя из формулы

""I-

приведенной в рассматриваемой работе, то для коэффициента

трения

двух сторон диска

получим значение са~ 0,0228

D

 

 

формуле (4.50) дает очень близ-

при -г- = 300. Подсчет по

кую величину

— 0,0229.

 

Заметим, что данные опытов Канаева по трению шерохо­

ватого

стального

диска о ртуть не согласуются с форму­

лой (4.50).

 

 

1) Канаев А.

А., Журнал техн, физики, т. XXIII, вып. 2, 1953.

§ 20]

ВЛИЯНИЕ

ОБДУВА ДИСКА ПРИ ТУРБУЛЕНТНОМ РЕЖИМЕ 87

Это

вполне

объяснимо горизонтальным расположением

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

Кокрэна (2.17)]; 2— для турбулентного режима при гладком диске [формула Кар­ мана (4.10)]; кружками указаны опытные данные Теодорсена и Регира для /?/£g==1205.

На рис. 26 представлен результат замеров Теодорсена и

Регирах) для диска,

покрытого песочной шероховатостью

с относительной высотой бугорка шероховатости k

1

Полученные значения

см несколько выше подсчитанных по

формуле (4.50). Возможно, что здесь сказывается влияние обода.

§ 20. Влияние обдува диска при турбулентном режиме

Применим для расчета2*) метод Кармана (§ 16). Для этого вначале выведем интегральные соотношения погранич­ ного слоя вблизи диска.

!) См. сноску на стр. 21.

2)Truckenbrodt Е., Turbulente Stromung an einer angeblasenen rotierender Scheibe, ZAMM, 1954, № 4/5,

88

ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА

[ГЛ. IV

Если

применить обычные соображения теории

погранич­

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

уравнениях

из

всех членов,

зависящих от

вязкости-, сохра­

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

на­

правлению нормали к плоскости диска; тогда получим:

 

 

 

dvr

и®.

 

 

 

dv

 

1

dz

1 др

)

 

 

 

 

> дгVr-^-----г

- +дг^г^—р=дг-- з-----р dr

I

 

4 5

vr

dv,„

 

v„vr

 

 

 

dv„

 

1

йт,

 

I

I

 

 

dr

H---- 1*----- h

г

dz

=-----3-2-.

 

J

 

 

r

1

r

'

 

 

p

dz

 

 

 

 

При этом

оставшиеся

члены,

зависящие

от вязкости,

мы

заменили, согласно (1.3)1), соответствующими производными от напряжений трения. В таком виде система уравнений по­

граничного слоя справедлива

как для ламинарного, так и

для турбулентного пограничного слоя.

Проинтегрируем первое из

уравнений (4.53) в направле­

нии,

перпендикулярном к диску, от поверхности диска (д = 0)

до

2 =

где о— толщина пограничного слоя. При этом

член

1 др

в соответствии с

уравнением Бернулли

 

 

 

I ?U2

, Р^2

 

 

 

----

1-= const,

в

котором

U — радиальная,

a W — осевая составляющие ско­

рости внешнего потока,

и формулами (2.29) можно записать

в виде

 

 

1

др__

., dU

Р

dr

dr '

Тогда получим, умножив еще на г:

h

h

h

f{rV^~rU^)dz+r fv^dz~ fv"dz = rzРr

0

 

0

0

 

(4.54)

 

 

 

 

 

где zr— радиальная

составляющая

напряжения трения на

стенке.

 

 

 

 

 

 

dv,

 

dvr

л

1) Пренебрегая

членом —по сравнению с

 

вблизи стенки.

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