Вихревая структура закрученных потоков, отрывных течений и следов тема диссертации и автореферата по ВАК РФ 01.02.05, доктор физико-математических наук Куйбин, Павел Анатольевич

• Куйбин, Павел Анатольевич
• доктор физико-математических наук
• 2003

Оглавление диссертации доктор физико-математических наук Куйбин, Павел Анатольевич

ВВЕДЕНИЕ Q £0 / (5О 'О'

ГЛАВА I. ВИНТОВАЯ ВИХРЕВАЯ НИТЬ J

1.1. Уравнения гидродинамики течений с винтовой симметрией

1.2. Гидродинамика течений с винтовыми вихревыми линиями

1.3.Поле скорости, индуцированное винтовой вихревой нитью

ГЛАВА II. РАЗВИТИЕ МОДЕЛЕЙ ВИХРЕВЫХ СТРУКТУР

2.1.Винтовая вихревая пелена

2.2. Колоннообразные винтовые вихри

2.3.Колоннообразный вихрь с винтообразными вихревыми линиями, движуагийея вдоль плоскости

2.4.Винтовой вихрь с ядром конечного размера

2.5. Структура поля завихренности в сферическом вихре Хикса

ГЛАЗА III. МЕТОД ДИСКРЕТНЫХ ВИХРЕВЫХ ЧАСТИЦ

3.1. МДВЧ для плоских течений в ограниченных односвязных областях

3.2. МДВЧ для осесимметричных течений без закрутки

ГЛАВА IV. СТРУКТУРА И НЕУСТОЙЧИВОСТЬ ДВУМЕРНЫХ

СДВИГОВЫХ ТЕЧЕНИЙ

4.1. Формирование разгонного вихря

4.2. Неустойчивость разгонного вихря

4.3. Неустойчивость следа за тонкой пластиной

ГЛАВА V. ОПИСАНИЕ И КЛАССИФИКАЦИЯ ВИХРЕВЫХ

СТРУКТУР В ЗАКРУЧЕННЫХ ПОТОКАХ

5.1. Параметры закрученного потока

5.2. Винтовая симметрия вихревых течений

5.3. Вихри с прямолинейной осью

5.4. Стационарные вихри винтовой формы

5.5. Прецессия вихревого ядра

5.6. О влиянии вихревой структуры течения на горение в закрученном потоке и энергоразделение в вихревой трубке Ранка

5.6.1. Форма пламени в закрученном потоке

5.6.2. Учет структуры течения в эффекте Ранка

ГЛАВА VI. ДВИЖЕНИЕ ВИХРЕЙ В ЦИЛИНДРИЧЕСКИХ ТРУБАХ

6.1. Уравнения движения вихревых частиц в круговой области

6.2. Прецессия прямолинейного вихря в трубе

6.3. Движение винтовых вихрей в трубе

6.3.1. Самоиндуцированное движение винтового вихря

6.3.2. Движение винтового вихря в трубе

6.3.3. Частота прецессии винтового вихря в трубе

6.3.4. Движение системы соосных винтовых вихрей в трубе

Введение диссертации (часть автореферата) на тему «Вихревая структура закрученных потоков, отрывных течений и следов»

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

На современном этапе исследований закрученных потоков стало ясно, что основную роль в них играют трехмерные вихревые структуры винтовой формы (см. обзор [С.В. Алексеенко, В.Л. Окулов, 1996]). В то же время до сих пор в огромном потоке публикаций преобладают работы, игнорирующие новые представления о сложной пространственной структуре сильнозакрученных потоков. Хотя для улучшения эффективности вихревых аппаратов увеличивают закрутку потока (при этом течение становится существенно не осесимметричным), в основной массе публикаций их продолжают рассматривать с позиций устаревших представлений. В создании расчетных методов здесь преобладают упрощенные осесимметричные модели, игнорирующие пространственную структуру реального течения. Многие экспериментальные исследования все еще ограничиваются лишь измерением осредненных характеристик потока, по которым невозможно получить представление о пространственной структуре реального течения. Продолжает удивлять живучесть устаревших (осредненных) подходов, несмотря на работы Faler & Leibovich [1977] и многих других исследователей по неустойчивости закрученных потоков и монографию А. Гупта [1987], определившую задачу изучения потери устойчивости и образующихся сложных вихревых структур как основную проблему дальнейшего развития вихревых технологий. Понятно, что без изучения трехмерной структуры сильно закрученного потока невозможно правильно оценить новые конструкторские решения и найти границы переходных зон работы для вихревых аппаратов, нельзя определить, достигнута ли максимальная эффективность в существующих технологиях, и создать теорию подобия сильно закрученных потоков, необходимую для решения задач масштабирования вихревых аппаратов.

При изучении трехмерных нестационарных течений важную роль играют модели трехмерных вихревых структур. Поле скорости, индуцированное прямолинейной вихревой нитью или вихревой нитью кольцевой формы описывается достаточно простыми формулами (см., например, [Г. Ламб, 1947]]. Эти простейшие модели послужили основой создания моделей колоннообразных вихрей и вихревых колец, а также для развития вихревых методов расчета плоских и осе-симметричных течений, обзор которых можно найти в книге С.М. Белоцерков-ского, М.И. Ништа [1978]. Третья по сложности после прямолинейной и кольцевой - винтовая вихревая нить оказалась изучена совершенно недостаточно. Обычное представление скорости через интеграл Био-Савара приводит к интегрированию осциллирующих функций на бесконечном интервале. J.C. Hardin [1982] нашел подходящее преобразование интегралов и получил аналитическое решение для поля скорости, индуцированного винтовой вихревой нитью в безграничном пространстве. Поле скорости при этом было записано через бесконечные ряды от модифицированных функций Бесселя. Дополнительная сложность возникает при попытке моделирования ограниченных течений, в частности, течений в трубах. В отличие от прямолинейных и кольцевых вихревых нитей, для которых наличие стенок можно учесть с помощью принципа отражения, для винтовой геометрии вихря последний не работает.

За последние годы был достигнут значительный прогресс в изучении винтообразных вихревых структур в ограниченных закрученных потоках. Теоретические исследования здесь базируются на точном решении для поля скорости, индуцированного простейшей структурой, - бесконечно тонкой винтовой вихревой нитью в цилиндрической трубе. Решение впервые получил В.Л. Окулов [1993], сведя задачу к решению дифференциального уравнения на функцию тока. Поле скорости, как и у Hardin [1982] представлено через бесконечные ряды от модифицированных функций Бесселя, но удовлетворяет условию непротекания на стенке трубы. Еще одно преимущество работы В.Л. Окулова заключается в прямом выделении особенностей в найденном решении. Сингулярные члены выражены через элементарные функции и содержат основную информацию о структуре течения.

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

Анализ экспериментальных данных по закрученным потокам [И.И. Смульский, 1992; А.А. Халатов, 1989; А.Н. Штым, 1985; А. Гупта и др., 1987; Escudier, 1984; Escudier et al., 1980, 1982 и др.] дает основание утверждать, что даже в случае осевой симметрии вихревые линии в таких течениях не прямолинейные, а винтообразные. Это означает, что профиль осевой скорости - неоднородный. Для аппроксимации опытных данных исследователи подбирают функции чаще всего интуитивно. В то же время, если предположить, что вихревые линии - канонические винтовые спирали, то удается описать обе компоненты скорости - окружную и осевую - через одну функцию, получаемую путем интегрирования радиального распределения осевой компоненты завихренности.

В диссертации построен целый класс моделей колоннообразных вихрей с винтовыми вихревыми линиями, пригодный для описания структуры различных закрученных течений. В результате обобщены известные модели вихрей - вихря Рэнкина [Г. Дамб, 1947], вихря с дробно-рациональной функцией распределения завихренности [Scully, 1975] и вихря Ламба (см. [Hopfmger & van Heijst, 1993].

Модель винтового вихря с ядром конечного размера давно привлекает внимание ученых-гидродинамиков. В первую очередь это связано с прикладными задачами: следом за гребным винтом [Н.Е. Жуковский, 1912, 1937а, 19376], пропеллером или гидротурбиной [Fanelli, 1989]. Однако факт конечного использовался лишь для оценки скорости самоиндуцированного движения вихря. Распределение скоростей в течении, индуцированном таким вихрем, до сих пор найдено не было. В данной работе сделан первый шаг на пути построения такой модели, а именно - найдены профили скоростей, осредненных по угловой координате. Если в эксперименте винтообразный вихрь совершает вращательное движение, то осреднение по углу будет соответствовать осреднению по времени, что позволяет проводить сопоставление опытных и модельных профилей. Тем самым по осредненным профилям скорости можно восстановить локальную пространственную структуру потока. В более полной постановке требуется решение задачи о сшивке вихревого и потенциального течений (при однородном распределении осевой завихренности в ядре) и определение формы ядра. В предложенной же модели рассмотрено круговое сечение ядра, что справедливо либо для тонких ядер, либо, наоборот, для слабоискривленных колоннообразных вихрей.

Теоретические подходы к оценке скорости самоиндуцированного движения винтовых вихревых струюур достаточно широко пред'ставлены в литературе. Однако все исследователи ограничиваются здесь изучением вихревых нитей в безграничном пространстве. Это объясняется тем, что для описания вихрей используется закон Био - Савара, который является удобным инструментом в неограниченных областях и становится мало пригодным для описания вихревых структур в присутствии границ (например, при описании закрученных потоков в трубах или рабочих участках вихревых аппаратов). Тем не менее, накопленный материал весьма полезен для дальнейших исследований. Кратко остановимся на основных результатах.

Известно [Дж. Бэтчелор, 1973; Kida, 1981; и др.], что искривленная вихревая нить в общем случае движется и деформируется не только под воздействием внешнего потока, но и за счет самоиндуцированного движения, причем в некоторых случаях искривленная нить не меняет своей формы. К примерам таких нитей относятся вихревое кольцо и винтовой вихрь, который может вращаться и двигаться поступательно вдоль своей оси. Изучение винтовых вихрей в несжимаемой идеальной жидкости имеет длинную историю. Впервые слабоискривлен-ный винтовой вихрь с большим шагом винта исследовал Lord Kelvin [1880]. Н.Е. Жуковский [1912], применяя закон Био-Савара, приближенно получил формулу для определения скорости самоиндуцированного вращения винтового вихря в безграничном пространстве. Позднее Levy & Forsdyke [1928] рассмотрели вихрь с малым шагом. Другие фундаментальные результаты по приближенному описанию динамики винтовых вихревых нитей в безграничном пространстве можно найти в целом ряде широко известных работ [Widnall, 1972; Moore & Saffman, 1972; Adebiyi, 1981; Ricca, 1994; и др.]. В последних перечисленных работах предложен, по-видимому, наиболее удачный подход, заключающийся в представлении скорости винтового вихря через скорость вихревого кольца эквивалентной кривизны. Но даже и в этом случае остается проблема вычисления интегралов от осциллирующих функций на бесконечном интервале.

Возможность решить задачу о самоиндуцированном движении винтового вихря и представить результат не в квадратурах, а в окончательной явной форме появилась после доработки методики прямого выделения особенностей из бесконечных рядов от модифицированных функций Бесселя [Kuibin & Okulov, 1998]. Вывод формул для скорости движения, как одиночного винтового вихря, так и их системы в безграничном пространстве и в цилиндрической трубе сделан в данной диссертации.

В вихревых аппаратах первое наблюдение винтообразного вихря связано с его самоиндуцированным вращением. Были обнаружены интенсивные низкочастотные окружные пульсации потока, которые назвали прецессией вихревого ядра. Правда, установить, что причиной прецессии является самоиндуцированное вращение винтового вихря, удалось не сразу. Наиболее детально это явление было изучено, пожалуй, для вихревых горелок. В монографии А. Гупта с соавторами [1987] проанализировано влияние прецессии вихревого ядра на горение. Однако то, что прецессия вызывается вращением винтового вихря (жгута), по-видимому, впервые было показано при изучении закрученных потоков в отсасывающих трубах гидротурбин. Идентификация винтовой формы пламен при горении в закрученном потоке была получена много позднее в работах Ishizuka [1984, 1989]. Теоретическое описание формы пламени в закрученном потоке в литературе до настоящего времени не встречалось.

Одно из первых наиболее полных описаний прецессирующего вихря с воспроизведением его винтовой структуры сделал R.C. Chanaud [1965]. Изучению прецессии вихревого ядра в вихревой камере с удлиненным рабочим участком посвящен цикл работ R.F. Guarga с соавторами [1985а, 1985b, 1985с, 1985d]. Они измеряли частоту пульсаций давления на стенке (частоту прецессии) и од-10'• новременно визуализировали трехмерную структуру вращающегося винтового^ вихря, порождающего эти пульсации. С.С. Кутателадзе, Э.П. Волчков, В.И. Терехов [1987] также изучали прецессирующий поток с одновременной визуализацией. Изучение трехмерной структуры течения в пылевых циклонах проведено Yazdabadi, Griffiths & Syred [1994] с помощью ЛДА. Однако в этой работе, как и в монографии А. Гупта с соавторами [1987], прецессия вихревого ядра не отождествляется с самоиндуцированным вращением винтообразной вихревой структуры. Несмотря на обнаружение связи между прецессией и возникновением рециркуляционной зоны (противотока вдоль оси) в циклоне, авторы не делают заключения о существовании винтового вихря (или вихрей), который является причиной, как первого, так и второго явления. В прозрачной модели гидроциклона винтовая структура потока была выявлена А.Г. Лопатиным [1991].

Перечисленные экспериментальные данные относятся к режиму с вращающимся винтовым вихрем. Работ по второй форме организации потока в вихревых аппаратах - с неподвижным винтовым вихрем - до недавнего времени практически не было. Данные о существовании таких режимов, конечно, имелись - это режимы с подавленной прецессией вихревого ядра [А. Гупта и др., 1987], режимы без пульсаций [С.С. Кутателадзе и др., 1987] и др. Однако явно структура потока здесь выделена не была. Kumar & Conover [1991, 1993] с помощью визуализации течения провели исследование трехмерного вращающегося потока в цилиндрическом циклоне с тангенциальными входом и выходом. Фотографическим путем ими показано существование устойчивой трехмерной структуры потока. Основной поток движется по спиральной траектории, окружающей вторичный поток, в котором имеется неподвижный одиночный вихрь. Ось вихря искривлена, но на ее винтообразный характер авторы явно не указывают. Установлено, что такая структура потока не меняется в диапазоне чисел Рейнольдса от 15000 до 60000. Alekseenko & Shtork [1992], проводя исследования на модели топки квадратного сечения, впервые в закрученном потоке, кроме достаточно хорошо изученного вращающего вихря, визуализировали стационарные (неподвижные) винтовые вихри. При вариации граничных условий они выявили режимы течения с неподвижными винтообразными вихревыми структурами (одиночными и двойными).

Концентрированные вихри винтовой формы весьма распространены в природе и технике. Помимо описанных выше вихрей, укажем далеко не полный их список: воронка в жидкости, истекающая из сосуда через донное отверстие [М. Ван-Дайк, 1986]; смерчи [Д. Сноу, 1984]; продольные вихри при обтекании дельтовидного крыла под большим углом атаки [Peyne et al., 1988]; продольные вихри в турбулентном пограничном слое [Kim et al., 1971] и в следе за обтекаемым телом на пластине [Tani et al., 1962]; система вихревых шнуров, образующаяся за несимметричной струей, вдуваемой в поток [Д.М. By и др., 1989]; система вихревых шнуров во вращающемся слое жидкости, подогреваемой снизу [Boubnov & Golitsyn, 1986]; вихревые шнуры в восходящем над закрученной жидкостью паре [В.А. Владимиров, 1977]; вихревые нити в модели турбулентности для сверхтекучего гелия [Р.Д. Доннели, 1989]; концевые вихри за винтом и пропеллером [Н.Е. Жуковский, 1912, 1937а, 19376]; вихревые нити при обтекании лункообразной каверны [Г.И. Кикнадзе и др., 1986]. Данная тематика широко освещена в литературе. Некоторые экспериментальные наблюдения винтовых вихрей приведены в хорошо известных обзорах [Сэффмэн, 2000; Widnall, 1972; Leibovich, 1984; Escudier, 1988]. Однако описанные в них результаты сконцентрированы в основном на явлении неустойчивости закрученных течений, подобном распаду вихря пузырьковой и спиральной форм. В общем объеме исследований вихрей винтовой формы более детально изучен случай вращающегося винтового вихря (или явления прецессии вихревого ядра). Исследования крупномасштабных вихревых структур, формирующихся в интенсивно закрученных потоках, практически отсутствуют. Хотя Takaki & Hussain [1984] и нашли стационарное решение в виде двух переплетенных винтовых вихрей с одинаковой циркуляцией, экспериментально получить неподвижную двойную винтовую структуру долгое время не удавалось. Этот факт говорит о том, что экспериментальное изучение вихревых структур в закрученных потоках - задача достаточно сложная. И только недавно в экспериментах С.В. Алексеенко и С.И. Шторка [1994]» эта структура была найдена.

Отдельно следует остановиться на исследованиях, посвященных эффекту энергоразделения в трубках Ранка. Дело в том, что, несмотря на интенсивное изучение этого явления в течение последних пятидесяти лет, исследователи не пришли к единому мнению о его механизме. Гертлеровские вихри [Stephan et al., 1983]; сжимаемость [Amitani et al., 1983; Д.К. Зайцев, E.M.-Смирнов, 1994]; турбулентный перенос тепловой энергии [Linderstrom-Lang, 1971]; акустические потоки [Kurosaka, 1982]; перестройка вынужденного вихря в свободный [А.П.Меркулов, 1969]; влияние вязкостных эффектов [В.И. Кузнецов, 1989; В.А. Баранов и др., 1994] и многое другое предлагалось как возможные механизмы.

В эффекте энергетического разделения, скорее всего, значительную роль играет прецессирующая винтообразная вихревая структура. А. Гупта с соавторами [1987] указывают на значительное влияние прецессии вихревого ядра на процесс энергоразделения. Н.А. Артамонов, Б.Ф. Абросимов и М.З. Максименко [1987] при изучении динамики течения в вихревой трубе сделали важное наблюдение о винтовой структуре потока в трубе Ранка. С.В. Лукачевым [1981] было дано объяснение регулярных пульсаций потока (типа прецессии) в вихревых трубах как следствие возникновения крупномасштабных вихревых структур. В своей последующей работе он предложил гипотезу о возможности существования двойных и тройных вихревых образований на основе анализа сдвига фаз между пульсациями давления, измеренными в двух диаметрально противоположных точках вихревой трубы [С.В. Лукачев, 1984]. Развитие этих представлений было продолжено в работе Ю.А. Кныша [1988]. Однако возникновение вихревых структур в вихревых трубах оба исследователя связывают не с неустойчивостью интенсивно закрученных потоков, а с неустойчивостью типа контактного разрыва на поверхности раздела периферийного подогретого газа и центрального холодного ядра потока. В связи с этим базовыми вихревыми структурами они считают вихревые кольца в виде тора, которые в реальном несимметричном потоке деформируются и сворачиваются в спиральные вихревые структуры. Точное представление о винтообразной форме вихревой структуры, образующейся в результате неустойчивости потока, в указанных работах отсутствует. Вместе с тем при исследовании структуры потока в другом аппарате - вихревом генераторе звука С.В. Лукачев [1988] связывает образование прецессирующего винтового вихря с распадом вихря. Наконец, явное подтверждение наличия винтообразных винтовых структур в вихревой трубе получено с помощью специальных методов визуализации в работе В.А. Арбузова с соавторами [1997].

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

При попытках сопоставления моделей вихрей с реальными течениями закономерно возникает вопрос о влиянии вязкости. Конечно, существуют области науки, где вязкость несущественна. Например, - при изучении сверхтекучего гелия. В этом случае применяют классическую модель вихревых нитей. Но, вследствие хаотического поведения нитей, рассматривают лишь статистические и макроскопические характеристики [Р.Д. Доннели, 1989; Nemirovskii, 1998а, 19986]. Однако и в обычной жидкости возможны условия, когда непосредственным влиянием вязкости можно пренебречь. Дж. Бэтчелор [1973] назвал их течениями эффективно невязкой жидкости. Такие условия могут возникать в течениях с достаточно большими числами Рейнольдса (Re > 103). Необходимо также, чтобы область завихренной жидкости находилась на достаточном удалении от твердых границ - вне пограничных слоев. Следует заметить, что применение моделей вихрей с ядром конечного сечения предполагает, что вязкость уже оказала свое действие на этапе формирования вихря. Размер ядра определяется Геометрическими параметрами установки-генератора и, в меньшей степени, числом Рейнольдса, т.е. скоростью потока и вязкостью жидкости [Maxworthy, 1977]. Еще один довод в пользу применения моделей вихрей заключается в том, что простыми зависимостями аппроксимируется поле завихренности, а сопоставление ведется по полю скорости, которое находится путем интегрирования. При этом повышается степень гладкости функций (например, кусочно-непрерывной завихренности соответствует непрерывное распределение скоростей, кусочно-линейной завихренности - гладкое и т.д.).

Новые возможности открываются перед исследователями в связи с развитием дифференциально-геометрического формализма, группового анализа и вариационных методов применительно к описанию течений несжимаемой жидкости. Так А.В. Кистовичем, Ю.Д. Чашечкиным [2000] предложен геометрический подход: линии тока рассматриваются как геодезические на интегральных поверхностях. В результате осуществляется переход от уравнений Эйлера к системе уравнений для параметров, описывающих геометрию интегральных поверхностей. Такой подход позволил авторам построить решение в форме совокупности модифицированных вихрей Рэнкина, которые в отличие от классического обладают ограниченными интегральными инвариантами (энергией, импульсом, моментом импульса, циркуляцией и спиральностью). Теми же авторами [А.В. Кистович, Ю.Д. Чашечкин, 2001] развит метод регулярного поиска дискретных симметрий дифференциальных уравнений, обладающий большими возможностями чем стандартный метод групп Ли и обобщающий упомянутый геометрический подход на другие классы дифференциальных уравнений.

Интересен подход, основанный на вариационном принципе максимума информационной энтропии, заимствованный из теории информации [Ю.Н. Григорьев, В.Б. Левинский, 1986; Ю.Н. Григорьев, В.Б. Левинский, В.П. Нартов, 1989]. Его применение позволило авторам построить модели цепочек вихрей, в которых распределение завихренности неоднородно, в отличие от моделей когерентных структур в сдвиговых слоях, построенных на базе вихрей Рэнкина или эллиптических вихрей Кирхгофа с равномерным распределением завихренности. Другое применение вариационных принципов - для построения методов дискретных вихревых частиц будет подробно рассмотрено ниже.

Классические модели элементарных вихревых структур - прямолинейная и кольцевая вихревые нити - явились основой построения методов расчета таких нестационарных явлений, как отрыв. Вихревые методы расчета развивались наряду с другими подходами к решению этой проблемы [П. Чжен, 1972-1973; Дж. Уильяме, 1979; Л.В. Гогиш и др., 1974; В.В. Сычев, 1987; В. Маккроски, 1977]. Первые попытки теоретического изучения отрыва сводились к построению стационарных моделей отрывных течений идеальной несжимаемой жидкости. В ходе их развития были сформулированы такие принципиальные положения теории отрыва, как гипотеза о конечности скорости на острых кромках обтекаемых тел (Н.Е. Жуковский, С.А. Чаплыгин, Кутга) и идея Н.Е. Жуковского о присоединенных вихрях [Н.Е. Жуковский, 1937а]. Эта идея дала возможность не только объ-^ яснить механизм образования подъемной силы крыла в идеальной жидкости, но и позволила создать удобный аппарат для численных расчетов.

Дальнейшее развитие теории отрыва связано с разработкой нестационарных моделей [В. Маккроски, 1977; С.М. Белоцерковский, М.И. Ништ, 1978; С.К. Бетяев, 1983]. Вследствие нелинейности уравнений гидродинамики при решении задач используются различные численные методы. В частности, в последнее время появляется все больше работ, в которых задачи о течении вязкой жидкости и газа при наличии отрывов исследуются путем прямого численного решения уравнений Навье-Стокса. Обзоры по этой теме содержатся в статьях [П. Чжен, 1972-1973; Дж. Уильяме, 1979; Л.В. Гогиш и др., 1974; Batchelor, 1956; Lissaman, 1983]. Прямое численное решение уравнений Навье-Стокса при больших числах Рейнольдса Re сопряжено, однако, с определенными трудностями. Во-первых, вследствие наличия малого коэффициента при старшей производной, пропорционального Re , появляются узкие области с большими локальными градиентами функций, что приводит к потере точности и устойчивости решения. Во-вторых, не во всех случаях известно существование стационарного устойчи-'т вого решения при достаточно больших числах Re. В-третьих, решение нестационарных задач требует существенных затрат машинного времени.

Поэтому в настоящее время все более широко используются более экономные и удобные в применении вихревые методы, основанные на лагранжевом подходе к описанию движения жидкости. Наиболее известным из вихревых методов является метод точечных вихрей (МТВ), применяемый для расчета плоских течений несжимаемой невязкой жидкости. Суть МТВ состоит в том, что поле завихренности аппроксимируется системой прямолинейных вихревых нитей. Впервые МТВ был применен Rosenhead [1931], который исследовал эволюцию вихревой пелены, моделируя ее точечными вихрями. С появлением ЭВМ этот метод начал широко использоваться для численного моделирования нестационарных отрывных течений [С.М. Белоцерковский, М.И. Ништ, 1978; К.П. Ильи

• чев, С.Н. Постоловский, 1972] и турбулентных струй [А.Б. Айрапетов, 1977,

С.М. Белоцерковский и др., 1985]. МТВ достаточно эффективен для класса задач, где требуется получить интегральные характеристики (например, силы, действующие на тело, размер отрывной зоны и т.п.), а в ряде случаев позволяет получать и правильную качественную информацию о течении. Применяют МТВ и в задачах с образованием пограничных слоев - как составную часть, пригодную для описания течения в зонах, где вязкие эффекты несущественны [С.М. Белоцерковский и др., 1988]. Что же касается более тонких характеристик, связанных с пульсационным характером течения, с развитием неустойчивости, исследованием структуры отрывной зоны, то здесь возникают значительные ошибки. Дело в том, что для моделирования указанных эффектов требуется увеличивать число точечных вихрей, а при этом происходит стохастизация их траекторий [Kuwahara & Takami, 1973; Birkhoff & Fisher, 1959], что связано с неограниченным ростом скоростей, индуцируемых вихрями при взаимном сближении. По этой же причине погрешность в определении поля скорости может быть очень значительной, даже, если движение точечных вихрей является регулярным [Beale & Maida, 1985]. Очевидно, также, что с помощью точечных вихрей нельзя сколь-нибудь удовлетворительно аппроксимировать физически интересные поля завихренности. Точечные вихри оказываются просто слишком сингулярными. В связи с этим в последние годы были развиты методы моделирования течений жидкости системой вихревых частиц, имеющих, в отличие от 5-образных точечных вихрей, конечные размеры [Kuwahara & Takami, 1973; Chorin, 1973; Chorin & Bernard, 1973; Moore, 1974; Leonard, 1980; Meng & Thomson, 1978]. Эти модели можно назвать феноменологическими, так как в них содержится ряд параметров, выбор которых никак не связан с физическими закономерностями задачи.

В работах А.Н. Веретенцева, В Я. Рудяка, Н.Н. Яненко [1982, 1985, 1986] развит вариационный принцип построения вихревых моделей, свободный от недостатков, характерных для упомянутых выше модели точечных вихрей и моделей вихревых частиц (vortex blobs). Более того, построенные модели сохраняют свойства консервативности, присущие континуальной среде, т. е. в них выполняются законы сохранения энергии, импульса и момента импульса. Исследование вопроса сходимости решений, полученных с помощью построенных моделей, к точным решениям уравнений Эйлера приводится в работах тех же авторов [1982, 1986].

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

Большинство упоминаемых здесь вихревых моделей строятся в рамках модели невязкой жидкости. Но отрыв потока всегда связан с вязким взаимодействием потока с поверхностью - с образованием пограничных слоев. В действительности же, при больших числах Рейнольдса важен сам факт наличия вязкости, а детали отрывного обтекания тел практически не зависят от величины вязкости среды и определяются инерционным взаимодействием в жидкостях и газах, описываемым нестационарным уравнением идеальной среды. Поэтому разрабатываемый в диссертации метод оказывается пригодным для описания отрывных течений при больших числах Рейнольдса.

Классической задачей в классе отрывных течений считается задача обтекания бесконечного клина, сформулированная еще Прандтлем [Prandtle, 1924] и моделирующая течение около тела (поверхность которого имеет изломы) в локальной области в окрестности излома. Существенный вклад в решении указанной проблемы внес А.А. Никольский [1957], получивший уравнение эволюции вихревой пелены, возникающей в результате отрыва, с использованием суммарной интенсивности вихревого следа в каждый момент времени в качестве ла-гранжевой координаты точки вихревой пелены, покидающей в этот момент кромку клина. В работе исследованы автомодельные режимы течения, когда потенциал внешнего (безотрывного) течения зависит от времени по степенному закону. Решение полученной таким образом автомодельной задачи применялось при исследовании некоторых стационарных пространственных задач для конических тел [А.А. Никольский, 1970; Ward, 1955] и при моделировании начальной стадии отрывного обтекания [JI.B. Гогиш, Г.Ю. Степанов, 1983; В.Ф. Молчанов,

1975]. Достаточно полное исследование разгонного вихря за бесконечным клином, основанное на численном решении автомодельной задачи (задача сведена к интегро-дифференциальному уравнению с одной переменной) проведено Pullin [1978]. Получена картина сворачивания в вихрь сходящей с клина вихревой пелены для различных углов раствора клина и в широком диапазоне показателей автомодельности. Проанализирована форма вихревой пелены в окрестности кромки и в ядре разгонного вихря. Вопрос же о рамках применимости автомодельного решения в конкретных задачах остается открытым.

В работе А.А. Никольского [1957] особо выделен случай вырожденного автомодельного течения, когда вся циркуляция стягивается в точечный вихрь. Реализация такого течения представляет интерес, как в численной модели, так и на практике: таким образом можно создавать вихри с наибольшей концентрацией завихренности. Актуальность их изучения диктуется, в частности, задачами переноса импульса силы, транспортировки легких примесей и др. [В .И. Меркулов, 1989].

Экспериментальные исследования показывают, что в процессе формирования разгонного вихря гладкая спиралевидная структура образуется не всегда. В альбоме течений жидкости и газа Ван-Дайка продемонстрировано развитие вторичных вихревых структур на фоне разгонного вихря (см. фото Ж№ 81, 82, 241 в [М. Ван-Дайк, 1986]). Для прямолинейного сдвигового течения (слой сдвига, слой смешения) механизмы образования подобных (когерентных) структур хорошо изучены (см. обзоры в [Но & Huerre, 1984; Hussain, 1986; Е.В. Власов, А.С. Гиневский, 1986]). В частности, экспериментально установлено, что образование когерентных структур (КС) в возбуждаемом слое смешения происходит на длине волны возбуждения. Результаты экспериментов подтверждаются также и численными расчетами [Acton, 1976; Ghoniem & Ng, 1987], где методом вихревых частиц моделировалось образование и последующее спаривание крупномасштабных вихревых структур в сдвиговом слое.

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

На первый взгляд, след за тонкой пластиной представляет собой частный случай слоя смешения - когда скорости с обеих сторон пластины одинаковы. Тем не менее, переносить результаты исследования слоя смешения на задачу о следе нельзя. Это связано, в первую очередь, с симметрией профиля средней скорости относительно оси следа, и, как следствие, с наличием двух мод неустойчивости различной четности. О существовании двух мод неустойчивости в течениях с симметричным профилем скорости сообщалось в работах Michalke, Schade [1963], Tatsumi, Kakutani [1958]. Там же было показано, что в следе наиболее неустойчивой является антисимметричная мода (с нечетным профилем пульсаций продольной скорости). На этом основании в подавляющем большинстве работ по неустойчивости следа, особенно в экспериментальных (см., например, Sato, Kuriki [1961], Mattingly, Criminale [1972]), второй модой - симметричной - пренебрегали. Вместе с тем в экспериментах обычно возбуждается возмущение, являющееся суперпозицией двух мод. В работах Wygnanski, Champagne, Marasli [1986], Marasli, Champagne, Wygnanski [1989] были предприняты попытки экспериментально генерировать отдельно каждую из мод. С другой стороны, ясно, что если на линейной стадии обе моды развиваются независимо, то на нелинейной - взаимодействие между ними будегг существенно влиять на характер развития течения.

Выше было показано, что основной механизм развития неустойчивости на нелинейной стадии в слое смешения - субгармонический. Из общих соображений подобный механизм представляется возможным и для следа [С.Я. Герценштейн, Ю.М. Штемлер, 1977]. В работе Kelly [1968] было однако показано, что хотя для симметричного профиля принципиально и возможно резонансное взаимодействие нейтральных возмущений симметричной и антисимг метричной мод, коэффициент связи мод, вычисленный в рамках слабонелинейной теории типа теории Стюарта-Ватсона, тождественно равен нулю. Полученный в цитируемой работе при многочисленных весьма сильных предположениях негативный результат вряд ли можно считать окончательным и, кроме того, он не дает ответа на вопрос о том, как протекает нелинейная стадия развития неустойчивости в следе.

Развитие возмущений в течении типа следа за телом (точнее в течении с профилем продольной скорости и = 1 - 0,7ехр(-0,9у )) изучалось ранее численно в [С.Я. Герценштейн и др., 1985]. Здесь был получен целый ряд интересных результатов. В частности, было установлено выделение длинноволновых составляющих спектра. Причем при нелинейном взаимодействии таких возмущений наибольшее нарастание наблюдается у их разностной составляющей. Важным представляется и обнаружение вторичной неустойчивости конечно-амплитудных волновых режимов к поперечным трехмерным возмущениям.

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

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

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

В первой главе выводятся уравнения гидродинамики течений с винтовой симметрией и рассматривается класс течений с винтовыми вихревыми линиями. На основе решения, полученного B.JI. Окуловым [1993], проводится анализ трехмерных течений, индуцированных вихревой нитью винтовой формы при различных значениях определяющих параметров.

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

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

В четвертой главе развитым методом решается задача об отрывном обтекании кромки полубесконечной пластины при различном нестационарном набегающем потоке. На основе известного автомодельного решения производится тестирование метода. Изучается неустойчивость разгонного вихря и процессы образования на его фоне вторичных вихревых структур. Исследуется нелинейная стадия развития неустойчивости в плоском следе за пластиной. Основное внимание здесь уделено нелинейному взаимодействию основного и субгармонического возмущений симметричной и антисимметричной мод.

Глава V посвящена описанию и классификации вихревых структур, образующихся в закрученных потоках. Предложены новые интегральные характеристики для описания закрученных потоков — шаг винтовой вихревой структуры и скорость ее чисто поступательного перемещения вдоль оси. Проведена проверка существования винтовой симметрии в закрученных течениях. Для вихрей с прямолинейной осью установлена возможность существования композиции винтовых вихрей. Изучены причины образования левовинтовых- и правовинтовых вихрей в вихревых камерах. Установлена взаимосвязь явления прецессии вихревого ядра с вращающимся винтовым вихрем. Проведен анализ влияния структуры потока на энергоразделение в вихревой трубке Ранка и на форму пламени в закрученном потоке.

В шестой главе изучается движение вихрей в цилиндрических трубах. На основе развитого метода дискретных вихревых частиц решается задача о неустойчивости положения прямолинейного вихря в трубе к возмущениям завихренности локализованным на окружности некоторого радиуса. Выводится формула для скорости самоиндуцированного движения винтового вихря в безграничном пространстве и в цилиндрической трубе. Получена формула для частоты прецессии винтового вихря в трубе. Исследовано движение системы соосных винтовых вихрей в трубе.

В заключении сформулированы основные результаты и выводы работы. В приложении I приводятся примеры вычисления характерных размеров вихревых частиц различными способами. В приложении II дано доказательство теоремы о законе изменения вихревого импульса при отрывном обтекании пластины, используемого для оценки погрешности метода в п. 4.1 четвертой главы.

На защиту выносится:

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

2. Новый класс точных решений уравнений Эйлера для вихрей с прямолинейной осью и распределением завихренности, обладающей винтовой симметрией. Модель колоннообразного вихря, движущегося вдоль плоскости, с винтообразной структурой вихревых линий. Модель винтового вихря с ядром конечного размера.

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

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

5. Новые интегральные характеристики, отвечающие за структуру закрученных потоков: шаг вихревых линий и переносная скорость, совпадающая со скоростью на оси потока. Выявление локальной винтовой симметрией в реальных закрученных потоках. Модель суперпозиции правостороннего и левостороннего винтовых вихрей (один находится в другом). Модель вихря с изменяющейся винтовой симметрией в пространстве (правый винтовой вихрь переходит в левый). Классификация вихревых структур в закрученных потоках.

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

Основные результаты диссертации опубликованы в 31 работе:

1. Веретенцев А.Н., Куйбин П.А., Меркулов В.И., Рудяк В.Я. О выводе уравнений движения дискретных вихревых частиц для осесимметричных течений // Изв. СО АН СССР. Сер. техн. наук. - 1986. - № 10. - Вып. 2. - С. 45-50.

2. Веретенцев А.Н., Куйбин П.А., Рудяк В.Я. Моделирование формирования вихря на острой кромке полубесконечной пластины // Изв. СО АН СССР. Сер. техн. наук. - 1988. -№ 7, вып. 2. - С. 21-25.

3. Veretentsev A.N., Kuibin P.A., Rudyak V.Ya. Coherent structures in artifically excited free shear flows // On turbulence: Proc. 5th EPS Liquid State Conf., Moscow, October 16-21,1989. - Institute for Problem in Mechanics, Moscow. -P. 198-201.

4. Веретенцев A.H., Гешев П.И., Куйбин П.А., Рудяк В.Я. О развитии метода вихревых частиц применительно к описанию отрывных течений // Ж. вычисл. мат. и мат. физики. - 1989. - Т. 29, № 6. - С. 878-887.

5. Kuibin P.A., Rudyak V.Ya., Veretentsev A.N. Processes of the instability development in separated flows behind the plate // Separated flows and jets: Proc. IUTAM Symp., Novosibirsk, USSR, July, 9-13, 1990. - Berlin: Springer-Verlag, 1991.-P. 747-750.

6. Куйбин ПА., Рудяк В.Я. Изучение эволюции возмущений в следе за пластиной // Изв. вузов. Авиационная техника. - 1992. - № 1. - С. 29-32.

7. Куйбин П.А., Рудяк В.Я. Развитие неустойчивости в следе за пластиной, размещенной параллельно потоку // Изв. АН СССР. МЖГ. - 1992. - № 1. - С. 2632.

8. Куйбин П.А. Исследование динамики и неустойчивости отрывных течений и следов // Дисс. канд. физ.-мат. наук, ИТ СО РАН, Новосибирск. - 1993. - 140 с.

9. Борисов А.А., Куйбин П.А., Окулов B.J1. Моделирование течения и конвективного энергоразделения в вихревых трубах // Изв. СО РАН. Сиб. физ.-техн. журн. - 1993. -№ 1.-С. 30-38.

Ю.Борисов А.А., Куйбин П.А., Окулов В.Л. Описание конвективного теплопере-носа в вихревой трубе // ДАН. - 1993. - Т. 331, № 1. - С. 28-31.

П.Борисов А.А., Куйбин П.А., Окулов В.Л. Спиновое и нормальное горение в закрученных потоках // Письма в ЖТФ. - 1993. - Т. 19, вып. 14. - С. 13-17.

12. Борисов А.А., Куйбин П.А., Окулов В.Л. К вопросу о горении газа в закрученном потоке // ФГВ. - 1993. -№ 5. - С. 107-109.

13.Borissov A.A., Kuibin P.A., Okulov V.L. Flame Shapes in Swirl Flow // Russian J. of Eng. Thermophys. - 1993. - Vol. 3, № 3. - P. 243-255.

14.Borissov A.A., Kuibin P.A., Okulov V.L. Convective Heat Transfer and its Action on the Ranque Effect in Vortex Tube // Experimantal and Numerical Flow Visualization: Proc. Int. Conf. FED Vol.172. Book no. H00849. 1993. ASME.

15. Куйбин П.А., Окулов В.Л. Определение частоты прецессии винтового вихря // Письма в ЖТФ. - 1994. - Т. 20, вып. 7. - С. 32-35.

16. Алексеенко С.В., Куйбин П.А., Окулов В.Л., Шторк С.И. Характеристики закрученных потоков с винтовой симметрией // Письма в ЖТФ. - 1994. - Т. 20, вып. 18.-С. 33-39.

17.Borissov A.A., Kuibin P.A., Okulov V.L. Calculation of Ranque Effect in Vortex Tube 11 Acta Mechanica. - 1994. - [Suppl.] № 4. - P. 289-295.

18.Borissov A. A., Kuibin P.A., Okulov V.L. The Helical Structure of Flame in a Swirl Flow // Proc. Zel'dovich Memorial "Combustion, detonation, shock waves". Vol. 2. - 1994.-P. 236-238.

19. Алексеенко C.B., Куйбин П.А., Окулов B.Jl., Шторк С.И. Стационарный вихрь с переменной винтовой симметрией // ДАН. - 1995. - Т. 345, № 5. - С. 611-614.

20. Alekseenko S.V., Kuibin Р.А., Okulov V.L., Shtork S.I. Stationary Vortex Structures in Intensively Swirling Flows // Proc. conf. "Numerical Methods in Laminar and Turbulent Flow". Vol.9. Part 1. - Atlanta, U.S.A. - 1995. P. 382-393.

21. Alekseenko S.V., Kuibin P.A., Okulov V.L., Shtork S.I. Large Scale Vortex Structures in Intensively Swirling Flows // Proc. conf. "Experimantal and Numerical Flow Visualization". FED Vol. 218. ASME. Hilton. U.S.A. - 1995. P. 181-188.

22. Alekseenko S.V., Kuibin P.A., Okulov V.L., Shtork S.I. Large Scale Vortex Structures in Intensively Swirling Flows // Proc. XXVI IAHR Conference, 11-15 September 1995, Vol. 2. - 1995. P. 90-95.

23. Куйбин П.А., Окулов В.Л. Одномерные решения для течений с винтовой симметрией // Теплофизика и аэромеханика. - 1996. - № 4. - С. 311-315.

24. Alekseenko S.V., Kuibin Р.А., Okulov V.L., Shtork S.I. (1997) Helical vortices in power engineering // N.V. Medvetskaya and R.S. Gromadskaya (eds.), The Physics of Heat Transfer in Boiling and Condensation, Institute for High Temperature, Russian Academy of Sciences, Moscow. P. 471-476.

25. Kuibin, P.A., Okulov, V.L. Self-induced motion and asymptotic expansion of the velocity field in the vicinity of helical vortex filament // Phys. of Fluids. - 1998. -Vol. 10, №3.-P. 607-614.

26.Kuibin, P.A., Okulov, V.L. Self-induced motion of helical vortex // Proc. of IUTAM-Symposium on Dynamics of Slender Vortices. Kluwer Ac. Pub. - 1998. -P. 55-62.

27. Alekseenko, S.V., Kuibin, P.A., Okulov, V.L., and Shtork, S.I. Theory of helical vortices // ibid. - P. 255-263.

28. Alekseenko, S.V., Kuibin, P.A., Okulov, V.L., and Shtork, S.I. Helical vortices in swirl flow // J. Fluid Mech. - 1999. - Vol. 382. - P. 195-243.

29. Kuibin, P.A. On motion of a double helical vortex in a cylindrical tube // Proc. IUTAM Symp. "Tubes, sheets and singularities in fluid dynamics". Kluwer. -2002.-P. 81-86.

30. Kuibin P.A., 2003. On structure of vorticity field in the spherical Hicks vortex // J. Eng. Thermophys. - Vol. 12, №. 3. - P. 22-29.

31. Алексеенко С.В., Куйбин П.А., Окулов B.JL, Введение в теорию концентрированных вихрей. - Новосибирск: Институт теплофизики им С.С. Кутателадзе СО РАН, 2003.

Основные материалы диссертации докладывались и обсуждались: на II Республиканской школе-семинаре "Актуальные вопросы теплофизики и физической гидрогазодинамики" (Алушта, 1987), на V Всесоюзной школе "Современные проблемы теплофизики" (Новосибирск, 1988), на Всесоюзном семинаре "Отрывные и струйные течения" (Новосибирск, 1988), на XI Всесоюзной школе по численным методам механики вязкой жидкости (Свердловск, 1988), на V конференции Европейского физического общества по турбулентности (Москва, 1989), на IV Всесоюзном симпозиуме "Методы дискретных особенностей в задачах математической физики" (Харьков, 1989), на V Школе по методам аэрофизических исследований (Абакан, 1989), на XII Всесоюзной конференции "Аэроупругость турбомашин" (Рига, 1989), на VI и VII (Москва, 1988 и 1990) Всесоюзных школах "Нелинейные задачи теории гидродинамической устойчивости", на Международном симпозиуме "Образование крупномасштабных структур в сплошной среде" (Пермь - Москва, 1990), на Международном симпозиуме по отрывным течениям и струям (Новосибирск, 1990), на VI и VII ASME конференциях (США, 1993 и 1995) «Экспериментальная и численная визуализация потоков», на международной конференции "Горение, детонация, ударные волны" (Москва, 1994), на международной конференции "Численные методы для ламинарных и турбулентных течений" (США, 1995), на XXVI конференции международной ассоциации исследователей-гидравликов (Швеция, 1995), на международной конференции "Физика и теплообмен при кипении и конденсации" (Москва, 1997), на IUTAM-симпозиуме "Моделирование и идентификация организованных структур в потоках" (Дания, 1997), на II и III конгрессах по индустриальной и прикладной математике (Новосибирск, 1996, 1998), на VI и VIII конференциях «Устойчивость течениий гомогенных и гетерогенных жидкостей» (Новосибирск, 1997, 2001), на IUTAM-симпозиуме по динамике тонких вихрей (Германия, 1998), на международном симпозиуме "Актуальные проблемы физической гидроаэродинамики" (Новосибирск, 1999), на X и XI международных конференциях "Потоки и структуры в жидккости" (С.-Петербург, 1999, Москва, 2001) на 4 Европейской конференции по механике жидкости (Нидерланды, 2000), на IUTAM-симпозиуме «Трубки, пелены и сингулярности в динамике жидкости» (Польша,

2001), на 426 коллоквиуме ЕВРОМЕХ и ЭРКОФТАК "Закрученные потоки" (Норвегия, 2001), на XXVI Сибирском теплофизическом семинаре (Новосибирск,

2002), а также на семинарах в Институте теплофизики им. С.С. Кутателадзе СО РАН, в Институте теоретической и прикладной механики СО РАН, в Санкт-Петербургском политехническом университете, в Институте проблем механики РАН, в Институте энергетики и механики жидкости (Датский технический университет, Копенгаген) и в Институте гидро- и термодинамики (Рурский университет, Бохум, Германия).

Заключение диссертации по теме «Механика жидкости, газа и плазмы», Куйбин, Павел Анатольевич

ЗАКЛЮЧЕНИЕ

По содержанию данной работы можно сформулировать следующие основные результаты и выводы:

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

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

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

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

Уточнена структура поля завихренности в сферическом вихре Хикса. В частности, установлено, что вихревые линии и линии тока, вьющиеся вокруг торооб-разных поверхностей тока, имеют конечный период.

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

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

2. С помощью обобщенного метода вихревых частиц систематически изучена задача о формировании автомодельного разгонного вихря на кромке полубесконечной пластины при различном законе изменения во времени скорости набегающего потока. В классе неавтомодельных задач проанализированы, в частности, течения, при которых концентрация завихренности в ядре вихря максимальна; описано течение, возникающее в режиме "разгон - торможение" (образующийся при этом грибовидный вихрь хорошо воспроизводит картину, наблюдаемую в эксперименте).

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

Впервые решена задача о развитии неустойчивости в разгонном вихре. Установлено, что возмущения различного типа преобразуются на кромке в возмущения завихренности. Эволюция последних приводит к образованию вторичных вихревых структур на фоне основного вихря. Картина течения, полученная в расчетах хорошо согласуется с наблюдавшейся в экспериментах. Найден оригинальный способ анализа спектральных характеристик развивающихся возмущений при сильно нестационарном основном потоке. В результате анализа обнаружена генерация субгармонических возмущений и активное ее резонансное взаимодействие с основной модой.

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

5. Введены новые интегральные характеристики, отвечающие за структуру закрученных потоков: шаг вихревых линий и переносная скорость, совпадающая со скоростью на оси потока. На основе анализа многочисленных экспериментальных данных установлено, что закрученные потоки обладают локальной винтовой симметрией. Продемонстрирована возможность реализации суперпозиции правостороннего и левостороннего винтовых вихрей (один находится в другом). На основе анализа влияния граничных условий на вид вихревой структуры предсказана возможность генерации вихря с изменяющейся винтовой симметрией в пространстве (правый винтовой вихрь переходит в левый). Стационарный вихрь такой формы позже был обнаружен в эксперименте. Проведена классификация вихревых структур в закрученных потоках по геометрическим признакам.

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

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

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

Список литературы диссертационного исследования доктор физико-математических наук Куйбин, Павел Анатольевич, 2003 год

1. Айрапетов А.Б., 1977. О статистических свойствах системы дискретных вихрей, моделирующих турбулентное струйное течение // Турбулентные течения. -М.: Наука,-С.188-195.

2. Алексеенко С.В., Куйбин П.А., Окулов В.Л., Шторк С.И., 1994. Характеристики закрученных потоков с винтовой симметрией // Письма в Журн. техн. физики. -Т. 20, Вып. 18.-С. 33-39.

3. Алексеенко С.В., Куйбин П.А., Окулов В.Л., Шторк С.И., 1995. Стационарный вихрь с переменной винтовой симметрией // Докл. РАН. Т. 345, № 5. - С. 611-614.

4. Алексеенко С.В., Окулов В.Л., 1996. Закрученные потоки в технических приложениях (обзор) // Теплофизика и аэромеханика. Т. 3, № 2. - С. 101-138.

5. Алексеенко С.В., Шторк С.И., 1994. Экспериментальное наблюдение взаимодействия вихревых нитей // Письма в Журн. эксперим. и теорет. физики. Т. 59, Вып. 11. - С. 746-750.

6. Алексеенко С.В., Шторк С.И., 1997. Бегущий распад вихря // Письма в Журн. техн. физики. Т. 23, Вып. 22. - С. 24-28.

7. Арбузов В.А., Дубнищев Ю.Н., Лебедев А.В., Правдина М.Х., Яворский Н.И., 1997. Наблюдение крупномасштабных структур в вихревой трубке и эффект Ранка // Письма в Журн. техн. физики. Т. 23. - Вып. 23. - С. 84-90.

8. Артамонов Н.А., Абросимов Б.Ф., Максименко М.З., 1987. Динамика струйных течений в вихревой трубе // ИФЖ. Т. 53, № 6. - С. 906-911.

9. Баранов В.А., Викульцев Ю.А., Рис В.В., Смирнов Е.М., 1994. Экспериментальное исследование влияния числа Рейнольдса на энергоразделение в вихревой трубе // Процессы горения и охрана окружающей среды: Материалы II Всерос. конф. Рыбинск. - С. 15-20.

10. Белоцерковский С.М., Гиневский А.С., 1995 Моделирование турбулентных струй и следов на основе метода дискретных вихрей. М.: Физматлит. - 368 с.

11. Белоцерковский С.М., Дворак А.В., Хлапов Н.В., 1985. Моделирование на ЭВМ плоских турбулентных струй // Докл.АН СССР. Т. 282. -N 3. - С. 542-545.

12. Белоцерковский С.М., Котовский В.Н., Ништ М.И., Федоров P.M., 1988. Математическое моделирование плоскопараллельного отрывного обтекания тел. М.: Наука. - 232 с.

13. Белоцерковский С.М., Ништ М.И., 1978. Отрывное и безотрывное обтекание тонких крыльев идеальной жидкостью. М.: Наука. - 351 с.

14. Бетяев С.К., 1983. Отрывные течения. Новосибирск, (Препринт //АН СССР. Сиб. отд-ние. Ин-т теор. и прикл. механики; N 14-83).

15. Бондаренко А.В., Завьялов П.С., 1979. К объяснению пространственной формы вихревого жгута, образующегося в отсасывающих трубах гидротурбин // Сб. Гидравлические машины. Харьков. Вып. 13.-С. 123-126.

16. Брутян М.А., Крапивский П.Л., 1984. Движение системы вихревых колец в несжимаемой жидкости // Прикл. математика и механика. Т. 48, Вып. 3. - С. 503-506.

17. Бэтчелор Дж., 1973. Введение в динамику жидкости. М.: Мир. - 758 с.

18. Ван-Дайк М., 1967. Методы возмущений в механике жидкости. М.: Мир. - 311 с.

19. Ван-Дайк М., 1986. Альбом течений жидкости газа. М.: Мир. - 182 с.

20. Васильев О.Ф., 1958. Основы механики винтовых и циркуляционных потоков. -М. Л.: Госэнергоиздат. - 144 с.

21. Веретенцев А.Н., Гешев П.И., Куйбин П.А., Рудяк В.Я., 1989. О развитии метода вихревых частиц применительно к описанию отрывных течений // Журн. вычисл. математики и мат. физики. Т. 29, № 6. - С. 878-887.

22. Веретенцев А.Н., Куйбин П.А., Меркулов В.И., Рудяк В.Я., 19866. О выводе уравнений движения дискретных вихревых частиц для осееимметричных течений // Изв. СО АН СССР. Сер. техн. наук. № 10. - Вып. 2. - С. 45-50.

23. Веретенцев А.Н., Рудяк В.Я., 1986. Динамика завихренности в двумерных течениях невязкой жидкости. Новосибирск. - 41 с. - (Препринт / АН СССР. Сиб. отд-ние. Ин-т теоретической и прикладной механики; № 4).

24. Веретенцев А.Н., Рудяк В.Я., 1987а. О процессах образования и эволюции вихревых структур в сдвиговых слоях // Изв. АН СССР. Механика жидкости игаза. № 1.-С. 31-37.

25. Веретенцев А.Н., Рудяк В .Я., 19876. О взаимодействии внешнего акустического поля со слоем смешения // Моделирование в механике. Новосибирск. - Т. 1/18, №6.-С. 14-19.

26. Веретенцев А.Н., Рудяк В.Я., Яненко Н.Н., 1982. Вариационный метод построения дискретных вихревых моделей. Новосибирск. - 41 с. -(Препринт / АН СССР. Сиб. отд-ние. Ин-т теоретической и прикладной механики; № 29).

27. Веретенцев А.Н., Рудяк В.Я., Яненко Н.Н., 1985. Моделирование двумерных течений невязкой жидкости вихревыми частицами конечных размеров // Проблемы динамики вязкой жидкости. Новосибирск, - С.69-733.

28. Веретенцев А.Н., Рудяк В.Я., Яненко Н.Н., 1986а. О построении дискретных вихревых моделей течений идеальной несжимаемой жидкости // Журн. вычисл. математики и мат. физики. Т. 26, № 1. - С. 103-113.

29. Владимиров В.А., 1977. Формирование вихревых шнуров из восходящих потоков над испаряющейся жидкостью // Докл. АН СССР. Т. 236, № 2. - С. 316-318.

30. Владимиров В.А., Луговцов Б.А., Тарасов В.Ф., 1980. Подавление турбулентности в ядрах концентрированных вихрей // Журн. прикл. механики и техн. физики. № 5. - С. 69-76.

31. Владимиров B.C., 1976. Обобщенные функции в математической физике. М.: Наука. - 512 с.

32. Власов Е.В., Гиневский А.С., 1986. Когерентные структуры в турбулентных струях и следах // Итоги науки и техники. ВИНИТИ. Сер. Механика жидкости и газа. -Т. 20. -С. 3-84.

33. Власов Е.В., Гиневский А.С., Каравосов Р.К., Франкфурт М.О., 1982. Пристеночные пульсации давления в зоне отрыва за двумерными препятствиями // Труды ЦАГИ. Вып. 2137. - С. 3-22.

34. By Д. М., Вэкили А.Д., Ю Ф.М., 1989. Взаимодействие асимметричной струи с поперечным потоком // Аэрокосмическая техника. № 5. - С. 58-67.

35. Гахов Ф.Д., 1977. Краевые задачи. М.: Наука, - 364 с.

36. Журн. прикл. механики и техн. физики. -№ 1. С. 167-170. Григорьев Ю.Н., Левинский В.Б., 1986. Вариационная модель организованной завихренности в плоских течениях // Журн. прикл. механики и техн. физики. -№5.-С. 60-68.

37. Жуковский Н.Е., 1912. Вихревая теория гребного вихря // Труды отделенияфизических наук Общества любителей естествознания. Т. 16, № 1. Жуковский Н.Е., 1937а. Полное собрание сочинений. Т. 5. Вихри. Теория крыла.

38. Авиация. Под ред. А.П. Котельникова. M.-JL: Гл. ред. авиац. лит. - 484 с.

39. Жуковский Н.Е., 19376. Полное собрание сочинений. Т. 6. Винты. Ветряки. Аэродинамическая труба. Под ред. В.П. Ветчинкина. M.-JI.: Гл. ред. авиац. лит. - 424 с.

40. Зайцев Д.К., Смирнов Е.М., 1994. Об условиях возникновения возвратного закрученного течения в неизотермическом закрученном потоке газа // Процессы горения и охрана окружающей среды: Материалы II Всерос. конф. -Рыбинск. С. 43-48.

41. Залманзон JI.A., 1973. Аэродинамические методы измерения входных параметров автоматических систем. М.: Наука. - 464 с.

42. Зельдович Я.Б., 1947. О теории детонации в газах // Журн. техн. физ. Т. 17, Вып. З.-С. 3-26.

43. Зобнин А.И., 1986. Исследование структуры вихревого следа за профилем с угловой кромкой в начальной стадии отрывного обтекания // Гидродинамика подводного крыла. Новосибирск. - С. 71-84.

44. Ильичев К.П., Постоловский С.Н., 1972. Расчет нестационарного отрывного обтекания тел потоком невязкой жидкости // Изв. АН СССР. Механика жидкости и газа. № 2. - С. 72-82.

45. Кикнадзе Г.И., Краснов Ю.К., Подымака Н.Ф., Хабенский В.Б., 1986. Самоорганизация вихревых структур при обтекании водой полусферической лунки // Докл. АН СССР. Т. 291, № 6. - С. 1315-1318.

46. Кистович А.В., Чашечкин Ю.Д., 2000. Геометрия спиральных вихрей в однородной идеальной жидкости // Докл. РАН. Т. 372, № 1. - С. 46-49.

47. Кистович А.В., Чашечкин Ю.Д., 2001. Регулярный метод поиска дискретных симметрий моделей физических процессов // Докл. РАН. Т. 380, № 6. - С. 757-760.

48. Кныш Ю.А., 1988. О механизме переноса энергии в вихревой трубе пульсирующими крупными вихрями // Вихревой эффект и его применение в технике: Материалы V Всесоюзн. научно-техн. конф. Куйбышев. - С. 71-74.

49. Кныш Ю.А., Урывский А.Ф., 1981. Теория взаимодействия вторичных вихревых структур в закрученных потоках жидкости // Изв. вузов. Авиационнаятехника. № 3. - С. 53-58. Кузнецов В.И., 1989. Вязкость и ее влияние на эффект Ранка // Изв. СО АН

50. Лейбович С., 1985. Устойчивость и разрушение вихрей: Современное состояние и перспективы исследований // Аэрокосмическая техника. Т. 3, № 4. - С. 162-181.

51. Маслоу С.А., 1984. Неустойчивости и переход в сдвиговых течениях // Гидродинамические неустойчивости и переход к турбулентности. М.: Мир. -С. 218-270.

52. Меркулов А.П., 1969. Вихревой эффект и его применение в технике. М.: Машиностроение. - 183 с.

53. Меркулов В.И., 1989. Гидродинамика знакомая и незнакомая. М.: Наука, - 136 с.

54. Молчанов В.Ф., 1975. Некоторые вопросы расчета течений с тангенциальными разрывами // Уч. зап. ЦАГИ. Т. 4, № 4. - С. 19-25.

55. Мурахтина Т.О., Окулов B.JL, 2000. Изменение топологии и симметрии поля завихренности при турбулентном распаде вихря // Письма в ЖТФ. Т. 26, № 10.-С. 66-72.

56. Найфе А., 1974. Методы возмущений. М.: Мир, - 455 с.

57. Никольский А.А., 1957. О "второй" форме движения идеальной жидкости около обтекаемого тела (исследование вихревых отрывных потоков) // Докл. АН СССР. Т. 116, № 2. - С. 193-196.

58. Никольский А.А., 1970. Законы подобия для трехмерного стационарного обтекания тел жидкостью и газом // Уч. зап. ЦАГИ. Т. 1. - N 1. - С. 1-7.

59. Окулов В.Л., 1993. Резонансные гидроакустические процессы в проточной части машин и агрегатов с интенсивной закруткой потока // Автореф. дисс. доктора, физ.-мат. наук. Новосибирск. - 34 с.

60. Окулов В.Л., Мартемьянов С.А., 2001. Эффект снижения массообменных процессов в закрученных потоках // Письма в ЖТФ. Т. 27, № 18. - С. 22-28.

61. Окулов В.Л., Соренсен Ж.Н., Войгт Л.К., 2002. Чередование право- и левовинтовых вихревых структур при увеличении интенсивности закрутки потока в цилиндрической каверне с вращающимися торцами // Письма в ЖТФ. Т. 28, № 2. - С. 37-44.

62. Пази Л.Г., 1968. Исследование периодических пульсаций давления в отсасывающих трубах гидротурбин // Тр. ВНИИ Гидромаш. Вып. XXXVII. -С. 151-166.

63. Пиралишвили Ш.А., Поляев В.М., Сергеев М.Н., 2000. Вихревой эффект.

64. Эксперимент, теория, технические решения / Под ред. А.И. Леонтьева. М.: УНПЦ «Энергомаш». -415 с.

65. Поликовский В.И., Перельман Р.Г., 1959. Воронкообразование в жидкости с открытой поверхностью. М.; - Л.: Госэнергоиздат, - 191 с.

66. Прудников А.П., Брычков Ю.А., Маричев О.И., 1981. Интегралы и ряды. Элементарные функции. М.: Наука.

67. Рабинович М.И., Сущик М.М., 1990. Регулярная и хаотическая динамика структур в течениях жидкости // Успехи физических наук. Т. 160, Вып. 1. - С. 3-64.

68. Сабельников В.А., Смирных Е.А., 1985. Численный расчет турбулентного течения на начальном участке плоского канала с острыми кромками методом дискретных вихрей // Уч. зап. ЦАГИ. Т. 16, № 4. - С. 59-64.

69. Седов Л.И., 1966. Плоские задачи гидродинамики и аэродинамики. М.: Наука, -287 с.ч

70. Седов Л.И., 1970. Механика сплошной среды. -М.: Наука. Т. 1. - 492 с.

71. Смульский И.И., 1992. Аэродинамика и процессы в вихревых камерах / Под ред. И.Р. Шрейбера. Новосибирск: ВО «Наука». Сибирская издательская фирма -301 с.

72. Сноу Д., 1984. Торнадо // В мире науки. С. 44-55.

73. Сычев В.В., 1987. Асимптотическая теория отрывных течений. -М.: Наука, 255 с.

74. Сэффмэн Ф.Дж., 2000. Динамика вихрей. М.: Научный Мир. - 376 с.

75. Уильяме Дж., 1979. Отрыв пограничного слоя несжимаемой жидкости // Механика (Новое в зарубежной науке, N 21). Вихревые движения жидкости. -М.: Мир.

76. Халатов А.А., 1989. Теория и практика закрученных потоков. Киев: Наукова думка.-180 с.

77. Хитрин Л.Н., 1957. Физика горения и взрыва. М.: МГУ.

78. Чжен П., 1972-1973. Отрывные течения. Т. 1 3. -М.: Мир,

79. Шторк С.И., 1994. Экспериментальное исследование вихревых структур в тангенциальных камерах // Дисс. канд. физ.-мат. наук. Новосибирск: ИТ СО РАН.-156 с.

80. Штым А.Н., 1985. Аэродинамика циклонно-вихревых камер. Владивосток:1. Дальневост. ун-т. 199 с.

81. Acarlal M.S., Smith C.R., 1987. A study of hairpin vortices in a laminar boundary layers // J. Fluid Mech. Vol. 175. - P. 1-83.

82. Acton E, 1980. A modelling of large eddies in an axisymmetric jet // J. Fluid Mech. -Vol. 98, № 1.-P. 1-31.

83. Acton E., 1976. The modelling of large eddies in a two-dimensional shear layer // J. Fluid Mech. Vol. 76. - P. 561-592.

84. Adebiyi A., 1981. On the existence of steady helical vortex tubes of small cross-section // Q. J. Mech. Appl. Math.-No. 34.-P. 153-177.

85. Albright L.F., Alexander L.G., 1957. Stabilization of flame by cyclone gas flow // Proc. VI Symp. on Combustion. N.Y. P. 464-472.

86. Alekseenko S.V., Kuibin P.A., Okulov V.L., Shtork S.I., 1999. Helical vortices in swirl flow // J. Fluid Mech. Vol. 382. - P. 195-243.

87. Alekseenko S.V., Shtork S.I., 1992. Swirling flow large-scale structures in a combustor model // Russ. J. Eng. Thermophys. Vol. 2, № 4. - P. 231-266.

88. Amitani Т., Adachi Т., Kato T.A., 1983. Study on Temperature Separation in a Large Vortex Tube // T. JSME. Vol. 49. - P. 877-884.

89. Arms R.J., Hama F.R., 1965. Localized-induction concept on a curved vortex and motion of an elliptic vortex ring // Phys. Fluids. Vol. 8. - P. 553-559.

90. Batchelor G.K., 1956. On steady laminar flow with closed streamlines at large Reynolds number // J. Fluid Mech. — Vol. 1. N 2. - P. 177-189.

91. Batchelor G.K., 1964. Axial flow in trailing line vortices // J. Fluid Mech. Vol. 20. -P. 645-658.

92. Beale J.Т., Maida A., 1985. High order accuracy vortex method with explicit velocity kernels // J. Comput. Phys.- Vol. 58. N 2. - P. 188-208.

93. Birkhoff G., Fisher J., 1959. Do vortex sheets roll up? //Rend. Circ. Mat. Palermo. -Ser. 2. Vol. 8. -P. 77-90.

94. Boersma J., Wood D.H., 1999. On the self-induced motion of a helical vortex // J. Fluid Mech. Vol. 384. - P. 263-280.

95. Boersma J., Yakubovich S.B., 1998. Solution to problem 97-18*: the asymptotic sum of a Kapteyn series // SIAM Rev. Vol. 40. - P. 986-990.

96. Borissov A.A., Kuibin P.A., Okulov V.L., 1993. Flame shapes in swirl flow // Russ. J.

97. Cantwell B.J., 1981. Organized motion in turbulent flow // Annu. Rev. Fluid Mech.

98. Vol. 13.-P. 457-515. Cassidy J.J., Falvey H.T., 1970. Observations of unsteady flow arising after vortexbreakdown // J. Fluid Mech. Vol. 41. - P. 727-736. Chanaud R.C., 1965. Observations of oscillatory in certain swirling flows // J. Fluid.

99. Mech.-Vol. 21.-P. 111-127. Chandrsuda C., Mehta R.D., Weir A.D., Bradshaw P., 1978. Effect of free-stream turbulence on large structure in turbulence mixing layers // J. Fluid. Mech. Vol. 85.-P. 693-704.

100. Chorin A.J., 1973. Numerical study of slightly viscous flow// J. Fluid Mech. Vol. 57. -P. 785-798.

101. Chorin A.J., Bernard P.S., 1973. Discretization of vortex sheet with an example of rollup // J. Comput. Phys. Vol. 13. - № 3. - P. 423^129.

102. Da Rios L.S., 1906. On the motion of an unbounded fluid with a vortex filament of tany shape (in italian) // Rend. Circ. Mat. Palermo. Vol. 22. - P. 117-135.

103. Dritschel D.G., 1991. Generalized helical Beltrami flows in hydrodynamics and magnetohydrodynamics // J. Fluid. Mech. Vol. 222. - P. 525-541.

104. Escudier M.P., 1984. Observations of the flow produced in a cylindrical container by a rotating endwall // Exps. Fluids. Vol. 2, № 4. - P. 189-196.

105. Escudier M.P., 1988 Vortex breakdown: observations and explanations // Prog. Aerospace Sci. Vol. 25. - P. 189-229.

106. Escudier M.P., Bornstein J., Maxworthy Т., 1982. The dynamics of confined vortices у

107. Proc. Roy. Soc. bond. A382. - P. 335-360.

108. Escudier M.P., Bornstein J., Zehnder N., 1980. Observations and LDA measurements of confined vortex flow // J. Fluid Mech. Vol. 98. - P. 49-63.

109. Evans R.A., Bloor M.I.G., 1977. The starting mechanism of wave-induced flow through a sharp-edged orifice // J. Fluid Mech. Vol. 82. - P. 115-128.

110. Faler J.H., Leibovich S., 1977. Disrupted states of vortex flow and vortex breakdown // Phys. Fluids. Vol. 20. - P. 1385-1400.

111. Fanelli M., 1989. The vortex rope in the draft tube of Francis turbines operating at partial load: a proposal for a mathematical model // J. Hydraulic Research. Vol. 27, № 6. - P. 769-807.

112. Feikema D., Chen R.-H., Driscol. J.F., 1990. Enhancement of flame blow-out limits by the use of swirl//Comb. & Flame.-Vol. 80.-P. 183-195.

113. Fiedler H.E., Fernholz H.H., 1990. On management and control of turbulent shear flows // Progr. Aerospace Sci. Vol. 27, № 4. - P. 305-387.

114. Freymuth P., Bank W., Palmer M., 1983. Visualization of accelerating flow around anairfoil at high angles of attack // Z. Flugwiss. Weltraumforsch. Bd. 7, Heft 6. - S. ^392.400.

115. Garg A.K., Leibovich S., 1979. Spectral characteristics of vortex breakdown flow fields // Phys. Fluids. Vol. 22. - P. 2053-2064.

116. Ghoniem A.F., Ng K.K., 1987. Numerical study of the dynamics of a forced shearlayer//Phys. Fluids. Vol. 30.-N 3.-P. 706-721. '

117. Guarga R.F., Rodal E.C., Solorio A., 1985b. Oscillations in confined swirling flow and its relation with the axial pressure gradient // Ibid. P. 14/1-14/9.

118. Guarga R.F., Sanchez A.H., Gracia J.S., Cafaggi A.F., 1985c. Oscillatory characteristics of swirling, confined, turbulent and non cavitating-flows // Ibid., P. 13/1-13/13.

119. Guarga R.F., Solorio A., Rodal E.C., 1985d. Analysis of some oscillatory phenomena associated with swirling and confined flows // Ibid., P. 15/1-15/12.

120. Hald O.H., 1979. Convergence of vortex methods for Euler's equations. II // SIAM J. Numer. Anal.-Vol. 16.-N 5. P. 726-755.

121. Hama F.R., 1962. Progressive deformation of a curved vortex filament by its own induction // Phys. Fluids. Vol. 5. - P. 1156-1162.

122. Hardin J.C., 1982. The velocity field induced by a helical vortex filament // Phys. Fluids. Vol. 25, № 11.-P. 1949-1952.

123. Hashimoto H., 1971. Theoretical study of swirling flow accompanied by helical cavity core in circular pipe // Rep. Inst. High Speed Mech. Vol. 23, № 228. - P. 61-97.

124. Hicks W.M., 1899. Researches in vortex motion-Part III. On spiral or gyrostatic vortex aggregates // Phil. Trans. Roy. Soc. Lond. A192. - P. 33-101.

125. Но C.M., Huerre P., 1984. Perturbed free shear layers //Ann. Rev. Fluid Mech. Palo Alto. Calif., - Vol. 16. - P. 365-424.

126. Hopfinger E.J., van Heijst G.J.F., 1993. Vortices in rotating fluids // Annu. Rev. Fluid Mech. Vol. 25. - P. 241-289.

127. Flames // Combustion and Flame. Vol. 75. - P. 367-379.

128. Kelly R.E., 1968. On the resonant interaction of neutral disturbances in two inviscid shear flows // J. Fluid Mech. Vol. 31. - P. 789-799.

129. Kelvin, Lord, 1880. Vibrations of a columnar vortex // Philos. Mag. Vol. 10. - P. 155-168.

130. Kida S., 1981. A vortex filament moving without change of form // J. Fluid Mech. -Vol. 112.-P. 397-409.

131. Kim H.T., Kline S.J., Reynolds W.C., 1971. The production of turbulence near a smooth wall in a turbulent boundary layer // J. Fluid Mech. Vol. 50. - P. 133-160.

132. Kuibin, P.A., 2002. On motion of a double helical vortex in a cylindrical tube // Proc. IUTAM Symp. "Tubes, sheets and singularities in fluid dynamics". Kluwer. - P. 81-86.

133. Kuibin P.A., 2003. On structure of vorticity field in the spherical Hicks vortex // J. Eng. Thermophys. Vol. 12, №. 3. - P. 22-29.

134. Kuibin P.A., Okulov V.L., 1998. Self-induced motion and asymptotic expansion of the velocity field in the vicinity of helical vortex filament // Phys. Fluids. Vol. 10, № 3.-P. 607-614.

135. Kuibin P.A., Rudyak V.Ya., Veretentsev A.N., 1991. Processes of the instability development in separated flows behind the plate // Separated flows and jets: Proc. IUTAM Symp., Novosibirsk, USSR, July, 9-13 1990. Berlin: Springer-Verlag. -P. 747-750.

136. Kumar R., Conover Т., 1993. Flow visualization studies of a swirling flow in a cylinder // Exp. Therm, and Fluid Sci. Vol. 7, No. 3. - P. 254-262.

137. Kurosaka M., 1982. Acoustic streaming in swirling flow and Rangue-Hilsch (vortex tube) effect // J. Fluid Mech. Vol. 124. - P. 139-172.

138. Palo Alto. Calif., Vol. 15. -P. 223-239. Liu J.T.C., 1989. Coherent structures in transitional and turbulent free shear flows //

139. Astronautica. ACTA. Vol. 17. - P. 363-374. Maxworthy, Т., 1977. Some experimental studies of vortex rings // J. Fluid Mech.

140. Vol. 81, Part 3.-P. 465—495 Maxworthy Т., Hopfinger E.J., Redekopp L.G., 1985. Wave motions on vortex cores //

141. Archiv. -Bd. 33, Heft l.-S. 1-23.

142. Milinazzo F., Saffman P.G., 1987. The calculation of large Reynolds number two-dimensional flow using discrete vortices with random walk // J. Comput. Phys. -Vol. 23,-N4.-P. 380-392.

143. Moffatt H.K., 1969. The degree of knottedness of tangled vortex lines // J. Fluid Mech. -Vol. 35.-P. 117-129.

144. Mollenkopf G., Raabe J., 1970. Measurements of velocity and pressure in the draft tube of a Francis turbine. B3. // Presented at 5-th IAHR Symposium, Stockholm, Sweden.

145. Moore D.W., 1974. A numerical study of the roll-up of a finite vortex sheet // J. Fluid Mech. Vol. 63,- P. 225-235.

146. Moore D.W., Saffman P.G., 1972. The motion of a vortex filament with axial flow // Phil. Trans. R. Soc. Lond. A272. - P. 403-429.

147. Murakami M., 1961. Vibration of water-turbine draft tubes // J. Eng. for Power (Trans. ASME, Ser. A.).-Vol. 83, № l.-P. 36^12.

148. Nemirovskii S.K., 1998a. Gaussian model of vortex tangle in superfluid turbulent Hell // Phys. Rev. B1. Vol. 57, № 10. - P. 5972-5985.

149. Nemirovskii S.K., 1998b. Superfluid mass current induced by chaotic vortex lines in Hell // Ibid. P. 5986-5994.

150. Nishioka M. Miyagi Т., 1978. Measurements of velocity distribution in the laminar wake of a flat plate // J. Fluid Mech. Vol. 84. - P. 705-715.

151. Okulov V.L., 1995. The velocity field induced by vortex filaments with cylindrical and conic supporting surface // Russ. J. Eng. Thermophys. Vol. 5, № 2. - P. 63-75.

152. Perry A.E., Lim T.T., 1978. Coherent structures in coflowing jets and wakes // J. Fluid. Mech. Vol. 88, № 3. - P. 451-463.

153. Peyne F.M., Ng T.T., Nelson R.S., 1988. Visualization and wake surveys of vortices flow over a delta wing // AIAA J. Vol. 26, No. 2. - P. 137-143.

154. Pierce D., 1961. Photographic evidence of the formation and growth of vorticity behind plates accelerated from rest in still air // J. Fluid Mech. Vol. 11. - P. 460464.

155. Pierrehumbert R.T., 1980. A family of steady, translating vortex pairs with distributedvorticity // J. Fluid Mech. Vol. 99. - P. 129-144.

156. Prandtle L., 1924. Uber die Entstehung von Wirbeln in der idealen Flussigkeit // Vortrage aus Hydro- und Aerodynamic. Berlin.

157. Pullin D.I., 1978. The large-scale structure of unsteady self-similar rolled-up vortex sheets // J. Fluid Mech. Vol. 88. - P. 401-430.

158. Pullin P.I., Perry A.E., 1980. Some flow visualization experiments on the starting vortex // J. Fluid Mech. Vol. 97. - P. 239-255.

159. Ricca R.L., 1994. The effect of torsion on the motion of a helical vortex filament // J. Fluid Mech. Vol. 273. - P. 241-259.

160. Ricca R.L., 1996. The contributions of Da Rios and Levi-Civita to asymptotic potential theory and vortex filament dynamics // Fluid Dynamics Research. Vol. 18.-P. 245-268.

161. Roberts K.V., Christiansen J.P., 1972. Topics in computational fluid michanics // Сотр. Phys. Commun. Vol. 3. - P. 14-32.

162. Rosenhead L., 1931. The formation of vortices from a surface of discontinuity // Proc. R. Soc. Lond. Vol.A134. - P. 170-192.

163. Rott N., 1956. Difraction of weak shock wave with vortex generation // J. Fluid Mech. -Vol. l.-P. 111-129.

164. Sakai Y., Ishizuka S., 1991. Structures of Tubular Flames // JSME, Ser. II. Vol. 34, No2.-P. 234-241.

165. Sato H., Kuriki K., 1961. The mechanism of transition in the wake of a thin flat plate placed parallel to a uniform flow // J. Fluid Mech. Vol. 11. - P. 321-342.

166. Scully M.P., 1975. Computation of helicopter rotor wake geometry and its influence on rotor harmonic airloads // Massachusetts Inst, of Technology Pub. ARSL TR 178-1. Cambridge.

167. Shimizu S., 1986. Discrete-vortex simulation of a two-dimensional turbulent impinging jet // Bulletin of JSME. Vol. 29, № 254. - P. 2440-2446.

168. Sozou C.J., Swithenbank J., 1969. Adiabatic transverse waves in a rotating fluid // J. Fluid. Mech. Vol. 38. - P. 657-671.

169. Stephan K., Lin S., Durst M. et al., 1983. An investigation of energy separation in a vortex tube // Int. J. Heat Mass Transfer. Vol. 26. - P. 341-348.

170. Takaki R., Hussain A.K.F.M., 1984. Dynamics of entangled vortex filaments // Phys.

171. Tung C., Ting L., 1967. Motion and decay of a vortex ring // Phys. Fluids. Vol. 10, № 5.-P. 901-910.

172. Ward G.N., 1955. Linearized theory of steady high speed flow. Cambridge. -Widnall S.E., 1972. The stability of helical vortex filament // J. Fluid Mech. - Vol. 54. -P. 641-663.