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

УДК 532.5, 533, 623.5

БОТ: 10.18698/2308-6033-2019-2-1852

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

© Н.В. Быков

МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

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

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

Введение. Разгон тел с помощью сжатого газа лежит в основе ряда технических устройств, к которым относятся лабораторные установки для аэробаллистических и ударно-прочностных исследований [1-4], спортивное, учебное-тренировочное и нелетальное оружие [5], пневматические линеметы [6], имитаторы запуска ракет с низкоскоростного носителя [7] и др. Указанный способ разгона отличают простота, безопасность, низкая стоимость и высокая воспроизводимость параметров [8], однако при этом в большинстве случаев скорости метаемых тел ограничены значениями 400...600 м/с [5-7]. В то же время имеются теоретические и экспериментальные исследования, показывающие возможность достижения скоростей более 1500 м/с для небольших масс (около 100 г) метаемых тел при использовании в качестве рабочего тела легких газов [3, 9, 10].

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

мулы, поскольку полное газодинамическое решение задачи даже в одномерной постановке довольно сложно в связи с наличием подвижной границы, которую представляет собой поршень [1, 2, 5, 6]. При этом следует различать две принципиально разные ситуации: для относительно тяжелых (по сравнению с массой сжатого газа) поршней проявление волновых эффектов в газе незначительно, поэтому ими можно пренебречь; в случае разгона относительно легких поршней наблюдается обратная картина — скорость движения тела сравнима со скоростью движения волн в газе, поэтому пренебрежение волновыми эффектами может существенно сказаться на результате расчета.

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

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

Постановка задачи и приближенные аналитические решения.

Схема разгона поршня с помощью сжатого газа показана на рис. 1. В начальный момент времени газ массой М занимает объем W0 (длину х0 в трубе диаметром й). Общая длина трубы вместе с газовой камерой — хт. Поршень массой т приводится в движение расширяющимся газом и разгоняется в трубе на длине хт - х0.

§ 1 1 Г

з_

-

-

Рис. 1. Схема разгона тела сжатым газом

Можно выполнить простую оценку достигаемой скорости, если записать условия баланса полной энергии системы «газ + поршень» в приближении однородной деформации (скорость газа вдоль оси канала трубы распределена линейно) при равномерном распределении массы газа по длине камеры:

1 + ^

тиР = МО _ Рту 2 к _ 1 к _ 1

где ц = M/m — отношение массы газа к массе поршня; up - скорость метаемого тела; p0 — давление газа в начальный момент времени;

1 xp

Pm (t) = — | p (x, t) dx — среднее давление в запоршневом простран-

xp о

стве в момент времени t; p(x, t) — истинное давление в сечении с координатой x в момент времени t; W = xpS, W0 = x0S — текущий и начальный объемы запоршневого пространства соответственно; S — площадь поперечного сечения трубы; xp — текущая координата поршня; k — показатель адиабаты газа.

Предполагая, что газ расширяется адиабатически, т. е. что

pmWk = p0Wok, можно получить

up =

dÇ p

dx \

2ц (l - Ç P9 )

(1)

0(9 +1)(1 + ц/3) '

где ир = ир / с0 — безразмерная скорость поршня; с0 = ^кЯ^ГО —

скорость звука в газе в начальный момент времени; — газовая постоянная; 20 — начальная температура газа; р = хр / х0 — безразмерная координата; т = С / х0 — безразмерное время; 0 = к - 1.

Приведенное уравнение позволяет сделать некоторые предельные оценки. Так, для сильного расширения газа р >> 1 и можно получить

и,

maxi

(i)= _2t

p \е (е+i)(i+ц/з)'

а в случае, если дополнительно ц >> 1, оценка зависит только от е:

max(2)

Up 1 ' =

9 (9 +1) '

Для холодных газов эта скорость составляет: для углекислого газа — 1064 м/с; для воздуха — 1136 м/с; для гелия — 2363 м/с; для водорода — 4308 м/с. Хотя эти значения скоростей являются теоретическим пределом, они позволяют быстро оценить возможность достижения требуемых в той или иной практической задаче скоростей на выбранном газе.

Уравнение (1) может быть проинтегрировано, в результате чего получится уравнение траектории поршня:

Çp F„ f I, ß, ,+ß, ÇP ß VÄßU

Г(1/2 + ß) V9 (9 + l)(l + ц/3)

6

гда ^ у) —

ж 1

Fh (аДу, у ) = 1+ £ -

п=1

п!

гипергеометрическая П (т + а)(т + в)

т=0

(т + у)

функция; уп; в = -1/0; Г(у) — гам-

ма-функция Эйлера.

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

