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

НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ май-июнь 2015 Том 15 № 3 ISSN 2226-1494 http://ntv.ifmo.ru/

SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS

May-June 2015

Vol. 15 No 3 ISSN 2226-1494

http://ntv.ifmo.ru/en

УДК 532.529

ИНТЕРФЕРЕНЦИЯ СКАЧКОВ УПЛОТНЕНИЯ ОДНОГО НАПРАВЛЕНИЯ

П.В. Булат3, П.В. Денисенкоь

a Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация b Университет Уорвика, Ковентри, CV4 7AL, Великобритания Адрес для переписки: pavelbulat@mail.ru Информация о статье

Поступила в редакцию 03.04.15, принята к печати 14.04.15

doi:10.17586/2226-1494-2015-15-3-500-508

Язык статьи - русский

Ссылка для цитирования: Булат П.В., Денисенко П.В. Интерференция скачков уплотнения одного направления // Научно-технический вестник информационных технологий, механики и оптики. 2015. Т. 15. № 3. С. 500-508.

Аннотация

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

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

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

скачок уплотнения, ударная волна, интерференция скачков уплотнения, скачки уплотнения одного направления, догоняющие скачки уплотнения.

INTERFERENCE OF UNIDIRECTIONAL SHOCK WAVES P.V. Bulat a, P.V. Denissenko b

a ITMO University, Saint Petersburg, 197101, Russian Federation b University of Warwick, Coventry, CV4 7AL, United Kingdom Corresponding author: pavelbulat@mail.ru Article info

Received 03.04.15, accepted 14.04.15 doi:10.17586/2226-1494-2015-15-3-500-508 Article in Russian

For citation: Bulat P.V., Denissenko P.V. Interference of unidirectional shock waves. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2015, vol.15, no. 3, pp. 500-508.

Abstract

Subject of study. We consider interference of unidirectional shock waves or, as they are called, catching up shock waves. The scope of work is to give a classification of the shock-wave structures that arise in this type of interaction of shock waves, and the area of their existence. Intersection of unidirectional shock waves results in arising of a shock-wave structure at the intersection point, which contains the main shock wave, tangential discontinuity and one more reflected gas-dynamic

discontinuity of unknown beforehand type. The problem of determining the type of reflected discontinuity is the main problem that one has to solve in the study of catching shock waves interference.

Main results. The paper presents the pictures of shock-wave structures arising at the interaction of catching up shock waves. The areas with a regular and irregular unidirectional interaction of shocks are described. Characteristic shock-wave structures are of greatest interest, where reflected gas-dynamic discontinuity degenerates into discontinuous characteristics. Such structures have a number of extreme properties. We have found the areas of existence for such shock-wave structures. There are also areas in which the steady-state solution is not available. The latter has determined revival of interest for the theoretical study of the problem, because the facts of sudden shock-wave structure destruction inside the air intake of supersonic aircrafts at high Mach numbers have been discovered.

Practical significance. The theory of interference for unidirectional shock waves and design procedure are usable in the design of supersonic air intakes. It is also relevant for application possibility investigation of catching up oblique shock waves to create overcompressed detonation in perspective detonation air-jet and rocket engines. Keywords

discontinuity, shock wave, interference of shocks, unidirectional shocks, catching up shock waves.

Введение

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

Наиболее просто представить себе взаимодействие газодинамических разрывов (ГДР) одного направления в одномерном случае, когда одна ударная волна догоняет другую (рис. 1, a). По аналогии с одномерным случаем пересекающиеся косые скачки уплотнения одного направления, т.е. разворачивающие поток в одном и том же направлении, называются «догоняющими» (рис. 1, б). Таких скачков в УВС может быть несколько: чем их больше, тем ближе такая конфигурация к центрированной изоэнтропической волне сжатия (рис. 1, в), которая является идеальной структурой для сжатия сверхзвукового потока и широко используется при проектировании гиперзвуковых воздухозаборников [1, 2]. Первые работы по изучению догоняющих скачков уплотнения (ДСУ), появились в 50-е и 60-е годы XX века и были связаны с появлением двух- и трехскачковых регулируемых воздухозаборников сверхзвуковых самолетов [3-5].

Di

D2

