Численное моделирование нестационарных течений реагирующего газа с явным выделением произвольного числа взаимодействующих разрывов тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Северина, Наталья Сергеевна

• Северина, Наталья Сергеевна
• кандидат физико-математических наук
• 2009

Оглавление диссертации кандидат физико-математических наук Северина, Наталья Сергеевна

СПИСОК ОБОЗНАЧЕНИЙ.

ВВЕДЕНИЕ.

Глава 1. Сеточно-характеристический метод для расчета квазиодномерных нестационарных течений реагирующего газа с явным выделением произвольного числа взаимодействующих разрывов.

1. 1 Математическая модель.

1. 2 Алгоритм решения начально-краевой задачи с явным выделением произвольного числа сильных и слабых разрывов.

1.2.1. Определение координат подвижных узлов.

1.2.2. Расчет точки в поле течения.

1.2.3. Расчет точки на контактном разрыве.

1.2.4. Расчет точки на ударной волне.

1.2.5. Расчет точки на границе области.

1.2.6. Расчет точки на характеристике.

1.2.7. Обработка взаимодействий.

1.3. Модельные задачи.

Выводы к главе 1.

Глава 2. Комплекс программ для задания, хранения, визуализации и анализа результатов численного моделирования одномерных нестационарных задач физической газовой динамики.

2.1. Назначение и структура комплекса.

2.2. Задание входных данных для расчетного модуля

2.3. Графическая обработка расчетных данных.

2.4. Особенности программы визуализация.

2.5. База данных термодинамических свойств индивидуальных веществ и кинетических механизмов.

Выводы к главе 2.

Глава 3. Моделирование влияния колебательной релаксации и электронного возбуждения на структуру ударной волны в кислороде.

3.1. Роль термической неравновесности в моделировании химической кинетики на примере диссоциации молекулярного кислорода за фронтом ударной волны.

3.1.1. Математическая модель.

3.1.2. Исследования термической неравновесности за фронтом стационарной ударной волны.

3.1.3. Термическая неравновесность при нестационарном распространении ударной волны.

3.2. Влияние электронного возбуждения на структуру ударной волны в кислороде.

3.2.1. Постановка задачи.

3.2.2. Математическая модель.

3.2.3. Стационарная ударная волна.

Выводы к главе 3.

Глава 4. Моделирование воспламенения водородо - кислородных смесей в ударной волне.

4.1. Задача об инициировании детонации за отражённой ударной волной.

4.2. Моделирование эксперимента В.А. Павлова в ударной трубе.

4.2.1. Экспериментальная установка.

4.2.2. Моделирование экспериментов в стационарной постановке.

4.2.3. Моделирование экспериментов в нестационарной постановке.

Выводы к главе 4.

Введение диссертации (часть автореферата) на тему «Численное моделирование нестационарных течений реагирующего газа с явным выделением произвольного числа взаимодействующих разрывов»

Разработка численных схем для расчёта разрывных решений законов сохранения газовой динамики - тема, интенсивно развивающаяся в течение последних десятилетий. Характерной особенностью течений невязкого газа около тел или в каналах сложной формы является взаимодействие ударных волн и других газодинамических разрывов. При численном решении подобных задач используются два основных типа разностных схем: 1) схемы сквозного счета, 2) схемы с выделением разрывов.

Оба подхода имеют определенные преимущества и недостатки. При использовании методов сквозного счета [22,29,33,38,82-83] все разрывы в исследуемом течении размазываются, т.е. представляются в виде конечных областей изменения параметров потока. Даже самые современные методы "сквозного счета" (TVD [39,40,77-81], UNO [41], ENO [42], WENO [43] и др.) "размазывают" ударную волну и контактные разрывы, особенно разделяющие газы с существенно различающимися свойствами, на 2-3 ячейки разностной сетки. Величина этого размазывания определяется размером расчетных ячеек. В схемах сквозного счета алгоритмы, как правило, весьма просты, единообразно описывают все поле течения и позволяют без серьезной модификации решать широкий круг задач. Однако точность таких схем не всегда бывает удовлетворительной. Кроме того, измельчение расчетной сетки, необходимое в области разрыва, приводит к значительному росту времени счета [29,39-40]. Схемы сквозного счета развиваются в направлении повышения порядка аппроксимации. Они используют отдельные черты схем с выделением разрывов для улучшения качества решения: измельчение сетки в области разрыва при сохранении ее в области непрерывного течения, увеличения порядка аппроксимации в области больших градиентов.

Методы явного выделения разрывов существенно повышают эффективность расчета. Основные принципы построения таких численных методов были разработаны еще 40 лет назад, причем примеры их применения подтвердили их высокую эффективность [5]. Из-за укоренившегося убеждения, что явное выделение поверхностей разрывов в потоке - это очень трудоемкий и плохо формализуемый процесс, число исследований, проведенных с использованием таких методов, крайне незначительно. Как правило, в них реализовано явное выделение ограниченного числа наиболее сильных ударных волн [24].

Связанное с выделением разрывов усложнение алгоритмов оправдывается не только ростом точности, но и сокращением времени счета. Отметим, что основное усложнение алгоритмов в этих подходах относится к прохождению разрыва. Поэтому смягчение основного недостатка, т.е. сложности алгоритма, следует искать на пути упрощения счета в окрестности разрыва.

