Математическое моделирование гидрогазодинамических процессов при движении двухфазной среды в кольцевых соплах Текст научной статьи по специальности «Физика»

УДК 519.6

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

А.Л. Карташев, М.А. Карташева

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

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

В современных ракетных двигателях ттаптли широкое применение металлизированные ракетные топлива, одним різ компонентов которых являются металлы либо их соединения, при горении которых образуются окислы металлов, находящиеся в продуктах сгорания топлива., как правило, в виде жидких капель ил и твердых частиц. .\1ассо ви и доля таких конденсированных продуктов сгорания в общем расходе продуктов сгорания через сопло двигателя может достигать ~ 40%.

В ракетных топливах в качестве добавок к основному горючему используются следующие металлы: алюминий (А1), бериллий (Ве), бор (В), магний (М$), цирконий {Ъх), а также соединения металлов - гидрид алюминия (АІН3) и ряд других соединений, в частности гидриды упомянутых выше металлов. Применение металлизированных топлив в ракетном двигателе имеет ряд преимуществ по сравнению с безметалытыми топливными составами: во-первых, позволяет существенно повысить плотность используемого ракетного топлива и разместить в одном и том же корпусе ракетного двигателя большее количество топлива, во-вторых, величина удельного импульса тяги металлизированного топлива значительно превышает аналогичную величину для топлив без металлических добавок, что обусловливается более высокой энергетической эффективностью химических реакций, протекающих с участием металлизированных добавок и более высокой температурой продуктов сгорания в камере ракетного двигателя.

В продуктах сгорания могут одновременно присутствовать как окислы металла, находящиеся в различных агрегатных состояниях: твердом, жидком и газообразном, так и окислы различных металлов (если в состав металлизированного ракетного топлива входит несколько металлов или их соединений), которые также могут находиться в различных агрегатных состояниях. Таким образом, течение продуктов сгорания в сопле ракетного двигателя представляет собой течение многокомпонентной (многофазной) среды.

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

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

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

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

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

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

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

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

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

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

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

- присутствием в мттогокомпоттетттттом потоке фазовых переходов и химических реакций.

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

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

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

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

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

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

В данной статье рассматривается осесимметричное движение двухфазной среды при исходных допущениях, представленных в работах [4, о] (впервые предложены Клигелем в работе: Kliegel ,J.R. One-dimensional Flow of Gas-Particle System. - IAS Paper, 1960, № 5, 20 p.):

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

- совокупность частрщ пррпшмается непрерывно распределенной по всему объему с плотностью ps, равной отношению массы частиц в рассматриваемом объеме к величине этого объема;

- течетше стацргопарпое, массовый расход газа pi частрщ вдоль по течетшто постоянны;

- давлетше создается только газом, влрштшем частрщ можно пренебречь;

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

- сртстема теплортзолрфовапа, теплообмен между частацамрт pi газом осуществляется только путем копвекцрш;

- объемом, затшмаемым частрщамрт, а также граврттацргоппымрт pi электрртческршрт сртламрт можно пренебречь;

- вязкрге сртлы проявляются только nppi взапмодействрш частрщ с газом, ускоретше частрщ обусловлено действргем ЭТР1Х же СР1Л;

- вследстврге высокой теплопроводпострт матерртала частрщ pix температура постоянна по всему объему частрщ;

- газ является ртдеальпым, хршртческрт не реагртругопщм, а теплоемкость газа pi частрщ постоянны.

Поставлена задача математртческого моделртроватшя гртдрогазодршамртческртх процессов nppi двртжепрш двухфазных моподртсперспых потоков в кольцевых соплах внешнего распшре-Т1Р1Я (в том чртсле с укороченным центральным телом) nppi комплексном учете всех основных особенностей рассматрртваемого потока.

Двухфазное моттодисперсттое течение в кольцевом сопле описывается системой уравнений осесимметричного течения двухфазной смеси в интегральной форме [1, 4], позволяющая вести «сквозной> счет без предварительного выделения разрывов в расчетной области. Система записана для плоскости сечения Х,У для фиксированной пощади а, ограниченной контуром Г:

д JJ pydxdy + J py (udy — vdx) = 0; o г

dt JJ psydxdy + J psy (Usdy — Vsdx) = 0;