Рис. 1. Интерференция разрывов одного направления: догоняющие одномерные ударные волны й. и й2, стрелками показано направление движения (а); косые догоняющие скачки уплотнения (б); центрированная изоэнтропическая волна сжатия (в). М - число Маха; Т - точка пересечения скачков уплотнения; И - отраженный разрыв; V.,, v2 - разрывные характеристики, ограничивающие центрированную волну сжатия, а., а2 - догоняющие скачки уплотнения; а3 - главный скачок уплотнения; т - тангенциальный разрыв; шс - волна сжатия. Стрелками показаны линии тока - левая ^ или правая ^ стрелка над обозначением газодинамического разрыва означает соответственно левое или правое

направление этого разрыва

а

в

Большой вклад в развитие теории ДСУ в выявление областей существования различных ударно -волновых структур с их участием внесли В.Н. Усков [6] и А.Л. Старых [7]. А.В. Омельченко и В.Н. Усковым было исследовано на экстремум поведение газодинамических переменных за УВС с косыми ДСУ [8]. В.Н. Усковым в полном объеме решена задача о взаимодействии одномерных ударных волн одного направления [9], в том числе, отвечающих некоторому критерию оптимальности [10]. К началу 2000-х годов казалось, что данное явление изучено практически полностью. Однако многие сведения оставались фрагментарными, например, области существования различных случаев интерференции ДСУ были исследованы только для некоторых чисел Маха. Отсутствовала и теория перехода от УВС одного типа к другому. Вернуться к ДСУ заставило выявление нестационарных процессов при обтекании головных частей с выдвинутой в поток иглой (рис. 2, а), а также изучение так называемого парадокса Неймана [11]. Суть его состоит в том, что при числах Маха, близких к единице, тройные конфигурации ударных волн теоретически существовать не могут [12], но, тем не менее, они наблюдаются визуально (рис. 2, б).

Рис. 2. Парадокс Неймана: догоняющие скачки при обтекании головной части летательного аппарата с выдвинутой в поток иглой (а); обтекание профиля трансзвуковым потоком (б)

Особенно остро данный вопрос встал при изучении истечения сверхзвуковых струй из сопел с числами Маха меньше 1,5 (рис. 3). Именно такие сопла используются сегодня на сверхзвуковых военных самолетах. Проблема заключается в том, что в данном случае отраженный скачок является, строго говоря, приходящим, а не исходящим. Следовательно, отражение от оси симметрии не может быть причиной его возникновения. Непонятны и механизмы трансформации УВС при изменении числа Маха в этом диапазоне. В то же время известно, что при небольших числах Маха сверхзвуковые струи характеризуются повышенным уровнем акустического излучения, что указывает на нестационарный характер течения. Все эти вопросы требуют детального изучения.

Рис. 3. Отражение скачков от оси симметрии сверхзвуковой струи, истекающей из сопла,

с числом Маха, близким к единице

Математическая модель

Регулярная интерференция. Если интенсивность скачков 1, 2 (рис. 1) и число Маха М таковы, что течение за вторым скачком - сверхзвуковое, то интерференция ДСУ носит регулярный характер. Образуется одна из трех возможных УВС, показанных на рис. 4. Основная задача при изучении регулярного взаимодействия ДСУ - определение типа отраженного разрыва: скачок (рис. 4, а), разрывная характеристика (рис. 4, б), центрированная волна разрежения (рис. 4, в). Задача может быть решена геометрически на плоскости поляр. Исходная поляра соответствует числу Маха набегающего потока. Если на этой поляре выбрать точку 1 такую, что угол поворота потока в и интенсивность скачка 3 будут соответственно

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

Л

Рис. 4. Ударно-волновые структуры и решения на плоскости поляр при регулярном взаимодействии догоняющих скачков уплотнения: случай, когда отраженный разрыв - скачок уплотнения а4 (а); случай, когда отраженный скачок - разрывная характеристика V,, (б); случай, когда отраженный разрыв - волна разрежения иц (в). 1 - точка на поляре, построенной для числа Маха М набегающего потока, соответствующая интенсивности скачка а1; 2 - точка на поляре, построенной для числа Маха потока за скачком а1, соответствующая интенсивности скачка а2; 3-4 - точка пересечения поляры, построенной по числу Маха потока за скачком а2, с полярой, построенной по числу М, соответствующая интенсивности главного скачка а3 и интенсивности отраженного разрыва. Л=1пЛ - натуральный логарифм интенсивности

скачка