Методы с явным выделением поверхностей разрывов позволяют легко строить адаптивные разностные сетки, состоящие как из неподвижных, так и подвижных узлов. Последние могут двигаться вдоль траекторий выделяемых сильных и слабых разрывов, траекторий газа (линий тока), подвижных границ и т.п. Это позволяет получать на сравнительно грубых сетках как общую картину течения, так и отдельные ее детали.

В литературе существуют два направления. Первое называется "скачок как граница". Вся область течения разбивается на подобласти, разделенные разрывами, либо границами другого типа. В различных подобластях используются различные сетки. Таким способом почти всегда решались задачи внешней аэродинамики [1-4] и решаются другие задачи [5]. Во втором направлении используется термин "плавающий скачок". При таком подходе разностная сетка фиксируется заранее. Ударные волны и другие разрывы могут двигаться по сетке. Специальный алгоритм должен следить за траекторией разрывов. Это второе направление представляется более универсальным. Оно развивалось в работах [6-11].

В [12] изложен способ построения разрывных нестационарных решений, обеспечивающий выделение всех, имеющих значение для данного расчета, разрывов решения и его производных. В предлагаемой схеме автоматически выполняется достаточное условие устойчивости, что гарантирует монотонность численного решения [15, 18]. Так же, имеет место аппроксимация законов сохранения на всём решении, что обеспечивает локализацию разрывов (отсутствие их «размазки »). По сложности алгоритм находится на уровне схем «сквозного» счёта, не использующих громоздкого сглаживания, и в этом смысле качественно отличается от предлагавшихся ранее способов расчёта с выделением разрыва.

Структура метода в некоторой степени повторяет известный прямой характеристический способ с той разницей, что вместо традиционных квазилинейных используются обобщённые характеристики [15]. Соответственно, роль характеристического преобразования квазилинейных уравнений выполняет распад произвольного разрыва для законов сохранения. Эта процедура используется не только при вычислении потоков и оценки допустимого шага по времени, но и для непосредственного построения сетки и решения. Она заменяет дивергентные законы сохранения системой обобщённо-характеристических соотношений. При этом в каждое соотношение входит только одна функция (плотность, импульс, энергия единицы объёма) и обобщённо-характеристические скорости, общие для всей системы уравнений (скорости обобщённых характеристик, возникающие при распаде произвольного разрыва).

В [13] разработан метод расчета течений невязкого газа с разрывами. Численная схема основана на явном выделении ударных волн и тангенциальных разрывов. Используется простой алгоритм выделения. В области непрерывного течения уравнения Эйлера решаются с помощью известной схемы Мак-Кормака, а прохождение ударных волн осуществляется явно с использованием условий Ренкина-Гюгонио.

Здесь в рамках концепции плавающего скачка предлагается новый подход к численному решению уравнений газовой динамики в окрестности разрывов. Разностная сетка не связана с разрывами, которые "плавают" по сетке в процессе счета. Главный принцип состоит в такой записи исходных уравнений, для которой при переходе через разрыв можно ограничиться разностной аппроксимацией производных непрерывных комплексов искомых функций. Недостающие переменные определяются из подходящих алгебраических соотношений. Движение разрыва по сетке рассчитывается известными методами (см., например, [14]). Можно предположить, что этот принцип является общим для любых газодинамических разрывов, т.е. таких, на которых выполняются определенные алгебраические соотношения между искомыми переменными.

Метод применяется к расчету стационарных сверхзвуковых двумерных течений с ударными волнами и тангенциальными разрывами. Допустимо обобщение на многомерные течения.

Процессы, протекающие во многих реальных установках, могут быть описаны в рамках квазиодномерного нестационарного приближения [23-29]. Построение экономичных схем расчета невозможно без серьезного анализа свойств и особенностей численно решаемых уравнений. Квазиодномерные нестационарные дифференциальные уравнения динамики течений реагирующего газа по типу относятся к гиперболическим и имеют три семейства характеристик [30-35]: звуковые характеристики и траектории газа. Вдоль них выполняются соотношения, которые записываются в полных дифференциалах. Соотношения, выполняемые вдоль звуковых характеристик, связывают дифференциалы скорости, давления газа и времени, коэффициенты при дифференциалах зависят от параметров, характеризующих газ. Вдоль траектории газа выполняются соотношения, число которых на единицу больше, чем количество компонент смеси. Первое из соотношений устанавливает связь между дифференциалами внутренней энергии, плотности газа и времени, причем коэффициенты при дифференциалах тоже зависят от параметров, характеризующих газ. Остальные - отвечают за изменение соответствующей концентрации вдоль траектории, которое может происходить за счет химических реакций. Таким образом, в каждой точке пространства приходится решать три уравнения газовой динамики и систему уравнений химической кинетики, число которых равно числу компонент газовой смеси. Реальные кинетические механизмы могут содержать более десяти компонент и несколько десятков реакций. Причем правые части уравнений химической кинетики содержат большое число экспонент, пропорциональное количеству элементарных стадий. Поэтому при решении уравнений реагирующей газовой динамики обычно девяносто процентов время тратится на интегрирование уравнений химической кинетики [29].

