Эффекты поступательной неравновесности во фронте ударной волны

1. Эффекты поступательной неравновесности во фронте ударной волны Великодный В.Ю., Московский физико – технический институт (технический университет), Долгопрудный, Россия E-mail: Vvelikodny@mail.ru

2. Постановка задачи Система кинетических уравнений Больцмана: Функция распределения: Граничные условия: Соотношения Рэнкина-Гюгонио:

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

4. Система уравнений (2)

5. Результаты (1) Сравнение скоростей и температур компонент смеси с данными, полученными в работе Куликов С.В. Поступательная неравновесность трехкомпонентного газа во фронте ударной волны// МЖГ. 1997. №4. С. 171-178 Ряд 1, 3, 5 – работа Куликова С.В. Ряд 2, 4, 6 – данная работа Ряд 1, 2, 3 – данная работа Ряд 4, 5, 6 – работа Куликова С.В. Параметры: M=4, потенциал –«твердые сферы»,

6. Результаты (2) Сравнение скоростей и температур компонент смеси с данными, полученными в работе Ruyev G.A., Fomin V.M., Shavaliev M. Sh. Shock – Wave Structure in a Ternary Disparate – Mass Mixture. : In. Rarefied Gas Dynamics. Weinheim, New York, Basel, Cambrige.: VCH Verlagsgesellshaft mbH., Ed. By Alfred E. Beylich. 1991. Vol. 1., pp. 183-190. Ряд 1, 2, 3 – данная работа Ряд 4, 5, 6 – работа Шавалиева М. Ряд 1, 2, 3 – данная работа Ряд 4, 5, 6 – работа Шавалиева М. Параметры: M= 3, потенциал –«твердые сферы»,

7. Результаты (3) Сравнение изменения профилей парциальных концентраций компонент смесиGmurczyk A. S., Tarczynski M., Walenta Z.A. Shock wave structure in the binary mixtures of gases with disparate molecular masses: In. Rarefied Gas Dynamics. Ряд 1, 2 – данная работа Ряд 3, 4 – эксперимент Параметры: M= 3,61, потенциал –

8. Результаты (4) Зависимость скачка температуры от числа Маха Ряд 1 – компонент 1, Ряд 2 – компонент 2 Параметры: потенциал –«твердые сферы»,

9. Результаты (5) Изменение параллельных составляющих температуркомпонент смеси в системе координат, связанной с собственной скоростью компонентов смеси Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

10. Результаты (6) Изменение параллельных составляющих температуркомпонент смеси в системе координат, связанной с центром масс компонент Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

11. Результаты (7) Изменение параллельных составляющих температур компонент смеси в системе координат связанной с центром масс компонент Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

12. Результаты (8) Изменение перпендикулярных составляющих температур компонент смеси Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

13. Результаты (9) Изменение тензора напряжений компонент смеси в собственных системах координат компонент Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

14. Результаты (10) Изменение тензора напряжений компонент смеси в системе координат центра масс смеси Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

15. Результаты (11) Изменение теплового потока компонент смеси в собственных системах координат компонент Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

16. Результаты (12) Изменение теплового потока компонент смеси в системах координат центра масс смеси Ряд 1 – компонент 1, Ряд 2 – компонент 2, Ряд 3 – компонент 3 Параметры: М = 3, потенциал –«твердые сферы»,

17. Заключение • Разработана методика решения системы кинетических уравнений Больцмана, основанная на обобщении методов Тамма – Мотт- Смита, позволяющая учитывать высокоскоростные хвосты функции распределения и описывать высокопороговые физико – химические процессы в ударных волнах при учете в правой столкновительной части кинетических уравнений неупругих членов. • 2. Проведено тестирование методики решения системы уравнений Больцмана для трехкомпонентных смесей путем сравнения с известными теоретическими расчетами на основе различных подходов и в вырожденном случае – с экспериментом ( для исчезающее малой концентрации третьего компонента). Проведена сравнительная характеристика результатов. • 3. Проведены расчеты высших моментов от функции распределения, параллельных и перпендикулярных парциальных температур, парциальных тепловых потоков, вязких напряжений. Проведен анализ поведения профилей этих величин во фронте ударной волны. • 4. Проведены расчеты различных характеристик, характеризующих поступательную неравновесность во фронте ударной волны в трехкомпонентной газовой смеси, при разных числах Маха и потенциалах взаимодействия частиц. Проведен анализ результатов.

18. Эффекты поступательной неравновесности

19. ЭФФЕКТЫ ПОСТУПАТЕЛЬНОЙ НЕРАВНОВЕСНОСИ

