Доказательство этой зависимости проведём на основе гидромеханических представлений

Вектор-градиент обозначается как grad r и вычисляется по формуле:

,

где - единичные векторы по направлению и вдоль координатных осей.

Проекция вектора grad r на некоторое направление определяет изменение плотностей в этом направлении:

где q - угол между направлениями и ; Cos a_I – направляющие косинусы вектора .

Наибольшее изменение плотности происходит в направлении, нормальном к поверхности (3.2.10).

3.3. ПОТОК ГИДРОМЕХАНИЧЕСКОЙ ХАРАКТЕРИСТИКИ ЧЕРЕЗ ПОВЕРХНОСТЬ

Законы механики формулируются для выделенных механических систем, или совокупностей физических тел. Для сплошной среды это жидкий объём, т.е. выделенный движущийся объём жидкости (текучего тела), сохраняющий при своём движении все составляющие его части (жидкие частицы). Это понятие соответствует лагранжеву методу описания движения текучих тел.

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

Рис.3.7. Поток скорости сквозь контрольную поверхность

Зафиксируем в пространстве, занятом движущейся жидкостью, поверхность А и выделим на этой поверхности около точки с координатами элементарную площадку dA(рис.3.7). Скорость жидкости в этой точке равна n - единичный вектор нормали к поверхности в этой же точке. Нормальная к поверхности составляющая скорости будет при этом равной . Объём жидкости, протекающей в единицу времени через площадку dA, равен .

В элементарном объёме dQ содержится dQ_B гидромеханической характеристики В, которая проносится жидкостью за единицу времени через площадку dA:

(3.3.1)

Поток Q_B гидромеханической характеристики В через контрольную поверхность А (количество характеристики, проносимое жидкостью за единицу времени через поверхность А) составляет

(3.3.2)

Поток гидромеханической характеристики В через контрольную поверхность единичной площади (подынтегральное выражение b ) называется плотностью потока гидромеханической характеристики.

Если принять , то из (3.3.1) следует, что при этом гидромеханической характеристикой является объём жидкости. Объём жидкости, протекающий через контрольную поверхность за единицу времени, или поток объёма жидкости называется расходом (объёмным расходом) Q:

. (3.3.3)

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

 

. (3.3.4)

Если положить , то получим поток Q_k кинетической энергии через контрольную поверхность:

 

. (3.3.5)

Положив , получим поток количества движения Q_I:

 

. (3.3.6)

 

В интегралах (3.3.3.)-(3.3.6) выражения представляют собой плотности потоков объёма, массы, кинетической энергии, количества движения, соответственно. Если в поле скорости u (или любой другой векторной величины) мысленно провести некоторую поверхность S и в каждой точке её задать нормаль , то для определения объёма жидкости, протекающей за единицу времени сквозь поверхность S, необходимо вычислить интеграл:

 

.

 

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

.

В декартовой системе координат дивергенция скорости вычисляется по формуле:

 

.

 

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

 

(3.3.7)

Это равенство называется формулой Гаусса.

3.4. ГИДРОМЕХАНИЧЕСКАЯ ИНТЕРПРЕТАЦИЯ ТЕОРЕМЫ ОСТРОГРАДСКОГО - ГАУССА

 

Зафиксируем неподвижную в пространстве контрольную поверхность А, ограничивающую контрольный объём V. Сквозь эту поверхность протекает жидкость со скоростью . Выделим на ней элементарную площадку dА. Единичный вектор нормали к площадке . Если воспользоваться ортами i, j, k, то

 

.

 

Обозначим модуль скорости ; по определению . Скалярное произведение двух векторов можно выразить через их проекции:

 

, (3.4.1)

а также через модули векторов и угол между ними,

 

(3.4.2)

где u_n - нормальная к поверхности dА составляющая скорости.

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

 

. (3.4.3)

 

Используя (3.4.3), запишем объёмный расход жидкости Q через поверхность dА:

 

. (3.4.4)

 

Согласно теореме Остроградского - Гаусса имеем

 

. (3.4.5)

Рис. 3.8. Определение расхода жидкости сквозь поверхность элементарного параллелепипеда

В данном случае - единичные векторы, направленные соответственно по касательной к L и по нормали к поверхности S.

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

.

На основании теоремы Стокса имеет место равенство

.

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

Из шести составляющих тензора, описывающего вращение жидкого объема вокруг полюса, только три отличаются друг от друга абсолютной величиной. Каждая из них определяет значение мгновенной угловой скорости вращения вокруг оси, параллельной одной из координатных осей. Эти угловые скорости w_x,w_y,w_z можно рассматривать как проекции на соответствующие координатные оси вектора w, определяющего угловую скорость вращения элементарного объема жидкости при его перемещении в трехмерном пространстве. Например, как было показано выше,

.

Обратим внимание на то, что буква z не входит в качестве индекса или координаты в правую часть равенства, определяющего w_z. Аналогично записываются выражения для w_z и w_у. Вектор w в матричной форме имеет вид