Исходя из вышесказанного, можно сделать вывод, что при расчете параметров газа в качестве сеточных линий в областях, занятых реагирующим газом, целесообразно использовать линии, совпадающие с траекториями газа. Такой выбор расчетной сетки позволяет существенно сократить количество интерполируемых величин, а также использовать высокоточные и высокоскоростные численные методы для решения систем обыкновенных дифференциальных уравнений химической кинетики [36,37].

Актуальность темы

Течения многокомпонентного газа при наличии химических превращений, колебательного и электронного возбуждения молекул реализуются во многих современных технологических устройствах и перспективных энергетических установках. В настоящее время с развитием вычислительной техники существенно возросла роль математического моделирования, как альтернативы и дополнения к физическому эксперименту. Математическое моделирование позволяет выявить детальную структуру течения, которую физически невозможно наблюдать в экспериментальных исследованиях. Разработка высокоточных и экономичных методов, позволяющих моделировать течения газа с физико-химическими превращениями при наличии в поле течения сильных разрывов (ударных волн, контактных разрывов), является актуальной практически важной задачей. Не смотря на развитие многомерных методик моделирования, одномерные модели сохранили свою важность, в частности они позволяют путем сравнения численных и экспериментальных результатов верифицировать модели, описывающие неравновесное протекание химических превращений в газовой фазе, которые в дальнейшем могут быть использованы при моделировании работы реальных перспективных установок, например детонационного двигателя.

Пели работы

Разработка вычислительных алгоритмов квазиодномерных нестационарных течений реагирующего газа с явным выделением взаимодействующих разрывов.

Разработка комплекса программ, позволяющего исследовать газодинамические течения с ударными и детонационными волнами, выделять их тонкие структуры, отслеживать распространение возмущений, определять места зарождения газодинамических разрывов.

Проведение широкомасштабного вычислительного эксперимента для получения высокоточных решений задач физической газовой динамики, а также с целью исследования неравновесных физических процессов, таких как: химическая кинетика, колебательная релаксация, детонация и другие.

В соответствии с вышесказанным в диссертации поставлены и решены следующие задачи:

1. Разработка вычислительных алгоритмов и комплекса программ для численного моделирования квазиодномерных нестационарных течений для моделирования многокомпонентного произвольного числа многокомпонентного реагирующего газа, с явным выделением произвольного числа взаимодействующих разрывов.

2. Разработка структуры хранения выходных данных и подсистемы визуализации результатов численного решения гиперболических двумерных разрывных задач математической физики на подвижных неструктурированных сетках, в том числе и для случая пересекающихся сеточных линий.

3. Исследование влияния колебательной релаксации и электронного возбуждения молекул на структуру ударной волны в кислороде.

4. Исследование процессов инициирования и развития детонации в ударных трубах. Проведение численного моделирования серии экспериментов по измерению задержек воспламенения сильно разбавленных аргоном водородо-кислородных смесей (3-5%) в ударных трубах за падающими ударными волнами.

Научная новизна

1. Разработаны вычислительные алгоритмы и программный комплекс, позволяющие решать одномерные нестационарные уравнения физической газовой динамики с точным выделением областей непрерывного и разрывного течений. Сеточными линиями в разработанной методике моделирования являются: траектории газа, ударные волны, контактные разрывы, характеристики и фиксированные точки.

2. Разработаны и апробированы алгоритмы, решающие задачи, которые возникают при пересечении сеточных линий друг с другом, в том числе решена задача о зарождении ударной волны, в результате пересечения характеристик одного семейства.

3. На примере расчета течения за сильной ударной волной в кислороде показана применимость разработанных методик к численному моделированию течений с колебательной релаксацией и электронным возбуждением.

Исследована зависимость длин зон температурной и химической неравновесности от числа Маха ударной волны.

4. На примере расчета нестационарного течения в ударной трубе, заполненной горючей смесью водород-кислород-аргон показана применимость разработанных методик к численному моделированию течений при наличии детонации. Получена детальная картина течения, которая включает: распад разрыва на границе инертный газ — горючая смесь, образование и распространение волны воспламенения, взаимодействие волны воспламенения с ударной волной, образование нестационарной пересжатой детонационной волны, переход от нестационарной к стационарной пересжатой детонационной волне. Исследована зависимость времени задержки воспламенения горючей смеси от начального давления и степени разбавления аргоном.

5. Разработаны алгоритмы и комплекс программ для прецизионной визуализации результатов численного моделирования на существенно нерегулярных сетках: сеточные линии могут пересекаться, исчезать, зарождаться в процессе расчета. Комплекс программ позволяет строить временные развертки течения, графики параметров течения в различные моменты времени, а также вдоль произвольных сеточных линий.

Апробация работы

Материалы диссертации докладывались и обсуждались

• На XV и XVI международных конференциях по вычислительной механике и современным прикладным программным системам (ВМСППС) (Алушта, 2007, 2009)

• На VII Международной конференции по неравновесным процессам в соплах и струях (NPNJ'2008) (Алушта, 2008)

• На семинаре НИИМЕХ МГУ по физической газовой динамике (Москва, 2009)

Публикации

Основное содержание диссертации отражено в 8 опубликованных работах [48-55].

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

Вторая глава работы посвящена разработке специализированных средств обработки результатов, их анализа и визуализации. Приведено описание комплекса программ, решающего следующие задачи:

1) Хранение, добавление и изменение информации о термодинамических свойствах индивидуальных веществ и кинетических механизмов.

