Слайд 1Эйлер, Ляпунов, Навье и Стокс
Слайд 2При обтекании цилиндра жидкостью, его поверхность должна быть линией тока,
а =const.
Для задания такой функции тока можно подобрать систему
источников и стоков.
Слайд 3х
цилиндр
Р
О
Q
R
источники
сток
Показать, что поверхность цилиндра является линией тока
Слайд 5Разложение комплексного
потенциала
в ряд
Слайд 6Всякая конечная, непрерывная и однозначная функция f(z), первая производная которой
конечна для всех точек области между двумя концентрическими окружностями, описанными
около начала координат, может быть разложена в ряд:
Слайд 7Если условия выполнены внутри круга, то
Если условия выполнены вне круга,
то
Если условия выполнены во всей плоскости ху
, вводя полярные координаты
, представляя комплексные постоянные
получим:
Слайд 10Если значение или
задано на концентрических окружностях, то это значение
может быть разложено в ряд Фурье по косинусам и синусам
кратного угла . Эти ряды должны быть эквивалентны рядам полученным выше. Приравнивая в отдельности коэффициенты при sin(n) и при cos(n) , можно получить уравнения для определения Pn, Qn, Rn, Sn.
Слайд 11U
2a
Цилиндр движется в неограниченной жидкости, которая на бесконечности находится в
покое
U
х
у
Предполагаем, что движение возникает из состояния покоя. Тогда оно остается
безвихревым, а потенциал скорости будет однозначным вне цилиндра
Слайд 12U
Нормальная составляющая скорости в жидкости на границе с цилиндром равна
для r=a
r
Слайд 13Дифференцируя потенциал скорости по r и приравнивая нормальные составляющие скорости
на границе цилиндра r=a, получаем:
Слайд 14Остальные коэффициенты равны нулю
Найти функцию тока
Слайд 15Так как циркуляция скорости
то - однозначна и может быть
представлена в виде ряда вместе с вне цилиндра
Слайд 18Полное решение задачи:
=const вдоль линии тока
а
Слайд 19у
Течение, возникающее при движении цилиндра с постоянной скоростью в неподвижной
жидкости, подобно течению, возникающему в жидкости при наличии дублета в
начале координат.
Слайд 20Задача
Найти
кинетическую
энергию
потока
Слайд 23Кинетическая энергия системы (цилиндр+жидкость),
отнесенная к единице длины цилиндра, равна
сумме энергии цилиндра ( 1/2 МU2)
и энергии жидкости (1/2 М выт.жидкU2)
Слайд 25Если при прямолинейном движении цилиндра на него действует внешняя сила
Х , отнесенная к единице длины, то уравнение изменения энергии
будет:
Слайд 26Со стороны жидкости на цилиндр действует сила
Она исчезает, если скорость
не зависит от времени.
Слайд 27Определим эту силу, используя интеграл движения, записанный через потенциал скорости
для нестационарного движения жидкости:
q - результирующая скорость:
Слайд 30Для того, чтобы найти силу, действующую на цилиндр со стороны
окружающей жидкости, надо проинтегрировать давление по границе цилиндра. Проекция на
ось х:
Слайд 31Если цилиндр движется в идеальной жидкости с постоянной во времени
скоростью, то со стороны жидкости на него не действует силы.
Слайд 32Обтекание
цилиндра
плоскопараллельным
потоком
Слайд 33-U
2a
Сообщим жидкости и цилиндру скорость -U. Тогда жидкость будет обтекать
неподвижный цилиндр
Пусть цилиндр движется со скоростью U, как показано на
рисунке.
Запишем комплексный потенциал для плоскопраллельного движения жидкости.
Слайд 34Потенциал скорости и функция тока для плоскопараллельного течения n=1
:
х
у
r
Слайд 36Является ли граница цилиндра линией тока?
Каковы линии тока?
Определим потенциал скорости
и функцию тока на границе цилиндра r=a
Слайд 38Ползущее течение
Обтекание круга в лотке Хил-Шоу
Слайд 39Подкраска позволяет увидеть линии тока в воде, текущей со скоростью
1 мм/с между двумя стеклянными пластинками, отстоящими на расстояние одного
миллиметра. Линии тока абсолютно симметричны, такая картина должна наблюдаться в идеальной жидкости (без вязкости)
Слайд 40На первый взгляд представляется парадоксальным, что наилучший способ получения безотрывной
картины плоского потенциального обтекания цилиндра характерного для идеальной жидкости, состоит
в переходе к противоположному крайнему случаю ползущего течения в узком зазоре, для которого влияние сил вязкости является определяющим. Фото D. Н. Регеgrine
Слайд 41Как
обтекает
цилиндр
реальная жидкость?
Слайд 42При этом числе Рейнольдса картина линий тока, очевидно, уже потеряла
ту симметрию передней и задней частей, которая имела место в
ползущем течении.
Обтекание кругового цилиндра при Re=1,54.
Слайд 43Однако поток сзади все же еще не оторвался. Отрыв начинается
примерно при Re = 5, хотя значение числа Рейнольдса начала
отрыва точно неизвестно. Линии тока визуализированы с помощью алюминиевого порошка в воде. Фото Sadatoshi Taneda
Слайд 44Обтекание кругового цилиндра при Re=9,6.
Произошел отрыв, и образовалась пара рециркуляционных
вихрей.
Слайд 45Цилиндр движется в бассейне с водой, содержащей алюминиевый порошок, и
подсвечивается световым ножом под свободной поверхностью. Экстраполяция результатов подобных экспериментов
на случай неограниченного потока указывает на возможность отрыва при Re=4-5, тогда как большинство численных расчетов дает Re=5-7. Фото Sadatoshi Taneda
Слайд 46
По мере увеличения скорости неподвижные вихри начинают вытягиваться в направлении
потока. Их длина линейно растет с ростом числа Рейнольдса, пока
значение Re не превысит 40
Обтекание кругового цилиндра при Re=13,1.
Слайд 47Расстояние вниз по потоку до центров вихрей также линейно возрастает
с ростом числа Рейнольдса.
Обтекание кругового цилиндра при Re=26.
Слайд 48При таком числе Рейнольдса рециркуляционный след простирается на целый диаметр
вниз по потоку, однако он полностью сохраняет свою стационарность, как
и в случае кругового цилиндра. Визуализация осуществляется тонким слоем сгущенного молока на шаре; молоко постепенно растворяется и уносится потоком воды. [Taneda, 1956b]
Обтекание шара при Re=104
Слайд 49Подкраска обнаруживает ламинарный пограничный слой, отрывающийся перед экватором, причем этот
слой остается ламинарным на длине, почти равной радиусу. Затем слой
становится неустойчивым и быстро превращается в турбулентный. Фото ONERA. [Werle, 1980]
Мгновенная картина потока при обтекании шара при Re=15 000
Слайд 51Придонный слой конечной толщины может быть остановлен этими силами
Силы действующие
на элементарный объем воды вблизи дна в замедляющемся потоке
х
Ffriction
u1p1
u2p2
y
Свободная поверхность
Слайд 52u=0
u=0
umax
umax
Fтр
Под действием силы трения и градиента давления происходит периодическая остановка
поверхностного слоя и образование вихрей
Слайд 53Вода обтекает цилиндр диаметром 1 см co скоростью 1,4 см/с.
Визуализация движения осуществляется так: частицы метятся белым коллоидным дымом, создаваемым
электролитическим способом и освещаются световым ножом. Видно, что по мере продвижения вниз по потоку на несколько диаметров ширина вихревой пелены возрастает. Фото Sadatoshi Taneda
Вихревая дорожка Кармана за круговым цилиндром при Re = 140.
Слайд 54Камера движется здесь вместе с вихрями, а не с цилиндром.
Структура линий тока весьма напоминает картину невязкого течения, рассчитанную Карманом.
Визуализация потока осуществляется с помощью частиц, плавающих на воде. Фото R. Wille, снимок взят из статьи [Werle, 1973]. Воспроизведено с соответствующего разрешения из Annual Review of Fluid Mechanics, Vol. 5, © 1973 by Annual Reviews Inc
Вихри Кармана в абсолютном движении
Слайд 55u=0
u=0
umax
umax
Fтр
u
x
y
В вязкой жидкости
Р1
Р2
Р1> Р2
Слайд 56Зависимость скорости потока от времени (внутри придонного d слоя) между
моментами вылета вихрей
T
Слайд 571. Почему не образуются вихри за цилиндром?
Подкраска позволяет увидеть линии
тока в воде, текущей со скоростью 1 мм/с между двумя
стеклянными пластинками, отстоящими на расстояние одного миллиметра. Линии тока абсолютно симметричны, такая картина должна наблюдаться в идеальной жидкости (без вязкости)
Слайд 58Над поверхностью воды существует горизонтальный поток воздуха, скорость потока уменьшается
в направлении движения. Как возникают волны на поверхности воды?
Слайд 59Ветер, Vвет>C
C
Р1> Р2
Р1
Р2
Усиление волн ветром
Слайд 61Для стационарного движения при отсутствии внешних сил
Для стационарного движения -
( ) выполняется вдоль данной
линии тока
Для безвихревого движения ( ) выплняется везде
Слайд 62Определить скорость на границе и компоненты силы, действующей со стороны
жидкости на цилиндр
u=0
u=0
umax
umax
Слайд 64Куда направлена сила, действующая на цилиндр в потоке жидкости со
сдвигом скорости?
u
u=0
u=0
u1
u2
x
y
Слайд 65u=0
u=0
umax
umax
Куда направлена сила, действующая на вращающийся цилиндр в однородном потоке
жидкости?
Слайд 66Для того, чтобы ряд, представляющий комплексный потенциал, соответствовал случаю произвольного
безвихревого движения между двумя концентрическими окружностями, к ряду
необходимо прибавить еще
слагаемое, описывающее влияние вихревых нитей
Слайд 67Если A=P+iQ, то соответствующие выражения для потенциала скорости и функции
тока будут:
Циклическая константа функции тока 2Р означает поток через внутреннюю
или внешнюю окружности.
Циклическая константа потенциала скорости 2Q означает циркуляцию по некоторой замкнутой кривой, заключающей начало координат.
Слайд 68Изучить
обтекание
цилиндра
потоком с циркуляцией
Задача
Слайд 69Цилиндр двигается с постоянной скоростью и вращается с циклической постоянной
.
Слайд 70Линии тока для
В этом случае точка, в которой скорость равна
нулю находится в жидкости
Слайд 72Для r = a
Сила действует на цилиндр вдоль вертикальной
оси
не зависит от радиуса
Слайд 73Всякое непрерывное движение жидкости, наполняющей неограниченное пространство и покоящейся в
бесконечности, можно рассматривать, как движение, вызванное соответствующим распределением источников и
вихрей с конечной плотностью.
Слайд 74Всякое непрерывное безвихревое циклическое и нециклическое движение несжимаемой жидкости, наполняющей
произвольную область, может рассматриваться как движение, вызванное некоторым распределением вихрей
по ограничивающей поверхности, которая отделяет область от остального неограниченного пространства. В случае области, простирающейся в бесконечность, это распределение относится к конечной части ограничивающей поверхности при условии, что жидкость покоится в бесконечности.