Авторефераты по всем темам  >>  Авторефераты по разным специальностям САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

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

ПАННЕВИЦ Оксана Владимировна

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

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

АВТОРЕФЕРАТ

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

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

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

Научный консультант: доктор физико-математических наук, профессор РЫДАЛЕВСКАЯ Мария Александровна

Официальные оппоненты: доктор физико-математических наук, профессор КУЗНЕЦОВ Михаил Михайлович (Московский государственный областной университет) доктор физико-математических наук, профессор ХАНТУЛЕВА Татьяна Александровна (Санкт-Петербургский государственный университет)

Ведущая организация: Балтийский государственный технический университет ВОЕНМЕХ им. Д.Ф. Устинова

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

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

М. Горького Санкт-Петербургского государственного университета по адресу: 199034, Санкт-Петербург, Университетская наб., 7/9.

Автореферат разослан У____Ф ___________ 20__ г.

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

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

.

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

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

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

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

Результаты, выносимые на защиту

:

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

2. Исследование разных стадий пространственно-однородной релаксации азота и кислорода.

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

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

Доказательство применимости метода Ньютона для определения этих параметров.

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

Апробация результатов. Результаты исследований неоднократно докладывались на семинарах кафедры гидроаэромеханики Санкт-Петербургского государственного университета; на Международных научных конференциях по механике Поляховские Чтения: Третьи (2003г), Четвертые (2006г), Пятые (2009);

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

Статьи [2,3,4] и доклады [5,6,7,9] из этого списка подготовлены вместе с научным руководителем М.А. Рыдалевской. В [2] соискателю принадлежит осуществление послойного исследования релаксационной зоны за прямым скачком уплотнения в диссоциирующем азоте (руководителю Ч постановка задачи о разбиении релаксационной зоны за фронтом ударной волны на несколько слоев). В [3] соискателю принадлежит доказательство применимости метода Ньютона для определения равновесного состояния газа за прямым скачком уплотнения и численная реализация этого метода (руководителю Ч постановка задачи и предложение перейти от экстенсивных к сопряженным интенсивным параметрам).

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

Структура и объем работы. Диссертация состоит из введения, 4-х глав и заключения. Каждой главе ставится в соответствие приложение, в котором дается численная иллюстрация полученных в главе результатов. Основной текст диссертации изложен на 65 страницах, список литературы включает наименований. Приложения, представленные на 49 страницах, содержат 1 таблицу и 95 рисунков.

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

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

Глава 1 посвящена кинетическому описанию высокотемпературных газов из гомоядерных молекул. Рассматриваются газы с вращательными и колебательными степенями свободы, процессами диссоциации и рекомбинации. Поступательные степени свободы атомов и молекул описываются квазиклассически. Уровни вращательной энергии заранее не конкретизируются, а вращательные статистические суммы рассчитываются в классическом пределе. Для описания колебательного движения атомов в молекуле используется модель ангармонического осциллятора (при этом учитывается, что колебательная энергия не должна превосходить энергию диссоциации). Кинетические уравнения для определения функций распределения свободных атомов и молекул в рассматриваемых условиях записываются в безразмерном виде [1,3] 1 D f1 = + , D f2 r = + , (1) 1 1 1 2 r 2 r 2 r где f1 и f2r - функции распределения атомов и молекул с фиксированными значениями вращательной 2r и колебательной 2 энергии ( r = 0, r, r - число возбужденных вращательных уровней, = 0,d, d - уровень колебательной энергии, соответствующий энергии диссоциации молекулы d2 ); D1 и D2r - дифференциальные операторы; J1 и J2r - интегральные операторы, соответствующие наиболее частым столкновениям, определяющим быстрые (по сравнению с газодинамическими) молекулярные процессы ( - среднее время между этими столкновениями, t*, t* - характерное время изменения макропараметров); J1 и J2r - интегральные операторы, соответствующие сравнительно редким столкновениям (эти столкновения характеризуют медленные молекулярные процессы); = / t*. Рассматриваемые течения соответствуют разным релаксационным режимам в зависимости от того, какие столкновения описывают операторы J1 и J2r.

Согласно известным данным (см., например, [2]), вероятность колебательных обменов возрастает (а среднее время убывает) более чем на порядок при уменьшении в 2 раза относительного дефекта резонанса колебательной энергии:

