Физико-математические модели и численные исследования трансформации массового спектра облачных капель за высокоскоростным летательным аппаратом тема диссертации и автореферата по ВАК РФ 01.02.05, кандидат физико-математических наук Здор, Александр Геннадьевич
- Специальность ВАК РФ01.02.05
- Количество страниц 174
Оглавление диссертации кандидат физико-математических наук Здор, Александр Геннадьевич
Список основных обозначений.
Введение.
1. СИСТЕМА УРАВНЕНИЙ ЭВОЛЮЦИИ
ПОЛИДИСПЕРСНОГО КОНДЕНСАТА.
1.1. Физические и прикладные аспекты моделирования многофазных течений.
1.2. Основные допущения и оценки их адекватности при построении моделей динамики диспергированного вещества.
1.3. Физические процессы в одиночной частице.
1.3.1. Фазовые переходы.
1.3.2. Аэродинамическое воздействие несущего потока на частицы.
1.4. Коагуляция и дробление частиц.
1.4.1 Методы описания эволюции полидисперсных включений.
1.4.2. Описание исходов парного взаимодействия капель.
1.4.3. Источниковые члены в уравнении эволюции функции распределения частиц по радиусам.
1.4.4. Уравнение эволюции массовой концентрации фракции полидисперсного конденсата.
1.4.5. Вывод уравнений эволюции импульса фракций.
1.4.6. Вывод уравнений эволюции полной удельной энергии фракций.
1.4.7 Система уравнений эволюции полидисперсного конденсата в случае конечного количества фракций.
1.5. Выводы по главе 1.
2. ВЗАИМОДЕЙСТВИЕ ДИСПЕРГИРОВАННОГО
ВЕЩЕСТВА С ТУРБУЛЕНТНОЙ НЕСУЩЕЙ СРЕДОЙ.
Введение.
2.1. Основные задачи и методы описания динамики частиц в турбулентных потоках.
2.2. Броуновское движение.
2.2.1. Эволюция одиночной броуновской частицы.
2.2.2. Частота взаимных столкновений в системе броуновских частиц.
2.3. Инерционные частицы в поле турбулентных пульсаций скорости несущей среды.
2.3.1. Эволюция одиночной частицы.
2.3.2. Оценка частоты взаимных столкновений капель, вовлечённых в пульсационное движение.
2.3.3. Диффузионные слагаемые в системе уравнений эволюции полидисперсного конденсата.
2.4. Выводы по главе 2.
3. РЕЗУЛЬТАТЫ ЧИСЛЕННОГО АНАЛИЗА.
Введение.
3.1. Течения газа.
3.1.1. Тестовые одномерные задачи.
3.1.2. Обтекание плоской пластины потоком сжимаемого вязкого теплопроводного газа под углом атаки.
3.2. Результаты расчета параметров полидисперсного конденсата вблизи поверхности и в следе за плоской пластиной.
3.2.1. Разные фракции в рамках одного режима «на бесконечности».
3.2.2. Влияние числа Маха набегающего потока.
3.2.3. Влияние угла атаки набегающего потока.
3.2.4. Влияние учета процессов коагуляции/дробления капель.:.
3.3 Выводы по главе 3.
Рекомендованный список диссертаций по специальности «Механика жидкости, газа и плазмы», 01.02.05 шифр ВАК
Обтекание тел потоком газовзвеси2005 год, доктор физико-математических наук Циркунов, Юрий Михайлович
Исследование турбулентной двухфазной неравновесной многокомпонентной полидисперсной неизотермической осесимметричной струи с фазовыми переходами, коагуляцией и дроблением капель1998 год, кандидат технических наук Парамонов, Андрей Николаевич
Математическое моделирование некоторых задач пограничного слоя в газовзвесях2005 год, кандидат физико-математических наук Забарин, Владимир Иванович
Многофазные турбулентные струйные течения в элементах тепловых двигателей, установках и устройствах различного назначения: Математическое моделирование, численное и экспериментальное исследования1999 год, доктор технических наук Зуев, Юрий Владимирович
Моделирование динамики парокапельных сред в процессе регазификации2019 год, кандидат наук Тукмакова Надежда Алексеевна
Введение диссертации (часть автореферата) на тему «Физико-математические модели и численные исследования трансформации массового спектра облачных капель за высокоскоростным летательным аппаратом»
Одним из важнейших направлений механики гетерогенных сред является исследование эволюции двухфазных потоков. Наличие дисперсных включений в течениях газа или жидкости не только влияет на количественные характеристики последних, но и существенно обогащает физику происходящих процессов. Исследования в этой области естествознания представляют научный интерес, а также имеют массу практических приложений. Среди них стоит упомянуть:
Метеорологию. Фактически все задачи метеорологии касаются физики и механики многофазных течений, где в качестве дисперсных включений выступают капли воды, льдистые частицы, пылинки. Изучение процессов зарождения и эволюции облаков актуальны в связи с проблемами локального регулирования климата (управление выпадением осадков), детектирования распространения атмосферных фронтов и предсказания погоды.
Экологию. Практический интерес представляет мониторинг состояния атмосферы с учётом влияния дисперсных включений природного (извержения вулканов, бури) и антропогенного происхождения. Наличие аэрозоля может существенно изменять скорости различных химических реакций, сопровождающиеся появлением и трансформациями активных соединений (озон, окислы азота, серы), влиять на процессы переноса солнечного и теплового излучения, что способно привести к нарушениям естественного температурного баланса.
Авиационную и ракетную технику. Здесь приложения особенно многочисленны. Исторически первыми появились задачи (актуальные и сейчас), связанные с моделированием динамики обледенения летательных аппаратов при движении в облаках и разработкой комплексов мер по его предотвращению или уменьшению. Сугубо практический интерес представляет изучение эрозионного и абразивного воздействия аэрозоля на обшивку ЛА. Наличие дисперсной фазы в набегающем потоке вызывает увеличение силовых и тепловых нагрузок на обтекаемые поверхности, особенно в затупленных головных частях фрагментов компоновки, что приводит к ускоренному разрушению покрытия. Также практическую значимость представляют оценки изменения аэродинамических характеристик планера самолёта при движении в атмосфере с высокой влажностью и большой удельной долей капель конденсата (взлёт и посадка в дождь). Отдельного внимания заслуживает проблема наблюдения за траекториями движения летательных аппаратов. Капли, входящие в состав атмосферного аэрозоля, при взаимодействии с создаваемыми ЛА газодинамическими неоднородностями участвуют в целом ряде физических процессов. Они могут испаряться в горячих приповерхностных зонах, дробиться при прохождении разрывов, коагулировать и снова конденсироваться в достаточно "холодных" областях течения. Все эти явления способны приводить к температурной и скоростной неравновесности фаз в потоке, а также существенно искажать массовый спектр аэрозоля. Возникающие при этом неоднородности могут быть хорошо заметными в некоторых диапазонах длин волн солнечного или сканирующего излучения. Подобные следы испарения визуально наблюдались при движении самолёта в низкой облачной дымке над поверхностью моря и имели характерную длину порядка нескольких десятков километров.
Машиностроение. Изучение физики двухфазных потоков, их эволюции и закономерностей взаимодействия с обтекаемыми поверхностями при наличии диспергированного вещества актуально в задачах проектировании оптимальных сопел разного назначения (реактивные двигатели, эжекторы), трактов двигателей, турбинной техники. Исследования в этой области играют важную роль при разработке методов и технологий нанесения покрытий путём напыления, упрочнения поверхностей деталей машин при ударном воздействии, а также управления движением примесей, что имеет -отношение к фильтрации и сепарации взвесей.
История исследования однофазных течений (газовая динамика и гидромеханика) насчитывает уже не одну сотню лет, и, несмотря на хорошую развитость моделей для большинства важных на практике случаев (здесь стоит упомянуть системы уравнений Эйлера, Навье-Стокса, Рейнольдса) по-прежнему изобилует массой нерешённых вопросов. Многофазные потоки по своим физическим свойствам гораздо сложнее, характеризуются большим числом степеней свободы, а, следовательно, и большим количеством необходимых характерных величин. Даже сравнительно простые синтетические модели кинетики полидисперсного конденсата включают порядка десятка различных безразмерных параметров подобия. Однако каждая из них необходимо содержит определённые способы описания двух важных аспектов:
1) Физических процессов, в которых участвует одиночная частица дисперсной фазы при взагшодействии с несущей средой. Во многих случаях дисперсная фаза в многофазном течении обменивается массой и теплом с несущей средой. В течение долгого времени эти процессы являются предметом теоретических и экспериментальных исследований различных научных школ. Как правило, эти явления определяются широкой совокупностью различных определяющих параметров: температурами капель и окружающей среды, её компонентным составом (предельными здесь выступают случаи испарения в собственный пар или в пассивный газ), скоростным скольжением частиц, числом Кнудсена их обтекания и другими величинами. Соответствующие зависимости для континуального и свободномолекулярного пределов достаточно хорошо изучены. Для промежуточных режимов существуют интерполирующие формулы.
Виды силового воздействия на частицу со стороны несущего потока также многочисленны. Как правило, принимается гипотеза об их независимости, что позволяет рассматривать разные эффекты по-отдельности. К наиболее существенным из них относится сила аэродинамического сопротивления. В литературе опубликованы аналитические зависимости коэффициента сопротивления от определяющих параметров для сферических частиц, прошедшие серьёзную экспериментальную проверку. Также существуют поправки к ним, учитывающие деформацию капель. При исследованиях эволюции дисперсных включений важную роль могут играть сила Сэфмана (существенная в областях сдвиговых течений) и сила Магнуса (при учёте вращения частиц). Необходимо упомянуть и стохастические воздействия несущей среды на частицы. Действительно, рассматриваемые потоки являются, как правило, турбулентными вследствие чего их параметры (скорость, температура, плотность) испытывают пульсации. Находясь в случайном поле флуктуаций скорости, частица испытывает диффузионные смещения. В случае малых радиусов частиц становится важным учет их броуновского движения. Исследованию этих явлений посвящено большое количество работ. При анализе неоднородных течений необходимо принимать во внимание форетические силы разной природы (турбофорез, термофорез). Также важно оценивать степень деформации и устойчивость капли, поскольку в ней могут развиваться колебания с незатухающими модами, приводящие к её разрушению. Рассматриваемые частицы в некоторых случаях могут нести электрический заряд, что приводит к необходимости учитывать силы их взаимодействия друг с другом и с внешними электромагнитными полями.
2) Эволюции ансамбля взаимодействующих частиц. Для моделирования эволюции системы дисперсных включений в несущем потоке в опубликованных работах предлагается несколько разных подходов. Все они подразумевают использование конечного количества параметров, характеризующих состояние частицы. В состав последних обычно включается её радиус, скорость центра масс, среднеобъёмная температура, кинетический момент вращения относительно центра масс. Для справедливости предположения о конечности количества параметров необходимо выполнение требований малости градиентов температур и скоростей вещества частицы, а также затухания вызванных взаимными столкновениями возмущений за промежуток времени между соударениями. Наиболее последовательные подходы подразумевают либо использование кинетического уравнения для функции распределения частиц системы по скоростям, пространственным координатам, радиусам, температурам, моментам импульса и другим представляющим интерес переменным, либо рассмотрение задачи о совместной эволюции параметров состояния капель вдоль их траекторий для достаточно больших совокупностей (до десятков миллионов «участников»). Их реализация на практике требует больших вычислительных затрат (например, использования методов прямого статистического моделирования). Несколько более экономичными и удобными в приложениях являются гидродинамические модели, в которых двухфазный поток рассматривается как многокомпонентная сплошная среда (газ и фракции диспергированного вещества). Допускается, что каждая из её составляющих может обладать собственными локальными значениями параметров состояния: температурой, скоростью, числовой и массовой концентрацией. Идентификация фракций конденсата осуществляется при помощи только одной величины — радиуса частицы, который также зависит от пространственных и временных координат. При учёте процессов коагуляции капель в результате взаимных столкновений необходимо адекватно описывать распределение массы, импульса, энергии и других параметров конденсата по фракциям. Широкое распространение для моделирования этих явлений получили непрерывный и кинетический подходы.
Основным предметом данной диссертации является эволюция и свойства двухфазных полидисперсных течений с малой объёмной долей конденсата около поверхности твёрдого тела и в его следе.
Первые модели, аналитические результаты и данные численных исследований движения частиц в потенциальных течениях несущей жидкости около сферы и цилиндра были обобщены в монографии Фукса Н.А. [89]. Вопросы мелкомасштабного движения вязкой несущей среды около частиц примеси рассматривались в монографиях Левина JT.M. [58] и ВолощукаВ.М. [16]. Родственные задачи в случае малых чисел Рейнольдса обтекания частиц примеси, а также их гидродинамическое взаимодействие подробно освещены в монографии ВолощукаВ.М. и Седунова Ю.С. [17]. Теоретические и экспериментальные данные о физических процессах в каплях конденсата (тепломассообмен с окружающей средой, силовое воздействие) при их движении в неоднородных газовых потоках систематизированы в работах Гилинского М. М. [22], Нигматулина Р.И. [66], Стасенко А.Л. [22, 76-81, 85], СтернинаЛ.Е. [86-88]. Наиболее известные исследовательские группы и отдельные ученые, занимавшиеся вопросами двухфазной газодинамики, проводили работы в Институте механики МГУ им. М.В. Ломоносова (Нигматулин Р.И., Стулов В.П., Осипцов А.Н., Шапиро Е.Г., Гилинский М.М. и др. [22, 66, 68, 69]), ЦИАМ им. П.И. Баранова (Крайко А.Н. [52]), ЦАГИ им. проф. Н.Е. Жуковского (Стасенко А.Л., Василевский Э.Б., АсмоловЕ.С. [1, 2, 12, 40, 76-85]), Ленинградском (Санкт-Петербургском) государственном университете (Матвеев С.К., Лашков В.А. и др. [55, 56, 61, 62]), Институте высоких температур АН СССР (РАН) (Полежаев Ю.В., Михатулин Д.С. и др. [10, 70, 71]), Томском государственном университете (Гришин A.M., Забарин В.И. и др. [25]), ЦНИИМаш (ВасинА.В. и др. [10]), БГТУ «Военмех» им. Д.Ф.Устинова (Циркунов Ю.М., Тарасова Н.В., Волков А.Н. и др. [14, 91-93]), ВЦ АН СССР (РАН) (Давыдов Ю.М., Толстых А.И. [26, 82, 83]).
Отдельным научным направлением при исследовании эволюции многофазных потоков является учёт влияния турбулентных пульсаций параметров несущей среды. В результате стохастических воздействий инерционные частицы дисперсной фазы также оказывается в той или иной мере вовлеченными в пульсационное движение,« что вызывает их диффузионные перемещения и влияет на частоту их взаимных столкновений. Изучением этих вопросов занимается целый ряд исследователей, существуют различные подходы к моделированию процессов. Наибольшее влияние на формирование взглядов автора диссертации в этой области оказали работы Деревича И.В. [29, 30, 99, 100] и Зайчика Л.И. [32-35].
В результате анализа опубликованных работ можно заключить, что теория многофазных потоков в настоящее время интенсивно развивается. Несмотря на многочисленные теоретические и экспериментальные работы, проведённые разными научными группами, можно отметить, что единый подход для математического моделирования рассматриваемых течений отсутствует, большинство исследований узко специализированы. Постановка опытов, данные которых позволяют тестировать, отлаживать и совершенствовать модели, крайне затруднительна как по техническим, так и методологическим причинам. В связи с этим, важная роль отводится численному эксперименту в рамках синтетических моделей эволюции полидисперсных потоков, принимающих в учёт сразу несколько физических процессов. Казалось бы, такой подход ведёт к резкому усложнению: при описании коагуляции/дробления, фазовых переходов и аэродинамического воздействия несущего газа на частицы основные зависимости содержат более десятка безразмерных критериев подобия. Основным в пользу предлагаемого подхода соображением служит то, что зачастую газодинамическая картина обтекания рассматриваемых тел крайне неоднородна, включая в себя области нагретого и «холодного» газа, а также, возможно, сильные разрывы. Следовательно, выделить априори доминирующий эффект (скажем, только фазовый переход или только коагуляцию) крайне затруднительно. На этих основаниях представляется целесообразным использовать модель, обогащённую учтёнными физическими процессами, а в численной реализации системы уравнений предусмотреть возможность расчёта каждого из них. При этом «главное» слагаемое должно выделяться автоматически в зависимости от локальных условий в точке течения.
Задачи, касающиеся эволюции полидисперсных течений и трансформации массового спектра конденсата вблизи поверхностей летательных аппаратов и в их следах, имеют большое прикладное значение. Их исследованию посвящено лишь небольшое количество публикаций, что подтверждает актуальность представляемой работы.
Цель диссертационной работы состояла в построении синтетической физико-математической модели эволюции полидисперсного двухфазного течения, создании численного кода и теоретических и методических численных исследованиях изменения дисперсности и динамики капель водного конденсата при обтекании модельной конфигурации ЛА высокоскоростным потоком.
Главными задачами работы являлись:
1. Анализ современного состояния исследований многофазных полидисперсных потоков.
2. Выбор математических моделей для описания физических явлений, происходящих с частицей конденсата: фазовых переходов, теплообмена, обмена импульсом с несущим газом.
3. Выбор модели эволюции ансамбля частиц дисперсной фазы, взаимодействующих с несущим газом и друг с другом.
4. Выбор, сравнительный анализ и дальнейшая разработка моделей взаимодействия дисперсной фазы с турбулентной несущей средой.
5. Построение синтетических моделей эволюции полидисперсного конденсата.
6. Создание численных кодов для расчёта течений газа и дисперсной фазы и отладка программ.
7. Методические численные исследования обтекания модельной конфигурации ЛА и его ближнего следа для разных режимов «на бесконечности» и разных комбинаций учитываемых эффектов.
8. Анализ численных результатов и определение диапазонов влияния параметров. Оценка важности учёта отдельных физических эффектов.
9. Оценки адекватности использованных физических моделей и результатов численного анализа.
Первая глава носит обзорный характер. В ней отмечен вклад различных научных школ в изучение физики явлений и построение математических моделей полидисперсных многофазных потоков. Определяется объект исследования: сверхзвуковое обтекание ЛА полидисперсным двухфазным потоком с совместным учётом ряда физических эффектов: силового воздействия течения на капли, фазовых переходов, коагуляции, дробления частиц под действием аэродинамических сил и при взаимных столкновениях. На основании уравнения Лиувилля приводится подробный вывод системы уравнений эволюции полидисперсного конденсата для гидродинамической модели в рамках непрерывного подхода.
Во второй главе обсуждаются результаты моделирования поведения инерционных дисперсных включений в поле турбулентных и молекулярных " флуктуаций скорости несущего потока. Относительное движение частиц конденсата и несущего газа можно условно разделить на две взаимосвязанные части: осреднённое скольжение фаз и случайные смещения. Как упоминалось выше, взаимодействие капель с пульсациями является одной из причин их взаимных соударений, а также вызывает диффузионные перемещения частиц. Информация о времени между столкновениями капель важна при анализе процессов их коагуляции. В представленной работе основной интерес концентрировался на анализе зависимости от свойств турбулентности флуктуационных параметров частиц: функции вовлечения в пульсационное движение, коэффициента диффузии и частоты взаимных столкновений дисперсных включений.
Третья глава посвящена численному анализу эволюции двухфазного полидисперсного потока при обтекании модельной конфигурации ЛА на основе предложенной синтетической модели. В качестве инструмента для проведения расчётов был использован метод крупных частиц Белоцерковского-Давыдова.
Неоспоримым достоинством метода является возможность проводить сквозной счёт, без предварительного выделения особенностей течения. Поскольку содержание дисперсной фазы в набегающем потоке принималось малым (отношение массовых концентраций дисперсного вещества и несущего газа порядка нескольких процентов), обратное влияние конденсата на газ не учитывалось. Поэтому вначале проводились исключительно газодинамические расчёты. Дальнейшие вычисления касались исследования эволюции трёхфракционного конденсата в поле течения около пластины. С целью установления влияния различных входных данных задачи выполнена серия параметрических расчётов. Проведено сравнение результатов численного анализа эволюции разных фракций в рамках одного режима «на бесконечности», изучено влияние параметров набегающего потока (числа Маха, угла атаки), показана важность учёта коагуляции/дробления капель при взаимных столкновениях при определении параметров конденсата вблизи обтекаемого тела и в следе. Автор защищает следующие положения:
• Синтетическую модель эволюции полидисперсных двухфазных течений, учитывающую ряд физических процессов: фазовые переходы, коагуляцию, дробление капель.
• Результаты верификации созданных численных кодов, позволяющих проводить исследования полидисперсных многофазных потоков в широком диапазоне параметров вблизи поверхностей тел и в следах.
• ' Результаты серии параметрических расчётов обтекания модельной конфигурации ЛА сверхзвуковым полидисперсным двухфазным потоком для случая частиц водного конденсата. Исследования закономерности влияния различных параметров (числа Маха набегающего потока, угла атаки, числовой концентрации дисперсной фазы в набегающем потоке) на картину течения.
• Результаты численных исследований влияния процессов коагуляции на трансформацию массового спектра капель в ударно-волновой картине обтекания ЛА.
• Предложенную (в соавторстве с Деревичем И.В.) замкнутую модель, позволяющую рассчитывать стохастические параметры инерционных частиц в однородной изотропной турбулентности в рамках подхода Лагранжа и сравнение получаемых по ней характеристик с результатами других моделей.
Результаты выполненных исследований по теме диссертации докладывались и получили положительную оценку на следующих конференциях и семинарах:
1. Научные конференции МФТИ: ХЬУ1 (2003 г.), ХЬУИ (2004 г.), ХЬУШ (2005 г.), ХЫХ (2006 г.), Ь (2007 г.), Ы (2008 г.).
2. Юбилейный международный семинар по струйным, отрывным и нестационарным течениям, Санкт-Петербург, 2004 г.
3. Международная научно-техническая конференция «Фундаментальные проблемы высокоскоростных течений», Жуковский, 2004 г.
4. Академические чтения по космонавтике («Королёвские чтения»), Москва: XXIX (2005 г.), XXX (2006 г.), XXXI (2007 г.), XXXII (2008 г.), XXXIII (2009 г.).
5. Научно-техническая школа-семинар «Современные проблемы аэрокосмической науки и техники» в рамках 2-ой Международной выставки авиационно-технического творчества молодёжи «ЮниМАКС-2005» и «Вузовская наука», Жуковский, 2005 г.
6. IX Международный симпозиум «Актуальные проблемы машиностроения и механики сплошных и сыпучих сред», Москва, 2006 г.
7. VI Международная конференция по неравновесным процессам в соплах и струях (№>N.1-2006), Санкт-Петербург, 2006 г.
8. XVI Школа-семинар молодых учёных и специалистов под руководством академика РАН А.И. Леонтьева, Санкт-Петербург, 2007 г.
9. Международный авиационно-космический научно-гуманитарный семинар имени С.М. Белоцерковского, Москва, 2007 г.
10. Научно-технический семинар НИО-8 ЦАГИ, Жуковский, 2007 г.
11. VI Минский Международный Форум по Тепломассообмену, Минск, 2008 г.
12. 7-ой семинар ЦАГИ-ОЫЕКА, Жуковский, 23-26 сентября 2008.
13. Семинар по аэродинамике ЦАГИ-ИТПМ, Жуковский-Новосибирск, 18 ноября 2008.
По результатам диссертационного исследования автором опубликовано 20 работ (5 статьей и 15 докладов), из них 11 без соавторов.
Автор глубоко признателен своему научному руководителю — доктору технических наук, профессору Стасенко А.Л., который оказал большое влияние на формирование научных интересов и взглядов соискателя, а также привлёк его внимание к рассмотренной задаче и принимал активное участие в обсуждении результатов исследований. Автор благодарен доктору технических наук, профессору Вышинскому В.В. за исчерпывающую информацию по ряду научных вопросов, представлявших интерес для автора, и за справедливые критические замечания по некоторым пунктам, улучшившие качество работы. Также автор хотел бы выразить благодарность коллеге по работе в ЦАГИ, кандидату технических наук, доценту Моллесон Г.В. за помощь в создании численных кодов.
Похожие диссертационные работы по специальности «Механика жидкости, газа и плазмы», 01.02.05 шифр ВАК
Нестационарное взаимодействие сверхзвуковых потоков газовзвеси с телами и преградами2008 год, кандидат физико-математических наук Семенов, Владимир Владимирович
Рассеяние частиц примеси при обтекании тел высокоскоростным потоком газовзвеси2008 год, кандидат физико-математических наук Панфилов, Сергей Владимирович
Тепломассообмен в двухфазных многокомпонентных турбулентных струйных течениях2006 год, кандидат технических наук Шустрова, Виктория Юрьевна
Динамика микро и макрообъёмов сброшенной с самолёта жидкости2010 год, кандидат технических наук Кудров, Максим Александрович
Численное моделирование особенностей течений идеального газа и двухфазных смесей газа с частицами2011 год, кандидат физико-математических наук Пьянков, Кирилл Сергеевич
Заключение диссертации по теме «Механика жидкости, газа и плазмы», Здор, Александр Геннадьевич
Основные результаты диссертационной работы состоят в следующем:
1. На основании анализа опубликованных разными авторами результатов, касающихся исследования отдельных процессов с каплями конденсата (фазовых переходов, коагуляции, дробления, обмена импульсом и энергией с окружающим газом), предложена синтетическая модель эволюции полидисперсных двухфазных течений, совместно учитывающая ряд физических явлений.
2. Созданы численные алгоритмы, позволяющие проводить расчёты обтекания тел полидисперсными многофазными потоками в широком диапазоне параметров вблизи их поверхностей и в следе.
3. Проведена серия параметрических расчётов обтекания модельной конфигурации ЛА сверхзвуковым полидисперсным двухфазным потоком для случая частиц водного конденсата. Исследованы закономерности влияния различных параметров на картину течения вблизи обтекаемого тела и в его следе.
4. На основании численных данных показана необходимость учёта процессов коагуляции капель в ударно-волновой картине обтекания ЛА.
5. Предложена (в соавторстве с Деревичем И.В.) замкнутая модель, позволяющая рассчитывать стохастические параметры инерционных частиц в однородной изотропной турбулентности в представлении Лагранжа. Проведено сравнение получаемых характеристик с данными, получаемыми по другим моделям.
В целом, в данной диссертации сформулирована задача об эволюции полидисперсного конденсата и исследовании трансформации его массового спектра для течения вблизи поверхности твердого тела и в его следе. Предложена синтетическая физико-математическая модель, учитывающая достаточно широкую совокупность явлений: обмен массой, теплом и импульсом фракций с несущим газом, дробление капель под действием аэродинамических сил, коагуляцию/дробление при их взаимных столкновениях, вовлечение частиц в турбулентность. Теоретический анализ различных аспектов проблемы, разработка математических и численных моделей, а также выполненные методические расчёты позволили получить представления об особенностях поведения конденсированной фазы при обтекании поверхности ЛА при различных сочетаниях размеров частиц и условий в набегающем потоке, установить роль некоторых определяющих параметров задачи, получить оценки порядков основных величин, характеризующих диспергированное вещество. Направления дальнейших исследований в рамках предложенной модели в этой области двухфазной аэродинамики, по мнению автора диссертации, могли бы состоять в следующем:
1. Расчёты зон смачивания и определение потоков массы конденсата к поверхности обтекаемого ЛА.
2. Расчёты поправок к аэродинамическим характеристикам профилей при высокой массовой доле конденсата в набегающем потоке (роль частиц в генерации и подавлении пульсаций скорости несущей среды).
3. Учёт обратного влияния дисперсной фазы на характеристики газа, в том числе и на параметры турбулентности.
4. Учёт возникающей в некоторых классах течений (пограничные слои, струи) анизотропии турбулентности, а также оценки влияния неоднородностей потока (ударные волны) на параметры стохастического движения капель.
Все перечисленные выше задачи являются актуальными как с научной, так и прикладной точки зрения.
ЗАКЛЮЧЕНИЕ
Список литературы диссертационного исследования кандидат физико-математических наук Здор, Александр Геннадьевич, 2008 год
1. Асмолов Е.С. О динамике сферической частицы в ламинарном пограничном слое // Изв. АН СССР, МЖГ. 1990, № 6, с. 91-96.
2. Асмолов Е.С. О движении дисперсной примеси в ламинарном пограничном слое на плоской пластине // Изв. РАН, МЖГ. 1992, № 1, с. 66-73.
3. БабухаЛ.Г., ШрайберА.А. Взаимодействие частиц полидисперсного вещества в двухфазных потоках. — Киев: Наук, думка, 1972.
4. Бахвалов Н., Жидков Н., Кобельков Г. Численные методы. — М. — С.-Пб.: Физматлит, Невский Диалект, Лаборатория базовых знаний, 2001.
5. Белоцерковский О.М., Давыдов Ю.М. Метод крупных частиц в газовой динамике. -М.: Наука, 1982.
6. Бусройд Р. Течение газа со взвешенными частицами. Пер. с англ. М.: Мир, 1975.
7. Вызова Н.Л., Иванов В.Н., ГаргерЕ.К. Турбулентность в пограничном слое атмосферы. Л.: Гидрометеоиздат, 1989.
8. Бэтчелор Дж. Теория однородной турбулентности. Изд. Иностранной литературы, М. 1955.
9. ВасенинИ.М., Архипов В.А., Бутов В.Г., Глазунов A.A., Трофимов В.Ф. Газовая динамика двухфазных течений в соплах. Издательство Томского университета, 1986.
10. Васин A.B., Михатулин Д.С., Полежаев Ю.В. О влиянии теплового состояния материалов на их эрозионную стойкость в запыленном газовом потоке // Изв. АН СССР, МЖГ. 1985, № 6, с.172-175.
11. Ван Кампен Н.Г. Стохастические процессы в физике и химии. М.: Высшая школа, 1990.
12. Василевский Э.Б., Осипцов А.Н., Чирихин A.B., Яковлева Л.В. Теплообмен на лобовой поверхности затупленного тела в высокоскоростном потоке, содержащем малоинерционные частицы // ИФЖ. 2001, т. 74, № 6, с. 29-37.
13. Васильков А.П. Окрестность критической точки затупленного тела в гиперзвуковом двухфазном потоке // Изв. АН СССР, МЖГ. 1975, № 5, с. 121-129.
14. Волков А.Н., Семёнов B.B., ЦиркуновЮ.М. Влияние моно- и полидисперсной примеси на течение и теплообмен при сверхзвуковом обтекании затупленного тела потоком газовзвеси // Математическое моделирование. 2004, т. 16, №7, с. 6-12.
15. Волков В.А. Исследование параметров двухфазной среды при сверхзвуковом обтекании затупленных тел // Труды XIX научной конф. МФТИ, 1973. Сер. Аэромех. и проц. управл. Долгопрудный: МФТИ, 1974, с. 14-21.
16. Волощук В.М. Введение в гидродинамику грубодисперсных аэрозолей. Л.: Гидрометеоиздат, 1971.
17. Волщук В.М., Седунов Ю.С. Процессы коагуляции в дисперсных системах. -Л.: Гидрометеоиздат, 1975.
18. ВороничИ.В. Сравнительный анализ группы численных методов газовой динамики: Учебное пособие. М.: МФТИ, 2007.
19. Гавин Л.Б., Шрайбер A.A. Турбулентные течения газа с частицами // Итоги науки и техники, серия МЖГ. т. 25. М.: ВИНИТИ, 1991.-е. 90-182.
20. Галкин B.C., Коган М.Н., Фридлендер О.Г. О некоторых кинетических эффектах в течениях сплошной среды. Изв. АН СССР, МЖГ, 1976, № 3, с. 13-21.
21. Гонор А.Л., РивкиндВ.Я. Динамика капли // Итоги науки и техники. Сер. Механика жидкости и газа. М.: Машиностроение, 1984.
22. Гилинский М.М., Стасенко А.Л. Сверхзвуковые газодисперсные струи. М.: Машиностроение 1990.
23. Горбачев Ю.Е., Лунькин Ю.П. Граничные условия в задаче о течении гетерогенной смеси // Письма в ЖТФ, 1980, т.6, Вып. 5, с. 299-301.
24. Горбачев Ю.Е., Круглов В.Ю. Расчет параметров течения двухфазной смеси при обтекании сферы с учетом столкновений частиц примеси между собой // Изв. АН СССР, МЖГ, 1989, № 4, с. 93-96.
25. Гришин A.M., Забарин В.И. Трение и теплообмен в двухфазном пограничном слое на пластине // ПМТФ, 1988, 1988, № 4, с. 78-86.
26. Давыдов Ю.М., Нигматулин Р.И. Расчет внешнего обтекания затупленных тел гетерогенным потоком газа с каплями или частицами // Докл. АН СССР, 1981, т. 259, № 1, с. 57-60.
27. Дрейк Р., Беккер Дж. Теплопередача от шара к разреженному газу в сверхзвуковом потоке // Вопросы ракетной техники. М.: Изд. иностр. ли-ры, 1953, вып. 2(14), с. 54-69.
28. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. -М.: Физматлит, 1963.
29. ДеревичИ.В. Столкновения частиц в турбулентном потоке. Известия РАН, МЖГ, 1996, №2, с. 104.
30. ДеревичИ.В., Здор А.Г. Замкнутая модель флуктуационного движения частиц в турбулентном потоке. Принято к опубликованию в Известия РАН, МЖГ, 2008.
31. Зайчик Л.И. Оценка времени между столкновениями дисперсных частиц в турбулентном потоке // ТВТ, 1998, т. 36, №3, с. 456 460.
32. Зайчик Л.И., Алипченков В.М., Петров О.Ф. Кластеризация заряженных частиц в условиях изотропной турбулентности // ТВТ, 2004, т. 42, № 6, с. 908-916.
33. Зайчик Л.И., Алипченков В.М. Кластеризация малоинерционных частиц в изотропной турбулентности // ТВТ, 2007, т. 45, № 1, с. 66-76.
34. Зайчик Л.И., Алипченков В.М. Столкновения частиц в турбулентном потоке // Изв. РАН, МЖГ, 2007, № 3, с. 94-109.
35. Здор А.Г. Эволюция массового спектра капель за пластиной, движущейся в облаке. Труды XLVII научной конференции МФТИ, часть VI, М.: Издательство МФТИ, 2004.
36. Здор А.Г. Динамика конденсированной фазы при движении плоской пластины в облаке. Труды XLVIII научной конференции МФТИ, часть VI, М.: Издательство МФТИ, 2005.
37. Здор А.Г. Эволюция массового спектра капель за телом, быстро движущимся в облаке. Материалы VI Международной конференции по неравновеснымпроцессам в соплах и струях (NPNJ 2006), 26 июня-1июля 2006 г., Санкт-Петербург. — М.: Вузовская книга, 2006.
38. Здор А.Г. Численное моделирование граничных условий взаимодействия JIA с атмосферными аэрозолями. Труды XLIX научной конференции МФТИ, часть VI, М.: Издательство МФТИ, 2006.
39. Здор А.Г., Миллер А.Б., Стасенко A.JI. Трансформация массового спектра капель облака, пересекаемого высокоскоростным крылатым JIA. Материалы XXX академических чтений по космонавтике, М.: Комиссия РАН, 2006.
40. Здор А.Г. Влияние на массовый спектр капель граничных условий на поверхности аппарата, движущегося в полидисперсной атмосфере. Материалы XXXI академических чтений по космонавтике, М.: Комиссия РАН, 2007.
41. Здор А.Г. Численные исследования моделей кинетики дисперсной фазы при обтекании высокоскоростного ДА. Материалы XXXII академических чтений по космонавтике, М.: Комиссия РАН, 2008.
42. Здор А.Г. Математическая модель и методические численные исследования эволюции атмосферного аэрозоля за пластиной под углом атаки. Труды ЦАГИ, выпуск 2676, 2008.
43. Здор А.Г. Физическая модель и численные исследования эволюции полидисперсного конденсата в следе за плоской пластиной, обтекаемой сверхзвуковым потоком. Принято к опубликованию в Учёные записки ЦАГИ, 2008.
44. Зилитинкевич С.С. Динамика пограничного слоя атмосферы. Ленинград: Гидрометеоиздат, 1970.
45. Карлсон Д., Хогланд Д. Сопротивление и теплопередача частиц в соплах ракетных двигателей. Ракетная техника и космонавтика, 1964, тю2, № 11, с. 104109.
46. Кашеваров A.B., Стасенко А.Л. Управление массовым составом и прозрачностью контрейла авиалайнера при помощи инжекции ионов в струи двигателей. Ученые записки ЦАГИ. 2006, т. XXVI, №3-4.
47. Кейн Э. Коэффициент лобового сопротивления шара при сверхзвуковой скорости и малых числах Яе. Вопросы ракетной техники. — М.: Изд-во иностранной литературы, 1953, вып. 2 (14), с. 70-86.
48. Коган М.Н. Динамика разреженного газа. М.: Наука, 1968.
49. Колмогоров А.Н. Уравнения турбулентного движения несжимаемой жидкости // Изв. АН СССР, сер. физ., 1942, т. 6, № 1-2.
50. Кошмаров Ю.А., Рыжов Ю.А. Прикладная динамика разреженного газа. М.: Машиностроение, 1977.
51. Крайко А.Н. К двухжидкостной модели течений газа с диспергированными в нем частицами // ПММ, 1982, т.6. Вып. 1, с. 96-106.
52. Лайхтман Д.Л. Физика пограничного слоя атмосферы. Л.: Гидрометеоиздат, 1970.
53. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. VI. Гидродинамика. -М.: «Наука», 1988.
54. ЛашковВ.А. Об экспериментальном определении коэффициентов восстановления скорости частиц потока газовзвеси при ударе о поверхность // ИФЖ. 1991, т.60, № 2, с. 197-203.
55. Лашков В.А. Аэродинамическое сопротивление цилиндра в двухфазном потоке // Изв. РАН, МЖГ. 1992, № 1, с. 123-129.
56. Лебедев П.Д., ЛеончикБ.И., Тыныбеков Е.К., МаякинВ.П., Лазарев В.Н. Исследование испарения капель в среде перегретого пара. ИФЖ, 1968, т. 16, № 4.
57. Левин Л.М. Исследования по физике грубодисперсных аэрозолей — М.: Изд-во АН СССР, 1961.
58. Леонтович М.А. Введение в термодинамику. Статистическая физика. М.: Наука, 1983.
59. Лойцянский Л.Г. Механика жидкости и газа. — М.: Наука, 1987.
60. Матвеев С.К. Модель газа из твердых частиц с учетом неупругих соударений //Изв. АН СССР, МЖГ. 1983, №6, с.12-16.
61. Матвеев С.К., Сеюкова Л.П. Обтекание сферы потоком газовзвеси // Динамика однородных и неоднородных сред. — Л.: Изд-во Ленингр. Ун-та, 1987. с. 16-23 («Газодинамика и теплообмен», Вып. 9).
62. Миллер А.Б., Моллесон Г.В., Стасенко A.JI. Механика и оптика сверхзвукового мелкодисперсного потока около освещаемой сферы // Учёные записки ЦАГИ, 2007, том XXXVIII, № 3-4, с. 92-101.
63. Монин A.C., Яглом A.M. Статистическая гидромеханика. Ч. 1. М.: Наука, 1965.
64. Найфэ А. Введение в методы возмущений. Пер. с англ. М.: Мир, 1984.
65. Нигматулин Р.И. Динамика многофазных сред. Т.1. М.: Наука, 1987.
66. Овсянников JI.B. Лекции по основам газовой динамики. М: Наука 1981.
67. Осипцов А.Н., Шапиро Е.Г. Влияние мелкодисперсной примеси на структуру пограничного слоя при гиперзвуковом обтекании затупленного тела // Изв. АН СССР, МЖГ, 1989, № 4, с. 85-92.
68. Осипцов А.Н., Шапиро Е.Г. Обтекание сферы запыленным газом с большой сверхзвуковой скоростью // Исследование газодинамики и теплообмена сложных течений однородных и многофазных сред / под ред. СтуловаВ.П. — М.: Изд-во Моск. Ун-та, 1990. с. 89-105.
69. Полежаев Ю.В. Проблемы теплообмена в запыленных потоках // Труды Второй Российской национальной конференции по теплообмену. Т. 1. Пленарные и общие проблемные доклады. Доклады на круглых столах. М.: Изд-во МЭИ, 1998. с. 64-69.
70. Полежаев Ю.В., Репин И.В., Михатулин Д.С. Теплообмен в сверхзвуковом гетерогенном потоке // ТВТ. 1992, т. 30, № 6, с. 1147-1153.
71. Резибуа П., Де Ленер М. Классическая кинетическая теория жидкостей и газов. Пер. с англ. М.: Мир, 1980.
72. Рёпке Г. Неравновесная статистическая механика. Пер. с нем. М.: Мир, 1990.
73. Салтанов Г.А. Неравновесные и нестационарные процессы в газодинамике. — М.: Наука, 1979.
74. Самарский A.A., Попов Ю.П. Разностные схемы газовой динамики. М.: Наука, 1975.
75. Сенковенко С.А., Стасенко А.Л. Релаксационные процессы в сверхзвуковых струях газа. -М.: Энергоатомиздат, 1985.
76. Стасенко A.JI. Газодисперсные течения в аэродинамике и летательной технике. Труды ЦАГИ, 1982, вып. 2133.
77. Стасенко А.Л. Модели динамики и теплообмена шаровых частиц в газодисперсных и капельных потоках// Труды ЦАГИ — 1983. Вып. 2220.
78. Стасенко А.Л. Модели дисперсных систем // Модели механики сплошной среды. Новосибирск: ИТПМ, 1983, с. 139-161.
79. Стасенко А.Л. Феноменология газодисперсных и парокапельных потоков с межфазным массообменом и лучистым переносом энергии. Труды ЦАГИ, 1994, вып. 2530, с. 3-27.
80. Стасенко А.Л. Проблемы авиационной экологии// Энергия, журн. Президиума РАН 1999. №7.
81. Стасенко А.Л., Толстых А.И., Широбоков Д.А. К моделированию оледенения самолёта: динамика капель и поверхность смачивания. РАН, Мат. моделирование, т. 13, №6, 2001.
82. Стасенко А.Л., Толстых А.И., Широбоков Д.А. Динамика деформируемых капель у поверхности крыла в вязком воздухе // Известия Академии наук, МЖГ, №5,2002, с. 180-190.
83. Стасенко А.Л., Здор А.Г., Миллер А.Б. Структура следа высотного летательного аппарата. Труды XLVI научной конференции МФТИ, часть VI, М.: Издательство МФТИ, 2003.
84. Стасенко А.Л. Физическая механика многофазных потоков. — М.: Издательство МФТИ, 2004.
85. Стернин Л.Е. Основы газодинамики двухфазных течений в соплах. М.: Машиностроение, 1974.
86. Стернин Л.Е., Маслов Б.Н., Шрайбер А.А., Подвысоцкий А.М. Двухфазные моно- и полидисперсные течения газа с частицами. М.: Машиностроение, 1980.
87. Стернин Л.Е., Шрайбер А.А. Многофазные течения газа с частицами. М.: Машиностроение 1994.
88. Фукс Н.А. Механика аэрозолей. М.: Изд-во АН СССР, 1955.
89. Хлопков Ю.И. Характеристики обтекания сферы в переходном режиме при сверх- и гиперзвуковых скоростях // Изв. АН СССР, МЖГ, 1981, № 3, с. 175-178.
90. ЦиркуновЮ.М. Влияние вязкого пограничного слоя на осаждение частиц при обтекании сферы газовзвесыо // Изв. АН СССР, МЖГ, 1982, № 1, с. 59-66.
91. Циркунов Ю.М., Тарасова Н.В. Влияние температуры преграды на осаждение тонкодисперсной примеси из сверхзвукового потока газовзвеси // ТВТ, 1992, №6, с. 1154-1162.
92. Циркунов Ю.М. Обтекание тел потоком газовзвеси. Докторская диссертация, С.-Пб., 2005.
93. Шафф С.А., Шамбре П.А. Течение разреженных газов // Основы газовой динамики. Изд-во иностранной литературы, 1963, с. 637-688.
94. Шидловский В.П. Введение в динамику разреженного газа. М.: Наука, 1965.
95. Яненко H.H., Солоухин Р.И., Папырин А.Н., Фомин В.Н. Сверхзвуковые двухфазные течения в условиях скоростной неравновесности частиц. -Новосибирск: Наука, 1980.
96. Aroesty J. Sphere drag in low-density supersonic flow // Raref. Gas Dynamics. — N.Y.-Lnd., Acad. Press, 1963, vol. 2, pp.261-277.
97. Coakley T.J. Turbulence modeling methods for the compressible Navier-Stokes equations // AIAA-83-1693.
98. Derevich I.V. Statistical modeling of particles relative motion in a turbulent gas flow // International Journal of Heat and Mass Transfer, No 49, 2006.
99. Derevich I.V. Coagulation kernel of particles in a turbulent gas flow // International Journal of Heat and Mass Transfer, No 50, 2007.
100. Henderson C.B. Drag coefficients of spheres in continuum and rarified flows. AIAA J., 1976, vol. 14, No 6, pp. 707-708.
101. Herner S.E. Fluid-dynamic drag. Published by the author, 1958.
102. Kavanau L.L. Heat transfer from spheres in a rarefied gas in subsonic flow // Trans. ASME. 1955, vol. 77, No 5, pp. 617-623.
103. Kogan M.N. Kinetic theory in aerothermodynamics. Prog. Aerospace Sei., Vol. 29, pp. 27-354, 1992. Printed in GB, Pergamon Press Ltd, 1993.
104. Mei R. An approximate expression for the shear lift force on a spherical particle on finite Reynolds number//Int. J. Multiphase Flow Vol. 18, No l,pp. 145-147, 1992.
105. Millikan R. A. The general law of fall of a small spherical body through a gas, and its bearing upon the nature of molecular reflection from surfaces// Phys. Rev. 1923. V. 22, No 1.
106. Nauman A. Luftwiederstand von Kugeln bei hohen Unterschallgeschwindigkeiten. Algem. Wärmetechnik, 1953, vol. 1, S. 217-221.
107. Osiptsov A.N. Mathematical modeling of dusty-gas boundary layers // Appl. Mech. Rev. 1997. vol. 50. No 6, pp. 357-370.
108. Pao Y.H. Structure of turbulent velocity and scalar fields at large wave numbers // Physics of Fluids, 1965, V.8, P. 1063-1075.
109. Sawford B.L. Reynolds number effects in Lagrangian stochastic models of turbulent dispersion// Phys. Fluids A 1991 - V. 3, No 6.
110. Sherman F.S. A survey of experimental results and methods for the transition regime of rarefied gas dynamics // Raref. Gas Dynamics. N.Y.-Lnd.: Acad. Press, 1963, vol. 2, pp. 228-260.
111. Sims W.H. Experimental sphere drag results in the near free molecular regime // Raref. Gas Dynamics. N.Y.-Lnd.: Acad. Press, 1969, vol. 1, pp. 751-756.
112. Stokes G.G. On the effect of internal friction of fluids on the motion of pendulums. Trans. Cambr. Phil., vol. 9, No 8.; Math. And Phys. Papers, Cambridge, 1901, vol. 3, pp. 1-141.
113. Corrsin S. Estimates of the relations between Eulerian and Lagrangian scales in large Reynolds number turbulens // J. Atmos. Sei., 1963, V. 5, No 2, pp. 115-119.
Обратите внимание, представленные выше научные тексты размещены для ознакомления и получены посредством распознавания оригинальных текстов диссертаций (OCR). В связи с чем, в них могут содержаться ошибки, связанные с несовершенством алгоритмов распознавания. В PDF файлах диссертаций и авторефератов, которые мы доставляем, подобных ошибок нет.