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

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

МЕХОНОШИНА Мария Андреевна

ПЕРЕНОС ТЕПЛА В СИЛЬНОНЕРАВНОВЕСНЫХ ТЕЧЕНИЯХ РЕАГИРУЮЩЕЙ СМЕСИ ГАЗОВ

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

АВТОРЕФЕРАТ

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

2 9 ИЮЛ 2015

Санкт-Петербург 2015

005571017

Работа выполнена в Санкт-Петербургском государственном университете.

Научный руководитель: доктор физико-математических наук,

доцент КУСТОВА Елена Владимировна

Официальные оппоненты: доктор физико-математических наук,

профессор Кузнецов Михаил Михайлович, Московский государственный областной универс профессор кафедры теоретической физики

кандидат физико-математических наук,

Бондарь Евгений Александрович,

Институт теоретической и прикладной механики

им. С.А. Христиановича СО РАН

зав. лабораторией вычислительной аэродинамик

Ведущая организация: Балтийский государственный технический

университет «Военмех» им. Д.Ф. Устинова (Санкт-Петербург)

Защита диссертации состоится " /»дстлф 2015 года в Л/ часов на заседании совета Д 212.232.30 по защите докторских и кандидатских диссертаций при Санкт-Петербургском государственном университете по адресу: 198504, Санкт-Петербург, Петродворец, Университетский пр., д. 28, математико-механический факультет, ауд. 405.

С диссертацией можно ознакомиться в Научной библиотеке Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская наб., д. 7/9 и на сайте http: / / spbu.ru/science/disser/soiskatelyu-uchjonoj-stepeni/dis-

list / details/14/563

«'//» иш

Автореферат разослан " '' " И"*71 2015 года.

Ученый секретарь диссертационного совета Е.В. Кустова

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

Актуальность темы. Изучение процессов переноса в высокотемпературных газах требуется для создания эффективной тепловой защиты космических аппаратов; описания процессов в ударных трубах, в рабочей части гиперзвуковых аэродинамических установок; использования в плазменно-химических технологиях. Поскольку при высоких скоростях и температурах экспериментальное определение транспортных свойств газа реализовать практически не удается, встает вопрос о корректном теоретическом моделировании диссипативных процессов в неравновесных условиях. Применение для этой цели хорошо апробированного аппарата кинетической теории представляется наиболее оптимальным. Кинетическая теория позволяет построить замкнутое описание течения и найти коэффициенты переноса при различных отклонениях от равновесия. Построение точных замкнутых моделей неравновесных высокотемпературных течений газов с внутренними степенями свободы является одной из главных задач данной работы.

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

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

Цель работы

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

2. Исследование коэффициентов переноса на основе точного кинетического подхода. Модификация формулы Эйкена для коэффициента теплопроводности в газах с учетом электронного возбуждения атомов и молекул.

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

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

Научная новизна

1. Учтено электронное возбуждение атомов и молекул в однотемпературном приближении, что впервые позволило обобщить формулу Эйкена на случай высокотемпературного неравновесного газа.

2. Проведены систематические расчеты коэффициентов теплопроводности компонентов воздуха N2, О2, N0, N,0 в широком диапазоне температур (200-20000 К); на основании анализа вклада возбужденных состояний атомов в интегралы столкновений установлены пределы применимости приближенных моделей расчета коэффициентов переноса в газе с электронным возбуждением.

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

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

5. Проведен систематический анализ вклада процессов теплопроводности, массовой диффузии, термодиффузии и диффузии колебательной

энергии в перенос тепла и массы в ударно-нагретых смесях N2^ и О2/О при различных начальных условиях; объяснены причины различного качественного поведения тепловых потоков в азоте и кислороде; обнаружен эффект компенсации тепловых потоков за счет различных диссипативных процессов.

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

Научная ценность и практическая значимость

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

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

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

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

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

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

2. Результаты расчета коэффициентов теплопроводности газов N2, О2, N0, Ы, О с электронным возбуждением в диапазоне температур 200 - 20000 К.