На этой вторичной поляре точка 2 определяет интенсивность и угол разворота потока на втором догоняющем скачке. Если вторая поляра расположена вне основной поляры, то отраженный разрыв будет скачком уплотнения, с параметрами в точке 3-4 на рис. 4, а. Если вторичная поляра проходит внутри основной, то отраженный разрыв - волна разрежения, интенсивность которой также определяется точкой 3-4 на плоскости поляр на рис. 4, в. Эти два случая разделяет характеристическая УВС (рис. 4, б), когда отраженный разрыв вырождается в разрывную характеристику. Эта УВС соответствует тройной конфигурации ударных волн типа ТК-3 [13].

Интенсивность отраженного разрыва определяется из следующих уравнений:

при 53 :р„(М,3А) + р,(М2,03А) = Р,(М,3,) + р,(М„32); (1)

при 0О3: (М,Л)-Рт(ММ2,03л) = Р,(М,3,) + р,(М,,32); (2)

где G = (3,32)-1; 31 е [1;35]; 32 е [1;321 ], М1 - число Маха за г'-м скачком уплотнения.

Характеристическая УВС. Найдем условия существования характеристической УВС, разделяющей случаи с отраженным разрывом - скачком уплотнения и отраженным разрывом - волной разрежения. Данная задача по технике решения напоминает анализ УВС, возникающих в центре волны сжатия [14].

Очевидно, что в точке на главной поляре, соответствующей характеристической УВС, главная и вторичные поляры должны иметь одинаковый угол наклона, т.е. порядок касания поляр - не ниже первого.

Найдем условия, при которых вторичная поляра располагается внутри основной поляры и вне ее. Вторичная поляра Р5 (МТ1) может находиться внутри поляры Р5 (М) или вне ее, в зависимости от значений производных Л' = (ЭЛ / Зр) в текущей точке 1 основной поляры и Л' = (ЭЛ / ЭР) в начале координат вторичной поляры, исходящей из этой точки. Так как (3 + е) [уМ2 - (3 -1)] - (1 + е3)(3 -1)

Л'(М.) =

[3(1 + е3)[(3 + е)(М2 -1) - (32 -1)]]''

(3)

в

где е = (у-1)/(у +1), у - показатель адиабаты, а Л' определяется формулой

М 3 3", (4)

Эр ^ Э3)

то, приравнивая производные (3) и (4), получаем кубическое уравнение относительно технической переменной х = 3{ - е , где 3 - интенсивность скачка в точке на поляре, в которой отраженный разрыв 4 вырождается в характеристику (рис. 4, б), со следующими коэффициентами:

X А„х" = 0:

и=0

А3 = 1; А2 =[(у -1)2 - 2]М2 + 2у - 3е; (5)

А = (2 - у)уМ4 + 2 [(у -1)е -1]М2 +1 + 4е + 3е2;

А0 = -е(1 + е + уМ2)2.

Такие точки называются характеристическими. На поляре может быть несколько характеристических точек. Уравнение (5), записанное относительно М2, является квадратным:

аМ4 + ЬМ2 + с = 0;

а = у [(2 - у)(3, + е) - уе] ; (6)

Ь =[(у-1)2 - 2](( + е)2 + 2[(у-1)е -1]] + е) - 2уе(1 + е);

с = (^ + е)3 + (2у - 3е)(^ + е) + (2е2 + 4е +1)(1 + е).

Дискриминант Д(у; М) уравнения (6) обращается в ноль при М = Мд (у). Если М < Мд , то на Ра (М) может иметь место лишь одна характеристическая точка /, если М > Мд , то из трех действительных корней уравнения необходимо отбросить те, которые не удовлетворяют неравенствам 1 < 3/ < 3т. Полагая в (6) = 1, получим граничные числа Маха, определяющие область существования характеристических точек на поляре:

(5 - 3у)М4 - 4(3 - у)М2 + 8 = 0. (7)

Решая уравнение (7), получим два корня:

МА = (а, -Ь,)1/2; М/2 = (а, + Ь,)1/2, (8)

где

а1 = 23-1) ; ь. = ^ЁЕ!. (9)

1 5 - 3у 1 5 - 3у

Приведенная математическая модель может быть использована для нахождения областей существования УВС с различными типами отраженных разрывов. С помощью (8)-(9) можно построить кривые М (у) и (М,у), соответствующие характеристическим УВС, они будут разделять области с различными типами отраженного разрыва.

Результаты и анализ