= 2 ' + 2 ' -2 -2 / 2 + 2. (2) ( ) 11 В [1] было предложено рассматривать вслед за поступательно-вращательной несколько последовательных стадий колебательной релаксации, когда к быстрым процессам относятся колебательные переходы, при которых величина не превосходит 1/8, 1/4 или 1/2, а затем (на стадии перехода в состояние равновесия) - все обмены колебательной энергией. Такой подход позволил теоретически обосновать установление некоторых молекулярных распределений по колебательным энергиям, которые ранее были обнаружены экспериментально.

Так как вероятность химической перестройки молекул также зависит от относительного дефекта резонанса внутренней энергии, в настоящей работе, по аналогии с (2), вводится в рассмотрение величина max - 2 / 2 (3) = 2 и предлагается рассматривать стадии релаксации, при которых к быстрым процессам относятся наряду с колебательными переходами столкновения, приводящие к распаду молекул или слиянию атомов, если и изменяются в одинаковых пределах. Расчеты, проведенные для ряда двухатомных газов, показали, что соответствующие ограничения сверху величины (3) выполнены, если колебательные уровни превосходят некоторое значение > 0. При 1/ 8 - > d, при 1/ 4 - d, а при 1/ 2 - < d.

Это позволило ввести в рассмотрение пять стадий релаксации:

1)поступательно-вращательную; 2)начальную колебательную (с учетом лишь тех колебательных переходов, при которых 1/ 8); 3)колебательно-химическую (с учетом колебательных обменов при 1/ 4 и процессов диссоциации с верхнего уровня, что соответствует лестничной модели); 4)еще одну колебательнохимическую (с учетом колебательных обменов при 1/ 2 и диссоциации с ряда верхних уровней d, что в каком-то смысле соответствует модели ТринораМаррона); 5)завершающую стадию перехода к локальному равновесию.

На каждой стадии релаксации ведущие столкновительные операторы J1 и J2r аннулируют распределения [1]:

3 m s2rm2 ( ) ( ) f1(0 ) = exp , f2(0) = exp (4) 1 r h h3 =0 2r = 0 ( ) Здесь 1() и - аддитивные инварианты столкновений, описываемых этими 2r операторами (разным стадиям релаксации соответствуют различные наборы инвариантов); - интенсивные макропараметры, которые могут зависеть от координат и времени (их число на каждой стадии релаксации равно числу аддитивных инвариантов столкновений, учитываемых на этой стадии); h - постоянная Планка; m1 и m2 - массы атомов и молекул; s2r -- статистический вес r -го вращательного уровня молекулы s2 = 1 ; поступательная энергия частиц ( ) представлена в собственных скоростях c = u - v, u - их скорость в неподвижной системе координат, v - газодинамическая скорость.

Функции распределения (4), нормированы на суммарные значения инвариантов в единице объема. При известных аддитивных инвариантах они позволяют исследовать колебательные заселенности n2 и химический состав газа (числа свободных атомов n1 и молекул n2 ).

В данной главе для каждой из пяти перечисленных выше стадий релаксации приведены аддитивные инварианты столкновений и исследованы температурные зависимости относительных колебательных заселенностей x2 = n2 / n20, а также относительных химических концентраций свободных атомов n1 / n и молекул n2 / n ( n = n1 + n2 ) в диссоциирующих азоте и кислороде. Результаты исследований приведены в Приложении 1 на 11 рис. При этом коэффициенты 0 (при ( ) аддитивных инвариантах , содержащих поступательную энергию mc2 / 2 ) и ( ) (при дополнительных энергетических инвариантах , которые являются обобщением инварианта Тринора и представляют собой кусочно-линейные аппроксимации колебательной энергии 2 ) представлены в виде:

=-1/ kT, = 1/ kT kT1, (5) ( )-1/ ( ) 0 где T - температура газа, а T1Ц так называемая температура 1-го колебательного уровня молекул [5].

Поведение относительных колебательных заселенностей и концентраций атомов и молекул смеси N + N2 на 4-й стадии релаксации иллюстрируют рис. 1 и 2, где приведены величины y2v = ln x2v /(21 / kT1 ) и относительные концентрации ( ) n1 / n, n2 / n при начальной плотности молекул N2, равной числу Лошмидта nL.

