Авторефераты по всем темам  >>  Авторефераты по разным специальностям

На правах рукописи

Порфирьев Денис Петрович

Моделирование структуры дозвуковых закрученных потоков в присутствии локализованных источников тепловыделения

01.02.05 - механика жидкости, газа и плазмы

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Казань - 2012

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

Официальные оппоненты:

Зарипов Шамиль Хузеевич, доктор физико-математических наук, профессор, федеральное государственное автономное образовательное учреждение Казанский (Приволжский) федеральный университет, заведующий кафедрой Моделирование экологических систем;

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

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

Защита состоится л27 декабря 2012 г. в 16.00 на заседании диссертационного совета Д 212.081.11 при ФГАОУ ВПО Казанский (Приволжский) федеральный университет по адресу: 420008, г. Казань, ул.

Кремлевская, 18, ауд. мех. 2.

С диссертацией можно ознакомиться в научной библиотеке им. Н.И.

обачевского ФГАОУ ВПО Казанский (Приволжский) федеральный университет по адресу: 420008, г. Казань, ул. Кремлевская, 18.

Автореферат разослан 23 ноября 2012 г.

Ученый секретарь диссертационного совета, к. ф-м.н., доц. Саченков А.А.

Общая характеристика работы

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

В тоже время, имеющиеся данные о совместном влиянии закрутки, турбулентности, теплообмена и разрядной плазмы на структуру течения весьма неполны. Существующие работы в этой области (В.Т.Волов, В.М.Шмелев, 1982; В.Т.Волов, Х.Д.Ламажапов, 1988; А.В.Казаков и др., 1998, 1999, 2003, 2005; С.Т.Суржиков, 2002, 2008; Е.С.Асмолов и др., 2005; Jeong H. Kim et al., 2003; A,Klimov, V.Bityurin et al., 2008; Моралев И.А., 2010) показывают, что наличие источников тепловыделения может заметно менять условия развития неустойчивости течений и области существования как различных типов течений, так и разрядных структур в вихревом потоке. Усложняет проблему тот факт, что с целью интенсификации рабочих процессов вихревые аппараты обычно эксплуатируются на режимах с повышенными массовыми расходами и сильной закруткой потока (Ю.А.Кныш, 1977, 1981, 1982, 1984; А.П.Меркулов, 1969; С.В.Лукачев 1977, 1981, 1984; В.К.Щукин, А.А.Халатов 1982;

И.И.Смульский 1992; А.Гупта, Д.Лилли, Н.Сайред, 1987; Ш.А.Пиралишвили, 2000, 2005; С.В.Алексеенко, П.А.Куйбин, В.Л.Окулов, 2003; С.И.Шторк и др., 2006; В.К.Ахметов, В.Я.Шкадов 2003, 2008, 2009; В.Т.Волов, 2006, 2011). При этом течение оказывается существенно трехмерным и, как правило, нестационарным. Поэтому комплексное рассмотрение данных факторов на процессы переноса и горения разряда в турбулентных закрученных потоках представляет весьма сложную и, в связи с практическими потребностями, актуальную задачу. Существует небольшое число работ, посвящнных влиянию плазменных зон на вихревые структуры и демонстрирующих заметную перестройку вихревых структур при относительно малых энерговкладах.

Детальные экспериментальные исследования вихревого потока с зонами тепловыделения, создаваемыми разрядами различного типа в камере с открытым концом были проделаны в ОИВТ РАН. В этих экспериментах наблюдались скачкообразные качественные перестройки структуры тлеющего и ВЧЕ-разряда в вихревых потоках (альфа-гамма переходы, контракция разряда, переход от филаментарного стримерного в режим стабильного протяжнного плазмоида и обратно, ряд других), а также параметров потока, в частности, образование зон противототока или, наоборот, сужения их вплоть до полного исчезновения. Наблюдалась сильная перестройка вихревого потока и светящихся областей, создаваемых тлеющим и ВЧЕ-разрядом в вихревой камере с открытым концом. Зафиксированные существенные изменения аксиальных и тангенциальных скоростей в разных разрядных условиях подтверждают гипотезу, что неравновесная плазма способна усиливать и разрушать вихревые образования. Эксперименты продемонстрировали качественное изменение структуры вихря в присутствие тепловыделения в области разрядной плазмы. В немногочисленных теоретических исследованиях показано, что подвод энергии в систему может привести к слабому увеличению угловой скорости вихря, наличие периодической накачки энергии может приводить к увеличению завихренности (Н.А.Винниченко, 2009).

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

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

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

