DOI: 10.15593/2224-9982/2021.64.10 УДК 519.9, 629.7
М.Ю. Егоров1, Д.М. Егоров2, С.М. Егоров2
1Пермский национальный исследовательский политехнический университет, Пермь, Россия Научно-исследовательский институт полимерных материалов, Пермь, Россия
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИКИ ВНУТРИКАМЕРНЫХ ПРОЦЕССОВ В МАРШЕВОМ РАКЕТНОМ ДВИГАТЕЛЕ НА ТВЕРДОМ ТОПЛИВЕ С УЧЕТОМ ПОЛЕТНЫХ ПЕРЕГРУЗОК. ЧАСТЬ 1. МЕТОДИКА РАСЧЕТА
Исследуется динамика переходных внутрикамерных процессов (внутренняя баллистика) маршевого РДТТ второй ступени крылатой ракеты с учетом в общем случае распределенных пространственно-трехмерных и изменяющихся во времени полетных перегрузок. Метод исследования - постановка вычислительного эксперимента. Рассматривается сопряженная постановка задачи, включающая в себя: нестационарное срабатывание воспламенительного устройства (скорость горения воспламенительного состава описывается на основе экспериментально-теоретического подхода, учитывающего догорание продуктов сгорания за корпусом воспламенительного устройства); прогрев, зажигание и последующее нестационарное и турбулентное горение заряда твердого топлива (используется квазигомогенная модель горения на основе уравнений теплопроводности и химической кинетики, записанных для конденсированной фазы (твердого топлива), с учетом влияния газовой фазы (факела) на процесс горения в конденсированной фазе; метод решения задачи - метод конечных разностей); нестационарное трехмерное гомогенно-гетерогенное четырехфазное течение воздуха и продуктов сгорания в камере сгорания, сопловом блоке и за сопловым блоком ракетного двигателя (используются подходы механики сплошных многофазных сред; базовая система уравнений - система вихревых дифференциальных уравнений газовой динамики; метод решения - многопараметрический класс разностных схем расщепления по физическим процессам метода Давыдова); разгерметизацию камеры сгорания РДТТ (уравнение движения заглушки соплового блока - второй закон Ньютона; метод решения задачи - метод Эйлера). Каждая из подзадач рассматривается во взаимосвязи и разрешается одновременно - на одном шаге по времени. Для решения сформулированной задачи разработан комплекс прикладных программ с использованием (для основного расчетного модуля) стандарта многопотоковой обработки информации OpenCL. Проведена проверка работоспособности программного продукта.
Ключевые слова: численное исследование, ракетный двигатель твердого топлива, горение твердого топлива, газовая динамика, полетные перегрузки, методика расчета.
M.Yu. Egorov 1, D.M. Egorov 2, S.M. Egorov 2
Perm National Research Polytechnic University, Perm, Russian Federation Research Institute of Polymer Materials, Perm, Russian Federation
NUMERICAL STUDY OF DYNAMICS INTRACHAMBER PROCESSES IN SOLID PROPELLANT SUSTAINER TAKING INTO ACCOUNT FLIGHT OVERLOADS.
PART 1. CALCULATION METHOD
The dynamics of transient in-chamber processes (internal ballistics) of the cruise missile's second-stage cruise missile propulsion system is studied, taking into account, in the general case, distributed spatially-three-dimensional and time-varying flight overloads. The research method is the formulation of a computational experiment. Be considered coupled formulation of the problem, including: - transient triggering of igniter device (the rate of combustion of the igniting composition is described on the basis of experimental and theoretical approach afterburn combustion products in case igniting device); preheating, ignition and subsequent unsteady and turbulent combustion of solid propellant charge (used quasi-homogeneous combustion model based on the equations of heat conduction and chemical kinetics recorded for a condensed phase (solid fuel), taking into account the influence of the gas phase (torch) on the process of combustion in the condensed phase; the method of solving the problem - finite difference method); - non-stationary three-dimensional homogeneous-heterogeneous four phase flow of air and products of combustion in the combustion chamber, the nozzle block and the block launchers rocket engine (used approaches of continuum mechanics of multiphase media; the basic system of equations system of vortex differential equations of gas dynamics solution method - a multi-parameter class of difference schemes of splitting into physical processes of the method Davydova); - depressurization of the combustion chamber of the SRB (equation of motion of the plug nozzle block - Newton's second law; the proposed solution method - Euler's
method). Each of the subtasks is considered in a relationship and resolved simultaneously - at one time step. To solve the formulated problem, a set of application programs has been developed using (for the main calculation module) the OpenCL multithreaded information processing standard. The performance of the software product was checked.
Keywords: numerical study, solid fuel rocket engine, burning of solid propellant, gas dynamics, flight overloads, calculation method.
Введение
Численное исследование переходных внутрикамерных процессов (так называемой внутренней баллистики) производится на этапе проектирования РДТТ [1-7]. Исследуются особенности процесса срабатывания воспламени-тельного устройства, прогрев, зажигание и горение основного заряда твердого топлива, газовая динамика в камере сгорания, в сопловом блоке, за сопловым блоком и другие вспомогательные задачи. Далее осуществляется конструирование ракетного двигателя, изготовление опытной партии РДТТ и его стендовая отработка. Результаты численных расчетов сопоставляются с данными стендовых испытаний РДТТ. После этого разработанное изделие (ракетный двигатель) в составе ракетного комплекса переходит на этап летных испытаний. Очевидно, что условия летных испытаний не в полной мере соответствуют условиям проектирования и стендовой отработки РДТТ, особенно тогда когда величина планируемых полетных перегрузок значительная (более 10g). В этом случае у разработчиков ракетной техники возникают вопросы, связанные как минимум с оценкой влияния перегрузок и как максимум с полным анализом влияния внешней аэродинамики и поля силы тяжести на внутри-камерные процессы и, соответственно, на выходные характеристики РДТТ (расход продуктов сгорания, тяга двигателя, время релаксации камеры сгорания и пр.).
Ранее в ряде работ одного из авторов статьи проводились исследования влияния осевой полетной перегрузки на динамику волнового процесса в камере сгорания РДТТ [8, 9 и др.].
Было установлено, что значительные осевые перегрузки влияют на амплитуду колебаний давления (по первой продольной моде) в камере сгорания. По мере увеличения уровня перегрузки наблюдается уменьшение амплитуды колебаний давления в камере сгорания и амплитуды колебаний тяги ракетного двигателя, хотя общий характер силового нагружения ракетной системы при перегрузках усложняется.
Предлагается продолжить исследование в этом важном для теории и практики направлении. Используя современные возможности вычислительной математики (в первую очередь численных методов), методов алгоритмирова-ния и программирования, методов обработки информации и высокопроизводительных вычислительных систем, можно перейти на более глубокий уровень понимания динамики рассматриваемых процессов.
Для этого в данной работе исследуется динамика переходных внутрикамерных процессов (внутренняя баллистика) маршевого РДТТ второй ступени крылатой ракеты с учетом распределенных пространственно-трехмерных и изменяющихся во времени полетных перегрузок. Принципиальная компоновочная схема ракетного двигателя показана на рис. 1.
Метод исследования - постановка прямого вычислительного эксперимента (основной расчетный метод - метод Давыдова). Рассматривается сопряженная постановка задачи, включающая в себя: нестационарное срабатывание воспламенительного устройства; прогрев, зажигание и последующее нестационарное и турбулентное горение заряда твердого топлива; нестационарное трехмерное гомогенно-гетерогенное четырехфазное течение
Рис. 1. Компоновочная схема РДТТ
воздуха и продуктов сгорания в камере сгорания, сопловом блоке и за сопловым блоком ракетного двигателя; разгерметизацию камеры сгорания РДТТ (вскрытие заглушки соплового блока). Каждая из подзадач рассматривается во взаимосвязи и разрешается одновременно -на одном шаге по времени.
Срабатывание воспламенительного устройства
Процесс срабатывания воспламенительного устройства (ВУ) с жестким несгораемым перфорированным корпусом и вкладным формованным зарядом пиротехнического воспла-менительного состава описывается на основе экспериментально-теоретического подхода. Путем решения обратной задачи внутренней баллистики для системы «ВУ - имитатор свободного объема камеры сгорания РДТТ» рассчитывается реальный газоприход от ВУ с учетом особенности горения заряда пиротехнического воспламенительного состава - догорания продуктов сгорания за корпусом ВУ. Скорость горения заряда ВУ ув напрямую связывается
с параметрами, определяющими процесс горения в корпусе ВУ (площадью перфорации корпуса ВУ площадью поверхности горения
заряда ВУ 5в и свободным объемом в корпусе
ВУ Ув), в виде зависимости ув = /(^кр,Бв,¥в).
Такого рода связи, не вникая в сущность сложного физико-химического процесса горения, можно найти в статистическом подходе, используя полиномиальные модели вида
V = А0 + А1Х1 + А2 X 2 + А3 XIX 2 +
+А4Х12 + А5 X 22 + А6 XIX 22
(1)
X1 =
Г F
- В. V 1
/ В2; X 2 =
^ - с
V. 1
/ с2
По полученной скорости горения заряда пиротехнического воспламенительного состава (1) определяется газоприход от ВУ в камеру сгорания РДТТ.
Прогрев, зажигание и горение заряда твердого топлива
Описание процесса нестационарного прогрева, воспламенения и последующего нестационарного и турбулентного горения заряда смесевого твердого топлива базируется на модели Мержанова - Дубовицкого-Соркина с учетом, в рамках подхода Горохова - Липанова -Русяка, влияния газовой фазы на процесс горения в конденсированной фазе (к-фазе) [1, 3, 4]. Используется так называемая квазигомогенная модель горения твердого топлива.
Будем рассматривать твердое топливо как твердое тело, к которому применимы известные уравнения теплопроводности и химической кинетики. Для удобства будем рассматривать эти уравнения в системе координат, связанной с поверхностью горения, направив ось от поверхности вглубь пороха. Считаем, что реакции в к-фазе удовлетворяют закону Арре-ниуса. Тогда в предположении 0-мерности порядка химических реакций данная система уравнений, описывающая процесс, имеет вид
дt
д 2Тк к я 2 дУк +к дУк
др ¥=ук др дУк .+Ф
Qk
(2)
где ф к (Тк ) = ^ ехР
V К0Тк у
а для нахождения коэффициентов полинома применить теорию планирования эксперимента. Здесь А0,...,А6,В1,В2,С1,С2 - коэффициенты полинома, зависящие от конкретного вида пиротехнического воспламенительного состава и интервалов варьирования конструктивных параметров ВУ. Более подробно методика расчета изложена в работе [10] и других работах.
До воспламенения необходимо положить Ук = 0. Условие воспламенения и последующего горения: Р^ = Р» = 1.
До воспламенения начальные и граничные условия для системы (2) имеют вид
t = 0, >к > 0, Тк = То, Р = 0;
дТ
t > 0, ^ = 0, -Т = ах( - Т,), Р<Р,; (3) дУк V У
t > 0, Ук = ®, Тк = Т0, р = 0.
После воспламенения граничные условия запишутся в виде
t > t„ yk = 0,
Tk = T, Р, = p.,
л t=q -Ck )pkvkT-[a -t bp;
(4)
t > t„ yk =<»,
Tk = To,
p = o.
Задача (2)-(4), дополненная замыкающими соотношениями, обезразмеренная и записанная в неравномерной (экспоненциальной) системе координат, интегрируется численным сеточным методом по явным и неявным конечно-разностным схемам [4, 11]. Неявные схемы разрешаются способом прогонки.
Температура по своду к-фазы, исключая поверхность горения пороха,
2 < I < (М -1). (5)
T/+1 = AT ++1 + Bt;
Температура на поверхности горения пороха:
- до воспламенения
Tj+1 = AjT/+1 + Bt;
- после воспламенения
T j+1 = A1.T2j+1 + B*.
Глубина превращения пороха:
- до воспламенения
(6)
(7)
Рj+1 = pj + Ata3Tj, 1 < i < (M -1); (8)
- после воспламенения
Р/+1 = Ai«Pi+1 + B,,„ 1 < i < (M -1).
(9)
Скорость горения порохового заряда vk определяется из формулы (9) с учетом граничных условий Р/+1 = Р1 = Р„, РM1 = pM = 0 итерационным способом (методом секущих).
В формулах (2)-(9) и далее по тексту приняты обозначения: T - температура, t -время, к - коэффициент температуропроводности, y - координата, Q - тепловой эффект реакции, c - удельная теплоемкость, в - глубина превращения пороха, Z - предэкспонент, E - энергия активации, R0 - универсальная газовая постоянная, А - коэффициент теплопроводности, ах - коэффициент теплоотдачи,
q - плотность теплового потока, a - коволюм, A, B - коэффициенты прогонки, At - шаг по времени, a3 - const, Т - безразмерная экспо-
ненциальная функция температуры. Символы: к - конденсированная фаза (твердое топливо), g - газ, 5 - поверхность горения, р - давление, 1 - номер расчетной точки по координате, ] -номер расчетной точки по времени, * - специальное значение.
По полученной скорости горения определяется газоприход от поверхности горения заряда твердого топлива в камеру сгорания РДТТ.
Газовая динамика в камере сгорания РДТТ
Для математического описания процесса течения в камере сгорания, сопловом блоке и за сопловым блоком РДТТ используем подходы механики сплошных многофазных сред, предложенные в работах Х.А. Рахматулина и Р.И. Нигматулина [12, 13]. Согласно этим подходам воздух, газообразные продукты сгорания пиротехнического воспламенительного состава и заряда смесевого твердого топлива назовем первой фазой, мелкодисперсные не-сгоревшие частицы (жидкие капли) в продуктах сгорания пиротехнического воспламени-тельного состава - второй фазой, мелкодисперсные несгоревшие частицы (жидкие капли) малого диаметра в продуктах сгорания заряда смесевого твердого топлива - третьей фазой, мелкодисперсные несгоревшие частицы (жидкие капли) большого диаметра в продуктах сгорания заряда смесевого твердого топлива -четвертой фазой. Первую, вторую, третью и четвертую фазы будем считать гомогенно-гетерогенной смесью со своими температурами и скоростями движения. В такой системе каждая фаза занимает часть объема смеси: а1, а2,а3,а4. Здесь ai =р; /р" (1 = 1...4), где ри - истинная плотность фазы. Движение этих фаз рассматривается как движение взаимопроникающих и взаимодействующих сред.
Тогда полная (нестационарная и трехмерная) система вихревых дифференциальных уравнений газовой динамики для сложного гомогенно-гетерогенного потока в камере сгорания, сопловом блоке и за сопловым блоком ракетного двигателя запишется в виде:
- уравнения неразрывности (сохранения массы)
^ + ¿ы (P1W1 ) = ^ + ^;
(Р2 ) = Ср2в';
Ф2
дГ
др*-+(Рз Wз ) = арз „;
р4и'
(10)
дх
др
др4 дt
(Р4) = Ор4ж;
д(р1ф)
дt ' ф = к, с ц, Х, а;
(рlФWl) = ФвОв + ФО'
^+^ (Р4^4 W4 )+а4 др
= Р4 ^ +тХ4 +иО 0(р4!4 )
дГ
= Р4 + тУ4 +
д(Р4 ^4) дt
= Р4 + Х,4
diу (Р4У4 )+а4~ =
ду
'р 4 и
(Р4^4) + а4 д- =
- уравнения сохранения импульса по осям координат
^+<ь- (ри,^) + а,|=
= РД г-х12-х13-х14 + и О + и О ;
~ 1 ^ х х х в гв и ги'
%> + ¿¿У (р.^ + а, £ =
ОГ ду
= Р^уё-хУ2-хуз-ху4 +УвО8в ;
0(Р1И1) л ( «7\ др
ОГ дг
= р#г-х12-х13-х14 + иО + и О ;
"1 го г г г в гв и ги'
р 4и '
- уравнения сохранения внутренней удельной энергии второй, третьей и четвертой фазы
0(р2—2)
ОГ
°(рз-л )
ОГ
0(р4 —4 )
ОГ
¿¿У (Р2—2^2 ) = Я? + £ + -р2вОр2в;
+ ¿¿У )= Я3 + Я13 + — рЗиОрЗи; (12)
■ ¿¿У (р4 ) = Чк + Я? + —р4иОр4и;
°(рГи^) + ¿¿У (р2и2^^2 ) + а 2 Ох = (11)
= р2 Кг+х1х2 +ивОр 2в;
+ ¿¿У (р2!2 W2 ) + а 2 °р =
дГ v 2 2 ^ 2 Оу
= Р2 Куё + Х!у2 +VвОp 2в;
^ + ¿¿У (Р2 И2 W2 ) + а 2 £
- уравнения сохранения полной удельной энергии смеси
0(р1Е1) , 0(р2Е2) , 0(р3Е3) , 0(р4Е4) +
дГ дГ дГ дГ
+ ШУ^Е^) + ¿¿У(р2Е^2) + ¿¿У(р3Е^3 ) + + ¿¿У (р4Е^4) + ¿¿У (а1 рЩ1) + ¿¿У (а2pW2) +
+
(азPWз)+(а4 ) =Кг Х(р,и,)
(13)
дГ
= Р2 ^ +х12 +ИвОр 2в;
дг
О(ри) + ¿¿У (рзUзWз ) + а3 £ =
= р3 Кг+х13 + и иО р3и;
^^ + ¿¿У (р3v3W3 ) + а3 °р =
дГ v 3 3 ъ! 3 Оу
= р3 Куг+хуз +УиОр 3 и;
+ ¿¿У (р3^3 ) + а3 =
дГ 3 ^ 3 дг
= р3 Кг+х13+иОр 3^;
+^ )+Кг )+Е^
+Е О + Е 2 О 2 + Е 3 О 3 + Е . О . ,
ги ги р 2в р 2в р3и р3и р 4и р4и'
где для декартовой системы координат
дх ду дг
5 = [р, , PlФ, Ри,-, , ри , р/] , рД-, а р ]; (14) / = (...4); ] = (2...4).
Для замыкания системы дифференциальных уравнений (10)-(14) используем уравнение состояния в виде
р = (к - 1)ри
(
Ег -
} 1
1 - ари
(15)
1=1
В уравнениях (10)-(15) и далее по тексту приняты обозначения: ср - удельная теплоемкость при постоянном давлении, Е - полная удельная энергия, J - внутренняя удельная энергия, О - расходно-приходный комплекс, g - ускорение свободного падения, к - показатель адиабаты, N - осевая перегрузка гравитационного поля силы тяжести, q - функция теплового межфазного взаимодействия; и, V, w - проекции вектора скорости по осям координат 0Х, 07, 02; W - вектор скорости, Ж - модуль вектора скорости; х, у, г - текущие координаты вдоль осей 0Х, 07, 02; ц - динамическая вязкость, х - функция силового межфазного взаимодействия. Символы: g - газ, к - конвективный, I - лучистый, р - частицы, ^ - заряд твердого топлива, х - вдоль оси 0Х, у - вдоль оси 07, г - вдоль оси 02.
Выражения (10)-(13) для расходно-при-ходных комплексов, функций силового и теплового межфазного взаимодействия, а также используемые в расчетах дополнительные соотношения подробно изложены в работе [4].
Система уравнений (10)—(15) с учетом дополнительных соотношений интегрируется численно с помощью метода Давыдова (метода крупных частиц) - метода постановки вычислительного эксперимента [3, 4, 8-10, 14-18]. Область интегрирования покрывается фиксированной в трехмерном пространстве (эйлеровой) равномерной (однородной и полностью изотропной) кубической расчетной сеткой с ячейками Ах х Ду х Дг. Для такой сетки геометрический центр ячейки совпадает с центром масс ячейки, что весьма существенно для получения точного результата. Значения целых чисел 1 (вдоль оси 0Х), ] (вдоль оси 07) и к (вдоль оси 02) обозначают геометрический центр ячейки, п - порядковый номер шага по времени.
Эйлеров этап метода. На этом этапе расчета изменяются величины, относящиеся к ячейке в целом, а исследуемая среда предполагается заторможенной (на этом этапе также учитывается вклад гравитационного поля силы тяжести):
PU
,5,j ,k
Л-0,5, j ,k
At
Ax
Pí
+ NxgAt; l = 1...4;
E" = Jn
m, j ,k m, j ,k
(u^ )2+()2
m = 2...4;
аналогично по другим направлениям V k, w" k ;(16)
En = e" — I
j ,k , j ,k
(E" — E" ))
V m.j,k m,j,k )
—I
a, P-n=.-,u, — а, р.^.-.и, At
Ax
P
a p--n=.,v, —a, P-■ At
Ay
P1'
al', j,k+0,5 ^+0,5^^+0,5 —al',j,k—0,5 P'j—0,5Wl',j,k—0,5 At
Az
p1',
+Nxg iRAJ
l=i l=4
+Nyg Z(P" „X J
At
n
pi., j ,k At
pi.,j,k
+Ng I (P" , j)),
l=i l=4
n , „"
n P',j,k + P'+i,j,k —n /1 ir- \
где P.'+0,5,j,k = ^^— ; и\.,,k = (i — alfa)x
хи" k + alfa • ¿I" k и т.д. В формулах (i6) alfa —
сеточный параметр, значение которого изменяется в пределах i,0-3,0.
Лагранжев этап метода. На данном этапе метода вычисляются эффекты переноса (потоковые комплексы), учитывающие обмен между ячейками при их перестройке на прежнюю эйлерову сетку:
- первая фаза вдоль оси 0Х
(P¿Ui )П+0,5,j,k =
( — beta )p" ,j,k + beta -p^ ^ j ^^, если U > 0;
i'+0 5 i k 7
(i — tetaKj + beta-pj ,
если U < 0;
i'+0,5,j ,k
£ = (i, k, Cp, Ц, A, a, ¿i, Vi, Wi, E );
2
m=2
l=i
аналогично по другим направлениям
п+1 п (р3и3 )+0,5,],к (р3и3 )-0,5,],к
(рМ )1]+0,5,к, (рМ)) ,к+0,5; (17) р3-,],к = р3.,],к--Д7-
ДГ -
вторая, третья и четвертая фазы (р ! )п -(р ! )п
_ УК3 3> 1,]+0,5,к УК3 3/ 1,]-0,5,к
„„ V ^ ]+0,5,к У ^ ]-0,5,к . (,о)
вдоль оси 0Х --:-ДГ - (18)
Ду
(р„5и„))
Р т1 ] к 5, ] ,кит,+05,] ,к , (Р3 и3 ). ,.+0_ -(р3 и3 ).
3 V],к+0,5 У^3 3],к-0,5 ДГ + О» ДГ;
" > 0; Дг р3и-,]к
¡+0,5,] ,к
я гп
Щ+1,]к , п+1 п (ри )П+0,5, ],к ~(р4и4У-
рП+1 =рп Уг 4 V¡+0,5,],к Уг4 4>¡-0,5,],к дг -
если ипт < 0; р4- ]к " р4] Дх
5=(1, ит , ^ , *т , - , еш ), т =2...4; — (р4!4 )) ]+0Дк - (рЛ )) ]-0,5,к Д __
Ду
аналогично по другим направлениям _
(р 5! ) ,(р ) и т.д. В форму- -(р4и X],к +0,5 -(ри],к-0,5 ДГ + о" дг■
т,]+0,5,к 'у1™^ ,к+0,5 д^ рк '
лах (17) beta - сеточный параметр, значение п
которого изменяется в пределах -0,2-0,0. п+1 = п рп,,к (рфи1 ^д.л - (р1фи1 Х-)] ДГ
п+1
(Р1Фи1 )п ],к +0,5 - (Рф1 )п ],к-0,5 ДГ
ф, =Ф. -^ ,
На лагранжевом этапе метода вычисля- 1 -] к 1]' р^ Дх р,
ются также приходные комплексы и функции
силового и теплового межфазного взаимодей- , . п , _ . п
ствия, входящие в уравнения (10)-(13), с уче- урф!1/I,]+0,5,к урфУI,._ДГ
том изменения параметров потока на эйлеро- Ду р вом этапе.
Заключительный этап метода. Здесь происходит перераспределение массы, импульса и энергии по пространству и опреде- / \ ДГ
~ + (фп О" +фп О" )——;
ляются окончательные поля эйлеровых пара- у"..к ^ик .к м/ рп+!
метров потока на фиксированной сетке в новый момент времени: ф = (к, ср,ц,Х,а);
- уравнения неразрывности (сохранения
массы) - уравнения сохранения импульса по
Дг Рп+\
рп+1 =рп (р1г^1 X+0,5,]к (р1и1 Х--0,5,]к ДГ -
Р1] Р1] Дх
осям координат
п+1 = ип Р",]к (Р1и1и1 Г+0,5,],к -(Р1и1и1 )0,5,],к
(р,!,) -(р,!1)п иГ1 = и? ^ -4 /1+0,5,],к 4 /¡—0ЛJ,k X
УН1 1 / 1,]+0,5,к Ур 1 Л,]-0,5,к Дг - 11,],к 11,]к ргп+1 Д ^
А "1, ] -к
Ду
(Р1#1 ),к+0,5 -(Р1#1
Дг 4 ~""" " п+г л,, „п+г
.0£ дг + (О; + О; ) ДГ; (Рги1!1 У,]+0,5,к -(Р1и1!1 Г,]-0,5,к ДГ
"¡J'k хрп+1 Ду рп+!
(р1и1и1)п -(р1и1и1)п д Г
/ _ чп , _ V, У1^1 1 ],к +0,5 У1^1 1 гЛ,],к-0,5
п+1 п (р2и;2 А+0,5, ] ,к (Р2и2 Л-0,5, ] ,к Дг Дг рп+!
2 =Р2-----—ДГ- ^ Р^к
Р2 = Р2 ] . Дх
(Р2!2 )п, ]+0,5,к -(Р2 !2 )п; ]-0,5,к дг -(х^" +хх3" +х"" + (19)
---—LJ--———ДГ — \ х1,] к х1.] к х1,] к) рп+г
Ду Р1,],к
(р2 и2 )п, ] к +0,5 (р2 и2 )п, ] к-0,5 дг + Оп ДГ; +(ип Оп + ип Оп );
--1-Дг-1-ДГ + Ор 2в,] к ДГ; +1Ч ] кОв,,]к + ии,] к°и,] к) рп+1 '
„ п+1 ~п 2 ,1 ,к
ио — и~) ,
2¡, 1 ,к 4, 1 ,к рП+1
2,1,
Р^, Лк (Р2и2и2 ) /+0,5,1,к-(р2и2и2 ) П-0,5, ;,к ^^ ) " ,к+0,5 -(Р2J2W2 ) п 0,5 Д
Дх
А (р2и2 Ч )П, 1+0,5,к -(Р2U2V2 Г-сд* Д
„п+1
рп,,
Ду р£
(р2и2^2 )П,1,к+0,5 -(р2и2 ^2 ^-0,5 Д
Дг
п+1
рп„к
-(( к + €,)
Дг Д
„п+1
Р2,
+Jn2в оп
, 1 ,к Д
рп+1 Р2вм,к р2вм,к рп+1
Р2м,к ^2,,1,к
J33+1 = J3
р3п,1,к (р3^3 )п+0,5,1,к -(р3-/3из)
/¿-0,5,1 ,к
3., 1,к 3,, 1,к рп+1 1 к
Дх
Д?
,1 ,к рп+1 в., 1 к р2в,,1 ,к рп+1
< 2, ,к У2, , к
, т-12п Д? . „,п г.
+х —— + и О
Д (р3^3 )п, 1+0,5,к -(р3 J3V3 )п, 1-0,5,к Д
., 1-0,5,к
п+1 _ ~п >,1 ,к
и^ — 7
3,,1 ,к 3,, 1 ,к рп+:1
рп (р3и3и3 )п+0 -(р3и3и3 )п_(
„п+1
р3
1 3,, 1 ,к
Ду
Ря
' 3 3 ^ х (р^) I ;,к+0,5 -(3-^3) I ;,к-0,5 Д?
Дх
Д (Р3^3 )п, 1+0,5,к - (Р3^3 )п, /-0,5,к Д?
', 1-0,5,к
„п+1 р3
1 3, I к
Ду
р3,
(+ < к)
Дг Д?
п+1
(20)
+Jn3w оп
, 1 ,к Д?
рп+1 Р3»',1 ,к Р3»',1 ,к рп+1 '
Р3.,1,к Р3.,1,к
(р3и3^3 X, 1 к+0,5 - (р3и3^3 ) п, 1 к-0,5 Д
¿, 1 к -
Дг
п+1
рд к
тп+1 _ тп ',1 ,к
^ 4,, 1 ,к = ^ 4,, 1 ,к рп+Т
рп,,к (р4 •/4и4) п+0,5,-(р4 •/4и4) ;
/ ¿-0,5,1 ,к
Дх
+ 13п Д + п Оп Д
+х„ —г- + и,., О
х'-1 ,к рп+1 ' "4,1 ,к рп+1 '
Г3, , к Р3, , к
Д? (р4 J4V4 )п, 1+0,5,к - (р4 J4V4 ¡I п, Д
¿, 1-0,5,к
„п+1
рп,,к
Ду
Ря
п+1 _ ~п 1 4.,1 ,к
и4,,1,к = Чм рп+Т
Р1к (р4и4и4 )п+0,5,лк - (р4и4и4 )^0,5,(4 J4^4 )п; к +0,5 - (р4J4)" ; к-0,5 Д?
Дх
х —
Д (Р4U4V4)" ;+0,5,к -(Р4U4V4)" ,0,5,к Д?
„п+1
Р4+1,к
Ду
„п+1
Р4,.„.
"(С +
Д
п+1
Р4,
+Jn4w оп
, 1 ,к Д?
рп+1 Р4»,,1 ,к р4»,,1 ,к рп+1 '
1 ,к
(Р4и4^4)п, +0,5 -(Р4и4^4) п,ик-0,5 Д?
¿, 1,к-0,5
Дг
п+1
- уравнения сохранения полной удельной энергии смеси
14п Д? + _ п оп Д? .
+Хх,, 1,к рп+1 + ,кОр4м>,, 1 ,к рп+1 ;
1 к
>., 1 к
рп т I
Еп+1 _ Еп '' 1 к ^ '
11,1 ,к ^,1 ,к рп+1 /_(
Р1,1,к т=2
(ртЕт )¡,,;■,k (ртЕт ).
1,к
Р1,
аналогично по другим направлениям v3+:к, :к;
- уравнения сохранения внутренней удельной энергии второй, третьей и четвертой фаз
Jn+1 = Jn Р2, 1 к (Р2J2U2 ),+0Д,;к (р2J2U2 ),_0,5„;,к
(р'Е'и')п+0,5,1,к "МЛ)п-0,5,а Д?
0,5,1,к
Дх
Р1,
2,1,к 2,,1к рп+1 1
Дх
-I -"
(рlElVl)) 1+0.-(р'E'V')п 1-0. Д?
Ду
р1,
Д? (Р2 J2V2 )п, 1 +0,5,к - (Р2 J2V2 )Д) -0,5,к Д
(Р'Е'#' )п,М+0,5 -(Р'Е'#' )) „ -0,5 Ы
„п+1
Р2,1к
Ду
Дг
Р2,
п+1
+ (21)
1=1
гЕп О" + Еп О" + Еп 2 О" 2 +Л
+ Е" О" + Е" О"
V р J к Р3 и к Р4 и к Р4 и
Аг
Для повышения точности вычислений (обеспечения условия полной консервативности) в схему метода при расчете давления по уравнению состояния (15) вводится поправка, обеспечивающая (уточняющая) баланс по внутренней удельной энергии смеси [4].
Постановка граничных условий. Для описания граничных условий вокруг расчетной области (камера сгорания, сопловой блок и за сопловым блоком - см. рис. 1) вводятся слои фиктивных ячеек (рис. 2). На закрытых поверхностях (поверхность горения заряда смесевого твердого топлива, стенки камеры сгорания, соплового блока и заглушки соплового блока) используются условия непротекания: нормальная к поверхности компонента вектора скорости из приграничных ячеек сносится в слой фиктивных ячеек с противоположным знаком, остальные параметры потока из приграничных ячеек сносятся в слой фиктивных ячеек без изменения. На открытых поверхностях (расчетная область за сопловыми блоками) используется экстраполяция параметров потока. Приход продуктов сгорания от пиротехнического заряда воспламенительного устройства осуществляется в расчетные ячейки, в которых находится воспламенительное устройство. Приход продуктов сгорания с поверхности горения заряда смесевого твердого топлива осуществляется в расчетные ячейки, центры масс которых расположены в потоке у поверхности горения.
На нерегулярных (несовпадающих с координатной сеткой) криволинейных границах расчетной области применяется предложенный Ю.М. Давыдовым аппарат дробных ячеек (см. рис. 2). Здесь используется процедура нормального отображения фиктивной ячейки относительно границы расчетной области в поток. Везде применяются расчетные формулы только для целых ячеек.
Для вычисления скалярных газодинамических параметров граничной (фиктивной)
ячейки, с учетом расщепления (параметры ф), используются зависимости
Ф* = Х(уь ф,);
IУЬ1 = 1; Ф = Р;,к,...,X,а,,р,,Е], (22)
г
] = 1.4.
Для вычисления векторных газодинамических параметров граничной (фиктивной) ячейки, согласно условиям непротекания, с учетом расщепления (параметры ф ), используются зависимости:
[V, {-Ши]ых + Шщ ]ых 2)];
**=Х (-{ + ^2 N 2)];
[V, (-ЩЫ; + Ж"2]Ы22)];
IV, = 1, ] = 1...4.
(23)
V
1 \ Л
сУ\ ЛГ
/ // \ а
А / 1 Тело
/
Рис. 2. Постановка граничных условий (схема)
В уравнениях (22), (23) приняты обозначения: Шл^Шщ^щ - проекции вектора скорости W на направляющие оси локальной системы координат (Шпх = 0); Ых, Ыу, Ы2 - проекции
нормали к плоскости отображения в базовой системе координат; Ых2, Ыу2, 2 - проекции
оси локальной системы координат в базовой системе координат. Символы: а - фиктивная ячейка, Ь - отраженная фиктивная ячейка.
Разгерметизация камеры сгорания
В силу конструктивных особенностей соплового блока рассматриваемого РДТТ, имеющего заглушку в виде тонкой алюминиевой
X
мембраны сферической формы, завальцованную в корпус сопла, вскрытие с непрогнозируемым разрушением на фрагменты и последующее движение фрагментов заглушки по каналу сопла рассматривать сложно. Ввиду этого разгерметизация камеры сгорания ракетного двигателя реализуется следующим упрощенным образом. До достижения определенного уровня давления продуктов сгорания в канале сопла (так называемого давления страгивания) заглушка сопла есть непроницаемая стенка. При превышении этого давления заглушка сопла мгновенно исчезает.
Комплекс прикладных программ
Для проведения численных расчетов на ЭВМ разработан комплекс прикладных программ, включающий в себя следующие основные модули:
- модуль расчета параметров, описывающих трехмерную постановку граничных условий на криволинейной образующей области интегрирования;
- основной расчетный модуль (main -модуль) - расчет газодинамического течения в камере сгорания, сопловом блоке и за сопловым блоком РДТТ с учетом срабатывания воспламенительного устройства, прогрева, воспламенения и последующего горения заряда смесевого твердого топлива, разгерметизации камеры сгорания;
- модуль визуализации полученной расчетной информации.
Комплекс прикладных программ написан в среде программирования Kdevelop 4.6 для рабочей станции высокой производительности с операционной системой Linux Mint 17.3x86-64 на алгоритмическом языке C/C++ с использованием (для основного расчетного модуля) стандарта многопотоковой обработки информации OpenCL [19-21]. Основная идея стандарта OpenCL состоит в реализации многопотокового выполнения кода по схеме «одна команда - много данных», т.е. одна операция одновременно применяется к большому массиву данных. Подобный подход позволяет в полной мере использовать вычислительный потенциал устройств с множеством относительно простых исполнительных модулей, например
таких, как современные графические процессоры, и существенно (на порядок и более) повысить производительность вычислений.
Проверка работоспособности программного продукта
Проверка работоспособности программного продукта проводилась путем анализа вычислительной устойчивости и точности численного решения, оценки сходимости расчета на различных сетках, а также сравнения результатов расчета с данными стендовых испытаний РДТТ.
Наиболее критичной с точки зрения устойчивости и точности численного решения является задача газодинамического течения. Здесь для анализа свойств разностных схем метода Давыдова (см. выражения (16)-(21)) использовался эвристический подход, основанный на рассмотрении параболической формы их дифференциальных приближений [4, 14]. В этом подходе оценивается знак коэффициентов диффузии ау у диссипативных членов дифференциальных приближений, содержащих частные производные второго порядка по пространственным переменным. Эти коэффициенты обычно группируются в виде матрицы - так называемой матрицы аппроксимационной вязкости. Положительность детерминанта, либо диагональных элементов, либо следа матрицы аппроксимаци-онной вязкости является условием вычислительной устойчивости исследуемой конечно-разностной схемы. Приближая значения этих параметров к нулю, повышаем точность вычислений.
Согласно трактовке лагранжева этапа метода Давыдова [14], каждая проекция конвективных членов исходной системы дифференциальных уравнений (10)-(13) на координатную ось может аппроксимироваться независимо. Это позволяет в рамках определенного приближения при анализе свойств разностной схемы ограничиваться одной координатой. Выпишем для одномерного однофазного аналога используемой выше конечно-разностной схемы метода (16)-(21) диагональные элементы матрицы аппроксимационной вязкости. Будем считать, что поток течет слева направо. Для противоположного направления потока достаточно поменять Дх на -Дх.
критерия устойчивости выбранной конечно-разностной схемы метода.
Проводилась также оценка сходимости численного решения на различных по величине расчетной ячейки разностных сетках. Оценка сходимости показала, что для обеспечения требуемого уровня сходимости численного решения (на момент времени начала процесса срабатывания РДТТ) необходимо иметь в расчетной области ~20 000 000 расчетных ячеек.
Дополнительно проводилось сравнение результатов численного расчета с результатами натурных стендовых испытаний по изменению во времени давления в камере сгорания в районе переднего днища и тяги РДТТ. Результаты расчета и результаты стендовых испытаний хорошо согласуются между собой.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и администрации Пермского края (грант № 19-41-590006-р-а) и отмечена дипломом 1-й степени XXI Всероссийской НТК «Аэрокосмическая техника, высокие технологии и инновации - 2020».
Библиографический список
1. Соркин Р.Е. Газотермодинамика ракетных двигателей на твердом топливе. - М.: Наука, 1967. - 368 с.
2. Алемасов В.Е., Дрегалин А.Ф., Тишин А.П. Теория ракетных двигателей. - М.: Машиностроение, 1980. - 533 с.
3. Численный эксперимент в теории РДТТ / А.М. Липанов, В.П. Бобрышев, А.В. Алиев [и др.]; под ред. А.М. Липанова; УИФ «Наука». - Екатеринбург, 1994. - 301 с.
4. Давыдов Ю.М., Егоров М.Ю. Численное моделирование нестационарных переходных процессов в активных и реактивных двигателях. - М.: НАПН, 1999. - 272 с.
5. Внутренняя баллистика РДТТ / А.В. Алиев, Г.Н. Амарантов [и др.]; под ред. А.М. Липанова и Ю.М. Милехина; РАРАН. - М.: Машиностроение, 2007. - 504 с.
6. Газовые течения в соплах энергоустановок / под ред. проф. В.Н. Емельянова. - М.: ФИЗМАТЛИТ, 2017. - 328 с.
7. Энергетика и внутренняя баллистика ракетных двигателей на твердом топливе / Ю.М. Миле-хин, Г.В. Бурский, Г.С. Лавров, В.С. Попов, Д.Н. Садовничий. - М.: Наука, 2018. - 359 с.
8. Давыдов Ю.М., Давыдова И.М., Егоров М.Ю. Влияние полетной перегрузки на неустойчивость рабочего процесса в камере сгорания ракетного двигателя на твердом топливе // Доклады академии наук. -2004. - Т. 398, № 2. - С. 194-197.
9. Davydov Yu.M., Davydova In.M., Egorov M.Yu. Flight overloading effect on the working process instability in the combustion chamber of a solid-fuel rocket engine // Doklady Physics. Moscow. Nauka / Interpe-riodica. - 2004. - Vol. 49, no. 9. - P. 527-529.
10. Численное моделирование внутрикамерных процессов при выходе на режим работы ракетного двигателя твердого топлива / Г.Н. Амарантов, М.Ю. Егоров, С.М. Егоров, Д.М. Егоров, В.И. Некрасов // Вычислительная механика сплошных сред. - 2010. - Т. 3, № 3. - С. 5-17.
11. Рихтмайер Р.Д., Мортон Х. Разностные методы решения краевых задач. - М.: Мир, 1972. - 420 с.
12. Рахматулин Х.А. Основы газовой динамики взаимопроникающих движений сплошных сред // Прикладная математика и механика. - 1956. - Т. XX, № 2. - С. 184-195.
Итак, имеем:
-^uAx-—u Ax2 +1 (p -u2)At; 2 4 x 2V p '
2 2 x 2 ^ 2
+ 2 (UPJJJX + UPjpP x + PjUx ) 2 +
1 л 1 л 2 1 ( 3 , , A 2
a22 = —puAx —pux Ax —I--beta I upx Ax +
(24)
2
1 (
+—
2
p Л
-3pu2 -ppp--Pj - u2pj At;
p J
гзз = "JpuAx -1 puxAx2 - 2(1 - beta)upxAx2 -
4 pjux Ax 2 - 2 (uPjjjx+uPpj p x )Ax 2+
1
+— 2
-pu + Pj
( P л
E -(1 - 2alfa ) pJ
At,
8п 8 р
гДе пх =—, Р/ = и тД 8х 8р8/
Условие положительности следа матрицы
аппроксимационной вязкости (ап + а22 + а33) > 0
из выражений (24) рассматривалось в качестве
13. Нигматулин Р.И. Динамика многофазных сред. - М.: Наука, 1987. - Ч. I. - 464 с.; ч. II. - 360 с.
14. Давыдов Ю.М. Крупных частиц метод // Математический энциклопедический словарь. - М.: Советская энциклопедия, 1988. - С. 303-304. (Математика. Большой энциклопедический словарь. - 2-е изд. -М.: Российская энциклопедия, 1996. - С. 303-304; Математика. Большой энциклопедический словарь. -3-е изд. - М.: Российская энциклопедия, 1998/2000. - С. 303-304.)
15. Давыдов Ю.М., Давыдова И.М., Егоров М.Ю. Неустойчивость рабочего процесса в двухкамерном ракетном двигателе на твердом топливе // Доклады академии наук. - 2011. - Т. 439, № 2. - С. 188-191.
16. Егоров М.Ю., Егоров Д.М. Численное исследование динамики внутрикамерных процессов при срабатывании бессоплового РДТТ // Известия вузов. Авиационная техника. - 2013. - № 1. - С. 51-54.
17. Егоров М.Ю. Численное исследование динамики внутрикамерных процессов при срабатывании специализированного РДТТ // Известия вузов. Авиационная техника. - 2017. - № 4. - С. 104-111.
18. Егоров М.Ю. Численное исследование динамики переходных внутрикамерных процессов при выходе на режим работы РДТТ особой компоновочной схемы // Известия вузов. Авиационная техника. -2019. - № 2. - С. 59-67.
19. Стахнов А.А. Linux. - СПб.: БХВ-Петербург, 2004. - 912 с.
20. Дерк Л. С и С++: справ. - М.: ВКК, 1997. - 592 с.
21. Programming Guide. AMD Accelerated Parallel Processing OpenCL. Advanced Micro Devices. -2011. - 294 p.
References
1. Sorkin R.E. Gazotermodinamika raketnykh dvigateley na tverdom toplive [Gas-Thermodynamics of solid-fuel rocket engines]. Moscow: Nauka, 1967, 368 p.
2. Alemasov V.E., Dregalin A.F., Tishin A.P. Teoriya raketnykh dvigateley [Theory of rocket engines]. Moscow: Mashinostroyeniye, 1980, 533 p.
3. Lipanov A.M. et al. Chislennyy eksperiment v teorii RDTT [Numerical experiment in the theory of solid propellant motors]. Yekaterinburg: Nauka, 1994, 301 p.
4. Davydov Yu.M., Egorov M.Yu. Chislennoye modelirovaniye nestatsionarnykh perekhodnykh prot-sessov v aktivnykh i reaktivnykh dvigatelyakh [Numerical modeling of non-stationary transient processes in active and jet engines]. Moscow: NAPN, 1999, 272 p.
5. A.V. Aliyev and others. Vnutrennyaya ballistika RDTT [Internal ballistics solid propellant motors]. Ed. by A.M. Lipanov and Y.M. Milekhina. Moscow: Mashinostroyeniye, 2007, 504 p.
6. V.N. Yemelyanova. Gazovyye techeniya v soplakh energoustanovok [Gas flows in the nozzles of power plants]. Moscow: FIZMATLIT, 2017, 328 p.
7. Milekhin Yu.M., Burskiy G.V., Lavrov G.S., Popov V.S., Sadovnichiy D.N. Energetika i vnutrennyaya ballistika raketnykh dvigateley na tverdom toplive [Power Engineering and internal ballistics of rocket engines on solid fuel]. Moscow: Nauka, 2018, 359 p.
8. Davydov Yu.M., Davydova I.M., Egorov M.Yu. Vliyaniye potetnoy peregruzki na neustoychivost rabochego protsessa v kamere sgoraniya raketnogo dvigatelya na tverdom toplive [The influence of flight overload on the instability of the working process in the combustion chamber of a solid-fuel rocket engine]. DOKLADY AKADEMII NAUK, 2004, vol. 398, No. 2, pp. 194-197.
9. Davydov Yu.M., Davydova In.M., Egorov M.Yu. Flight Overloading Effect on the Working Process Instability in the Combustion Chamber of a Solid-Fuel Rocket Engine. Doklady Physics. Moscow. Nauka /Interpe-riodica. 2004, no. 9, Vol. 49, pp. 527-529.
10. Amarantov G.N., Egorov M.Yu., Egorov S.M., Egorov D.M., Nekrasov V.I. Chislennoye modelirovaniye vnutrikamernykh protsessov pri vykhode na rezhim raboty raketnogo dvigatelya ^rdogo topliva [Numerical modeling of intra-chamber processes when entering the operating mode of a solid fuel rocket engine]. Computational continuum mechanics, 2010, vol. 3, no. 3, pp. 5-17.
11. RikhtmayYer R.D., Morton Kh. Raznostnyye metody resheniya krayevykh zadach [Difference methods for solving boundary value problems]. Moscow: Mir, 1972, 420 p.
12. Rakhmatulin Kh.A. Osnovy gazovoy dinamiki vzaimopronikayushchikh dvizheniy sploshnykh sred. PMM [Fundamentals of gas dynamics of interpenetrating motions of continuous media]. Journal of Applied Mathematics and Mechanics, 1956, vol. XX, no. 2, pp. 184-195.
13. Nigmatulin R.I. Dynamics of multiphase media. M.: Nauka, 1987, part I. 464 p.; part II. 360 p.
14. Davydov Yu. M. [Krupnykh chastits metod]. Large particle method. - In the book: Mathematical Encyclopedia. Moscow: Soviet Encyclopedia, 1988, pp. 303-304. Mathematics. A large encyclopedic dictionary.
Ed-e 2-e-M: Russian Encyclopedia, 1996, p. 303-304. Mathematics. A large encyclopedic dictionary. Ed-e 3-e-M: Russian Encyclopedia, 1998/2000, pp. 303-304.
15. Davydov Yu.M., Davydova I.M., Egorov M.Yu. Neustoychivost rabochego protsessa v dvukhkamer-nom raketnom dvigatele na tverdom toplive [Instability of the working process in a two-chamber rocket engine on solid fuel]. Reports of the Russian Academy of Sciences, 2011, vol. 439, no. 2, pp. 188-191.
16. Egorov M.Yu., Egorov D.M. Chislennoye issledovaniye dinamiki vnutrikamernykh protsessov pri srabatyvanii bessoplovogo RDTT [Numerical study of the dynamics of intra-chamber processes in the event of a popleless solid fuel rocket engine]. Russian Aeronautics, 2013, no. 1, pp. 51-54.
17. Egorov M.Yu. Chislennoye issledovaniye dinamiki vnutrikamernykh protsessov pri srabatyvanii spetsializirovannogo RDTT [Numerical study of the dynamics of in-chamber processes when a specialized SRM is triggered]. Aviatsionnaya Tekhnika, 2017, no. 4, pp. 104-111.
18. Egorov M.Yu. Chislennoye issledovaniye dinamiki perekhodnykh vnutrikamernykh protsessov pri vykhode na rezhim raboty RDTT osoboy komponovochnoy skhemy [Numerical study of the dynamics of transient in-chamber processes when entering the operating mode of the special layout circuit RTTT]. Aviatsionnaya Tekhnika, 2019, no. 2, pp. 59-67.
19. Stahnov A.A Linux. St.Peterburg.: BKHV-Peterburg, 2004, 912 p.
20. Derk L. C and C++. Guide. Moscow: VKK, 1997, 592 p.
21. Programming Guide. AMD Accelerated Parallel Processing OpenCL. 2011. Advanced Micro Devices. 294 p.
Об авторах
Егоров Михаил Юрьевич (Пермь, Россия) - доктор физико-математических наук, профессор, профессор кафедры «Высшая математика» ФГБОУ ВО ПНИПУ (614990, г. Пермь, Комсомольский пр., д. 29, e-mail:_egorov-m-j@yandex.ru).
Егоров Дмитрий Михайлович (Пермь, Россия) - кандидат технических наук, первый заместитель генерального директора - главный конструктор АО НИИ полимерных материалов (614113, г. Пермь, ул. Чистопольская, д. 16, e-mail: egorovdimitriy@mail.ru).
Егоров Сергей Михайлович (Пермь, Россия) - кандидат физико-математических наук, начальник расчетного отдела АО НИИ полимерных материалов (614113, г. Пермь, ул. Чистопольская, д. 16, e-mail: know_nothing@bk.ru).
About the authors
Mikhail Yu. Egorov (Perm, Russian Federation) - Doctor of Physical and Mathematical Sciences, Professor of Higher Mathematics Department, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation; e-mail: egorov-m-j@yandex.ru).
Dmitry M. Egorov (Perm, Russian Federation) - Candidate of Technical Sciences, First Deputy General Director - Chief Designer, JSC "Scientific-Research Institute of Polymeric Materials" (16, Chistopolskaya st., Perm, 614113, Russian Federation; e-mail: egorovdimitriy@mail.ru).
Sergey M. Egorov (Perm, Russian Federation) - Candidate of Physical and Mathematical Sciences, Head of the Calculation Department, JSC "Scientific-Research Institute of Polymeric Materials" (16, Chistopolskaya st., Perm, 614113, Russian Federation; e-mail: know_nothing@bk.ru).
Получено 16.02.2021