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

УДК 629.762

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

© А.В. Плюснин МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

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

Ключевые слова: летательный аппарат, подводный выброс, кольцевое пространство, втекание воды

Введение. Задачи исследования и математического моделирования нестационарных взаимодействий газовой и жидкой сред представляют интерес для разработки энергетических установок, технологических процессов, создания, эксплуатации и применения различных технических систем. В этих приложениях наиболее широко распространены явления, связанные с динамикой парогазовых включений в жидкости [1-6], кавитацией [7, 8], струйными взаимодействиями [9, 10], испарением и конденсацией [11-13]. Оригинальным решением при экспериментальном определении параметров ударного проникания тел в воду является моделирование сжимаемости воды с помощью специально приготовленной смеси несущей жидкости с пузырьками газа [14]. Стоит отметить, что необычные явления, похожие на перечисленные, встречаются и в живой природе [15, 16]. Дополнительным стимулом к физическому исследованию этих процессов и развитию средств их математического и компьютерного моделирования [17, 18] служит присущее этим явлениям фрактальное изящество [19, 20], обусловленное динамикой свободных границ жидкости.

Тематика подводного старта летательных аппаратов (ЛА) содержит множество задач экспериментального определения и расчета параметров нестационарных газожидкостных течений и их воздействий на деформируемые конструкции ЛА и подводного носителя [21-23]. «Минометная» схема [24] подводного газодинамического выброса (рис. 1) имеет свои особенности [25]. Рассмотрим процесс, который может реализоваться при разгерметизации кольцевого пространства (КП)

между корпусом ЛА и внутренней стенкой пускового контейнера (ПК). Если не предпринимать специальных мер, таких, например, как предварительный наддув КП до давления, примерно равного наружному гидростатическому давлению рн = ратм + рв gH [26, 27] (где

ратм — атмосферное давление; рв — плотность воды; g — ускорение силы тяжести; Н — глубина погружения по верхнему срезу ПК), то под действием перепада давлений рн/ратм вода начнет быстро заполнять КП, сжимая находящийся там газ. Повышение давления в КП, которое при этом может реализоваться, способно привести к разрушению корпуса ЛА. Иногда это явление называют «гидроударом», хотя оно и не вполне соответствует процессу, наблюдаемому при гидравлическом ударе в трубах [28].

1

■Ратм ,

Рш®

Рис. 1. Схема подводного газодинамического выброса летательного аппарата (ЛА):

1 — поверхность воды; 2 — верхний срез контейнера; 3 — ЛА; 4 — начальное положение ЛА; 5 — задонное пространство (ЗП) контейнера; 6 — энергоустройство, создающее избыточное давление

рЗП (/); 7 — кормовое уплотнение ЛА; 8 — кольцевое пространство

контейнера; 9 — начальное давление р0 до разгерметизации;

10 — переднее уплотнение контейнера после разрушения; 9 — угол наклона контейнера к горизонту

Для проектных целей представляет интерес определение максимального давления р^, реализующегося в КП при втекании воды

под действием перепада давлений рн /р0, где р0 — начальное давление газа в КП. Случай ратм < р0 < рн соответствует частичному

предварительному наддуву КП. Для того чтобы гарантировать необходимый уровень давления, наддув осуществляют, подавая газ непосредственно в КП [26]. При этом учитывается также начальный свободный объем ЗП, куда газ из КП способен вполне свободно перетекать под кормовым уплотнением ЛА (см. рис. 1). Состав среды для задачи, рассматриваемой в данной статье, не имеет значения, хотя на практике во избежание вторичного догорания в ЗП продуктов сгорания энергоустройства газодинамического выброса [24, 25, 29] предпочтительнее производить наддув газом, не содержащим кислород (азотом) [26].

В данной работе процесс втекания воды в КП исследуется для наихудшего случая р0 = ратм в рамках одномерной расчетной схемы [25, 30], построенной на балансовых соотношениях массы и импульса для столба воды в КП. Основное внимание уделено рассмотрению внешней гидродинамической задаче и оценке влияния на величину Ртах упругости стенок ПК и корпуса ЛА. При этом движение ЛА в ПК считается не зависящим от дополнительного сопротивления, возникающего вследствие втекания воды в КП. На рис. 2 приведены графики пути Ь (1), скорости V (2), ускорения А (4) ЛА и давления

в ЗП Рзп в функции времени X, использованные в данной работе

при проведении расчетов. Учет связанности задачи не составляет труда. В практических приложениях к расчетной схеме стоит лишь добавить уравнения, описывающие движение ЛА и процессы в ЗП [23, 29, 31]. Угол $ наклона ПК к горизонту (см. рис. 1) принят вертикальным.

0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 Время г

Рис. 2. Заданные параметры движения летательного аппарата в пусковом контейнере: 1 — путь Ь; 2 — скорость V; 3 — давление в задонном пространстве Рзп

4 — ускорение А

Предполагается, что в ненагруженном состоянии боковые поверхности ПК и ЛА представляют собой цилиндры толщиной

сПК Г)ПК г>ПК сЛА Г)ЛА г>ЛА ^^^ л л „ О пПК г>ЛА

о = Квнеш - К и о = К - кввут соответственно. Здесь к , Лвнут —

знут

внутренние, а кв^КШ, КЛА — внешние радиусы ПК и ЛА. Учет прогибов

стенок ПК и ЛА под действием избыточного давления среды учитывается в упрощенной постановке, мотивированной работой [28]. Однако в отличие от работы [28] вода как снаружи ПК, так и в КП предполагается несжимаемой. Такой подход позволяет избежать решения задачи в частных производных, и его можно расценивать как один из начальных этапов исследования и математического моделирования в рамках методологической стратегии, изложенной в работе [32]. Развитием этого подхода является описание деформаций стенок ПК и ЛА на основе динамической теории оболочек [33, 34]. Если учесть сжимаемость воды [28, 35], получим постановку задачи, естественную с точки зрения общей теории механики сплошных сред [36].

Кормовое уплотнение ЛА (см. рис. 1) считаем непроницаемым для среды в обе стороны. Это означает, что горячие продукты сгорания из ЗП не могут повлиять на термодинамическое состояние газового объема в КП, а газ из КП не может перетекать в ЗП под действием положительного перепада давлений, что вполне согласуется с начальной стадией анализа физических процессов в рамках стратегии [32].

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

Внешняя гидродинамическая задача. Для составления математического описания процесса втекания воды в КП воспользуемся идеализированной схемой, приведенной на рис. 3. Считаем, что до начала этого процесса жидкость занимает верхнее полупространство над

плоскостью х = Ь

ПК

отождествляемой

с наружной поверхностью корпуса подводного носителя ЛА. На бесконечности жидкость покоится и имеет давление,

Рис. 3. Идеализация втекания равное давлению рн на уровне верхнего

воды в к°лъцев°м пр°странстве среза пускового контейнера (ВСПК) в

постоянного поперечного ~ л

исходной физической постановке задачи

сечения

(см. рис. 1). Весомость жидкости учитывается только для столба воды, проникшего в КП.

Под действием начального перепада давлений рн /р0 формируется течение воды, направленное внутрь ПК. Скорость втекания воды в КП Увтек и давление рвтек, которое она приобретает на уровне ВСПК, требуется определить в ходе решения задачи. Наряду с величиной Увтек будем использовать объемный расход воды внутрь КП