0 4 8 12 16 20 24 28 32 36 40 v -------- y2v Рис. 1. Относительные колебательные заселенности молекул N2 на стадии 4 (колебательно-химической релаксации).

Кривая 1 соответствует T1 /T = 1; 2 - T1 /T = 0,1; 3 - T1 /T = 8.

1,1,1,1,0,0,4' 1' 0,0,3' 0 T 500 3000 5500 8000 10500 13000 15500 180Рис. 2. Относительные концентрации для атомов N и молекул Nна стадии 4 (колебательно-химической релаксации).

Кривые 1 и 1Т соответствуют отношению T1 /T = 1;

3 и 3Т - T1 /T = 10 ; 4 и 4Т - T1 /T = 0,5.

Расчеты показали, что относительные колебательные заселенности x2 на стадиях релаксации 2Ц4 могут сильно отличаться от больцмановских. Величины y2v = ln x2v /(21 / kT1 ) зависят лишь от отношения T1 / T. Отношение T1 /T = ( ) соответствует состоянию равновесия. При T1 /T < 1 величины y2v меньше равновесных, но очень близки к ним. При T1 /T > 1 их поведение на нижних колебательных уровнях совпадает с поведением распределения Тринора [5], на верхних - с поведением распределения Больцмана при температуре T, а на средних - их значения тем больше, чем больше T1 / T. На стадии 4, если T1 / T > 5, наблюдается характерный платообразный участок, о котором много говорилось в литературе. В настоящей работе показано, что наличие атомарной компоненты (учет диссоциации и рекомбинации) оказывает влияние лишь на абсолютные значения колебательных заселенностей и не сказываются на значениях относительных заселенностей x2. На стадиях 2Ц4 последние совпадают с соответствующими заселенностями в химически однородном газе, приведенными в монографии [1].

Из рис. 2 видно, что увеличение T1 / T, которое оказывает существенное влияние на возбуждение колебаний, влияет на относительные концентрации химических компонентов (свободных атомов и молекул) и приводит к повышению степени диссоциации.

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

Изначально газ находился в состоянии равновесия с температурой T0. Затем газ резко вывели из состояния равновесия за счет какого-то краткосрочного воздействия таким образом, что в газе установилось максвелл-больцмановское распределение по поступательным и вращательным степеням свободы с температурой T при сохранении его химического состава и колебательных заселенностей молекул. Его неравновесное начальное состояние соответствует завершению 1-й стадии RT-релаксации с функциями распределения 3 ( n1(0) m1 m1c2 n20) s2rm2 m2c2 / 2 + 2r f1(1) = exp exp, (6) - , f2(1) = - Z1 T h3 2kT r Z2 RT T h3 kT ( ) ( ) где Z1 T и Z2 RT T - статистические суммы атомов по поступательным и ( ) ( ) ( молекул по поступательным и вращательным степеням свободы, n1(0) и n20) соответствуют начальному состоянию газа (распределениям (4) в равновесных условиях). После прекращения внешнего воздействия на систему она будет релаксировать к состоянию равновесия.

Согласно предложенной в первой главе схеме разделения на стадии релаксационного процесса, в настоящей главе исследуются состояния газа на этапах завершения каждой стадии. При этом на каждом этапе в газе устанавливаются функции распределения (4). В показателях экспонент содержатся аддитивные инварианты столкновений, которые описываются операторами J1 и J2r на данной стадии. Неизвестные коэффициенты при этих инвариантах определяются из условий сохранения 0 0 1 1 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) f1 1 dc + f2 vr dc = f1 1 dc + f2vr dc. (7) 2 vr 2vr v,r v,r В левых частях равенств (7) стоят функции распределения (4) с неизвестными коэффициентами , в правых частях - известные функции (6). На каждой стадии релаксации число условий сохранения совпадает с числом инвариантов учитываемых столкновений. Поэтому на каждой стадии число уравнений (7) равно числу неизвестных параметров . Системы алгебраических уравнений, к которым сводятся соотношения (7) после вычисления соответствующих интегралов, решались методом Ньютона.

В настоящей главе проведено последовательное исследование разных стадий пространственно-однородной релаксации смесей N2 + N и O2 + O (в частности, поведение температур T и T1, введенных в (5)). Результаты соответствующих исследований приведены в Приложении 2 на рис.12Ц25.