Классификация УВС, возникающих при нерегулярном взаимодействии ДСУ. Алгоритм решения задачи об интерференции ДСУ в случае их регулярного взаимодействия достаточно очевиден. Нужно только проанализировать расположение на главной поляре точки, соответствующей вторичной поляре, относительно характеристических точек (рис. 4). Для этого используют уравнения (5)-(8). Гораздо сложнее обстоит дело, если решение уравнений (1)-(2) отсутствует. На приходящих скачках образуются тройные точки и интерференция становится нерегулярной. На рис. 5 рассмотрены случаи, когда течение за скачком 2 - дозвуковое. УВС, показанная на рис. 5, а, относится к тройным конфигурациям ТК-3. Точ -ка 2, соответствующая второму скачку, расположена на вторичной поляре выше звуковой точки 5, в которой число Маха за скачком равно 1, а интенсивность скачка Второй догоняющий скачок искривляется. Его интенсивность изменяется от параметров, соответствующих точке 2, до параметров в точке 3-4, в которой давление и направления векторов скорости по обе стороны тангенциального разрыва должны быть одинаковыми.

В зависимости от числа М и интенсивностей скачков 1 и 2 второй скачок может искривляться настолько сильно, что поменяет свое направление (рис. 5, б). В этом случае УВС относится к классу тройных конфигураций ТК-2, а точка 3-4, соответствующая параметрам в тройной точке Т, лежит на левой ветви вторичной поляры. Если течение за скачком 4 в тройной точке Т - сверхзвуковое, а за вторым приходящим скачком - дозвуковое, то образуется УВС, показанная на рис. 5, в. Механизм перехода от УВС,

изображенной на рис. 5, б, в УВС на рис. 5, в, и обратно, не изучен.

Л 3-4/Т-\ 2

Л

М2<1

М2<1

^

СТ2 б

Л

3-4

5-6 2

М2<1

Рис. 5. Ударно-волновые структуры и решения на плоскости поляр при нерегулярном взаимодействии догоняющих скачков уплотнения при дозвуковом течении за вторым приходящим скачком: случай, соответствующий образованию ТК-3 (а); случай, когда скачок а2 искривляется настолько сильно, что меняет свое направление (б); маховское взаимодействие догоняющих скачков (в); б - звуковая точка

на поляре, соответствующая М=1 за скачком а2

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

При увеличении числа М течение за скачком 4 может стать сверхзвуковым. Тогда для скачка 2 скачок 4 становится встречным, и образуется УВС, характерная для интерференции встречных скачков (рис. 6, б).

Исходящий из точки Т скачок 4 является встречным для скачка 2 и приходящим для УВС, состоящей из двух тройных точек 5 и 8. При дальнейшем увеличении М взаимодействие скачков становится регулярным, точки 5 и 8 сливаются в одну.

3-4

___т

2<1

Т

а б

Рис. 6. Ударно-волновые структуры и решения на плоскости поляр при нерегулярном взаимодействии догоняющих скачков уплотнения при сверхзвуковом течении за вторым приходящим скачком

Области существования характеристических УВС. Анализ производных Л' и Л' в различных областях существования характеристических точек Л позволяет сделать вывод, что при

Т

Т

а

в

"л' > л;

г 1 г 1 (10)

_ ¿1 4 Л М Л }

начало первой поляры Ра (М^) проходит внутри Ра (М) (рис. 4, в)

"л' < л;

г п г п (11)

_Ф>МЛ\

вне основной поляры (рис. 4, а). Таким образом, условия (10)-(11) определяют тип отраженного разрыва в УВС с ДСУ при интенсивностях приходящего скачка уплотнения, близких к характеристическим.

