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

Механика жидкости и газа. В 4 ч. Ч. 2

.pdf
Скачиваний:
2
Добавлен:
24.11.2025
Размер:
2.15 Mб
Скачать
Рис. 1.6
Рис. 1.5

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

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

y

 

 

 

 

углубленному пониманию. Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

рассмотрим также этот путь.

 

 

 

 

 

 

 

 

 

 

 

Пусть жидкая частица вращается

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вокруг оси z с угловой скоростью z

 

 

 

 

 

r

M

 

 

 

 

 

 

 

 

 

 

 

(рис. 1.5). Запишем выражение для

 

 

 

 

 

 

 

 

 

x

ротора скорости в проекциях на оси

 

 

 

 

 

z

координат (рис. 1.6). Имеем

 

 

 

 

 

 

 

 

 

 

z

rot xu uz uy ;

y z

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

uy

 

 

 

 

 

ux

 

 

uz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot yu

 

;

 

 

 

 

 

 

 

 

 

 

 

z

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

ux

 

 

 

 

rot z

 

 

 

uy

 

 

ux

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

x x Рассмотрим точку M на жидкой частице.

Линейная скорость этой частицы u z r . Запишем выражения для проекций скоростей на оси координат:

ux z y ;

uy z x ;

uz 0 ,

11

откуда находим

uy z ;x

ux z .y

Таким образом

rot zu uy ux 2 z .

x y

Аналогично для двух других компонент

rot xu 2 x и rot yu 2 y .

Либо в векторной форме

12 rot u ,

что полностью совпадает с (1.5).

Движение, при котором rot u 0 , называют вихревым, если же rot u 0 – безвихревым либо потенциальным. Сказанное означает, что если течение вихревое, то движение жидких частиц происходит с вращением.

12

2. КИНЕМАТИКА ВИХРЕВОГО ДВИЖЕНИЯ

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

12 rot u.

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

2x 2y 2z .

Движение, при котором величина вихря скорости не равна нулю, т. е. rot u 0 , называют вихревым. Если же rot u 0 , то движение безвихревое (потенциальное).

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

dx dy dz .

x y z

Вихревая трубка – аналог трубки (поверхности) тока, т. е. это поверхность, образованная вихревыми линиями, прове-

13

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

Интенсивность вихря

Понятие интенсивности вихря достаточно абстрактно и вводится чисто математически. Напомним, что потоком векторного поля называют интеграл вида

u ndA .

A

Поскольку вихрь скорости (ротор) есть вектор, то вместо u можно подставить rot u , что и приводит к понятию интенсивности вихря, т. е. интенсивность вихря – это поток вектора вихря скорости:

i rot u ndA.

A

Эту формулу, используя очевидное соотношение rot u n rotn u , можно переписать как

 

 

 

 

i rotnudA .

 

A

 

 

 

 

Имея в виду, что

rot

 

2 (рис. 2.1),

 

u

 

можем записать

 

 

 

Рис. 2.1

i 2 ndA 2 n dA.

A

 

 

A

 

 

 

14

Используя формулу Гаусса–Остроградского и переходя от интеграла по поверхности к интегралу по объему, получим

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

x

 

 

 

 

z

i 2 n dA 2 div dV 2

x

y

z

dV .

A

V

V

 

 

 

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

 

 

 

 

1

 

u

 

u

 

 

 

 

x

 

 

 

 

z

 

y

;

 

 

 

 

 

 

2

 

y

z

 

 

 

 

 

 

 

 

 

y

1

 

u

x

u

 

 

 

 

 

 

 

z ;

 

 

 

 

 

 

 

2

 

z

x

 

 

 

 

 

 

1

 

u

y

u

 

 

 

 

z

 

 

 

 

 

 

x .

 

 

 

 

 

 

 

2

 

x

y

 

 

 

 

 

 

 

 

 

Получим

 

 

 

2

 

2

 

2

 

2

 

2

1

 

uz

uy

ux

uz

uy

 

 

 

 

 

 

 

2

 

y x

x z

z y

x y

x z

 

Следовательно, можно записать

n dA 0.

A

 

 

2

 

0 .

 

ux

 

 

 

 

 

 

y z

 

(2.1)

15

Применим (2.1) к вихревому шнуру (рис. 2.2). На боковой поверхности n 0 , так как вектор направлен по касательной к поверхности. Поэтому можем записать

n1dA1 n2 dA2 0 ;

 

A1

A2

 

n1dA1 n2 dA2 .

 

A1

A2

 

Если допустить, что в пределах сечения

 

n const , то

 

Рис. 2.2

n1A1

n2 A2 .

Либо в общем случае

 

 

A cons t ,

(2.2)

