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

2982

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
2.73 Mб
Скачать

Рис. 4.4

Как видим, в обоих случаях с увеличением числа Re убывает коэффициент сопротивления. Но в турбулентном пограничном слое коэффициент сопротивления убывает значительно медленнее, чем в ламинарном пограничном слое.

Как показывают опыты, более точно коэффициент сопротивления пластины в случае турбулентного пограничного слоя выражается формулой

Cx тр

 

0,074 .

 

 

Re0,2

 

 

 

В пределах 106 Re 109 часто используется следующая

формула:

 

 

 

 

Cx тр

0,445

 

 

 

 

 

lg Re 2,58

 

Как видим, теория достаточно хорошо подтверждается экспериментом, и тем самым оправдывается введенная гипотеза об аналогии между турбулентными течениями по трубам и в пограничном слое плоской пластины.

61

4.5. Расчет смешанного пограничного слоя для плоской пластины

Характер пограничного слоя, образующегося при обтекании потоком жидкости плоской пластины, существенно зависит от режима обтекания, определяемого числом Re.

При сравнительно небольших числах Re вдоль всей пластины образуется ламинарный пограничный слой, расчет которого рассмотрен.

При очень больших значениях числа Re практически вдоль всей пластины образуется турбулентный пограничный слой.

Этот случай также рассмотрен.

При числах Re от Re 105 до Re 5 106 в начале пластины образуется ламинарный пограничный слой, который, начиная с некоторого места пластинки, разрушается и переходит в чисто турбулентный пограничный слой.

Рассмотрим случай, когда на передней части пластины ОА образуется ламинарный пограничный слой, а за ним на участке AB – турбулентный.

Этот случай впервые рассмотрен Прандтлем.

Рис. 4.5

Сделаем два упрощающих предположения:

1. Переход от ламинарного пограничного слоя к турбулентному происходит мгновенно в точке А; расстояние ОА обозначим

через X кр .

62

2. Изменение толщины турбулентного слоя, распределение скоростей и касательных напряжений в нем, аналогично тому, которое было бы, если бы турбулентный слой начинался не от точки А, а от передней кромки О.

Обозначим через X тр силу трения всей пластины длинной в

предположении, что пограничный слой на всем ее протяжении турбулентный.

Тогда, чтобы получить силу трения X тр для смешанного слоя,

нужно из Xтр вычесть силу трения наружного участка слоя ОА,

считая его турбулентным, и прибавить к полученной разности силу трения этого же участка ОА, считая его ламинарным, т.е.

Xтр Xтр Xтр ОА Xтр ОА ,

где X тр - сила трения ламинарного участка ОА.

Силу трения переднего участка пластины шириной b при ламинарном пограничном слое можно записать в виде:

 

 

V2

 

 

C

 

0

X

 

b,

x тр

2

тр

 

 

 

 

 

 

 

 

а при турбулентном слое в виде

 

 

 

 

 

 

 

 

 

 

V2

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

0

 

 

X

 

b .

 

 

 

 

 

 

x тр

2

 

 

тр

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Их разность равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2

X тр

X тр ОА

 

X тр ОА

 

 

 

 

С х тр

0

X кр b

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2

 

 

 

 

V

2

 

 

 

 

 

 

 

C х тр

0

X кр b

 

 

0

 

С х тр C х тр X кр b

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V 2

Разделим обе части полученного равенства на

 

 

0

b .

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда для изменения C х тр

пластины из-за наличия ламинарного

участка будем иметь

63

Reкр

 

 

 

 

С х тр

C х тр

X кр V0

 

C х тр

С х тр

C х тр

X кр

 

 

A

,

 

 

 

 

 

 

 

 

 

 

V0

 

 

 

Re

 

 

 

 

 

 

 

 

где через А обозначена следующая величина:

A

Сх тр

Cх тр

XкрV0

Сх тр

Cх тр

Re кр

 

и число

Re кр

Xкр V0

– называется критическим.

 

 

 

Опытным путем установлено, что для гладких пластин величина А=1700. При увеличении шероховатости пластины или степени турбулентности набегающего потока уменьшаются критическое число и величина А. Для шероховатых пластин А=300.

Таким образом 300 A 1700.

Итак, для величины коэффициента трения плоской пластины в случае смешанного пограничного слоя получаем следующее выражение:

Cх тр

 

0,074

 

A ,

 

Re0,2

 

Re

 

 

 

или

 

 

 

 

 

 

 

 

Cх тр

0,445

 

 

 

A .

lg Re 2,58

 

 