Q(t) = увтек2™к, где Е^ — площадь сечения КП на уровне ВСПК.

Заметим, что эти величины отрицательны при фактическом втекании воды внутрь КП.

Данную задачу и родственные ей [25, 37] обычно рассматривают с использованием двух модельных гидродинамических течений. Одно из них — течение от точечного источника/стока [38], непосредственно определяемое фундаментальным решением уравнения Лапласа в №3 [39]. В видоизмененной форме оно описывается уравнением Релея и характеризует динамические колебания сферического газового пузырька, находящегося в жидкости [2, 40]. Другое течение определяется решением Ламба задачи о вытекании жидкости из полупространства через круговое отверстие [41].

Решение Ламба в большей степени соответствует рассматриваемой задаче. Имея в виду определенную работу с ним, дадим его вывод, следуя решению задачи о потенциале электростатического поля, созданного диском [42], имеющей тождественную математическую формулировку. Двойственной по отношению к рассматриваемой задаче является задача об осесимметричном безотрывном обтекании диска [41] или, что эквивалентно, линеаризованная задача о плоском ударе диска по поверхности идеальной несжимаемой жидкости [35]. Замкнутую форму решения этих задач обычно приводят только для границы области определения, ограничиваясь для внутренних точек интегральным представлением, которое неудобно для вычислений [41, 42]. Однако благодаря формулам [43] все решение можно представить в простой замкнутой форме [30, 44].

В постановке, близкой к рис. 3, пока игнорируя присутствие ЛА,

_ д2ф 1 дф д 2ф п будем искать решение уравнения Лапласа —- н----1--- = 0 в по-

дг г дг дг

лупространстве г = х - ЬПК > 0, подчиняющееся условиям ф| дПК = С(X),

-л

— = 0 при г = 0 и ф| г ^ = 0 на бесконечности (где г — рассто-

дх г >ЯПК " ^

яние от оси ПК). Применяя преобразование Ф (р, г) =

= |ф(r, z) J0(rp)rdr Ханкеля нулевого порядка [42], для изображе-

0

ния Ф (р, z) получаем уравнение d2Ф/dz2 -р2Ф = 0, решением которого с заданным поведением на бесконечности является функция Ф (р, z) = A (p)e-pz. Применяя к ней обратное преобразование, имеем:

ад

ф (r, z) = JA (p)e-pzJo(rp) p dp. (1)

0

Подстановка выражения (1) в граничные условия приводит к системе дуальных интегральных уравнений

ад

| A(p) Jo(rp) p dp = C(t), r < RП

0 ад

J A(p) Jo( rp) p2dp = 0, r > Rn

0

2C(t) sin (R™p)

откуда A(p) = ——-^-2-- [45]. Ее подстановка в выражение (1)

к p

и формулы [43] дают:

Ф(г, z) = 2C(t) Jp-1e-pz sin (RnKp) J0(rp) dp = ^^arcsin —; (2) ni v ' к A

£(r,z)=-Je-pzSin(RnKp)J0(rp)dp=-^Ш^-;,-' , (3)

dz к 0 к Ay 'A[+>

где A = 0,5(A(-)+ A(+)), A(±)^z2 +(r±RПК)2.

Для получения предельных значений выражений (2), (3) при z = + 0 предварительно найдем предельные значения A( ) = r - R

A(+) = r+RПК. Теперь, если r < RПК, получаем:

A(-)= RПК - r, A=RПК, ф=arcsin 1 = C (t), dф=0,

к dr

2 2 2 2 - r2 - z2

ПК

с"ф_ 2C (t) ^¡(R ™)2 -r' = 2C (t)

ft к (R пк - r )(RПК + r) xJR^)- r 2

2

- r"

Если r >R , имеем:

Л-) г>пк а 2C(t) . Rn

A ' = r-R , A=r, ф=—— arcsin —

к r

5ф=- 2 R ПКС (t) 5ф=0 ör к r^r2 -(RnK )2 ' & '

Интегрирование по площади входного сечения дает объем-

dz

ный расход

RnK г 2 1-/2 Q(t) = -4С(t) j r (Rпк ) -r2 dr = -4C(t)R

ПК

Q(t) 4R

входного сечения ПК выражаются формулами

т. е. С(t) = - , „ПК . В итоге параметры внешнего течения в точках

„ext ^W ext р. ext ÜW i л\

ф =-4RK, vr =0, vz =-/ . 2 . (4)

4R 2KRn^(RnK )2 - r2

Полученное математическое решение задачи Ламба справедливо при любом знаке Q(t), но физический смысл оно имеет лишь в случае истечения жидкости из полупространства. При втекании жидкости в полупространство реальная картина течения становится совершенно иной [46]. Вместо расталкивания жидкости во все стороны, что имеет место при течении от точечного источника, втекающая в полупространство струя в результате вязкого взаимодействия отдает свой импульс окружающей жидкости. Модельная задача для такого течения должна ставиться в соответствии с теорией затопленных струй [9]. Применительно к рассматриваемым в данной статье задачам достаточно считать, что при вытеснении жидкости из ПК