В соответствии с поставленной целью определены основные задачи диссертации:

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

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

- Изучить влияние источников тепловыделения на устойчивость стационарных вихревых профилей и закрученных потоков различного типа.

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

Научная новизна полученных результатов состоит в следующем:

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

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

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

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

Основные положения, выносимые на защиту:

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

2. Результаты моделирования газодинамических условий возникновения переходов между коронной и шнуровой формами одноэлектродных ВЧЕ - разрядов. Механизм формирования зон свечения от разрядных областей в вихревых потоках в воздухе при атмосферном давлении.

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

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

Связь с государственными программами.

Работы по теме диссертации выполнялись в соответствии с планами фундаментальных научно-исследовательских работ по программам:

аналитическая ведомственная целевая программа Развитие научного потенциала высшей школы (2009-2010 годы), проекты 1.1.11, 1.2.08, 2.1.1/309, 2.1.1/13492 грантами ФЦП Научные и научно-педагогические кадры инновационной России на 2009-2013 гг. №№ П2315, П2450, 14.740.11.0999, 14.В37.21.0767, государственного задания Минобрнауки РФ на 2012-2014 гг, № 2.560.2011.

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

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

Апробация результатов работы. Основные результаты диссертации докладывались и обсуждались на семинарах Казанского (Приволжского) федерального университета, Самарского государственного аэрокосмического университета им. С.П. Королева, Самарского филиала физического института им. П.Н.Лебедева РАН и были представлены на следующих Всероссийских и Международных конференциях: 6th - 11th International Workshop on MagnetoPlasma Aerodynamics (Москва, 2007 - 2012), V - VIII Международных междисциплинарных научных конференциях Курдюмовские чтения. Идеи синергетики в естественных науках (Тверь, 2009-2012), II и III Международных конференциях по математической физике и ее приложениям (Самара, 2010, 2012), Х Королвских чтений (Самара, 2009), Международной конференции с элементами научной школы для молодежи Перспективные информационные технологии для авиации и космоса (ПИТ-2010), (Самара, 2010), Международной конференции Физика высокочастотных разрядов, (Казань, 2011), 9th Euromech Fluid Mechanics Conference (EFMC9) (Roma, Italy, 2012). Результаты, вошедшие в диссертацию, отмечены премией победителя конкурса Молодой учный Самарской области 2011 г (номинация аспирант, Механика).

Основные публикации. По теме диссертации опубликовано 14 работ, в том числе 6 статей в изданиях, рекомендованных ВАК РФ, 8 трудов Международных и Всероссийских конференций.

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

Структура и объем диссертации. Диссертация состоит из введения, трех глав, заключения, списка литературы (150 наименований), изложена на 1страницах, содержит 57 рисунков.

Основное содержание работы

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

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

В разделе 1.1 проведена постановка задачи. При моделировании использовались уравнения Навье-Стокса, записанные в цилиндрической системе координат и дополненные уравнением энергии. Поставлены граничные условия, проведено обоснование выбранной модели турбулентности СпалартаАллмараса (P.R.Spalart, M.L.Shur, 1997). С помощью программы Gambit 2.4.создана расчтная геометрия завихрителя и трубы (рисунок 1).

Геометрия завихрителя моделировалась в соответствии с экспериментальными данными. Число элементов расчетной сетки менялось от 80 до 240 тысяч, сетка сгущалась вблизи стенок трубы для устойчивого счета течений газа с большими градиентами и минимизации влияния численной диффузии.

Рисунок 1. Геометрия завихрителя и трубы В разделе 1.2 исследовались условия существования прямых потоков и потоков с рециркуляционными зонами, а также найдены пространственные распределения основных параметров потока.

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

