Численное профилирование плоских и осесимметричных сверхзвуковых сопел и каналов для моделирования газодинамических течений тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

ВВЕДЕНИЕ.

ГЛАВА I. ПРОФИЛИРОВАНИЕ СВЕРХЗВУКОВЫХ! КАНАЛОВ ДЛЯ ТЕЧЕНИЙ ИДЕАЛЬНОГО ГАЗА НА ОСНОВЕ РЕШЕНИЯ СМЕШАННЫХ КРАЕВЫХ ЗАДАЧ.

§ I. Постановка обратной задачи профилирования и схема ее решения.

§ 2. Построение модификации сеточно-характеристического метода по слоям ^ = const . gj

§ 3. Результаты расчетов.

§ Формулировка двухграничных смешанных краевых задач профилирования.

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

ГЛАВА П. ПРОФИЛИРОВАНИЕ СВЕРХЗВУКОВЫХ КАНАЛОВ ДЛЯ НЕАДИАБАТИЧЕСКИХ ТЕЧЕНИЙ И С УЧЕТОМ НЕРАВНОВЕСНОГО ПРОТЕКАНИЯ ХИМИЧЕСКИХ РЕАКЦИЙ В ПОТОКЕ.

§ I. Профилирование каналов с учетом энерговыделения в потоке.

§ 2. Метод численного профилирования для течений с неравновесными химическими реакциями. 7д

§ 3. Результаты расчетов для течений реагирующей смеси газов, содержащих атомы С, Н} О, л/. qq

ВЫВОДЫ. юз

ГЛАВА Ш. ПРОФИЛИРОВАНИЕ СВЕРХЗВУКОВЫХ КАНАЛОВ С ЗАДАННЫМИ

НЕИЗЭНТРО ПИЧЕСКИМИ ПАРАМЕТРАМИ НА ВЫХОДЕ.

§ I. Постановка и метод решения задачи профилирования каналов с гладкими параметрами на выходе. Ю

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

§ 3. Численное профилирование каналов с разрывным распределением энтропии на выходе.

ВЫВОДЫ.

ч ется задача о построении (профилировании) плоского иди осесимметричного сверхзвукового канала с некоторыми априори заданными свойствами. Решение подобной задачи необходимо при построении рабочих трактов экспериментальных установок, моделирующих реальные течения, при разработке сопел реактивных двигателей, аэродинамических труб, сопел газодинамических лаверов (ГДЛ), переходных каналов различного класса, воздухозаборников и целого ряда других газодинамических устройств*

Уравнения газовой динамики являются нелинейными уравнения* ми в частных производных /I/» аналитические решения их возможны лишь в некоторых частных случаях* Кроме того, в практических задачах необходимо учитывать реальные процессы в газе: химические реакции, излучение, конденсацию и т.д., что значительно усложняет как постановку задачи, так и создание эффективны!; методов ее решения* В связи с этим основным методом профилирования является проведение численных расчетов на ЭВМ*

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

В зависимости от расположения R- относительно потока задачу профилирования южно подразделить на два класса:

I) Я = Г - является сечением потока, обыкновенно его полагают выходным;

2) R. = L - линия тока, в качестве которой в большинстве случаев выбирается либо ось (плоскость) симметрии, либо образующая центрального тела.

Первый алгоритм профилирования осесимметричных сверхзвуковых сопел с равномерными изэнтропическими параметрами в выходном сечении Г % основанный на использовании члененного метода характеристик, был предложен в /2/. Так как требуемый поток даже в классе плоских и осесимметричных течений можно получить бесчисленным количеством способов, то необходим дополнительный критерий выбора. В /2/ за основу был взят принцип быстрейшего разгона потока. Это приводит к появлению точек излома контура, если не наложено ограничение на кривизну поверхности канала. Звуковая линия предполагается прямолинейной, что дает возможность отдельно профилировать дозвуковую и сверхзвуковую части сопла, так как при этом линия перехода является одновременно характеристикой первого и второго семейств /3/. Профилирование начинается с расчета обтекания заданного в виде угловой точки закритического разгонного участка сопла. Для расчета характеристики второго семейства центрированной волны, близкой к звуковой линии, используется представление решения в виде ряда по характеристической переменной* Находится характеристика второго семейства, ограничивающая область свободного расширения потока, и дополняется справа равномерной характеристикой первого семейства* Последующее решение задачи Гурса с данными на этих характеристиках завершает определение контура сверхзвуковой части сопла.