О

dt JJ Puydxdy + JJ psfxydxdy + J y ((p + pu2) dy — puvdx) = 0;

o o г

dt JJ Pvydxdy + JJ y (p. fy — dxdy + J y (puvdy — (p + pv2) dx) = 0;

o o Г

dt JJ p.u.ydxdy — JJ psfxydxdy + J yp.u. (u.dy — v.dx) = 0; (1)

o o Г

Jt JJ psvsydxdy — JJ psfxydxdy + J ypsvs (usdy — vsdx) = 0;

o o Г

py (Vf + e) dxdy + JJ psy V ■ f + ^ dxdy + J py (Vf + ^ (udy — vdx) = 0;

o г

JJ psesydxdy — JJ psqydxdy + J psesy (usdy — vsdx) = 0,

d_ dt

а а Г

где p и ps - плотности газа и <газа> частиц; p - давление газа; и и v - осевая и радиальная компоненты вектора скорости газа V; us и vs - осевая и радиальная компоненты вектора скорости частиц Vs; x и y - осевая и радиадьная координаты; e = e (p,T)- удельная внутренняя энергия газа; es = es (Ts) = свTs - удельная внутренняя энергия вещества частиц. св - удельная теплоемкость вещества частиц; p = pRT - уравнение состояния совершенного газа; t - время; f = Q = •

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

Расчет поля течетшя осуществляется методом устаповлепрш с ртспользоватшем вычртслрт-тельпого алгоррттма, построенного па основе разностной схемы С.К. Годунова - В.П. Колгапа [1, 6, 7]. Гратшчпые условрш сртстемы уравпепрш прртяты следугопщмрт:

- на жестких стенках принято условие непротекания vn = 0;

- на входной дозвуковой границе: энтропия S=const, полная энтальпия Ho = const; распределение угла наклона вектора скорости 0* = 0* (x,y).

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

Параметры частиц на входной границе us, vs, Ts совпадают с соответствующими пара-метрамрт газовой фазы.

Алгоритм расчета двухфазных моттодисперсттых потоков в кольцевых соплах реализован в модулях расчета параметров движения двухфазной моттодисперспой смеси в до-, трансзвуковой области течения и в сверхзвуковой области течения, который отличается особенностями вычислительного алгоритма, связанными с гиперболичностью системы уравнений во всей рассматриваемой области течения, и использованием, в связи с этим, вместо схемы С.К. Годунова - В.П. Колгатта модифицированной маршевой схемы М.Я. Иванова -

А.Н. Крайко - Н.В. Михайлова.

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

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

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

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

- ширина минимального сечения Н = 0, 2;

- угол наклона минимального сечения к оси сопла 0* = 40°:

длина сопла = 1, 65;

- радиус выходного сечения сопла Ка = 1, 57.

0 X

Рис. 1. Схема кольцевого сопла для параметрических расчетов

В качестве конденсированной фазы рассмотрены сферические частицы диаметром d = 1 — 10 мкм, массовая доля которых в двухфазной смеси г изменялась от 0,05 до 0,5.

Параметры газовой фазы приняты следующими:

- давление газа па входе в сопло р0 = 10 МП а (~ 100 атм):

- температура газа на входе в сопло Т0 = 3000 К:

- показатель изоэнтропы газа к = 1, 25;

- коэффициент динамической вязкости ц = 9, 5 • 10_5 (кг/м-с).

Число Прандтля Рг принято равным 0,75, характерный размер для расчета числа Рейнольдса Ее принят равным ширине щели кольцевого критического сечения и составляет ~ 0,05 м.

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

- зависимость Е и О от диаметра частиц ^ ^^^^^^^ганном значении г (рис. 2):

- зависимость Е и О от массовой доли г при фиксированном значении d (рис. 3):

- зависимость Е и О от диаметра частиц ^ ^^^^^^^ганном значении п (рис. 4).

Рис. 2. Зависимость Я и С от диаметра частиц Н при г = 0, 3

0.1 0.2 0.3 0.4 0.5 0.6 0.7 г 0.1 02 03 0.4 0.5 Об о.7 г

Рис. 3. Зависимость Я и С от массовой доли конденсированной фазы г при Н = 5 мкм

Анализ полученных результатов позволяет сделать следующие выводы. Установлена величина диаметра частиц конденсированной фазы, при котором значения тяги и расхода газа достигают своего максимального значения при фиксированной доле конденсированной фазы в двухфазном потоке, для данной конфигурации сопла и принятых параметрах рабочего тела эта величина составила ~ 5 мкм. Зависимости тяги и расхода от г при фиксированном диаметре частиц (в данном случае рассмотрены частицы с диаметром d = 5 мкм) и зависимости тяги и расхода от диаметра частиц d при постоянной концентрации частиц конденсированной фазы, в отличие от предыдущего случая, носят монотонный характер. При этом обнаружена существенная особенность этих зависимостей: слабое изменение тяги

в. кг/с

9.9

9.3

9.7

9.6

/

А мкм

Рис. 4. Зависимость Я и С от диаметра частиц Н при п = сопв1

и расхода при изменении диаметра частиц от й ~ 1 мкм до й ~ 3 мкм (при п = const) и существенный рост тяги и расхода при увеличении диаметра частиц свыпте 3 мкм. И аналогичный характер зависимости тяги и расхода от массовой доли конденсированной фазы г (при й = 5 мкм): слабый рост исследуемых параметров при изменении массовой доли от г ~ 0,1 до г ~ 0, 5, и существенный рост тяги и расхода при увеличении г свыше 0,5.

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

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

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

Необходимо отметить еще одну особенность численного моделирования течений в кольцевых соплах при наличии в рабочем теле конденсированной фазы, которая была обнаружена в ходе параметрических исследований. Для проведения численных исследований в настоящем исследовании была выбрана конфигурация полноразмерного кольцевого сопла сметанного расширения с внешней обечайкой и центральным телом. Рассмотренную конфигурацию можно назвать «гладкой:». Математическое моделирование течений в других, «негладких:», конфигурациях кольцевых сопел, показало наличие значительных вычислительных трудностей, возникающих при расчете характеристик таких конфигураций [1]. В этом случае при исследовании характеристик кольцевого сопла внешнего расширения с укороченным центральным телом за торцем центрального тела происходит пересечение траекторий частитт, конденсированной фазы. Для разрешения возникающей особенности можтто рассмотреть возникающую в потоке поверхность (в рассматриваемом осесимметричном случае она представляет собой «шнур»), которая имеет название «пелена» и модели расчета которой представлены в работе [2]. Решение, получаемое тта пелетте частитт,, будет влиять только тта параметры доттттой области, которые определяют величину тяги, создаваемой торцем укороченного центрального тела. Учитывая, что величина тяги торца невелика по

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

Литература

1. Карташев, А.Л. Математическое моделирование течений в кольцевых соплах: монография / А.Л. Карташев, М.А. Карташева. - Челябинск: Издательский центр ЮУрГУ, 2011. - 158 с.

2. Глазунов, A.A. Исследование неравновесных двухфазных течений в осесимметричных соплах Лаваля / A.A. Глазунов, А.Д. Гычков /'/ Изв. АН CCCF. МЖГ. - 1977. - № 6. -С. 86-91.

3. Васенин, И.М. Численное решение задачи о течении смеси газа и частиц в осесимметричном сопле Лаваля ./' И.М. Васенин, А.Д. Гычков /'/' Изв. АН СССГ. МЖГ. - 1973. -.№ 5. - С. 178-181.

4. Двухфазные мотто- и полидисперсттые течения газа с частицами / Л.Е. Стерттитт, Б.Н. Маслов, A.A. Шрайбер и др. - М.: Машиностроение, 1980. - 176 с.

5. Стерттитт, Л.Е. Основы газодинамики двухфазных течений в соплах / Л.Е. Стерттитт. -М.: Машиностроение, 1974. - 212 с.

6. Численное решение многомерных задач газовой динамики / С.К. Годунов, A.B. Забродит!, М.Я. Иванов и др. - М.: Наука, 1976. - 400 с.

7. Колгатт, В.П. Применение принципа минимальных значений производной к построению коттечтторазттостттых схем для расчета разрывных решений газовой динамики /

В.П. Колгатт // Уч. зап. ЦАГИ. - 1972. Т.З, № 6. - С. 68-77.

Александр Леонидович Карташев, доктор технических наук, профессор, кафедра «Летательные аппараты и автоматические установки:», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), al_kartashev@mail.ru.

Марина Анатольевна Карташева, кандидат технических наук, доцент, кафедра «Лета»

верситет (г. Челябинск, Российская Федерация), ma_kartasheva@mail.ru.

Bulletin of the South Ural State University.

Series «Mathematical Modelling, Programming & Computer Software»,

2013, vol. 6, no. 1, pp. 25-33.

MSC 35Q30, 76T10, 76T15

The Mathematical Modelling of Hydro-, Gas-Dynamic Processes at Flow of Two-Phase Medium in Annular Nozzles

A.L. Kartashev, South Ural State University, Chelyabinsk, Russian Federation, al_kartashev@mail.ru,

M.A. Kartasheva, South Ural State University, Chelyabinsk, Russian Federation, rna kartasheva@mail.ru

The flow of two-phase medium in annular nozzles is considered. The mathematical model for analysis of hydro-, gas-dynamic processes is introduced. The computational field of annular nozzles and system of equations for mathematical modelling of two-phase monodisperse flows in integral form are presented. The solving of system equations carry out by procedure of assignment with using scheme Godunov - Kolgan. The results of mathematical modelling of hydro-, gas-dynamic processes (dependences of propulsion of annular nozzle and flow rate of working medium from parameters of two-phase medium), and carries at flow of two-phase medium in annular nozzle with short-cut central body are presented.

Keywords: mathematical modelling, annular nozzle, hydro- gas-dynamic processes, system of equations, propulsion of annular nozzle, flow rate of working medium, parameters of two-phase medium.

References

1. Kartashev A.L., Kartasheva M.A. Maternaticheskoe rnodelirovanie techeniy v koltsevykh soplakh [Mathematical Simulation of Flows in Annular Nozzles]. Chelyabinsk, Yuzhno-Uralskiy gos. Univ., 2011. 158 p. (in Russian)

2. Glazunov A.A., Rychkov A.D. Research of Nonequilibrium Two-Phase Flows in Axisyinmetric Laval’s Nozlles [Issledovanie neravnovesnykh dvukhfaznykh techeniy v osesiinmetrichnykh soplakh LavalyaJ. Izvestiya AN SSSR. MZHG, 1977, no. 6, pp. 86-91.

3. Vasenin I.М., Rychkov A.D. Computational Solution of Problem of Flow Mixture of Gas and Particles in Axisyinmetric Laval’s Nozlle [Chislennoe reshenie zadachi о techenii srriesi gasa i chastits v osesiininetrichnom sople LavalyaJ. Izvestiya AN SSSB.. MZHG, 1977, no. 6, pp. 86-91.

4. Sternin L.E., Maslov B.N., Shrayber A.A., Podvysot.skiy A.M. Dvukhfasnye mono- i polidispersnye techeniy gasa s chastitsarni [Two-Phase Mono- and Polydisperse Flows of Gas with Particles]. Moscow, Machinostroenie, 1980. 176 pp. (in Russian)

5. Sternin L.E. Osnovy gasodinarniki dvukhfasnykh techeniy v soplakh [Fundamentals of Gas-Dynamic of Two-Phase Flows in NozzlesJ. Moscow, Machinostroenie, 1974. 212 p.

6. Godunov S.K., Zabrodin A.V., Ivanov M.Ya., Kraiko A.N., Prokopov G.P. Chislennoe reshenie rnnogornernykh zadach gasovoy dinarniki [Computational Solution of Multidimension Problem of GasdynamicJ. Moscow, Nauka, 1976. 400 p.

7. Kolgan V.P. Primenenie printsipa ininimalnykh znacheniy proizvodnoy к postroeniyu konechnoraznostnykh skhein dlya raschet.a razryvnykh resheniy gazovoy dinarniki [Application of Principle Minimal Values of Derivative for Construction Finite-Difference Schemes for Calculation of Discontinuous Solutions of Gas-DynamicJ. Uch. Zap. TSAGI, 1972, vol. 3, no. 6, pp. 68-77.

Поступила в редакцию 24 декабря 2012 г.