Труды МАИ. Выпуск № 95
http://trudymai.ru/
УДК 533.6.011
Условия возникновения и величина эффекта высокоскоростного перехлёста в ударно-сжатой смеси газов
Кузнецов М.М.*, Кулешова Ю.Д.**, Решетникова Ю.Г.***, Смотрова Л.В.****
Московский государственный областной университет, ул. Радио, 10 а, Москва, 105005, Россия *е-mail: kuznets-omn@yandex. ru **е-mail:juliaybogdanova@mail. ru
***e-mail: gau1972@mail.ru ^^е-mail: lilysmotrova@mail.ru
Аннотация
Даны строгая формулировка и строгое обоснование необходимых и достаточных условий высокоскоростной поступательной неравновесности в аналитической бимодальной модели ударной волны для бинарной смеси газов. Этот вид неравновесности наиболее существенно проявляется в так называемом эффекте «перехлёста». В этом эффекте, как показано в численных исследованиях, концентрация пар молекул с большими относительными скоростями по сравнению с их тепловыми внутри фронта ударной волны значительно превосходит её значение в поступательном равновесии за волной. Приведены также аналитические оценки этого эффекта, учитывающие наличие внутренних степеней свободы у молекул тяжелого компонента в рэлеевской смеси газов.
Ключевые слова: кинетический, уравнение, неравновесный, смесь газов, ударная волна.
1. Введение
Изучение процессов, протекающих в ударных волнах, весьма актуально. В большинстве работ, посвященных данной тематике, рассматривается течения однокомпонентных газов [1]. В предыдущих работах авторов [2,3], в которых были сформулированы аналогичные теоремы для однокомпонентного (простого) газа, отмечалось, что распространение их на случай бинарной ударно-сжатой смеси газов сопряжено с преодолением ряда существенных трудностей.
Главной из них является то, что в смесях газов бимодальное распределение по скоростям молекул не может быть безоговорочно применено для расчёта структуры ударной волны, как это делалось в простом газе. Как правило, область применимости классического бимодального распределения Тамма - Мотт-Смита с постоянными параметрами в виде аппроксимационных макроскопических скоростей и кинетических температур групп молекул в сверхзвуковых и дозвуковых «крыльях» этого распределения ограничены малыми значениями концентраций одного из компонентов смеси [4]. Только в этом случае удаётся сохранить все преимущества бимодальной аппроксимации парциальной функции распределения, как и в случае ударных волн в простом газе и получить простое аналитическое решение задачи.
В силу этого получение необходимых и достаточных условий высокоскоростной поступательной неравновесности для произвольных, в общем
случае, значений концентрации компонентов смеси газов становится значительно более сложным.
Тем не менее авторами ранее удалось выделить основной масштабный фактор, от которого зависит главный порядок величины рассматриваемого эффекта [2].
В настоящей работе дано обоснование и строгая формулировка необходимых и достаточных условий высокоскоростной поступательной неравновесности для произвольного состава ударно-сжатых смесей газов с бимодальным распределением молекул. Дана так же численная оценка рассматриваемого эффекта в зависимости от определяющих факторов задачи.
2. Теоремы о необходимых и достаточных условиях эффекта «перехлеста» в эволюции функции распределения пар молекул бинарной смеси газов внутри фронта ударной волны
В бинарных смесях газов с функциями распределения для легкого и тяжелого компонентов можно составить три функции распределения пар молекул по модулю относительной скорости [5]. Этими функциями являются: 0(//) - функция распределения пар внутри легкого компонента, вторая - 0(1к) - функция пар легкий-тежелый компонент, третья - 0(М) - функция пар тяжелый - тяжелый компонент. Усреднение данных функций, умноженных на то или иное сечение химической реакции (зависящей так же от модуля относительной скорости) дает константу скорости этой реакции.
Функции пар молекул получаются на основе интегрирования произведения двух одночастичных функций для компонентов бинарной смеси. Для
формулировки искомых теорем, используется аппроксимация Тамм-Мотт-Смита. Согласно этой аппроксимации каждая из одночастичных функций представляет из себя линейную комбинацию (суперпозицию) «холодного»/«горячего» максвелловского распределения в ударной волне, на входе («холодного»), на выходе («горячего») соответственно.
Аналитические исследования, проведенные к настоящему моменту, показали, что постоянство по толщине ударной волны вспомогательных макроскопических параметров, играющих роль температуры и скорости в «холодном» и «горячем» крыле Тамм - Мотт-Смитовской аппроксимации в общем случае невозможно в рамках данной аппроксимации [4]. Исключением является случай малой концентрации «легкого» (газ Лоренца) или «тяжелого» (газ Рэлея) компонента, когда предположение о постоянстве частичных макроскоростей и концентраций, входящих в «холодные»/ «горячие» части Тамм - Мотт-Смитовской аппроксимации функции распределения каждого компонента остается справедливым (при исследовании структуры ударной волны).
Заметим, что в нашем исследовании эффекта высокоскоростного перехлеста функции распределения пар молекул интересен как раз случай не только малых значений концентраций тяжелого компонента по сравнению с легким, но также и его еще более частный вариант, когда mh намного превосходит Ш/. Именно в этих условиях, как показывают численные исследования, эффект перехлеста оказывается наиболее сильным [5].
Возвращаясь к общему случаю произвольных соотношений концентраций компонентов бинарной смеси и их масс, заметим, что в настоящее время, в
численных аналитических исследованиях структуры ударной волны в бинарных смесях на основе распределения Тамм - Мотт-Смита наиболее часто используются три схемы. Общей для всех трех схем является переменность весовых множителей в виде концентраций «холодного» и «горячего» крыла Тамм - Мотт-Смитовского распределения для каждого компонента смеси. Переменность же остальных двух аппроксимационных параметров: макроскоростей пучков молекул и их кинетических температур различна в этих трех схемах.
В наиболее общем случае переменными по толщине волны являются аппроксимационные кинетические температуры и скорости как в «холодном» так и в «горячем» крыльях одночастичных функций распределения.
Однако, из-за большого числа искомых неизвестных макропараметров (в общем случае, не менее восьми) иногда задачу существенно упрощают, сохраняя переменность минимального числа аппроксимационных макропараметров: концентрации, макроскорости и кинетической температуры только в «горячем» или только в «холодном» крыле Тамм - Мотт-Смитовского распределения.
Рассматриваемые теоремы о необходимых и достаточных условиях эффекта высокоскоростного перехлеста относятся к случаю переменности аппроксимационных макропараметров только в «горячем» крыле бимодальных распределений смесей газов, что и составляет главную трудность при их формулировке. Численная же реализация соответствующего эффекта высокоскоростного перехлеста рассматривается для частного случая малых концентраций тяжелого компонента (рэлеевская смесь), где этот эффект наиболее заметен.
Для удобного сопоставления рассматриваемых теорем в случаях простого газа и бинарной смеси приведем краткую сводку результатов, полученных ранее для простого газа [2].
2.1. Простой (однокомпонентный) газ.
Воспользуемся аппроксимацией Тамма - Мотт-Смита для одночастичной функции распределения ¥(х,с) и функции распределения пар молекул б(х ,
Р (х,с) = хР0(с) + ( 1- х) Р (с) (1)
С=х 2^о + 2х ( 1 -х) + ( 1- х) 2 (2)
где *= х = (3)
где б = —, б = —, в1 = 1 , в0 1 = — , О0, 01; 001 - соответственно
«холодная» (перед волной), «горячая» (за волной) и «перекрестная» моды распределений;
¥1 - «холодное» и «горячее» распределения перед и за волной;
Г т ^
Рг (с) = 1 ч 2лкТ] ,
ехр
т(с - и )2
2 кТ
(4)
т - масса молекулы; Т/, щ - вспомогательные аппроксимационные скорости, температуры и концентрации в холодном (/=0) и горячем крыле Тамма-Мотт-Смитовского распределения (/=1) волной, к - постоянная Больцмана, (с-и) -собственная скорость молекулы. Параметры п - переменны по ширине ударной волны, мг-, Т1 - совпадают с их значениями в состояниях равновесия перед (1=0) и за волной (1=1).
У 2
Распределения О0 и О} являются максвелловскими функциями по относительным скоростям g :
Ог (8) =
Г т\
V кТ J
8 ехр
т8 4кТ
(5)
Перекрестная мода имеет вид [6]:
О01( 8) =
т
У2
2жк (Т0 + Т)
8 \ ехр и
т( ^ - иУ
2к (Т0 + Т)
ехр
т( 8 + и)2 ' (То + Т)
(6)
Макроскопические параметры, входящие в соотношения (1) - (3), связаны законами сохранения потоков массы, импульса и энергии в сечениях I = 0 (перед
волной) и I = 1 (за волной):
Т
у = 1 + то(1 -^о2) Т
^ = ^ = ^о = е(1 + т-1)
П и0
и = и0 - щ = и0 (1-е0)
Здесь т0 =е(1 -е) М02, £ =
к-1 х+1
х - отношение удельных теплоемкостей
при постоянном давлении ср и объеме сг, я = — = ( 1 + £) ( 1 — £) 1 , М0 - число
и
Маха перед волной, М0 = —, а0 - скорость звука перед волной, а0 =
а
хкТп
т
Выражение (3) позволяет сформулировать следующие теоремы о «перехлесте» сверхскоростной поступательной неравновесности в бимодальной
ударной волне.
1
Теорема 1. (Необходимость и достаточность) Для высокоскоростного превышения («перехлеста», О > 1) величины поступательно неравновесной функции распределения пар молекул внутри фронта ударной волны (0<х<хь<1) над соответствующей равновесной величиной за волной необходимо и достаточно, чтобы величина перекрестной моды От удовлетворяла соотношению
2Оо1 > 1 + Оо (7)
Теорема 2. (о максимуме относительного «перехлеста» ( 6тах — 1 ) ) Величина сверхскоростного превышения (О > 1) в бимодальном однокомпонентном газе при выполнении соотношения (4) достигает своего
максимального значения
(д — 1) „ = ( 2 1 — 1 — 00 ф 2 (8)
здесь хь - значение переменной х, при котором внутри волны ( 0 < / < 1 ) «обнуляется» эффект перехлеста как и в состоянии «горячего» равновесия за волной при х=0 [2].
Хь = 2(^01 - 1)(2С01 - 1 - С0) ^
Справедливость утверждений обеих теорем непосредственно следует из выражения (2), рассматриваемого как квадратное уравнение относительно параметра х и анализа его дискриминанта на положительную определенность.
2.2. Бинарная смесь газа.
Для обобщения предыдущих результатов, полученных для простого газа, на случай смеси газов соотношение (2) удобно записать, используя значения переменнойх= Хъ в следующем виде:
(С-1 )= х ( 2 С01 -1-(Го) (Хь-Х) (10)
Видно, что в интервале (0<х<хь<1) условие (7) и эффект перехлеста (6-1 ^ > 0 следует друг из друга.
Рассмотренные выше теоремы о необходимых и достаточных условиях эффекта высокоскоростного перехлеста в простом газе переносятся, практически дословно) на случай бинарной смеси, когда переменными по ширине ударной волны являются все вспомогательные параметры и1, Т1, п1 в «горячем» крыле Тамм-Мотт-Смиттовского распределения. В «холодном» же крыле этого распределения переменной величиной является только концентрация молекул п0. Далее, для перехода к смеси газов, остаётся только произвести надлежащую замену формул.
Так, формулы (1), (2), (3) заменяются, соответственно, на следующие:
Р 00 = х ООр^ + (1 - х(у)) рТ (11)
С(а,ю = х(а)х(Юс^+х(а\1 - - х(а)) Х^С™ + (1 -
х (а) (1-Х(Ю) б(а,Ю (12)
Дг) __п
(г)
^ (г) , (г) (13)
Здесь: у - сорт молекул, у=а,Р (в дальнейшем под сортом а будем понимать легкий газ, причем а=1, под в - тяжелый газ, Р=И.
Функции ^ \ ^ ' ), ^ ' ) совпадают в простом газе с функциями
с{а,Р)
( о ' (1' ( о 1 ' (1 о ' при а = //, (7 (а'Р) =
Аналогично формулу (10) следует заменить на:
{—Х + Хь) (14)
где Х = Х (а)+Х (//), »7 =
Хь =
Заметим, что формула (14) непрерывно переходит в формулу (10), поскольку при переходе от смеси к простому газу будет:
у(а) _ у(/?) г, — 1 г(а,Ю-г(а,/3) _ я (а,/?) _ ^ л(а,/?) _ 1
X = х » Л = "о 1 _С/1 0 - Ь 0 1' Ь0 = Ь0' Ч = 1 .
Как упоминалось ранее, в простом газе необходимые и достаточные условия эффекта «перехлеста» следовали из формулы (10). Подобно этому, аналогичные условия в случае смеси газов будут следовать из формулы (14).
Более того, текст теорем 1 и 2 сохраняется прежним, а формулы (7) и (8) заменяются на следующие:
(е^+> Ф^ + (15)
Ф(аР) - С.; (а'Ю)тах = щ ( 1 - ф^ + С^ - С^ - С^) Й (16)
Заметим однако, что в случае бимодальной смеси газов эффект перехлеста в
ударной волне состоит в выполнении локального неравенства: б(а 'Р) > £ а'Р ,
поскольку величина и. зависит от координаты внутри ударной волны. В простом же газе аналогичное неравенство будет следующим:
(б - 1 ) >0. В общем
случае, нельзя быть уверенным в том, что всегда £ / а'Р > 1 . Однако в рэлеевской смеси газов, когда концентрация легкого компонента п(1) значительно превосходит
концентрацию тяжелого компонента п(ъ) всегда будет £ / а 'Р = 1 [4].
Для рэлеевской смеси не трудно рассчитать также максимум эффекта высокоскоростного перехлеста, аналогично случаю простого газа. Например, для максимального значения функция распределения пар тяжелого компонента смеси
£ (Кк), когда отношение массы молекулы тяжелого компонента тъ к массе молекулы легкого т1 равно двум, можно составить следующую таблицу.
Таблица 1
Максимум "перехлеста" Отах в сверхзвуковом потоке
Оа8 А (А2), без учета колебательных степеней свободы молекулы (А2), с учетом колебательных степеней свободы молекулы (А3), с учетом вращательных и колебательных степеней свободы молекул С8 Н16
У 5/3 7/5 9/7 7/6 22/21
s 1/4 1/6 1/8 1/13 1/43
G 1,86 2.8 11.7 648/392 6.48-1016
В этой таблице величина s задавалась в качестве параметра для молекул газов с различным количеством атомов: одноатомных (А), двухатомных (A2) , трехатомных (A3) , и многоатомных (типа ). Рассмотрены отдельно случаи
s = 1 (у=1,4) и s = 1 |/ = 9j, соответствующие отсутствию или наличию
возбужденных колебательных степеней свободы у двухатомных газов (A2) . Данные значения параметров у и s задавались параметрически. Справедливость такого параметрического задания, т.е. завершенность процессов релаксации внутри фронта ударной волны, показана для вращательных степеней свободы в случаях: (A2) и (A3) в работе [7].
Работа поддержана грантом РФФИ№17-07-00-945А.
Библиографический список
1. Рыжов Ю.А., Никитченко Ю.А., Парамонов И.В. Численное исследование гиперзвукового обтекания острой кромки на основе модели Навье-Стокса-Фурье // Труды МАИ. 2012. № 55. URL: http://trudy.mai.ru/published.php?ID=30027
2. Кузнецов М.М., Матвеев С.В., Молоствин Е.В., Решетникова Ю.Г., Смотрова Л.В. Высокоскоростная поступательная неравновесность смеси газов в аналитической модели ударной волны // Физико-химическая кинетика в газовой динамике. 2016. Т.17. №1. URL: http://chemphys.edu.ru/issuse/2016- 17-1/articles/613/
3. Кузнецов М.М., Кулешова Ю.Д., Смотрова Л.В, Решетникова Ю.Г. О максимуме эффекта высокоскоростной поступательной неравновесности в ударной волне // Вестник МГОУ. Физика-математика. 2016. № 3. С. 84-95.
4. Oberai M.M. A Mott-Smith distribution to describe the structure of a plane shock wave in binary mixture // Phys Fluids.1966. Vol. 9. P.1634-1678.
5. Генич А.П., Куликов С.В., Манелис Г.Б., Черешнев С.Л. Распределение молекулярных скоростей во фронте ударной волны в газовых смесях // Механика жидкости и газа. 1990. № 2. С. 144-150.
6. Куликов С.В., Терновая О.Н., Черешнев С.Л. Специфика поступательной неравновесности во фронте ударной волны в однокомпонентном газе // Химическая физика. 1993. Т. 12. № 3. С. 340-342.
7. Ching Shen, Zhenhua Hu, Wanquan Wu, Xiaoyan Xu. Shock Waves in Gas Mixtures with Internal Energy Relaxation. Rarefied Gas Dynamics // Proceedings of the 17th International Symposium on Rarefied Gas Dynamics July 8-14, 1990, Aachen, Germany. P. 247-254.