Re

 

 

 

Из изложенного следует, что сопротивление трения пластины будет тем меньше, чем больше длинна ламинарного слоя, т.е. чем дальше точка А перехода ламинарного слоя в турбулентный отстоит от передней кромки пластины.

Рассматривая обтекание плоской пластины как первое приближение к обтеканию крыльевого профиля, приходим к выводу что с целью уменьшения сопротивления трения на профиле выгодно иметь возможно больший участок ламинарного пограничного слоя. Иначе говоря, точку А перехода ламинарного пограничного слоя в турбулентный выгодно отодвигать как можно дальше назад по потоку. По этому принципу сконструированы специальные крыльевые профили, так называемые ламиниризованные профили, которые имеют значительно меньшее сопротивление.

64

4.6. Пограничный слой на криволинейной поверхности

При обтекании криволинейной поверхности скорость V0 на

внешней границе пограничного слоя будет величиной переменной, зависящей от координаты х.

Как следует из уравнения Бернулли, давление на внешней границе пограничного слоя, а стало быть, и в пограничном слое также будет функцией координаты х.

Рассмотрим поток, обтекающий криволинейную поверхность, например, профиль крыла (рис. 4.6).

Обозначим через V и P соответственно

скорость и давление набегающего потока вдали от профиля крыла (на бесконечности).

На верхней поверхности профиля скорость вначале возрастает (до некоторой точки М), а затем убывает.

На основании уравнения

Рис. 4.6

Бернулли давление, наоборот, будет вначале убывать, а затем возрастать.

В точке М скорость будет иметь максимальное значение, а давление – минимальное.

Следовательно, частицы жидкости в пограничном слое около рассматриваемой криволинейной поверхности будут двигаться при

наличии градиента давления dPdx как отрицательного, так и положительного по знаку.

Этот факт существенно отличает пограничный слой около криволинейной поверхности от пограничного слоя вдоль плоской

пластинки, где dPdx 0 .

Учитывая эту особенность пограничного слоя около криволинейной поверхности, можно выяснить причины отрыва

65

потока от обтекаемого тела и образования вихрей, срывающихся с обтекаемой поверхности.

При течении вязкой жидкости касательная и нормальная составляющие скорости в точках поверхности обращаются в нуль.

На некотором небольшом удалении от поверхности течение жидкости мало отличается от течения идеальной жидкости, при котором нормальная составляющая скорости на поверхности тела обращается в нуль, а касательная отлична от нуля.

Имея это в виду, приходим к выводу, что изменение касательной составляющей скорости вдоль нормали к поверхности тела должно иметь вид, изображенный на рис. 4.6 (точка А).

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

Кроме этого, было установлено, что градиент давления dPdx для случая криволинейной поверхности отличен от нуля и,

что важно, в кормовой части тела перепад давления направлен в сторону, противоположную основному движению жидкости

dPdx 0 , так как в этой части давление возрастает по потоку.

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

Это дополнительное притормаживание, очевидно, особенно сильно скажется на частицах, движущихся в непосредственной близости от тела, так как их кинематическая энергия в сравнении с другими частицами и так уменьшена силами вязкости.

В результате касательная скорость ряда частиц обратится в нуль на некоторой линии SD (рис. 4.6), а в области DSE касательная скорость изменит знак, т.е. возникает возвратное движение жидкости в пограничном слое.

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

Два столкнувшихся слоя жидкости, сорвавшись с поверхности тела, свертываются и образуют вихри.

66

Точка поверхности тела, начиная от которой поток срывается с обтекаемого тела, называется точкой отрыва пограничного слоя. Условно принимают, что

слева от точки отрыва S производная:

dVx

tg 0,

dy

y 0

 

справа от точки S производная

dVx

tg 0,

dy

y 0

 

в самой точке S производная

dVx

tg 0.

dy

y 0

 

Таким образом, точка S отрыва пограничного слоя характеризуется равенством нулю производной

dVx

0.

dy

y 0

 

Отсюда следует, что и напряжение трения

 

dVx

0

0

dy

 

 

y 0

 

 

в точке отрыва.

Полученный результат хорошо подтверждается экспериментом.

Если в поток поместить цилиндр, то можно обнаружить срывающиеся с него вихри, образующие вихревую дорожку.

Если пограничный слой отсасывать внутрь цилиндра, то срыв вихрей можно устранить.

Таким образом, теория пограничного слоя позволяет управлять пограничным слоем.

67