Среди методов нахождения приближенных решений баллистической задачи Лагранжа наибольшее распространение получило точное решение газодинамической задачи в области простой волны [1, 2]:

ир =-

р к -1

к-1\

1 + -

к +1 2к

-цт

к+1

(2)

§ р = 1 +

к -1

к

т + — Ц

1 + -

к +1 2к

цт

к+1

(3)

Решение справедливо до момента ?а (или тА в безразмерном варианте), соответствующего возвращению отраженной от дна канала волны на поршень (рис. 2) [2]:

к +1

т а = 2 +--ц .

А 2к

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

С практической точки зрения гораздо важнее времени та расстояние хА - х0, которое успевает пройти поршень за это время (в литературе не приводится). Подставляя выражение для та в формулу (3), получаем

£рА = 1+

к -1

„ к +1 ц к

2 +-- + -

2 к

ц

\\

1 + -

к+1ц| к+1ц

2 к

2 + V 2 к

\\

к+1

/У

Для различных значений показателя адиабаты это функция одной переменной — ц. Графики этой функции при трех значениях к показаны на рис. 3.

Рис. 3. Зависимость положения поршня от значения ц в момент времени та:

1 — углекислый газ (к = 1,3); 2 — водород, азот, воздух (к = 1,4); 3 — гелий (к = 1,67)

Отсюда следует важный вывод: при ц >> 1 решение в области простой волны в большинстве практических случаев справедливо до момента вылета поршня из канала ствола. Основное неудобство решения в форме (2), (3) состоит в том, что оно не показывает связи между скоростью и положением поршня, поэтому момент времени, соответствующий определенной координате, можно найти, только решая уравнение (3) численно.

Формула (2) позволяет получить другую предельную оценку скорости поршня при ц>>1:

и

тах(3) _

2

к -1

Это известная в газовой динамике формула предельной скорости нестационарного истечения газа в вакуум, в соответствии с которой предел ц>>1 в сущности можно воспринимать как отсутствие порш-

ня, поэтому и

тах(1) < итах(2) < ^;шах(3)

< и.

Один из наиболее популярных в зарубежных исследованиях метод решения основан на приближении центральной волны разрежения, создаваемой движущимся поршнем [11, 12]. В этом случае уравнение движения поршня можно записать в виде

йир йт

= и,

йир

к -1

2к к-1

и,

Это уравнение может быть проинтегрировано относительно и, тогда

ь р=1+

2к / ц

к +1

к -1 к-1

к +1Л к -1

1--ир

1

к -1

к+1 к-1

+1

ир

(4)

Неудобство этой формулы заключается в ее неявности. Легко показать, что формулы (2) и (3) представляют собой параметрическую форму зависимости (4), поэтому расчеты по этим формулам полностью совпадают, однако в обоих случаях необходимо решать численно нелинейные уравнения.

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

Математическая модель и численный метод решения. Учесть волновые явления в газе позволяет одномерная газодинамическая модель, которая в векторной безразмерной форме может быть записана как

дц дГ п -± + — = 0. дт д§

(5)

где

г

ги

V

г(е + и2 /2)

ги

о + ги

2

V

ги (е + и2 /2)

+ ои

/

Здесь введены следующие безразмерные переменные:

х и Жо

§ = —, и = —, о = р—2 хо Со тсо

г = р-

Жо

е =

т

где 8 — внутренняя энергия единицы массы газа.

л

8

С

Система (5) дополняется уравнением состояния идеального газа в безразмерной форме а = 9ге и уравнением движения поршня

йир

йт

а,

(6)

где нижний индекс р характеризует значение соответствующих функций при х = хр.

Начальные условия (при т = 0) имеют следующий вид: и = 0, г = ц, а = ц / к. На левой границе (дно трубы, В = 0) ставится граничное условие непротекания, на правой границе (подвижный поршень, В = ВР) — равенства скоростей газа и поршня, при этом скорость поршня определяется из (6).

Как видно, отношение масс газа и поршня (ц) входит в начальные условия. Таким образом, решение задачи является двухпараметриче-ским и зависит от параметров к и ц. При условии задания конкретного газа остается один параметр — ц.

Численное решение системы (5) осуществляется с помощью двухшаговой схемы типа «предиктор-корректор» на подвижной сетке, отслеживающей положение поршня:

г И \

и+1/2

Чг+1/2

лви

и+1

Чг+1/2

ЛВ

лв

и+1/2

чИ+1/2 -

2ЛВ'

{Гш - $)

ЛВ

и+1

Чг+1/2 -

ЛВ'

(^ и+1/2 — ^ и+1/2 )