3. Результаты анализа вклада различных диссипативных процессов в тепловой поток за фронтом ударной волны в смесях N2/Ы, О2/О] эффект компенсации потоков за счет теплопроводности, массовой диффузии и диффузии колебательной энергии ведет к существенному

уменьшению полного теплового потока; термодиффузия не вносит заметного вклада в перенос тепла.

4. Сравнение потоков тепла в ударно-нагретых смесях Л^/ЛГ, Ог/О в однотемпературном и поуровневом приближениях выявило важную роль неравновесной колебательной кинетики в переносе тепла за ударными волнами.

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

Апробация работы. Основные результаты диссертации докладывались на следующих Всероссийских и международных конференциях: Международная конференция по механике "Шестые Поляховские чтения "(Санкт-Петербург, 2012); XXIII Всероссийский семинар с международным участием по струйным, отрывным и нестационарным течениям (Томск, 2012); 28-ой международный симпозиум по динамике разреженного газа (Испания, 2012); 5-я Европейская конференция по астронавтике и космическим наукам (Германия, 2013); IX Международная конференция по Неравновесным процессам в соплах и струях, №Ш'2014 (Алушта, 2014); 29-ый международный симпозиум по динамике разреженного газа (Китай, 2014); 8-я Всероссийская школа-семинар "Аэротермодинамика и физическая механика классических и квантовых систем" (Москва, 2014); Международная конференция по механике "Седьмые Поляховские чтения"(Санкт-Петербург, 2015). Результаты также докладывались на научных семинарах кафедры гидроаэромеханики Санкт-Петербургского государственного университета.

Публикации. Основные результаты диссертации опубликованы в работах [1]-[15], из них шесть ([1]-[6]) в журналах, входящих в перечень рецензируемых научных журналов, рекомендованных ВАК. Список работ приведен в конце автореферата. В работах [1,11] Кустовой Е.В. принадлежит общая постановка задачи и обобщение метода Энскога-Чепмена; автору принадлежит разработка алгоритма расчетов релаксационного давления и анализ результатов. В работах [2,9] Кустовой Е.В. принадлежит общая постановка задачи и обобщение метода Энскога-Чепмена; Макаркину Д.В. принадлежит вывод, расчет и анализ дифференциальных сечений столкновений; автору принадлежит разработка алгоритма численного расчета и анализ вклада релаксационного давления в диагональные члены тензора напряжений. В работах [3,5,12,14,15] Кустовой Е.В. принадлежит общая постановка задачи, разработка методов кинетической теории; Истомин В.А.

выполнил расчет и анализ коэффициентов сдвиговой и объемной вязкости, чисел Прандтля и Стокса; автору принадлежит разработка и реализация алгоритма программной части вычисления коэффициентов теплопроводности, определение фактора /¿nf в модифицированной поправке Эйке-на, анализ результатов. В работах [4,13] Кустовой Е.В. и Нагнибеда Е.А. принадлежит идея исследования, постановка задачи и обобщение метода Энскога-Чепмена, Куновой О.В. разработан алгоритм и произведен расчет макропараметров за фронтом ударной волны, автору принадлежит расчет коэффициентов теплопроводности и диффузии за фронтом ударной волны, расчет и анализ потоков диффузии и потоков тепла вблизи фронта. В работе [7] Куновой О.В. принадлежит расчет макропараметров за фронтом ударной волны в поуровневом приближении, автору принадлежит расчет и анализ потоков диффузии и потоков тепла за фронтом ударной волны. В работе [6] автору принадлежит расчет потоков диффузии и потоков тепла у поверхности затупленного тела в гиперзвуковом потоке; соавторами проведен численный расчет поля течения в двумерной задаче. В работе [8] Крылову A.A. и Пашкову В.А. принадлежит описание экспериментальных установок и постановка эксперимента, Кустовой Е.В. - теоретическая модель течения, автору принадлежит проведение расчетов коэффициентов переноса.