Расчеты показали, что на этапах начальной колебательной и колебательнохимической релаксации температура газа T мало меняется и остается близкой к температуре T (независимо от температуры T0 ). Температура первого колебательного уровня T1 возрастает с увеличением T0. При достаточно больших значениях T0 температура T1 возрастает с увеличением T0 тем больше, чем ниже температура T (см. рис. 3). В равновесной ситуации, как и следовало ожидать, температура газа T тем выше, чем больше температуры T0 и T.

T160140120100806040201000 2000 3000 4000 5000 6000 7000 8000 9000 10000 11000 120TРис. 3. Температура T1 на стадии 4 (колебательно-химической релаксации) в газовой смеси N + N2.

Кривые 1, 2, 3 соответствуют значениямT : 1000 K, 6000 K, 11000 K.

Глава 3 посвящена описанию течений диссоциирующего двухатомного газа на разных стадиях релаксации. Газодинамические уравнения выводятся из систем кинетических уравнений (1). В ситуации, когда при вычислении потоковых членов на каждой стадии релаксации используются функции распределения (4), эти уравнения соответствуют приближению идеальной жидкости и имеют вид dv = F - р, (8) dt % d de %% = = -( ) e + p div v + e, (9) dt dt d =- div v +, = 1,...,. (10) dt Здесь F - сила, действующая на единицу массы газа; - массовая плотность % газа; p - давление; e - энергия единицы объема газа, включающая наряду с поступательной энергией молекул ту часть внутренней энергии, которая обменивается с поступательной при столкновениях, описываемых операторами J1 и J2r ; -- суммарные значения аддитивных инвариантов этих операторов;

% % e и - изменение параметров e и ( =1,...,) за счет столкновений, определяемых операторами J1 и J2r. Часть уравнений (10) соответствует сохранению числа частиц, которые не изменяются при столкновениях, которые на данной стадии релаксации описываются операторами J1 и J2r. Следствием этих уравнений является уравнение неразрывности.

В настоящей главе исследовались свойства уравнений (8)Ц(10) для разных релаксационных режимов, перечисленных выше. Особое внимание уделялось изоэнтропийным течениям, когда медленные процессы % заморожены ( e и равны нулю).

Одним из наиболее важных газодинамических параметров является скорость звука. Поведение этого параметра в разных газах с физикохимическими процессами и его зависимость от интенсивных параметров в ситуациях, когда для функций распределения справедливы выражения, подобные (4), исследовалось в [4]. Там было показано, что во всех этих случаях сохраняется традиционное выражение a2 = p / для квадрата скорости звука, и выведена аналитическая формула, определяющая зависимость коэффициента от параметров .

В настоящей главе формула для была конкретизирована для каждой стадии релаксации рассматриваемых газов. На основе этих конкретизаций были проведены расчеты коэффициента и квадрата скорости звука для смесей N + N2 и O + O2 в широком диапазоне температур и давлений (при этом использовались соотношения (5)). Результаты соответствующих исследований приведены в Приложении 3 на рис. 26-51.

Расчеты показали, что значение коэффициента существенно зависит от того, описываются внутренние степени свободы классически или считаются квантованными. На 1-й стадии релаксации, когда наряду с поступательными учитываются вращательные степени свободы, описываемые классически, значения зависят лишь от состава газа и соответствуют известному выражению = 1 + 2n /(3n1 + 5n2 ).

На стадиях 2Ц5, когда подключается учет квантованных степеней свободы молекул (колебательных переходов, распада и рекомбинации), значения зависят от плотности и температуры газа, а в неравновесных условиях (на стадиях 2Ц4) еще и от температуры 1-го колебательного уровня T1.

При сравнительно больших фиксированных значениях температуры T1, коэффициент монотонно возрастает с увеличением температуры T.

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

Поведение коэффициента в равновесных условиях при разных начальных плотностях газа проиллюстрировано на рис. 4.

1, 1, 1 5 4 1, 1, 1, 1, 1, 1 T 5 0 0 2 0 0 0 3 5 0 0 5 0 0 0 6 5 0 0 8 0 0 0 9 5 0 0 1 1 0 0 0 1 2 5 0 0 1 4 0 0 0 1 5 5 0 0 1 7 0 0 0 1 8 5 0 Рис. 4. Температурная зависимость в смеси N + Nв равновесных условиях при разной начальной плотности газа.