20. Выводы: • объяснить «аномально» малые времена индукции цепной реакции H2 +O2 в случае инициации этой смеси слабой ударной волной, наблюдаемые в ряде экспериментов [6,8,9], • неравновесное излучение в смесях CS2, SO2/Kr, Xe (компоненты близкие по массе) [11], • « аномально» высокая ионизация Ar во фронте ударной волны умеренной интенсивности (М=4.4-7), описанная в работах [12,13]. • В результате теоретических исследований показана возможность влияния эффекта поступательной неравновесности на время задержки ( индукции) цепных реакций на рис. 1 приведено сравнение с экспериментом [8]. • Анализировались так же и другие эксперименты и теоретические модели [9]. Например, время индукции цепной реакции H2 /O2 + инертные примеси в диапазоне температуры за фронтом ударной волны 720-780К и давлении 10-100 атм. в присутствии радикалов O (1%) по модели [9] совпадает со случаем присутствия 0,2% молекулы O2(b1Σg+) в состоянии с возбужденным электроном по расчетам данной работы. Однако для получения радикала О необходимо потратить энергию ~6 ev, а для получения O2(b1Σg+) – только ~1.6 эВ. Спрашивается откуда взять такое количество радикалов O при столь малых интенсивностях ударной волны? С позиций эффектов поступательной неравновесности во фронте слабой ударной волны можно объяснить «аномальную» задержку воспламенения смеси H2 + O2 + инертные примеси , наблюдаемую в работе [9].

21. литература • 1. Bird G.A. Collisdion rates and collisional energy distributions within shock waves//Proc. Of Tenth. Int. Shoch Tuybe Symp. Japan. 1975. P 284-287. • 2. КолышкинИ.Н., ЭндерА.Я., ЭндерИ.А. Разложение по максвеллианам в задачах монотонной и немонотоннй релаксации: Москва, Наука Молекулярная газодинамика и механика неоднородных сред.1990. С 75-82.. • 3. Velikodniy V.Yu., Emel’yanov A.V., Eremin A.V. Nonadiabatic iodine molecule excitation in the shock wave translational nonequilibrium a region .// Zhurnal Technicheskoi Fiziki. 1999. V.89. №10. P. 23-33. • 4. Velikodniy V.Yu. Decision of Boltzmann equation for task of structure of shock wave front// XII International conference on the computing mechanics and modern applied program systems, theses of reports. V.1, 2003. P.144-145. • 5. Velikodniy V.Yu. The effect of translational nonequilibrium in a shock wave on initiation of the H2 +Cl2 reaction// Chem.Phys. Reports 2001. Vol.19(4). P.787-792. • 6. Divakov O.G., Velikodnyi V.Yu., Eremin A. V., Ziborov V. S. Nonequilibrium ignition of H2/O2 deluter mixtures in the weak shock wave front// IV Int. Schol-Seminar. Nonequilibrium processes and their Application. 1998. Minsk. Belarus. P. 155-160. • 7. Velikodniy V. Yu., Bityurin V.A. The effect of translational nonequilibrium on the kinetic of physicochemical conversion in the shock wave front// Chem. Phys. Reports. 1997. Vol.16(9). P. 1521-1531. • 8. Divakov O. G., Eremin A.V., Ziborov V.S., Fortov V.E. Non-equilibrium ignition oxygen -hydrogen mixtures in front of a weak shock wave // Dokl. RAS. 2000. V.373, № 4. p.487-490. • 9. Gel'fand B.E., Medvedev S.P., Homik S.V., Popov O.E., Kusharin A.J., Agafonov G.L. Spontaneous ignition of mixtures hydrogen-oxygen at high initial pressure // Dokl. RAS, 1996, v.349, №4. p.482-485 • 11. Genich A.P., Kulikov S.V., Manelis G.B., Chereshnev S.L. Thermophysiucs of translational relaxation in shock waves in Gases // Sov. Tech. Rev. B Therm. Phys. 1992. v.4. P. 1-69. • 12. Bazhenova T.V., Gvozdeva L.G., Lobastov J.S., Naboko I.M., Nemkov R.G., Predvoditeleva O.A. Shock waves in real gases. M.: Nauka, 1968, p.189 • 13. Kozlov P.V., Losev S.A., Romanenko J.V. Translational non–equlibrium in front of shock wave in a mixtuer of argon and helium // Letters in Zhurnal Technicheskoi Fiziki 2000, v.26. Issue 22. p.69-75. • 14. Bird G.A. The velocity distribution function within a shock wave// J. Fluid Mech. 1967. vol. 30. Part.3. P. 479 –487. • 15. Kogan M. N. Rarefied Gas Dynamics. M.: Nauka. 1967. 440 p.

22. КЛАСТЕРНЫЙ ТЕРМОЯДЕРНЫЙ СИНТЕЗ (КРИТИЧЕСКИЙ АНАЛИЗ ПУБЛИКАЦИЙ)