На основе этого метода в классе сопел с равномерным выходом были спрофилированы осесимметричные без центрального тела /4-6/, кольцевые /7-9/, плоские симметричные /Ю/ и плоские несимметричные /11-13/ каналы* Обобщением метода /2/ явилась разработка алгоритмов, позволяющих профилировать сопла с заданными неравномерными изэнтропическими параметрами на выходе /14-20/. В отличие от /2J здесь замыкающая характеристика первого семейства не является равномерной, а определяется из решения задачи Коши с неравномерными данными в выходном сечении. В работах /14, 16/ данный метод был применен для построения сопел модельных газодинамических установок. В /15, 17/ подробно рассмотрен вопрос о профилировании кольцевых сопел, в частности, задача "адаптации11, т.е. согласования входного до- и трансзвукового участка с искомым сверхзвуковым. В /18-20/ задача профилирования решалась сеточно-характеристическим методом при больших неоднородно стях сверхзвукового потока в выходном сечении Г •

Обобщением рассмотренных выше постановок для профилирования каналов с заданными изэнтропическими параметрами на Г является предложенная в /17/ схема построения плоского сверхзвукового канала с заданным неиззнтропическим потоком на выходе при изэнтропическом течении на входе» Требуемое распределение энтропии 5 , которое предполагается достаточно гладким, создается специально выстраиваемой ударной волной, расположенной на входе в искомый канал* Однако, в предложенной схеме возникает проблема согласования входного потока с заданным выходным, решение которой в ооесимметричном случае достаточно сложно*

Эффективным методом профилирования с граничными условиями на известной линии тока L , подучившим широкое распространение, является численный метод У.Г.Пирумова /21/, сводившийся к решению задачи КОши, которая в дозвуковой области интерпретируется как условно корректная* Метод позволяет по единому алгоритму строить до-, транс- и сверхзвуковые участки сопел. На основе его получены сопла с равномерным выходом и заданным распределением давления (скорости) на оси (образующей центрального тела), и » течения в которых не содержат скачков уплотнения. Метод достаточно просто обобщается на течения с неравновесными физивю-хи-мическими превращениями, на пространственные течения в каналах сложной формы /22/. Однако, он накладывает жесткие ограничения на классы профилируемых каналов и на характер реализуемых в них течений. Так, некорректность задачи Коши в дозвуковой» эллиптической области обусловливает необходимость задания на L аналитических функций, специального выбора узлов разностной сетки и приводит к большим вычислительным трудностям при построении достаточно коротких сопел.

Детальный сравнительный анализ методов /2/ и /21/ профилирования сверхзвуковых сопел с равномерным потоком на выходе проведен в работе /23/.

Отметим, что постановка задачи профилирования с граничными условиями на априори известной линии тока L является не единственно возможной. Так, в ряде работ распределение некоторых параметров задается на неизвестной образующей искомого канала. Например, такой подход использован в /24, 25/ для профилирования дозвуковых сопел, где ревается краевая задача для уравнения Чаплыгина в плоскости годографа. С помощью этого метода удается построить короткие безударные и безотрывные сопла. Однако, этот метод не допускает простого обобщения на случай непотенциальных течений газа. В работе /26/ предложен метод профилирования дозвуковых и трансзвуковых каналов, течение в которых может отличаться от потенциального, а распределение скорости (давления) может задаваться как на оси, так и на искомой стенке канала. Метод основан на расщеплении исходной системы газодинамических уравнений на подсистемы более низкого порядка* Основная идея расщепления связана с выделением эллиптического или смешанного уравнения второго порядка, для численного решения которого используется метод блочной верхней релаксации.