Кривые 1Ч5 соответствуют плотностям: nL ; 10 nL ; 100 nL ; 0,1 nL ; 0,01 nL.

На стадии 5 (в равновесных условиях) зависимость коэффициента от температуры T не является монотонной. При температурах T < 500K его значение близко к 1,4. Повышение температуры приводит к возбуждению колебательных степеней свободы и уменьшению . При дальнейшем повышении температуры подключаются процессы диссоциации, и уменьшается доля внутренней энергии, приходящаяся на одну частицу. Коэффициент начинает расти, приближаясь к 1,6(6) при температуре, соответствующей полной диссоциации газа. Расчеты показали, что при увеличении степени разреженности газа (уменьшении начального числа молекул в единице объема) убывание и возрастание начинаются при более низких температурах, все кривые как бы сдвигаются влево. Поведение коэффициента в равновесных условиях иллюстрирует рис. Глава 4 посвящена изучению ударных волн, возникающих в высокотемпературном двухатомном газе. Обычно ударная волна рассматривается как узкий (по сравнению с размерами обтекаемого тела) переходный слой от одного состояния равновесия к другому. В силу динамических законов сохранения, параметры на ударной волне удовлетворяют некоторым соотношениям. Их называют условиями динамической совместности. Обычно используются законы сохранения импульса, полной энергии и массовой плотности.

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

Наше рассмотрение ограничивается прямыми скачками уплотнения, возникающими в равновесных газовых смесях. Так как конкретные расчеты проводятся для прямых скачков уплотнения, соответствующие условия совместности записываются в симметричном виде %% p + e p + e ++ +-- += +, 22 (11) ++ ( ) ( - ) (12) = , = 1,..., , + + p = + p, (13) + + + - - где знаками л- и л+ обозначены значения параметров до скачка и после (на границе, завершающей рассматриваемую зону за фронтом волны), параметры соответствуют суммарным значениям аддитивных инвариантов столкновений, учитываемых в масштабах выделенной зоны. До ударной волны параметры считаются известными, а на границе каждой из зон они выражаются через квазистационарные функции распределения (4), сформировавшиеся в соответствующих зонах. Поэтому в левые части условий совместности (11)Ц(13) + + + входят неизвестные интенсивные параметры 0,1,...,, и эти условия можно рассматривать как системы алгебраических уравнений относительно + + + неизвестных ,1,..., и +. На границе каждой релаксационной зоны число 0 уравнений (11) - (13) равно числу неизвестных.

Уравнения (11)Ц(13) являются сложными трансцендентными уравнениями и могут быть записаны в виде + + F +,0,..., = 0, = 0, +1. (14) () Для решения подобных уравнений обычно используют метод Ньютона. Для его применения необходимо найти область значений искомых параметров, в + + которой якобиан перехода D F0, F1,..., F+1 / D +,0,..., отличен от нуля. В () ( ) работе показано, что при любых значениях этих параметров якобиан системы (14) отличен от нуля и обращается в ноль только тогда, когда скорость газа на ударной волне равна скорости звука. Из газовой динамики известно, что это невозможно. Тем самым доказана применимость метода Ньютона для решения систем (11) - (13).

В высокотемпературных двухатомных газах с процессами диссоциации и рекомбинации в ударных волнах можно выделить пять зон, соответствующих пяти стадиям релаксации, которые рассмотрены в первой главе. В настоящей главе рассмотрены релаксационные зоны за фронтом прямых скачков уплотнения. На границах каждой из таких зон выписаны условия совместности (11)Ц(13), приведены схемы численной реализации метода Ньютона для + + + определения неизвестных интенсивных параметров ,1,..., и скорости + 0 и формулы для расчета значений всех наиболее важных газодинамических параметров.

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

Расчеты показали, что скорость газа + на границе любой из зон за скачком меньше скорости до скачка (как и должно быть, согласно теореме Цемплена); температура T+ > T- и возрастает с ростом T- ; скорость звука a+ > a- и возрастает с ростом a- ; число Маха М+ < M- и убывает с ростом М.

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

Зависимость скорости звука a+ на границах разных релаксационных зон от скорости звука a в равновесном набегающем потоке показана на рис. 5.

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