23. ПРИНЦИПИАЛЬНАЯ СХЕМА ДВИГАТЕЛЬНОЙ УСТАНОВКИ И ЭЛЕМЕНТОВ ЕЕ КОНСТРУКЦИИ

24. collimatedbeam of neutrons lolock - guillotine - blasting blasting charge charge ice porousiceice nozz mechanism of blasting chargesupply Cf252 Bbarrel of СХЕМА РАБОТЫ УСТРОЙСТВА ПРИ СРАБЫТАВАНИИ ЗАРЯДА

25. Проведение предварительных испытаний Взрыв проволочки из молибдена. Одно из тестовых испытаний. Крепления для проволочки

26. equipment for production of porous ice water shield people, load equipment for production of blasting charge engine ПРИНЦИПИАЛЬНАЯ СХЕМА КОСМИЧЕСКОГО КОРАБЛЯ

27. ВЫВОДЫ • Создание детонационного космического двигателя позволит в 5 10 раз сократить время полета к Марсу, Юпитеру, Ио при существенно большей полезной нагрузке (по сравнению с ЖРД). • Важное преимущество детонационного ядерного двигателя перед ядерными двигателями других конструкций состоит в том, что в данном случае не так остро стоит проблема сброса избыточного тепла в космос. Последнее позволяет существенно уменьшить вес конструкции и удешевить стоимость космического корабля • Предварительные расчеты показали, что можно получить тягу P = 400 т и удельный импульс J = 5400 с. Это существенным образом превышает параметры для ЖРД по J и электрореактивных ядерных двигателей по P • Важным преимуществом данного детонационного ядерного двигателя перед ЖРД и ядерными двигателями других конструкций является то, что вместо водорода и кислорода для ЖРД, водорода для целого ряда ядерных двигателей традиционных конструкций в нем в качестве вещества - движетеля предполагается использовать пористый лед. Это связано с тем, что в детонационном двигателе такой конструкции удельный импульс Isp~ ()-1/4, а не Isp~()-1/2 (где  - молекулярный вес) как в ЖРД и ядерных двигателях традиционных конструкций [9]. Поэтому вместо жидкого водорода, который имеет малую плотность =0.172 г/см3 возможно использовать пористый лед, который можно готовить прямо на борту летательного аппарата. • Возможен вариант двигателя для корабля с целью перевозки людей с более высокой скоростью, но с меньшей полезной нагрузкой, и для перевозки грузов с меньшей скоростью, но с большей полезной нагрузкой. Крейсерская скорость в первом случае составляет V = 55 -120 км/c во втором V = 40-25 км/c/

29. Расчет развала молекулы бензола

30. Ударная волна в воде

31. Рис.3. Изменение отношения длины свободного пробега к толщине ударной волны в зависимости от числа Маха в плотном газе (b=0,8). Сравнение различных подходов.

32. Рис.5. Изменение логарифма отношения констант скоростей реакции разрыва двойной связи в молекуле бензола, посчитанных по формуле (18) - кривая 2, к константе скорости за зоной поступательной неравновесности - кривая 1 и на основе закона Аррениуса .

33. ВЫВОДЫ • Получено решение задачи о структуре ударной волны в реагирующем газе большой плотности и жидкости. • Получено аналитическое выражение для константы скорости физико-химических превращений в зоне поступательной неравновесности во фронте ударной волны для плотных газов и жидкостей. • Показано, что аррениусовская кинетика не пригодна для описания процессов физико-химических превращений в зоне поступательной неравновесности во фронте ударной волны. • Предложен новый метод, обобщающий метод Мотт-Смита для решения задачи о структуре ударной волны в реагирующих плотных газах и жидкостях. • Наиболее вероятным механизмом, по которому может реализоваться разрыв двойной связи C=C в молекуле , является механизм сильных однократных неадиабатических суперстолкновений. • В зависимости от молекулярной структуры, величины порогов энергии активации, интенсивности ударных волн эффекты поступательной неравновесности ( существенное превышение выхода продуктов реакции, наработанных в зоне поступательной неравновесности, над выходом продуктов за этой зоной) могут проявляться или отсутствовать. • Работа поддержана Российским Фондом Фундаментальных Исследований-грант N 01-03-32729 a.

34. Внешний вид установки для обжатия микропористой среды при взрыве проволочки кольцевой конфигурации

35. Гамма-спектры неактивированного и активированного индия Совмещённый гамма-спектр. (Индий ненаактивирован – спектры фона и индия практически совпадают.) 28.10.2005. Совмещённые гамма-спектры индия– синяя кривая и фона лиловая кривая. Жёлтым цветом выделена энергетическая область расположения фотопика Е=417кэВ. Индий активирован.

36. Результаты расчета потока нейтронов f – плотность потока нейтронов, Ф=4πR2f* - полный поток «в четыре пи», где R- расстояние от мишени до эпицентра реакции