Проведена классификация режимов течения на плоскости (Qt, Qax).

Соответствующая параметрическая диаграмма приведена на рисунке 2.

Рисунок 2. Теоретически рассчитанные области существования различных типов закрученного течения в трубе. Зоны 1, 2 и 3 - зоны развитого противотока, промежуточная зона и зона прямого течения соответственно. Штриховые линии - верхняя (U) и нижняя (D) границы промежуточной зоны. Радиус трубы R ~ 1.9 см, длина L = 40 см В области 1, соответствующей большим степеням закрутки потока, в приосевой зоне существует развитый противоток, возникающий у открытого конца трубки и занимающий значительную часть е длины, а в области 3, соответствующей слабо закрученным потокам, имеет место режим прямого течения. В области 2 между кривыми U и D имеет место промежуточный режим.

Получено, что переход от режима прямого течения к течению с противотоками является мягким и является непрерывным при переходе параметра Z = Qt/(Qax + Qt) при его увеличении через область критических значений, что соответствует известным результатам (А.Гупта, Д.Лилли, Н.Сайред, 1987; С.И.Шторк и др., 2006). Показано, что практически во всей рассматриваемой области параметров зона противотока возникает у открытого конца трубки. При этом пристеночная зона прямого течения и приосевая зона противотока довольно резко отделены друг от друга узкой застойной зоной с малой аксиальной компонентой скорости течения (рисунок 3). В области малых массовых расходов Qt и Qax, кроме зоны противотока, расположенной у выхода, может возникать область противотока в окрестности завихрителя.

Рисунок 3. Радиальные распределения аксиальной скорости течения (z = 33 см) при Qt= 4 г/с, Qax= 0 г/с, N = 0 Вт. Пунктирная линия Ч режим ламинарного течения, сплошная - режим турбулентного течения.

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

Расчеты показали, что в трубках, где скорости течения достаточно велики (Re > 104), присутствие источников нагрева, с мощностями N 1 кВт слабо меняет условия возникновения рециркуляционных зон и диаграмма режимов потока практически совпадает с приведенной на рисунке 2. Показано, что рассчитанные профили качественно совпадают с профилями компонент тангенциальной и аксиальной скоростей для частично закрученных потоков, полученными экспериментально при наличии тлеющих и ВЧЕ - разрядов в воздухе при атмосферных давлениях со средней мощностью N = 0.2 - 1 кВт (А.И.Климов, И.А.Моралев и др., 2010). При наличии локализованного источника тепловыделения в трубке возникают приосевые области высокой температуры. Если осевой и тангенциальный расходы газа соответствуют зоне 1, то нагретая зона формируется в области, расположенной вверх по потоку, а ее длина полностью определяется длиной зоны противотока, если осевой и тангенциальный расходы газа соответствуют зоне 3, то нагретый воздух сносится потоком в сторону открытого конца трубки, рисунок 4. Если осевой и тангенциальный расходы газа находятся в зоне 2, то структура температурного поля определяется соотношением между длиной зоны противотока и расположением поперечно-локализованного источника тепловыделения.

а б Рисунок 4. Изолинии температуры закрученного потока в воздухе, массовый расход:

а) 2 г/с, б) 9 г/с. Поток распространяется слева направо. Положение источника указано решеткой При малых скоростях течения (Re < 104) наличие источника нагрева может приводить к разрушению зоны противотока. При наличии локализованного источника тепловыделения приосевая зона пониженных скоростей сильно перегревается и плотность в ней резко падает. В результате, на границе с этой областью также резко растт радиальная составляющая прямого потока, направленная к оси трубы по градиенту температуры vr1 / vr0 0 / 1 T1 / T0, где 0,T0 и 1,T1 - плотности и температуры газа в холодной и нагретой областях, соответственно. Следствием этого является преломление линий тока к оси трубки. В результате проникновения закрученных периферийных областей потока в осевую область тангенциальная скорость вблизи оси трубы увеличивается по сравнению со случаем без источника. Осевая скорость также растт. Эффективность конвективного охлаждения растт с ростом расхода газа. При малых расходах газа (или больших энерговкладах) приосевая зона нагревается сильнее. В результате преломление линий тока происходит под большим углом, и основной поток протекает через локальную зону источника к выходному отверстию, пробивая и разрушая зону противотока. На аналогичный эффект влияния разрядной области на закрученный поток ранее было указано в работе (Ю.П.Райзер, 1980).

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

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

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

