Теоретические основы инженерного конструирования бустерных фрикционно-вихревых турбин приводов насосов аварийной подпитки для АЭС Текст научной статьи по специальности «Механика и машиностроение»

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

УДК 621.039:621.311

Погосов Олексій Юрійович, д-р техн. наук, професор, Одеський національний політехнічний університет, м. Одеса, Україна

Дерев’янко Ольга Володимирівна, аспірант, Одеський національний політехнічний університет, м. Одеса, Україна. ov derevyanko@mail.ru

ТЕОРЕТИЧНІ ОСНОВИ ІНЖЕНЕРНОГО КОНСТРУЮВАННЯ БУСТЕРНИХ ФРИКЦІЙНО-ВИХРОВИХ ТУРБІН ПРИВОДІВ НАСОСІВ АВАРІЙНОГО ПІДЖИВЛЕННЯ ДЛЯ АЭС

Розглянуто теоретично обгрунтований підхід до проектування вдосконалених бустерних фрикційно-вихрових турбін, заснований на аналізі динаміки робочого тіла в фрикційно-вихрових турбінах приводних агрегатів методами математичної фізики. На основі теоретичних розробок представлений варіант вирішення завдання інженерного конструювання комбінованих турбоприводу насосів аварійного підживлення обладнання АЕС.

Ключові слова: безпека атомної енергетики, турбоприводу насосів аварійного підживлення для АЕС, фрикційно-вихрова турбіна.

Погосов Алексей Юрьевич, д-р техн. наук, профессор, Одесский национальный политехнический университет, г. Одесса, Украина, pogosov@ukr.net

Деревянко Ольга Владимировна, аспирант, Одесский национальный политехнический университет, г. Одесса, Украина, ov derevyanko@mail.ru

ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ИНЖЕНЕРНОГО КОНСТРУИРОВАНИЯ БУСТЕРНЫХ ФРИКЦИОННОВИХРЕВЫХ ТУРБИН ПРИВОДОВ НАСОСОВ АВАРИЙНОЙ ПОДПИТКИ ДЛЯ АЭС

Рассмотрен теоретически обоснованный подход к проектированию усовершенствованных бустерных фрикционно-вихревых турбин, основанный на анализе динамики рабочего тела в фрикционно-вихревых турбинах приводных агрегатов методами математической физики. На основе теоретических разработок представлен вариант решения задачи инженерного конструирования комбинированных турбоприводов насосов аварийной подпитки оборудования АЭС.

Ключевые слова: безопасность атомной энергетики, турбоприводы насосов аварийной подпитки для АЭС, фрикционно-вихревая турбина.

Pogosov Aleksandr Yurevich , Dr. Eng. Sc, professor, Odessa national polytechnic university, Odessa, Ukraine. Derevianko Olga Vladimirovna, post-graduate student, Odessa national polytechnic university, Odessa, Ukraine.

ov derevyanko@mail.ru

THEORETICAL FOUNDATIONS OF ENGINEERING DESIGN BOOSTER FRICTION-VORTEX TURBINE DRIVES OF PUMPS FOR EMERGENCY IRRIGATION FOR NPP

Considered theoretically grounded approach to the design of improved booster friction-vortex turbines, based on the analysis of the dynamics of the working fluid in the friction-vortex turbine drive units methods of mathematical physics. Based on the theoretical developments presented variant of the solution of engineering design combined turbine drive emergency feed water pumps NPP equipment.

Keywords: safety of nuclear energy, turbine pumps for nuclear power plant emergency makeup, friction-vortex turbine.

Введение

При аварийном электрообесточивании энергоблоков АЭС обеспечение их безопасности во многом определяется надежным восполнением дефицита водных сред при помощи насосов подпитки. Рациональное конструирование фрикционно-вихревых турбин (турбин трения) в качестве бустерных конструктивных элементов комбинированных турбоприводов насосов аварийной подпитки тепломассообменного оборудования АЭС требует предварительного аналитического рассмотрения соответствующих физических процессов для обоснованных инженерных решений при профилировании рабочих колес роторов этих агрегатов [1].