2) Быстрое и однозначное задание начальных данных.

3) Извлечение из расчетных данных информации о каждой функции и каждом параметре, интерполяция данных.

4)Визуализация расчетных данных.

В третьей главе настоящей работы рассматриваются процессы термически неравновесной химической кинетики в ударных волнах при отсутствии равновесия между поступательными и колебательными степенями свободы молекул. Представлены уравнения, описывающие термическую неравновесность. Проведено математическое моделирование неравновесных процессов за фронтом УВ как в стационарном, так и в нестационарном режимах. Приводятся результаты оценки влияния колебательной релаксации на структуру УВ, а также проводится оценка влияния электронного возбуждения молекул на структуру УВ в кислороде. Приводятся результаты оценки влияния электронного возбуждения на распределение параметров течения за стационарной ударной волной.

Четвертая глава посвящена вопросам численного моделирования газовой детонации. Рассмотрена задача об инициировании детонации за отраженной ударной волной в ударной трубе в горючей смеси водорода с кислородом. Проведено численное моделирование серии экспериментов по измерению задержек воспламенения сильно разбавленных аргоном водородо-кислородных смесей (3-5%) в ударных трубах за падающими ударными волнами. Моделирование представлено в стационарной и нестационарной постановках. Проведено сравнение полученных результатов с экспериментальными данными Павлова В.А. (НИИ Механики МГУ).

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Северина, Наталья Сергеевна

ВЫВОДЫ К ГЛАВЕ 4

Разработаны и оттестированы методика математического моделирования и комплекс программ для расчета детонационных процессов в многокомпонентных газовых смесях в ударных трубах.

Исследовано влияние числа Маха падающей УВ на параметры течения за отраженной УВ в горючей смеси водорода с кислородом. Расчетным путем получено, что чем больше число Маха падающей УВ, тем быстрее происходит воспламенение горючей смеси на торце УТ и, как следствие, более быстрое образование в горючей смеси висячей ударной волны, оказывающей существенное влияние на развитие детонации.

При моделировании экспериментов в ударной трубе исследована зависимость времени задержки воспламенения горючей смеси от начального давления и степени разбавления аргоном. Сравнение измеренных в ходе эксперимента и рассчитанных численно времен задержек воспламенения сильно разбавленных аргоном водородо-кислородных смесей (3-5%) в ударных трубах за стационарными падающими ударными волнами свидетельствует об их удовлетворительном согласии.

Выявлены режимы существования и отсутствия стационарных детонационных волн в сильно разбавленных аргоном горючих смесях водорода с кислородом.

Расчетным путем получена детальная картина процесса инициирования детонационных волн за падающей УВ в ударных трубах. Картина течения включает: воспламенение горючей смеси у контактного разрыва, образование волны горения, взаимодействие волны горения и головной ударной волны и, как следствие, образование нестационарной пересжатой детонационной волны, переход детонационной волны к стационарному режиму распространения.

ЗАКЛЮЧЕНИЕ

В диссертационной работе получены следующие результаты: Разработана физико-математическая модель, вычислительные алгоритмы и комплекс программ для моделирования квазиодномерных нестационарных течений многокомпонентного невязкого нетеплопроводного реагирующего газа.

Создан вариант сеточно-характеристического метода для расчета нестационарных реагирующих течений произвольных многокомпонентных газовых смесей, с выделением произвольного числа сильных и слабых разрывов параметров течения. Особенностью реализации является то, что все разрывы (ударные волны, контактные разрывы, характеристики, ограничивающие веера волн разряжения) являются сеточными линиями, при этом отслеживаются все возможные взаимодействия между ними, также расчетными узлами являются траектории частиц газа и характеристики.

Разработаны оригинальные формат выходных данных и система визуализации результатов численного моделирования одномерных нестационарных течений многокомпонентного реагирующего газа с произвольным числом взаимодействующих линий разрыва. Особенностью системы визуализации является то, что она ориентирована на работу с данными, соответствующими подвижным сеточным линиям. Причем сеточные линии могут пересекаться между собой, число их может увеличиваться и уменьшаться в процессе расчетов.

Комплекс программ визуализации полностью независим от расчетного модуля, его можно использовать для хранения и визуализации результатов численного решения произвольной физической задачи, описываемой двумерной гиперболической системой дифференциальных уравнений в частных производных. Чтобы использовать программу визуализации разработчику необходимо предусмотреть сохранение данных в файлы описанного выше формата, задать конфигурационные файлы, описывающие структуру расчетных данных, что позволяет программе анализа распознавать их и "настроиться" на произвольное физическое содержание.

Разработаны вычислительные алгоритмы и комплекс программ, позволяющие моделировать тонкую структуру ударных волн в кислороде в ударных трубах в стационарной и нестационарной постановках с учетом колебательной и химической неравновесности, а также электронного возбуждения.

Расчетным путем обнаружено, что при распространении сильных ударных волн в молекулярном газе наблюдается значительное нарушение термического равновесия по колебательным степеням свободы молекул, как в стационарных, так и в нестационарных режимах течения, которое происходит непосредственно за фронтом ударной волны. Это способствует неравновесному возбуждению электронных состояний атомов и молекул и образованию значительных радиационных потоков в ударно-нагретых газах в дополнение к поступательной неравновесности во фронте ударной волны.