Здесь индексы и и г отвечают временным и координатным узлам

И А V И

разностной сетки соответственно; т — шаг по времени; Лс, — шаг по координате на и-м временном слое. Сетка растягивается равномерно, что приводит к появлению множителя перед скобкой, который отражает деформацию ячеек сетки за один шаг по времени.

Потоки на границах ячеек (интерфейсах) определяют решением задачи о распаде разрыва методом ЛИ8М+, модифицированным для подвижной сетки:

Г 0 л

с

Гг = 2 (МГ (Ф г+1/2 + Фг-1/2)

)- Мг (Фг+1/2 -Фг_1/2)) +

аг

аи

где сг — скорость звука на г-м интерфейсе; Мгг — число Маха относительно г-го интерфейса; Ф = (г, ги, гН )т; Н — полная энтальпия

единицы массы газа; иг — скорость перемещения г-го интерфейса. Число Маха на интерфейсе находят по формуле

мг = /р+ (мг_1/2)+/р- (мг+1/2),

где

Мг

_ и/±1/2 - иI . „ _ с/+1/2 + с1-1/2 .

/±1/2 - ; с/ - I ;

с 2

' т (м ± 1М1), М > 1;

±}(М± 1)2 (1 + 4в(М +1)2), |М| < 1. Давление на интерфейсе определяют с помощью выражения

