CC BY
24
УДК 532.529:518.5
БО!: 10.14529/ттр170404
ОСЛАБЛЕНИЕ СФЕРИЧЕСКИХ УДАРНЫХ ВОЛН В ГЕТЕРОГЕННЫХ СРЕДАХ
Н.Л. Клиначева, Ю.М. Ковалев
Южно-Уральский государственный университет, г. Челябинск, Российская Федерация
В данной работе приводится анализ численного исследования взаимодействия сферических ударных волн (УВ) с механическими и тепловыми неоднородностями, расположенными на плоской поверхности. Результаты расчетов показали, что наличие на плоской поверхности гетерогенного слоя приводит к ослаблению падающей УВ внутри этого слоя в зависимости от объемной доли конденсированной фазы. Достаточно большая объемная доля твердой фазы в гетерогенном слое приводит к деформации головного фронта УВ и сильное ее торможение, в результате чего отражения от стенки практически не происходит и ножка Маха не образуется, однако при а2 > 5 • 10-2 - она образуется над поверхностью гетерогенного слоя. Увеличение температуры газа в слое приводит к возрастанию скорости (УВ), в то время как амплитуда ее снижается. Приведенные в работе расчеты показали, что увеличение температуры газа в слое способствует увеличению момента силы, действующей на решетку из твердых частиц. При высоких температурах газа в гетерогенном слое, несмотря на силовое и тепловое взаимодействие газа и частиц, происходит образование предвестника.
Ключевые слова: математическое моделирование; уравнения Эйлера; газовзвесь; ударные волны; число Куранта.
Введение
Исследование поведения ударных волн (УВ) в неоднородных средах представляет как практический, так и теоретический интерес. Это связано с тем, что в настоящее время являются актуальными проблемы обеспечения безопасности при изучении быстропротекающих процессов, сопровождающихся образованием УВ в результате взрыва в атмосфере различных тел (например, метеоритов). В последние годы активность вокруг проблемы «Астероидно-кометная опасность> резко возросла. Вопросы возможных столкновений Земли с астероидами и кометами и их вероятных последствий начали обсуждаться на уровне правительств, национальных и международных организаций и ведомств. При входе космического тела в плотные слои атмосферы происходит резкое торможение тел, в ударном слое создается давление порядка тысячи атмосфер, возникают большие лучистые и тепловые потоки. Внутри тела поя в-ляются значительные массовые силы инерции; под действием высоких механических и термических напряжений происходит полное или частичное разрушение тела с выделением энергии. Все это приводит к образованию сферической УВ, которая взаимодействует с механическими и тепловыми неоднородностями вблизи поверхности Земли.
К настоящему времени проведено достаточно большое количество исследований процессов взаимодействия плоской У В с гетерогенными слоями, посвященных различным аспектам данной проблемы: ослабление УВ при прохождении их через решетки или слои запыленного газа [1-7]; инвариантность к преобразованиям Галилея законов сохранения, описывающих движение многокомпонентных сред [8, 9], и т.д.
В работах [10. 11] было установлено, что при взаимодействии мелкомасштабных сферических ударных волн? инициированных лазерным пробоем, с абсолютно твердой поверхностью в присутствии теплового слоя происходит: частичное запирание УВ внутри слоя, образование предвестника, формирование двух тройных точек. Если изучению закономерностей взаимодействия плоских УВ посвящено достаточно большое количество исследований, то процессы взаимодействия сферических УВ с гетерогенным слоем на плоской поверхности в настоящее время изучены еще недостаточно. несмотря на значительный практический интерес.
Целью ДаННОИ работы является исследование особенностей взаимодействия сферической УВ с гетерогенным слоем, расположенным на плоской поверхности.
1. Постановка задачи
Рассмотрим задачу о взаимодействии сферической УВ с гетерогенным слоем твердых частиц в следующей постановке! пусть в результате взрыва над поверхностью Земли образовался сферический объем сильно нагретого газа. Давление внутри этого объема ПОСТОЯННО Р = со^ (рис. 1). Центр объема находится на ос и симметрии Ог цилиндрической системы координат. В качестве расчетной области выбирается плоскость, левая граница которой совпадает с осью Ог, а нижняя - с осью Ог. К поверхности прилегает гетерогенный слой высотой И. Для простоты проведения расчетов и анализа результатов газ будем считать идеальным во всей расчетной области. Реализация расчетов с более сложными уравнениями состояния [12] в области повышенного давления не представляет принципиальных трудностей. Для изучения движения УВ в гетерогенном слое принимаются обычные для «замороженной> газовзвеси допущения. Согласно этим допущениям гетерогенный слой можно рассматривать как неде-формированную решетку, заполненную газом. Узлами решетки являются твердые частицы, а влиянием связей пренебрегаем. Обоснование модели кзамороженнойтъ газовзвеси для описания законов взаимодействия с преградами УВ в слоях запыленного газа и решеток детально обсуждалась в работах [2. 3, 6].
к
\ 1
\R
1У
П
H Ш
Рис. 1. Схема расчетной области в начальный момент времени
При достаточно сильном взрыве на начальном этапе развития процесса в гетерогенном слое значительную роль играет излучение, но в данной постановке задачи его влияние не рассматривается. Это обусловлено тем. что в рассматриваемой задаче
«высвечивание> области энерговыделения уже произошло, и на момент входа УВ в гетерогенный слой роль излучения несущественна.
В качестве математических моделей изучаемого физического процесса выбраны система уравнений Эйлера, записанная в цилиндрической системе координат в области чистого газа, модель «замороженной> газовзвеси, которые замыкаются уравнениями состояния фаз. С целью получения обобщенных результатов система законов сохранения была приведена к безразмерному виду. Характерными параметрами для
обезразмеривания были выбраны: высота гетерогенного слоя Н, плотность при нор-
0^02 2 мальных условиях р0, скоростной напор р0 • а2 и квадрат скорости звука а0 в невозмущенном газе соответственно. безразмерная скорость газа определялась отношением скорости газа к скорости звука а0, а характерное время значением Ь0 = Н/а0. Вид остальных безразмерных параметров определен соотношениями (1):
г = Нг, г = Нг и = а0и, V = а0у
Р = Р°Р, Р = Р0а0Р Ь = Ь0Ь, Е
(1)
а0Е.
Для выбранных параметров система законов сохранения и уравнения состояния в безразмерной форме имеют следующий вид:
д (ргг) + д (риг) + д (р1 Vlr) = 0 дЬ дг дг
д (риг) д (ргщщг) д (ргииг) д (гагР)
-тт.--1--~--1--~-= а\Р-----гЕг аг,
дЬ дг дг дг
д (ргVlr) д (рlVlUlг) д (ргVlVlr) д (гагР) „
-^--1--~--1--~-=--~--гГх аг,
дЬ дг дг дг
д (ргЕгг) д (ргщЕгг) д (риЕгг) д (щагРг) д Рг)
-^--1--~--1--~--1--~--1--~-= -гЯаг,
дЬ дг дг дг дг
др2&2 ~ 0 Р2
= Я, W2 = 0, р° = — = СО^
дЬ а2
Р Е + Ж о рг .
ег = ~0Г<-ГТ, Ег = ег + -—, рг = —, аг = 1 - а2.
(2)
р0 (7 - 1У г г' 2 ' - аг
При записи системы законов сохранения (2) черта над безразмерными переменными для удобства была опущена. Здесь г и г - радиальная и осевая координаты; Ь -время, и и V - компоненты вектора скорости Ж; Р - давление; ег, Ег, р0 удельная внутренняя и полная энергии газа, истинная плотность; 7 - показатель адиабаты газа; а2 - объемная доля частиц; ср - удельная теплоемкость при постоянном давлении твердой фазы; Г и Я - интенсивности силового и теплового взаимодействия газа с
узлами решетки:
^ = 0, ЧЬСааР^тБка, = 0, 756^?!
24 4 —
С? = С? (Ее) Ф (М) ф (а), С? = — + -= + 0, 4, Ее = р°| - ,
Ее V Ее
К \ ! 1АТ { р о о — Н ^
^ = 6а2Л + ( в1 - —еП ргЕво , Яв? = р1 а°—, ¿Д? = — ,
Мп = 2ежр (-М) + 0, 459Яе°'55Рг°'33, Рг = ,
А1
где С? - коэффициент аэродинамического сопротивления сферы; С? - коэффициент аэродинамического сопротивления частиц в несжимаемом безграничном потоке газа; М, Яе, N11 и 1'г числа Маха, Рейнольдса, Нуссельта и Прандтля; (I диаметр частиц; А1; а1 - теплопроводность, вязкость и местная скорость звука газа; а° -характерная скорость звука в невозмущенном газе.
Система дифференциальных уравнений (2) - (4) решалась методом крупных частиц [13, 14] и его модификациями [15, 16]. В начальный момент времени параметры имеют следующие значения: в зонах II и III газ находится при нормальных условиях Р° = 105 Па и плотность р° = 1, 21 кг/м3. В 1 зоне давление равно 2500Р°, высота взрыва равна 7 Н, а давление за У В у поверхности гетерогенного слоя составляло около 5 атм. Теплоемкость газа с1 =716 Дж/(кг • К), теплопроводность А1 = 26,2 • 105 Дж/(м- К), вязкоеть /л1 = 18, 5 • 105(Па•с), местная старость звука а1=341 м/с, характерная скорость звука ао=287 м/с. Для твердой фазы диаметр частиц принимался
р2 = 7800 / 3 2 •
емная доля частиц изменялась от а2 = 10_4 до а2 = 5 • 10_2. Расчеты проводились на сетке размером 112x72 узла.
2. Основные результаты
При значениях а2 < 1 • 10_3 сферическая УВ, падающая на гетерогенный слой частично отражается от его поверхности, в результате образуется ударная волна, которая движется к центру взрыва (рис. 2). При значениях а2 < 5 • 10_4 отражения УВ от верхней границы гетерогенного слоя практически нет, поскольку наличие неодно-родностей не оказывает существенного влияния на картину течения. В этом случае УВ проходя внутрь слоя достигает стенки и отражается от нее. Происходит переход от регулярного отражения к нерегулярному и образуется ножка Маха. Фронт ударной волны остается близким к сферическому.
Увеличение объемной доли частиц в слое (а2 > 5 • 10_3) приводит к изменению формы головного фронта УВ. Часть УВ отражается от его поверхности, при этом отраженная волна движется в направлении к центру взрыва (рис. 3), а часть проходит в слой, где испытывает торможение. Интенсивность УВ в слое становится слабой и она практически гасится внутри него.
Дальнейшее увеличение объемной доли частиц до значений а2 > 5 • 10_2 приводит к сильной деформации головного фронта УВ внутри решетки. Ножка Маха вблизи стенки не образуется, но она формируется на верхней границе гетерогенного слоя (рис. 4).
Рис. 4. Изолинии давлений а2 =5 • 10 2
Зависимость давления от времени в эпицентре взрыва представлена на рис. 5: кривая 1 - а2 = 0 кривая 2 - а2 = 0, 0001, кривая 3 а2 = 0, 001 и кривая 4 - а2 = 0, 005. Видно, что с ростом а2 растет и время, за которое УВ достигает поверхности, а ее амплитуда падает. С увеличением объемной доли частиц в слое максимальное давление в эпицентре взрыва падает (рис. 6), а время, за которое достигается максимум давления, линейно растет с ростом а2 (рис. 7).
Распределения давления вдоль стенки представлены на рис. 8. Видно, что наличие частиц приводит к снижению амплитуды падающей сферической УВ и более медленному падению давления за фронтом УВ. Для больших значений объемной доли УВ в слое не имеет четко выраженного фронта.
Рассмотрим теперь влияние температуры газа в слое на динамику течения. В данном случае характерными параметрами задачи являются объемная доля твердой фазы и относительная плотность газа в теплом слое w. Первый параметр оказывает влияние на характер отражения УВ от верхней границы гетерогенного слоя, а второй на образование предвестника.
На рис. 9 представлены поля давлений для трех значений параметра w на один и тот же момент времени. Скорость УВ, движущейся по слою, тем больше, чем вьтттте температура газа в слое. Уменьшение температуры газа в гетерогенном слое приводит к деформации фронта ударной волны, движущейся по этому слою. Ножка Маха растет с увеличением прогрева газа в гетерогенном слое. При значении параметра w 0.1.
р
Рис. 5. Зависимость давления в эпицентре взрыва от времени
t
О 85 —
О 0 002 0004 0 006
Рис. 7. Зависимость времени достижения максимальногого давления в эпицентре взрыва от объемной доли частиц
О О 002 0 004 0 006
Рис. 6. Зависимость времени достижения максимальногого давления в эпицентре взрыва от объемной доли частиц
0 12 3 4
Рис. 8. Распределение давления вдоль по— верхности
несмотря на наличие силового и теплового взаимодействия газа и частиц. УВ в слое ускоряется настолько, что обгоняет головной фронт, что. в свою очередь, приводит к образованию предвестника.
На рис. 10. 11 представлены зависимости давлений в эпицентре взрыва от времени и от координаты на момент времени Ь 0.76: кривая 1 ш 0.1. кривая 2 ш 0.5, кривая 3 ш 1. Рост температуры приводит к снижению интенсивности УВ. а ее скорость увеличивается.
Анализ силового воздействия сферической УВ на наземные ооъекты мо^кно провести. рассчитав значение момента силы, действующей на решетку:
н
М (г,г) = У (4)
о
w = 0,1
w = 0, 5
w = 1
Рис. 9. Изолинии давлений
Рис. 10. Зависимость давления от време- Рис. 11. Распеделение давления вдоль по~ ни в эпицентре взрыва верхности
В данной работе представлено распределение максимального по времени момента силы, действующей на решетку вдоль стенки (рис. 12):
Mmax (r) = max M (r, t) 0<t<t'
(5)
где - время окончания расчета по координате т.
Увеличение температуры газа в слое при фиксированной объемной доле способствует увеличению максимального момента силы.
Выводы
В результате проведения численных исследований было показано:
1. Наличие гетерогенного слоя на поверхности приводит к ослаблению падающей УВ внутри этого слоя.
2. При достаточно большой объемной доле твердой фазы происходит деформация головного фронта УВ и сильное ее торможение, в результате чего отражения от стенки практически не происходит и ножка Маха не образуется, однако при а2 > 5 • 10_2
она образуется над поверхностью гетерогенного слоя.
3. С ростом температуры газа в гетерогенном слое скорость ударной волны воз-растает^ а интенсивность падает.
4. Момент силы, действующей на решетку, тем больше, чем вьтттте температура газа в гетерогенном слое.
5. При высоких температурах газа в гетерогенном слое, несмотря на силовое и тепловое взаимодействие газа и частиц, происходит образование предвестника.
Литература
1. Иваттдаев, А.И. Влияние экранирующих слоев газовзвеси тта отражение ударных воли / А.И. Иваттдаев, А.Г. Кутуптев /'/ Прикладная механика и техническая физика. - 1985. -№ 21. - С. 115-120.
2. Кругликов, Б.С. Ослабление воздушных ударных волтт слоями занылеттттого газа и решетками / Б.С. Кругликов, А.Г. Кутуптев /'/ Прикладная механика и техническая физика. - 1988. - № 1. - С. 51-57.
3. Кругликов, Б.С. Ослабление ударных волтт экранирующими решетками / Б.С. Кругликов, А.Г. Кутуптев /'/ Физика горения и взрыва. - 1988. - № 1. - С. 115-117.
4. Газовая динамика многофазных сред. Ударные и детоттациоттттые волтты в газовзвесях / А.И. Иваттдаев, А.Г. Кутуптев, Р.И. Нигматулитт /'/ Итоги ттауки и техники. Серия: Механика жидкости и газа; Всероссийский институт научной и технической информации. - 1981. - Т. 16. - С. 209-287.
5. Кутушев, А.Г. Ослабление воздушных ударных волн слоями запыленного газа и решетками / А.Г. Кутушев, Д.А. Рудаков // Прикладная механика и техническая физика. -1993. - № 25. - С. 25-31.
6. Ковалев, Ю.М. Ослабление воздушных ударных волн системой решеток / Ю.М. Ковалев, А.Ю. Черемохов // Вопросы атомной науки и техники. Серия: Математическое моделирование физических процессов. - 1997. - № 3. - С. 39-43.
7. Беляев, П.Е. Влияние экранирующего слоя газовзвеси на силовое воздействие ударной волны на жесткую стенку / П.Е. Беляев, Н.Л. Клиначева // Вестник ЮУрГУ. Серия: Математика. Механика. Физика. - 2016. - Т. 8, № 4.- С. 49-55.
8. Клиначева, Н.Л. Модифицированная математическая модель «замороженной» газовзвеси / Н.Л. Клиначева, Ю.М. Ковалев, В.Ф. Куропатенко // Инженерно-физический журнал. - 2014. - Т. 87, № 6.- С. 1398-1403.
9. Ковалев, Ю.М. Математический анализ уравнений сохранения двухфазных смесей / Ю.М. Ковалев, Е.А. Ковалева // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2014. - Т. 7, № 2.- С. 29-37.
10. Андрущенко, В.А. Отражение сферической ударной волны от плоскости при наличии на ней слоя нагретого газа / В.А. Андрущенко, М.В. Мещеряков, Л.А. Чудов // Механика жидкости и газа. - 1989. - № 24. - С. 141-147.
11. Андрущенко, В.А. Взаимодействие сферических ударных волн с приповерхностными тепловыми газовыми неоднородностями / В.А. Андрущенко, М.В. Мещеряков // Физика горения и взрыва. - 1990. - № 3. - С. 77-82.
12. Антонов, В.А. Моделирование взрыва шнурового заряда в пологе леса при отсутствии пожара / В.А. Антонов, А.М. Гришин, Ю.М. Ковалев, Л.Ю. Наймушина // Физика горения и взрыва. - 1993. - № 4. - С. 115-123.
13. Белоцерковский, О.М. Метод крупных частиц в газовой динамике / О.М. Белоцерков-ский, Ю.М. Давыдов. - М.: Наука, 1982.
14. Белоцерковский, О.М. Нестационарный «метод крупных частиц» для газодинамических расчетов / О.М. Белоцерковский, Ю.М. Давыдов // Прикладная механика и техническая физика. - 1971. - № 21. - С. 182-207.
15. Гришин, Ю.А. Повышение устойчивости вычислительного алгоритма метода крупных частиц / Ю.А. Гришин, В.А. Зенкин // Наука и образование. - 2011. - № 9. - 5 с. - URL: http://old.technomag.edu.ru/doc/221488.html
16. Ковалев, Ю.М. Анализ некоторых модификаций метода крупных частиц на примере исследования течений газовзвесей / Ю.М. Ковалев, Е.А. Ковалева, Е.Е. Пигасов // Вестник ЮУрГУ. Серия: Математика. Механика. Физика. - 2015. - Т. 7, № 3.- С. 71-77.
Наталия Леонидовна Клиначева, кандидат физико-математических наук, доцент, кафедра «Вычислительная механика:», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Юрий Михайлович Ковалев, доктор физико-математических наук, профессор, кафедра «Вычислительная механика», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Поступила в редакцию 25 августа 2017 г.
MSC 76L05 DOI: 10.14529/mmp170404
THE ATTENUATION OF SPHERICAL SHOCK WAVES IN HETEROGENEOUS MEDIA
N.L. Klinacheva, Y.M. Kovalev,
South Ural State University, Chelyabinsk, Russian Federation E-mail: [email protected], [email protected]
This article represents the analysis of the numerical investigation of the interaction of spherical shock waves (SW) with mechanical and thermal inhomogeneities located on a flat surface. The results of the calculations showed that the presence of the mechanical and termal heterogeneity on a flat surface of the heterogeneous layer leads to attenuation of the incident UV within this layer, depending on the volume fraction of the condensed phase. A sufficiently large volume fraction of the solid phase in the heterogeneous layer leads to the deformation of the head front to the SW and strong braking, with the result that reflections from the wall hardly occurs and the leg of the Mach is not formed, however, when a2 > 5 • 10-2 - it is formed over the surface of the heterogeneous layer. With the increase of the gas temperature in the heterogeneous layer the velocity (SW) increases while the intensity of its amplitude falls. The calculations given in this particular work showed that the increasing of the gas temperature in the layer promotes the increasing of the moment force acting on the lattice of solid substances. At high temperatures of the gas in the heterogeneous layer, despite power and heat interaction of gas and substances the precursor is formed.
Keywords: numerical procedure; mathematical model; conservation laws; Courant number.
References
1. Ivandaev A.I, Kutushev A.G. [The Effect of the Shielding Layers of the Gas Suspension on the Reflection of Shock Waves]. Applied Mechanics and Technical Physics, 1985, no. 21, pp. 115-120. (in Russian)
2. Kruglikov B.S., Kugushev A.G. [The Attenuation of Air Shock Waves by Layers of Dusty Gas and Lattices]. Applied Mechanics and Technical Physics, 1988, no. 1, pp. 51-57. (in Russian)
3. Kruglikov B.S., Kutushev A.G. [The Weakening Of Shock Waves Louvers]. Physics of Combustion and Explosion, 1988, no. 1, pp. 115-117. (in Russian)
4. Ivandaev A.I., Kutushev A.G., Nigmatulin R.I. [Gas Dynamics of Multiphase Media. Shock and Detonation Waves in Gas Suspensions]. Itogi nauki i tehniki. Seriya: Mehanika zhidkosti i gaza, 1981, vol. 16, pp. 209-287. (in Russian)
5. Kutushev A.G., Rudakov D.A. [Numerical Investigation of the Effect of a Shock Wave on an Obstacle Screened by a Layer of Porous Powdery Medium]. Applied Mechanics and Technical Physics, 1993, no. 25, pp. 25-31. (in Russian)
6. Kovalev Yu.M., Cheremohov A.Yu. [The Attenuation of Air Shock Waves by a System of Grids]. Voprosy atomnoy nauki i tekhniki. Series: Mathematical Modelling of Physical Processes, 1997, no. 3, pp. 39-43. (in Russian)
7. Belyaev P.E., Klinacheva N.L. Impact of Gas Suspension Shielding Layer on The Force Effect of Shock Waves on a Rigid Wall. Bulletin of the South Ural State University. Series: Mathematics. Mechanics. Physics, 2016, vol. 8, no. 4, pp. 49-55. DOI: 10.14529/mmphl60406 (in Russian)
8. Klinacheva N.L., Kovalev Yu.M., Kuropatenko V.F. Modified Mathematical Model of a "Frozen" Gas Suspension. Journal of Engineering Physics and Termophysics, 2014, vol. 87, no. 6, pp. 1456-1462. DOI: 10.1007/sl0891-014-1150-x
9. Kovalev Yu.M., Kovaleva E.A. [Mathematical Analysis of the Conservation Equations of Two-Phase Mixtures]. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2014, vol. 7, no. 2, pp. 29-37. DOI: 10.14529/mmpl40202 (in Russian)
10. Andrushhenko V.A., Meshcheryakov M.V., Chudov L.A. [The Reflection of Spherical Shock Wave from the Plane When It Has a Layer of Heated GasJ. Fluid Mechanics, 1990, no. 3, pp. 141-147. (in Russian)
11. Andrushhenko V.A., Meshcheryakov M.V. [The Interaction of Spherical Shock Waves with Near-Surface Thermal Gas InhomogeneitiesJ. Physics of Combustion and Explosion, 1990, no. 1, pp. 77-82. (in Russian)
12. Antonov V.A., Grishin A.M., Kovalev Yu.M., Naymushina L.Yu. [Modelling of Explosion of a Cord Charge in the Forest Canopy in the Absence of FireJ. Physics of Combustion and Explosion, 1993, no. 4, pp. 115-123. (in Russian)
13. Belocerkovskiy O.M., Davydov Yu.M. Metod krupnyh chastic v gazovoj dinamike [The Method of Large Particles in Gas Dynamics], Moscow, Nauka, 1982.
14. Belocerkovskiy O.M., Davydov Yu.M. [Unsteady "Method of Large Particles" for Gas Dynamics Calculations]. Applied Mechanics and Technical Physics, 1971, no. 21, pp. 182-207. (in Russian)
15. Grishin Yu.A., Zenkin V.A. Improving the Stability of the Numerical Algorithm of the Method of Large Particles. Journal of Science and Education, 2011, no. 13, 5 p. Available at: http://old.technomag.edu.ru/doc/221488.html
16. Kovalev Yu.M., Kovaleva E.A., Pigasov E.E. The Analysis of Some Modifications of the Large-Particles Method on the Basis of Research of Gas-Suspension Currents. Bulletin of the South Ural State University. Series: Mathematics. Mechanics. Physics, 2015, vol. 7, no. 3, pp. 71-77. (in Russian)
Received August 25, 2017