Расчеты параметров течения за стационарной ударной волной в молекулярном кислороде проведены с учетом процессов колебательной релаксации и кинетических процессов образования и тушения двух нижних синглетных электронных состояний молекулярного кислорода. Учет ЭВ влияет на результаты расчета состава смеси: происходит наработка кислорода в состояниях а и Ь, и удлиняется релаксационная зона за фронтом ударной волны.

Разработаны и оттестированы методика математического моделирования и комплекс программ для расчета детонационных процессов в многокомпонентных газовых смесях в ударных трубах.

Исследовано влияние числа Маха падающей УВ на параметры течения за отраженной УВ в горючей смеси водорода с кислородом. Расчетным путем получено, что чем больше число Маха падающей УВ, тем быстрее происходит воспламенение горючей смеси на торце УТ и, как следствие, более быстрое образование в горючей смеси висячей ударной волны, оказывающей существенное влияние на развитие детонации.

При моделировании экспериментов в ударной трубе исследована зависимость времени задержки воспламенения горючей смеси от начального давления и степени разбавления аргоном. Сравнение измеренных в ходе эксперимента и рассчитанных численно времен задержек воспламенения сильно разбавленных аргоном водородо-кислородных смесей (3-5%) в ударных трубах за стационарными падающими ударными волнами свидетельствует об их удовлетворительном согласии.

Выявлены режимы существования и отсутствия стационарных детонационных волн в сильно разбавленных аргоном горючих смесях водорода с кислородом.

Расчетным путем получена детальная картина процесса инициирования детонационных волн за падающей УВ в ударных трубах. Картина течения включает: воспламенение горючей смеси у контактного разрыва, образование волны горения, взаимодействие волны горения и головной ударной волны и, как следствие, образование нестационарной пересжатой детонационной волны, переход детонационной волны к стационарному режиму распространения.

Список литературы диссертационного исследования кандидат физико-математических наук Северина, Наталья Сергеевна, 2009 год

1. Гипинский С.М., Теленин Г.Ф., Тиняков Г.П. Метод расчета сверхзвукового обтекания затупленных тел с отошедшей ударной волной // Изв.АН СССР. Мех. и маш, 1964, № 4, С.9-28.

2. Белоцерковский О.М. и др. Обтекание затупленных тел сверхзвуковым потоком газа. Теоретические и экспериментальные исследования. М.: ВЦ АН СССР, 1966, С.340.

3. Любимов А.Н., Русанов В.В. Течения газа около тупых тел. 4.1. М.: Наука, 1970, С.282.

4. Бутаков В.А. Метод расчета течений идеального газа с выделением ударных волн. Примеры течений около тел с протоком // Исследование газодинамики и теплообмена сложных течений однородных и многофазных сред. М: Изд-во МГУ, 1990, С.20-25.

5. Моретти Дж. К вопросу о выделении скачка //Численное решение задач гидромеханики. М.: Мир, 1977, С.55-63.

6. Белоцерковский О.М., Грудницкий В.Г., Рыгалин В.Н. Выделение разрывов при расчете одномерных нестационарных течений газа // ДАН СССР, 1983, Т.272, № 1, С.49-52.

7. О.А. Азарова, В.В. Власов, ВТ. Грудницкий, В.Н. Рыгалин. Расчёты одномерных газодинамических течений с выделением ударных волн и контактных разрывов. // Сб. Алгоритмы для численного исследования разрывных течений, 1993, с.56-80.

8. О.А. Азарова, В.В. Власов, В.Г. Грудницкий, Н.А. Попов, В.Н. Рыгалин. Разностная схема на минимальном шаблоне и её применение в алгоритмах выделения разрывов. // Сб. Алгоритмы для численного исследования разрывных течений, 1993, с.9-56.

9. Крайко А.Н., Макаров В.Е., Тгшляева Н.И К численному построению ударных волн //ЖВМ и МФ, 1980, Т.20, № 3, С.716-723.

10. Капитанец П.И., Лунев В.В. Ударно-характеристический метод решения разрывных газодинамических задач //ЖВМ МФ, 1988, Т. 28, № 6, С.945-950.

11. Кенцер Ч. Дискретизация граничных условий на движущихся разрывах //Численные методы в механике жидкости. М.: Мир, 1973, С.62-72.

12. Грудницкий В,Г. Прямой обобщенно-характеристический метод для расчета разрывных решений законов сохранения газовой динамики. // Матем. Моделирование, 2004, т. 16., №1, с.90-96.

13. Стулов В.Я. Расчет невязких течений газа с внутренними разрывами. // Матем. Моделирование, 1992, т.4., №12, с.97-105.

14. Кутлер П., Уорминг Р., Ломаке Ю. Расчет обтекания транспортного космического корабля с использованием нецентральных разностных схем //Ракетная техника и космонавтика, 1973, T.l 1, №2, С.86-97.

15. В.Г. Грудницкий. Обобщённые характеристики для систем уравнений Эйлера и их применение к конструированию численных схем. // Матем. Моделирование, 1992, т.4, №12, с.45-48.

16. Котельников В.А., Волков В.А., Гидаспов В.Ю., Котельников М.В., Хохлов