О/ - Ш +(МГ-1/2 ) о/-1/2 + Ш-(М/*+1/2 ) о/+1/2 ,

где

(M ± |м|) /2M, |М| >1;

(M± 1)2 (2 + M) / 4± аМ (М +1)2 ]), \и\ < 1.

Ша± (М) -

Значения констант: а = 3/16, р = 1/8 [14].

Адекватность одномерного газодинамического описания показана в работах [7, 9, 10], где проведено сравнение результатов расчета вылета поршней разной массы из баллистических установок на сжатом газе с экспериментальными данными и расчетами по двумерной теории. Расчеты по приведенной модели с использованием предложенного численного метода также дают хорошее согласие с результатами указанных работ.

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

Рис. 4. Зависимость ир от £ и ц при к = 1,4

Видно, что возрастание скорости при ц > 5 быстро замедляется. Кроме того, при числе объемов расширения газа более = 2-3 также не наблюдается значительного возрастания скорости — очевидно, что эффективность разгона по мере расширения быстро падает.

Для оценки практических областей использования полученной зависимости проведен расчет достигаемых скоростей и необходимых объемов газового резервуара для разгона тела массой 1 г в трубе диаметром 5 мм при давлении в резервуаре р0 = 50 МПа для различных холодных газов (Г0 = 300 К) при = 4 и ц = 5.

В таблице указаны параметры газа, необходимые для решения прямой задачи, а также соответствующая приведенная длина газового резервуара х0, общая длина трубы хт, соответствующая значению = 4, и достигаемая дульная скорость метаемого тела, рассчитанная различными методами: ир1 — путем численного решения системы уравнений (5), (6); ир2 — путем численного решения системы уравнений (2), (3), описывающих решение в области простой волны, и ир3 — по формуле (1).

Таблица

Оценка геометрических размеров трубы и скоростей метаемых элементов для различных газов

Газ к ^, Дж/(кг • К) С0 х0 х т иР2

м/с м м/с

Углекислый газ 1,30 189 271,5 0,30 1,20 488,7 491,4 491,1

Воздух 1,40 287 347,2 0,46 1,84 590,2 583,3 586,2

Гелий 1,67 2078 1020,3 3,36 13,44 1428,4 1435,1 1455,7

Водород 1,40 4125 1316,2 6,68 26,72 2237,5 2211,3 2222,1

Данные таблицы демонстрируют резкое различие между обычными газами, такими как воздух и углекислый газ, и легкими — водородом и гелием, что особенно заметно проявляется в необходимых объемах резервуара и, как следствие, длинах трубы. Отсюда можно заключить, что на практике для легких газов в одноступенчатых баллистических установках (без предварительного сжатия газа поршнем, как это делается в двухступенчатых установках) режимы работы при значениях ц >>1 слабо достижимы.

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

решением (5), которое обозначим ипр, и аналитическим решением (1), которое обозначим иар:

д _ п а

1\и р — и р и р .

Зависимость Аир от и ц при к = 1,4 показана на рис. 5.

1 О

Рис. 5. Зависимость Аир от В и ц при к = 1,4

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

Рис. 6. Зависимость Аир от ц при В = 4 и к = 1,4

Из приведенной зависимости следует, что выходная скорость, рассчитанная по формуле (1), может довольно существенно отличаться от полученной численно при ц > 10 в сторону занижения. В то же время при 0 < ц < 5 различие в результатах расчета по формуле (1) и газодинамического расчета не превышает ±0,03с0 (см. таблицу).

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

ческую энергию метаемого тела. Энергоэффективность любой системы разгона тел может быть оценена с помощью КПД, равного в данном случае [1]:

_ е тир _ е тир _ е (е+1) ир _ е (е+1) 2 П _ 2 ~ 2 МЁГ _ ц с2 _ ц ир.

Зависимость КПД от В и ц показана на рис. 7.

Из рисунка ясно, что КПД может быть достаточно большим при малых значениях ц (п - 1 и менее), однако с увеличением ц быстро снижается. Это объясняется, в частности, тем, что при расширении газа энергия тратится не только на разгон метаемого элемента, но и на ускорение самого газа. Кроме того, при больших значениях ц большая часть газа вообще не успевает поучаствовать в разгоне.

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

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

Л

10 1

Рис. 7. Зависимость КПД от В и ц при к = 1,4

женная от дна канала трубы волна давления не успевает догнать поршень.

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

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

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

Работа выполнена при частичной финансовой поддержке РФФИ (проекты № 16-38-00948 мол-а и 16-29-09596 офи-м).

ЛИТЕРАТУРА

[1] Seiler F., Igra O., ed. Hypervelocity Launchers. Springer, 2016, 300 p.

[2] Златин Н.А., Красильщиков А.П., Мишин Г.И., Попов Н.Н. Баллистические установки и их применение в экспериментальных исследованиях. Москва, Наука, 1974, 344 с.

[3] Павленко А.В., Балабин С.И., Козелков О.Е., Казаков Д.Н. Одноступенчатая газовая пушка для изучения динамических свойств конструкционных материалов в диапазоне до 40 ГПа. Приборы и техника эксперимента, 2013, № 4, с. 122-124.

[4] Jones D.R., Chapman D.J., Eakins D.E. A gas gun based technique for studying the role of temperature in dynamic fracture and fragmentation. J. Appl. Phys., 2013, vol. 114, p. 173508.

[5] Horak V., Do Duc L., Vitek R., Beer S., Mai Q.H. Prediction of the Air Gun Performance. Advances in Military Technology, 2014, vol. 9 (1), pp. 31-44.

[6] Григорьев В.В., Исаков С.Н., Петров Р.Л., Юркин С.В. Газодинамическое исследование пневматического линемета. ЖТФ, 2006, т. 76, вып. 3, c. 75-80.

[7] Johnston A., Krishnamoorthy L.V. A Numerical Simulation of Gas Gun Performance. Defence Science and Technology Organisation, 2008, DST0-TN-080.

[8] Plassard F., Mespoulet J., Hereil P. Analysis of a single stage compressed gas launcher behaviour: from breech opening to sabot separation. 8th European LS-DYNA Users Conference, 2011, vol. 8, pp. 1-11.

[9] Moradi A., Ahmadikia H. One-Dimensional and Axisymmetric Numerical Simulation of a Single-Stage Gas Gun. Adv. Theor. Appl. Mech., 2011, vol. 4 (3), pp. 101-111.

[10] Moradi A., Khodadadiyan S. Study of Real Gas Behaviour in a Single-Stage Gas Gun. Int. J. of Mechanical andMechatronicsEng., 2011, vol. 5 (6), pp. 948-952.

[11] Seigel A.E. The theory of high speed guns. NATO AGARD, 314 p.

[12] Gardiner P.A., Egawa Y., Watanabe K. Performance evaluation of single stage diaphragmless vertical gas gun for nitrogen and helium gas propellants. Mechanical Engineering Journal, 2016, vol. 3 (6), pp. 16-00273.

[13] Denny M. Gas gun dynamics. Eur. J. Phys., 2013, vol. 34, pp. 1327-1336.

[14] Luo H., Baum J.D., Lohner R. On the computation of multi-material flows using ALE formulation. J. of Comp. Phys., 2004, vol. 194, pp. 304-328.

Статья поступила в редакцию 08.02.2019

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

http://dx.doi.org/10.18698/2308-6033-2019-2-1852

Быков Никита Валерьевич — канд. техн. наук, доцент кафедры «Ракетные и импульсные системы», заместитель заведующего кафедрой по научной работе, заведующий отделом «Ракетные и ствольные системы» Научно-исследовательского института специального машиностроения МГТУ им. Н.Э. Баумана. Автор более 60 научных работ в области исследования внутрикамерных процессов, газовой динамики, внутренней баллистики и методов синтеза технических систем. e-mail: bykov@bmstu.ru

Comparative analysis of numerical and approximate analytical solutions to the problem of piston acceleration in a pipe by compressed gas

© N.V. Bykov

Bauman Moscow State Technical University, Moscow, 105005, Russia

The paper focuses on mathematical and numerical simulation of piston acceleration by compressed gas in a guide pipe, also known as Lagrange ballistic problem. We consider two approaches to an approximate analytical solution of the problem: the uniform deformation approximation, within which a new exact dimensionless solution for the piston trajectory is obtained, and a solution in the simple wave region, for which new estimates of applicability limits are obtained. The wave motion of a gas is described with the numerical solution of gas dynamics dimensionless equations on a moving mesh. We analyze the working areas of approximate analytical solutions comparing the latter with the numerical solutions. Furthermore, we demonstrate the calculation of a model example, which compares the results of calculations for various models and estimates the variation range in dimensionless parameters from the point of view of practical application. The dimensionless representation of the solutions allows us to summarize the results obtained. Based on the efficiency analysis, we deal with the issue of energy performance of piston acceleration by compressed gas. Findings of research show that with an increase in the length of the pipe and the ratio between the mass of compressed gas and the mass of the piston, the efficiency decreases.

Keywords: gas-dynamic acceleration, Lagrange ballistic problem, acceleration methods, uniform deformation, hypervelocity launcher, interior ballistics, ballistic installations

REFERENCES

[1] Seiler F., Igra O., ed. Hypervelocity Launchers. Springer, 2016, 300 p.

[2] Zlatin N.A., Krasilschikov A.P., Mishin G.I., Popov N.N. Ballisticheskie ustanovki i ikh primenenie v eksperimentalnykh issledovaniyakh [Ballistic ranges and their application in experimental research]. Moscow, Nauka Publ., 1974, 344 p.

[3] Pavlenko A.V., Balabin S.I., Kozelkov O.E., Kazakov D.N. Pribory i tekhnika eksperimenta — Instruments and Experimental Techniques, 2013, no. 56 (4), pp. 482-484.

[4] Jones D.R., Chapman D.J., Eakins D.E. A gas gun based technique for studying the role of temperature in dynamic fracture and fragmentation. J. Appl. Phys., 2013, vol. 114, p. 173508.

[5] Horak V., Do Duc L., Vitek R., Beer S., Mai Q.H. Prediction of the Air Gun Performance. Advances in Military Technology, 2014, vol. 9 (1), p. 31-44.

[6] Grigor'ev V.V., Isakov S.N., Petrov R.L., Yurkin S.V. Zhurnal tekhnicheskoy fiziki — Technical Physics. The Russian Journal of Applied Physics, 2006, vol. 51, no. 3, pp. 367 372.

[7] Johnston A., Krishnamoorthy L.V. A Numerical Simulation of Gas Gun Performance. Defence Science and Technology Organisation, 2008, DST0-TN-080.

[8] Plassard F., Mespoulet J., Hereil P. Analysis of a single stage compressed gas launcher behaviour: from breech opening to sabot separation. 8th European LS-DYNA Users Conference, 2011, vol. 8, pp. 1-11.

[9] Moradi A., Ahmadikia H. One-Dimensional and Axisymmetric Numerical Simulation of a Single-Stage Gas Gun. Adv. Theor. Appl. Mech., 2011, vol. 4 (3), pp. 101-111.

[10] Moradi A., Khodadadiyan S. Study of Real Gas Behavior in a Single-Stage Gas Gun. Int. J. of Mechanical and Mechatronics Eng., 2011, vol. 5 (6), pp. 948-952.

[11] Seigel A.E. The theory of high speed guns. NATO AGARD, 314 p.

[12] Gardiner P.A., Egawa Y., Watanabe K. Performance evaluation of single stage diaphragmless vertical gas gun for nitrogen and helium gas propellants. Mechanical Engineering Journal, 2016, vol. 3 (6), pp. 16-00273.

[13] Denny M. Gas gun dynamics. Eur. J. Phys., 2013, vol. 34, pp. 1327-1336.

[14] Luo H., Baum J.D., Lohner R. On the computation of multi-material flows using ALE formulation. J. of Comp. Phys., 2004, vol. 194, pp. 304-328.

Bykov N.V., Cand. Sc. (Eng.), Assoc. Professor, Department of Missile and Kinetic Warfare Systems, Deputy Head of Department for Scientific Work, Head of Department of Missile and Barrel Systems, Institute of Special Engineering, Bauman Moscow State Technical University. Author of over 60 scientific papers in the fields of gas dynamics, interior ballistics and optimization techniques in engineering systems. e-mail: bykov@bmstu.ru