Для у=5/4 в интервале М е [м^,М^ ^ находится только одна характеристическая точкаа при М > М^ - две: Л и/2. При у=5/4 существует область Мд < М < М^ , в которой поляры сжатия имеют все

три характеристические точки (рис. 7).

В диапазоне М > М^ нижние точки Л отсутствуют для любых значений у . Особые числа Маха

сведены в таблицу.

Y 1,1 1,25 1,4 1,67

M д 1,62 2 2,46 4,125

MA 1,302 1,265 1,245 1,224

Mf J2 1,666 2 2,54 да

Таблица. Особые числа Маха при разных показателях адиабаты

График, отражающий область существования характеристических УВС, показан на рис. 7. Хорошо видно, что при различных числах Маха на ударной поляре может существовать от одной до трех характеристических точек.

М

4

2

III / /

II у/ /

Mf2 >> /У

M7l / I

3

1 1,1 1,2 1,3 1,4 1,5 1,6 у

Рис. 7. Зависимость характерных чисел Маха от показателя адиабаты. I - одна характеристическая точка £,, II - две характеристические точки f и f2, III - три характеристические точки f.,, f2, f3

На рис. 8 показаны зависимости от числа М интенсивности скачка 1 в системе догоняющих скачков, при которой образуются характеристические УВС. Именно такие УВС являются оптимальными по критерию степени восстановления полного давления. Во врезке «А» показаны в увеличенном масштабе кривые, соответствующие характеристическим точкам f1. Пунктирная дискриминантная кривая, соответствующая равенству нулю дискриминанта в уравнении (6), разделяет характеристические кривые на две ветви: верхнюю f и нижнюю f3.

Числа Маха Mк ограничивают слева область существования характеристических УВСf и f3. При у<1,2 дискриминантная кривая уходит в область J<1, следовательно, существует только верхняя ветвь f характеристической кривой. Таким образом, в зависимости от числа Маха и показателя адиабаты характеристические УВС могут возникать в одной, двух и трех точках на ударной поляре.

Рис. 8. Характеристические точки на ударной поляре. Сплошные линии - характеристические кривые, соответствующие точкам f,, f2, f3; пунктирная линия - дискриминантная кривая, соответствующая условию Д=0

Заключение

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

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

Литература

1. Bulat M.P., Bulat P.V. The analysis centric isentropic compression waves // World Applied Sciences Journal. 2013. V. 27. N 8. P. 1023-1026. doi: 10.5829/idosi.wasj.2013.27.08.13719

2. Усков В.Н., Булат П.В. О задаче проектирования идеального диффузора для сжатия сверхзвукового потока // Фундаментальные исследования. 2012. № 6-1. C. 178-184.

3. Вуст В. К теории разветвления скачков уплотнения / Газовая динамика. М.: Изд-во иностр. лит., 1950. С. 131-143.

4. Веккен Ф. Предельные положения вилкообразных скачков уплотнения / В сб. Механика. 1950. № 4. С. 24-34.

5. Росляков Г.С. Взаимодействие плоских скачков уплотнения одного направления // Численные методы в газовой динамике. М.: МГУ, 1965. С. 28-51.

6. Адрианов А.Л., Старых А.Л., Усков В.Н. Интерференция стационарных газодинамических разрывов. Новосибирск: Наука, 1995. 180 c.

7. Росляков Г.С., Старых А.Л., Усков В.Н. Интерференция стационарных скачков уплотнения одного направления // Изв. АН СССР. Механика жидкости и газа. 1987. № 4. С. 143-152.

8. Омельченко А.В., Усков В.Н. Оптимальные догоняющие скачки уплотнения с ограничениями на суммарный угол поворота потока // Прикладная механика и техническая физика. 1999. Т. 40. № 4. С. 99-108.

9. Усков В.Н. Бегущие одномерные волны. СПб.: БГТУ «ВОЕНМЕХ», 2000. 224 с.

10. Усков В.Н. Оптимальные бегущие по потоку газа одномерные волны // XV Сессия Международной школы по моделям механики сплошной среды. 2000. С. 63-78.

11. Vasilev E., Elperin T., Ben-Dor G. Analytical reconsideration of the von Neumann paradox in the reflection of a shock wave over a wedge // Physics of Fluids. 2008. V. 20. N 4. Art. 046101. doi: 10.1063/1.2896286

12. Усков В.Н., Чернышев М.В. Особые и экстремальные тройные конфигурации скачков уплотнения // Прикладная механика и техническая физика. 2006. Т. 47. № 4. С. 39-53.

13. Усков В.Н., Чернышев М.В. Тройные конфигурации стационарных ударных волн // XII Всероссийская научно-практическая конференция «Актуальные проблемы защиты и безопасности». Санкт-Петербург, 2009.

14. Bulat P.V., Bulat M.P. Discontinuity of gas-dynamic variables in the center of the compression wave // Research Journal of Applied Sciences, Engineering and Technology. 2014. V. 8. N 23. P. 2343-2349.

Булат Павел Викторович — кандидат физико-математических наук, кандидат экономических

наук, руководитель научной лаборатории, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, pavelbulat@mail.ru

Денисенко Петр Викторович — кандидат физико-математических наук, научный сотрудник,

Университет Уорвика, Ковентри, СУ4 7ЛЬ, Великобритания, p.denissenko@gmail.com

Pavel V Bulat

Petr V Denissenko

PhD in physics and mathematics, PhD in economics, Scientific Supervisor of International Laboratory, ITMO University, Saint Petersburg, 197101, Russian Federation, pavelbulat@mail.ru PhD, researcher, University of Warwick, Coventry, CV4 7AL, United Kingdom, p.denissenko@gmail.com