Расчет инициации горения в осесимметричных каналах на основе кинетической модели многостадийной химической реакции водородо-воздушной смеси

Многочисленные исследования посвящены воспламенению водородо-воздушной смеси и переходу горения в детонацию. Преимущество водорода как моторного топлива заключается в детонационном топливном цикле, более энергетически выгодном, по сравнению с обычным топливным циклом. В связи с этим задача конструирования водородного детонационного двигателя является крайне актуальной. В лабораторных экспериментах воспламенение водорода наблюдается инициацией электрическим разрядом или как результат сжатия в каналах с сужениями. Большое число работ посвящено численному моделированию течений реагирующих смесей в каналах с сужениями, а также инициации воспламенения при возникновении горячего пятна (как модели электрического разряда). В настоящей работе проведена модификация алгоритма, предложенного автором ранее для расчета инициации горения водородо-воздушных газовых смесей, заключающаяся в использовании строго противопоточного варианта разностной схемы Чакраварти–Ошера, позволившей устранить в расчетах осцилляции концентраций компонент смеси. По этой модификации разностного алгоритма проведены расчеты воспламенения водородо-воздушных смесей в каналах с сужениями и угловыми точками, в окрестности которых происходит инициация горения. Для расчета уравнений кинетики использовались два алгоритма, основанные: первый — на решении полной системы уравнений кинетики и второй — на теории ветвящихся цепных реакций Н. Н. Семенова. Проведенные расчеты демонстрируют пригодность алгоритма ветвящихся цепных реакций для расчета задач инициации горения и детонации водородо-воздушных смесей в полной постановке многостадийных химических реакций.