Структура и объем диссертации. Диссертация состоит из введения, 3 глав, заключения, приложения и списка литературы из 117 наименований. Общий объем диссертации составляет 112 страниц, включая 31 рисунок и 8 таблиц.

Содержание работы

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

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

Пусть тгар, tsi — характерные времена быстрых и медленных физико-химических процессов, происходящих при столкновениях частиц. Тогда, при условии тгар « ts[ ~ в (в — характерное время изменения макроско-

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

^ + и = (1)

£ = тгар1тв1 ~ Тгар/в « 1 — малый параметр, — соответственно

интегральные операторы быстрых и медленных процессов, /с/ — функция распределения частиц сорта с с внутренним состоянием I = (п, п, г, 3 — соответственно уровни электронной, колебательной и вращательной энергии, ис — скорость частиц сорта с. Приближенное решение уравнений (1) строится в виде обобщенного ряда Энскога-Чепмена по параметру с:

и (г, = £ £ТГс1 (и, р\, Чрх, У2рА, . . .) • (2)

г

Пространственная и временная зависимость коэффициентов этого ряда определяется макропараметрами газа р\ (г, ¿) и их градиентами всех порядков.

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

Р= (Р-Рге1) 1-2ф- СУ -VI, (3)

ч = -Л'УТ - р БТЛЫ + £ (+ {еЧ)го1 + ^ + еЛ пыУы,(4)

сг сг '

здесь р — давление, рте1 — релаксационное давление, г] — коэффициент сдвиговой вязкости, С — коэффициент объемной вязкости, Б — тензор скоростей сдвига, V — макроскопическая скорость газа, I — единичный тензор, Т — температура газа, А' — коэффициент теплопроводности поступательных и вращательных степеней свободы, (¡^ — диффузионная термодинамическая сила для частицы каждого сорта с на каждом колебательном уровне г, & — постоянная Больцмана, (еа)гог — осредненная вращательная энергия, е\ — колебательная энергия частиц сорта с, ес — энергия образования частиц сорта с, п„ — заселенность г-го колебательного уровня,

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

Va = -x;Dcidkddk - DTciVInT, (5)

dk

Dcidk — коэффициент диффузии для каждого химического сорта и колебательного уровня.

В однотемпературном приближении выражение для тензора напряжений формально совпадает с (3). Поток полной энергии имеет более простой вид:

q = —A'VT - р DTdc + ^ PchcVc, (6)

с с

hc — удельная энтальпия частиц сорта с. Поток характеризуется только градиентом температуры и градиентами концентраций химических компонентов, а от градиентов заселенностей колебательных уровней не зависит. Коэффициент теплопроводности А' в однотемпературном приближении описывает перенос всех видов внутренней энергии и определяется сечениями всех упругих и неупругих столкновений без химических реакций, коэффициенты диффузии Dai и термодиффузии Dtc в однотемпературном приближении не зависят от колебательных состояний.

Для определения коэффициентов переноса, входящих в выражения (3)-(5), с помощью метода Энскога-Чепмена из интегральных уравнений для функции распределения первого порядка выводятся замкнутые системы линейных алгебраических уравнений. Для решения этих систем необходим расчет интегральных скобок и интегралов столкновений (так называемых Г2^-интеграл ов).

В диссертации проанализированы модели для расчета О^-интегралов в различных диапазонах условий. Проведена оценка влияния размера возбужденных атомов на интегралы столкновений. Показано, что при температурах ниже 20000 К размер электронно возбужденных атомов слабо влияет на интегралы столкновений. Однако учет размера возбужденных атомов необходим при более высоких температурах.

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

В п. 2.1 коэффициент теплопроводности рассчитывался в газах N, О, N2, О2, NO с помощью точной кинетической теории с учетом возбужденных колебательных, вращательных и электронных степеней свободы, а также по приближенным формулам Эйкена и Гиршфельдера. В 1913 на осно-

ве представлений о длине свободного пробега молекул Е.Эйкен предложил разделять Л на две части (коэффициенты теплопроводности поступательных и внутренних степеней свободы):

А = Atr + A;nt = (Cc,tr/tr + Сс,int/int) Vi

здесь /tr = §, /и = 1 — безразмерные факторы, cCitr, cCiint — удельные теплоемкости поступательных и внутренних степеней свободы частиц. Дж.Гиршфельдер на основе метода Энскога-Чепмена предложил модификацию формулы Эйкена: /¡nt = ^ = 1.328, р - плотность, D - коэффициент самодиффузии.

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

Рис. 1. Фактор Эйкена для различных сортов как функция от Т

На основе точной кинетической теории были оценены факторы в формуле Эйкена. Среднее значение для /1г равно | во всем диапазоне рассматриваемых температур. На рис.1 представлены значения фактора внутренних степеней свободы /ы для всех рассмотренных химических сортов. Для молекул при низких температурах полученный фактор совпадает со значением, предложенным Гиршфельдером. При температурах выше 5000 К среднее значение фактора для молекул стремится к константе /¡^ = 1.40. В среднем диапазоне температур значение фактора зависит от химического

сорта и температуры и может быть аппроксимировано простыми формулами.

Для атомов без учета электронного возбуждения /¡п1 = 0. При учете электронных степеней свободы в настоящей работе обнаружено, что при температурах выше 3000 К среднее значение /¡п1 для атомов также стремится к константе, несколько большей, чем для молекул, /¡пЬ = 1.42. При низких температурах теплоемкость внутренних степеней свободы атомов близка к нулю, поэтому значение /¡п4 не влияет на теплопроводность, и можно использовать предложенную константу во всем диапазоне температур.

Таким образом, для инженерных расчетов коэффициентов теплопроводности мы можем предложить значения фактора /¡пЬ для всего рассмотренного диапазона температур, от 200 до 20000 К. Значения /¡п4 приведены в таблице 1.

Таблица 1. Значение поправки /¡^ для формулы Эйкена.

с Т < 20000 К

N 1.42

О 1.42

с Т < 1000 к 1000 < Т < 5000 К 5000 < Т < 20000 К

N0 1.328 -2 • ИГТ2 + 3 • 10-Т + 1.2988 1.40

о2 1.328 -3 ■ 10"9Г2 + 3 • 10"5Г + 1.3436 1.40

N2 1.328 -8 • 10"9Т2 + 6 • 10"5Т + 1.2783 1.40

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

Оценен вклад релаксационного давления в диагональные члены тензора напряжений. Показано, что без учета электронного возбуждения вклад релаксационного давления не превышает 8%, в то время как при учете электронного возбуждения может достигать 14% при концентрации атомов азота, равной п^/п = 0.5. Интересно отметить, что знак релаксационного давления определяется направлением протекания химической реакции.

В п. 2.3 исследованы скорости диффузии и поток тепла в бинарной смеси газов за фронтом ударной волны в однотемпературном приближении

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

бинарной смеси атомов и молекул:

^ + + = (7)

аъ

^ + ПаУ • V + У • (ПаУа) = -2(8) аъ

Р"ТГ + V • Р = 0, (9)

аъ

р^ + У- Ч + Р:Уу = 0, (Ю)

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

Ус = У™В + У с = а,т (11)

-уШЭ _ _ Г) А _ Т-) н » т — ^таиа пи

дгМ!) _ _ Г) ^ _ п Л

