книги из ГПНТБ / Дорфман Л.А. Гидродинамическое сопротивление и теплоотдача вращающихся тел
.pdf§ 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) Пренебрегая |
членом —по сравнению с |
|
вблизи стенки. |