Для профилирования плоских и осесимметричных каналов с заданным распределением давления вдоль искомых стенок для до- и трансзвуковых течений в /27/ применен метод установления. Стенки канала на участках, где задано распределение давления* предполагаются непроницаемыми и гибкими. В начальный момент их форма определяется в соответствии с приближенной оценкой (например, по одномерной теории) и в процессе счета изменяется до достижения стационарного состояния. Метод достаточно просто обобщается на трехмерный случай, однако использование метода установления требует значительных затрат времени ЭВМ.

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

Однако в этих работах профилирование ограничивалось рассмотрением каналов, обеспечивающих либо достаточно гладкие изэнт-ропические параметры на Г , либо аналитическое распределение давления вдоль L при равномерных параметрах на выходе. А предложенная в /17/ схема профилирования каналов с заданными неиз-энтропическиии параметрами на выходе является довольно сложной, в связи с этим в этой работе отсутствуют численные решения даже в классе плоских каналов. Разработанные обобщения данных методов профилирования на течение невязкого нетеплопроводного rasa с химической релаксацией также рассматривали течения с гладкими распределениями параметров на L и Г . в настоящее время отсутствуют работы» посвященные решению обратных задач профилирования каналов, реализующих заданные разрывные параметры как на L , так и на Г с учетом реальных свойств потока: неравновесного протекания химических реакций) наличия высокознтропийннх слоев, ударных волн и источников энергии» Решению этих задач посвящена данная диссертационная работа» целью настоящей работы является:

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

- разработка численного алгоритма сеточно-характеристиче-ского метода по сдоям, образованным линиями тока, ориентированного на решение обратных задач профилирования широкого класса каналов с учетом отмеченных выше реальных свойств потока;

- разработка метода профилирования каналов с заданными параметрами на выходе на основе решения смешанных краевых задач с граничными условиями на заданной характеристике и поперечно-ориентированной относительно потока линии;

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

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

- позволяет достаточно просто учитывать физические особенности течения, определять положение ударных волн и других разрывов; 1

- обеспечивает регулярность разностной сетки;

- позволяет рассчитывать области с разрывами небольшой интенсивности сквозным образом;

- является устойчивым.

Диссертация состоит из трех глав.

В первой главе диссертации для течений идеального газа разработана новая постановка обратной задачи профилирования сверхзвуковых каналов, сводящая ее к трехграничной смешанной краевой задаче с новыми граничными условиями в выходном сечении Г , накладываемыми на распределения давления Р или угла наклона вектора скорости W к оси ос . На граничные условия вдоль образующей центрального тела ^ и на Г не накладываются требования их непрерывности. Для численного решения задачи профилирования каналов разработана новая модификация сеточ-но-характеристического метода по слоям ^ = con.it ( f функция тока) /28-30/, в которой в отличие от известной модификации по слоям х- = costst интегрирование ведется поперек потока слоями, образованными линиями тока. Проведены исследования возможностей численной схемы, и на ее основе разработан алгоритм расчета с выделением ударных волн в неоднородном ^потоке. Впервые при решении обратных задач получены сверхзвуковые каналы, реализующие заданные разрывные параметры как на образующей центрального тела, так и в выходном сечении /30/.

В §§ 4-5 главы I исследуется вопрос о построении сверхзвуковых каналов на основе решения новых смешанных краевых задач профилирования с граничными условиями на характеристике и в сечении выхода /31/. Использование двухграничных смешанных задач позволило построить более широкий класс сверхзвуковых каналов в отличие от известного, получаемого на основе решения задачи Гурса /2, 14-17/. Так, построены каналы минимальной длины, реализующие течение от источника и заданные разрывные параметры в выходном сечении, которые обеспечиваются фокусировкой характеристик сжатия. Приводятся результаты параметрических расчетов по профилированию каналов с различными распределениями давления Р или угла г.