V™ = -ДгсУ1пТ, с = а,т

V™ — соответственно вклады массовой диффузии и термодиффузии.

4 = + + (12)

здесь и ^^ потоки тепла за счет теплопроводности (поток Фу-

рье) поступательных и всех внутренних степеней свободы, массовой диффузии и термодиффузии, которые определяются следующим образом:

чнс = —А'УТ,

ЧМ° = РшЬтУ^ + РаК (13)

ЧТ° = -р(0Ттйт + 0ТаАа) + ртКУ™ + раНаУ1°.

Апт, -Ста) Аш, Агт, Е>Та ~ бинарные коэффициенты диффузии и коэффициенты термодиффузии, с!то, с1а - диффузионные термодинамические силы, Нт, Ъ,а — удельные энтальпии молекул и атомов. Такое представление удобно для анализа вклада различных диссипативных процессов в перенос тепла.

Система уравнений (7) - (10) решалась в одномерной стационарной постановке для смесей ЛГ2/N и О2/О. В набегающем потоке брались следую-

щие условия: Т0 = 271 К, р0 = 100 Па, М0 = 10; 15; 18, в невозмущенном потоке атомы отсутствовали. Параметры газа за фронтом ударной волны вычислялись с помощью соотношений Рэнкина-Гюгонио при условии замороженного химического состава смеси и полностью возбужденных внутренних степеней свободы во фронте. Показано, что кинетические и диссипа-тивные процессы протекают более эффективно при высоких числах Маха. Исследована роль теплопроводности, термодиффузии и массовой диффузии в переносе тепла; обнаружено, что поведение полного потока тепла в основном определяется массовой диффузией. Все процессы в смеси О^/О протекают заметно быстрее, чем в Л^/Л/", это ведет к более высоким значениям градиентов макропараметров и, следовательно, к более интенсивному переносу массы и энергии.