В соответствии с результатами, полученными в разделе 1.3, при условиях, соответствующих зоне 1, происходит формирование расположенной вверх по потоку от электрода нагретой приосевой зоны. В этой зоне изменяется баланс заряженных частиц, который определяется конкуренцией процессов ионизации и диссоциативной рекомбинации, как в электроположительных газах, а не трехтельным прилипанием (Ю.П.Райзер, М.Н.Шнейдер, И.А.Яценко, 1995). Это приводит к сильному падению погонного сопротивления канала и формированию визуально однородной волны ионизации (А.И.Климов, И.А.Моралев и др., 2010). Волна ионизации останавливается в неоднородном поле на некотором расстоянии от электрода и формирует шнуровой разряд (рисунок 5).

Рисунок 5. Теоретически рассчитанные области существования различных типов закрученного течения в трубе и области существования шнурового (выше заштрихованной зоны) и филаментарного (ниже заштрихованной зоны) ВЧЕ-разряда (Климов А.И., Моралев И.А. и др., 2010). Зоны 1, 2 и 3 соответствуют рисунку 2. Радиус трубы R = 1.9 см, длина L = 40 см, мощность источника N = 700 Вт Проведена оценка длины шнура и показано хорошее соответствие экспериментальным данным. В области прямого течения поток тепла направлен в сторону открытого конца трубки. Конвективный прогрев разрядной области и соответствующее изменение ионизационного баланса не происходит. Поэтому условия формирования разряда в зоне, расположенной вверх по потоку от электрода, слабо отличаются от имеющих место в отсутствие потока и формируется высокочастотная кистевая корона. В промежуточном диапазоне значений Z одновременно существуют обе формы разряда.

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

В разделе 2.2 получено и проанализировано дисперсионное соотношение для малых возмущений модельного стационарного невязкого несжимаемого потока с распределением тангенциальной скорости, соответствующей вихрю Ренкина и различными кусочно-непрерывными распределениями аксиальной скорости потока и плотности газа в радиально ограниченной области. Анализ устойчивости с учетом скачка плотности газа проведен при условиях, допускающих существование приосевой зоны противотока. В предельном случае однородного потока дисперсионное соотношение переходит в полученные в (T. Loiseleux, J. M. Chomaz, P. Huerre, 1998; F. Gallaire, J.-M.

Chomaz, 2003), а для закрученного потока типа следа со спутным потоком и скачком плотности - в полученные в (M.S.Uberoi, C.Y. Chow, J.P. Narain, 1972;

D.W.Lim, L.G.Redekopp, 1998).

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

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

В третьей главе проведено численное исследование структуры нестационарного трехмерного вихревого течения в дозвуковом потоке в круглой трубе с объемным подводом энергии. В разделе 3.1 проведн обзор результатов работ, посвященных исследованию крупномасштабных вихревых структур в закрученных потоках. В разделе 3.2 приведены физическая и математическая постановка задачи. Геометрия завихрителя и трубки соответствовали приведенным в разделе 1.1. Число элементов расчетной сетки менялось 2 до 6 млн. У стенок трубы для уменьшения влияния численной диффузии при помощи алгоритма Map была сгенерирована структурированная сетка с минимальным радиальным размером ячейки 0.05 мм. В центральной области использовался алгоритм Pave. В обоих случаях проводился контроль качества сетки. Расчеты показали, что сетка была выполнена с запасом точности: уже при достижении количества ячеек в 3 млн. дальнейшее ее измельчение не приводило к изменению результатов расчетов более чем на 12%. Для расчета конвективных слагаемых применялась противопоточная разностная схема QUICK. Сходимость контролировалась как стандартным мониторингом невязок, так и проверкой массового расхода газа. Проверялось также постоянство массового расхода в зависимости от времени: на выбранных промежутках расчта оно достигалось с точность до ~10 -7г/с.

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

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