До последнего времени физика вихревого движения потоков несепарированного влажного пара, выработка которого технологически свойственна АЭС (и который при аварийном электрообесточивании энергоблока может быть использован в таких устройствах в качестве рабочего тела) не была описана применительно к таким инженерным задачам [2].

18

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

Основная часть

Целью исследования является разработка аналитически обоснованного подхода к проектированию усовершенствованных бустерных фрикционно-вихревых турбин и использование его результатов для инженерного конструирования турбоприводов насосов аварийной подпитки оборудования АЭС.

Для достижения поставленной цели, прежде всего проведем анализ динамики рабочего тела в фрикционно-вихревых турбинах приводных агрегатов методами математической физики, как это принято в классическом подходе А. В. Борисова и С. М. Рамоданова при решении задач по теории вихрей, затем представим вариант инженерного решения целевой задачи. В качестве исходного положения примем, что начало организованного вихревого движения в цилиндрическом корпусе турбины всегда связано с предшествующим безвихревым движением среды, поступающей через сопловой аппарат. При безвихревом движении рабочего тела на входе в рабочее колесо, сила R и момент силы M, действующие на твердое тело - в данном случае ротор турбины в ортоаксиальном плоском сечении -традиционно, в задачах математической физики, определяются формулами Седова [3]:

Здесь X, Y - проекции силы на неподвижные оси, р - плотность жидкости.

При этом момент силы вычисляется относительно начала неподвижной системы координат. В нашем случае эти формулы нуждаються в дальнейшем уточнении. Следуя подходу, описаному в [4], далее рассмотрим геометрию формирования вихря рабочего тела, представленную на рис.1. Пусть, Q и С2 - концентрические окружности, определяющие

пределы вихревого движения, центр которых совпадает с точкой, вокруг которой в данный момент концентрируется вихрь.

Рис. 1. Физико-геометрическое представление движения и взаимодействия рабочего тела с ротором фрикционно-вихревой турбины

Пусть радиус r2 окружности С2 настолько велик и радиус r1 окружности С1

настолько мал, что эти окружности не имеют общих точек с контуром С формирования дополнительных вихрей [5]. Пользуясь принятым физико-геометрическим представлением,

(1)

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

19

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

проведем разрезы AB и A1B1 так, чтобы в области D, ограниченной контурами C, C1, C2, AB и A1B1, потенциал обтекания ротора был однозначной функцией. Импульс фрагмента рабочего тела, находящейся в области D, представим в комплексной форме:

Р = Px + iPy =p\^~dT+ *р\^&dT , (2)

D Sx D дУ

где Px, Py - проекции импульса на неподвижные оси. Применяя формулу Грина и формулу интегрирования по частям, получим:

Р = Рх + iPy = pI (pdy + ipI (pdx = ip I zdp. (3)

dD dD dD

На двух сторонах разрезов AB и A1B1 значение p отличаются на постоянную величину Г и Г * - Г , следовательно, интегралы, соответствующие сторонам разрезов, упраздняются:

Р = Px + iPy = ipI zdp-ip zdp-ipzdp . (4)

C2 C1 C

Для фрагмента рабочего тела, находящегося в области D, может быть записано выражение для изменения импульса фрагмента рабочего тела (фрагмента вихря):

dP_

dt

ipd | zdp-i pd-1zdp-i pd-1zdp = -RC -RCi + RC2,

dt

dt

(5)

где RC , RC1 , RC2 - главные векторы сил, вызываемых давлением рабочего тела, действующие на контурыC, C1, C2, соответственно. Можно принять, что разность:

RC2- ipdp I zdP

dt

C2

стремится к нулю при r2 ^ да. Действительно, давление определяется по формуле Коши-Лагранжа [6]:

Тогда

Учитывая, что:

p = -pSP-pA + F (,). St 2

RC2 = i| pdz = -i| p^S-dz -ipI~dz .

C2 C2 St C2 2

d I zdp=I dtzdp+1 zdd~=I (v+ioy Ур-1 id+(V+V )dz,

d

, dp

dt c c dt c dt c

dp

dt

(6)

(7)

(8)

в результате получаем:

Rc2 -ipd I zdp = p I^xVydx -I

2 2 ^ f

Vx -v„

-dy

V C2

+ ip

IVxVdy +I

V2-V2 -

----dx

V C2

(9)

20

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

г і ^

В силу того, что v = O —т при r2 ^ да, приходим к выводу, что оба интеграла в правой

V r J

части стремятся к нулю при r1 ^да. Найдем предел разности:

d г

RC - ip— I zdp при r ^ 0.

dt c

C1

Поместим начало выбранной условно неподвижной системы координат в центр окружности С1. При фиксированном t справедливо представление [7]:

(

Vx + ivy =

У

\ (

2— х2 + у2

+ Vx

+ І

-Г

2— х2 + у2

+ vx

(10)

Здесь Vx, vy - аналитические в окрестности начала координат функции переменных х у . Обозначим Vx (0,0) = u0, vy (0,0) = v0 . Согласно (9) имеем:

Rq - ipd I zdp = ——2 I (-V0у + u(}x)dx - I (у + vQx) dy

-Гр

dt.

2—r

V C2

іГр

2—r2

| (-Vу + U0x)dx - | ((у + U0x)у

C2

\

(11)

V C2

+ ...

Троеточие в (11) указывает на наличие слагаемых, стремящихся к нулю при (x, у) ^ 0, которыми можно пренебречь. Таким образом, искомый предел равен Гpv0 - іГри0. Переходя в (5) к предельным случаям r2 ^+да и r1 ^ 0, получаем выражение для силы, действующей на контур C :

RC = ip—— |zdp - Гpv0 + іГpu0.

dt

(12)

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

Обратимся теперь к определению момента силы, действующего на контур C . Кинетический момент (момент импульса) фрагмента вихря рабочего тела, находящейся в области D, равен [8]:

K=р\\

d(xP) д(ур)

ду dx

\ r 2 r 2 r 2

dr = -p| p(dx + ydy) = PI —dp-pj —dp-pj —dp . (13)

Линии действия сил давления на контурах С1 и С2 проходят через начало координат, поэтому их момент относительно начала координат равен нулю. Момент силы, действующий на контур С, имеет вид M = -dK/dt. Исследуем производные по времени от интегралов в правой части (13). Равенство предела нулю:

r

У d

lim — р

2 ^+да dt

(14)

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

21

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

следует из того, что модуль скорости фрагмента рабочего тела убывает на бесконечности как —. Для интеграла по контуру C1 имеем:

dtp\L2d(P = p\ d r—\d<p = p\(xul + yoyox )dy .

d ^p >

dt' C 2

C1

C1 dt I 2 J

(15)

2

Можно понять, что с учетом (11), что правая часть (15) обращается в ноль при r1 ^ 0. Продифференцировав (13) и перейдя к пределам r2 ^+да и r1 ^ 0, получим выражение для момента силы, действующего на контур C:

м = pdL[ Д d„ (16)

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

Получим интеграл для силы, входящий в выражение (12). Конформное отображение

внешности единичной окружности C плоскости С на внешность C имеет вид:

z - f (С)-kC +ko + КС 1 +... , (17)

> 0. С учетом этого равенства, контурный интеграл в (12) может быть преобразован следующим образом:

, dz где k ----

dC

Jzd?e - J( + ко + +...)I т—(g(С-h))-arg(-^l1) + -2— arg^

ггч Г*

2a

-Г

2ai

-Г

J

2ai

kC + k0 + k1C

СГ"

f (9-к'

- dc-J

kC + k0 + ktC 1 + ... dC

c-Vh

2a

Г fkC + k0 + kC 1 + ••• d^ - (18)

-J-

2—

C-1 h

2ai

2—k l

1J

— a—k0 --rf (l1 ) + r*k0 2a

C

Формально заменяя в формуле (18) р на (рв, мы фактически определяем ту часть силы,

действующей со стороны фрагмента рабочего тела на ротор турбины, которая обусловлена присутствием вихря и соответствующей циркуляцией [9].

Таким образом, согласно (12) и (18), для силы, обусловленной циркуляцией и наличием вихря, получаем следующее выражение:

R - 1Г*Р^К irPdtf (ll) - ГР°~0 +іГPU0 (19)

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

22

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

исполнению ротора турбины, как это буде показано далее.

Опираясь на аналитические результаты выполненного математического моделирования, позволившего уяснить качественную физическую картину

рассматриваемого явления, далее исследуем возможности инженерной реализации рациональных конструкций турбин. Очевидно, что возможность решения инженерной задачи для «плоского случая» в замкнутом пространстве, ограниченном цилиндрическим корпусом, в котором должен размещаться ротор фрикционно-вихревой бустерной турбины, технически связана с формированием рабочего колеса в виде пакета плоско-параллельных дисков.

Если рабочее колесо ротора, собранно по типу рабочего колеса турбины Теслы, то для организации плоского вихревого движения среды в междисковом пространстве от периферии к центру каждый диск должен иметь в аксиальной части отверстия для выхода отработавшего рабочего тела, как показано на рис. 2. При таком инженерном решении организованный вихрь движения частиц среды геометрически может представлять собой пук скручивающихся спиралей [10]. Организация вихря технически может осуществляться за счет соответствующих спиральных проточек. Проанализируем приемлемую форму спиральных кривых.

Рис. 2. Инженерный эскиз и объемное изображение спирального профилирования элементов дисковой фрикционно-вихревой турбины

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

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

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

23

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

Рис. 3. Спиральная траектория движения элементов нерасширяющегося рабочего тела в фрикционно-вихревой турбине

Расстояние между двумя последовательными витками спирали, учет которого важен при разработке пространственных геометрических форм рабочего колеса ротора, определяется в плоском сечении формулой:

l = а (ф + 2ж) - ар = 2аж , (20)

где -полярный угол;

а - названный выше коэффициент пропорциональности (величина, являющаяся некоторой константой для конкретной спиральной кривой, параметр спирали).

Исходя из изложенного, может быть предложено инженерное обобщение рассмотренного плоского случая на случай трехмерного пространственного вихреобразования для усовершенствованных бездисковых безлопастных турбинах трения [11] с рационально организованным (спиральным) движением рабочего тела в проточной части.

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

Конструкция предусматривает корпус с круговой цилиндрической камерой, который является статором турбины. В корпусе установлен вал вращения с рабочим колесом, образующие ротор турбины. Геометрическая ось цилиндрической камеры 1 совпадает с геометрической осью вала 2. Входной патрубок 3 для впуска рабочего тела в цилиндрическую камеру расположен на боковой поверхности корпуса 4 и содержит сопло 5. Для выхода рабочего тела служит канал внутри вала 6, на котором крепится рабочее колесо турбины. Рабочее колесо образовано скручено-изогнутым прямоугольным листовым элементом 7 в виде концентрически охватывающих вал витков с контактной спиральной поверхностью трения. Линия крепления рабочего колеса параллельна сквозному проему 8, сделанному в стенке трубчатого вала 6. Фрикционная листовая поверхность (поверхность трения) рабочего колеса содержит параллельные друг другу дистанционирующие ограничительные буртики 9, которые изготовлены так, что имеют профиль в форме плоской

24

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

спирали, размеры и форма которой повторяет конфигурацию профиля спиральной свертки листового элемента ротора.

Рис. 4. Бустерная фрикционо-вихревая турбина тубопривода с дополнительным эффектом использования энергии рабочего тела

Плоскости спиралей буртиков параллельны между собой и перпендикулярны геометрической оси вала, геометрические центры этих спиралей принадлежат этой оси. Буртики являются дистанционирующими элементами (сохраняющими неизменную расстояние между витками, как того требует спираль архимеда) и придают конструкции рабочего колеса жесткость. В полости трубчатого вала, на его внутренних стенках, расположены лопастные элементы 10, которые воспринимают рабочее тело, отработавшее на поверхности трения рабочего колеса, поступившее в прорезь вала.

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

Исследования ориентировано на реализацию идей, предложенных профессором А. В. Королёвым в части использования фрикционно-вихревых турбин в качестве бустерных композиционных элементов турбоприводов.

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ

25

АЛЬТЕРНАТИВНАЯ ЭНЕРГЕТИКА

Выводы

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

Список литературы

1. Деревянко О. В., Королёв А. В., Погосов А. Ю. О выявлении латентных факторов, требующих

превентивного срабатывания специализированной системы подпитки водных сред ЯЭУ //Енергетика та електрифікація, 2014, № 4 (368), С. 44-50.

2. Деревянко О. В., Королёв А. В., Погосов А. Ю. Предаварийные физические процессы и надежный теплоотвод в ядерных энергоустановках. Монография: Одесса, Наука и техника, 2014, - 268 с.

3. Седов Л. И. Плоские задачи гидродинамики и аэродинамики. 3-е изд. М.: Наука,1980. - 488 с.

4. Волков К. Н., Емельянов В. Н. Вычислительные технологии в задачах механики жидкости и газа. -Санкт-Петербург: Лань, 2013. - 468 с.

5. Борисов А. В., Мамаев И. С., Васькина А. В. Новые относительные равновесия в системе трех точечных вихрей в круговой области и их устойчивость //Нелинейная динамика. 2011. Т. 7. № 1. - С. 119-138.

6. Борисов А. В., Мамаев И. С. Математические методы динамики вихревых структур. М.- Ижевск: НИЦ

«РХД», Инст. компьютерн. исслед., 2005. - 368 с.

7. Куракин Л. Г. Об устойчивости томсоновских вихревых конфигураций внутри круговой области // Нелинейная динамика, 2009, т. 5, № 3, С. 295-317.

8. Hampton M., Moeckel R. Finiteness of stationary configurations of the four-vortex problem // Trans.Amer. Math. Soc., 2009, vol. 361, P. 1317-1332.

9. А. В. Борисов, И. С.Мамаев Математические методы динамики вихревых структур. - Москва-Ижевск: НИЦ "Регулярная и хаотическая динамика", 2005, 368 с.

10. Савелов А. А. Плоские кривые: Систематика, свойства, применения. Справочное руководство / Под ред. А. П. Нордена. Изд. 3-е. - М.: Книжный дом «ЛИБРОКОМ», 2010. - 296 с.

11. Патент Украины на полезную модель по заявке № u 2014 05622, решение о выдаче патента от 15.09.2014.

Referenses

1. Деревянко О. В., Королёв А. В., Погосов А. Ю. О выявлении латентных факторов, требующих превентивного срабатывания специализированной системы подпитки водных сред ЯЭУ //Енергетика та електрифікація, 2014, № 4 (368), С. 44-50.

2. Деревянко О. В., Королёв А. В., Погосов А. Ю. Предаварийные физические процессы и надежный теплоотвод в ядерных энергоустановках. Монография: Одесса, Наука и техника, 2014, - 268 с.

3. Седов Л. И. Плоские задачи гидродинамики и аэродинамики. 3-е изд. М.: Наука,1980. - 488 с.

4. Волков К. Н., Емельянов В. Н. Вычислительные технологии в задачах механики жидкости и газа. -Санкт-Петербург: Лань, 2013. - 468 с.

5. Борисов А. В., Мамаев И. С., Васькина А. В. Новые относительные равновесия в системе трех точечных вихрей в круговой области и их устойчивость //Нелинейная динамика. 2011. Т. 7. № 1. - С. 119-138.

6. Борисов А. В., Мамаев И. С. Математические методы динамики вихревых структур. М.- Ижевск: НИЦ <РХД>, Инст. компьютерн. исслед., 2005. - 368 с.

7. Куракин Л. Г. Об устойчивости томсоновских вихревых конфигураций внутри круговой области // Нелинейная динамика, 2009, т. 5, № 3, С. 295-317.

8. Hampton M., Moeckel R. Finiteness of stationary configurations of the four-vortex problem // Trans.Amer. Math. Soc., 2009, vol. 361, P. 1317-1332.

9. А. В. Борисов, И. С.Мамаев Математические методы динамики вихревых структур. - Москва-Ижевск: НИЦ "Регулярная и хаотическая динамика", 2005, 368 с.

10. Савелов А. А. Плоские кривые: Систематика, свойства, применения. Справочное руководство / Под ред. А. П. Нордена. Изд. 3-е. - М.: Книжный дом «ЛИБРОКОМ», 2010. - 296 с.

11. Патент Украины на полезную модель по заявке № u 2014 05622, решение о выдаче патента от 15.09.2014.

26

№1 0 (129) 2014 ЭНЕРГОСБЕРЕЖЕНИЕ • ЭНЕРГЕТИКА • ЭНЕРГОАУДИТ