В третьей главе изучено влияние поуровневой кинетики на диффузию и перенос тепла за фронтом ударной волны в бинарных смесях N2^ и О2/О (электронное возбуждение не учитывалось). Решалась система уравнений за фронтом ударной волны для заселенностей колебательных уровней щ, числовых плотностей атомов па, скорости у и температуры Т, поток предполагался одномерным и стационарным, газ — невязким:

¿(ущ) _ ЫЬг пМзе-гес

= ДГГ + Я?"-™0, г = 0 ..Ь, (14)

йх (15)

йх й{уПа)

р0у0 = ру, (16)

Л)«о+Ро = РЪ2+Р, (17)

2 2

+ ^ = Л + (18)

здесь х — расстояние от фронта ударной волны, ро, г>о, ко — плотность, скорость и энтальпия газа в набегающем потоке, Ь — число возбужденных колебательных уровней. Правые части уравнений содержат поуров-невые коэффициенты скорости УТ и УУ переходов колебательной энергии, диссоциации и рекомбинации. Вероятности колебательных переходов вычислялись на основе обобщенной БЗН-теории для ангармонических осцилляторов, для описания диссоциации использовалась модель Тринора-Маррона, модифицированная для поуровневого приближения. Для исследования процессов переноса применялся приближенный подход: макропараметры и их производные, полученные в результате решения уравнений (14)—(18), подставлялись в выражения для потоковых членов, полученных в первом приближении метода Энскога-Чепмена.

Характерной особенностью диффузии и теплопереноса в поуровневом

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

Уг = УР + V?10 + (19)

Уа = У™ + УГ- (20)

Здесь, У-1'-0, V™ Vа° и У?УЕ - соответственно вклад массовой

диффузии, термодиффузии и диффузии колебательной энергии:

утМО _ _Г) ^ _ п я

УР = -1?г^1пТ, (21)

V = У1п—,

\ "гот и та / пт

= - АЛ - Оаайа, V™ = -ВТаЧ 1п Т,

