Механика жидкости и газа. В 4 ч. Ч. 2
.pdf
этом и заключается гидромеханический смысл равенства нулю дивергенции вектора скорости.
Полученную выше связь между поступательной и вращательной скоростями жидкой частицы можно получить и более коротким путем, представляющим определенный интерес. Разные подходы к одному и тому же вопросу способствуют его
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
Понятие об интенсивности вихря является весьма важным, но, к сожалению, непосредственное определение этой величины экспериментальным путем связано с непреодолимыми трудностями. Кроме того, если пытаться распространить это понятие на вихри конечных размеров, то по аналогии со средней скоростью пришлось бы вводить понятие о средней угловой скорости, что связано с определенными трудностями чисто математического характера. Поэтому гидромеханики избрали другой путь, заменив это понятие другим, более удобным для целей практики.
Циркуляция скорости
Для введения понятия о циркуляции скорости воспользуемся методикой, предложенной Н.Я. Фабрикантом. Несомненным ее преимуществом является то, что она позволяет ввести понятие циркуляции не формально математически, а исходя из достаточно простых и ясных физических предпосылок.
Рассмотрим крыловой профиль, находящийся в равномерном потоке воздуха. Как известно, на профиль в этом случае будет действовать подъемная сила 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