Во второй главе диссертации сформулирована и численно решена задача о профилировании каналов с заданным распределением параметров на известной линии тока /ив сечении выхода Г с учетом наличия в поле сверхзвукового течения невязкого совершенного газа системы объемно распределенных стоков и источников энергии. С помощью послойной модификации метода характеристик проведено численное построение сверхзвуковых каналов с различными способами задания функций теплоподвода (отвода) и показано, что наличие тепловыделения может приводить к значительной деформации образующих каналов. Во втором и третьем параграфах главы П предложено обобщение рассмотренного в главе I метода профилирования, а также численной схемы сеточно-характеристиче-екого метода по слоям ^ = const ддя течений химически неравновесной или замороженной снеси невязких нетеплопроводных совершенных газов, содержащих атомы С , Н , О , А/ /32/'. Приводятся результаты расчета сверхзвуковых каналов для неравновесных и замороженных течений при различных (в том числе и разрывных) граничных условиях, результаты сравнения с известными решениями из /22/.

В главе I для течений идеального газа разработан метод профилирования каналов с заданными неизэнтропическими параметрами в сечении выхода Г . Метод позволяет строить плоские и осесимметричные каналы с заданным распределением энтропии S (su) (иди полного давления р0 (V) ), а также числа Маха М (V) либо угла наклона i? ( у ) вектора скорости W к оси ос на Г и вдоль профилируемых стенок. Требуемое распределение S(V) обеспечивается специально выстраиваемой ударной водной, расположенной на входе в искомый канал. На основе предложенного метода построены плоские и осесимметричные воздухозаборники, преобразующие входное сверхзвуковое течение в практически однородный по Л/ ( у ) и лоток с близкой к звуковой скоростью и с заданным распределением & на выходе. Спрофилированы контуры плоских и осесимметричных каналов, обеспечивающих сжатие и разворот потока до А/ s 1,05 и г7~= 0° и различные распределения £ на Г . Построен плоский канал, моделирующий в сечении выхода Г вихревое течение у боковой поверхности плоского затупленного тела, обтекаемого гиперзвуковым потоком.

В последнем параграфе главы Ш впервые предложен метод профилирования каналов, реализующих заданное разрывное распределение энтропии 5 в выходном сечении. Разрыв 3 ( V ) обеспечивается с помощью новой конфигурации взаимодействующих разрывов: ударной волны и центрированного веера сжатия. В связи с этим проведено численное исследование конфигурации разрывов в широком диапазоне числа Маха невозмущенного набегающего потока: 1,2 ^ Но < Ю. Выявлены области существования конфигураций с отраженный веером разрежения и ударной годной, вычислена интенсивность отраженных разрывов. Исследованная конфигурация применена к расчету плоского сверхзвукового канала с центральным телом с разрывными распределениями S ( ^ ) и числа М ( ) в выходном сечении.

Основные результаты диссертационной работы опубликованы в /18-20, 30-32/. По материалам, содержащимся в диссертации, сделаны доклады: на научных семинарах под руководством акад. Г.И. Петрова (НИЙВЦ 'МГУ, 1980) и акад. О.М.Белоцерковского (МФТИ,

1983 г.), на научно-технических конференциях "Проблемы газовой динамики двигателей и силовых установок" (ЦЙАМ им.П.И.Баранова, 1981 г.) и Ленинградского механического института (Ленинград,

1984 г»), на школах-семинарах "Многомерные задачи механики сплошной среды" (ВЦ 00 АН СССР, институт прикладной математики АН СССР и институт гидродинамики 00 АН СССР, Красноярск, 1982 г.) и "Высокотемпературная газодинамика, ударные трубы и ударные волны (институт тепло- и массообмена иы.А.В.Лыкова АН БССР, институт механики МГУ, Минск, 1983 г.), на УШ Всесоюзном семинаре по комплексам программ математической физики (институт теоретической и прикладной механики 00 АН СССР, институт кибернетики АН УССР, Ташкент, 1983 г.).

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