ч=чНС + (1МП + чТП + чОУЕ (22)

qБУЕ _ П0Т0К за счет диффузии колебательной энергии:

Ч°УЕ = £ (|*Г + (е')гог + щУГЕ- (23)

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

Сначала рассмотрен случай равновесного набегающего потока (То = Ту, Ту — колебательная температура ); проведена оценка влияния начальных условий на перенос тепла за ударной волной при следующих начальных условиях: Т0 = 271 К, р0 = 100 Па, пт = р0/кТ0, па = 0, М0 = 10,15,18. На рис. 2 приведены тепловые потоки за фронтом ударной волны как функции расстояния х от фронта для смесей -Л^/ЛГ и О2/О. При Мз = 10 поток тепла близок к нулю, так как градиенты всех макропараметров малы. В смеси N2^ при Мо = 15 полный поток тепла быстро меняется вблизи фронта ударной волны, затем стремится к нулю при установлении равновесного состояния. Самый интересный случай соответствует смеси Л/г/ЛГ при Мо = 18 и смеси О2/О при Мэ = 15 и Мо = 18. Полный поток тепла ведет себя немонотонно, так как основной вклад в полный поток тепла в зависимости от х вносят различные процессы. Показано, что для более высоких чисел Маха кинетические и диссипативные процессы происходят эффективнее.

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

д, кВ/м2

X, см

(а)

(Ь)

м0=ю

М0=15 М„=18

Рис. 2. Полный поток тепла за ударной волной как функция от х. 2о = Т„. (а) - Л^/ЛГ, (Ь) - Ог/О.

(а)

(Ь)

Рис. 3. Вклад различных процессов в поток тепла за фронтом ударной волны как функция от х. Т0 = Т,. Л/о = 18. (а) - ЛГ2/ДГ, (Ь) - 02/0.

потока при Мо = 18. Видно, что термодиффузия практически не влияет на поток тепла за фронтом ударной волны. В смеси Л^/ЛГ (рис. 3(а)) поток Фурье и поток за счет диффузии колебательной энергии вносят вклад одного порядка сразу за фронтом ударной волны. Их абсолютные значения велики при х < 0.1 см. Знаки этих потоков противоположны, поэтому происходит компенсация, и значение полного потока оказывается значительно меньше. Вклад потока за счет диффузии колебательной энергии оказывается больше, что ведет к тому, что знак полного потока тепла противоположен знаку потока Фурье. Это объясняет, почему полный поток тепла - отрицательный. Влияние массовой диффузии становится значительным лишь при х > 0.05 см из-за задержки реакции диссоциации. Это приводит к повышению полного потока тепла при х > 0.05 см. Поэтому при Мо = 18 поведение полного потока тепла является немонотонным. При более низких числах Маха диссоциация молекул азота недостаточна, чтобы оказать

эффект на полный поток тепла. В смеси 02/0 все процессы протекают быстрее и градиенты газодинамических параметров больше по сравнению со случаем N2/./V, что приводит к большим значениям потока; диссоциация в смеси О2/О начинается раньше и протекает активнее, поэтому влияние массовой диффузии на полный поток тепла сильнее, чем в смеси N2/N.

Далее рассмотрен случай неравновесного набегающего потока, когда колебательная температура % = 8000 К. В этом случае (рис. 4) в смеси Л^/ЛГ полный поток тепла меняет свой знак и становится положительным. Для смеси О2/0 поток не меняет знак, но его абсолютное значение становится на порядок больше по сравнению со случаем равновесного набегающего потока. Чтобы понять поведение полного потока тепла, проанализируем вклад теплопроводности, термодиффузии, массовой диффузии и диффузии за счет переноса колебательной энергии в q. Для смеси азота (рис. 5(а)) влияние массовой диффузии на полный поток тепла мало, а знак д определяется дВ]/Е. В случае сильной неравновесности (То = 271 К, Ту = 8000 К), полный тепловой поток, поток Фурье и поток за счет диффузии колебательной энергии меняют знаки, так как процессы возбуждения заменяются деактивацией, а температура газа возрастает с ростом х. В случае кислорода (рис. 5(Ь)) начальная колебательная неравновесность ускоряет диссоциацию, поэтому основной вклад в полный поток тепла вносит массовая диффузия. Несмотря на то, что поток Фурье и поток за счет диффузии колебательной энергии меняют знаки, знак полного потока тепла в случае сильной неравновесности определяется цмв, абсолютное значение которого значительно выше, чем значения и цОУЕ. Поэтому полный тепловой поток в смеси О2/О остается отрицательным.

(а) (Ь)

Рис. 4. Полный поток тепла смеси за ударной волной как функция от х. Мо = 10. (а) - ЛГ2/ЛГ, (Ь) •

Ог/О.

кВ/м

,Т,-Т,-271

-*-ч"',Т,=Т1=271 п-ц^. т_=яс;ю

о - Т¥=8000

л-ч"°,т,=8000

-«-Ч~.Т,=Т,=271

- □ - ч**, Т,=8000 -О-ч'^'.ТевООО

- Л - Ч**, Т *8000

(Ь)

Рис. 5. Вклад различных процессов в поток тепла за фронтом ударной волны как функция от х.

М0 = 10. (а) - Nг|N, (Ь) - 02/0.

ч, кВ/м

-БТБ ■ 1-Т

0.3

д, кВ/м

0.0

-ЭТЭ

---1-Т

0.1

Рис. 6,

х, см х, см

(а) (Ь)

Полный поток тепла за ударной волной как функция от х. М0 = 18. (а) - /Л', (Ь) - О2/О.

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

Заключение содержит список основных результатов, полученных в работе.

1. В диссертации представлены результаты расчета коэффициентов пе-

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

2. На основе однотемпературного приближения метода Энскога-Чепмена рассчитан коэффициент теплопроводности молекул N2, О2, N0 с возбужденными колебательными, вращательными и электронными степенями свободы и электронно возбужденных атомов Аг, О. Предложена модификация поправки Эйкена для коэффициента теплопроводности в температурном диапазоне 200-20000 К, учитывающая электронное возбуждение молекул и атомов. Проведена оценка влияния размера возбужденных атомов на интегралы столкновений. Установлены пределы применимости приближенных моделей расчета коэффициентов теплопроводности при учете электронного возбуждения.

3. Изучено влияние химического состава, температуры и учета электронного возбуждения атомов и молекул на релаксационное давление в смеси N2/И при температурах 3000-10000 К; выявлена зависимость величины релаксационного давления от типа реакции и степени отклонения от равновесия, оценен вклад релаксационного давления в диагональные члены тензора напряжений.

4. Диффузия и поток тепла при течении бинарных смесей Л^/ЛГ и О2/О за фронтом ударной волны исследованы для различных чисел Маха в набегающем потоке (М0 = 10; 15; 18) с использованием различных кинетических моделей. Проведено сравнение результатов, полученных при использовании поуровневого и однотемпературного подходов. Показано, что в однотемпературном приближении сразу за фронтом ударной волны наблюдается существенное завышение числовых плотностей атомов, что ведет к переоценке вклада массовой диффузии и, как следствие, завышенным значениям полного потока тепла.

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

6. Показано, что в случае неравновесного набегающего потока при Мо = 10 в смеси N2^ с начальным возбуждением колебательных степеней свободы тепловой поток меняет знак; в смеси О2/О происходит заметный рост теплового потока.

Работы автора по теме диссертации

Публикации в журналах, рекомендованных ВАК:

[1] Кустова Е.В., Мехоношина М.А. Релаксационное давление в смеси N2/N с учетом неравновесной реакции диссоциации // ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА. СЕРИЯ 1: МАТЕМАТИКА, МЕХАНИКА, АСТРОНОМИЯ, 2012. № 1. С. 86-95

[2] Kustova E.V., Makarkin D.V, Mekhonoshina M.A. Normal Mean Stress in Non-Equilibrium Viscous N2/N Flows with Dissociation and Electronic Excitation // AIP Conference Proceedings, 2012. V. 1501, P. 1086-1093.

[3] Istomin V.A., Kustova E. V, Mekhonoshina M.A. Eucken correction in high-temperature gases with electronic excitation // Journal of Chemical Physics, 2014. V. 140, № 18. P. 184311.

[4] Kunova O.V., Kustova E.V., Mekhonoshina M.A., Nagnibeda E.A. The influence of state-to-state kinetics on diffusion and heat transfer behind shock waves // AIP Conference Proceedings, 2014. V. 1628, P. 1202-1209.

[5] Istomin V.A., Kustova E. V., Mekhonoshina M.A. Validity of Eucken formula and Stokes' viscosity relation in high-temperature electronically excited gases // AIP Conference Proceedings, 2014. V. 1628, P. 12291236.

[6] Josyula E., Burt J., Kustova E., Vedula P., Mekhonoshina M. State-to-State Kinetic Modeling of Dissociating and Radiating Hypersonic Flows // AIAA Aerospace Sciences Meeting, 2015. AIAA 2015-0475.

Другие публикации:

[7] Кунова О.В., Мехоношина М.А. Влияние неравновесной кинетики на теплоперенос и диффузию за фронтом ударной волны // ФИЗИКО-ХИМИЧЕСКАЯ КИНЕТИКА В ГАЗОВОЙ ДИНАМИКЕ, 2014. Т. 15. http://chemphys.edu.ru/issues/2014-15-l/articles/106/

[8] Kustova Е., Krylov A., Lashkov V., Mekhonoshina М. Experimental Facilities and Modelling for Rarefied Aerodynamics // RTO-EN-AVT-194 - Models and Computational Methods for Rarefied Flows. — von Karman Institute, Rhode St. Gense, Belgium: NATO Science and Technology Organization, 2011. 624, 20 стр.

[9] Kustova E.V., Makarkin D.V, Mekhonoshina M.A. Normal Mean Stress in Non-equilibrium Viscous N2/N Flows with Dissociation and Electronic Excitation // 28th Int. Symp. on Rarefied Gas Dynamics, Book of Abstract, Zaragoza, Spain, July 9-13, 2012 — Zaragoza, 2012. P. 134-135

[10] Мехоношина М.А. Влияние неравновесной реакции диссоциации на нормальные напряжение в смеси N2-N с электронным возбуждением // XXIII семинар по струйным, отрывным и нестационарным течениям: сборник трудов / под ред. Г.В. Кузнецова и др.; Томский политехнический университет, г. Томск, 2012. С. 430, С. 250-251.

[11] Кустова Е.В., Мехоношина М.А. Релаксационное давление в смеси N2/N с неравновесной реакцией диссоциации-рекомбинации и электронным возбуждением // Шестые Поляховские чтения: Тезисы докладов Международной научной конференции по механике, Санкт-Петербург, 31 января-3 февраля 2012г. г. Санкт-Петербург, 2012. С. 304

[12] Istomin V.A., Kustova E.V., Mekhonoshina М.А. Heat flux and stress tensor in plasma flow with electronic excitation // Proceedings of 5th European Conference for Aeronautics and Space Science. June 2013, Munich, Germany, 2013.

[13] Kunova O.V., Kustova E.V., Mekhonoshina M.A., Nagnibeda E.A. The influence of state-to-state kinetics on diffusion and heat transfer behind shock waves // Abstract book of 29th International Symposium on Rarefied Gas Dynamics. Xi'an, China, 2014. - P. 287.

[14] Istomin V.A., Kustova E.V., Mekhonoshina M.A. Validity of Eucken formula and Stokes viscosity relation in high-temperature electronically excited gases // Abstract book of 29th International Symposium on Rarefied Gas Dynamics. Xi'an, China, 2014. P. 283-284.

[15] Истомин В.А., Кустова E.B., Мехоношина М.А. Теплопроводность высокотемпературных газов с электронным возбуждением: уточнение поправок Эйкена и Гиршфельдера // Материалы X Международной конференции по Неравновесным процессам в соплах и струях (NPNJ'2014), 2014.

Подписано в печать 10.07.2015. Формат 60 х 84 V16. Бумага офсетная. Гарнитура Times. Печать цифровая. Усл. печ. л. 1,00. Тираж 100 экз. Заказ № 026.

Отпечатано в Издательстве ВВМ. 198095, Санкт-Петербург, ул. Шзецова, 41.