Показано, что течение, соответствующее области 3 на рисунках 2 и 5 при переходе параметром Z через область критических значений становится нестационарным и характеризуется возникновением так называемого прецессирующего вихревого ядра (ПВЯ), что соответствует известным результатам (С.В. Алексеенко, П.А. Куйбин, В.Л. Окулов, 2003; С.И.Шторк и др., 2006). Исследована зависимость частоты прецессии от массового расхода газа, параметра кручения и мощности источника нагрева. Частота прецессии определялась путем построения спектра осевой и тангенциальной скоростей, а также плотности и давления. Для этого снималась зависимость параметров потока от времени в точках, расположенных на разном удалении от оси системы. К полученные данным применялось быстрое преобразование Фурье, реализованное в пакете Matlab, и строились спектры колебаний, по которым можно судить о частоте прецессии. Характерный спектр колебаний давления приведен на рисунке 6. Обработка данных численного эксперимента показала, что частоте прецессии соответствуют спектральные пики с максимальной амплитудой и минимальной частотой. Для определения номера азимутальной моды, неустойчивость которой приводит к развитию ПВЯ, был построен спектр колебаний тангенциальной скорости течения на оси трубки.

Рисунок 6. Спектр колебаний статического давления при значении массового расхода Qt = 4 г/с, Qax = 1 г/с, мощности 100 Вт Этот результат показывает, что наблюдаемая прецессия вихревого ядра связана с развитием неустойчивости моды с азимутальным числом |m| = 1, так как только эта мода возмущает течение на оси потока, а остальные, в силу симметрии, не возмущают (С.И.Шторк и др., 2005, 2006). Кроме того, расчеты показывают, что направление вращения ядра совпадает с направлением вращения основного потока, но ось винта закручена против вращения. Это соответствует левовинтовым вихрям m = 1. Расчеты структуры потока при фиксированных мощности источника и параметре крутки показывают, что частота прецессии растет с ростом массового расхода газа (рисунок 7).

Рисунок 7. Зависимость частоты прецессии от суммарного расхода газа Q при мощности источника N = 700 Вт и расходе газа через аксиальный подвод Qax =0 г/с Этот результат качественно совпадает с известными результатами как экспериментальных, так и теоретических исследований (А.Гупта, Д.Лилли, Н.Сайред, 1987; С.В. Алексеенко, П.А. Куйбин, В.Л. Окулов, 2003; С.И.Шторк и др., 2006; В.Т.Волов, 2006, 2011). Расчеты также показали, частота прецессии растет с ростом мощности источника тепла (рисунок 8). Результат получен впервые и качественно коррелирует с зависимостью частоты изгибной моды от скачка плотности газа, полученной при аналогичных условиях в главе 2 и с зависимостям частоты прецессии от массового расхода газовой фазы для газожидкостных потоков (А.А.Дектерев, А.А.Гаврилов, А.А.Дектерев, 2011).

Рисунок 8. Зависимость частоты прецессии от мощности источника N. Массовый расход через тангенциальный и аксиальный подводы Qt = 4 г/с и Qax =1 г/с соответственно.

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

В заключении перечислены основные результаты, полученные при выполнении данной диссертационной работы:

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

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

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

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

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

Публикации по теме диссертации Статьи в изданиях, рекомендуемых ВАК России 1. Завершинский, И.П. Эволюция вихря Рэнкина в газе с источником тепловыделения [Текст] / И.П.Завершинский, А.И.Климов, Н.Е.Молевич, Д.П.Порфирьев // Письма в ЖТФ. - 2009. - Т 35. - № 7. - С. 106-111.

2. Завершинский, И.П. Об особенности свечения поперечного разряда в зависимости от массового расхода газа в вихревой камере [Текст] / И.П.Завершинский, А.И.Климов, В.Г.Макарян, Н.Е.Молевич, И.А. Моралев, Д.П.Порфирьев // Письма в ЖТФ. - 2009. - Т.35. - № 24. - С. 59-67.