ВЫВОДЫ

I. Для течений идеального газа разработаны постановка и метод профилирования сверхзвуковых плоских и осесимметричных каналов, реализующих заданное (в общем случае) разрывное распределение энтропии $ (или полного давления), а также распределения статического давления р (или числа Маха) или угла if наклона вектора скорости к оси X в сечении выхода и вдоль профилируемых стенок,

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

3. Для профилирования каналов с заданным разрывным распределением S на выходе предложена конфигурация взаимодействующих разрывов - ударной волны и центрированного веера сжатия. Проведено численное исследование данной конфигурации в широком диапазоне числа Мо (1»2^ Мс ^ Ю) невозмущенного набегающего потока - выявлены области существования решения с отраженными ударной волной и волной разрежения, вычислена интенсивность отраженных разрывов. Для Мо = 3 спрофилирован плоский канал, реализующий разрывное распределение S в сечении выхода.

ЗАКЛЮЧЕНИЕ

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

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

2. Для решения обратных задач профилирования разработан алгоритм сеточно-характеристического метода по слоям <р =

- сок.si { - функция тока), на основе которого интегрирование уравнений газовой динамики производится в направлении поперечном потоку. Метод используется как для расчета течений без выделения, так и с выделением ударных волн в потоке.

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

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

5. Для течений идеального газа предложен метод профилирования сверхзвуковых плоских и осесимметричных каналов, обеспечивающих заданные распределения энтропии 3 (полного давления), а также р или угла $ в выходном сечении и вдоль профилируемых стенок канала. Проведено численное профилирование каналов такого типа для различных распределений $ , f> и & на выходе: плоских и осесимметричных воздухозаборников, переходных и моделирующих каналов.

6. Впервые для течений идеального газа предложен метод профилирования каналов, реализующих заданное разрывное распределение S на выходе, который основан на использовании новой конфигурации взаимодействующих разрывов - ударной волны и центрированного веера сжатия. Проведено численное исследование данной конфигурации в диапазоне числа Маха Мо ( 1,2 ^ Но ^ Ю) невозмущенного набегающего потока. Для М0 = 3 построен плоский канал, реализующий разрывное распределение 5 на выходе.

1. Мизес Р. Математическая теория течений сжимаемой жидкости. - M.S Мир, 1961. - 588 с.

2. Кацкова О.Н., Шмыглевский Ю.Д. Осесимметричное сверхзвуковое течение свободно расширяющегося газа с плоской переходной поверхностью. В сб.: Вычислительная математика. М,, 1957, №2, с.45-89

3. Рыжов О.С. Исследование трансзвуковых течений в соплах Лаваля. М.: Изд-во Вычисл. центра АН СССР, 1965. - 238 с.

4. Кацкова О.Н., Крайно А.Н. Расчет осесимметричных изэн-тропических течений реального газа. Ж. вычислит, матем. и ма-тем. физ., 1962, т.2, Ш I, с.125-132

5. Пирумов У.Г., Рубцов В.А., Суворова В.Н. Расчет осесимметричных сопел с учетом равновесных физико-химических превращений. В кн.: Численные методы в газовой динамике. М., 1963, вып.2, с.48-60

6. Кацкова О.Н. Расчет равновесных течений газа в сверхзвуковых соплах. М.: Изд-во Вычисл. центра АН СССР, 1964. - 61 с.

7. Пирумов У.Г., Рубцов В.А. Расчет осесимметричных сверхзвуковых кольцевых сопел. Изв. АН СССР. Сер. механ. и маши-ностр., 1961, №6, с.15-25

8. Денисова Н.В. Численный расчет профилированных сверхзвуковых сопел с центральным телом. Тр. ЦАГИ, 1974, вып.1571.-53с.

9. Верховский В.П., Денисова Н.В., Межиров И.И. Расчет сверхзвуковых кольцевых профилированных сопел. Уч. зап. ЦАГИ, 1976, т.7, № 3, с.Ю8-ПЗ