5. ОБТЕКАНИЕ ТЕЛ ПОТОКОМ ИДЕАЛЬНОЙ ЖИДКОСТИ И ГАЗА

5.1.Уравнение Эйлера

Впотоке идеальной жидкости выделим элементарный

параллелепипед (см. рисунок). Пусть R x, y, z – единичная массовая сила, а X ,Y, Z – проекции единичной массовой силы. На параллелепипед будут действовать инерционные силы:

dV

dxdydz;

dVy

dxdydz;

dV

dxdydz.

x

 

z

 

 

 

 

dt

 

dt

 

dt

 

Так как жидкость идеальная, то из всех поверхностных сил будут действовать только силы давления. Пусть давление в точке A

P x, y, z .

Тогда

ввиду

малости граней проходящих

через точку А, и не зависимости давления от ориентации площадки, на эти грани будет тоже действовать давление Р.

Очевидно,

 

что

P P x

dx, y, z .

1

 

 

Разлагая

P1

в ряд Тейлора,

получаем:

 

 

P P x dx, y, z P x, y, z

 

dP

dx.

 

 

 

 

1

 

dx

 

 

Аналогично

 

 

 

 

 

P P x, y dy, z P x, y, z

dP

 

dy;

 

2

dy

 

P P x, y, z dz P x, y, z

dP

dz.

 

3

dz

 

68

 

 

 

 

 

Согласно принципу Даламбера сумма проекций на любую ось всех сил, действующих на параллелепипед, равна нулю.

В результате в проекции на ось х получим

 

 

 

 

 

 

 

dVx

 

 

 

 

 

 

 

 

dP

 

 

 

 

 

 

X dxdydz

 

 

 

 

 

 

 

 

 

dxdydz

 

 

 

 

 

 

 

dxdydz

0,

 

 

 

dt

dx

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dVx

 

 

1 dP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

dx

 

 

 

 

 

 

Аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dVy

 

 

1 dP

 

 

 

 

dV

 

 

1 dP

 

 

 

 

Y

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

z

 

Z

 

 

 

 

.

 

dt

 

 

 

 

dy

 

 

 

 

 

 

dt

 

 

 

dz

Уравнения Эйлера содержат 5 неизвестных: V , V , V , , P.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

z

К этим уравнениям можно добавить уравнение

неразрывности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

dV

 

dVy

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

z

 

0.

 

 

 

 

 

 

dt

 

 

 

 

 

 

dx

 

dy

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пятое уравнение получить из принципов механики нельзя. Его можно получить, исходя из физических свойств жидкости или газа. Например, const для несжимаемой жидкости. Для

сжимаемой жидкости P C

и

Cp

– для сжимаемого газа,

CV

 

 

 

где Cp и C V – удельные теплоемкости при постоянном давлении и объеме.

5.2. Уравнение Бернулли

Дифференциальные уравнения движения жидкости в форме Эйлера в общем виде не интегрируются.

Только в частных случаях, когда движение жидкости потенциальное или установившееся, можно проинтегрировать уравнения движения жидкости.

Рассмотрим случай установившегося движения. Представим, что

69

1 dP

 

dD

;

1 dP

 

dD

;

1 dP

 

dD

,

 

 

dx

 

dx

 

 

dy

 

dy

 

 

dz

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

где D x, y, z– искомая функция. Уравнения Эйлера имеют вид:

dV

 

1 dP

dVy

 

1 dP

dV

 

1 dP

x

X

 

 

 

;

 

Y

 

 

 

;

z

Z

 

 

 

.

dx

 

 

dx

dy

 

 

dy

dz

 

 

dz

Так как движение жидкости установившееся, то траектории и линии тока совпадают.

За время dt частица жидкости пройдет по траектории путь ds , так что

ds Vdt .

Спроектируем элементарное перемещение частицы вдоль линии тока на оси координат.

Тогда получим

dx Vx dt; dy Vy dt; dz Vz dt.

Умножим уравнения Эйлера на элементарные перемещения и сложим.

В результате получим:

Vx dVx VydVy Vz dVz

 

 

Xdx Ydy Zdz

dD

dx

dD

dy

dD

dz.

 

 

 

 

dx

dy

dz

Справа в скобках стоят полные дифференциалы, т.е.

Xdx Ydy Zdz d ,

и уравнение можно переписать в виде:

 

d

 

V 2

 

d

dD,

 

 

2

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

d

V 2

 

D

0,

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

V 2

С .

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Причем,

так

 

как

 

элементарные

перемещения

рассматривались

вдоль

линии

 

тока, то и

интегрирование

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]