17. A.В. Математическое моделирование обтекания тел слабоионизованной столкновительной плазмой. // М.: Изд-во МАИ, 2007. 144 е.: ил.

18. Волков В.А., Гидаспов В.Ю., Пиру мое У.Г., Стрельцов В.Ю. Численное моделирование течений реагирующих газокапельных и газовых смесей в экспериментах по воспламенению метанолаУ/ ТВТ, 1998, N 3.

19. В.Г. Грудницкий. Достаточное условие устойчивости многомерного расчёта нестационарных разрывных решений уравнений Эйлера. // Матем. Моделирование, 2000, т. 12, №1, с.65-77.

20. В.Ю.Гидаспов, И.Э.Иванов, ИА.Крюков, И.М.Набоко, В.А.Петухов,

21. B.Ю.Стрельцов. Исследование распространения волн горения и детонации вкумулирующем объеме// Математическое моделирование, 2004, т. 16, № 6, с. 118-122.

22. В.Ю. Гидаспов. Вычислительный алгоритм решения задачи о распаде произвольного разрыва в равновесно-реагирующем газе// Математическое моделирование, 2006, т. 18, н. 8, стр. 64-76.

23. Численные методы в газовой динамике. Под редакцией Павлова Б.М., Рослякова Г.С., Чудова JI.A. М.: Изд-во МГУ, 1965.-304 с.

24. Волков В.Ф. Расчет взаимодействия двух пластин при обтекании сверхзвуковым потоком невязкого газа //Аэродинамическая интерференция при обтекании пространственных тел. Новосибирск: ИТПМ СО АН СССР, 1980, СЗ-19.

25. Баженова Т. В., Гвоздева Л.Г. Нестационарные взаимодействия ударных волн. М., Наука, 1977, 274 с.

26. Демьянов Ю. А., Секриеру Г. В., Игогиин А. И., Киреев В. Т., Пинский B.JI. Одномерные нестационарные течения реального газа. Кишинев, Штиница, 1980, 188 с.

27. Денисов Е.Т. Кинетика гомогенных химических реакций. М., Высшая школа, 1978, 367 с.

28. Денисов Ю. Н. Газодинамика детонационных структур. М., Машиностроение, 1989, 176 с.

29. Детонация конденсированных и газовых систем. М., Наука, 986, 318 с.

30. Нетлетон М. Детонация в газах. Пер. с англ. М.: Мир, 1989, 280

31. Э. Оран, Дж. Борис. Численное моделирование реагирующих потоков. М.: Мир, 1990, 661 с.

32. Основы газовой динамики: пер. с англ. / ред. Г. Эммонс. М., ИИЛ, 1963, 702 с.

33. Пирумов У. Г. Обратная задача теории сопла М: Машиностроение, 1988, 240 с.

34. Пирумов У. Г, Росляков Г. С. Течение газа в соплах. -М, МГУ, 1978, 352 с.

35. Пирумов У. Г., Росляков Г. С. Газовая динамика сопел. -М.: Наука, Гл. ред. физ.-мат. лит., 1990. 368 с.

36. К.П.Станюкович. Неустановившиеся движения сплошной среды. Москва, Наука, 1971, 856 с.

37. Черный Г.Г. Газовая динамика. М.: Наука, 1988, 424 с.

38. Современные численные методы решения обыкновенных дифференциальных уравнений. / Пер, с англ. //Под ред. Дж. Холл и Дж. Уатт. -М., Мир, 1979.

39. Волков В. А., Мусин В. Р., Прохоров М. Б. Численный метод решения системы конечно-разностных уравнений химической кинетики. // Тематический сборник научных трудов. МАИ. 1988.

40. Годунов С. К., Забродин А. В., Иванов М. Я., Крайко А. П., Прокопов Г. П. Численное решение многомерных задач газовой динамики. Под ред. С. К. Годунова. М.: Наука, 1976,400 с.

41. Oran E.S., Boris J.P. Weak and Strong Ignition: II, Sensitivity of the Hydrogen-Oxyden System, Comb. Flame, 48, 145-161, 1982.

42. Oran E.S., Young T.R., Boris J.P., Coxen A. Weak and Strong Ignition: I, Numerical Simulation of Shock Tube Experiments, Comb. Flame, 48, 135-148, 1982.

43. Harten A., Osher S. Uniformly high-order accurate nonoscillatory schemes. // SIAM Journal Numer. Anal., 1987, № 2, pp. 279-309.

44. Harten A. Preliminary results on the extension of ENO schemes to two-dimensional problems, in Nonlinear Hyperbolic Problems (St. Etienne, 1986), C. Carasso, P.-A. Raviart and D. Serre (Eds.), Lect. Notes in Math., pp. 23-40,Springer, Berlin, 1987.

45. Liu X.-D., Osher S., Chan T. Weighted essentially non-oscillatory schemes, // Journal of Computation Physics, 1994, №1, pp.200-212.

46. Bauer P. A., Dabora E. K, Manson N. Chronology of early research on detonation wave // Dynamics of Detonations and Explosions: Detonations / A. L.

47. Gelfand В. E., Frolov S. M., Nettleton M. A. Gaseous detonations — a selective review // Prog. Energy Combust. Sci. 1991.V. 17. P. 327-371.

48. Васильев А. А., Митрофанов В. В., Топчиян М. Е. Детонационные волны в газах // Физика горения и взрыва. 1987. Т. 23, № 5.С. 109-131.