10. Ю. Верховский В.П. Численное исследование течений в плоских соплах с угловой точкой при числах М = 3 5 и нерасчетных показателях адиабаты. Тр. ЦАГИ, 1976, вып.1802, с.3-32

11. Затолока В.В., Зудов В.Н., Шумский В.В. Расчетный анализ плоских несимметричных сопел при сверхзвуковой скорости на входе. Изв. СО АН СССР. Сер. техн. наук, 1972, вып.З, № 13, с.42-46,

12. Рылов А.И, Решение вариационной задачи о профилировании сопла, реализующего равномерный сверхзвуковой поток. Изв. АН СССР. Сер. механ. жидкости и газа, 1974, № 14, с.85-92

13. Рылов А.И. К построению сверхзвуковых несимметричных сопл. Изв. АН СССР. Сер. механ. жидкости и газа, 1977, № 3, с.87-94

14. Киреев В.И. О моделировании газовых струй, истекающих из реальных сопел Лаваля. Уч. зап. ЦАГИ, 1976, т.7, fe 2» с. 143-155

15. Иванов И.Я., Крайко А.Н. Профилирование двумерных и пространственных сопел и расчет течений в них. В кн.: Газодинамические лазеры и лазерная фотохимия. М.: Изд-во Моск. ун-та, 1978, с.52-61

16. Киреев В.И., Пирумов У.Г. О профилировании сопел модельных систем. Изв. АН СССР. Сер. механ. жидкости и газа, 1979,1. JI2 2, с. 187-189

17. Крайко А.Н., Шеломовский В.В. О профилировании плоских и осесимметричных сопел и каналов, реализующих заданный сверхзвуковой поток в сечении выхода. Изв. АН СССР. Сер. механ. жидкости и газа, 1981, № 2, с.94-Ю2

18. Войновский А.С., Киреев В.И. Численное решение задачи.Коши для внутренних течений. Гос. фонд алгоритмов и программ СССР, рев. № П004516, 1979

19. Войновский А.С., Киреев В.И. Профилирование сверхзвуковых сопел послойным методом характеристик. М., 1980. - 12 с. ВИНИТИ, деп. рук. 1980, Ш 2279

20. Войновский А.С., Киреев В.И. Обратная задача профилирования сверхзвуковых сопел и каналов численными методами. -Изв. АН СССР. Сер. механ. жидкости и газа, 1981, № 4, с.176 (семинары НИ ВЦ МГУ).

21. Пирумов У.Г. Расчет течения в сопле Лаваля. Изв. АН СССР. Сер. механ. жидкости и газа, 1967, N2 5, с.Ю-22

22. Пирумов У.Г., Росляков Г.С. Течения газа в соплах. -М.: Изд-во Моск. ун-та, 1978. 351 с.

23. Крайко А.Н., Шеломовский В.В. Сравнение двух методов профилирования контуров сверхзвуковых частей сопел, реализующих равномерный поток. Уч. зап. ЦАГИ, 1981, т.12, № 4, с.46-53

24. Подсыпанина Н.А., Шифрин Э.Г. Об одном методе профилирования коротких сопел. Изв. АН СССР. Сер. механ. жидкостии газа, 1975, № I, с.54-58

25. Подсыпанина Н.А. Использование плоскости годографа при профилировании численным методом осесимметричного сопла Лаваля.- Изв. АН СССР. Сер. механ. жидкости и газа, 1977, № I, с.164-168

26. Осипов И.Л. Численный метод построения двумерных сопл.- Изв. АН СССР. Сер. механ. жидкости и газа, 1979, № 2, с.179-183

27. Дзаннетти Л. Решение обратной задачи для внутренних течений методом установления.-Ракетная техника и космонавтика, 1980, т.18, № 7, с.167-172

28. Киреев В.И. Разностные методы решения обратной задачи профилирования сверхзвуковых сопел. М., 1980, - 21 с. ВИНИТИ, деп. рук. 1980, № 3302

