Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
3367.pdf
Скачиваний:
8
Добавлен:
15.11.2022
Размер:
4.53 Mб
Скачать

4.2. Кинематика жидкой частицы (движение жидкой частицы в общем виде)

Поставим перед собой задачу.

 

 

 

 

 

 

 

 

 

Пусть известна скорость жидкости

 

в точке

 

 

, оп-

ределенной радиус-вектором .

 

 

 

 

 

 

 

 

 

 

 

 

Требуется найти скорость

в некоторой другой точке,

отстоящей от

 

на некотором малом расстоянии

 

. Как

это сделать?

 

 

 

 

 

 

 

 

 

 

 

 

 

Применяя формально правила дифференциального

 

 

=

+

=

+

.

 

 

 

(4.7)

 

 

 

 

 

 

 

 

 

 

 

 

исчисления, можно записать равенство

 

 

 

 

 

 

 

 

 

Замечаем,

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

 

 

представляет собой

 

 

 

отношение двух бесконечно малых векторов

 

 

и .

 

 

 

Раньше мы договорились, что

отношение двух век-

 

 

 

 

 

 

 

 

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

Отыскание этого тензора требует специальных рас-

суждений.

( )

Раньше мы получили (см. формула 3.47), что тензор можно рассматривать как производную от вектор-

ной функции по векторному аргументу .

Рассмотрим более детально тензор

 

, в кото-

 

 

 

 

ром под вектором

 

будем понимать скорость( .)Для этого

разложим тензор

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

с равенством

( )

(3.11) на симметричную и антисимметричную часть:

69

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ( )

 

=

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

+

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

+

 

 

 

1

 

 

 

 

 

 

2

 

 

+

 

 

 

 

+

 

1

+

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

21

21

 

 

. (4.8)

 

 

 

 

 

 

 

 

 

21

 

 

 

1

0

 

 

 

 

21

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

0

 

 

 

Выясним, какую физическую трактовку можно дать этим двум составляющим тензора.

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

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

В результате рассматриваемая частица изменит свою форму и через некоторый промежуток времени δt приобретет вид, показанный на рисунке штриховой линией (рассматриваем сечение этой частицы плоскостью ).

70

Рис. 5. Геометрическое изображение частицы

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

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

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

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

Примечание: такой подход к деформации частицы основывается на переносе уже известных представлений из механики и сопротивления материалов в гидродинамику и распространении их на движение жидкой частицы.

71

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

Рис. 6. Последовательное превращение частицы

Вообразим, что первоначальный контур частицы MNKL за интервал времени δt перешел в контур MABD путем трехи последовательных превращений.

1. Вначале частица растянулась или сжалась вдоль осей .

2.Затем она повернулась как единое целое на угол δγ.

3.Затем произошел её сдвиг (сплющивание) на удвоенный угол δθ.

Замечание: выбор таких трех последовательных

превращений частицы условен. Никаких оснований для

72

этого нет, т.е. мы выдвигаем гипотезу, что растяжение, поворот и сдвиг деформируют частицу и т.о. переводят контур MNKL в контур MABD.

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

одновременно.

 

 

, ,

точку внутри исходной частицы с

Возьмем любую

координатами

 

по отношению к точке М.

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

 

Посмотрим, какой величине будет пропорциональна

скорость удаления этой точки от точки M в

направлении

 

 

 

 

оси x вследствие расширения или сжатия частицы.

 

 

 

 

 

 

Если скорость в точке M принять равной

 

,

то

смещение выбранной (рассматриваемой) точки по

отноше-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нию к точке M за время δt оказывается равной величине:

 

а скорость

δx =

+

 

,

 

 

 

 

 

 

удаления или приближения к точке М вдоль ко-

при

0

, т.е.

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

ординатной оси x определяется пределом отношения

 

Любая точка

= lim→0 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частицы будет приближаться или удаляться от точки M вдоль оси x со скоростью, пропорцио-

нальной величине .

Очевидно, что, повторяя аналогичные рассуждения, найдем коэффициенты пропорциональности скоростей смещения по осям y и z в виде множителей.

73

и .