Рвтек (t)|q((),о = PH . (5)

В решении (2)-(4) компоненты скорости у краев отверстия не

(ext ext \ ( ext\2

v ^ да, vz ^ да , а квадрат скорости (v ) =

r ^R™ +0 r ^R™ -0 ! V '

= (vezxt) является неинтегрируемым по площади входного сечения. Решение модельной задачи можно усовершенствовать, если распре-

ех( дфех

деление скорости =- на границе 2 = 0 считать ограниченным,

дг

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

г < Rn

f vf (r,0)dS = 2л Г ^^— (г,0)Ыг = 1. J— J dz

Решение Фех (г, г) («единичный потенциал») соответствующей задачи Неймана для уравнения Лапласа в полупространстве 2 > 0 имеет вид [39]:

ф- (г, z)=-2-г mvt (p )ds (p),

2л 0

z=0

где гмР — расстояние между точкой «интегрирования» Р(г cos ф, г sin ф,0) и точкой «наблюдения» M(г,0, z) (указаны декартовы координаты точек, ф — полярный угол цилиндрической системы координат). Учитывая характер распределения v^ и формулы для вычисления потенциала простого слоя [37, 47], для точек плоскости z = 0 имеем:

„ПК

2

ф (г,0)=-- [ — к (к 2хх{ (г,0)ёг,

к 0 г + г

где К (к2) — полный эллиптический интеграл первого рода в форме

2 4г г

Лежандра, k =--. Результаты вычислений для ряда распреде-

(г + г)2

лений уе,,х'(г,0)| кПК, включая формулы (4), приведены на рис. 4.

Кривая 4 соответствует задаче втекания воды в ПК через кольцевой зазор между корпусом ЛА и стенкой ПК. Расчеты проводились для значений КПК = 0,5, 5 = КПК - КЛА = 0,02. Потенциал

Ф^ (г, г, I) = )Ф ^ (г, г ), (6)

соответствующий этому расчету, будем использовать в дальнейшем. Присутствие корпуса ЛА учитывалось схематично условием

уе2х (г ,0) & т = 0. Для уточнения требуется серия расчетов, учитывающих выдвижение ЛА из ПК и изменение текущей скорости ЛА [37, 47, 48].

ПК

0

Расстояние от оси пускового контейнера, г/Л1™

Рис. 4. Поведение единичных потенциалов Фпри г=0:

— вычисления решения Ламба по аналитической формуле

Для внешнего течения выполняется интеграл Коши — Лагранжа

„ext i^.ext\2 ^ext

++pe dt 2

(r, Z )

= PH.

Рв

(7)

Для входного сечения КП из формул (6), (7) получаем:

|2 P

Q (t )Фext + 0,502(t )|grad Фe

В результате осреднения по сечению имеем:

Рн

2л с

РВХ(t)J Pext(r,0)rdr = Рн-«в

Рв Q

вх £) ПК ^вх

втек яЛ

R

-Рв

Рв0 2

2 (ЕКП f

у втек j

(8)

а„

2л

R

ПК КП

RTP

• J Фext (r, 0)r dr,

2л

v2 vKn

/уКП \ _ втек Rm l^BreK / втек R™

J |grad Ф^ (r,0)|2 rdr. (9)

Из формулы (6) следует, что размерность Ф^ обратна длине, поэтому коэффициенты авх, Рвх безразмерны. Результаты вычислений

этих коэффициентов по формулам (9) приведены на рис. 4.

Характер течения во входном сечении ПК, показанный на рис. 3, практически невозможно реализовать в маловязкой жидкости. В пре-

ПК

R

деле идеальной жидкости всякое безотрывное течение, обтекающее выпуклый угол, должно иметь в вершине угла бесконечную скорость [38, 42]. Решение Ламба соответствует этому случаю: при аналитическом продолжении потенциала (2) в полупространство г < 0 получится безотрывное обтекание жидкостью выпуклого угла величиной 2л.

Согласно терминологии гидравлики процесс втекания жидкости в ПК можно охарактеризовать как течение с внезапным сужением потока [49, 50]. Выделяют два основных режима таких течений: «безотрывный» (рис. 5, а) и с полным отрывом потока от стенок (рис. 5, б). В дальнейшем исследовании ограничимся рассмотрением первого режима с локальной зоной отрыва потока от стенок ПК.

Рис. 5. Схемы «безотрывного» (а) и с полным отрывом потока от стенок (б) режимов течения реальной жидкости с внезапным

сужением потока

Нестационарное течение (6) в момент времени ? = т совпадает со стационарным течением, имеющим потенциал ф^)(г,2) = б(х)х х ф<пХ (г, г) и удовлетворяющим интегралу Бернулли

'(у")2 , Р1

ext стац

(r,z)

pH

Рв

(10)

Сравнивая формулы (7) и (10), получаем:

дф^ p

at

+-

(11)

(r, z)

ext ext ext

^ г нест г г стац'

Для стационарного течения на длине входного участка имеет место потеря давления

рвхг _ ры . с рвQ

/вх. стац /вх. стац ^ вх / \ 2 *

2 (Е™ Г

\ втек /

Здесь рВХ стац — давление в месте присоединения струи к стенкам трубы при стационарном течении (см. рис. 5, а); коэффициент сопротивления входа Свх «0,5 [49].

Нестационарная составляющая давления, будучи связанной с потенциалом формулой (11), имеет более высокую степень гладкости. Ограничившись предположением р'^ нест _ Р^ нест и полагая, что

Рвтек = РвГ _ Р^ стац + К! нест , и3 формулЫ (8) имеем:

Рвтек (0 _ Рн вх ^ "Р вх ^,2, Р вх _ Рвх +С вх . (12)

К 2 (ЕКПк)

Одномерная схема расчета втекания воды в кольцевое пространство с недеформируемыми стенками. Пусть в момент X _ t0

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

^ _ Рв^втек (ОЕКПк , П (X, t) _ УвТек (t) _ V. (t) _ (13)

Здесь тв (X) _рв (р™ - хгв )ЕКП — текущая масса воды внутри КП; Ув (х, t) — скорость частиц воды в сечении с координатой хгв (X) < х < Р™ в момент времени t; ЕКП _ Е^ — площадь поперечного сечения КП.

Производная по времени от импульса Кв (X) _ тв (t )угв (t) _

_ рвvгв (РПК - хгв) ЕКП, находящегося в КП столба воды равна сумме

действующих сил и скорости переноса импульса через входное сечение

дКв__Г/ ГПК „ \ „ ,.2

- = рЕКП Г(РПК - х ) -V2

г^в I у гв у гв п

дх

_ ЕКП [ Рг (X) - Рвтек (х) ] - Рв - тв * + РСОпр, (14)

где агв _ 1>гв _ Хгв — ускорение границы воды.

Давление газа в объеме Ог 0) = £КП [ хгв - Щ) - х0 ] определяем как

( О Т

Рг О = Ро , (15)

где Ог0 = Ог (70) — начальный объем газа; п — показатель политропы. Величина ^сопр = ^С™ + ^.^Лр складывается из сил вязкого трения,

действующих на столб воды со стороны стенок ПК и ЛА. Оценим эти силы, используя формулы гидравлики для стационарных течений в трубах [49] и учитывая, что сопротивление вязкого трения определяется величиной и направлением скорости н жидкости относительно стенки. В результате получаем соотношение [30, 49]

^[х-кяпкv2 sgn^ +^лаяла ( -г)2 sgn( -г)]. (16)

Коэффициенты сопротивления трения X ПК, ЛА вычисляем из уравнения

1

Х=-

1 + Ъ^ (Яе>/Х) +c1lg А

а

Числовые значения коэффициентов а1, Ъ1, с1 [49] зависят от диапазона значений, принимаемых выражением А , где А=—А--

"гидр

относительная шероховатость стенки; А — средняя высота выступов

рв

шероховатости, м; Яе=—в-—--число Рейнольдса; — динамике

4 £КП

ческая вязкость воды; "гидр = ЦЯпк + Яла ) = 2 ( - Я" ) = 2§ш —

гидравлический диаметр КП.

Из формулы (14) получаем выражение

= (Рг - Рн )/Рв + 0, 5 РВхУг2в - (( - Хгв ) 8 + ^СОпр/(рв£КП )

агв 1ПК - Хгв -«вх £ЯПК , ( )

являющееся основой для составления системы обыкновенных дифференциальных уравнений (ОДУ) Хгв = угв , Угв = агв, легко интегрируемой

методом Рунге — Кутты. Результаты вычислений давления рг (^) для

начальных данных хгв (/0) = £ПК, Угв (¿0) = 0 и двух предельных процессов сжатия газа, изотермического (п = 1) и адиабатического (п = 1,4), приведены на рис. 6.

Время I

Рис. 6. Давление в кольцевом пространстве при втекании воды рг (1, 2), рвтек (4) без учета деформаций стенок ПК и ЛА и рг (3),

рвтек (5) с учетом деформаций; изотермический (_, — ,_,-)

и адиабатический (___) процессы сжатия газа

Результаты расчетов для случаев п = 1 и п = 1,4 сильно различаются. Реальный процесс предположительно протекает ближе к изотермическим условиям, поскольку сжатие газа в КП происходит в контакте с холодными стенками ПК и ЛА и границей воды. С одной стороны, ввиду отрыва и завихрения втекающей жидкости у ВСПК граница может быть нерегулярной и возможно проникновение в газовый объем струек и капель жидкости, что вызовет интенсификацию теплообмена. С другой стороны, стадия интенсивного сжатия газового объема кратковременна, поэтому требуется более детальное рассмотрение энергетического баланса для газового объема в КП. Причем при решении практических задач целесообразно учитывать и податливость стенок КП, поскольку, как видно из рис. 6 (кривая 3), учет деформаций стенок КП существенно влияет на динамику процесса сжатия газового объема.

Учет податливости стенок кольцевого пространства. Если учесть деформации стенок, площадь поперечного сечения КП будет определяться выражением

■^КП / г>ПК . ПК \2 / пЛА ЛА\2 л о\

£ = пук + иг ) -пук + иг ) , (18)

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

^КП ^

ж -) = 0.

в /

х в , (19)

В рамках упрощенного подхода, мотивированного работой [28], будем учитывать податливость стенок КП, подобно методу плоских сечений [48], исходя из решения классической задачи Ламе об определении напряжений и перемещений в упругой трубе, нагруженной внутренним и внешним давлением [34, 51].

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

Р2 = 0 / р1 = рКП, а на стенке ЛА — р2 = рКП / Д = 0, где рКП — местное избыточное давление газа или воды в КП. Считая, что толщины стенок малы по сравнению с их радиусами, в соответствии с решением задачи Ламе имеем:

ПК

Р КП Я

ПК

(я™ Г-(ЯПКУ

у внеш у у у

кПК + 2дПК

(ЯПК )2 (ЯеПКш )2

кПК + ДПК

Д

РКП (ЯПК )2

' 45ПК

РКП (ЯПК)2 1-(^ПК)2

ПК+ДПК

и

ЛА

(к

РКП Я ЛА

)ДПК

45г

Е П

(20)

(ЯЛА )2 -(ЯЛА )2

х кЛА + 2Д ЛА Х(к ЛА + д ЛА )дЛ

(ЯЛА )2 ((ЛА )2

кЛА+ДЛА

Д

РКП (ЯЛА )2

" 4§ла

РКП (ЯЛА)2 1-(ЛА)2

45

ЛА

Е

ЛА

где к, д — упругие константы Ламе; Е — модуль Юнга; V — коэффициент Пуассона.

Руководствуясь приведенными формулами, положим:

ПК

= к

ПК

РКП,

ЛА __ т ЛА — — _

иг = к РКП, РКП =

Р КП

Г

2

2

атм

, „ _ р„ (ДПК )21-(у ПК )2 ла _ Р„ (ЛА ) '-(у т )2 (21)

4§пк Е ПК ' 45 "А ЕЛА '

Для многих материалов у_0,25 -0,30, что дает у2 «1, тогда соотношения (20) соответствуют формуле в работе [28]. Однако в широком диапазоне значений у такое предположение может оказаться неточным [34, 52], поэтому будем придерживаться соотношений (20), (21).

Из уравнения баланса импульса для КП между жесткими стенками постоянного радиуса (см. рис. 3) следует, что дрв/дх не зависит от х, т. е. в каждый момент времени давление вдоль столба воды изменяется линейно. При учете упругости стенок КП сохраним эту зависимость в качестве аппроксимации. Тогда, согласно формулам (21), радиальные координаты стенок ПК и ЛА также будут изменяться по линейному закону, т. е. в каждый момент времени столб воды будет занимать кольцевой объем между усеченными конусами величиной

Ов (( + а2Ъ2 + Ъ22 - а2 - сЬ - Ь ) = Лв®0, (22)

где

¿в _ ¿ПК - Хгв _-¿гв ,

«1,2 _ ДЛА + <ПК (Хгв, 0 _ ДЛА+С Рг, (23)

7 пЛА ЛА / 7-ПК Л г,ЛА_7ЛА —

ъ1,2 _ ДПК + иг,ПК ( , * ) = ДПК + кПК Рвтек .

Из формул (18), (20), (21) получаем:

ЕКП (х, *) _ л (ДПК + кПК ркп )2 - л (ЛА - к - ркп )2,

ЕКП _ЕКП(Хгв,0 _*(( -«12), ЕКПк _ЕКП (,*)л(Ъ2 -Ь2). (24)

Уравнение баланса массы (13) для столба воды перепишем как

^_-у е™ . (25)

1, втек втек V /

м

Затем, дифференцируя выражения (22), (23), приходим к соотношению

0^гв + ЕККПк ^втек + К «г рг,* + КвСвтек Рвтек, * _ 0 (26)

где

«г _| [кПК ( + Ъ2)+кЛА (2а + Ъ1)], Рг,* - Рг;

«втек [ ^ («2 + 2¿2 ) + кЛА («1 + 2b1 )] , p

= Р в

втек, t Jr втек '

«г Рг, + «тек Рвтек, , ='® 0 = ®0,, .

Объемный расход воды внутрь КП определяется выражением О = увтек(,^ (,). Учитывая формулы (23), соотношение (12) запишем в виде:

d(v EKI1)- b2

втек втек у

í

а„

ph - Рв

Л

- °,5P вх v

,2

вх втек

(27)

■"вх V Рв у

Из соотношения (15) при п = 1 получаем:

Рг,,^г + 2Рг,,Ог + Ог (1 + Рг ) = 0. (28)

Объем газа в КП вычисляем по модифицированной формуле

О г (,) = 2КП [ ХгВ - Щ) - Х0 ]^ЕКП К. Согласно формулам (24), имеем:

2 КП = «2 Рг^ £ ™ = «3 Рг2,+ «12 ^

«12 = 2 л (кПК«2 + к ЛА«), «3 = 2л Г(кПК )2 - (к ЛА )2

Тогда соотношение (28) принимает вид

2КП (1 + Рг )гв +ГОг + «12 К (1 + Рг )] %.,=-2Рг,,Ог,,-«4 (1 + Рг), (29)

где Ог,, = «^г Рг,, Кг,, = Угв - Г (,); «4 = «3Йг Рг2, + 2^12 Кг,, Рг,, -2™ Ж 0.

Наконец, аналогичное формуле (14) уравнение баланса импульса позволяет записать:

dt

-ZKnp -Z™ p -р v2 ZKn -р Q gsinF

г .тг втек ^ втек г^в втек втек г^в в& с

КП сопр 5

^-Рв J Vb2ш dx-рв

гПК -ь -л

(в))-J X |(вZK1 )

(30)

Используя соотношение (19), получаем:

ЛуКП 7

рв^ = РПКу £КПП - х V £КП + Г х — —х=РПКV Е™ + — Г х£КП—х.

г в в втек втек гв гв г I ^ втек втек » I

dt

dt

Поскольку Пв = | ЕКПшх, то | хЕКПдх = | (х - хгв )ЕКПдх + хгвПв. Ин-

хгв хгв хгв

ьш

теграл | (х - хгв )ЕКПдх — это статический момент инерции усечен-

хгв

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

¿ПК 2

| (х- хгв)ЕКПдх = ^[>, + Ъ )2 + 2Ь2 ]^=2 = И2®,

хгв

Следовательно

Р-Х = ^втек ЕКППк + ^гв Пв + д (¿вЧ),

ш

1 шк ,-,КП , ш / ,-,КП \ „ шп ш2 ,, 2

в1 =_ Угв ^Е-к + ¿в щ; (втек Е^ ) + П ^ + ^ ^ + (ИЧ) :

= п V -2у У ЕКП + 2

в гв гв втек втек

Ъ2И

Г Рн - Р Л ш2

авх ^ Рв У

&

Вычислим последнее слагаемое:

ш (^в2®1) = - 2ИвУгв® 1 + ¿в2® 1,; , С® 1 = ® 1,; = агрг; + ¿втекРвтек, ; ,

ш

ат =26 [кПК (¿2+Ъ2)+кЛА (¿1 + Ъ1) ], ¿гвтек [ кПК (¿2+3Ъ2)+кЛА (¿1 + 3Ъ1)],

¿г = к(рг; + рвтек, ; ), 4тек = к(рг; + 3рвтек, ; ), к = 12-1 ^ ® 1 = ¿грг; + ¿¿втекрвтек,; + 4рг,; + 4текрвтек, ; = ¿¿грг; + ¿втекрвтек, ; + ®2,

®2 = к [(( + рвтек, ; ) ) + (( + 3рвтек,;) рвтек,;] = к (( + рвтек,; )2 + 2рв2тек,;

ш2 _ _

—- (^в2®1) = 2У>1 - 2Ив®1Угв - 4ИвУгв®1,; + № рг; + ¿¿втек ,?втек, ; +®2).

В результате подстановки в формулу (30) получим: (П - 2к®1 )у + ё И2 р . + ё И2 р . =

V в в 1/ гв г в ^гвтек в ^втек,;

ПК

у™ п _уКП п рКП ЬЬ (п _ п ~ Л

^г /Л ^втек/Уек О О , сопР и2г 1в Ун /Уек л -,,2 --"в6 81П --0,5РвхУ]

вх втек

Рв Рв авх ^ Рв У

+ 0(2^гв _ ^втек ) _ 2У>! + 4ЬвУгв©и _ %©2- (31)

Для вычисления силы сопротивления оставляем прежнюю формулу (16).

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

шаге интегрирования системы ОДУ. Дифференцируя соотношение (26) еще раз, имеем:

_ © 01>гв + ЬвагРгЛ + Ьв«втекРвтек, г = _ 4кЬв (Р1л + Рг,гРвтек, г + Р^тек, г ) +

+ 2у.

в©0,г _ Ь2/«вх (р_1 (Рн _ Рвтек ) _ 0 5Рвх^тек ) , (32)

поскольку

««гРгЛ + «втекРвтекл = 2к [(2Рг,г + РвтекД )рг,г + (Рг,г + 2рвтек,г )Рвтек, г ] = = 4кк(Р1< + РгдРвтекЛ + Рвтек, г )-

В момент вскрытия КП (г0 = 0) вода неподвижна, но ускорение частиц и давление в ней неоднородны из-за условия несжимаемости. Из

р _ р рпк

выражений (12), (17) имеем: «гв (г0) =—-1 -КП , Рвтек (г0) = Ро = Ратм —

авхРв У

т. е. стенки ПК и ЛА деформируются непрерывно. Матрица СЛАУ, составленная для старших производных, при г = г0 вырождена, поскольку Ьв (г0) = 0. По этой причине первые несколько расчетных шагов выполнялись по уравнению (17), затем вычисления велись по уравнениям (27), (29), (31), (32), а при 0 > 0 вместо соотношения (27) использовалось условие (5). Начальные условия для ОДУ были приведены ранее. При дальнейшем расчете с использованием других формул новые начальные условия согласовывались с результатами предыдущих вычислений. Необходимо, однако, обратить внимание на соблюдение соотношений (15) и (25), из-за дополнительного дифференцирования уже не являющихся автоматическим следствием решаемых ОДУ. Интегрирование систем ОДУ проводилось численно методом Рунге — Кутты, результаты расчета параметров рг (г) и

Рвтек (г) приведены на рис. 6. Данные графиков служат подтверждением того, что учет упругости стенок КП значительно снижает прогнозируемый максимум давления от втекания воды в КП.

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

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

литература

[1] Нигматулин Р.И. Динамика многофазных сред. В 2 ч. Москва, Наука, 1987, ч. 1, 464 с.; ч. 2, 360 с.

[2] Накоряков В.Е., Покусаев Б.Г., Шрейбер И.Р. Волновая динамика газо-и парожидкостных сред. Москва, Энергоатомиздат, 1990, 248 с.

[3] Яковлев Ю.С. Гидродинамика взрыва. Ленинград, Судпромгиз, 1961, 316 с.

[4] Приходько Н.А., Сирый В.С. Расчет подъема пузырька горячего газа в вязком неравномерно нагретом слое жидкости. Гидромеханика, 1980, вып. 43, с. 85-93.

[5] Cole R.H. Underwater explosions. Princeton University Press, 1948, 437 p.

[6] Leighton T.G. The acoustic bubble. London, Academic Press, 1994, 613 p.

[7] Грумондз В.Т., Журавлев Ю.Ф., Парышев Э.В., Соколянский В.П., Шоры-гин В.П. Гидродинамика и динамика высокоскоростного движения тел в жидкости. Москва, Наука, 2013, 574 с.

[8] Рождественский В.В. Кавитация. Ленинград, Судостроение, 1977, 247 с.

[9] Абрамович Г.Н. Теория турбулентных струй. Москва, Наука, 1984, 716 с.

[10] Якимов Ю.Л. Сборник статей. Москва, Изд-во Московского университета, 2013, 536 с.

[11] Кутателадзе С.С., Накоряков В.Е. Тепломассобмен и волны в газожидкостных системах. Новосибирск, Наука, 1984, 301 с.

[12] Алемасов В.Е. Теория ракетных двигателей. Москва, Оборонгиз, 1962, 476 с.

[13] Ваулин С. Д., Кириллов В.В., Феофилактов В.И. Математическая модель газодинамических процессов в низкотемпературном газогенераторе с камерой охлаждения. Ракетно-космическая техника. Сер. XIV, 2004, вып. 1 (50), ч. I, с. 181-190.

[14] Якимов Ю.Л., Ерошин В. А., Романенков Н.И. Моделирование движения тела в воде с учетом ее сжимаемости. Некоторые вопросы механики сплошной среды. Москва, Изд-во Московского университета, 1978, с. 29-33.

[15] Lohse D., Schmitz B., Versluis M. Snapping shrimp make flashing bubbles. Nature, 2001, no. 413, pp. 477-478.

[16] Scorer R.S. Environmental aerodynamics. New York, Halsted Press, 1978, 551 p.

[17] Александров А.А., Димитриенко Ю.И. Математическое и компьютерное моделирование — основа современных инженерных наук. Приветствие

к выходу первого номера журнала «Математическое моделирование и численные методы». Математическое моделирование и численные методы,

2014, № 1 (1), с. 3-4.

[18] Perot B., Nallapati R. A moving unstructured staggered mesh method for the simulation of incompressible free-surface flows. Journal of Computational Physics, 2003, vol. 184, no. 1, pp. 192-214.

[19] Peregrine D.H. The fascination of fluid mechanics. Journal of Fluid Mechanics, May 1981, vol. 106, pp. 91-119.

[20] Jürgens H., Peitgen H.-O., Saupe D. The Language of Fractals. Scientific American, 1990, no. 263, pp. 60-67.

[21] Апальков Ю.В., Мант Д.И., Мант С. Д. Отечественные баллистические ракеты морского базирования и их носители. СПб., Галея, 2006, 216 с.

[22] Морские стратегические ракетные комплексы. Москва, Военный парад, 2011, 268 с.

[23] Дегтярь В.Г., Пегов В.И. Гидродинамика подводного старта ракет. Москва, Машиностроение, 2009, 448 с.

[24] Ефремов Г.А., Страхов А.Н., Минасбеков Д.А., Горлашкин А.А., Плюснин А.В., Соколов П.М., Бондаренко Л.А., Говоров В.В. Отработка газодинамики подводного старта на наземном газодинамическом стенде предприятия. Ракетные комплексы и ракетно-космические системы — проектирование, экспериментальная отработка, летные испытания, эксплуатация. Труды секции 22 имени академика В.Н. Челомея XXXVIII Академических чтений по космонавтике. Реутов, 2014, с. 65-74.

[25] Плюснин А.В., Бондаренко Л.А., Сабиров Ю.Р. Анализ газогидродинамических процессов и методов их расчета на основе опыта предприятия в отработке подводного минометного старта. Ракетные комплексы и ракетно-космические системы — проектирование, экспериментальная отработка, летные испытания, эксплуатация. Труды секции 22 имени академика В.Н. Челомея XXXIX Академических чтений по космонавтике. Реутов,

2015, с. 74-83.

[26] Дергачев А.А., Бондаренко Л. А., Сабиров Ю.Р., Лобзов Н.Н., Плюснин А.В. Способ старта ракет с подводной лодки, надводных кораблей и наземных носителей из незатопленной пусковой установки и пусковая установка для его реализации. Пат. Российская Федерация № 2536961, Опубл. 27.12.2014, бюл. № 36.

[27] Плюснин А.В. Моделирование параметров наддува свободного пространства контейнера при газодинамическом выбросе летательного аппарата с учетом свойств реального газа. Математическое моделирование и численные методы, 2016, № 3 (11), c. 53-78.

[28] Жуковский Н.Е. О гидравлическом ударе в водопроводных трубах. Москва, Либроком, 2011, 104 с.

[29] Плюснин А.В. Учет эффекта вторичного догорания при расчетах систем газодинамического выброса летательного аппарата. Математическое моделирование и численные методы, 2014, № 3, с. 55-73.

[30] Плюснин А.В. Исследование зависимости максимальной расчетной оценки давления «гидроудара» при подводном газодинамическом выбросе ЛА от выбора физической модели явления. Сборник тезисов XLI Академических чтений по космонавтике, посвященных памяти академика С.П. Королёва и других выдающихся отечественных ученых — пионеров освоения космического пространства. Москва, МГТУ им. Н.Э. Баумана, 2017, с. 516.

[31] Плюснин А.В., Бондаренко Л. А. Способы крупномасштабного моделирования систем газодинамического выброса. Вестник МГТУ им. Н. Э. Баума-

на. Сер. Естественные науки, 2012. Спец. выпуск № 4. «Математическое моделирование», с. 111-122.

[32] Зарубин В.С., Кувыркин Г.Н. Особенности математического моделирования технических устройств. Математическое моделирование и численные методы, 2014, № 1, с. 5-17.

[33] Григолюк Э.И., Горшков А.Г. Взаимодействие упругих конструкций с жидкостью (удар и погружение). Ленинград, Судостроение, 1976, 200 с.

[34] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 4. Основы механики твердых сред. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2013, 624 с.

[35] Ерошин В.А., Плюснин А.В., Романенков Н.И., Созоненко Ю.А., Якимов Ю.Л. О влиянии атмосферы на величину гидродинамических сил при плоском ударе диска о поверхность сжимаемой жидкости. Известия Академии наук СССР. Механика жидкости и газа, 1984, № 3, с. 15-20.

[36] Димитриенко Ю.И. Механика сплошной среды. В 4 т. Т. 2. Универсальные законы механики и электродинамики сплошных сред. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2011, 560 с.

[37] Плюснин А.В. Расчет нестационарной гидродинамической нагрузки на крышку пускового устройства при подводном выбросе ЛА избыточным давлением газа. Инженерный журнал: наука и инновации, 2013, вып. 7. URL: http://engjournal.ru/ catalog/hmodel/aero/847.html (дата обращения 03.10.2017).

[38] Лойцянский Л.Г. Механика жидкости и газа. Москва, Наука, 1987, 840 с.

[39] Соболев С.Л. Уравнения математической физики. Москва, Наука, 1966, 444 с.

[40] Седов Л.И. Механика сплошной среды. В 2 т. Т. 2. Санкт-Петербург, Лань, 2004, 560 с.

[41] Lamb H. Hydrodynamics. Cambridge, Cambridge University Press, 1932, 928 p.

[42] Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. Москва, Наука, 1987, 688 с.

[43] Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. В 3 т. Т. 3. Специальные функции. Москва, Наука, 1983, 752 с.

[44] Плюснин А.В., Бондаренко Л.А., Сабиров Ю.Р. Расчет нестационарной гидродинамической нагрузки, действующей на деформируемую панель оперения ЛА. Ракетные комплексы и ракетно-космические системы — проектирование, экспериментальная отработка, летные испытания, эксплуатация. Труды секции 22 имени академика В.Н. Челомея XL Академических чтений по космонавтике. Реутов, 2016, с. 185-196.

[45] Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. В 3 т. Т. 2. Функции Бесселя, функции параболического цилиндра, ортогональные многочлены. Москва, Наука, 1974, 296 с.

[46] Batchelor G.K. An introduction to fluid dynamics. Cambridge, Cambridge University Press, 2000, 615 p.

[47] Плюснин А.В. Расчет коэффициентов внутреннего и внешнего нестационарного взаимодействия корпуса ЛА с жидкостью методом граничных элементов. Математическое моделирование и численные методы, 2014, вып. 2, с. 77-100.

[48] Плюснин А.В., Доденко И.А. Исследование точности метода распределенной присоединенной массы при расчете нестационарной поперечной нагрузки на деформируемый корпус ЛА при подводном выбросе. Инженерный журнал: наука и инновации, 2013, вып. 7.

URL: http://engjournal.ru/catalog/mathmodel/aero/841.html (дата обращения 03.10.2017).

[49] Идельчик И.Е. Справочник по гидравлическим сопротивлениям. Москва, Машиностроение, 1992, 672 с.

[50] Сиов Б.Н. Истечение жидкости через насадки. Москва, Машиностроение, 1968, 140 с.

[51] Тимошенко С.П., Гудьер Дж. Теория упругости. Москва, Наука, 1979, 560 с.

[52] Полилов А.Н. Этюды по механике композитов. Москва, Физматлит, 2016, 320 с.

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

Ссылку на эту статью просим оформлять следующим образом: Плюснин А.В. Математическое моделирование процесса втекания воды в кольцевое прстранство контейнера при подводном газодинамическом выбросе летательного аппарата. Математическое моделирование и численные методы, 2017, № 2, с. 39-64.

Плюснин Андрей Владимирович окончил Днепропетровский государственный университет по специальности «Гидроаэродинамика» и аспирантуру кафедры «Волновая и газовая динамика» МГУ им. М.В. Ломоносова. Канд. физ.-мат. наук, доцент-совместитель кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана, заместитель начальника отдела АО «ВПК «НПО машиностроения». Автор порядка 30 работ по нестационарным задачам гидрогазодинамики и теории упругости. e-mail: andrey.plusnin@gmail.com

Mathematical simulation of the process of water entering the annular space of a canister during submarine gas-driven aircraft ejection

© A.V. Plyusnin Bauman Moscow State Technical University, Moscow, 105005, Russia

The study deals with a one-dimensional analytical model for computing the loads on the body of an aircraft caused by water entering the annular space of a launch canister. We used potential theory to solve the "external" hydrodynamic problem. Solving Lamé equations for the static case accounts for the strain in the walls of the aircraft and the launch canister.

Keywords: aircraft, submarine ejection, annular space, water influx

references

[1] Nigmatulin R.I. Dinamika mnogofaznykh sred [Dynamics of multi-phase media]. In 2 vols. Moscow, Nauka Publ., 1987, vol. 1, 464 p.; vol. 2, 360 p.

[2] Nakoryakov V.E., Pokusaev B.G., Shreyber I.R. Volnovaya dinamika gazo-i parozhidkostnykh sred [Wave dynamics of gas-liquid and vapour-liquid media]. Moscow, Energoatomizdat Publ., 1990, 248 p.

[3] Yakovlev Yu.S. Gidrodinamika vzryva [Hydrodynamics of explosion]. Leningrad, State Publishing House for Shipbuilding Literature, 1961, 316 p.

[4] Prikhodko N.A., Siryy V.S. Gidromekhanika — Hydromechanics, 1980, no. 43, pp. 85-93.

[5] Cole R.H. Underwater explosions. Princeton University Press, 1948, 437 p.

[6] Leighton T.G. The acoustic bubble. London, Academic Press, 1994, 613 p.

[7] Grumondz V.T., Zhuravlev Yu.F., Paryshev E.V., Sokolyanskiy V.P., Shorygin V.P. Gidrodinamika i dinamika vysokoskorostnogo dvizheniya tel v zhidkosti [Hydrodynamics and dynamics of high-speed motion of bodies in fluid]. Moscow, Nauka Publ., 2013, 574 p.

[8] Rozhdestvenskiy V.V. Kavitatsiya [Cavitation]. Leningrad, Sudostroenie Publ., 1977, 247 p.

[9] Abramovich G.N. Teoriya turbulentnykh struy [Theory of turbulent jets]. Moscow, Nauka Publ., 1984, 716 p.

[10] Yakimov Yu.L. Sbornik statey [Collected publications]. Moscow, Moscow State University Publ., 2013, 536 p.

[11] Kutateladze S.S., Nakoryakov V.E. Teplomassobmen i volny v gazozhidkostnykh sistemakh [Heat and mass transfer and waves in gas-liquid systems]. Novosibirsk, Nauka Publ., 1984, 301 p.

[12] Alemasov V.E. Teoriya raketnykh dvigateley [Theory of rocket engines]. Moscow, State Publishing House for the Defense Industry, 1962, 476 p.

[13] Vaulin S.D., Kirillov V.V., Feofilaktov V.I. Matematicheskaya model gazo-dinamicheskikh protsessov v nizkotemperaturnom gazogeneratore s kameroy okhlazhdeniya [Mathematical model of gas-dynamic processes in a low-temperature gas generator with a cooling chamber]. Raketno-kosmicheskaya tekhnika. Ser. XIV [Aerospace technology. Series 14], 2004, no. 1 (50), vol. 1, pp. 181-190.

[14] Yakimov Yu.L., Eroshin V.A., Romanenkov N.I. Modelirovanie dvizheniya tela v vode s uchetom ee szhimaemosti [Simulating motion of a body in water, considering its compressibility]. Nekotorye voprosy mekhaniki sploshnoy sredy [Certain problems in continuum mechanics]. Moscow, Moscow State University Publ., 1978, pp. 29-33.

[15] Lohse D., Schmitz B., Versluis M. Nature, 2001, no. 413, pp. 477-478.

[16] Scorer R.S. Environmental aerodynamics. New York, Halsted Press, 1978, 551 p.

[17] Aleksandrov A.A., Dimitrienko Yu.I. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 1 (1), pp. 3-4.

[18] Perot B., Nallapati R. Journal of Computational Physics, 2003, vol. 184, no. 1, pp. 192-214.

[19] Peregrine D.H. Journal ofFluid Mechanics, May 1981, vol. 106, pp. 91-119.

[20] Jürgens H., Peitgen H.-O., Saupe D. Scientific American, 1990, no. 263, pp. 60-67.

[21] Apalkov Yu.V., Mant D.I., Mant S.D. Otechestvennyye ballisticheskiye rakety morskogo bazirovaniya i ikh nositeli [National sea-launched ballistic missiles and their carriers]. St. Petersburg, Galeya Print Publ., 2006, 216 p.

[22] Morskiye strategicheskiye raketnye kompleksy [Marine strategic missile system]. Moscow, Voyennyy parad Publ., 2011, 268 p.

[23] Degtyar V.G., Pegov V.I. Gidrodinamika podvodnogo starta raket [Hydrodynamics of the submarine missile takeoff]. Moscow, Mashinostroeniye Publ., 2009, 448 p.

[24] Efremov G.A., Strakhov A.N., Minasbekov D.A., Gorlashkin A.A., Plyusnin A.V., Sokolov P.M., Bondarenko L.A., Govorov V.V. Otrabotka gazodinamiki podvodnogo starta na nazemnom gazodinamicheskom stende predpriyatiya [Underwater launch gas dynamic tests on the ground gas dynamics ground test bench of the enterprise]. Raketnye kompleksy i raketno-kosmicheskie sistemy — proektirovanie, eksperimentalnaya otrabotka, letnye ispytaniya, ekspluatatsiya. Trudy sektsii 22 im. akad. V.N. Chelomeya XXXVIII Akademicheskikh chteniy po kosmonavtike [Rocket and space-rocket systems — designing, experimental

A.B. nmcnuH

tests, flight tests, exploitation. Proc. of the 22nd section acad. V.N. Chelomei of

XXXVIII Academic Space Technology Readings]. Reutov, 2014, pp. 65-74.

[25] Plyusnin A. V., Bondarenko L.A., Sabirov Yu.R. Analiz gazogidrodinamich-eskikh protsessov i metodov ikh rascheta na osnove opyta predpriyatiya v otra-botke podvodnogo minometnogo starta [Analysis of gas and hydro dynamic processes and their calculating methods on the basis of the enterprise experiments in underwater mortar launch tests]. Raketnyye kompleksy i raketno-kosmicheskiye sistemy — proyektirovaniye, eksperimentalnaya otrabotka, letnyye ispytaniya, ekspluatatsiya. Trudy sektsii 22 im. akad. V.N. Chelomeya

XXXIX Akademicheskikh chteniy po kosmonavtike [Rocket and space-rocket systems — designing, experimental tests, flight tests, exploitation. Proceedings of the 22nd section named after acad. V.N. Chelomei of the XXXIX Academic Space Technology Readings]. Reutov, 2015, pp. 74-83.

[26] Dergachev A.A., Bondarenko L.A., Sabirov Yu.R., Lobzov N.N., Plyusnin A.V. Sposob starta raket s podvodnoy lodki, nadvodnykh korabley i nazemnykh nos-iteley iz nezatoplennoy puskovoy ustanovki i puskovaya ustanovka dlya ego real-izatsii [A way for launching rockets from a submarine, surface combatants and ground vehicles using a non-submerged launcher and a launcher to implement it]. Patent RU2536961, bulletin no. 36, Dec 27th, 2014.

[27] Plyusnin A.V. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2016, no. 3 (11), pp. 53-78.

[28] Zhukovskiy N.E. O gidravlicheskom udare v vodoprovodnykh trubakh [On water hammer in water plumbing pipes]. Moscow, LIBROKOM Publ., 2011, 104 p.

[29] Plyusnin A.V. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 3 (3), pp. 55-73.

[30] Plyusnin A.V. Issledovanie zavisimosti maksimalnoy raschetnoy otsenki davleni-ya "gidroudara" pri podvodnom gazodinamicheskom vybrose LA ot vybora fizi-cheskoy modeli yavleniya [Investigating how the maximum estimated value of fluid hammer pressure during submarine gas-driven aircraft ejection depends on the physical model selected to describe the phenomenon]. Sbornik tezisov XLI Akademicheskih chteniy po kosmonavtike, posvycshennykh pamyati akademika V.N. Chelomeya i drugikh vydayushcihsya otechestvennykh uchenykh — pionerov osvoeniya kosmicheskogo prostranstva [Abstracts XLI Academic readings on cosmonautics, dedicated to the memory of academician S.P. Korolev and other outstanding domestic scientists — pioneers of space exploration]. Moscow, Bau-man Moscow State Technical University, 2017, p. 516.

[31] Plyusnin A.V., Bondarenko L.A. Vestnik MGTU im. N.E. Baumana. Ser. Yestestvennyye nauki. Spets. vypusk "Matematicheskoye modelirovaniye" — Herald of Bauman Moscow State Technical University, Series Natural Sciences. Special iss. "MathematicalModeling", 2012, no. 4, pp. 111-122.

[32] Zarubin V.S., Kuvyrkin G.N. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 1, pp. 5-17.

[33] Grigolyuk E.I., Gorshkov A.G. Vzaimodeystvie uprugikh konstruktsiy s zhidkostyu (udar i pogruzhenie) [Elastic structure interaction with fluids (impact and submerging)]. Leningrad, Sudostroenie Publ., 1976, 200 p.

[34] Dimitrienko Yu.I. Mekhanika sploshnoy sredy. V 4 tomakh. Tom 4. Osnovy mekaniki tverdykh sred [Continuum mechanics. In 4 vols. Vol. 4. Fundamentals of mechanics of solid mechanics]. Moscow, BMSTU Publ., 2013, 624 p.

[35] Eroshin V.A., Plyusnin A.V., Romanenkov N.I., Sozonenko Yu.A., Yakimov Yu.L. Izvestiya Akademii nauk SSSR. Mekhanika zhidkosti i gaza — Fluid Dynamics, 1984, no. 3, pp. 15-20.

[36] Dimitrienko Yu.I. Mekhanika sploshnoy sredy V 4 tomakh. Tom 2. Universalnye zakony mekhaniki i elektrodinamiki sploshnykh sred [Mechanics of Continua. In 4 vols. Vol. 2. Universal Laws of Mechanics and Electrodynamics of Continuous Media]. Moscow, BMSTU Publ., 2011, 560 p.

[37] Plyusnin A.V. Inzhenernyy zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovations, 2013, no. 7 (19), Available at: http://engjournal.ru/ catalog/hmodel/aero/847.html (дата обращения 03.10.2017).

[38] Loytsyanskiy L.G. Mekhanika zhidkosti i gaza [Fluid mechanics]. Moscow, Nauka Publ., 1987, 840 p.

[39] Sobolev S.L. Uravneniya matematicheskoy fiziki [Equations of mathematical physics]. Moscow, Nauka Publ., 1966, 444 p.

[40] Sedov L.I. Mekhanika sploshnoy sredy [Continuum mechanics]. In 2 vols. Vol. 2. St. Petersburg, Lan Publ., 2004, 560 p.

[41] Lamb H. Hydrodynamics. Cambridge, Cambridge University Press, 1932, 928 p.

[42] Lavrentyev M.A., Shabat B.V. Metody teorii funktsii kompleksnogo peremen-nogo [Complex variable theory methods]. Moscow, Nauka Publ., 1987, 688 p.

[43] Prudnikov A.P., Brychkov Yu.A., Marichev O.I. Integraly i ryady [Integrals and series]. In 3 vols. Vol. 3. Spetsialnye funktsii [Special functions]. Moscow, Nauka Publ., 1983, 752 p.

[44] Plyusnin A.V., Bondarenko L.A., Sabirov Yu.R. Raschet nestatsionarnoy gidro-dinamicheskoy nagruzki, deystvuyushchey na deformiruemuyu panel opereniya LA [Computing unsteady hydrodynamic loads on a deformable empennage panel]. Raketnye kompleksy i raketno-kosmicheskie sistemy — proektirovanie, ek-sperimentalnaya otrabotka, letnye ispytaniya, ekspluatatsiya: Trudy sektsii 22 imeni akademika V.N Chelomeya XL Akademicheskikh chteniy po kosmonavtike [Rocket and aerospace systems — design, experiment-based refinement, flight testing, operation: Proc. of the section 22 of V.N. Chelomey 40th Academic readings in cosmonautics]. Reutov, 2016, pp. 185-196.

[45] Bateman H., Erdelyi A. Higher Transcendental Functions. In 3 vols. Vol. 2. McGraw-Hill Book Company, 1953, 396 p. [In Russ.: Bateman H., Erdely A. Vysshie transtsendentnye funktsii. In 3 vols. Vol. 2. Funktsii Besselya, funktsii parabolicheskogo tsilindra, ortogonalnye mnogochleny [Bessel functions, functions of the parabolic cylinder, orthogonal polynomials]. Moscow, Nauka Publ., 1974, 296 p.].

[46] Batchelor G.K. An introduction to fluid dynamics. Cambridge, Cambridge University Press, 2000, 615 p.

[47] Plyusnin A.V. Matematicheskoe modelirovanie i chislennye metody — Mathematical Modeling and Computational Methods, 2014, no. 2, pp. 77-100.

[48] Plyusnin A.V., Dodenko I.A. Inzhenernyy zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2013, no. 7. Available at: http://engjournal.ru/catalog/mathmodel/aero/841.html (accessed October 3, 2017).

[49] Idelchik I.E. Spravochnik po gidravlicheskim soprotivleniyam [Handbook on friction loss values]. Moscow, Mashinostroenie Publ., 1992, 672 p.

[50] Siov B.N. Istechenie zhidkosti cherez nasadki [Fluid flow through attachments]. Moscow, Mashinostroenie Publ., 1968, 140 p.

[51] Timoshenko S.P., Goodier J.N. Theory of Elasticity. McGraw-Hill Education, 1970, 608 p. [In Russ.: Timoshenko S.P., Goodier J.N. Teoriya uprugosti. Moscow, Nauka Publ., 1979, 560 p.].

[52] Polilov A.N. Etyudy po mekhanike kompozitov [Studies in composite mechanics]. Moscow, FIZMATLIT Publ., 2016, 320 p.

Plyusnin A. V., graduated from Dnepropetrovsk State University, majoring in Aerohy-drodynamics and obtained a degree in the Department of Wave and Gas Dynamics, Lo-monosov Moscow State University. Cand. Sci. (Phys. & Math.), Assoc. Professor of the Computational Mathematics and Mathematical Physics Department at the Bauman Moscow State Technical University, Deputy Director of the Department in at JSC "MIC "Mashinostroeniya". Author of approximately 30 publications in the field of nonstation-ary problems of hydrodynamics, gasdynamics and theory of elasticity. e-mail: andrey.plusnin@gmail.com