CC BY
64
DOI: 10.15593/2224-9982/2015.42.07 УДК 621.454.3: 519.63
Р.В. Мормуль1' 2, А.В. Порубов1, М.Ю. Егоров2
1 ПАО «Научно-производственное объединение „Искра"», Пермь, Россия 2 Пермский национальный исследовательский политехнический университет,
Пермь, Россия
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ТЕРМОВИБРАЦИОННОГО
КОНВЕКТИВНОГО ТЕЧЕНИЯ ПРОДУКТОВ СГОРАНИЯ ТВЕРДОГО ТОПЛИВА В УСЛОВИЯХ МИКРОГРАВИТАЦИИ
Проанализирована и исследована структура осредненного и пульсационного движения продуктов сгорания твердого топлива при реализации термовибрационного конвективного режима течения. Изучены механизмы трансформации структуры течения в зависимости от ориентации векторов напряженности гравитационного, вибрационного поля и значения силы Кориолиса. Исследована динамика конвективных течений в широком интервале частот и параметры задачи, определяющие интенсивность конвекции. Проведен анализ степени влияния параметров управления конвективными процессами неизотермических газообразных продуктов сгорания старших ступеней твердотопливных ракетных двигателей на распределение характеристик тепломассо-переноса в переменных по величине и направлению силовых полях в условиях реальной невесомости и микрогравитации. Исследованы параметрические резонансные конвективные течения при высокочастотных модуляциях поля силы тяжести. В рамках численного моделирования решены уравнения тепловой вибрационной конвекции неравномерно нагретого газа для любых взаимных ориентаций силы тяжести и оси вибраций в полости, имеющей форму бесконечного горизонтального цилиндрического слоя для произвольных значений его кривизны при наличии однородного и неоднородного тепловыделения. Показано, что изменение взаимной ориентации указанных векторов в плоскости, перпендикулярной к слою, позволяет контролировать структуру осредненных течений.
Представленные в работе численные результаты имеют фундаментальное значение для понимания общих закономерностей теплопередачи в термогазодинамических системах с объемным энерговыделением, находящихся в модулированном поле внешних сил в условиях орбитального полета. Основная часть численных исследований служит теоретическим базисом при подготовке и проведении экспериментов по гидромеханике невесомости, их лабораторному моделированию в теплообменных устройствах.
Ключевые слова: конвекция, вибрации, твердое топливо, тепловыделение, невесомость, численное моделирование, микрогравитация, тепломассоперенос, метод конечных разностей, теплопроводность.
R.V. Mormul1, 2, A.V. Porubov1, M.Yu. Egorov2
PJSC Research and Production Association "Iskra", Perm, Russian Federation 2 Perm National Research Polytechnic University, Perm, Russian Federation
NUMERICAL SIMULATION OF THERMAL VIBRATION CONVECTION OF COMBUSTION PRODUCTS OF SOLID PROPELLANT UNDER MICROGRAVITY
The article analyzes and studies the averaging and pulsating motion of solid propellant combustion products under thermal and vibration conditions of the convective flow. The transformation mechanisms of the flow structure depending on the orientation of the vectors of the gravitational field and the vibration values of the Coriolis force are studied. The dynamics of convective flows in a wide frequency range and parameters of the problem, determining the intensity of convection are investigated. The influence analysis of control parameters of non-isothermal processes of convective combustion gases of upper stages of solid rocket motors on the distribution of heat and mass transfer characteristics in varying magnitude and direction of the force fields in a real weightlessness and microgravity is conducted. Parametric resonance convection flows in the high-frequency modulation of the gravity field are studied. The equations of thermal and vibration convection of non-uniform heated gas are solved by numerical simulation for any mutual orientation of gravity force and vibration axis in a cavity having the form of endless horizontal cylindrical layer for arbitrary values of its curvature in the presence of homogeneous and non-homogeneous heat. It is shown that the change in the mutual orientation of these vectors in a plane perpendicular to the layer allows controlling the average flow structure.
The presented numerical results are of fundamental importance for the understanding the general laws of heat transfer in thermal gas dynamic systems with volumetric energy emission, which are in the modulated field of external forces under conditions of orbital flight. The main part of numerical studies serves the theoretical basis for preparation and carrying out experiments in zero-gravity hydromechanics, their laboratory simulation in heat exchangers.
Keywords: convection, vibration, solid propellant, heat, weightlessness, numerical simulation, microgravity, heat and mass transfer, finite difference method, the thermal conductivity.
На сегодняшний день большой цикл исследований посвящен изучению конвективных движений жидкостей и газов под воздействием разного рода вибраций. Параметрический характер вибрационного воздействия служит причиной того, что осредненное течение отличается от соответствующего течения без вибрации. Это отличие особенно отчетливо проявляется в предельном случае отсутствия статического поля силы тяжести (невесомость), когда одна лишь вибрация вызывает регулярное осредненное течение. Конвекция, состоящая из осредненной и колебательной компонент, может условно рассматриваться как комбинированное течение, в котором колебательная компонента играет роль вынужденного течения. Переносные тепловые и массовые процессы в неоднородно нагретых жидких и газообразных средах могут играть положительную роль во многих технологических
циклах при перемешивании компонентов топлива космических аппаратов. Серия экспериментальных исследований термовибрационной конвекции на борту орбитальной станции «Мир» представлена в работах [1, 2]. Однако вопрос влияния внутренних источников тепла на процессы конвективного теплообмена продуктов сгорания твердого топлива в условиях микрогравитации является малоизученным.
В первом приближении на примере цилиндрического слоя исследована структура тепловой вибрационной конвекции продуктов сгорания твердого топлива при наличии однородного и неоднородного тепловыделения, закономерности теплопередачи в условиях невесомости. Процесс внутреннего тепловыделения обусловлен протеканием в реагирующей среде химических экзотермических реакций нулевого порядка [3].
Рассмотрим полость, имеющую форму бесконечного цилиндрического слоя, изображенного на рис. 1. Полость окружена массивом, температуропроводность которого много больше температуропроводности х продуктов сгорания твердого топлива. Коэффициент линейного расширения газообразных продуктов течения b и их кинематическую вязкость V будем считать постоянными.
На рис. 1 введены следующие обозначения: (x, y, z) - декартова
система координат; n, k - единичные векторы нормали и вибраций, составляющие с осью x углы а, р соответственно; g - напряженность гравитационного поля; Q0 - циклическая частота вращения цилиндров. Рассмотрим цилиндрическую систему координат (r, 3, z), связанную с декартовой системой соотношениями x = r cos 3, y = r sin 3, z = z.
В настоящей статье моделируются продольные высокочастотные гармонические модуляции напряженности гравитационного поля (инерционные силы). Для рассмотрения поведения газообразных про-
Рис. 1. Геометрия задачи термовибрационной конвекции в цилиндрическом слое
дуктов течения используются уравнения вибрационной конвекции [4], хорошо обоснованные как теоретически [5-8], так и экспериментально [9, 10].
Согласно принципу суперпозиции скорость, температура и давление в системе координат, связанной с полостью, представляются в виде сумм двух слагаемых V + Т + С и р + 5, где V, Т ир - осред-ненные величины, медленно меняющиеся во времени, а С и 5 -пульсационные величины, осциллирующие с частотой ю.
Конвективное движение продуктов сгорания твердого топлива сопровождается тепловыделением, обусловленным протеканием экзотермической реакции согласно закону Аррениуса:
б =
КТ
(1)
где к - константа скорости химической реакции горения твердого топлива; б0 - начальная мощность внутренних источников тепла; Е -
энергия активации химической реакции; К - универсальная газовая постоянная; Т - температура газа.
Систему уравнений Зеньковской-Симоненко при математическом моделировании термовибрационной конвекции продуктов сгорания твердого топлива запишем в безразмерной форме следующим образом:
+ ф -УГ = Дф + ОУТ х п + £ (УЖ •к) хУТ + Та (фх-),
дТ - - 1 Гг ! ^ + V-УТ = —ДТ + — ехр дл Рг Рг I
Д\- + ф = о, ДГ + УТ х к = 0,
V = Уху, Ж = УхГ, У •V = о, У-Ж = о.
1
1 + ВТ
(2)
Е
В системе уравнений (2) введены:
1) безразмерные критерии подобия: О = ^--число Грасго-
V X
фа; 8 = — 2
кЕх ехр
{а.ЫЮ{) 113 Л
V ^ у
КБ 20(
вибрационное число Грасгофа;
V
Бг =---— - параметр Франк-Каменецкого; Рг = -— число
К X
4П2 б 4
Прандтля; Та =-^--число Тейлора;
V 2
2) геометрические параметры модели: а - радиус внутреннего цилиндра; А = а + Б - радиус внешнего цилиндра; Б - величина зазора;
3) соленоидальные вектор-функции: V - поле осредненной скорости конвективного течения; Ж - амплитуда скорости пульсационно-го течения; у, ф - функция тока и вихря скорости осредненного течения; Р - функция тока пульсационного движения газообразных продуктов сгорания;
4) дифференциальные операторы цилиндрической системы коорди-
Л - г д „ Д д2 1 5 1 д2 д2 нат: у= е —- - оператор Гамильтона; Д = —- +--+ ——- +—- -
д^г дг2 г дг г2 дд2 дг2
оператор Лапласа.
Компоненты единичных векторов нормали п и вибраций к в цилиндрической системе координат имеют вид
п = {соб (а-д) , б1п (а-д) ,0}, к = {соб (р-3) , б1п (р-3) ,0}. (3)
Из условия прилипания осредненной скорости течения, условия непротекания пульсационной составляющей скорости газа формируются граничные условия для функций тока:
У = = Р I г=а, а+1 = 0. (4)
дг
Границы полости предполагаются изотермическими:
ТI г=а, а+Б = 0. (5)
Численное решение системы уравнений (2), (3) с учетом граничных условий (4), (5) получено методом конечных разностей [11]. На рис. 2 представлены профили распределения равновесной температуры газа.
Рис. 2. Распределение равновесной температуры в газе (Рг «1) для различных значений числа Франк-Каменецкого: 1 - Гг = 0,4; В = 0,6; 2 - Гг = 0,6; В = 0;
3 - Гг = 0,7; В = 0,2; 4 - Гг = 0,8; В = 0,2; 5 - Гг = 1; В ^то
Максимум равновесной температуры в газе достигается в случае 5, соответствующем однородному внутреннему тепловыделению.
Рассмотрим предельные случаи реализации конвективных режимов течения газа.
1. Естественная (свободная) тепловая конвекция при отсутствии вибраций (О Ф 0,8 = 0).
Вариации значения радиуса внутреннего цилиндра а приводят к различным модификациям структур вихревых течений. Из условия ду дг
асг = 0,758, при котором возникает тот или иной вид вихревого течения газообразных продуктов топлива. При значениях а < асг формируется двухвихревая структура течения, а при значениях а > асг (рис. 3)
течение приобретает двухэтажную структуру по радиусу и в слое образуется четыре вихря. На рис. 3 случаи 1а, 2а, 3а соответствуют распределению ее радиальной части для различных отношений внутреннего радиуса к толщине слоя, а случаи 1Ь, 2Ь, 3Ь - ее полярному распределению относительно радиальной и угловой координат.
г=ас
= 0 оценено критическое значение внутреннего радиуса
Рис. 3. Распределение функции тока при естественном конвективном режиме:
1а - а = асг = 0,758; 2а - а = 1,2; 3а - а = 2,0; О = 1000, 8 = Та = 0, Б = 0,1 м
2. Вибрационная тепловая конвекция в цилиндрическом слое при невесомости.
Численное решение исходной задачи (2)-(5) с условиями О = 0, 8 Ф 0 показано на рис. 4. При увеличении интенсивности вибрационного поля структура течения в слое слабо меняется, при этом сохраняется одноэтажная четырехвихревая структура осредненного течения вне зависимости от отношения внутреннего радиуса цилиндра к величине зазора. При этом осуществляется значительная интенсификация осредненного движения, которая приводит к изменению температурного поля в полости и перераспределению теплового потока на поверхности цилиндров. Конвективное движение газа в условиях микрогравитации обеспечивается достаточно малыми значениями О « 0.
Структура течения газа при термовибрационном конвективном движении продуктов сгорания твердого топлива численно исследована при выполнении условий О Ф 0, 8 Ф 0. Карты распределения функции тока осредненного движения при переходных процессах свободной и вибрационной конвекции представлены на рис. 5. При больших значениях 8 « 500 и достаточно малом значении О « 0,1 (микрогравитация) структура течения принимает вид, как и при вибрационной конвекции, а при малой величине параметра 8 «1 влияние вибрационного поля незначительно и при этом сохраняется структура течения, как и в случае свободной конвекции. Такой механизм трансформации вих-
ревых структур течений объясняется преобладанием скорости осред-ненного движения вибрационного конвективного режима течения над скоростью осредненного течения естественного режима.
Рис. 4. Распределение функции тока осредненного течения вибрационного конвективного режима при невесомости: 1а - а = асг = 0,758; 2а - а = 1,2; 3а - а = 2,0;
О = Та = 0, 8 = 100, Б = 0,1 м
а б
Рис. 5. Распределение функции тока осредненного течения при термовибрационном конвективном режиме: 1 - 8 = 10; 2 - 8 = 117; 3 - 8 = 182; 4 - 8 = 224, Та = 0, а = 0,2 м, Б = 0,1 м; а - в = 100, а = р = 7л/6; б - О = 80, а= 7к/6, р = тс/3
Таким образом, при увеличении параметра 8 происходит относительное уменьшение интенсивности внешних вихрей и увеличение интенсивности внутренних. Данный процесс сопровождается смещением центров вихрей. При совместном воздействии напряженности вибра-
ционного и гравитационного поля осевая симметрия течения газа сохраняется при их коллинеарном или взаимно ортогональном расположении (рис. 5, а), в случае произвольной ориентации - осевая симметрия течения в слое нарушается (рис. 5, б), при этом формирование ярусных структур в слое невозможно.
В рамках численного эксперимента проводилось исследование структуры тепловой вибрационной конвекции при невесомости с учетом влияния силы Кориолиса. Результаты численного моделирования по распределению изолиний функции тока осредненного течения представлены на рис. 6.
Увеличивая вклад влияния силы Кориолиса ^Ког ~Та~ О2 на
протекание вибрационных конвективных процессов в условиях невесомости, т.е. повышая циклическую частоту вращения цилиндров, че-тырехвихревая структура вибрационной конвекции трансформируется в восьмивихревую, состоящую из двух ярусов.
Интенсивность внутренних вихрей, обусловленная действием вибрационных конвективных сил, ослабевает по мере увеличения числа Тейлора, вихри разных знаков меняются местами (поворачиваются), при этом осевая симметрия течения в слое не нарушается. Постепенно формируется пограничный слой, в который «сжимается» течение, наблюдаемое при малых частотах вибраций, а его место занимает течение похожей структуры, но другого знака.
Рис. 6. Распределение функции тока осредненного течения при вибрационном конвективном режиме в условиях невесомости под воздействием силы Кориолиса: 1 -
Та = 1; 2 - Та = 4; 3 - Та = 11; 4 - Та = 18; 5 - Та = 23; 6 - Та = 28; 7 - Та = 33; 8 - Та = 35; 8 = 70, а = 0,2, Б = 0,1, О = 0
Для описания закономерностей теплопередачи на поверхности цилиндров введем в рассмотрение локальные тепловые потоки: ^(9) -
для поверхности внутреннего цилиндра и ^(9) - для поверхности
внешнего. Будем определять их на некотором малом расстоянии 5 от стенок слоя, на котором конвективная часть теплового потока пренебрежимо мала по сравнению с теплопроводной:
^ ^^ , ^(9) = .
При значении тепловых потоков ^ (9) > 0 или ^ (9) < 0 через
окрестность точки на поверхности внутреннего (внешнего) цилиндра с координатой & происходит передача тепла от газа к соответствующему цилиндру, а при значении ^ (9) < 0, ^(9) > 0 - от цилиндра к газу.
Общие тепловые потоки через поверхность цилиндров определяются как
2л 2л
N1 = (а + 5) | ^ (9)49, N2 = (А -5) | ^ (9)49.
0 0
На рис. 7 представлена зависимость отношения общих тепловых
1 А - а
потоков внутреннего и внешнего цилиндров от кривизны 1 =-
а + А
цилиндрического слоя в случае однородного тепловыделения (см. рис. 2, случай 5).
Рис. 7. Зависимость отношения общих тепловых потоков внутреннего и внешнего цилиндров от кривизны слоя
Пульсационные характеристики среды определяются по осред-ненным формулам:
^ = -О-1>/28 (Ж Т) эт О1, Л = у[288ЖсовО,
юБ 2
где О =-= 81 - безразмерная частота колебаний (параметр Стокса).
V
Карты распределения функции тока пульсационного течения представлены на рис. 8, где:
1 - свободный конвективный режим: О = 1000; 8 = Та = 0; Г = 0,45; Бг = 0,6; а = 90°;
2 - вибрационный конвективный режим: О = Та = 0; 8 = 500; Г = 0,45; Бг = 0,6; р = 90°;
3 - термовибрационный конвективный режим: О = 100; 8 = 1000; Г = 0,4; Бг = 0,6; а = 90°; р = 180°;
4 - термовибрационный конвективный режим: О = 100; 8 = 1500; Г = 0,6; Бг = 0,6; а = 90°; р = 180°.
Рис. 8. Карты функции тока пульсационного течения
Достоверность результатов распределения параметров тепломас-сопереноса при реализации свободной и вибрационной тепловой конвекции в частном случае однородного тепловыделения подтверждается
данными численного расчета, представленного в теоретических работах [12, 13], и данными орбитальных опытов [2], позволяющих расширить представление о термовибрационных эффектах в условиях реальной невесомости.
Использование методики численного интегрирования системы дифференциальных уравнений (2) с помощью конечно-разностных аппроксимаций позволяет воспроизвести термовибрационные конвективные режимы течения, подобные обнаруженным в ходе космических экспериментов.
Библиографический список
1. Полежаев В.И., Соболева Е.Б. Тепловая гравитационная и вибрационная конвекция околокритического газа в условиях микрогравитации // Известия РАН. Механика жидкости и газа. - 2000. - № 3. -С. 70-80.
2. Исследование околокритической жидкости в условиях микрогравитации: эксперименты на станции «Мир» и численное моделирование / А.В. Зюзгин, А.И. Иванов, В.И. Полежаев, Г.Ф. Путин, Е.Б. Соболева // Космонавтика и ракетостроение. - 2000. - № 19. - С. 56-63.
3. Франк-Каменецкий Д. А. Диффузия и теплопередача в химической кинетике. - М.: Наука, 1987. - 502 с.
4. Зеньковская С.М. Исследование конвекции в слое жидкости при наличии вибрационных сил // Известия АН СССР. Механика жидкости и газа. - 1968. - № 1. - С. 55-58.
5. Повицкий А. С., Любин Л.Я. Основы динамики и тепломассообмена жидкостей и газов в невесомости. - М.: Машиностроение, 1972. - 252 с.
6. Бабский В.Г., Копачевский Н.Д., Мышкис Л.Д. Гидродинамика невесомости. - М.: Наука, 1976. - 504 с.
7. Бурдэ Г.И. Численное исследование конвекции, возникающей в модулированном поле внешних сил // Известия АН СССР. Механика жидкости и газа. - 1970. - № 2. - С. 196-201.
8. Гершуни Г.З., Любимов Д.В. Термовибрационная конвекция, 1998. - 358 с.
9. Иванова А.А., Козлов В.Г. Экспериментальное исследование тепломассопереноса в условиях вибрационной конвекции // 3-й Всесо-юз. семинар по гидромеханике и тепломассопереносу в невесомости. -Черноголовка, 1984.
10. Острах С. Роль конвекции в технологических процессах, проводимых в условиях микрогравитации // Космическая технология / под ред. Л. Стега. - М.: Мир, 1980. - С. 9-27.
11. Ладыженская О. А. Метод конечных разностей в теории уравнений c частными производными // Успехи математических наук. -1957. - Т. 12, № 5(77) - С. 123-148.
12. Roschina N.A., Uvarov A.V., Osipov A.I. Natural convection in an annulus between coaxial horizontal cylinders with internal heat generation // International Journal of Heat Mass Transfer. - 2005. - Vol. 48. -Р. 4518-4525.
13. Шарифулин А.Н. Численное исследование конечно-амплитудной вибрационной конвекции в цилиндрическом слое // Конвективные течения / Перм. пед. ин-т. - Пермь, 1983. - С. 102-107.
References
1. Polezhaev V.I., Soboleva E.B. Teplovaya gravitatsionnaya i vibratsionnaya konvektsiya okolokriticheskogo gaza v usloviyakh mikro-gravitatsii [Thermal gravitational convection and vibration of near-critical gas in microgravity]. Izvestiya Rossiyskoy akademii nauk. Mekhanika zhidkosti i gaza, 2000, no. 3, pp. 70-80.
2. Zyuzgin A.V., Ivanov A.I., Polezhaev V.I., Putin G.F., Soboleva E.B. Issledovanie okolokriticheskoy zhidkosti v usloviyakh mikrogra-vitatsii: eksperimenty na stantsii "Mir" i chislennoe modelirovanie [Investigation of near-critical fluids in microgravity: experiments at the station "Mir" and numerical simulation]. Kosmonavtika i raketostroenie, 2000, no. 19, pp. 56-63.
3. Frank-Kamenetskiy D.A. Diffuziya i teploperedacha v khimiches-koy kinetike [Diffusion and heat transfer in chemical kinetics]. Moscow: Nauka, 1987. 502 p.
4. Zenkovskaya S.M. Issledovanie konvektsii v sloe zhidkosti pri nalichii vibratsionnykh sil [Study of convection in a liquid layer with vibrating forces]. Izvestiya Akademii nauk SSSR. Mekhanika zhidkosti i gaza, 1968, no. 1, pp. 55-58.
5. Povitskiy A.S., Lyubin L.Ya. Osnovy dinamiki i teplomassoobmena zhidkostey i gazov v nevesomosti [Fundamentals of dynamics and heat and mass transfer fluids in weightlessness]. Moscow: Mashinostroenie, 1972. 252 p.
6. Babskiy V.G., Kopachevskiy N.D., Myshkis L.D. Gidrodinamika nevesomosti [Hydrodynamics of weightlessness]. Moscow: Nauka, 1976. 504 p.
7. Burde G.I. Chislennoe issledovanie konvektsii, voznikayushchey v modulirovannom pole vneshnikh sil [Numerical study of convection occurring in the modulated field of external forces]. Izvestiya Akademii nauk SSSR Mekhanika zhidkosti i gaza, 1970, no. 2, pp. 196-201.
8. Gershuni G.Z., Lyubimov D.V. Termovibratsionnaya konvektsiya [Thermal and vibration convection]. 1998. 358 p.
9. Ivanova A.A., Kozlov V.G. Eksperimentalnoe issledovanie teplomassoperenosa v usloviyakh vibratsionnoy konvektsii [Experimental study of heat and mass transfer in a vibration convection]. III Vsesoyuznyy seminar po gidromekhanike i teplomassoperenosu v nevesomosti. Cherno-golovka, 1984.
10. Ostrakh S. Rol konvektsii v tekhnologicheskikh protsessakh, provodimykh v usloviyakh mikrogravitatsii. [The role of convection processes conducted in microgravity]. Kosmicheskaya tekhnologiya. Ed. L. Steg. Moscow: Mir, 1980, pp. 9-27.
11. Ladyzhenskaya O.A. Metod konechnykh raznostey v teorii uravneniy c chastnymi proizvodnymi [Finite difference method in the theory of partial derivative equations]. Uspekhi matematicheskikh nauk, 1957, vol. 12, no. 5(77), pp. 123-148.
12. Roschina N.A., Uvarov A.V., Osipov A.I. Natural convection in an annulus between coaxial horizontal cylinders with internal heat generation. International Journal of Heat Mass Transfer, 2005, vol. 48, pp. 4518-4525.
13. Sharifulin A.N. Chislennoe issledovanie konechno-amplitudnoy vibratsionnoy konvektsii v tsilindricheskom sloe [Numerical investigation of the vibration amplitude convection in a cylindrical layer]. Konvektivnye techeniya. Permskiy pedagogicheskiy institut, 1983, pp. 102-107.
Об авторах
Мормуль Роман Викторович (Пермь, Россия) - инженер-конструктор 2-й категории, ПАО «Научно-производственное объединение „Искра"» (614038, г. Пермь, ул. Академика Веденеева, д. 28), аспирант кафедры «Высшая математика» ФГБОУ ВПО ПНИПУ (614990, г. Пермь, Комсомольский пр., д. 29, e-mail: rmormul@yandex.ru).
Порубов Андрей Венальевич (Пермь, Росия) - инженер-конструктор 3-й категории, ПАО «Научно-производственное объединение „Искра"» (614038, г. Пермь, ул. Академика Веденеева, д. 28, e-mail: leo2701@bk.ru).
Егоров Михаил Юрьевич (Пермь, Россия) - доктор физико-математических наук, профессор кафедры «Высшая математика» ФГБОУ ВПО ПНИПУ (614990, г. Пермь, Комсомольский пр., д. 29, e-mail: egorov-m-j@yandex.ru).
About the authors
Roman V. Mormul (Perm, Russian Federation) - Design engineer, PJSC Research and Production Association "Iskra" (28, Academica Vedeneeva st., Perm, 614038, Russian Federation), Doctoral Student, Department of Higher Mathematics, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: rmormul@yandex.ru).
Andrey V. Porubov (Perm, Russian Federation) - Design engineer, PJSC Research and Production Association "Iskra" (28, Academica Vedeneeva st., Perm, 614038, Russian Federation, e-mail: leo2701@bk.ru).
Mikhail Yu. Egorov (Perm, Russian Federation) - Doctor of Physical and Mathematical Sciences, Professor, Department of Higher Mathematics, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: egorov-m-j@yandex.ru).
Получено 15.05.2015