29. Киреев В.И. Профилирование сверхзвуковых сопел на основе решения обратной задачи. М., 1980. - 15 с. ВИНИТИ, деп. рук. 1980, №2278

30. Войновский А.С., Киреев В.И. О смешанных краевых задачах профилирования сверхзвуковых сопел и каналов. Изв. АН СССР. Сер. механ. жидкости и газа, 1983, № 4, с.112-118

31. Войновский А.С., Киреев В.И. Численное профилирование плоских и осесимметричных сверхзвуковых каналов для неравновесных течений, реализующих неравномерные параметры на выходе. -Письма в ЖТФ, 1983, т.9, вып.16, с.Ю15-Ю18

32. Пчёлкина Л.В., Солодкин В.К. Корректировка влияния пограничного слоя на течение в соплах с изломом образующей. В кн.: Численные методы в газовой динамике. М., 1965, вып.4, с. ЮЗ-Ю7

33. Запрянов З.Д., Миносцев В.Б. Метод расчета пространственного обтекания тел сверхзвуковым потоком газа. Изв. АН СССР. Сер. механ. и машиностр., 1964, № 5, с.53-58

34. Чушкин П.И. Метод характеристик для пространственных сверхзвуковых течений. М.: Изд-во Вычисл. центра АН СССР, 1968. - 122 с.

35. Магомедов К.М., Холодов А.С. О построении разностных схем для уравнений гиперболического типа на основе характеристических соотношений. Ж. вычислит, матем. и матем. физ., 1969, т.9, № 2, с.373-386

36. Ашратов Э.А., Волконская Т.Г., Росляков Г.С., Усков В.И.

37. Исследование сверхзвуковых течений газа в струях. В кн.: Некоторое применение метода сеток в газовой динамике. М., 1974, с.241-407ричных сверхзвуковых течений при наличии необратимых процессов. М.: Изд-во вычисл. центра АН СССР, 1964. - 43 с.

38. J, Soc. Off*. is. £ лГ^ f. Я5-1М

39. Елисеев C.H. Модифицированный метод характеристик для расчета двумерных сверхзвуковых течений газа с выделением разрывов. Тр. ЦАГИ, 1983, вып.2199, с.З-Збнов В.В. Пространственное обтекание гладких тел идеальным газом.- М.: Наука, 1964. 505 с.

40. Овсянников A.M., Пирумов У.Г., Плетнева Е.М., Росляков Г.С. Атлас плоских сопел. М.: Изд-во Моск. ун-та, 1976.- 108 с.

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

42. Пирумов У.Г. Исследование двухслойных течений в сверхзвуковых осесимметричных соплах. Изв. АН СССР. Сер. механ. жидкости и газа, 1970, № 4, с.76-81

43. Абрамович Г.Н. Прикладная газовая динамика. М.: Наука, 1976. - 888 с.

44. Кочин Н.Е., Кибель И.Л., Розе Н.В. Теоретическая гидромеханика. М.: Физматгиз, 1963. - 728 с.

45. Кацкова О.Н., Крайко А.Н. Расчет плоских и осесиммет

46. Бабенко К.И., Воскресенский Г.П., Любимов А.Н., Руса<

47. Коган М.Н., Кучеров А.Н., Михайлов В.В., Фонарев А.С, Плоские течения газа при слабом подводе энергии. Изв. АН СССР Сер. механ. жидкости и газа, 1978, № 5, с.95-Ю2

48. Шахнов И.Ф., Фрост В.А., Обтекание плоской пластины сверхзвуковым потоком идеального газа при наличии объемного тепловыделения. Тр. МФТИ, 1961, № 7, с.ЦО-123

49. Кривцов В.М., Осипов И.Л. Численное решение задачи о профилировании дозвуковой части сопла Лаваля для течений излучающего газа. Ж. вычислит, матем. и матем. физ., 1981, т.21, № 5, с.1347-1352

50. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1973. - 848 с.

51. Годунов С.К. Уравнения математической физики. М.: Наука, 1979. - 392 с.

