Приближение поуровневой кинетики является одним из наиболее точных методов описания сильнонеравновесных течений газов с внутренними степенями свободы и химическими реакциями [1-9]. Для моделирования вязких течений необходим расчет коэффициентов переноса, входящих в выражения для тензора напряжений и теплового потока. Алгоритм расчета коэффициентов переноса в поуров- невом приближении впервые был предложен в работе [6]. Одно из предположений, лежащих в основе алгоритма, состоит в том, что сечения упругих столкновений считаются независящими от колебательного состояния молекулы. Это предположение значительно упрощает расчеты поуровневых коэффициентов переноса, однако в [6] не приводится строгого обоснования пределов его применимости. Более того, известно, что сечение упругого столкновения возбужденных молекул заметно увеличивается [10-13], что может влиять на интегралы столкновений и коэффициенты переноса.
В работе [14] был предложен простой метод для получения приближенных оценок влияния молекулярного диаметра на коэффициент сдвиговой вязкости. Однако для расчета диаметра колебательно возбужденной молекулы использовалась модель [10], существенно переоценивающая размер молекул на верхних колебательных уровнях. В связи с этим расчеты приходилось обрывать при достаточно низких значениях колебательной энергии, что заметно ухудшает точность оценок. Целью настоящей работы является оценка вклада возбужденных состояний в коэффициент сдвиговой вязкости компонентов воздуха на основании более точных моделей потенциалов с учетом всех возможных колебательных состояний. Исследуется отношение коэффициента вязкости к соответствующему коэффициенту в основном колебательном состоянии для различных значений температуры, равновесных и неравновесных распределений.
Проведенный в работе расчет диаметров молекул N2, O2, NO для разных колебательно-вращательных состояний по трем моделям потенциала позволил определить, что простая модель Канга—Кунца дает экспоненциальный рост диаметра молекулы для уровней выше 10, поэтому ее применение целесообразно только при низких температурах, а модели Тица—Хуа и Морзе дают близкие значения диаметров. При этом вкладом вращательного возбуждения в диаметр рассмотренных молекул можно пренебречь.
Анализ результатов расчета отношения поуровневого коэффициента сдвиговой вязкости к коэффициенту вязкости для газа, состоящего из невозбужденных молекул, показал, что во всех рассмотренных случаях, включая сильнонеравновесные, эффект увеличения молекулы практически не влияет на вязкость, в среднем отклонение составляет лишь 1-2%. Некоторый эффект наблюдается только при очень высоких температурах для близких к равновесным колебательных распределений. Однако такие условия на практике реализовать достаточно сложно, поскольку в околорав- новесном высокотемпературном течении доля недиссоциирующих молекул пренебрежимо мала. Таким образом, нам удалось для широкого диапазона условий показать справедливость предположения о том, что при расчете поуровневых коэффициентов переноса зависимость сечения упругого столкновения от колебательного уровня можно не учитывать. Это позволяет обоснованно применять упрощенные алгоритмы расчета коэффициентов переноса.
1. Montroll E., Shuler K. Studies in nonequilibrium rate processes. I. The relaxation of a system of harmonic oscillators //J. Chem. Phys. 1957. Vol. 26. P. 454-464.
2. Adamovich I., Macheret S., Rich J., Treanor C. Vibrational relaxation and dissociation behind shock waves // AIAA Journal. 1995. Vol. 33. P. 1064-1075.
3. Armenise I., Capitelli M., Colonna G., Gorse C. Nonequilibrium vibrational kinetics in the boundary layer of re-entering bodies //J. Thermophys. Heat Transfer. 1996. Vol. 10. P. 397-405.
4. Candler G., Olejniczak J., Harrold D. Detailed simulation of nitrogen dissociation in stagnation regions // Phys. Fluids. 1997. Vol. 9. P. 2108-2117.
5. Giordano D., Bellucci V., Colonna G., Capitelli M., Armenise I., Bruno C. Vibrationally relaxing flow of N past an infinite cylinder //J. Thermophys. Heat Transfer. 1997. Vol. 11. P. 27-35.
6. Kustova E. V., Nagnibeda E.A. Transport properties of a reacting gas mixture with strong vibrational and chemical nonequilibrium // Chem. Phys. 1998. Vol. 233. P. 57-75.
7. Colonna G., Capitelli M., Tuttafesta M., Giordano D. Non-arrhenius NO formation rate in onedimensional nozzle airflow //J. Thermophys. Heat Transfer. 1999. Vol. 13. P. 372—375.
8. Capitelli M., Ferreira C., Gordiets B., Osipov A. Plasma Kinetics in Atmospheric Gases. Vol. 31 of Springer series on atomic, optical and plasma physics, Springer-Verlag, Berlin, 2000.
9. Kim J., Boyd I. State-resolved master equation analysis of thermochemical nonequilibrium of nitrogen // Chem. Phys. 2013. Vol. 415. P. 237-246.
10. Kang S.H., Kunc J.A. Molecular diameters in high-temperature gases //J. Phys. Chem. 1991. Vol. 95. P. 6971-6973.
11. Gorbachev Yu. E., Gordillo-Vazques F. J., Kunc J. A. Diameters of rotationally and vibrationally excited diatomic molecules // Physica A. 1997. Vol. 247. P. 108-120.
12. Gordillo-Vazquez F.J., Kunc J.A. Rotational—vibrational levels of diatomic molecules represented by the Tietz—Hua rotating oscillator //J. Phys. Chem. 1997. Vol. 101. P. 1595-1602.
13. Gordillo-Vazquez F. J., Kunc J. A. Radial probability of atoms in diatomic molecules represented by the rotating Morse and Tietz—Hua oscillators // J. Mol. Structure. 1998. Vol. 425. P. 263-270.
14. Kustova E. V., Kremer G. M. Effect of molecular diameters on state-to-state transport properties: the shear viscosity coefficient // Chem. Phys. Lett. 2015. Vol. 636. P. 84-89.
15. Kunova O., Kustova E., Mekhonoshina M., Nagnibeda E. Non-equilibrium kinetics, diffusion and heat transfer in shock heated flows of N2/N and O2/O mixtures // Chem. Phys. 2015. Vol. 463. P. 70-81.
...