Кроме растяжения или сжатия, на скорости удаления выбранной точки от точки M по координате x скажутся деформация сдвига и поворот частицы в пространстве как твердого тела.

Очевидно, что эти добавки скорости пропорцио-

нальны быстроте изменения углов

 

 

и .

 

 

 

частицы

Непосредственно

из схемыδθ

деформации

 

 

 

 

 

 

Выразим

 

 

δa = δj + δθδβ; δβ = δj δβδθ

 

 

вытекают такие соотношения:

 

 

 

 

 

.

 

(*)

 

 

 

 

 

 

δθ

 

 

 

 

 

 

 

.

Сами углы

 

и = найдем2 ; δθ,

=

2

 

 

 

 

 

 

 

углы

 

и

δa+

 

 

 

 

δa−

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рассматривая вращение

граней MN и ML

относительно точки M.

 

 

 

 

 

δa δβ

 

 

 

 

 

 

 

 

 

 

 

Угловые смещения, как известно, равны линейным

смещениям, деленным на радиус поворота.

 

 

 

 

Линейное смещение точки N относительно точки M

а угловое –

 

 

 

 

 

=

 

δt,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вдоль оси x за время δt составит величину:

 

 

 

 

Так же

 

 

 

 

 

=

 

 

δt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соответственно:

 

.

 

 

 

 

 

 

 

 

Знак минус δβ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

находится угол:

 

 

 

 

 

 

 

 

 

здесь поставлен потому, что точка L движется относительно точки M в отрицательном направлении по оси Y.

74

Подставляя найденные значения углов в формулы

(*), деля интервал времени

к нулю, получаем скорости

скоса частицы в плоскостиxyδиt вращения её вокруг оси z, т.е.

Скос:

=

lim

 

=

1

 

 

 

.

 

 

2

+

Вращение:

→0

 

 

 

1

 

 

 

= lim

 

=

 

2

 

.

Знак минус

→0

 

 

означает вращение частицы

 

 

перед

 

вокруг оси z по часовой стрелке, т.к. мы уславливаемся

считать положительным направлением вращения против часовой стрелки (от оси x к оси y).

Скорость смещения любой точки внутри элементар-

ного кубика относительно точки M в направлении оси x

 

1

 

1

 

 

.

оказывается пропорциональна величинам:

 

Мы

; 2 +

 

; 2

+

рассмотрим только деформацию и вращение в плоскости xy.

Легко убедиться, что смещение по оси x вызывается аналогичными превращениями с частицей и в плоскости xz, а в плоскости zy коэффициенты пропорциональности

равны нулю.

, 21 + , + ;

(**)

 

 

Поэтому полный набор коэффициентов такой:

(***)

 

0, 21 + , 21 ;

75

они

 

( )

 

 

Сравним найденные коэффициенты с первой стро-

 

тензора

 

в формуле (4.8) и сразу увидим, что

кой

тождественно

совпадают с его элементами. Причем

группа (**) дает верхнюю строку симметричного тензора, а группа (***) – антисимметричную.

Очевидно, что если провести аналогичные рассуждения для отыскания коэффициентов, определяющих сме-

щение точек

жидкой частицы относительно её полюса

( )

 

 

(точки M) по осям y и z, то получим весь набор элементов

тензораПри этом.

первая группа (**) опишет деформацион-

ное движение частицы, а вторая группа (***) – вращательное.

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

деформационным и вращательным движением её частичек,

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

что выражается симметричной и кососимметричной частя-

ми тензора

 

.

 

 

 

 

 

 

 

 

 

 

Симметричная часть тензора

 

поэтому назы-

def , аТ.о.

 

 

 

 

 

деформации

 

 

 

 

вается тензором

 

скоростей

 

( )

 

 

 

 

 

 

кососимметричная – вихрем скорости –

 

 

 

.

и с учетом (4

.=

( )

 

= def +

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

(4.9)

= +

 

= + def +

.

(4.10)

 

 

7) получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (4.10) представляет собой важное для гидродинамики утверждение, известное под названием теорема Коши-Гельмгольца: движение любой точки жидкой частицы можно рассматривать как результат сложения поступательного движения по траектории вместе с некоторой начальной точкой (точкой M в рассмотренном приме-

76

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