49. Гусева Н.С. Сравнение различных методов расчета АД. // Тезисы докладов XI Международной конференции по вычислительной механике и современным прикладным программным системам, Москва-Истра, 2-6 июля 2001 г. — М.: Изд-во МАИ, 2001, с. 165-166.

50. Гидаспов В.Ю., Пирумов У.Г., Северина Н.С. Математическое моделирование квазиодномерных нестационарных течений реагирующего газа с произвольным числом взаимодействующих разрывов. // Вестник МАИ. М.: Изд-во МАИ, 2008. -т.15, №5, с. 83-94

51. Гидаспов В.Ю., Лосев С.А., Северина Н.С., Ярыгина В.Н. Влияние электронного возбуждения на структуру ударной волны в кислороде. // Вестник МАИ. М.:Изд-во МАИ, 2009. т.16, №2, с. 93-100

52. Гидаспов В.Ю., Лосев С. А., Северина Н.С, Роль термической неравновесности в моделировании химической кинетики на примере диссоциации молекулярного кислорода за фронтом ударной волны. // Математическое моделирование. 2009, т.21, №9, с. 3-15.

53. Физико-химические процессы в газовой динамике. Справочник. Том первый. Динамика физико-химических процессов в газе и плазме (второе издание). М.: Изд-во Научный мир, 2007, 400 с.

54. Лосев С.А., Ковач Э.А, Погосбекян М.Ю., Сергиевская А.Л. Моделирование физико-химических процессов в сильных ударных волнах. // Электронный журнал «Физико-химическая кинетика в газовой динамике» http://www.chemphys.edu.ru Том 1,2003-12-24-001.

55. Моделирование колебательной релаксации двухатомных молекул при высоких температурах. Отчет НИИМеханики МГУ № 4736. стр.46-54, 2004.

56. Park С. Nonequilibrium Hypersonic Aerothermodynamics. Wiley. NY., 1990.

57. Шаталов О.П. Исследование колебательного возбуждения, дезактивации и диссоциации молекул кислорода в сильных ударных волнах и охлаждающихся потоках. Дисс. канд. физ.-мат. наук. М. 1973.

58. Шаталов О.П. О диссоциации молекулярного кислорода в отсутствие колебательного равновесия. Физика горения и взрыва, 1973, № 5, стр.699-703.

59. Куликов С.В., Смирнов A.JI., Терновая О.Н. Влияние поступательной неравновесности смеси 02 и инертных газов на одноступенчатую диссоциацию С>2 во фронте ударной волны. // Химическая физика. 2000, том 19, № 12, стр.53-61.

60. Горелов В.А., Комаров В.Н., Кузнецов М.М., Юмашев В.Л. О влиянии поступательной неравновесности на скорость молекулярной диссоциации в гиперзвуковой ударной волне. // Журнал Прикладной Механики и Технической Физики. 2001, том 42, № 2, стр.42-51.

61. Великодный В.Ю., Битюрын В.А. Влияние эффектов поступательной неравновесности на кинетику физико-химических превращений во фронте ударной волны. // Письма в Журнал Технической Физики. 1996, том 22, вып.4, стр.39-45.

62. Нагнибеда Е.А., Kycmoea Е.В. Кинетическая теория процессов переноса и релаксации в потоках неравновесных реагирующих газов. // СПб.: Изд-во СПб гос. ун-та, 2003, 272 с.

63. Ступоченко Е.В., Лосев С.А., Осипов А.И. Релаксационные процессы в ударных волнах. // М.: Наука, 1965, 474 с.

64. А.В.Виноградов, В.А Волков, В.Ю. Гидаспов. Метод расчета нестационарных реагирующих течений с явным выделением поверхностей сильных и слабых разрывов. // Математическое моделирование, том 8, N 3, 1996, с. 79-90.

65. Aliat A., Kustova E.V., Chikhaoui А. II Chemical Physics,2005, Vol.314, pp.3 7-47.

66. A.M. Старик, H.C. Титова. //Кинетика и катализ. 20036 Т.43, №1, с.1-12.

67. Физико-химические процессы в газовой динамике. Справочник. Том второй. Физико-химическая кинетика и термодинамика. М. Научно-изд. Центр механики. 2002.

68. Л.В. Гурвич, И.В. Вещ, В.А. Медведев и др. Термодинамические свойства индивидуальных веществ: Справочное издание в 4-х т. //. М.: Наука, 1978.

69. Павлов В.А., http://www.chemphvs.edu.ги/, 2009 г, том 8.

70. Ибрагимова Л.Б., Смехов Г.Д., Шаталов О.П., http://www. chemphvs. edu. ги/, 2009 г, том 8.

71. Roe P.L. Approximate Riemann Solvers, Parameter Vectors and Difference Schemes, J. Сотр. Phys.,43,357-372,1981.

72. Roe P. Approximate Riemann Solvers, Parameter Vectors and Different Schemes//J. Comput. Phys., 1981, v.43,N2, 357-372.

73. Schuster A. Philos. Trans. Roy. Soc. London A, 1893, vol. 152, p. 21.

74. Sod G. A' Survey of Several Finite Difference Methods Systems of Nonlinear Hyperbolic Conservation Zaws.- J. Comput. Phis,v. 27, 1978, pp. 1-31.