т. е. это своеобразное «уравнение неразрывности» в интегральной форме для завихренности. Полученный результат носит название теоремы Гельмгольца о вихрях (второй теоремы Гельмгольца), которую можно сформулировать следующим образом: интенсивность вихревого шнура на всей его протяженности остается постоянной. Из выражения (2.2) следует и другой весьма важный вывод. Поскольку произведение A

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

16

Рис. 2.3

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

Циркуляция скорости

Для введения понятия о циркуляции скорости воспользуемся методикой, предложенной Н.Я. Фабрикантом. Несомненным ее преимуществом является то, что она позволяет ввести понятие циркуляции не формально математически, а исходя из достаточно простых и ясных физических предпосылок.

Рассмотрим крыловой профиль, находящийся в равномерном потоке воздуха. Как известно, на профиль в этом случае будет действовать подъемная сила F (рис. 2.3). Физически наличие этой силы можно объяснить лишь тем, что давление

под профилем p1 больше, а давление над профилем p2 меньше, чем давление на некотором удалении от него (давление невозмущенного потока), которое обозначают p . Это позволяет утверждать, что под крыловым профилем скорость u1 u , а над ним u2 u . В данном случае u – скорость невозмущенного потока.

Теперь из скоростей u1 и u2 вычтем скорость u , т. е. получим разности u1 u и u2 u . Это действие приводит к

17

понятию возмущенного потока, т. е. движения, которое возни-

кает в среде из-за того, что в нее внесено инородное тело. По

существу, это реакция потока, в данном случае обусловленная

тем, что в ней появился крыловой профиль. Установим теперь

направление потоков возмущения. Под профилем u1 u , и он

направлен против скорости u , над профилем соответственно

наоборот. В результате появляется циркуляционный поток,

направленный по часовой стрелке, как это показано на рис. 2.3.

Чтобы характеризовать этот поток количественно, вводится

понятие циркуляции скорости по замкнутому контуру.

dl

 

u

Рассмотрим замкнутый контур C, по-

 

казанный на рис. 2.4. Пусть в произ-

 

 

 

 

 

 

 

 

вольной точке M скорость равна u . Со-

 

M

 

ставим скалярное произведение u dl ,

 

 

 

где dl – направленный элемент дуги.

L

 

 

Циркуляцией скорости называют кон-

 

 

турный интеграл вида:

Рис. 2.4

 

u dl .

Обратим внимание на структуру этого соотношения. Оно

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

гда говорят, что циркуляция – это своеобразная «работа» век-

тора скорости. Имея в виду, что u ux,uy,uz и dl dx, dy, dz , по

правилу скалярного произведения получим

uxdx uydy uzdz .

Для плоского течения:

uxdx uydy .

(2.3)

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

18

ность вихря. Действительно, из формулы (2.3) следует, что для определения циркуляции достаточно знать проекции скорости, нахождение которых не связано с существенными трудностями. Однако при этом пока открытым остается вопрос о том, существует ли связь между циркуляцией скорости и интенсивностью вихря.

Теорема Стокса

В движущейся жидкости рассмотрим вихревое поле и выделим в нем малый замкнутый контур со сторонами dx и dy (рис. 2.5). Пусть в начале координат скорости будут ux и u y .

 

 

z

 

 

 

 

Запишем выражение для эле-

 

 

 

 

 

 

 

 

ментарной циркуляции по

 

 

O

 

uy

C

этому контуру, имея в виду,

 

 

 

 

что поток двумерный:

 

ux

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

ux

 

ux dy

d uxdx uydy .

 

 

 

 

 

 

y

 

x

u

 

 

uy

B

 

 

Рассмотрим контур OABC.

 

y

x

dx

 

 

 

 

 

 

 

 

 

Если вдоль OA скорость ux ,

 

 

 

 

 

Рис. 2.5

 

 

 

 

 

 

 

 

 

то вдоль CB ее приращение

 

 

 

 

 

 

 

 

составит

ux

dy , и аналогично вдоль AB

uy

dx . Это следу-

y

x

ет из выражения для полного дифференциала скорости, например:

dux ux dx ux dy .

x y

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

d u

 

 

 

 

u

 

 

 

 

 

u

 

 

 

dy .

x

dx u

y

 

y dx dy u

x

 

x dy dx u

y

 

 

 

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

19

Раскрывая скобки и выполнив сокращения, получаем

uy d

x

 

 

 

 

dA.

ux dxdy 2

 

y

 

z

 

 

 

 

 

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

2 ndA i .

A

Это и есть формула Стокса, показывающая, что циркуляция по произвольному контуру равна сумме интенсивностей (напряжений) вихрей, пронизывающих поверхность, натянутую на контур.

20