CC BY
20
УДК 517.9:541.12+544.4
МЕТОДИКА И РЕЗУЛЬТАТЫ ЧИСЛЕННОЙ ПРОВЕРКИ ГИПЕРБОЛИЧЕСКОЙ ПРИРОДЫ АТТРАКТОРОВ ДЛЯ РЕДУЦИРОВАННЫХ МОДЕЛЕЙ РАСПРЕДЕЛЕННЫХ СИСТЕМ
В. П. Круглов
Саратовский государственный университет им. Н.Г. Чернышевского Институт радиотехники и электроники им. В.А. Котельникова РАН, Саратовский филиал
Метод проверки гиперболической природы хаотических аттракторов, основанный на анализе статистики распределения углов между подпространствами устойчивых и неустойчивых направлений, применяется к редуцированным конечномерным моделям нескольких распределенных систем, сконструированных на основе модификации уравнений Свифта-Хохенберга и модели брюсселятора, а также к задаче о параметрическом возбуждении стоячих волн модулированным сигналом накачки.
Ключевые слова: Однородно гиперболический аттрактор, структурная устойчивость, показатель Ляпунова.
Введение. Описание метода
Аттракторы, состоящие только из седловых траекторий, устойчивые и неустойчивые многообразия которых могут пересекаться лишь трансверсально, называют однородно гиперболическими [1-5]. Из трансверсальности многообразий следует свойство структурной устойчивости: устройство их фазового пространства и характер динамики нечувствительны к изменению параметров и функций, входящих в уравнения. Примеры физически реализуемых систем с однородно гиперболическими аттракторами были предложены недавно в работах С.П. Кузнецова и соавторов [6-12], причем некоторые из них относятся к категории распределенных систем и описываются уравнениями в частных производных. Существует необходимость проверки гиперболичности аттракторов этих систем на базе математически обоснованных методик.
Известный критерий конусов [1-5, 7, 13] позволяет подтвердить гиперболичность аттрактора с высокой достоверностью, но его использование требует значительных вычислительных ресурсов и, что наиболее существенно, затруднительным представляется его распространение на случай систем высокой размерности фазового пространстве, в том числе на распределенные системы.
Альтернативный метод состоит в анализе распределений углов между устойчивыми и неустойчивыми многообразиями на аттракторе [2, 14].
Предположим, что необходимо проверить гиперболичность аттрактора некоторого отображения M : xra+1 = M(xn), имеющего размерность фазового пространства к. Если аттрактор гиперболический, то в любой его точке пространство векторов малых возмущений представляется как прямая сумма двух инвариантных подпространств, являющихся касательными к устойчивому и неустойчивому многообразиям [5].
Пусть неустойчивое подпространство одномерным (во всех рассмотренных ниже примерах это так), тогда устойчивое подпространство имеет размерность к — 1.
Неустойчивый вектор в точках произвольной траектории на аттракторе можно получить, выполняя многократные итерации отображения M совместно с линеаризованным отображением DM : Sx„+i = DM(xn)5x„. С другой стороны, для той же траектории в обратном времени можно провести многократные итерации набора из к — 1 линеаризованных отображений DM, подвергая полученные на каждой итерации векторы процедуре ортогонализации и нормировки. Таким образом, для каждой точки траектории находим касательный вектор к неустойчивому многообразию и ортонормированный базис гиперповерхности, касательной к устойчивому многообразию. Это позволяет определить угол между устойчивым и неустойчивым подпространствами. Отсутствие углов, близких к нулю, свидетельствует о гиперболической природе аттрактора. С другой стороны, обнаружение близких к нулю углов говорит о нарушении гиперболичности.
Хотя описанный тест не может претендовать на строгое обоснование гиперболичности аттрактора, он позволяет на практике различать гиперболические и негиперболические аттракторы. В описанной выше форме, однако, использование этого метода вызывает трудности для систем высокой размерности, поскольку требует выполнения решения большого числа уравнений в вариациях. Благодаря недавно открытым алгоритмам вычисления ковариантных ляпуновских векторов [15] появилась возможность модифицировать этот метод [16]. Как и в исходной версии метода, для определения вектора, принадлежащего неустойчивому подпространству, в прямом времени осуществляются итерации линеаризованного отображения DM. Для итераций в обратном времени вдоль той же траектории вводится сопряженное линеаризованное отображение DMT. Для него матрица Якоби, отвечающая линеаризованному отображению DM, транспонируется и берется с обратным знаком. В результате в коде итераций в каждой точке траектории получается вектор, ортогональный гиперповерхности, касательной к устойчивому многообразию.
Опишем алгоритм более детально, имея в виду его применение к аттракторам в сечении Пуанкаре систем обыкновенных дифференциальных уравнений. На первом этапе решаем численно систему уравнений X = F (x, t) совместно с уравнениями в вариациях Sx = F7 (x, t) Sx, нормируя вектор возмущения на каждом шаге, пока не установится режим, соответствующий хаотическому аттрактору, а вектор возмущения будет принадлежать неустойчивому подпространству. На втором этапе, продолжая вычисления, на каждой итерации отображения Пуанкаре записываем в файл значения динамических переменных x = {xi,X2,... ,Xk} и компонент вектора возмущения Sx = {Sxi, Sx2,..., Sxk}. Этот участок траектории должен быть достаточно длинным. На третьем этапе продолжаем решать уравнения и записы-
вать в файл значения динамических переменных на каждой итерации отображения Пуанкаре. На четвертом этапе движемся вдоль записанной в файл опорной траектории в обратном времени и численно решаем сопряженные уравнения в вариациях SU = — [F7 (x, t)]T Su. Когда процедура доходит до участка траектории, где записывался вектор возмущения в прямом времени, начинаем записывать вектор возмущения сопряженной системы на каждой итерации отображения Пуанкаре. В результате для точек на длинном отрезке траектории получаем набор векторов {Sxn}, касательных к неустойчивому многообразию, и набор векторов {Sun}, ортогональных к устойчивому многообразию. В каждой точке можно найти угол |Зп € [0, п/2] между векторами Sxn и Sun по теореме косинусов. Соответственно, угол между устойчивым и неустойчивым подпространствами в этой точке ап = п/2 — |Зп.
Если система дифференциальных уравнений автономная, на каждом шаге итераций из векторов возмущения необходимо исключить составляющую, отвечающую одномерному нейтральному подпространству. Для этого из вектора возмущения вычитается вектор F, направленный вдоль траектории, компоненты которого отвечают функциям в правых частях уравнений, с выбранным надлежащим образом коэффициентом (скалярное произведение вектора возмущения и нормированного на единицу вектора F).
В данной работе описанная методика проверки гиперболичности применяется к редуцированным конечномерным моделям распределенных систем, которые предположительно имеют хаотические аттракторы типа Смейла-Вильямса.
1. Исследуемые системы и их конечномерные модели
1.1. Автономная система на основе уравнения Свифта-Хохенберга. В статье [8] предложен пример автономной распределенной системы на основе модификации уравнения Свифта-Хохенберга [13], описывающей возникновение и затухание в одномерной среде структур, пространственная фаза которых на определенном характерном периоде претерпевает удвоение (с обратным знаком). Модель описывается системой уравнений в частных производных
dtu + (1 + д%)2и = (ц + u2 — v2/5)u + ev cos 3x,
2 2 (1) dtv = (u2 — 1)v + u2,
с периодическими граничными условиями
u(x,t) = u(x + L,t), v(x,t) = v(x + L,t). (2)
Первое уравнение системы (1) представляет собой модифицированное уравнение Свифта-Хохенберга. Функция e cos 3x в последнем слагаемом правой части описывает пространственную неоднородность среды. Второе уравнение - вспомогательное и не содержит производных по координате.
В соответствии с граничными условиями (2) функция, описывающая паттерн в фиксированный момент времени, периодическая по координате x, поэтому переменные u и v можно разложить в ряд Фурье.
Качественно функционирование системы можно описать следующим образом. Когда величина ц + u2 — v2/5, входящая в первое уравнение системы (2), стано-
вится больше нуля, переменная и нарастает, при этом основной вклад в ее величину вносит первая пространственная гармоника с фазой ф: и & \U1\ cos(x + ф). Когда величина и2 — 1, входящая во второе уравнение, становится больше нуля, нарастает переменная v, при этом возбуждаются нулевая и вторая пространственные гармоники, и за счет присутствия слагаемого и2, второй гармонике передается удвоенная фаза: \V2\ cos(2x + 2ф). Произведение второй гармоники переменной v и функции, описывающей пространственную неоднородность среды, обеспечивает передачу первой гармонике переменной и удвоенной фазы с обратным знаком: cos(2x + 2ф) cos 3x = cos(x — 2ф)/2 + ... Далее переменные и и v убывают, и остаются близкими к нулю, пока величина ц + u2 — v2/5 снова не станет положительной. Таким образом, на каждой стадии возбуждения фазы пространственных спектральных компонент претерпевают удвоение с изменением знака: фга+1 = — 2фга. Это -растягивающее отображение окружности, демонстрирующее хаос с показателем Ляпунова, равным ln2 ~ 0.6931.
Система уравнений (1)-(2) решалась численно с помощью явно-неявной разностной схемы. На рис. 1 представлена пространственно-временная зависимость переменных и и v при значениях параметров, указанных в подписи к рисунку. Колебания переменной и в среде представляют собой длительные стадии роста возмущений, чередующиеся с резкими спадами почти до нуля. Средний период между последовательными сбросами возбуждения в данном режиме составил т = 50.37. Можно видеть, что фазовый сдвиг от паттерна к паттерну меняется нерегулярно.
Поскольку функционирование распределенной системы основано на взаимодействии пространственных гармоник, приближенно динамику можно описать конечномерной системой обыкновенных дифференциальных уравнений, в качестве переменных которой выступают несколько наиболее значимых фурье-компонент переменных и и v. Представим решение уравнений (1)-(2) в виде усеченных рядов Фурье
и = Cieix + C*le-ix, v = w + C2e2ix + C%e-2ix
(3)
Здесь Ci - комплексная амплитуда первой гармоники переменной и; C2 - комплексная амплитуда второй гармоники переменной v; w - вещественная амплитуда нулевой гармоники переменной v. Подставляя соотношения (3) в уравнения (1), умножая
Рис. 1. Пространственно-временные реализации переменных и (а) и v (б) при значениях параметров ц = 0.03, е = 0.03, L = 2я, At = 0.001, Ax = L/64 ~ 0.098
первое из получившихся выражений на e гх, а второе на e 2гх или 1, и усредняя их в каждом случае по пространственному периоду 2п, приходим к системе уравнений
1 2 2 1
Ci = (ц - 5w2 + 3|Ci|2 - 5|C2|2)Ci - 5CvClw + 2eC2>,
C2 = (2|Ci|2 - 1)C2 + (w + 1 )C2, (4)
w = -w + 2Re(C2C1*2) + 2|Ci|2(w + 1).
Система уравнений (4) решалась численно с помощью метода Рунге-Кутты четвертого порядка при ц = 0.03, е = 0.03. Для построения отображения Пуанкаре секущая поверхность задавалась условием S = | C 11 - 1 = 0, причем отслеживаются значения переменных в моменты времени, когда пересечение этой поверхности сопровождается уменьшением абсолютной величины переменной C .
На рис. 2, а изображена диаграмма для фазы переменной C1. Как видно из диаграммы, динамика фазы в конечномерной системе описывается растягивающим отображением окружности. На рис. 2, б изображена проекция аттрактора в сечении Пуанкаре на плоскость действительной и мнимой частей переменной C . Предположительно, аттрактор в сечении Пуанкаре представляет собой соленоид Смейла-Вильямса.
Для аттрактора потоковой системы вычислялись показатели Ляпунова по алгоритму Бенеттина с ортогонализацией по Граму-Шмидту. При значениях параметров ц = 0.03, е = 0.03 показатели Ляпунова составили X =
= {0.0123, 0, -0.892, -0.951, -0.972}. Поскольку средний период между последовательными прохождениями сечения Пуанкаре равен т = 52.61, старший показатель Ляпунова для отображения Пуанкаре Лi = Xiт = 0.65, что приблизительно соответствует показателю Ляпунова растягивающего отображения окружности ln2. Остальные показатели Ляпунова отображения Пуанкаре отрицательные (нулевой показатель, полученный для потоковой системы, у отображения отсутствует). Таким образом, элемент объема в фазовом пространстве отображения Пуанкаре испытывает за одну итерацию растяжение по одному направлению и сжатие по трем направлениям. Это соответствует построению гиперболического аттрактора типа Смейла-Вильямса (растяжение происходит по угловой переменной) в четырехмерном фазовом пространстве.
Рис. 2. Итерационная диаграмма для фазы переменной С\ (а) и проекция аттрактора в сечении Пуанкаре на плоскость действительной и мнимой частей переменной С\ (б). Параметры ц = 0.03, е = 0.03
Используя вместо (3) подстановки с большим числом гармоник, можно построить более точные конечномерные модели с гиперболическими аттракторами. В частности, используя вместо (3) соотношения
u = Cieix + C*e-ix,
(5)
v = w + C2e2ix + C*e-2ix + C\e4ix + C* e-4ix,
где C4 - комплексная амплитуда четвертой гармоники переменной v, на основе уравнений (2), получим систему
Ci = (ц -1 w2 + 3|Ci|2 - 21C2I2 - 5 |C4|2)Ci - 2 C2C*,w-
- 2 C*C* C4 + 1 e(C* + C4),
C2 = (2|Ci|2 - 1)C2 + C4C*2 + (w + 1 )C2, (6)
C4 = (2|Ci|2 - 1)C4 + C2C2,
w = -w + 2Re(C2 C*2) + 2|Ci|2(w + 1).
(В подстановке (5) не учитываются гармоники переменной u, следующие за первой, так как их амплитуды оказываются намного меньше, чем амплитуда первой гармоники.)
Система уравнений (6) решалась численно методом Рунге-Кутты четвертого порядка. Сечение Пуанкаре определялось так же, как и для системы (4).
Динамика фазы (рис. 3, а) приближенно описывается растягивающим отображением окружности. Преобразование фазы в системе (6) по сравнению с системой (4) оказывается сглаженным, благодаря учету дополнительной гармоники. На рис. 3, б изображена проекция аттрактора в сечении Пуанкаре на плоскость действительной и мнимой частей переменной C1. Предположительно, аттрактор в сечении Пуанкаре представляет собой соленоид Смейла-Вильямса.
Показатели Ляпунова для аттрактора при значениях параметров ц = 0.03, e = 0.03 составили А = {0.0135, 0, -0.841, -0.934, -0.955, -0.96, -0.9639}. Средний период между последовательными сечениями т = 49.6, поэтому старший показатель Ляпунова отображения Пуанкаре Лl = А^т = 0.67, то есть близок
Рис. 3. Итерационная диаграмма для фазы переменной C (а) и проекция аттрактора в сечении Пуанкаре на плоскость действительной и мнимой частей переменной С\ (б). Параметры ц = 0.03, е = 0.03
к величине ln2, отвечающей равномерно растягивающему отображению окружности. Остальные показатели Ляпунова отрицательны. Это соответствует построению гиперболического аттрактора типа Смейла-Вильямса (если учесть, что растяжение происходит по угловой переменной) в шестимерном фазовом пространстве.
1.2. Система на основе брюсселятора. В работе [9] предложена модель неавтономной распределенной системы на основе системы типа брюсселятор [17, 18] с аттрактором Смейла-Вильямса. Модель описывается системой уравнений в частных производных
dt u = (A — u)(1 + е cos3k0x) — Bu + u2v + y(t)odxxu, dt v = Bu — u2v + y(t)dxxv,
(7)
где
( l, nT <t < nT + T/2,
Y(t) = <
[1/4, nT + T/2 <t < (n + 1)T.
Граничные условия - периодические:
u(x,t) = u(x + L,t), v(x,t) = v(x + L,t), (8)
где L - пространственный период. Фактор (1 + е cos 3kox) описывает пространственную неоднородность затухания в среде. Функция y(t) описывает периодическое переключение коэффициентов диффузии. В соответствии с периодическими граничными условиями, функции u и v периодические по координате и могут быть представлены рядами Фурье.
Система функционирует следующим образом. На первой половине периода модуляции (nT < t < nT + T/2) наиболее значимый вклад вносит первая пространственная гармоника, которая характеризуется некоторой фазой ф: cos (k0x + ф). На второй половине периода модуляции (nT + T/2 <t < (n + 1)T), вторая пространственная гармоника cos (2k0x + 2ф) нарастает, а первая гармоника затухает. На следующем периоде первая гармоника возбуждается под воздействием возмущения, обусловленного произведением второй гармоники u и функции (1 + е cos 3k0x), учитывающей пространственную неоднородность. Этим обеспечивается передача первой гармонике вновь возникающего паттерна удвоенной пространственной фазы (с обратным знаком), что следует из соотношения cos (2k0x + 2ф) cos 3k0x = = 1/2 cos (k0x — 2ф) + .... Этот процесс повторяется на каждом периоде. Таким образом, пространственные фазы на последовательных стадиях формирования волновых структур эволюционируют в соответствии с растягивающим отображением окружности фга+1 = — 2фп + const. Это хаотическое отображение характеризуется положительным показателем Ляпунова Л = ln 2 = 0.693 ...
Уравнения (7)-(8) решались численно при значениях параметров A = 2, B = 4.1, о = 0.25, е = 0.03. Пространственный период был задан L = 2п, а период модуляции T = 32п. При таком выборе параметров волновое число, отвечающее максимальной неустойчивости, составляет k0 = 4A2/о = 2. На рис. 4 показаны
Рис. 4. Пространственно-временные реализации переменных и (а) и v (б) при значениях параметров A = 2, B = 4.1, а = 0.25, е = 0.03, L = 2я, T = 32я, ко = 2
пространственно-временные реализации, полученные в численных расчетах для переменных и и v. Можно видеть, что фазовый сдвиг от паттерна к паттерну меняется нерегулярным образом.
Так как работа системы основана на взаимодействии пространственных гармоник переменных и и v, можно свести ее описание к конечномерной системе обыкновенных дифференциальных уравнений для амплитуд наиболее значительных мод -первых слагаемых разложенного в ряды Фурье возмущения в окрестности неподвижной точки (A, В/А)системы (7)-(8). Минимальный набор укороченных уравнений, который адекватно описывает динамику, можно получить с помощью подстановки
u(x, t) = A + x\ cos 2x + ui sin 2x + x2 cos 4x + u2 sin 4x, v(x, t) = B/A + y1 cos 2x + v1 sin2x + y2 cos4x + v2 sin4x.
Коэффициенты Фурье x1, u1, x2, u2 и y1, v1, y2, v2 - это новые переменные, зависящие только от времени. Подставляя эти разложения в (7), после усреднения по пространственному периоду 2п получим следующую систему уравнений:
x 1 = -x1 — ex2/2 — Bx1 — 4ax1y(t) + F1, u 1 = —u1 + eu2/2 — Bu1 — 4au1y(t) + G1, x 2 = —x2 — 6x1/2 — Bx2 — 16ax2Y(t) + F2, u 2 = —u2 + eu1/2 — Bu,2 — 16au2y(t) + G2,
(9)
y1 = Bx1 — 4y1Y(t) — F1, v1 = Bu1 — 4v1y(t) — G1, y2 = Bx2 — 16y2 Y(t) — F2, v>2 = Bu2 — 16v2Y(t) — G2,
Fi = (2B + mvi/2 + U2V2 + Х2У2 + 3xiyi/4)xi + (A2 + Ж2/2 + u2/2 + u2/4)y1 +
+A(uV + U2V1 + Х1У2 + Х2У1) + BA-1 (U1U2 + X1X2),
G1 = (2B + Х1У1/2 + U2V2 + Х2У2 + 3u1V1/4)u1 + (A2 + x2/2 + u2/2 + x2/4)v1 +
+A(X1V2 + ^У1 - UlУ2 - V1x2) + BA-1 (x1U2 - Ux),
F2 = (2B + U1V1 + x^1 + U2V2/2 + 3X2У2/4)X2 + (A2 + x22 + u1/2 + u2/4)у2+
+A(xlУl - U1V1) + BA_1(x1 - u2)/2,
G2 = (2B + U1V1 + x1y1 + x2y2 / 2 + 3U2V2/4)U2 + (A2 + x2/2 + U1/2 + x2/4)V2 +
+A(x1 V1 + u1 y1) + BA-1x1u1.
(10)
На рис. 5, а можно видеть, что динамика фазы описывается растягивающим отображением окружности. На рис. 5, б изображен аттрактор стробоскопического отображения в проекции на плоскость переменных x1 и U1. Визуально аттрактор можно отнести к аттракторам типа Смейла-Вильямса.
Показатели Ляпунова стробоскопического отображения Пуанкаре составили Л = {0.672, -2.49, -138.14,-148.25, -351.65, 362.86, -1263.72, -1264.64}. Отметим, что старший показатель положителен и близок к ln 2, а остальные - отрицательны, что соответствует аттрактору типа Смейла-Вильямса в восьмимерном фазовом пространстве стробоскопического отображения Пуанкаре.
T = 32я
1.3. Параметрическое возбуждение стоячих волн. В работах [6, 7, 12]
предложена система, в которой осуществляется параметрическое возбуждение длинноволновых и коротковолновых паттернов стоячих волн благодаря соответствующим образом заданной модуляции накачки. Примером может служить нелинейная струна. Уравнение в частных производных имеет вид
р (x) d2u = — (а + u2) dtu — yu + G (t) d%u. (11)
Граничные условия - периодические:
u(x,t) = u(x + L,t). (12)
Слагаемое —(а + u2)dtu обеспечивает нелинейную диссипацию, слагаемое —yu подавляет возмущения с волновым числом к = 0. Присутствие нелинейной диссипации необходимо для стабилизации параметрической неустойчивости на некотором уровне амплитуды колебаний, однако для рассматриваемой здесь системы весьма существенным является еще и то обстоятельство, что кубическая нелинейность обеспечивает возникновение третьей и других нечетных гармоник при колебательноволновых движениях. Функция р (x) = 1+ е sin 4kox описывает линейную плотность струны, при малом значении е распределение массы вдоль струны - слабо неоднородное. Функция G (t) = 1 + acos2(2nt/T) sin2mot + bsin2(2nt/T) sin6mot описывает натяжение струны, величины а и b неотрицательные, a + b < 1. Частота mo и волновое число ко равны друг другу. Длину L примем равной целому числу длин волн с волновым числом ко: L = 2nN/ко.
Предположим, что при накачке на частоте 2mo в системе возбуждена стоячая волна на частоте mo с волновым числом ко, у которой расположение узлов и пучностей характеризуется фазовой постоянной 0, так что в грубом приближении u ~ cos m0t sin (h0x + 0).
Амплитуда волны стабилизируется на конечном уровне благодаря нелинейной диссипации, и из-за ее присутствия колебательно-волновое движение будет иметь также составляющую на третьей гармонике. Пространственная фаза гармоники на частоте 3mo будет, очевидно, даваться величиной 30.
На следующем этапе временной эволюции амплитуда накачки на частоте 2mo становится малой, и колебания на частоте mo затухают. Амплитуда накачки на частоте 6mo, с другой стороны, с какого-то момента оказывается достаточной для развития параметрической неустойчивости стоячей волны с частотой 3mo и волновым числом 3ко. Эта волна формируется из начального возмущения, отвечающего третьей гармонике волны, существовавшей на предыдущей стадии, так что пространственная фаза будет определяться величиной 30.
Далее вновь наступает стадия интенсивной накачки на частоте 2mo. Теперь затравочное колебательно-волновое возбуждение с волновым числом ко обеспечивается комбинацией возмущения u ~ cos3mot sin (3^x + 30), оставшегося от предыдущей стадии. Поскольку эта комбинация содержит слагаемое, пропорциональное cos mot sin (^x — 30), новая фаза получается из предыдущей посредством растягивающего отображения окружности: 0ra+i = —30га+const. Это хаотическое отображение характеризуется положительным показателем Ляпунова Л = ln3 ~ 1.098. При
наличии сжатия в фазовом пространстве по остальным направлениям, для отображения, описывающего изменение состояния системы за период модуляции накачки, на этом основании можно предположить присутствие аттрактора типа соленоида Смейла-Вильямса.
В статье [1] рассмотрены также варианты модели, где для фазы имеет место отображение окружности с растяжением в 5, 7, 9, 11 раз.
Система уравнений (11)-(12) решалась численно с использованием явной схемы типа «крест». На рис. 6 представлена пространственно-временная диаграмма, визуализирующая эволюцию огибающей для паттернов стоячих волн; форма струны показана в моменты, кратные периоду быстрых колебаний. Параметры приведены в подписи к рисунку. Как можно видеть, в системе поочередно возбуждаются длинноволновые и коротковолновые структуры, притом пространственная фаза их меняется от одного периода модуляции накачки к другому, претерпевая растягивающее отображение окружности. Фаза колебаний стоячих волн во времени, с другой стороны, жестко привязана к фазе накачки.
Поскольку параметрическое возбуждение в рассмотренной системе имеет место на волновых числах ко и 3ko, для сведения к конечномерной системе обыкновенных дифференциальных уравнений, естественно искать решение уравнений (11)-(12) в виде
Рис. 6. Пространственно-временные реализации переменной и при значениях параметров ш0 = 2п, к0 = 2п, L =1, a = 0.4, b = 0.2, е = 0.2, а = 0.4, Y = 0.03, T = 60
u(x, t) = y1(t) cosk0x + v1(t) sink0x + y3(t) cos 3k0x + v3(t) sin3k0x. (13)
Подстановка этого выражения в (11)-(12) и разложение обеих частей уравнения в ряд Фурье с учетом гармоник того же порядка приводят к системе уравнений для коэффициентов
.. Fi - еСз/2 ..
У1 = 1 - £2/4 ’ v1
Gi - sF:i/2 ..
1 - е2/4 ’ У3
F3 - eGi/2 1 - е2 / 4 ’
V3
G3 - eFi/2 1 - е2/4 ’
(14)
где
F1 = -ш0 (1 + a cos2 T sin2m0t + b sin2 П sin6rn0t)y1 - yy1 - ay1-
-2 [(3y2/2 + y1y3 + у3 + v2/2 + v1v3 + vDy1 + 2 (y2 + 4y1y3 - v2)y3+
+V1V1 - y3v1 + У1 v^v 1 + (y1 v1 + 2yv)v3],
( 2 nt 1 2 nt \
G1 = -шо{ 1 + a cos t sin2wot + b sin t sin6wot)v1 - yv1 - a v 1-
-2 [(3v2/2 - vrv3 + v§ + y1/2 - У1У3 + yl)v 1 + 2 (y2 + 4vrv3 - v1)v3+
+ (yv - V3v1 + ym)У1 + (-yv + 2y3v1)y3],
Рис. 7. Итерационная диаграмма для фазы ф = arctan(uivi/yi) (а) и проекция аттрактора в сечении Пуанкаре на плоскость переменных х1 и и1 (б). Параметры ю0 = 2я, a = 0.4, b = 0.2, е = 0.2,
а = 0.4, у = 0.03, T = 60
F3
-9юо(1 + a cos2 T sin 2mot + b sin2 T sin6rnot)y3 - ууз - ауз-
1 [(y2 + 4yiy3 - vjOyi + 2(-yiVi + 2уз^1 )vi+
+(2y2 + 2v2 + 3y2 + vl)y3 + 2y3V3V3],
G3 = -9ю0(1 + a cos2 T sin2m0t + b sin2 T sin6m0t)v3 - yv3 - аV3-
-4 [2(y i v i + 2yi v3)y/i + (y2 + 4v iv3 - v2)v i + 2y3v3y3+
+(2y 2 + 2v2 + yf + 3v2)v)3].
На рис. 7, б изображен аттрактор стробоскопического отображения системы (14) в проекции на плоскость yi и vi. Визуально его можно отнести к аттракторам типа Смейла-Вильямса. На рис. 7, а приведена итерационная диаграмма для фаз ф = arctan(v i/y i). Можно видеть, что динамика фазы описывается растягивающим отображением окружности.
С помощью алгоритма Бенеттина был вычислен спектр показателей Ляпунова аттрактора стробоскопического отображения Пуанкаре: Л = {1.099, -22.1, -30.44, -40.78, -43.61, -46.43, -56.67, -69.97}. Отметим, что старший показатель положителен и близок к величине ln3, а остальные отрицательные. Это соответствует аттрактору типа Смейла-Вильямса в восьмимерном фазовом пространстве стробоскопического отображения Пуанкаре.
2. Результаты численной проверки гиперболической природы аттракторов
Для всех конечномерных моделей систем, описанных в предыдущем разделе, были получены гистограммы распределений углов между устойчивыми и неустойчивыми многообразиями аттракторов в сечении Пуанкаре.
На рис. 8 представлены гистограммы распределений углов для отображений Пуанкаре пятимерной и семимерной модели модифицированного уравнения Свифта-Хохенберга. Оба распределения дистанцированы от нуля, что подтверждает гиперболическую природу аттракторов. Распределение для пятимерной аппроксимации
более широкое и состоит из нескольких «частей». Более узкое распределение для семимерной аппроксимации означает, что аттрактор устроен примерно одинаково во всех областях. Семимерная модель лучше отражает свойства исходной распределенной системы, чем пятимерная.
На рис. 9 приводится гистограмма распределения углов между устойчивыми и неустойчивыми направлениями для аттрактора конечномерной модели модифицированной системы на основе брюсселятора. Нулевые углы в распределении отсутствуют, что подтверждает предположение о гиперболичности аттрактора. Отметим, что распределение сосредоточено в очень узком диапазоне углов.
На рис. 10 представлена гистограмма распределения углов для конечномерной модели струны с параметрическим возбуждением. Распределение дистанцировано от нуля, что подтверждает гиперболическую природу аттрактора системы.
Заключение
В данной работе к трем редуцированным конечномерным моделям распределенных систем, предложенных как примеры, характеризующиеся присутствием хаотических аттракторов типа Смейла-Вильямса, была применена предложенная в [12] модифицированная методика проверки гиперболичности. Тест состоит в построении распределения углов между подпространствами устойчивых и неустойчивых направлений на хаотическом аттракторе и позволяет отличить режимы с гиперболическим и негиперболическим аттракторами. Благодаря использованию нового алгоритма [14], метод требует относительно небольшого объема вычислений, поскольку и устойчивое и неустойчивое подпространство задаются лишь одним
Рис. 8. Гистограмма распределения углов между устойчивым и неустойчивым направлениями на аттракторе в сечении Пуанкаре пятимерной (а) и семимерной (б) модели модифицированного уравнения Свифта-Хохенберга. Параметры ц = 0.03, е = 0.03
2500л---------------------------------------
о-и------------------------------------
0 а к/2
Рис. 9. Гистограмма распределения углов между устойчивым и неустойчивым многообразиями аттрактора в сечении Пуанкаре конечномерной модели модифицированной системы типа брюсселятор. Параметры A = 2, B = 4.1, а = 0.25, е = 0.03, T = 32я
Рис. 10. Гистограмма распределения углов между устойчивым и неустойчивым направлениями на аттракторе в сечении Пуанкаре конечномерной модели струны с параметрическим возбуждением при ш0 = 2я, a = 0.4, b = 0.2, е = 0.2, а = 0.4, Y = 0.03, T = 60
вектором каждое. Это делает удобным применение теста к системам произвольной размерности. В работе методика применена только к системам с одномерным неустойчивым подпространством, однако ее можно применять и к системам с произвольно большими размерностями устойчивых и неустойчивых подпространств [12]. Для всех рассмотренных примеров - системы на основе модификации уравнений Свифта-Хохенберга и модели брюсселятора, а также системы с параметрическим возбуждением стоячих волн - тест подтверждает предположение об однородно гиперболической природе аттракторов в отображении Пуанкаре.
Автор благодарен профессору С.П. Кузнецову за обсуждение и помощь.
Работа выполнена при поддержке гранта РФФИ№14-02-31162 и гранта поддержки ведущих научных школ НШ-1726.2014.2.
Библиографический список
1. Кузнецов С.П. Динамический хаос и гиперболические аттракторы: от математики к физике. Москва; Ижевск: Институт компьютерных исследований, 2013. 488 с.
2. Кузнецов С.П.Динамический хаос и однородно гиперболические аттракторы: от математики к физике // Успехи физических наук. 2011. T.181, № 2. С. 121.
3. Синай Я.Г. Стохастичность динамических систем // В кн. Нелинейные волны. М.: Наука, 1979. С. 192.
4. Anishchenko V.S., Kopeikin A.S., KurthsJ., Vadivasova T.E., Strelkova G.I. Studying hyperbolicity in chaotic systems // Physics Letters A. 2000. Vol. 270. P. 301.
5. KatokA., Hasselblatt B. Introduction to the Modern Theory of Dynamical Systems. Cambridge: Cambridge University Press, 1996.
6. Круглов В.П., Кузнецов А.С., Кузнецов С.П.Гиперболический хаос в системах с параметрическим возбуждением паттернов стоячих волн // Нелинейная динамика. Т.10, №3. 2014. c. 265-277.
7. Isaeva O.B., Kuznetsov A.S., Kuznetsov S.P. Hyperbolic chaos of standing wave patterns generated parametrically by a modulated pump source // Phys. Rev. E. 2013. Vol.87. 040901.
8. Kruglov V.P., Kuznetsov S.P., Pikovsky A. Attractor of Smale-Williams type in an autonomous distributed system // Regular and Chaotic Dynamics. 2014. Vol. 19, № 4. P. 483.
9. Kruglov V.P. Attractor of Smale-Williams type in modified Brusselator model // Book of Abstracts. International Conference <Nonlinear Dynamics of Deterministic and Stochastic Systems: Unraveling Сomplexity> 19-23 May 2014, Saratov: Saratov State University. 2014. P. 26.
10. Kuptsov P. V, Kuznetsov S.P., Pikovsky A. Hyperbolic chaos of Turing patterns// Phys. Rev. Lett. 2012. Vol. 108. 194101.
11. Kuznetsov S.P. Example of a physical system with a hyperbolic attractor of the Smale-Williams type // Phys. Rev. Lett. 2005. Vol. 95. 144101.
12. Kuznetsov S.P. Some mechanical systems manifesting robust chaos // Nonlinear Dynamics and Mobile Robotics. 2013. Vol. 1, № 1. P. 3.
13. Cross M.C. and Hohenberg P.C. Pattern formation outside of equilibrium // Rev. Mod. Phys. 1993. Vol. 65, № 1. 1993. P. 851.
14. Lai Y.-C., Grebogi C., Yorke J.A., Kan I. How often are chaotic saddles nonhyperbolic? // Nonlinearity. 1993. Vol. 6. P. 779.
15. Kuptsov P.V., Parlitz U. Theory and computation of covariant Lyapunov vectors // Journal of nonlinear science. 2012. Vol. 22, № 5. P. 727.
16. Kuptsov P V Fast numerical test of hyperbolic chaos // Phys. Rev. E. 2012. Vol. 85, № 1. 015203.
17. Гленсдорф П., Пригожин И. Термодинамическая теория структуры, устойчивости и флуктуаций. М.: Мир, 1973. 280 с.
18. Kuznetsov S.P., Mosekilde E., Dewel G., Borckmans P Absolute and convective instabilities in a one-dimensional Brusselator flow model // The Journal of chemical physics. 1997. Vol. 106, № 18. P. 7609.
Поступила в редакцию 25.12.2014
TECHNIQUE AND RESULTS OF NUMERICAL TEST FOR HYPERBOLIC NATURE OF ATTRACTORS FOR REDUCED MODELS OF DISTRIBUTED SYSTEMS
V. P. Kruglov
Saratov State University
V.A. Kotelnikov Institute of Radio Engineering and Electronics of RAS, Saratov Branch
A test of hyperbolic nature of chaotic attractors, based on an analysis of statistics distribution of angles between stable and unstable subspaces, is applied to reduced finitedimensional models of distributed systems which are the modifications of the Swift-Hohenberg equation and Brusselator model, as well as to the problem of parametric excitation of standing waves by the modulated pump.
Keywords: Uniformly hyperbolic attractor, structural stability, Lyapunov exponents.
Круглов Вячеслав Павлович - родился в 1990 году. В 2007 году окончил среднюю школу № 8 г. Саратова. В 2012 окончил факультет нелинейных процессов Саратовского государственного университета им. Н.Г. Чернышевского.В настоящее время является аспирантом ФНП СГУ и работает в Саратовском филиале Института радиотехники и электроники им. В.А. Котельникова РАН, в должности стажера-исследователя. Участвовал в научных конференциях «ХАОС», «Нелинейные волны», «Нелинейные дни для молодых», «Наноэлектроника, нанофотоника и нелинейная физика», «Гамильтонова динамика, неавтономные системы и структуры в уравнениях с частными производными», «Динамика, бифуркации и странные аттракторы». Автор 5 статей в российских и международных научных журналах.
410012 Саратов, ул. Астраханская, 83
Саратовский государственный университет им. Н.Г. Чернышевского 410019, Саратов, ул. Зеленая, 38
Институт радиотехники и электроники им. В.А. Котельникова РАН, Саратовский филиал E-mail: [email protected]
