Комплексное прямое Монте-Карло и кинетическое моделирование неравновесных процессов во фронте ударной волны.



Complex direct monte-carlo and kinetic modelling of nonequilibrium processes in shock wave front

The aggregate of nonequilibrium processes described peaks of radiation and free electrons concentration experimentally
observed in the front of weak shock waves in inert gases with small admixture of heavy iron pentacarbonyl
molecules Fe(CO)5 are studied in this work.
Analysis of pair collision energy distribution function in the zone of transitional relaxation of shock wave was performed
using direct Monte-Carlo simulation method. It was shown that distribution function exceeding in compare
with equilibrium state behind the shock wave front may amount 104 – 105. Subsequent analysis of non-adiabatic
processes probabilities had shown that this effect can’t describe experimental data.
For analysis of kinetic of formation of excited and ionized iron clusters forming behind the shock wave front a numerical
modeling was performed using developed detailed kinetic scheme. Results of calculation are in good agreement
with the experimental dependences.

Александр Викторович Еремин, Александр Всеславович Дракон, С. В. Куликов

Том 7, 2008 год



В работе рассмотрена совокупность неравновесных процессов, описывающих пики излучения и концентрации свободных электронов, наблюдаемые во фронте слабых ударных волн в инертных газах, содержащих
малую примесь тяжелых молекул пентакарбонила железа Fe(CO)5.
Анализ функции распределения энергий парных соударений в зоне поступательной релаксации ударной вол-
ны проведен с использованием методов прямого статистического моделирования Монте-Карло. Показано,
что перехлест функции распределения по отношению к равновесию за фронтом может достигать 104 – 105.
Последующий учет вероятностей неадиабатических переходов при однократных соударениях показал, что
данный эффект не позволяет описать экспериментальные данные.
Для анализа кинетики образования возбужденных и ионизованных кластеров железа, формирующихся за фронтом волны, разработана детальная кинетическая схема процесса и проведено моделирование в программе
ChemKin. Результаты расчетов хорошо согласуются с закономерностями, наблюдаемыми в экспериментах.

Александр Викторович Еремин, Александр Всеславович Дракон, С. В. Куликов

Том 7, 2008 год



1. Зельдович Я.Б., Генич А.П., Манелис Г.Б. // Доклады Академии Наук СССР, 1979, т. 248, №2, c. 349-351.
2. Вагнер Х.Г., Деппе Й., Емельянов А.В., Еремин А.В., Заслонко И.С., Шумова В.В. // ДАН. Химия, 2001, Т. 379, № 1, С. 63−68.
3. Великодный В.Ю., Емельянов А.В., Еремин А.В. // ЖТФ. 1999. Т. 69. Вып. 10. С. 23–33.
4. Генич А.П., Куликов С.В., Манелис Г.Б., Сериков В.В., Яницкий В.Е. // Ж. вычисл. матем. и матем. физ., 1986, 26, №12, С. 1839−1854.
5. Никитин Е.Е. Теория элементарных атомно-молекулярных процессов в газах. -М.: Химия, 1970.
6. Г.Ф. Карабаджак и др. (ЦНИИМаш) Докторская дисс. Май 2007 г.
7. Giesen A., Kovalik A., Roth P. // Phase Transitions, 2004, Vol. 77, №. 1–2, pp. 115–129.
8. Eremin A., Gurentsov E., Schulz C. // NEPCAP-3, Sochi, June 2007.
9. Смирнов В.Н. // Кинетика и Катализ, 1993, т. 34, №4, с. 591-598.
10. Chen J.L., Wang C.S., K.A. Jackson, M.R. Pederson // Phys. Review B., 1991, V.44, №12, pp. 6558−6561.