52. Колмогоров В.Ф. Численное решение обратной задачи теории сопла Лаваля применительно к двумерным неравновесным течениям совершенного газа. Изв. АН СССР. Сер. механ. жидкостии газа, 1974, № 2, с.136-142

53. Кацкова О.Н., Крайко А.И. Расчет плоских и осесиммет-ричных сверхзвуковых течений при наличии необратимых процессов. Прикл. механ. и техн. физика, 1963, № 4, С.П6-П8

54. Камзолов В.Н., Пирумов У.Г. Расчет неравновесных течений в соплах. Изв. АН СССР. Сер. механ. жидкости и газа, 1966,6, с.25-33

55. В^о^ъ jt-J'F. CLtL^a-ctz^tssbcs о^p. Z16-222

56. Кацкова O.H., Чушкин П.И. Пространственные сверхзвуковые течения газа с неравновесными процессами. Ж. вычислит.матем. и матем. физ., 1968, т.8, й 5, с.Ю49-Ю62

57. Киреев В.И. Пирумов У.Г. Расчет стационарных сверхзвуковых течений с неравновесными химическими реакциями. S. вычислит. матем. и матем. физ., 1980, т.20, № I, с.182-190

58. Термодинамические свойства индивидуальных веществ. Справочник в 2-х т. Под ред. акад. В.П.Глушко. T.I М.: Наука, 1962. - 875 с.

59. Кондратьев В.Н. Константы скорости газофазных реакций. М.: Наука, 1971. - 351 с.

60. Лунев В.В. Гиперзвуковая аэродинамика. М.: Машиностроение, 1975. - 327 с.

61. Кларк Д., Макчесни М. Динамика реальных газов. М.: Мир, 1967. - 566 с.

62. Елисеев С.Н. Таблицы контуров центральных тел осесим-метричных воздухозаборников с изэнтропическим торможением потока для чисел М 1,8 - 5. - Тр. ЦАГИ, 1972, вып.1409. - 46 с.

63. Ландау Л.Д., Лифшиц Е.М. Механика сплошных сред. М.: Гостехиздат, 1954. - 796 с.

64. Wezsc. А. ЖесъСе iZ&'cdcoAtiu^ty-ty^. -Л. oCe^SrAc /u/tfa&^tfa^scA. 194S; AO.

65. Курант P., Фридрихе К. Сверхзвуковое течение и ударные волны. М.; Изд-во Йн. лит., 1950. - 427 с.

66. Вуст В. К теории разветвленных скачков уплотнения. -В сб.: Газовая динамика, М., 1950, с.131-143

67. Шмыглевский Ю.Д. Некоторые вариационные задачи газовой динамики. М.: Изд-во Вычисл. центра АН СССР, 1963. - 141 с.

68. Росляков Г.С. Взаимодействие плоских скачков одного направления. В сб.: Численные методы в газовой динамике. М., Изд-во Моск. ун-та, 1965, с.28-51

69. Усков В.Н. Интерференция стационарных газодинамических разрывов. В сб.: Сверхзвуковые газовые струи. Новосибирск, 1983, с.22-45

70. Тугазаков Р.Я. Точное решение задачи взаимодействия движущегося со сверхзвуковой скоростью клина с границей раздела двух сред. Прикл. механ. и техн. физика, 1983, № 5, с.94-97

71. Землянский Б.А., Лесин А.Б., Лунев В.В., Шманенкова Г.А. Об интерференции косых скачков уплотнения одного семейства в гиперзвуковом потоке. Изв. АН СССР. Сер. механ. жидкости и газа, 1982, № 5, с.134-138

72. Крайко А.Н. Некоторые вопросы геометрической акустики одномерных нестационарных и двумерных стационарных течений. -Изв. АН СССР. Сер. механ. жидкости и газа, 1967, № 5, с.104-109

73. Черный Г.Г. Течения газа с большой сверхзвуковой скоростью. М.: Физматгиз, 1959. - 220 с.