75. Van Leer B. Lect. Notes Phys. 1982. v. 170. pp. 507-512 .

76. Иванов И.Э., Крюков HA. Квазимонотонный метод повышенного порядка точности для расчета внутренних и струйных течений невязкого газа. // Математическое моделирование, т.8, №6, 1996, с. 47-55.

77. Дремин А.Н., Савров С.Д., Трофимов В.А., Шведов К. К. Детонационные волны в конденсированных средах / М.: Наука, 1970.

78. Воронин Д. В., Ждан С. А. Об одномерной неустойчивости детонационных волн в распылах // Физика горения и взрыва. 1986. Т. 22,№ 4. С. 92-98.

79. Зельдович Я. Б., КомпанеецА. С. Теория детонации. М.: Гостехиздат, 1955.

80. Николаев Ю.А., Топчиян М.Е. Расчет равновесных течений в детонационных волнах в газах // Физика горения и взрыва. 1977.Т. 13, № 3. С. 393-404.

81. Левин В. А., Марков В. В. Возникновение детонации при концентрированном подводе энергии // Физика горения и взрыва. 1975.Т. 11, № 4. С. 623-633.

82. Щетинков Е. С. Физика горения газов. М.:Наука, 1965.

83. Николаев Ю. А. Обобщенная модель кинетики химических реакций в водородокислородных газовых смесях // Физика горения и взрыва. 1994. Т. 30, № 1. С. 66-72.

84. Азатян В. В., Бакланов Д. И., Гвоздева Л. Г. и др. Ингибирование развившейся детонации водородно-воздушных смесей //Докл. АН. 2001. Т. 376, № 1.С. 55-58. '

85. Николаев Ю. А. Приближенное моделирование, модель кинетики и калорическое уравнение состояния химически реагирующих газовых смесей при высоких температурах //Физика горения и взрыва. 2001. Т. 37, № 1.С. 6-15.

86. Николаев Ю. А., Фомин П. А. О расчете равновесных течений химически реагирующих газов // Физика горения и взрыва. 1982.Т. 18, № 1. С. 66-72.

87. Фомин П; А., Трогрок А. В. Приближенный расчет изоэнтропы химически равновесного газа // Физика горения и взрыва. 1995.Т. 31, № 4. С. 59-62.

88. Николаев Ю. А. Распространение плоских пересжатых детонационных волн в газах//Физика горения и взрыва. 1995. Т. 31, № 5.С. 101-113.

89. Химия горения: Пер. с англ / Под ред.У. Гардинера-мл. М.: Мир, 1988.

90. Николаев Ю.А., ЗакД.В. Квазиодномерная модель самоподдерживающейся многофронтовой газовой детонации с учетом потерь и турбулентности // Физика горения и взрыва.1989. Т. 25, № 2. С. 103-112.

91. Николаев Ю. А., Гапонов О. 77. О пределах детонации в газах // Физика горения ивзрыва. 1995. Т. 31, № 3. С. 139-145.

92. Бедарев И. А., Федоров А. В. Сравнительный анализ трех математических моделей воспламенения водорода // Физика горения и взрыва. 2006, Т. 42, № 1, с. 26-33.

93. Гриб А. А., 1941 — цитируется по: К. П. Станюкович. Неустановившиеся движения сплошной среды. М.: Гос. изд-во техн.-теор. лит.,1955. С. 364.

94. Ждан С. А. Расчет взрыва газовых смесей с учетом сдвига химического равновесияпродуктов // Физика горения и взрыва. 1983.Т. 19, № 1. С. 131-135.

95. Астапов Н. С, Николаев Ю. А., Улъяницкий В. Ю. О параметрах детонации водородно-кислородных и водородно-воздушных смесей при высокой начальной плотности // Физика горения и взрыва. 1984.Т. 20, № 1. С. 98-105.

96. Николаев Ю. А., Васильев А. А., Улъяницкий В.Ю. Газовая детонация и ее применение в технике и технологиях (обзор). // Физика горения и взрыва, 2003, т. 39, №4.

97. Котельников В.А., Волков В.А., Гидаспов В.Ю., Котельников М.В., Хохлов А.В. Математическое моделирование обтекания тел слабоионизованной столкновительной плазмой. // М.: Изд-во МАИ, 2007. 144 е.: ил.

98. Адрианов А.Л., У сков В.Н., Старых A.JI. Интерференция стационарных газодинамических разрывов. Новосибирск: ВО «Наука». Сибирская издательская фирма, 1995. - 180 с.

99. Гладилин A.M., Григоров А.И., Сагидуллин Г.Г. Детонационные процессы вдвухфазных средах. М.: Недра, 1991. - 144 е.: ил.

100. Нетлетон М. Детонация в газах. Пер. с англ. М.:Мир, 1989. - 280 с.:илл.

101. Shu, С. —W. & Osher, S. 1989 Efficient implementation of essentially non-oscillatory shock-capturing schemes II. J. Сотр. Phys. 83, 32-78.

102. Dohi Т., Hayashi A.K., Ogawa S. Numerical simulations of layered detonation using a detailed hydrogen chemistry /21 si International Symposium on Shock Waves, Great Keppel Island, Australia, July 20-25, 1997, Paper 2660.