. (3.5.1)

В векторном анализе вместо w рассматривают вектор 2w, который обозначают rot u и называют вихрем вектора u (или вихрем скорости):

. (3.5.2)

Рис. 3.10. Сдвиговое продольно-однородное течение (течение Куэтта)

Следует обратить внимание, что ранее рассматривалось лишь пере­мещение объема за бесконечно малые интервалы времени Dt. Если рас­сматривать конечные интервалы времени, то, например, после поворота с угловой скоростью w за время Dt элементарный объем жидкости займет другие точки в поле скорости, и за следующий отрезок времени его угло­вая скорость может измениться. В результате даже при наличии во всех точках скорости и rot u ¹ 0 (такие поля называются вихревыми) мо­жет оказаться, что при перемещении на конечные расстояния эле­ментарного объема жидкости поворот его вокруг полюса не будет вос­приниматься как вращение в обычном смысле. Например, при сдвиговом продольно-однородном течении (течении Куэтта), представленном на рис. 3.10 с полем скорости , имеем, что только , а все остальные составляющие тензора grad u равны нулю.

Отметим, что, хотя мгновенная угловая скорость этого объема отлична от нуля,

,

при перемещении объема abсd на большое расстояние его поворот вокруг точки а не воспринимается как вращение.

В том случае, когда все проекции скорости могут быть определены одной функцией j (х₁, х₂, х₃, t) в виде , то есть = grad j, то говорят, что поле скоростей потенциальное, а функция j - потенциал скорости.

Проекция скорости v_l на любое направление l определяется производной dj/d l.

Необходимым и достаточным условием существования потенциальных течений являются равенства ():

.

Следовательно, безвихревое течение жидкости потенциально.

3.6. ПОЛЯ В ГИДРОДИНАМИКЕ

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

В общем случае поле является пространственным (трёхмерным), иногда задачу упрощают, и рассматривают двумерные (плоские) или одномерные поля. В этом случае полагают, что физические величины зависят от одной или двух пространственных координат.

Если физические величины не зависят от времени, то поле называют стационарным, в противном случае - нестационарным.

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

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

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

Основное свойство скалярной функции а(х₁,х₂,х₃) состоит в том, что её численное значение не меняется при преобразовании координат.

Если перейти от старой х₁,х₂,х₃ к новой х^¢₁,х^¢₂,х^¢₃ системе координат, то значения плотности или температуры в фиксированной точке пространства, естественно, не изменяются:

а (х^¢₁,х^¢₂,х^¢₃) = а (х₁,х₂,х₃).

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

Например, поле скоростей жидкости. Вектор в пространстве трёх измерений может быть задан тремя компонентами:

а₁(х₁,х₂,х₃), а₂(х₁,х₂,х₃), а₃(х₁,х₂,х₃),

то есть, тремя функциями от трёх переменных. Это можно записать в виде матрицы-столбца:

а ÜÞ

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

Пусть l _ij - направляющий косинус оси x^¢_j относительно оси x_i (i = 1,2,3; j = 1,2,3). Вычислим проекции того же вектора на новые оси координат:

a^¢₁ = l ₁₁ a₁ + l ₂₁ a₂ + l_{3 1} a₃;

a^¢₂ = l ₂₁ a₁ + l ₂₂ a₂ + l₂₃ a₃;

a^¢₃ = l ₃₁ a₁ + l ₃₂ a₂ + l₃₃ a₃.

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

То есть,сам вектор не меняется в новых координатах, а меняются его компоненты.

Это выражение можно представить в индексной форме записи как сумму:

Или ещё более короткой

При такой записи пользуются двумя правилами:

1. Соглашение о суммировании. По индексу, встречающемуся дважды (немой индекс), производят суммирование от 1 до 3.

2. Соглашение о ранге. Индекс, встречающийся один раз (свободный индекс), пробегает значения от 1 до 3.

Таким образом, уравнение с одним свободным индексом означает запись трёх уравнений.

Помимо скалярных и векторных полей в механике сплошной среды рассматриваются ещё и тензорные поля.

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

Описание происходит в прямолинейных (декартовых) системах координат. Координаты обозначаем х₁, х₂, х₃, единичные векторы по осям - i₁,i₂, i₃.

Предположим, что в результате вращения осей координат как единого целого вокруг начала координат, мы перешли к новой системе координат Ox₁´x₂´x₃´. Обозначим косинус угла между осями x_i и x´_k старой и новой системы a_ik = cos (x´_i^x_k). Для удобства пользования дальнейшими формулами приводим таблицу.

Таблица 3.1

Теперь перейдём к определению тензора. Пусть каждому направлению соответствует вектор (не обязательно коллинеарный n).

Направлениям осей соответствуют векторы , разложение которых опишем подробно:

(3.6.1)

Если векторы для любого направления выражаются лишь через 3 вектора согласно формуле

, (3.6.2)

то множество векторов образует тензор Т.