Условные обозначения
Введение
1. Колебательная энергия
2. Заселенность колебательных уровней
2.1 Триноровское распределение
2.2 Неравновесное больцмановское распределение
2.3 Равновесное больцмановское распределение
2.4 Составное распределение
3. Поуровневый коэффициент скорости диссоциации
3.1 Равновесный коэффициент скорости диссоциации
3.2 Поуровневый неравновесный фактор
3.3 Поуровневый коэффициент скорости диссоциации
4. Осредненный коэффициент скорости диссоциации
4.1 Осредненный коэффициент скорости диссоциации, вычисленные с использованием неравновесных распределений
4.2 Модель Парка
4.3 Модифицированная формула для осредненного коэффициента скорости диссоциации
5. Результаты
Заключение
Список литературы
Для решения целого спектра современных научно-технических задач, таких как исследование течений в высокоэнтальпийных установках, процессов в активной среде молекулярных лазеров, очистка загрязненной атмосферы, необходимо изучение многокомпонентных смесей реагирующих газов в условиях далеких от равновесия. Отличительной особенностью моделирования таких течений является то, что система уравнений для макропараметров потока включает в себя уравнения релаксации наряду с уравнениями сохранения. Уравнения релаксации содержат коэффициенты скорости физико-химических процессов. Актуальность данной проблемы заключается в том, что коэффициенты скорости этих процессов необходимы для вычисления параметров потока с учетом колебательной и химической кинетики, а также для решения уравнений неравновесной аэродинамики.
Большое количество работ посвящено изучению взаимного влияния распределений молекул по скоростям и уровням внутренней энергии и химических реакций. Исследования И.Пригожина [13] были одними из первых в этом вопросе, за этими исследованиями вышли работы Р.Презента[12] и Дж. Людвига и М.Хейля [1]. В работе [2] было рассмотрено влияние неравновесных распределений на коэффициенты скорости диссоциации, а в работе[17] на скорости обменных реакций в воздухе.
Значительный прогресс, который был достигнут в последние годы в исследовании детальной поуровневой кинетики неравновесных процессов в потоках реагирующих газов, привел к появлению более точных моделей для скорости неравновесных реакций. Основываясь на заселенностях колебательных уровней, полученных из решения уравнений поуровневой кинетики, в работах [6,7,9] были найдены коэффициенты скорости реакции диссоциации. Они значительно отличаются как от равновесных коэффициентов скорости диссоциации, описываемых законом Аррениуса, так и от многотемпературных коэффициентов. В работе Е.В.Кустовой, Е.А.Нагнибедаи соавторов[8] был рассмотрен коэффициент скорости диссоциации в расширяющихся потоках газа с учетом колебательно-химической кинетики. Сравнение результатов траекторных расчетов с моделью Тринора-Маррона для скоростей диссоциации обсуждается в [5].
Целью данной работы было нахождение коэффициентов скоростей диссоциации на основе неравновесных распределений и на основе моделей Парка, Аррениуса и модели [14] в широком спектре температур, а также последующий анализ и сравнение полученных результатов. Для этого в работе приводятся результаты вычислений распределения Больцмана (равновесного и неравновесного), распределения Тринора и составного распределения. Отдельно рассматриваются коэффициенты скорости диссоциации, вычисленные на основе этих распределений, с целью определения влияния на них выбора параметра неравновесности U. Также рассмотрены модели Парка и Аррениуса для того, чтобы определить влияние выбора различных моделей реакций на коэффициент скорости диссоциации.
В работе рассматривается пятикомпонентная воздушная смесьN_2/O_2/ NO/N/O . Принимаются во внимание происходящие в смеси реакции диссоциации:
N_2+M⟷N+N+M (1)
O_2+M⟷O+O+M (2)
NO+M⟷N+O+M (3)
(M=N_2,O_2,NO,N,O-партнер по столкновению)
и обменные реакции:
N_2+O⟷NO+N (5)
O_2+N⟷NO+O (6)
Осредненный коэффициент скорости диссоциации (1-3) определяется следующим соотношением[1]:
k_diss^((M))=1/n_mol ∑_i^L▒〖n_i k_(i,diss)^((M)) 〗 (7)
M =N_2,O_2, NO
где k_(i,diss)^((M)) – поуровневые коэффициенты скорости диссоциации,n_i- заселенности i – го уровня,-количество колебательных уровней молекулы,n_mol-общее число молекул в единице объема.
В данной работе рассматриваются коэффициенты скорости диссоциацииk_diss^((M)), вычисленные на основе неравновесного распределения Больцмана, распределения Тринора и составного распределения. Все вышеперечисленные распределения для заселенностей колебательных уровней зависят от колебательной энергии молекулы и от температуры газа.
В работе приведены результаты расчетов поуровневых факторов неравновесности и коэффициентов скорости диссоциации для молекул N_2,O_2, NO. При этом поуровневые коэффициенты скорости усреднялись по разным колебательным распределениям: триноровскому, больцмановскому и составному. Проведен анализ зависимости поуровневого коэффициента скорости диссоциации от выбора параметра U. Выполнены расчеты равновесного коэффициента скорости диссоциации при использовании различных данных для констант в формуле Аррениуса, рекомендованных в современной литературе. Рассмотрена интуитивная модель Парка длядвухтемпературного коэффициента скорости диссоциации. Проведены расчеты двухтемпературного коэффициента скорости диссоциации, полученного на основе распределения Тринора с учетом новой модели [14], учитывающей зависимость параметра модели от колебательных уровней. Показано влияние уровневых и двухтемпературных моделей диссоциации и неравновесных распределений на осредненные коэффициенты скорости диссоциации.
Полученные результатыпоказали, чтовыбор параметровA и nв законе Аррениуса и параметра Uв модели Тринора-Маррона оказываетзначительное влияние на результаты вычислений коэффициентов скорости диссоциации. Наиболее значительное различие коэффициентов скорости диссоциации, найденных на основе разных распределений, в рассмотренных в работе условиях обнаружено при более низкой температуре газа, что соответствовало значительному колебательному возбуждению (Т
[1]Нагнибеда Е.А., Кустова Е.В. Кинетическая теория процессов переноса и релаксации в потоках неравновесных реагирующих газов. СПб.: Изд-во СПб ун-та, 2003.
[2] Нагнибеда Е.А., Синицын К.А.,Базылевич С.С. Коэффициенты скорости диссоциации в колебательно-неравновесном газе // Вестник СПБГУ. Сер.1,2006,вып.3
[3]Физико-химические процессы в газовой динамике / под редакцией Черного Г.Г. и Лосева С.А. / М.:Науч. мир. 1995. Т. 1., 2002. Т. 2.
[4] Неравновесные физико-химические процессы в аэродинамике / М.: 1972. 344с
[5] M. Capitelli, F. Esposito, E. Kustova, E. Nagnibeda, Chem. Phys. Lett. 330 (2000) 207
[6] CapitelliM., ArmeniseI., Gorse C. State-to-state approach in the kinetics of air components under re-entry conditions// J.Thermophys. Heat Transfer.1997.Vol.11.N 4.570-578.
[7] Colonna G., Tuttafesta M., Capitelli M. e.a. Influence of dissociation rates on the state-to-state vibrational kinetics in nozzle expansions// Rarefied Gas Dynamics/Ed.R. Brun. 1999. Vol.2 .P. 281-288
[8] E. Kustova, E. Nagnibeda, T. Alexandrova, A. Chikhaoui, Chem. Phys. Lett. 377 (5) (2003) 663.
[9] Kustova E.V., Nagnibeda E.A. The influence of the state-to-state distributions behind shock wave on the dissociation rates//Proc. of the 22d Internat.Symp.on Shock Waves/ Ed. G.J. Ball, R. Hillier, G.T. Roberts. Vol. 1. Univ. of Southhampton, UK.2000.P.783-788
[10] Internal Workshop on Radiation of High Temperature Gases in Atmospheric Entry – Part II// 30 Sep.- 10 Oct. 2004 Porquerolles, France
[11] Ludwig G., Heil M. Boundary layer theory with dissociation an ionization // Adv. in Appl.Mech.1960.Vol.6.P.39-118
[12] Present R.D. Chapman-Enskog method in chemical kinetics// J.Chem.Phys.1960.Vol.48.P.4875-4877
[13]Prigogine I.,Xhrouet E., On the perturbation of Maxwell distribution function by chemical reaction in gases// Physica. 1949.Vol.1.P1079-1085
[14] E. Kustova, E. Nagnibeda, G. Oblapenko , A. Savelev, I. Sharafutdinov, Chemical Physics.464 (2016) 1–13.
[15] Marrone P.V., Treanor C.E. Chemical relaxation with preferential dissociation from excited vibrational levels// Phys. Fluids. 1963.Vol. 6. N 9.P. 1215-1221.
[16] Capitelli M., Ferreira.,Gordiets B.F. e.a. Plasma kinetics in atmospheric gases . Berlin: Springer-Verlag , 2000. 302p.
[17] O.V. KunovaO.V. , NagnibedaE.A. State-to-state description of reacting air flows behind shock waves//Chem. Phys. Lett.441 (2014) 66–76.