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

2019 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 6 (64). Вып. 1

МАТЕМАТИКА. МЕХАНИКА. АСТРОНОМИЯ

МЕХАНИКА

УДК 533.6.011 МБС 82С40

Время релаксации вращательной энергии колебательно возбужденных молекул*

А. И. Бенина, Е. В. Кустова

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

Для цитирования: Бенина А. И., Кустова Е. В. Время релаксации вращательной энергии колебательно возбужденных молекул // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2019. Т. 6(64). Вып. 1. С. 118-130. https://doi.org /10.21638/11701^рЬи01.2019.109

В работе изучается влияние колебательного уровня молекулы на время релаксации вращательной энергии в приближении поуровневой кинетики. Вращательные уровни молекул описываются моделью нежесткого ротатора; взаимодействие молекул описывается моделью переменной мягкой сферы. По этой модели рассчитываются сечения столкновений N2-^ 02-0, N0-0 для разных колебательных и вращательных уровней молекул. Методами кинетической теории неравновесных процессов вводится время релаксации вращательной энергии для каждого колебательного уровня. В широком диапазоне температур проводится численный расчет времен релаксации и сравнение с временем релаксации, полученным по известной формуле Паркера. Анализируется влияние различных многоквантовых вращательных переходов на точность расчета времени вращательной релаксации; показана сходимость решения при увеличении максимально возможного числа квантов, переданных при переходе. Установлено, что колебательное состояние молекулы заметно влияет на время релаксации вращательной энергии в поуровневом приближении; при расчете поуровневых коэффициентов переноса применение формулы Паркера может приводить к заметной ошибке. В од-нотемпературном приближении при умеренных температурах формула Паркера дает удовлетворительное согласие с результатами, полученными при усреднении поуровне-вых времен релаксации с распределением Больцмана по колебательной энергии. Ключевые слова: время релаксации, вращательная энергия, возбужденное колебательное состояние, поуровневое приближение.

* Работа выполнена при поддержке СПбГУ (тема НИР 6.37.206.2016) и РФФИ (грант №18-0100493).

