Слайд 1Эйлер, Ляпунов, Навье и Стокс
Слайд 2Поле скорости по
заданному полю
вихрей и
расхождения скорости
Слайд 3Компоненты ротора скорости:
Вихревые линии - линии, направление которых совпадает всюду
с мгновенной осью вращения жидкости.
Слайд 4Дифференциальное уравнение вихревых линий
Если через каждую точку малой замкнутой кривой
провести соответствующую вихревую линию, то получим трубку, которая называется вихревой
трубкой. Жидкость внутри трубки образует вихревую нить или просто вихрь.
Слайд 5 Задача
Заданы распределения вихря и дивергенции скорости в любой точке
жидкости, нормальная составляющая скорости на поверхности, ограничивающей данный объем жидкости
Слайд 6Предположения
Жидкость заполняет все пространство,
находится в покое на бесконечности.
Заданы вихрь
скорости Ω расхождение (дивергенция) скорости , равные 0 вне объема
.
Область может быть разложена на конечное число частей, в которых и Ω равномерно непрерывны, так же как и их частные производные.
Слайд 75. На поверхностях разрыва нормальная составляющая Ω остается непрерывной
Искомую скорость
u будем рассматривать как сумму двух слагаемых: одно определяется расхождением
скорости, и ее вихрь равен нулю, а второе – ротором скорости, а ее дивергенция равна нулю.
Слайд 8Скорость задается уравнениями:
существует
Слайд 9Для определения u1 получаем уравнение Пуассона
Предположим, что функция Θ равна
0 всюду, кроме очень малой окрестности τ0 начала координат, причем
Слайд 10По теореме Гаусса
Т.е. объемный интеграл от расхождения вектора скорости равен
потоку вектора скорости через поверхность S0, ограничивающую объем τ0. Поток
должен равняться единице в силу предположения
Слайд 11Пусть τ→0. Получаем картину течения от точечного источника в начале
координат интенсивности 1. В силу симметрии потенциал скорости φ –
функция только r
φ всюду, кроме начала координат, удовлетворяет уравнению Лапласа
Слайд 12φ – функция только r, в сферических координатах нет зависимости
от широты и долготы.
Интегрируем:
Слайд 13Произвольная постоянная определяется из условия, что поток скорости через произвольную
сферу с центром в начале координат равен 1. На сфере
нормальная составляющая скорости имеет постоянное значение:
а площадь поверхности сферы равна 4πr2
Слайд 14Получаем:
Произвольная постоянная не влияет на величину скорости течения и может
быть отброшена
Слайд 15Предположим, что интенсивность источника имеет значение q
В этом случае
Слайд 16Если задано распределение источников Θ(ξ, η, ζ)
Где r расстояние от
точки N(x,y,z), где ищем поле скорости до точки М(ξ,η,ζ), где
расположен источник. Интегрировать надо по (ξ,η,ζ).
Слайд 18Определим вектор u2 для системы
Удовлетворим второму уравнению, если положим
где А
– векторный потенциал. Тогда
Можно доказать тождество:
Слайд 19Получаем уравнение:
(Можно считать не нарушая общности)
Слайд 21Одна вихревая нить
Несжимаемая
жидкость, покоящаяся
на бесконечности
Слайд 22Применяя к вихревой трубке свойство
учитывая, что боковые поверхности трубки –
есть вихревые линии, т.е. параллельны ротору скорости, получаем для суммарного
потока вихря:
ds
Слайд 27ds
x,y,z
Составляющие единичного вектора k
Составляющие единичного вектора m
Слайд 28Вклад в величину скорости в точке (x,y,z) от элемента вихревой
трубки ds определяется выражением:
Слайд 29Электродинамика:
Сила, действующая на магнитный полюс в точке (x,y,z) от элемента
проводника ds, по которому течет ток (Био и Савара)
Слайд 30Прямолинейные
вихри
плоское движение
несжимаемая жидкость
Слайд 31Пусть движение происходит в плоскости х, у. Вихревые линии являются
прямыми, параллельными оси z. ξ, η, ζ – координаты точек
вихревой трубки, ds – элемент дуги трубки
Слайд 32Скорость жидкости в точке (х, у) определяется:
Считая ξ, η, (х,
у, z) постоянными, интегрируем
Слайд 33Достаточно рассматривать движение на плоскости 0xy, причем вместо вихревой нити
точку пересечения ее с плоскостью 0xy. Будем называть ее точечным
вихрем. Под влиянием такого вихря частицы жидкости двигаются по окружностям, центром которых является вихрь. Положительным соответствует движение против часовой стрелки.
Вследствие симметрии движения центр вихря не будет смещаться.
Слайд 34Найти комплексный потенциал для точечного вихря.
Слайд 36 2 вихря
Две параллельные прямые вихревые нити в точках z1
и z2. Показать, что нити всегда сохраняют одинаковое расстояние друг
от друга и вращаются с постоянной угловой скоростью вокруг общего центра С. Пусть циркуляция одной нити равна 1 , а второй 2 .
r
Запишите комплексный потенциал для 2 нитей.
Слайд 37Комплексная сопряженная скорость
Слайд 38Скорость первого вихря в точке z1 (сам на себя не
действует)
Скорость второго вихря в точке z2
Отделяем мнимые и действительные части.
Обозначим расстояние между вихрями
Слайд 40Получаем интеграл движения центра инерции системы двух вихрей
Точка С (центр
инерции) с координатами
Остается неподвижной во все время движения
Слайд 41Вычитаем из первого третье уравнение, из второго – четвертое.
Слайд 42Интегрируем и получаем:
Расстояние между вихрями не меняется в процессе перемещения
Слайд 43A
B
x
С
Два вихря, вращающиеся в разных направлениях, находятся в начальный момент
времени на вещественной оси.
Найти скорости вихрей, расстояние до центра
системы, угловую скорость вращения вихрей.
Слайд 44Вихревая нить с координатами
и циркуляцией сообщает жидкости в точке
(х, у) скорость, компоненты которой равны
х
(х, у)
r
Слайд 46Угловая скорость вращения вихрей вокруг центра С
Где будет точка С,
если вихри вращаются в одном направлении?
Слайд 47x
r
Два вихря, вращающиеся в разных направлениях, находятся в начальный момент
времени на вещественной оси.
Найти скорости вихрей, расстояние до центра
системы.
Слайд 48Скорость первого вихря в точке z1 (сам на себя не
действует)
Скорость второго вихря в точке z2
Отделяем мнимые и действительные части.
Обозначим расстояние между вихрями
Слайд 49Отделяем вещественную часть от мнимой :
Вихри перемещаются с постоянной скоростью,
перпендикулярно прямой, соединяющей вихри в положительном направлении оси у
Слайд 50Пример 1
Одна вихревая нить в точке (х,у), циркуляция скорости внутри
бесконечно малого сечения имеет постоянное значение.
Найти центр системы
Слайд 51х
у
r
Центр одиночного вихря не смещается во времени.
Слайд 52Пример 2
Внутри круга радиуса а жидкость вращается как твердое тело
с угловой скоростью . Определить скорость вне круга.
Слайд 53х
у
r
На границе вихря скорость равна a.
а
Слайд 54Как будут двигаться 2 вихря радиуса а, если они имеют
циркуляцию разного знака, но одинаковую по модулю? Вихри вращаются как
твердое тело.
Пример 3
А
а
В
а
А
В
Слайд 55Движение двух вихрей с противоположным направлением вращения и
С
А
АС=
V1=V2=
Вихри двигаются по
прямой с одинаковой скоростью
а
В
Слайд 57а
Куда
будет двигаться
вихрь?
Пример 4
Слайд 58Найти скорость перемещения вихря у твердой стенки
Слайд 59Так как скорость жидкости во всех точках плоскости симметрии направлена
по касательной, то можно предположить, что эта плоскость образует твердую
границу для жидкости. Таким образом систему «вихрь у твердой границы» можно смоделировать системой, представляющей собой пару вихрей. Вихрь у твердой границы будет перемещаться с постоянной скоростью
а
2а
Слайд 60Пример 5
В точке с координатами х,у находится прямая вихревая нить.
Жидкость ограничена твердыми стенками, образующими прямой угол вдоль осей координат.
Найти траекторию вихря.
Слайд 62Так как циркуляция скорости вихря постоянна, то течение стационарно. В
этом случае траектория вихря совпадает с линией тока, проходящей через
точку (x,y)
Для того, чтобы описать влияние стенки, расположенной вдоль мнимой оси, поместим вихрь той же интенсивности, но с противоположным знаком циркуляции в точку (-x,y ) (отразим вихрь)
Слайд 64Для того, чтобы описать влияние стенки, расположенной вдоль действительной оси,
поместим вихри той же интенсивности, но с противоположными знаками циркуляции
в точки (-x,-y ) и (x,-y )
Слайд 65x,y
-x,y
x,-y
-x,-y
Надо найти потенциал в точке (x,y)
Слайд 67x,y
-x,y
x,-y
-x,-y
Траектория вихря в углу, образованном твердыми стенками
Слайд 78расположенных в точках z1,…, zn и имеющих интенсивности 1, …,
n
Комплексный потенциал
Комплексная скорость
Записать скорость вихря в точке zp (k=p)
n
вихрей
Слайд 79Уравнение движения вихря в точке zp (k=p)
Умножая на p и
суммируя по переменной p от 1 до n, получаем в
правой части
Так как слагаемому с (zp-zk) соответствует слагаемое с (zk-zp)
Слайд 80Следовательно
Получаем
Отделяем мнимую и действительную части
, то получаем интегралы
движения центров инерции (следствие 1)
Уравнение движения вихря в точке zp (k=p)
Умножая на p zp и суммируя по переменной p от 1 до n, получаем
Отделить мнимую и действительную части
Слайд 83Сумма моментов количеств движения масс p относительно начала координат не
меняется со временем
(1)
(2)
(2)
Слайд 85Сумма моментов инерции масс p относительно начала координат не меняется
со временем
Слайд 86Умножая выражения скорости на
и суммируя по переменной p от 1
до n, получаем еще один интеграл
Следствие 3
Слайд 87А
а
В
Такие вихри называются «парой вихрей». Они являются плоской аналогией вихревому
кольцу и обладает многими свойствами последнего.
Слайд 89круговая вихревая нить
радиуса а лежит в плоскости х,у
Слайд 90Координаты точки М на вихревой нити в декартовой системе (x,
y, z)
Расстояние от точки М на вихревой нити до точки
N(x,y,z)
Слайд 91(1)
Запишем компоненты скорости через векторный потенциал
Слайд 93Для векторного потенциала в точке x,y,z
Слайд 94Для плоскости x,y при θ=0
Вследствие симметрии задачи это справедливо для
всех углов θ
Слайд 96Подставляем выражения векторного потенциала в (1) и находим компоненты скорости
Слайд 97ОПРЕДЕЛИТЬ объем жидкости, протекающей через круг с радиусом ρ, лежащий
в плоскости перпендикулярной оси z с центром на оси z
, через переменную А
Слайд 98Объем жидкости, протекающей через круг с радиусом ρ, лежащим в
плоскости перпендикулярной оси z с центром на оси z
Слайд 99ОПРЕДЕЛИТЬ
скорость перемещения круговой вихревой нити
Слайд 101Вихревой слой
0
а
S
S1
x
y
u1
u
Найти ротор скорости
Слайд 102Внутри слоя = const
Разрыв скорости может быть смоделирован цепочкой
вихревых нитей
Слайд 103Вихревой цилиндрический слой
x
r
Интенсивность вихрей на элементе дуги rd
Записать комплексный потенциал
и комплексную скорость для точки z
Слайд 104Комплексный потенциал цилиндрического слоя в точке z
Комплексная скорость
Делим подинтегральное выражение
на z, умножаем на
Слайд 105z рассматривается как постоянная при интегрировании
Слайд 106Если z внутри окружности, т.е.
При изменении от 0 до
2 аргумент разности
меняется на 2, так как вектор обходит начало координат. Тогда
z
r
0
Слайд 107Если z вне окружности, т.е.
При изменении от 0 до
2 конец вектора
обходит замкнутый путь не содержащий начала координат. Тогда
z
r
0
Записать комплексную скорость внутри и вне цилиндра
Слайд 108Комплексная скорость
движения нет внутри цилиндра
вне цилиндра движение описывается воздействием вихревой
нити интенсивности 2r, расположенной в начале координат
Слайд 110U
Записать функцию тока для вихревой пары в потоке жидкости, имеющего
постоянную скоростью U. Скорость потока равна по модулю скорости пары,
но противоположна по направлению.
№1
Слайд 112Функция тока для течения, имеющего постоянную скорость вдоль оси ординат
Слайд 113Относительные линии тока вихревой пары
О
на оси у и на овале
О (жирная линия).
Слайд 114Жидкость внутри овала О движется вместе с вихрями.
Жидкость вне
овала О обтекает этот овал как твердый цилиндр.
Полуоси овала О
приблизительно имеют длину 2.09а и 1.73а.
Слайд 115Найти в ортогональной системе координат уравнения линий тока для
1)
случая двух вихрей одинаковой интенсивности
2) пары вихрей
№2
Слайд 119При обтекании цилиндра жидкостью, его поверхность должна быть линией тока,
а =const.
Пусть в центре цилиндра имеется сток, а на
оси 0х лежат два источника: один внутри цилиндра, другой вне цилиндра, Q – инверсия точки Р, мощность одинакова. Показать, что
х
Р
О
Q
R
источники
сток
поверхность цилиндра является линией тока
Слайд 120При обтекании цилиндра жидкостью, его поверхность должна быть линией тока,
а =const.
Пусть в центре цилиндра имеется сток, а на
оси 0х лежат два источника: один внутри цилиндра, другой вне цилиндра, Q – инверсия точки Р, мощность одинакова. Показать, что
х
Р
О
Q
R
источники
сток
поверхность цилиндра является линией тока
Слайд 122х
цилиндр
Р
О
Q
R
источники
сток
Так как,
То треугольники ORQ и ORP подобны
Углы ORQ и RPO
равны
Слайд 123Вихревой цилиндрический слой
x
r
Интенсивность вихрей на элементе дуги rd
Найти координаты центра
системы
№4
Слайд 124Интенсивность твердотельно вращающихся вихрей такова, что тангенциальные составляющие скорости в
точках соприкосновения равны
Найти координаты центра системы
№5
Слайд 126U
Записать функцию тока для вихревой пары в потоке жидкости, имеющего
постоянную скоростью U. Скорость потока равна по модулю скорости пары,
но противоположна по направлению.
№1
Слайд 127Найти в ортогональной системе координат уравнения линий тока для
1)
случая двух вихрей одинаковой интенсивности
2) пары вихрей
№2