В настоящее время продолжает развиваться космическая и авиационная промышленность. При этом она нуждается в разработке все более совершенных и экономически выгодных технических устройствах. Важнейшим из таких устройств является реактивный двигатель, создающий необходимую для движения силу тяги посредством преобразования потенциальной энергии топлива в кинетическую энергию реактивной струи рабочего тела. Любой реактивный двигатель должен иметь, по крайней мере, две составные части: камеру сгорания и реактивное сопло. Камера сгорания («химический реактор») в нем происходит освобождение химической энергии топлива и её преобразование в тепловую энергию газов. Реактивное сопло («газовый туннель») в котором тепловая энергия газов переходит в их кинетическую энергию, когда из сопла газы вытекают наружу с большой скоростью, создавая тем самым реактивную тягу.
Исторически первые применения сопел связаны с паровыми турбинами, в которых они были использованы в 1890 г. шведским инженером Густафом де Лавалем [1]. В 1925 г. сопло впервые было использовано в ракетном двигателе. Сопло Лаваля - это устройство, которое служит для ускорения рабочего тела до скоростей, превышающих скорость звука. Рабочими телами в соплах Лаваля являются продукты сгорания жидких или твердых топлив. Сопла широко используется в технике: паровых и газовых турбинах, ракетных и воздушно¬реактивных двигателях, газодинамических лазерах, магнитно -газодинамических установках, аэродинамических установках и газодинамических стендах, при создании молекулярных пучков, химической технологии, струйных аппаратах, расходомерах, при дутьевых процессах и многих др.
Исследование характеристик прорабатываемых сопел проводят при помощи экспериментальных исследований и (или) математического моделирования. Проведение экспериментальных исследований связано со значительными расходами, в связи с этим предварительную проработку технических решений осуществляют на основе математического моделирования. При проектировании ракетных двигателей необходимо учитывать влияние параметров окружающей среды на энергетические характеристики сопла и истекающую струю. В общем случае данная задача носит пространственный характер.
При старте ракет-носителей процесс истечения сверхзвуковых струй из сопел ракетного двигателя сопровождается излучением акустических шумов с высокими уровнями звуковых давлений. Они создают значительные акустические нагрузки на различные части корпуса ракеты -носителя и обтекатель, под которым установлена полезная нагрузка. Проникая сквозь обшивку, акустическое излучение пагубно влияет на выводимый груз. Распространяясь в окружающей среде, акустическое излучение создает дополнительные нагрузки на стартовые сооружения [2].
Повышение требований к подобным задачам влечет за собой повышение требований к точности и качеству методов моделирования. В связи с этим разрабатываются более совершенные численные методы, внедряются методы параллельных вычислений, немалое внимание уделяется также методам обработки результатов и визуализации. В конечном виде эти разработки представляются в виде прикладного программного обеспечения, которое может распространяться под различными лицензиями, работать под различными операционными системами, распространяться в виде исходного кода [3, 4], библиотек разработки или в виде самодостаточных программ .
Актуальность исследования сверхзвуковых струйных течений обусловлена их широким применением в ракетно-космической технике. Целью настоящей работы является реализация методики расчета и исследование течения газа в соплах и струях с использованием пакета прикладных программ OpenFOAM
Реализована методика расчета течения газа в соплах и струях с использованием пакета прикладных программ OpenFOAM.
Проведено сравнение результатов расчетов с экспериментальными данными других авторов, получено хорошее совпадение. Совпадение наблюдается в положении диска Маха, скачков сжатия, отраженных скачков, «висячего» скачка уплотнения и «соплового» скачка.
Проведены параметрические исследования для одиночной сверхзвуковой струи истекающей в затопленное пространство при степенях нерасчетности n=0.57, 0.7, 1.0, 1.5, 2.0. Показаны ударно-волновые структуры струй и
распределения давления в поперечном сечении струи. Получено, что с ростом степени нерасчетности увеличивается диаметр струи.
Проведено исследование акустических характеристик сверхзвуковой струи в одной точке на удалении 0.21 метра от среза сопла. Показано, что максимальная амплитуда колебаний достигается при степени нерасчетности n=2.0 на частоте 175 Гц, а уровень максимального звукового давления составляет
1. Дорофеев А. А. Основы теории тепловых ракетных двигателей: (Общая теория ракетных двигателей): Учебник для авиа - и ракетостроительных специальностей вузов. - М.: Изд-во МГТУ им. Н. Э. Баумана, 1999. - 415 с.
2. Струйные и нестационарные течения в газовой динамике /В. Н. Глазнев, В. И. Запрягаев, В. Н. Усков и др. ; Под ред. С. А. Гапонова, А. А. Маслова; Рос. акад. наук, Сиб. отд. , Балтийский гос. техн. ун-т "Военмех" им. Д. Ф. устинова
3. OpenFOAM // Официальный сайт [Электронный ресурс] / User Guide - http://www.openfoam.com/docs/user (дата обращения 13.11.2020).
4. OpenFOAM // Официальный сайт [Электронный ресурс] / OpenFOAM C++ Documentation - http://www.openfoam.com/docs/cpp (дата обращения 13.11.2020).
5. Теория ракетных двигателей: Учебник для студентов высших технических учебных заведений / В. Е. Алемасов, А. Ф. Дрегалин, А. П. Тишин; Под ред. В. П. Глушко. - М.: Машиностроение, 1989. - 464 с.
6. Абрамович Г. Н. Прикладная газовая динамика. В 2 ч. Ч. 1: Учеб. руководство: Для втузов. - 5-е. изд., перераб. и доп. - М.: Наука. Гл. ред. физ.- мат. лит., 1991. - 600 с.
7. Смирнов Е. М. Течения вязкой жидкости и модели турбулентности: методы расчета турбулентных течений / Е. М. Смирнов, А. В. Гарбарук // Конспект лекций. - 2010.
8. Toro E.F. Restoration of the contact surface in the HLL-Riemann solver / E.F. Toro, M. Spruce, W. Speares // Shock Waves. - 1994. - V. 4. - P. 25-34
9. Белов И.А. Моделирование турбулентных течений : учеб. пособие / И.А. Белов, С.А. Исаев. - СПб.: Балт. гос. тех. ун-т, 2001. - 108 с.
10. Menter F.R. Two-equation eddy-viscosity turbulence models for engineering applications / F.R. Menter // AIAA Journal. - 1994. - V. 32, № 8. - P. 1598-1605
11. Годунов С.К. Разностный метод численного расчета разрывных решений уравнений гидродинамики / С.К. Годунов // Матем. сб.. - 1959. - Т. 47, № 3. - С. 271-306
12. Salome // Официальный сайт [Электронный ресурс] / User Section - URL http://www.salome-platform.org/user-section (дата обращения 10.11.2020).
13. Математическое моделирование акустических излучений при взаимодействии сверхзвуковой струи с плоской преградой: / С.А. Николин, Г.И. Сокол. Вестник Днепровского национального университета. 2018. Т. 26, №4, С. 73-80.
14. Запрягаев В.И. расчетно-экспериментальное исследование
газодинамической структуры структуры сверхзвуковой перерасширенной струи / И. Н. Кавун, С. Г. Кундасев // ПМТФ. 2010. № 4, с. 58—64.
15. Годунов С.К. Численное решение многомерных задач газовой динамики / С.К. Годунов, А.В. Забродин, М.Я. Иванов, А.Н. Крайко, Г.П. Прокопов. - М.: Наука, 1976. - 400 с.