(¡5 Санкт-Петербургский государственный университет, 2019

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

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

Tel < Trot < Tvibr < Treact ~ в. (1)

Здесь Tei, Trot, Tvibr, Treact — соответственно времена релаксации поступательных, вращательных, колебательных степеней свободы и характерное время химических реакций, в — характерное время изменения макропараметров газа. Данное соотношение времен подтверждается экспериментально [2].

При использовании строгих методов кинетической теории время релаксации вращательной энергии при условии (1) должно вводиться для каждого колебательно возбужденного состояния молекулы [1]. Однако экспериментальные данные для таких времен релаксации в настоящее время отсутствуют. На практике при вычислении времени релаксации вращательной энергии обычно пользуются формулой Паркера [3], полученной в предположении, что Trot не зависит от колебательного спектра молекулы. Данное предположение заметно упрощает расчет времени релаксации и, следовательно, коэффициентов переноса в поуровневом приближении.

Количественная оценка влияния учета колебательного состояния молекулы на время релаксации вращательной энергии и коэффициенты переноса в поуровневом приближении обобщенного метода Энскога—Чепмена до настоящего времени не проводилась. Отсутствие подобного анализа не позволяет оценить точность расчета коэффициентов переноса, особенно коэффициента объемной вязкости [4, 5]. Цель настоящей работы — исследовать времена вращательной релаксации для столкновений молекул с атомами с учетом влияния колебательного возбуждения на сечения наиболее вероятных упругих столкновений и оценить пределы применимости упрощенных моделей.

2. Приближение поуровневой кинетики. Введем макропараметры, характеризующие течение неравновесной реагирующей смеси газов в поуровневом приближении [1]: nci(r,t) —заселенность колебательного уровня i молекул сорта c (r — координата, t — время), v(r, t) — макроскопическая скорость газа, T(r, t) — температура (вместо температуры может рассматриваться полная энергия смеси в расчете на единицу массы U(r,t)). Массовая плотность смеси газов р = ^ci mcnci, mc — масса частицы сорта c.

Система уравнений для макропараметров пС^(т,г), ,о(г,г), и(г,г) записывается в виде [1]

dnc ~dt dv „ „

dU ^ „ ^ p— + V • q + P : Vv

dt

+ nciV ■ v + V ■ ( nci Vci

0.

Rc

, L, i = 0,

I Lc

(2)

(3)

(4)

Здесь Ь — число химических сортов в смеси; ЬС — число возбужденных колебательных уровней молекул сорта с; — скорость изменения заселенностей колебательных уровней за счет медленных процессов (обменов колебательной энергией и химических реакций), Р — тензор напряжений, — скорость диффузии молекул сорта с, находящихся на колебательном уровне г, д — поток полной энергии. Выражения для потоковых членов в поуровневом приближении приведены в [1].

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

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

а.ег(г) Ег(г) - ег°

dt

Trot

(5)

где Ег(г) —энергия вращательных степеней свободы, отнесенная к единице объема, Е° —равновесное значение этой энергии. Тогда время релаксации тгс^ можно записать в виде

E0

Trot

dET(t) dt

(6)

Er=0

Широко используемая формула Паркера [3] была получена при тех же предположениях, что и уравнение (5). В соответствии с этой формулой число столкновений ССЫ, необходимое для установления равновесия по вращательным степеням свободы молекул сорта с, равно

С г

rotc

Fc(T )'

где функция ЕС(Т) имеет вид

Fc(T ) = 1 +

т3/2

кт \-1/2 +

( кТ\

1

( кТ\

Т+2HtJ +'3/4tJ-

2

2

Здесь к —постоянная Больцмана; Е* —глубина потенциальной ямы в потенциале взаимодействия; величины ^ определяются из соответствия теоретических и экспериментальных значений. Формула для времени вращательной релаксации молекул сорта с тГы выглядит следующим образом [1]:

4 ртгЫ

СГ = (7)

п п

•с

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

Определим оператор осреднения в поуровневом приближении [1]:

1/2 сг

,rp \ 1/2 „сг„аь r

) £ z^ Iexp(-^2 - ef - (8)

jlj' l' ci dk

Здесь j, l —вращательные уровни молекул сортов c и d до столкновения, j', l' — после столкновения, — статистическая сумма вращательных степеней свободы молекул сорта c на колебательном уровне i; scl = 2j + 1 — статистический вес, который характеризует вырождение состояния молекулы с вращательной энергией £сг; mcd = mcmd/(mc + md) —приведенная масса частиц сортов c и d; Fcij —некоторая функция, зависящая от скорости, химического сорта c и уровней внутренней энергии i, j; y = (mcd/(2kT)1/2g), j = £cl/(kT) — безразмерные относительная скорость

и вращательная энергия; ajJ—сечение наиболее вероятных упругих столкновений и столкновений с переходами вращательной энергии; d2il — телесный угол, в который попадает вектор относительной скорости после столкновения. Выражение для вращательной энергии j имеет вид

£d

-fc = B?j(j + l)-D?j2(j + l)2, (9)

где h — постоянная Планка, c — скорость света, а спектроскопические постоянные Bc, Dc приводятся в [6]. Видим, что в общем случае вращательная энергия зависит от колебательного состояния молекулы i. Для модели жесткого ротатора спектроскопические постоянные равны своим равновесным значениям, Bc = Bc, Dc = Dc, и вращательная энергия не зависит от колебательного состояния.

Предположим, что вероятность столкновений, при которых происходит одновременное изменение вращательной и колебательной энергии, мала. Это позволяет определить время вращательной релаксации rc°dk частиц химического сорта c на колебательном уровне i при столкновении с частицами сорта d, находящихся на k-м колебательном уровне [1]:

i = -^-(^Т (10)

1 cidk mccrot , ci

Здесь n = c i nci — числовая плотность смеси газа; безразмерный дефект резонанса вращательной энергии при столкновении ^EcStk определяется следующим образом [1]:

Ascotk = + K£dt,

0.7 -

3 0.6 -

О

"С 0.5 -

ий 0.4 -

0.3 -

0.2 -

0.1 -

0 -0

-N0, Т= 1300 К N0, Т= 5300 К

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

Д£ТО1 = V-Д£Го1

к Г ' кТ '

В поуровневом приближении вводится удельная вращательная теплоемкость сГо1,сг молекул сорта с, находящихся на колебательном уровне г [1]. Для модели жесткого ротатора вращательная энергия не зависит от колебательного состояния, поэтому сГоь,сг = сГ01;,с для любого г. Обычно теплоемкость вращательных степеней свободы принимается постоянной, сГо11,с = к/тс. Для более реалистичной модели нежесткого ротатора удельная теплоемкость вращательных степеней свободы при постоянном объеме V зависит от температуры и колебательного уровня:

_ I дЕТо^ст,

V дТ / V

удельная вращательная энергия ЕГо1сг молекул сорта с на колебательном уровне г вводится выражением

1 ес М I £

Е^сг —У, ехР ( ) ' & = Е >Т 'Ч Т-

с^сг

На рис. 1 приведены отношения сГо1сг/сГо1с как функции колебательного уровня для компонентов воздуха при фиксированных температурах. Хорошо заметно, что для нижних уровней отношение близко к единице. Однако, чем выше температура газа, тем быстрее теплоемкость уменьшается с ростом колебательного уровня. Для проверки расчета была также вычислена теплоемкость сГо1,с, полученная с помощью осреднения теплоемкости сГо1;,сг с больцмановским распределением заселенностей псг по колебательной энергии; результат в точности совпал с сГо1,с, рассчитанной методами термодинамики. В случае отклонения распределения псг от больцмановского возможно заметное влияние модели ротатора на коэффициенты переноса; этот факт отмечается ив [4].

Если расписать подробно время вращательной релаксации (10) с помощью оператора осреднения (8), подставив в него выражения для изменения вращательной

£ак

£ак

энергии (9), то получим выражение, в котором заданы все величины, кроме сечения столкновения • Вид сечения столкновения зависит от модели взаимодействия

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

4. Статистическая модель неупругого сечения столкновения. Реалистичные модели для сечений упругих и неупругих столкновений необходимы, в первую очередь, при моделировании течений разреженного газа методом прямого статистического моделирования (ПСМ). Для упругих молекулярных столкновений Бердом была предложена модель переменной твердой сферы (VHS) [7], основанная на обратном степенном потенциале. Модель твердых сфер проста, но не дает корректных значений коэффициентов вязкости и диффузии для такого потенциала. Чтобы снять это ограничение, Коура и Мацумото [8] ввели модель переменной мягкой сферы (VSS). Рассмотрим подробнее модель VSS.

В методах ПСМ для описания молекулярных столкновений вводят полные сечения а^ (ij, е) неупругих столкновений:

= [С(E) - Zj,i (E)]aei (е). (11)

Здесь E — полная энергия столкновения, е —относительная энергия до столкновения, ае1(е) —полное упругое сечение. Функциональная форма Zj,i(E) определяется следующим образом [9]:

Zji i(E) = C (E)csfsf(E - eci - efk)(1-u), (12)

Z(E) аппроксимируется с хорошей точностью выражением

Z(E) = (4 - u)AEr(E)/E. (13)

Функция C(E) в равенстве (12) специально определяется для вращательных неупругих столкновений с использованием классической функции энергии Паркера [3]. Для каждого потока молекул C(E) принимает вид

C (E) = 2AEr(E )/Gr (E). (14)

Функция вращательной энергии, приобретенной за счет столкновений молекул, определяется следующим образом [5]:

AER(e) = [е/2 + (2n/3)(eE*c)1/2 + (п2/4 + 2)EC]/ZT- (15)

Функция Gr(E), определяемая суммированием по дискретному вращательному спектру, может быть аппроксимирована интегрированием по непрерывному вращательному спектру, поскольку значение характеристической вращательной температуры 6R мало. Для упругого сечения VSS, классическая форма функции GR(E) принимает вид

Gr(E) = [c/[2(2 - ш)(3 - ш)(4 - u)(k0R f]}E-. (16)

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

aei (е) = се-и, (17)

где значения постоянных c и ш приведены в [10] для различных газов.

1

Рис. 2. Сечения столкновений как функции вращательного уровня ] при фиксированной относительной энергии и различных колебательных уровнях.

е, К

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

Параметры для столкновений N2— N О2 — О и N0-0 приведены в таблице.

Использованные константы

Столкновение ш с, м" лоо Е*, Дж вн, К

N2^ 0.275 240.3 • КГ™ 20.3898 1.15127 • Ю--1 2.86256

о2-о 0.199 114.7- КГ™ 20.8478 1.22331 • Ю--1 2.05706

N0-0 0.217 135.4 • КГ"и 20.3898 1.17846 • 10~-х 2.39326

На рис. 2 приведена зависимость сечения столкновения молекулы азота с атомом азота, рассчитанного по модели VSS, от вращательного уровня при фиксированной относительной энергии до столкновения и различных колебательных уровнях. Видно, что с ростом вращательного уровня сечение уменьшается, причем на нулевом колебательном уровне примерно на 13%, а на 40-м — на 3%. Также можно отметить, что сечение уменьшается с ростом колебательного уровня. На рис. 3 представлено изменение сечения столкновения в зависимости от относительной энергии до столкновения (выраженной в Кельвинах) при фиксированном вращательном уровне для разных колебательных уровней. С ростом относительной энергии сечение уменьшается.

5. Расчет времени релаксации вращательной энергии. Для упрощения расчетов будем считать, что партнер по столкновению — атом. Тогда формула (10) сводится к более простой:

сг тссго%,сг

Интегральный оператор осреднения (8) и сечения столкновений в данном случае не зависят от уровней внутренней энергии партнера по столкновению к, I. В результате

1 4 кп ( кТ \1/2

тссго^сг

г™* тсСГо1,сг I 27Гтсс1.)

Х Е 111 - ^) V ехР("Г - ) (7, X, е) 8ШХ ,1\,1 <*7. (19)

Здесь использовано соотношение cPí! = sinxdxde, х — угол рассеяния, е — азимутальный угол рассеяния. Для модели VSS сечение не зависит от е, и интегрирование по этому углу дает 2п. Кроме того, сечение, определенное формулой (11), уже осред-нено по углам рассеяния, то есть

= 2п j acd, ij Ь,х) sin Xdx-Это позволяет переписать время релаксации в виде

1 4 кп f кТ \

Tlf rriccI0t d V 2imicd I

1/2

EÄt (£?-£Г)2^£Г / 73e-<2iaE)-QAE)Wi(eyb- (20)

зз'

Здесь учтено выражение (11) и из-под интеграла вынесены все величины, не зависящие от скорости 7.

Рассмотрим отдельно интеграл, входящий в это выражение. Подставим в него равенства (13), (12) и (17). Тогда интеграл представляется в виде

3__

Y e 1

(4 - _ c(E)cs?(E -

E j j

се-ш dq. (21)

J7min

Учитываем, что полная энергия столкновения имеет вид

E = е + Sci + SC,

и принимаем во внимание, что е = (mCdg2)/2, y = (mcd/(2kT))1/2g. Подставляем эти соотношения в выражения для AEr(E), C(E) и Gr(E) и в результате получаем окончательное выражение для времени релаксации вращательной энергии для модели VSS.

Пусть при столкновении происходят многоквантовые переходы вращательной энергии j' = j ± Aj, где Aj = 2, 4, 6,... для гомоядерных молекул и Aj = 1, 2, 3,... для гетероядерных. Тогда сумма в формуле (20) имеет вид

NroM-1 j+Aj<Nrot,i-1

£= £ £ , (22)

jj' j=0 j-Aj>0

где Nrot,i — количество вращательных уровней у молекулы на колебательном уровне i.

Закон сохранения энергии записывается следующим образом:

Y2 + sci + SC = y2 + Ep + SC,

где Sp — безразмерная колебательная энергия молекул сорта c на уровне i. В быстром процессе (поступательно-вращательная релаксация) колебательная энергия не изменяется. Перенесем SjC'i в левую часть

Y2 + Sci -Sp = Y2.

х

х

Y

оо

-Д >2

-Д ¡=6

-д >10

- Рагкег

0 2000 4000 6000 8000 10000

Т, К

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

Правая часть уравнения неотрицательна. Следовательно левая тоже должна быть неотрицательна. Если < ] ,то Е' — £' > 0, т. е. левая часть всегда неотрицательна. Следовательно, в этом случае относительная скорость может быть любой, и тогда 7тт = 0. Если ]' > ], то Е' — £' < 0 и тогда появляется условие на 7:

В этом случае 7тт = у (£р — £/).

Так как интеграл (21) аналитически взять не удается, он рассчитывается численно методом квадратур Гаусса—Кронрода по 15 точкам. Следует также отметить, что расчет сумм (22) при больших Дj требует больших вычислительных затрат, поскольку число вращательных уровней в молекулах велико. Так, для N2 имеем N■0^0 = 197, для 02 — Жго-^о = 168, для N0 — ^^о = 176. С ростом колебательного уровня число возможных вращательных состояний уменьшается.

6. Результаты. Ниже представлены результаты расчета времен вращательной релаксации, зависящих от колебательного состояния, в компонентах воздуха N2, О2, N0.

Сначала оценим влияние числа возможных многоквантовых обменов на время релаксации. На рис. 4 представлены времена релаксации, рассчитанные численно по формуле (20) для столкновения молекулы азота, находящейся на 20-м колебательном уровне, с атомом азота. В качестве максимальных значений Дj брались 2, 6 и 10. Видно, что чем больше вращательных переходов мы учитываем, тем ближе решение к решению Паркера. При дальнейшем увеличении Дj точность расчета повышается незначительно, но при этом решение требует больших вычислительных ресурсов. Дальнейший анализ проводится при Дj = 10.

На рис. 5 сравниваются времена релаксации вращательной энергии для столкновений N2—N и О2—О, рассчитанные для разных колебательных уровней, с временем релаксации, полученным по формуле Паркера. Можно видеть, что на разных колебательных уровнях время релаксации ведет себя по-разному. Для верхних уровней зависимость от температуры может быть немонотонной. Такое поведение

10

10

10

10

4000 6000

Т, К

4000 6000

Т, К

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

10

10

10

10

10

0

2000

8000

10000

0

2000

объясняется прямой зависимостью времени релаксации от удельной вращательной теплоемкости (см. рис. 1). Для нижних уровней теплоемкость практически постоянна, но резко убывает с ростом г и Т. Убывание времени вращательной релаксации для верхних колебательных состояний обусловлено аналогичным поведением сго^Сг. Для столкновения N0—0 график выглядит аналогично.

Интересно сравнить время релаксации вращательной энергии, рассчитанное по формуле Паркера, с т(Гак, которое получено осреднением по колебательной энергии поуровневого времени релаксации (18). Для этого используем распределение Больц-мана с температурой газа Т:

Сго^с _ 1 \ ' ^ > у ! (^3)

Здесь ИрЬг — статистическая сумма колебательных степеней свободы молекулы сорта с, ЕС —безразмерная колебательная энергия. На рис. 6 и 7 можно видеть, что при

4000 6000

Т, К

УЗЗ N0-0 Рагкег N ^

-Рагкег N0-0 8000 10000

Рис энергии, ням, как

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

5

4.5

4

3.5

1 3

я

-4: 2.5

2

1.5

0.5

0

4000 6000

Т, К

Рис. 7. Отношение усредненного времени релаксации к времени релаксации по Паркеру как функция температуры.

10

10

0

2000

0

150

I

100

20 30 40 7=9300 К

1200 1000 800 600

з >

>- 400

200 0

50

0

Т < 10000 К время релаксации, усредненное с распределением Больцмана по колебательным уровням, достаточно близко ко времени релаксации, рассчитанному по формуле Паркера. Из рис. 7 видно, что отношение времен релаксации растет почти линейно с температурой; расчеты в более широком диапазоне температур показали, что при Т = 30000 К это отношение достигает 17.6 для столкновения N2 — N 14.1 для 02—О и 7.9 для N0—0. Поскольку формула Паркера получена для однотем-пературного приближения и умеренных температур, данный результат можно рассматривать как проверку достоверности модели при условии слабого отклонения от равновесия по внутренним степеням свободы. При сильных отклонениях распределений от больцмановских усредненное время релаксации может заметно отличаться от времени, рассчитанного по формуле Паркера.

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

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

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

Литература

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

2. Ступоченко Е. В., Лосев С. А., Осипов А. И. Релаксационные процессы в ударных волнах. М.: Наука, 1965. 484 с.

3. Parker J. G. Rotational and vibrational relaxation in diatomic gases // Phys. Fluids. 1959. Vol. 2. P. 449-461.

4. Kustova E., Mekhonoshina M., Oblapenko G. On the applicability of simplified state-to-state models of transport coefficients // Chem. Phys. Lett. 2017. Vol. 686. P. 161-166.

5. Kremer G. M., Kunova O., Kustova E., Oblapenko G. The influence of vibrational state-resolved transport coefficients on the wave propagation in diatomic gases // Physica A. 2018. Vol. 490. P. 92-113.

6. Huber K. P., Herzberg G. Molecular Spectra and Molecular Structure. IV. Constants of Diatomic Molecules. New York: Van Nostrand, 1979.

7. Bird G. A. Monte-Carlo simulation in an engineering context // Prog. Astronaut. Aeronaut. 1981. Vol. 74. P. 239-255.

8. Koura K., Matsumoto H. Variable soft sphere molecular model for inverse-power-law or Lennard-Jones potential // Phys. Fluids A. 1991. Vol. 3. P. 2459-2465.

9. Koura K. Statistical inelastic cross-section model for the Monte Carlo simulation of molecules with discrete internal energy // Phys. Fluids A. 1992. Vol. 4. P. 1782-1788.

10. Koura K., Matsumoto H. Variable soft sphere molecular model for air species // Phys. Fluids A. 1992. Vol. 4. P. 1083-1085.

Статья поступила в редакцию 10 июля 2018 г.;

после доработки 15 августа 2018 г.; рекомендована в печать 27 сентября 2018 г.

Контактная информация:

Бечина Анна Ильинична — студент; anya_bechina@mail.ru

Кустова Елена Владимировна — д-р физ.-мат. наук, проф.; e.kustova@spbu.ru

Rotational energy relaxation time of vibrationally excited molecules

A. I. Bechina, E. V. Kustova

St. Petersburg State University, Universitetskaya nab., 7-9, St. Petersburg, 199034, Russian Federation

For citation: Bechina A. I., Kustova E. V. Rotational energy relaxation time of vibrationally excited molecules. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2019, vol. 6(64), issue 1, pp. 118-130. https://doi.org/10.21638/11701/spbu01.2019.109 (In Russian)

The influence of the molecular vibrational state on the rotational energy relaxation time is studied in the state-to-state approach. The rotational levels of the molecules are described by the nonrigid rotator model; their interaction is described by the variable soft sphere model. According to this model, collision cross sections for N2-N, O2-O, NO-O interactions are calculated for different vibrational and rotational levels of molecules. On the basis

of the kinetic theory of nonequilibrium processes, the rotational relaxation time is introduced for each vibrational level. Relaxation times are calculated numerically in the wide temperature range, and a comparison with the relaxation time obtained by the well-known Parker formula is carried out. The influence of various multi-quantum rotational transitions on the accuracy of the calculation of the rotational relaxation time is analyzed; the convergence of the solution is shown with increasing the maximum possible number of quanta transferred during the transition. It is shown that the vibrational state of the molecule significantly affects the rotational energy relaxation time in the state-to-state approach; when calculating the state-specific transport coefficients, the use of the Parker formula can lead to a noticeable error. In the one-temperature approximation for moderate temperatures, the Parker formula provides a satisfactory agreement with the results obtained by averaging the state-resolved relaxation times with the Boltzmann distribution over vibrational energy.

Keywords: relaxation time, rotational energy, vibrationally excited state, state-to-state approach.

References

1. Nagnibeda E. A., Kustova E.V., Kinetic Theory of Relaxation and Transport Processes in Nonequilibrium Reacting Flows (Saint Petersburg University Press, Saint Petersburg, 2003).

2. Stupochenko E. V., Losev S. A., Osipov A. I., Relaxation Processes in Shock Waves (Nauka Publ., Moscow, 1965).

3. Parker J.G., "Rotational and vibrational relaxation in diatomic gases", Phys. Fluids 2, 449—461 (1959).

4. Kustova E., Mekhonoshina M., Oblapenko G., "On the applicability of simplified state-to-state models of transport coefficients", Chem. Phys. Lett. 686, 161—166 (2017).

5. Kremer G. M., Kunova O., Kustova E., Oblapenko G., "The influence of vibrational state-resolved transport coefficients on the wave propagation in diatomic gases", Physica A 490, 92—113 (2018).

6. Huber K. P., Herzberg G., Molecular Spectra and Molecular Structure. IV. Constants of Diatomic Molecules (Van Nostrand, New York, 1979).

7. Bird G. A., "Monte-Carlo simulation in an engineering context", Prog. Astronaut. Aeronaut. 74, 239-255 (1981).

8. Koura K., Matsumoto H., "Variable soft sphere molecular model for inverse-power-law or Lennard-Jones potential", Phys. Fluids A 3, 2459-2465 (1991).

9. Koura K., "Statistical inelastic cross-section model for the Monte Carlo simulation of molecules with discrete internal energy", Phys. Fluids A 4, 1782-1788 (1992).

10. Koura K., Matsumoto H., "Variable soft sphere molecular model for air species", Phys. Fluids A 4, 1083-1085 (1992).

Received: July 10, 2018 Revised: August 15, 2018 Accepted: September 27, 2018

Author's information:

Anna I. Bechina — anya_bechina@mail.ru Elena V. Kustova — e.kustova@spbu.ru