T1+ /T+ a+ 28262422201816141210a_ 455 765 981 1157 1297 1405 1554 1788 20Рис. 5. Зависимость скорости звука a+ от a в диссоциирующем азоте на границах разных релаксационных зон.

Кривые 2 - 5 соответствуют границам зон 2 - 5.

T1+ /T+ 1,1 4 0,0,0,0,T_ 500 1500 2500 3500 4500 5500 6500 7500 8500 95Рис. 6. Зависимость отношения температур T1+ / T+ от T- в смеси N2 + N на стадии 4 (колебательно-химической релаксации).

Кривая 1 соответствует значению = 1000 м/с;

2 - = 1500 м/с; 3 - = 2000 м/с; 4 - = 2500 м/с; 5 - = 3000 м/с.

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

Публикации автора по теме диссертации:

Статьи в рецензируемых журналах и изданиях:

1. Жаркова О.В. Газодинамика диссоциирующего двухатомного газа на разных стадиях релаксации // Вестн. С.-Петерб. ун-та. Сер. 1, 2005. Вып. 4.. 87С. 8791.

2. Жаркова О.В., Рыдалевская М.А. Моделирование структуры ударных волн в диссоциирующем двухатомном газе // Вестн. С.-Петерб. ун-та. Сер. 1, 2007. Вып. 4. С. 7-12.

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

3. Жаркова О.В., Рыдалевская М.А. Состояние диссоциирующего азота за прямым скачком уплотнения // Аэродинамика / Под ред. Р.Н.Мирошина. СПб., 2002. С. 70-78.

4. Жаркова О.В., Рыдалевская М.А. Релаксация двухатомного газа с диссоциацией и рекомбинацией // Аэродинамика / Под ред. Р.Н.Мирошина.

СПб., 2003. С. 94-113.

5. Жаркова О.В., Рыдалевская М.А. Прямые скачки уплотнения в диссоциирующем газе // Тезисы докладов Международной научной конференции по механике Третьи Поляховские Чтения. СПб., 2003. С. 138.

6. Жаркова О.В., Рыдалевская М.А. Интегралы движения и скорость звука изоэнтропийных течений диссоциирующего газа // Тезисы докладов Международной научной конференции по механике Четвертые Поляховские Чтения. СПб., 2006. С. 137.

7. Жаркова О.В., Рыдалевская М.А. Интегралы движения и скорость звука изоэнтропийных течений диссоциирующего газа // Избр. труды Четвертых Поляховских Чтений. СПб., 2006. С. 347-357.

8. Панневиц О.В. (Жаркова О.В.) Пространственно однородная релаксация в изолированном объеме диссоциирующего двухатомного газа // Тезисы докладов Всероссийского семинара по аэрогидродинамике. СПб., 2008. С. 47.

9. Панневиц О.В., Рыдалевская М.А. Разные стадии пространственно однородной релаксации диссоциирующего двухатомного газа // Тезисы докладов Международной научной конференции по механике Пятые Поляховские Чтения. СПб., 2009. С. 140.

Список цитируемой литературы:

1. Рыдалевская М.А. Статистические и кинетические модели в физикохимической газодинамике. СПб., 2003. 248с.

2. Осипов А.И., Уваров А.В. Кинетические и газодинамические процессы в неравновесной молекулярной физике // Успехи физ. наук. 1992. Т. 162, №11. С.

1-42.

3. Нагнибеда Е.А., Кустова Е.В. Кинетическая теория процессов переноса и релаксации в потоках неравновесных реагирующих газов. СПб., 2003. 272с.

4. Rydalevskaya M.A. Kinetic foundation of nonextensive gas dynamics // Proceeding of 24-th Internation symposium on Parefied Gas Dinamics / Ed/ M/ Capitelli Melville, NY, 2005. P. 1073-1079.

5. Treanor C.E., Rich J.W., Rehm R.G. Vibrational relaxation of anharmonic oscillators with exchange-dominated collisions // J. Chem. Phys. 1968. Vol. 48. N 4.

P. 1798-1807.

6. Рыдалевская М.А. Газодинамические скачки в релаксирующем газе // Прикл.

Механика и техн. Физика. 1995. Т. 36. № 3. С. 92-97.

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