3. Завершинский, И.П. Численное моделирование вихревого потока в присутствии поперечного разряда постоянного тока [Текст] / И.П.Завершинский, А.И.Климов, В.Г.Макарян, Н.Е.Молевич, И.А. Моралев, Д.П.Порфирьев//Теплофизика высоких температур. - 2010. - № (дополнительный). - С. 157-161.

4. Завершинский И.П. Исследование газодинамики закрученного потока в присутствии объемного разряда постоянного тока [Текст] / Завершинский И.П., Макарян В.Г., Молевич Н.Е., Порфирьев Д.П. // Физическое образование в вузах. - 2011. - Т.17. - №1. - С.П19.

5. Завершинский, И.П. Структура высокочастотного емкостного разряда в закрученном потоке воздуха при атмосферном давлении [Текст] / И.П.Завершинский, А.И.Климов, В.Г.Макарян, Н.Е.Молевич, И.А. Моралев, Д.П.Порфирьев // Письма в ЖТФ. - 2011. - Т.37. - № 23. - С. 68-75.

6. Troshkin, N.V. Acoustical instability of gas flow with heat release [Текст] / Troshkin N.V., Molevich N.E., Porfiriev D.P. // Acta Acustica uni Acustica. - 2011 - V. 97. -Suppl.1. - P. S76.

В других изданиях 7. Завершинский, И.П. Vortex structures in excited molecular gas-plasma media [Текст] / И.П.Завершинский, А.И.Климов, Н.Е.Молевич, Порфирьев Д.П.

// 8th Workshop on Magnetoplasma Aerodynamics. Moscow, 2009.

8. Завершинский, И.П. Устойчивость вихря Рэнкина в газе с источником тепловыделения [Текст] / И.П.Завершинский, Н.Е.Молевич, Д.П.Порфирьев // Материалы международной междисциплинарной научной конференции VI Курдюмовские чтения. Идеи синергетики в естественных науках, Тверь: ТГУ, 2010. С.91-92.

9. Завершинский И.П. Численное моделирование структуры вихревого потока в присутствии разряда. Перспективные информационные технологии для авиации и космоса (ПИТ-2010) [Текст] / Завершинский И.П., Макарян В.Г., Молевич Н.Е., Порфирьев Д.П. // Труды Международной конференции с элементами научной школы для молодежи. Самара, 2010. С. 569-573.

10. Завершинский И.П. Факельный разряд в закрученном потоке воздуха [Текст] / Завершинский И.П., Макарян В.Г., Молевич Н.Е., Порфирьев Д.П. // Труды Международной конференции Физика высокочастотных разрядов, 5-апреля 2011 г., Казань. С.103-105.

11. Завершинский И.П. Structures of the single electrode RF gas discharge in swirling air flow [Текст] / Завершинский И.П., Макарян В.Г., Молевич Н.Е., Порфирьев Д.П. // 10th Int. Workshop on Magneto-Plasma Aerodynamics, March 22 - 24, 2011, Moscow. С.62-12. Makaryan V.G., Molevich N.E., Porfiriev D.P., Sugak S.S., Zavershinsky I.P. Nonstationary structure of swirling flows in open tube channel with heat sources // 11 Int.Workshop on Magneto-Plasma Aerodynamics,April 9 Ц13, 2012, M.,P. 6263.

13. Завершинский И.П., Молевич Н.Е., Порфирьев Д.П., Сугак С.С.

Численное моделирование структуры вихревых потоков в присутствии газовых разрядов // Материалы III международной конференции Математическая физика и ее приложения. Самара. 2012. C. 136-137.

14. Molevich N.E., Porfiriev D.P., Sugak S.S., Zavershinsky I.P. Precessing vortex core in a wake with counterflow in presence of heat source [Текст] // Proc.

Euromech Fluid Mechanics Conference (EFMC9), Roma, Italy, 2012. P. 201-202.

Подписано в печать 20 ноября 2012 г. Формат 60х48/16. Тираж 100 экз.

Отпечатано с готового оригинал-макета заказчика.

443086, г. Самара, Московское шоссе, 34, СГАУ.

Авторефераты по всем темам  >>  Авторефераты по разным специальностям