1
w= +i 2
Условия Коши-Римана для и iy х - это условия аналитичности функции w= +i комплексного аргумента z=x+iy. Это значит, что эта функция имеет определенную производную u v 3
D = x(x 0, y 0)D x + y(x 0, y 0)D y +x (x, y); D = x(x 0, y 0)D x+ y(x 0, y 0)D y+h (x, y); V(x, y)=x (x, y)+ih (x, y) 4
5
iy х - комплексная скорость Функция w= +i называется комплексным потенциалом, однозначно определяет плоское течение 6
Задан комплексный потенциал Определить линии тока, линии равного потенциала, поле скорости 7
8
y Линия тока х 9
10
v a u b 11
iy a х b 12
вдоль линии тока Линии тока определены как: iy = a co ns t х b 13
14
15
y z х Линии х=0 и у=0 являются частями одной и той же линии тока =0. 16
y на линии равного потенциала х 17
1. Уравнение Фридмана и Гельмгольца 2. Можно ли вызвать движение неподвижной идеальной жидкости, если приложить постоянное давление по всему объему? 3. Написать уравнение неразрывности через потенциал скорости, раскрыв полную производную по времени 4. Исследовать течение жидкости, которое описывается комплексным потенциалом 18
19
Запишем z в полярных координатах 20
21
y х 22
y х 23
Значение скорости становится бесконечно большим в начале координат. Исключаем эту точку, окружив произвольной замкнутой линией. Это особая точка для комплексного потенциала и для комплексной скорости. 24
25
26
Линии = 0 и представляют собой части одной и той же линии тока. Если положить , то получим безвихревое движение, при котором граница состоит из двух твердых наклонных стенок 27
Решение: На линии тока =const 28
Компоненты скорости в направлении r и в перпендикулярном направлении равны: 29
30
31
= co t ns На линии равного потенциала 32
В начале координат особая точка для скорости и комплексного потенциала. Исключаем эту точку, окружив произвольной замкнутой линией. 33
Линии равной скорости концентрические окружности. В начале координат - особая точка. Поле скорости непрерывно и конечно во всех точках кроме одной. Начало должно быть исключено с помощью описанной вокруг него замкнутой кривой. 34
Гидродинамическое истолкование потенциала = co t ns Источник из начала координат. Определить поток массы в единицу времени. 35
Поток массы (мощность источника), проходящий через окружность радиуса r в единицу времени, определяется произведением радиальной скорости на длину окружности Мощность источника 2 не зависит от радиуса (сохранение массы). Если <0, то определяет сток. 36
Если в точках плоскости z=a 1, a 2, …an находятся источники или стоки мощностью 2 1, 2 2, …, 2 n k – модуль, k – аргумент числа (z – ak) 37
Линии тока r = const концентрические окружности. 38
Движение во всех точках безвихревое кроме начала кординат, где. Эту точку назовем вихревой. Потенциал скорости - неоднозначная функция координат на плоскости ОХУ. При обходе вокруг начала координат величина потенциала меняется на. Скорость: Циркуляция скорости: Интенсивность вихревой точки 39
Движение безвихревое везде, кроме начала координат Вихревая точка, расположенная в начале координат z=0 создает плоское движение, определяемое комплексным потенциалом 40
1. Исследовать течение жидкости, которое описывается комплексным потенциалом 2. a, b - вещественные 41
Река имеет участок, где оба берега представляют собой сегмент вложенных окружностей с единым центром. Показать, что у берега А скорость течения больше, а уровень ниже, чем у берега В. Считать движение установившимся и безвихревым. z Уравнение неразрывности A B r 42
Проекции скорости на оси координат Так как движение безвихревое rotu=0 Из уравнения неразрывности Из уравнения Бернулли 43
z r 1 -а 1 r 2 а 2 44
=const на линии тока 45
z r 1 -а 1 r 2 а 2 46
=const источник -а сток а Источник и сток равной мощности 47
Пусть , т. е. источник и сток расположены в начале координат. Такую пару называют дублетом. 48
=const -а а 49
При обтекании цилиндра жидкостью, его поверхность должна быть линией тока, а =const. Почему? Так как нормальная составляющая скорости равна нулю (нет протекания) Для задания такой функции тока можно подобрать систему источников и стоков. 50
R цилиндр О сток Q Р источники х 51
-3 а -2 а а О а 2 а 3 а Источники одинаковой мощности на равном расстоянии 52
Эти же формулы имеют место для источника, помещенного на одинаковом расстоянии от двух параллельных твердых стенок, имеющих координаты y х 53
54
Всякая конечная, непрерывная и однозначная функция f(z), первая производная которой конечна для всех точек области между двумя концентрическими окружностями, описанными около начала , может быть разложена в ряд: 55
Если условия выполнены внутри круга, то Если условия выполнены вне круга, то Если условия выполнены во всей плоскости ху 56
Полагая координаты постоянные получим: , вводя полярные , представляя комплексные 57
Если значение или задано на концентрических окружностях, то это значение может быть разложено в ряд Фурье по косинусам и синусам кратного . Эти ряды должны быть эквивалентны рядам полученным выше. Приравнивая в отдельности коэффициенты при sin(n ) и при cos(n ) , можно получить уравнения для определения Pn, Qn, Rn, Sn. 58
2 a U у U Так как Цилиндр движется в неограниченной жидкости, которая на бесконечности находится в покое Предполагаем, что движение возникает из состояния покоя. Тогда оно остается х безвихревым, а потенциал скорости будет однозначным вне цилиндра то - однозначна и может быть представлена в виде ряда 59 вместе с вне цилиндра
r U Нормальная составляющая скорости в жидкости на границе с цилиндром равна для r=a Так как Pn=0, Qn =0, то дифференцируя потенциал скорости по r и приравнивая нормальные составляющие скорости на границе цилиндра r=a, получаем: 60
Остальные коэффициенты равны нулю Полное решение задачи: =const вдоль линии тока 61
у 2 a U Течение, возникающее при движении цилиндра с постоянной скоростью в неподвижной жидкости, подобно течению, возникающему в жидкости при наличии дублета в начале координат. 62
63
Кинетическая энергия жидкости: Кинетическая энергия системы (цилиндр+жидкость), отнесенная к единице длины цилиндра, равна сумме энергии цилиндра ( 1/2 МU 2) и энергии жидкости (1/2 М выт. жидк. U 2) 64
65
Если прямолинейном движении цилиндра на него действует внешняя сила Х , отнесенная к единице длины, то уравнение изменения энергии будет: Со стороны жидкости на цилиндр действует сила Она исчезает, если скорость не зависит от времени. 66
Определим эту силу используя интеграл движения, записанный через потенциал скорости для нестационарного движения жидкости: q - результирующая скорость: 67
Для того, чтобы найти силу, действующую на цилиндр со стороны окружающей жидкости, надо проинтегрировать давление по границе цилиндра. Проекция на ось х: 68
Если цилиндр движется в идеальной жидкости с постоянной во времени скоростью, то со стороны жидкости на него не действует силы. 69
70
2 a 2 a Пусть цилиндр движется со скоростью U, как показано на рисунке. U -U Сообщим жидкости и цилиндру скорость -U. Тогда жидкость будет обтекать неподвижный цилиндр Запишем комплексный потенциал для плоскопраллельного движения жидкости. 71
Потенциал скорости и функция тока для плоскопараллельного течения n=1 у : r х 72
2 a 2 a U -U Надо прибавить 73
Определим потенциал скорости и функцию тока на границе цилиндра r=a Является ли граница цилиндра линией тока? Каковы линии тока? 74
Обтекание круга в лотке Хил-Шоу Подкраска позволяет увидеть линии тока в воде, текущей со скоростью 1 мм/с между двумя стеклянными пластинками, отстоящими на расстояние одного миллиметра. Линии тока абсолютно симметричны, такая картина должна наблюдаться в идеальной жидкости (без вязкости) 75
На первый взгляд представляется парадоксальным, что наилучший способ получения безотрывной картины плоского потенциального обтекания цилиндра характерного для идеальной жидкости, состоит в переходе к противоположному крайнему случаю ползущего течения в узком зазоре, для которого влияние сил вязкости является определяющим. Фото D. Н. Регеgrine 76
77
Для стационарного движения при отсутствии внешних сил Для стационарного движения -( ) выполняется вдоль данной линии тока Для безвихревого движения -( ) выплняется везде 78
Какая сила действует на тело в потоке жидкости, если движение безвихревое? Скорость на границе: y umax=2 U u u=0 x umax 0 79
Куда направлена сила, действующая на тело в потоке жидкости со сдвигом скорости? y u 1 u u=0 x u 2 80
Какая сила действует на вращающееся тело в потоке жидкости, если движение безвихревое? y umax u u=0 x umax 81
82
Обтекание кругового цилиндра при Re=1, 54. При этом числе Рейнольдса картина линий тока, очевидно, уже потеряла ту симметрию передней и задней частей, которая имела место в ползущем 83 течении.
Однако поток сзади все же еще не оторвался. Отрыв начинается примерно при Re = 5, хотя значение числа Рейнольдса начала отрыва точно неизвестно. Линии тока визуализированы с помощью алюминиевого порошка в воде. Фото Sadatoshi Taneda 84
Обтекание кругового цилиндра при Re=9, 6. Произошел отрыв, и рециркуляционных вихрей. образовалась пара 85
Цилиндр движется в бассейне с водой, содержащей алюминиевый порошок, и подсвечивается световым ножом под свободной поверхностью. Экстраполяция результатов подобных экспериментов на случай неограниченного потока указывает на возможность отрыва при Re=4 -5, тогда как большинство численных расчетов дает Re=5 -7. Фото Sadatoshi Taneda 86
Обтекание кругового цилиндра при Re=13, 1. По мере увеличения скорости неподвижные вихри начинают вытягиваться в направлении потока. Их длина линейно растет с ростом числа Рейнольдса, пока значение 87 Re не превысит 40
Обтекание кругового цилиндра при Re=26. Расстояние вниз по потоку до центров вихрей также линейно возрастает с ростом числа 88 Рейнольдса.
Обтекание шара при Re=104 При таком числе Рейнольдса рециркуляционный след простирается на целый диаметр вниз по потоку, однако он полностью сохраняет свою стационарность, как и в случае кругового цилиндра. Визуализация осуществляется тонким слоем сгущенного молока на шаре; молоко постепенно растворяется и уносится 89 потоком воды. [Taneda, 1956 b]
Мгновенная картина потока при обтекании шара при Re=15 000 Подкраска обнаруживает ламинарный пограничный слой, отрывающийся перед экватором, причем этот слой остается ламинарным на длине, почти равной радиусу. Затем слой становится неустойчивым и быстро 90 превращается в турбулентный. Фото ONERA. [Werle, 1980]
91
Силы действующие на элементарный объем воды вблизи дна в замедляющемся потоке Если и p -давление u 1>u 2 y рхн бодная пове Сво u 1 p 1 ость u u 2 p 2 х Придонный слой конечной толщины может быть остановлен этими силами Ffriction 92
umax Fтр u=0 umax Под действием силы трения и градиента давления происходит периодическая остановка u=0 поверхностного слоя и образование вихрей 93
Вихревая дорожка Кармана за круговым цилиндром при Re = 140. Вода обтекает цилиндр диаметром 1 см co скоростью 1, 4 см/с. Визуализация движения осуществляется так: частицы метятся белым коллоидным дымом, создаваемым электролитическим способом и освещаются световым ножом. Видно, что по мере продвижения вниз по потоку на несколько диаметров ширина вихревой пелены возрастает. Фото Sadatoshi 94 Taneda
Вихри Кармана в абсолютном движении Камера движется здесь вместе с вихрями, а не с цилиндром. Структура линий тока весьма напоминает картину невязкого течения, рассчитанную Карманом. Визуализация потока осуществляется с помощью частиц, плавающих на воде. Фото R. Wille, снимок взят из статьи [Werle, 1973]. Воспроизведено с соответствующего разрешения из Annual Review of Fluid Mechanics, Vol. 5, © 1973 by Annual Reviews Inc 95
Зависимость скорости потока от времени (внутри придонного d слоя) между моментами вылета вихрей T 96
y В вязкой жидкости u y u umax x x Fтр Р 1 Р 2 u=0 Р 1> Р 2 umax 97
Ветер, Vвет>C Р 1 C Р 2 Р 1> Р 2 Усиление волн ветром 98
99
Для того, чтобы ряд, представляющий комплексный потенциал, соответствовал случаю произвольного безвихревого движения между двумя концентрическими окружностями, к ряду Необходимо прибавить еще член 100
Если A=P+i. Q, то соответствующие выражения для потенциала скорости и функции тока будут: Циклическая константа функции тока 2 Р означает поток через внутреннюю или внешнюю окружности. Циклическая константа потенциала скорости 2 Q означает циркуляцию по некоторой замкнутой кривой, заключающей начало. 101
2 a U Цилиндр двигается с постоянной скоростью и вращается с циклической постоянной . Граничные условия будут выполнены, если предположить, что На одной стороне цилиндра скорость больше, а на другой - меньше. Возникает разность давлений. Для поддержания горизонтального движения надо приложить силу в вертикальном 102 направлении.
Линии тока для В этом случае точка, в которой скорость равна нулю находится в жидкости 103
104
2 a -U Обтекание вращающегося цилиндра Для r = a Сила действует на цилиндр вдоль вертикальной оси не зависит от радиуса 105
Всякое непрерывное движение жидкости, наполняющей неограниченное пространство и покоящейся в бесконечности, можно рассматривать, как движение, вызванное соответствующим распределением источников и вихрей с конечной плотностью. 106
Всякое непрерывное безвихревое циклическое и нециклическое движение несжимаемой жидкости, наполняющей произвольную область, может рассматриваться как движение, вызванное некоторым распределением вихрей по ограничивающей поверхности, которая отделяет область от остального неограниченного пространства. В случае области, простирающейся в бесконечность, это распределение относится к конечной части ограничивающей поверхности при условии, что жидкость покоится в бесконечности. 107