Вычислительные технологии
Том 11, часть 2, Специальный выпуск, 2006
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ СТАЦИОНАРНЫХ ДЕТОНАЦИОННЫХ ВОЛН*
А. В. Троцюк, А.Н. Кудрявцев, М. С. Иванов Институт теоретической и прикладной механики им. С.А. Кристиановича СО РАН, Новосибирск, Россия e-mail: trotsyuk@itam.nsc.ru
Numerical simulations have been performed to study regular and Mach reflections of oblique shock waves in a steady supersonic flow of a homogeneous combustible gas mixture for different inflow Mach numbers Min. The main objective of the present study is to investigate numerically the feasibility of “stationary” detonation at the Mach reflection of oblique shock waves. It has been found that for some values of Min there may exist a standing Mach stem with a system of unsteady transverse waves on its front.
Введение
Одной из наиболее интригующих возможностей при разработке новых перспективных прямоточных воздушно-реактивных двигателей является использование камеры сгорания со стационарной (стабилизированной) детонационной волной (ДВ). Этот подход имеет существенные преимущества, однако нестабильная природа взаимодействия зоны реакции и ведущей ударной волны в ДВ осложняет его практическую реализацию. Быстрая скорость горения или переработка вещества, обычно в десятки тысяч раз большая, чем во фронте пламени, может предоставить ряд преимуществ для детонационных реактивных двигателей, например их большую компактность и эффективность. Детонационные двигатели должны теоретически давать наивысшую эффективность в превращении потенциальной химической энергии топлива в кинетическую энергию продуктов сгорания [1].
Предложения по использованию детонационной волны для создания тяги появляются в литературе в 1940-х годах. Уже тогда были выдвинуты две различных концепции. Использование стационарных (стабилизированных) ДВ предложено в работах [1, 2]. В работе [3] перемежающийся (или пульсирующий) режим с нестационарными ДВ осуществлен в реактивном устройстве. Значительное число экспериментальных, теоретических и численных исследований этих двух концепций проведено в течение прошедших лет, и в настоящее время наблюдается новый интерес к данной проблеме [4, 5].
* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 05-03-32131) и Программ фундаментальных исследований РАН (№ 17 и № 19, 2003 г.).
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2006.
К настоящему времени очень мало известно о стационарных детонационных волнах. Экспериментальные доказательства их существования немногочисленны и противоречивы [4], а их теория далека от завершения. К сожалению, большинство известных экспериментальных результатов может быть должным образом интерпретировано скорее как горение за ударными волнами (УВ), чем собственно явление детонации.
Основной целью данного исследования является численное изучение осуществимости “стационарной” детонации при маховском отражении (М11) наклонных ударных волн, генерируемых системой из двух клиньев, в стационарном сверхзвуковом потоке водороднокислородной смеси. По сравнению с отражением наклонных ударных волн в инертной среде [6-8], взаимодействие наклонных ударных и детонационных волн между собой или с твердыми поверхностями еще недостаточно полно исследовано.
В настоящих исследованиях использован принцип расщепления по физическим процессам, предложенный в работах Н.Н. Яненко, а именно при нахождении потоков консервативных переменных через грани контрольного объема при использовании конечнообъемной схемы для моделирования химически реагирующих течений применялся алгоритм приближенного решения задачи Римана для химически замороженного газа. Учет протекания химических реакций производился на следующем этапе, когда вычислялись новые значения переменных на верхнем временном слое.
1. Формулировка проблемы и физическая модель
Моделирование проведено для течения горючей смеси 2Н2 + 02 при начальном давлении и температуре р0 = 0.2 бар, Т0 = 298.15 К и различных числах Маха набегающего потока М[п. Рассматривалось двумерное течение смеси над телом сжатия, которое представляло собой клин с углом в и длиной поверхности сжатия ш, в канале с высотой входной части Н и выходной части К (рис. 1,а). В отличие от моделирования отражения УВ в инертной среде, в случае химически реагирующих потоков существуют дополнительная характерная скорость, а именно скорость детонационной волны в режиме Чепмена — Жуге (С-Л), и дополнительный характерный линейный размер. Поперечный размер детонационной ячейки а0 часто используется в качестве этого размера. Другой часто используемой длиной является расчетная толщина зоны реакции идеализированной одномерной стационарной детонационной волны.
О 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.Э 0.9 1 1.1 х, см 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 х, см
а б
Рис. 1. Численная шлирен-визуализация поля течения при М^п = 7.0, в = 23.5° и Н/т = 0.4: а — регулярное отражение; б — стационарное маховское отражение.
Для указанной смеси и в данных начальных условиях мы имеем Дел = 2757 м/с и Мел = 5.13. Формирование стационарно стоящей в потоке ДВ возможно лишь при соблюдении условия, что число Маха набегающего потока больше Мел- При невыполнении этого условия, в случае возникновения детонации в канале, детонационная волна будет распространяться вверх по потоку, поэтому образование стационарной конфигурации невозможно.
Течение химически реагирующей газовой смеси описывалось нестационарными двумерными уравнениями сохранения массы, импульса и энергии в представлении Эйлера. Вязкость, молекулярная диффузия и эффекты теплопроводности не рассматривались.
Химические превращения в газовой смеси описывались двухстадийной кинетической моделью (включающей в себя индукционную стадию и стадию тепловыделения) [9]. Индукционная стадия моделировалась согласно экспериментальной кинетике [10]. Последующий этап тепловыделения описывался, следуя модели химических реакций при высоких температурах [11, 12] совместно с калорическим уравнением состояния [12], согласованным с выбранной кинетической моделью по второму закону термодинамики. Отметим также, что данный подход позволяет учесть сильное изменение в ходе химических реакций значений средней молярной массы смеси, тепловыделения, удельных теплоемкостей и их отношения. Система определяющих уравнений замыкалась термическим уравнением состояния идеального газа. Применение данных кинетических моделей описано более подробно в работах [13, 14].
На различных сторонах вычислительной области ставились следующие граничные условия (см. рис. 1): условие равномерного сверхзвукового набегающего потока задано на левой границе области; условие сверхзвукового истечения определено на правой границе; условие непротекания положено на верхней границе, рассматривающейся в качестве жесткой стенки; нижняя граница считается линией симметрии, что в случае уравнений Эйлера эквивалентно идеальной твердой стенке, поэтому условие непротекания определено и для этой границы.
2. Численный метод
Полученная гиперболическая система уравнений решалась численно с использованием конечно-объемной схемы и \11’5С1. ТУБ-иптерполяции четвертого порядка точности [15]. Для приближенного решения задачи Римана использовался современный НЬЬС-алго-ритм [16]. Применение данного алгоритма для случая химически реагирующего газа стало возможным при помощи “метода релаксации энергии” [17]. Этот метод позволяет избегать трудностей, возникающих при численном решении задачи Римана для среды со сложным нелинейным уравнением состояния (включая случай с переменным показателем адиабаты 7). Интегрирование по времени осуществлялось со вторым порядком точности путем использования недавно разработанных аддитивных явно-неявных методов Рунге — Кут-
= 0 . 2 . . . 0 . 3
При моделировании использовалась адаптированная к границам расчетной области четырехугольная сетка с числом ячеек около Мх = 600 и N = 200 в направлении х и у соответственно. Для проверки независимости получаемых численных результатов от пространственного разрешения по сетки некоторые вычисления повторены с удвоенными значениями ячеек по каждой координате. Обнаружено, что результаты расчетов, выполненные на грубой и мелкой сетках, очень хорошо согласуются между собой.
3. Результаты расчетов
В работе проведен анализ регулярного (ЕЕ) и маховского отражений наклонных ударных волн с использованием метода ударных и детонационных поляр. Для получения детонационной поляры использована модель детонационной волны с мгновенной химической реакцией на фронте для некой модельной газовой смеси с 70 = 1.4 Тел = 1.2 и Мел = 5.13. Показано, что существование как регулярного, так и маховского типа отражений возможно при одинаковых параметрах течения. Эта проблема исследована численно в последующем двумерном моделировании течения.
На рис. 1,о приведена численная шлирен-визуализация [19] регулярного отражения наклонной ударной волны при угле клипа в = 23.5° и числе Маха набегающего потока Мщ = 7.0. Структура течения в данном случае аналогична структуре регулярного отражения в инертной среде [6-8], за исключением того, что за отраженной волной имеется зона тепловыделения. Высота входной части канала Н = а0 = 0.64 см, т. е. выбрана равной размеру детонационной ячейки в свободно распространяющейся нестационарной многофронтовой детонационной волне [13].
Стационарное маховское отражение получено для чисел Мщ = 7.0 и Мщ = 8.0. В этих случаях ножка Маха (МН) представляет собой гладкую ударную волну, подобную У В в инертной среде, с последующей за ней зоной индукции и зоной тепловыделения (см. рис. 1,6). Поперечные волны на фронте МН отсутствуют, ножка является стационарной и стоящей в потоке пережатой детонационной волной. Исследовано влияние величины тепловыделения в зоне реакции на величину (высоту) ножки Маха (рис. 2, о). Для этого проведено моделирование течения модельной смеси с константами Ер, меньшими по величине, чем значение константы Ер для нормальной смеси 2Н2 + 02. Эта константа входит в модель кинетики тепловыделения [11, 12], используемую нами, и определяет максимальное значение тепловой энергии, которая может выделиться в зоне реакции; Ер = 0 соответствует инертной среде. Высота ножки Маха растет с увеличением тепловыделения
в
нормальной смеси.
На рис. 3, а приведен случай маховского отражения для Мщ = 5.5. Структура течения за ножкой Маха подобна структуре многофронтовой (ячеистой) свободно распространяющейся детонационной волны [13]. На фронте маховской ножки имеются нестационарные поперечные волны, которые периодически отражаются от плоскости симметрии и области течения в окрестности тройной точки. Движение этих поперечных волн вдоль ножки Маха приводит к существенному изменению ее формы.
В отличие от маховского отражения в инертной среде, ножка Маха для случая химически реагирующего потока оказалась нестационарной при Мщ = 5.5. После формирования маховского отражения детонационный фронт двигается вверх или вниз по потоку. Наши исследования показали, что при фиксированных значениях Мщ и в отношение К/ш является управляющим параметром. Моделирование проведено при в = 25°. При величине К/ш вплоть до 0.7 ножка Маха всегда двигалась вперед против потока вплоть до входной границы. В этом случае в расчетах получается газодинамическая картина течения, аналогичная той, что имеет место в незапущенном сверхзвуковом диффузоре с непрерыв-
К/ш
критическое значение, то ножка Маха, возникнув, затем перемещалась вниз по потоку и исчезала, и мы фиксировали регулярное отражение в качестве финальной стационарной
К/ш = 0.8
Рис. 2. Высота ножки Маха Нмэ при М^ = 7.0 и Н/т = 0.4 как функция: а — величины тепловыделения в модельной смеси (для в = 24.5°); б — угла клина (для нормальной смеси).
Рис. 3. Численная шлирен-визуализация поля течения при М^ = 5.5, в = 25° и Н/т = 0.6: а — нестационарная многофронтовая ножка Маха при Н = ао] б — стационарная ножка Маха с подавленной многофронтовой структурой при Н = ао/3.
Рис. 4. Численная шлирен-визуализация поля течения со стационарной многофронтовой ножкой Маха при Min = 6.0 и в = 25°: а — h/w = 0.5; 6 — h/w = 0.6.
кое значение h/w находится между 0.7 и 0.8. Чем ближе значение h/w к этой критической величине, тем с меньшей скоростью перемещается ножка Маха вверх или вниз по потоку. Основываясь на этих данных, мы считаем, что ножка Маха находится в состоянии
h/w
Дальнейшие исследования показали, что стационарная ножка Маха может существовать с системой нестационарных поперечных волн на своем фронте, показанной на рис. 3, а. Результаты моделирования для Min = 6.0 приведены па рис. 4. Несмотря па наличие поперечных волн, их движение вдоль фронта ножки и периодическое изменение его формы, ножка Маха в среднем представляет собой стационарно стоящую волну, которая слабо осциллирует около некоторой точки равновесия. Сравнивая структуру фронта ножки на рис. 3, о и рис. 4, можно заметить, что во втором случае поперечные волны менее выражены и имеют меньшую длину в направлении осп x (особенно на рис. 4, б). Это объясняется большей степенью пережатия (т. е. Min/Mcj) детонационной волны (ножки Маха) для Min = 6.0. Предположено, что существование интенсивных и протяженных поперечных волн па фронте ножки для Min = 5.5 приводит к нестационарное™ маховского отражения.
Для проверки справедливости этого предположения моделировалось маховское отражение для двух каналов с одинаковыми в и h/w и одинаковыми параметрами течения, но различными значениями H. Для случая H = а0 (см. рис. 3, а) получалась нестационарная многофронтовая маховская ножка, которая двигалась вверх по потоку до входной границы. Однако для H = а0/3 (см. рис. 3,6) структура течения оказалась качественно отличной. Стационарная и неподвижная ножка имела гладкий фронт, поперечные волны отсутствовали полностью (см. рис. 4,6). Изменение пространственного масштаба течения привело к формированию стационарной ножки Маха. Характерный размер течения (высота ножки Маха) во втором случае оказался меньше, чем характерный внутренний масштаб в детонационной волне (т. е. поперечный размер ячейки на фронте пережатой детонационной волны). Это и явилось причиной подавления развития неустойчивости фронта ножки и как следствие стационарности структуры течения в целом. Таким образом, исследован эффект интенсивности поперечных волн на ножке Маха на стационарность маховского режима отражения.
Влияние начальных условий на конечные стационарные ударно-волновые конфигурации, получающиеся в расчетах, а также результаты исследования перехода регулярного режима отражения к маховскому под воздействием локальных возмущений в натекающем потоке описаны в [14].
Нами также изучено течение водородно-воздушных смесей. Размер и геометрия вычислительной области выбраны соответствующими тестовой секции высокотемпературной аэродинамической трубы [20] с числом Маха стационарного потока Min = 3.0. В этих расчетах в = 23°, h/w = 0.8655, h = 5.25 см и H = 7.62 см. Эти величины определены с максимально возможной точностью по фотографиям и схемам из [20]. Для бедной водородно-воздушной смеси 0.4H2+4.772Air (коэффициент избытка воздуха 5, Mcj = 3.04) число Маха набегающего потока и параметры точки торможения выбраны идентичными условиям эксперимента [20], а именно Min = 3.15, p00 = 7.485 бар, T00 = 1033.15 К. Для выбранных параметров потока и размеров вычислительной области не было зарегистрировано химических реакций как за отраженной волной, так и за ножкой Маха в численных экспериментах по исследованию как регулярного, так и маховского отражения. Такой же результат получен и в работах других исследователей. Однако при пятикратном увеличении линейных размеров обнаружен очень интересный режим маховского отражения с
а
б в
Рис. 5. Осциллирующая ножка Маха в бедной водородно-воздушной смеси 0.4H2 + 4.772Air для в = 23°, h/w = 0.8655 и Min = 3.15: а — траектория ножки Маха, частота пульсаций около 840 Гц; численная шлирен-визуализация поля течения при позиции фронта: б— в точке 1; в — в точке 2. Здесь tw = w/c0, где с0 — скорость звука в набегающем потоке.
сильно переменной по высоте (пульсирующей) ножкой Маха (рис. 5). Осцилляции положения маховской ножки вверх и вниз по потоку имеют значительную амплитуду и являются строго периодическими во времени. Несмотря на эти осцилляции, фронт ножки Маха всегда остается гладким, какие-либо поперечные волны или возмущения на нем не были зарегистрированы. Этот режим очень напоминает режим одномерных пульсаций детонационной волны, свободно распространяющейся в узком капилляре (околокритический режим галопирующей детонации).
Список литературы
[1] Зельдович Я.Б. К вопросу об энергетическом использовании детонационного горения // Журн. техн. физики. 1940. Т. 10, № 17. С. 1453-1461.
[2] Roy М.М. Propulsion par Statoreacteur a Detonation // Comptes Rendus de TAcademie des Sciences de Paris. 1946. Vol. 222. Feb. 1. P. 31-32.
[3] Hoffmann N. Reaction Propulsion by Intermittent Detonative Combustion // German Ministry of Supply. A1152365 Volkenrode Translation. 1940.
[4] Shepherd J.S. Detonation waves and propulsion // Combustion in High-Speed Flows / Eds T. Jackson, J. Buckmaster, S. Kumar. Dordrecht; Boston; London: Kluwer Acad. Publ., 1994. P. 373-420.
[5] SlSLlAN J.P. Detonation-wave ramjets // Scramjet Propulsion. Progress in Astronautics and Aeronautics, AIAA, Washington, DC. 2000. Vol. 189. P. 823-889.
[6] Ivanov M.S., Gimelshein S.F., Beylich A.E. Hysteresis effect in stationary reflection of shock waves // Phys. of Fluids. 1995. Vol. 7, N 4. P. 685-687.
[7] Ivanov M.S., Vandromme D., Fomin V.M. et al. Transition between regular and Mach reflection of shock waves: new numerical and experimental results // Shock Waves. 2001. Vol. 11, N 3. P. 199-207.
[8] Ben-Dor G., Ivanov M.S., Vasilev E.I., Elperin T. Hysteresis processes in the regular reflection ^ mach reflection transition in steady flows // Progress in Aerospace Sci. 2002. Vol. 38. P. 347-387.
[9] Левин B.A., Коробейников В.П. Сильный взрыв в горючей смеси газов // Изв. АН СССР. Механика жидкости и газа. 1969. № 6. С. 48-51.
[10] White D.R. Density induction times in very lean mixtures of D2, H2, C2H2, and C2H4 with O2 // 11th Symp. (Intern.) on Combustion, Proc., Berkeley, CA, 1966. P. 147-154.
[11] Николаев Ю.А., Фомин П.А. О расчете равновесных течений химически реагирующих газов // Физика горения и взрыва. 1982. Т. 18, № 1. С. 66-72.
[12] Николаев Ю.А., Зак Д.В. Согласование моделей химических реакций со вторым началом термодинамики // Физика горения и взрыва. 1988. Т. 24, № 4. С. 87-90.
[13] ТРОЦЮК А.В. Численное моделирование структуры двумерной газовой детонации смеси H2 — O2 — Ar // Физика горения и взрыва. 1999. Т. 35, № 5. С. 93-103.
[14] Trotsyuk A.V., Kudryavtsev A.N., Ivanov M.S. Numerical simulation of reflection of detonation and shock waves in supersonic flows // Confined Detonations and Pulsed Detonation Engines / Eds G. Roy, S. Frolov, R. Santoro, S. Tsyganov. Moscow: Torus Press, 2003. P. 73-86.
[15] Yamamoto S., Daiguji H. Higher-order-accurate upwind schemes for solving the compressible Euler and Navier — Stokes equations // Comput. Fluids. 1993. Vol. 22, N 2/3. P. 259-270.
[16] Batten P., Leschziner M.A., Goldberg U.C. Average-state Jacobians and implicit methods for compressible viscous and turbulent flows // J. Comput. Phys. 1997. Vol. 137. P. 38-78.
[17] COQUEL F., Perthame B. Relaxation of energy and approximate Riemann Solvers for general pressure laws in fluid dynamics // SIAM J. Numer. Anal. 1998. Vol. 35, N 6. P. 2223-2249.
[18] Shen J.W., Zhong X. Semi-implicit Runge — Kutta schemes for non-autonomous differential equations in reactive flow computations // AIAA Paper. 1996. N 96-1969.
[19] Quirk J.J. A contribution to the great Riemann solver debate // Intern. J. Numer. Methods in Fluids. 1994. Vol. 18. P. 555-574.
[20] Gross R.A., Chinitz W. A study of supersonic combustion // J. of the Aero/Space Sciences. 1960. Vol. 7, N 7. P. 517-524.
Поступила в редакцию 13 апреля 2006 г.