CC BY
11Нелинейная динамика. 2016. Т. 12. № 2. С. 197-209. Полнотекстовая версия в свободном доступе http://nd.ics.org.ru Б01: 10.20537/па1602003
ОРИГИНАЛЬНЫЕ СТАТЬИ
УДК: 53.043 М8С 2010: 34Е05
Химерные режимы в кольце элементов с локальным однонаправленным нелинейным взаимодействием
И. А. Шепелев, Т. Е. Вадивасова
Сложные пространственные структуры, состоящие из стационарных областей с когерентным и некогерентным поведением соседних элементов, названные химерами, вызывают в последнее время большой интерес исследователей. С существованием подобных структур связан ряд еще не решенных проблем. Одна из таких проблем касается характера взаимодействия элементов ансамбля, при котором возможно наблюдать устойчивые химерные структуры. До последнего времени считалось, что важнейшим условием существования химеры является нелокальный характер взаимодействия. Однако возможно, что это предположение не вполне справедливо. При особом выборе локальных связей химеры могут быть реализованы (например, химера в ансамблях с инерционной локальной связью). Для описания связи вводится дополнительная переменная, задаваемая линейным дифференциальным уравнением. Существование виртуальных химер в осцилляторах с запаздывающей обратной связью также позволяет предположить возможность получения химер в кольце локально связанных осцилляторов с однонаправленным взаимодействием, которое является безынерционным, но носит нелинейный характер. Такое предположение основано на качественной аналогии, которая прослеживается в поведении системы с запаздывающей обратной связью и кольца соответствующих осцилляторов с локальной однонаправленной связью.
В данной работе за основу принимается система с запаздывающей обратной связью, в которой существует виртуальная химера, и построен распределенный аналог, представляющий собой кольцо осцилляторов с однонаправленной нелинейной локальной связью.
Ключевые слова: осциллятор с запаздывающей обратной связью, распределенная система, пространственная структура, химера, динамический хаос
Получено 03 марта 2016 года После доработки 21 апреля 2016 года
Шепелев И. А. благодарит за финансовую поддержку РНФ, грант 16-12-1075, в рамках технического задания которого выполнены разделы 1-3 настоящей статьи. Вадивасова Т. Е. благодарит за финансовую поддержку РФФИ, грант №14-52-12002.
Шепелев Игорь Александрович igor_sar@li.ru
Вадивасова Татьяна Евгеньевна vadivasovate@yandex .ги
Саратовский государственный университет им. Н.Г.Чернышевского 410012, Россия, г. Саратов, ул. Астраханская, д. 83
Введение
В последние годы большой интерес исследователей вызывают сложные пространственные структуры, получившие название химер [2, 13, 23]. Химеры возникают в осцилляторных ансамблях с различной динамикой элементов, как регулярной [3, 22, 29, 31], так и хаотической [19-21], и представляют собой чередующиеся в пространстве кластеры с когерентным и некогерентным поведением соседних осцилляторов. Одним из важнейших условий возникновения химер принято считать нелокальный характер взаимодействия элементов: каждый осциллятор непосредственно связан с целой группой соседей. В большинстве работ, посвященных химерам, исследуются модели ансамблей, учитывающие связь каждого элемента с несколькими соседями, то есть нелокальный характер связи. Имеются также результаты, свидетельствующие о существовании химер в ансамблях с глобальным взаимодействием [28].
Вопрос о возможности существования химерных структур в ансамблях с локальной связью остается открытым. В большинстве случаев при уменьшении радиуса взаимодействия химеры исчезают. Однако можно предположить, что при некоторых видах локального взаимодействия химероподобные структуры все-таки могут существовать. Так, в [14] приведены примеры химер в ансамблях с локальной инерционной связью. Связь задается специальной переменной, которая описывается линейным дифференциальным уравнением. Возможность метастабильных амплитудных химер в ансамбле осцилляторов Стюарта-Ландау с локальной связью отмечалась в [31]. Однако случай локальной связи не был специально исследован.
Особый случай химер имеет место в системах с запаздывающей обратной связью [15, 16]. Это так называемые виртуальные химеры. Они представляют собой режим перемежаемости во времени, когда на периоде запаздывания укладывается несколько интервалов с регулярным и нерегулярным поведением, причем чередование интервалов происходит упорядо-ченно, практически повторяясь через период запаздывания. Период запаздывания в таких системах может рассматриваться как условное пространство [1, 15, 16].
Системы с запаздыванием составляют особый класс распределенных динамических систем, для которых мгновенное состояние в момент времени Ь задается набором функций, определенных на интервале [Ь; ¿+ТУ, где Т^ — период запаздывания. Модели с запаздыванием широко используются в механике [18], задачах управления [5, 24-26], экологии [7, 9, 17], нейродинамике [6, 8, 27] и ряде других областей. Осциллятор с запаздывающей обратной связью при некоторых условиях демонстрирует поведение, подобное поведению распределенной системы типа кольца осцилляторов с локальным однонаправленным взаимодействием и может в определенной степени рассматриваться как упрощенная модель такой распределенной системы [4, 10, 11].
Существование виртуальных химер в системах с запаздыванием и отмеченная аналогия позволяют предположить, что химеры могут возникать в кольце локально связанных осцилляторов, по крайней мере в случае однонаправленной связи. Если это утверждение верно, то по крайней мере определенный класс химерных структур может быть реализован в ансамблях с локальным взаимодействием элементов. В данной работе приводятся результаты численного моделирования, подтверждающие сделанное предположение. Рассматривается ансамбль осцилляторов типа системы Икеды без запаздывания, соединенных в кольцо посредством нелинейных однонаправленных локальных связей. Показано, что данная распределенная система ведет себя во многом аналогично модели с запаздыванием, исследованной в [15, 16], и демонстрирует химеры, имеющие большое сходство с описанными
в [15, 16] виртуальными химерами. В работе также рассматривается эволюция химер при вариации различных параметров исследуемой системы.
Статья организована следующим образом: в разделе 1 описывается исследуемая система и методы численного моделирования, в разделе 2 описывается эволюция системы, показывается существование режимов химер и производится их анализ, в разделе 3 рассматривается явление мультистабильности, в разделе 4 рассматриваются диаграммы режимов исследуемой системы на плоскостях различных управляющих параметров. В заключении на основании полученных данных делаются определенные выводы.
1. Исследуемая система и методы численных экспериментов
Виртуальная химера, численно исследованная в [15] и реализованная экспериментально в устройстве с оптоволоконной запаздывающей обратной связью [16], может быть исследована с помощью простой математической модели, представляющий собой уравнение линейного осциллятора с запаздывающей нелинейной обратной связью — модифицированный осциллятор Икеды:
= ~5°У ~х + kfix(T ~ Т<г)Ь , dt (1.1)
dt '
где т — время в исходных единицах, T¿ — время запаздывания, ео и §о — параметры осциллятора, k — коэффициент обратной связи, f (x) — нелинейная функция, характеризующая преобразование в цепи с запаздывающей обратной связью. Если нормировать время на
период запаздывания t = то можно переписать уравнения системы в виде
Td
(1.2)
с!у _ М ~ Х'
где е = 5 = Т^£о, у = 7^-. Функция /(ж) может быть задана в виде Т d Т d
№ = 7--.1 (1-3)
1 + т 8ш2(ж + Фо)
в, т, Фо — параметры нелинейной обратной связи.
Модель, исследуемая в настоящей работе, получается заменой цепочки запаздывающей обратной связи на цепочку идентичных осцилляторов, замкнутую в кольцо. Связь элементов является локальной, однонаправленной и задается нелинейной функцией (1.3). Соответствующие уравнения имеют вид
= ~5У3 ~ Х3 +
^ (1.4)
сМ ~Хз'
1 + m sin2(xj-i + Ф0)
—2>
Здесь ] = — номер осциллятора (дискретная пространственная координата), N —
число осцилляторов в кольце. В проведенных исследованиях рассматривалось кольцо из N = 300 элементов. Параметры связи фиксировались постоянными: в = 4, т = 4.7,
Модель (1.4) может быть физически реализована, например, в виде кольца из полупроводниковых лазеров с интерферометрами Фабри-Перо, характеризующиеся нелинейностью вида (1.3), связанных оптоволоконными линиями. Длина этих линий предполагается достаточно малой, чтобы запаздыванием сигнала в линии можно было пренебречь.
Система (1.4) численно интегрировалась методом Гюна [30]. Шаг интегрирования фиксировался равным Н = 0.0004. По полученным массивам данных строились пространственно-временные графики режимов, на которых разные цвета соответствуют различным уровням значений переменной х. Кроме того, для диагностики режимов анализировались колебания во времени в фиксированной точке пространства и мгновенные пространственные профили. Начальные условия выбирались случайным образом или задавались периодическими с добавлением случайной компоненты. Реализация какого-либо режима со случайных начальных условий позволяет сделать вывод о его устойчивом характере. При построении диаграмм режимов на плоскости управляющих параметров отдельно отслеживался каждый устойчивый режим. Важно отметить, что для исследуемой системы характерен очень длительный процесс установления (>10 000), поэтому при построении пространственно-временных графиков использовались данные, полученные после интегрирования на достаточно больших интервалах времени.
2. Эволюция режимов с ростом параметра связи и возникновение химероподобной структуры
При вариации параметров е, 5 и к в кольце (1.4) наблюдаются различные пространственно-временные структуры. Зафиксируем значения е = 0.01, 5 = 0.01 и рассмотрим, как меняются динамические режимы в кольце (1.4) с ростом связи к.
При очень слабой связи в системе отсутствуют какие-либо колебания и она находится в устойчивом однородном равновесном состоянии х = 0, у = 0. С ростом к появляется колебательный режим, представляющий собой перемежаемость во времени хаотического поведения и режима, близкого к однородному равновесию (рис. 1). При режиме, близком к равновесию, некоторые пространственные колебания имеют место, но их амплитуда очень мала (рис. 1Ь). В области хаотической динамики, напротив, в пространстве наблюдаются хаотические колебания с достаточно большой амплитудой (рис. 1с). Аналогичное поведение имело место в системе (1.2) и было названо бризером [15]. Колебания во времени в указанном режиме представляют собой периодические переключения между двумя метастабиль-ными состояниями, в одном из которых наблюдается временной хаос (рис. 1ё).
При достижении коэффициентом связи некоторого значения в системе (1.4) происходит резкий переход к режиму устойчивой бегущей волны (рис. 2ё). Особенность режима бегущей волны, представленного на рисунке 2, состоит в том, что часть элементов кольца находится в одном и том же метастабильном состоянии покоя (темные области на пространственно-временном графике), а часть — в состоянии пространственно-временного хаоса (светлые полосы). Данный режим соответствует виртуальной химере, наблюдаемой в системе с запаздыванием (1.2) [15, 16]. Применительно к пространственно распределенной системе (1.4) данный режим может быть назван «вращающейся химерой».
Фо = 0.4.
100 150 200 250 300
3
300
0 50 100 150 200 250 300
г
Рис. 1. Режим перемежаемости во времени: (а) пространственно-временной график колебаний; (Ь) пространственный профиль в режиме, близком к однородному равновесию; (е) пространственный профиль в хаотической фазе временной перемежаемости; (^ реализация колебаний во времени."
Границы областей когерентного и некогерентного поведения вращаются по кольцу с некоторой постоянной скоростью. Стационарную картину областей можно получить, перейдя в движущуюся систему координат. На рисунках 2а, Ь приведены пространственно-временные графики в режиме вращающейся химеры для неподвижных и движущихся координат. Из-за конечного числа элементов кольца не удается точно задать скорость вращения химеры и полностью избавится от дрейфа границ областей: границы областей на рисунке 2 не являются строго вертикальными. Однако примененная процедура позволяет наблюдать химерную структуру на пространственно-временном графике для достаточно большого интервала времени. Как видно из рисунка 2Ь, пространственная структура почти не меняется во времени, то есть можно говорить о стабильном во времени химерном состоянии. Форма соответствующих колебаний в пространстве и во времени приведена на рисунках 2с, ^
Важной чертой отмеченного режима, указывающей на то, что этот режим действительно можно называть химерой, является независимость от начальных условий. На рисунке 3 приведены два различных начальных распределения переменных х^ и соответствующие им установившиеся пространственно-временные структуры. Можно видеть, что с разных начальных условий система выходит на один и тот же установившийся режим «вращающейся химеры». Можно отметить, что для разных начальных распределений время установления различается: так, при периодическом распределении оно составляет примерно 6 000-8 000 единиц, а для случайного начального профиля — уже порядка 30 000.
аДля читателя печатной версии: здесь и далее полноцветные версии рисунков см. в эл. версии статьи — http://nd.ies.org.ru/nd1602003/
3 t
Рис. 2. Режим «вращающейся химеры» в исследуемой системе при г = 0.01, 5 = 0.41, k = 0.456: (a) пространственно-временной график в неподвижной системе координат; (b) пространственно-временной график в движущейся системе координат; (с) мгновенный пространственный профиль; (d) форма колебаний во времени в фиксированной точке пространства. Время установления tset = = 40 000, начальные условия — случайные.
Рис. 3. Химерные состояния для разных начальных условий при г = 0.01, 5 = 0.01, к = 0.51: начальные распределения х^ (а, Ь) и соответствующие им пространственно-временные графики установившегося режима «вращающейся химеры» (е,^.
Химероподобные структуры существуют в некоторой конечной области значений параметра связи. При дальнейшем увеличении связи между элементами химеры разрушаются (рис. 4). Разрушение химер происходит постепенно.
В области некогерентности возникают кластеры с когерентным поведением, которые существуют на конечных интервалах времени (рис. 4а, Ь). С ростом к таких кластеров появляется все больше, а время их жизни уменьшается (рис. 4с). При некотором достаточно большом значении к в системе вновь устанавливается пространственно-временной хаос (рис. 4ё).
Следует отметить, что механизм разрушения виртуальной химеры в осцилляторе с запаздывающей обратной связью при увеличении коэффициента обратной связи является совершенно аналогичным [15]. Он также наблюдается в экспериментах с лазером и электронной запаздывающей обратной связью [16].
О 50 100 150 200 250 300
3
50 100 150 200 250 300
3
0 50 100 150 200 250 300
3
Рис. 4. Разрушение химеры при е = 0.01, 5 = 0.01 с ростом значений параметра связи: (а) к = 0.519; (Ь) к = 0.5258; (е) к = 0.56; к = 0.8.
3. Мультистабильность
Независимость режима вращающейся химеры от начальных условий, проиллюстрированная на рисунке 3, не является абсолютной. Во многих случаях один и тот же режим наблюдается при изменении начальных условий только в определенных пределах. В то же время существует множество начальных условий, которые приводят к установлению другого режима (например, химеры с другим числом голов (кластеров некогерентности) на длине системы), то есть имеет место явление мультистабильности.
Пример двух химерных структур, сосуществующих при одних и тех же значениях параметров, приведен на рисунке 5.
Рис. 5. Два режима химер, наблюдаемые при е = 0.01, 5 = 0.41, к = 0.456 в зависимости от начальных условий.
4. Диаграммы режимов
на плоскостях управляющих параметров
Были исследованы области существования различных пространственно-временных структур при изменении параметров системы (1.4) и построены диаграммы режимов на плоскостях различных управляющих параметров. Они приведены на рисунках 6-8. Различным цветам на диаграммах соответствуют различные режимы. Там же приведены пространственно-временные графики, иллюстрирующие характер режима, наблюдаемого в той или иной области. Границы областей определены достаточно приближенно, так как их детальное построение в силу длительности процессов установления, мультистабильности и сложности режимов представляется затруднительным. Тем не менее, приведенные диаграммы дают адекватное представление об эволюции режимов при изменении того или иного параметра.
На рисунке 6 приведены области различных динамических режимов на плоскости параметров (к,е) при 5 = 0.01. В светлой области (область 1) наблюдается колебательный режим, представляющий собой перемежаемость во времени между хаотическим режимом
0.485 0.49 0.495 0.5 0.505 0.51 0.515 0.52 к
Рис. 6. Диаграмма режимов в системе (1.4) на плоскости управляющих параметров (к,е) при 5 = 0.01. Цветами выделены области разных динамических режимов. Характер режимов проиллюстрирован вставками, представляющими пространственно-временные диаграммы для переменной х(£).
и режимом, близким к однородному равновесию (см. рис. 1). Две области более темных тонов (2.1 и 2.2, синяя и голубая) соответствуют существованию химеры с одним кластером некогерентности. Меньшая по размеру и более светлая область 2.2 соответствует случаю максимально устойчивой химеры, когда она реализуется при случайных начальных условиях. В более широкой темной области 2.1 для получения того же режима химеры требуются особые начальные условия. При переходе в светлую область справа (область 3) химера разрушается. Граница между областями является достаточно условной, поскольку разрушение химеры с ростом параметра к происходит постепенно (см. рис. 4).
Для значения 5 = 0.01 явление мультистабильности не наблюдается. Однако оно возникает с ростом 5, что можно видеть на диаграмме режимов, приведенной на рисунке 7. Управляющими параметрами теперь являются к и 5, а значение е постоянно: е = 0.01. Здесь в областях 2.1-2.4, обозначенных градациями тона (от светло- до темно-синего) существуют режимы химер, причем с переходом в более темные области (слева направо 2.1 ^ 2.4) возникают новые химеры с большим числом некогерентных кластеров и сохраняются уже существующие химеры.
Таким образом, в системе развивается мультистабильность: в самой темной из отмеченных областей 2.4 сосуществуют химеры с числом кластеров некогерентности от одного до четырех. Далее, с переходом в светлую область 3, происходит разрушение химер. Разрушение химерных структур сопровождается длительными переходными процессами, причем длительность этих процессов различна для химер с различным числом кластеров.
1 2.1 2.2 2.3 2.4 3
0.35 0.3 0.25 0.2 0.15 0.1 0.05
0 0.3 0.35 0.4 0.45 0.5 к
Рис. 7. Диаграмма режимов в системе (1.4) на плоскости управляющих параметров (к, 5) при е = 0.01.
Меняя параметры е и 5 при постоянном коэффициенте связи (было выбрано значение к = 0.505), также можно наблюдать возникновение вращающихся химер и развитие мультистабильности (рис. 8). Если при малых значениях е и 5 (область 1) имеет место колебательный режим, представляющий собой перемежаемость во времени между хаотическим и близким к однородному режимами (см. рис. 1), то при увеличении сразу обоих параметров происходит возникновение всё новых химер (перекрывающиеся области 2.1-2.5).
Сравнивая диаграмму режимов на рисунке 8 с аналогичной диаграммой, полученной экспериментально для лазера с запаздывающей обратной связью [16], можно увидеть их несоответствие: имеет место как бы зеркальное отражение диаграмм. В модели (1.4) области мультистабильности химерных режимов наблюдаются при больших значениях е, а на диа-
1 2.1 2.2 2.3 2.4
0
0 0.005 0.01 0.015 0.02 0.025
Рис. 8. Диаграмма режимов в системе (1.4) на плоскости управляющих параметров (к, 5) при е = 0.01.
грамме из [16], напротив, при малых значениях е. Это несоответствие легко объяснить. Параметры рассматриваемых систем неидентичны и могут быть связаны с нелинейным преобразованием. Так, в модели (1.4) аналогом времени запаздывания Та является период Тг вращения волны вокруг кольца осцилляторов. Он зависит от параметров системы.
Расчеты показали, что с ростом е период Та в системе (1.4) увеличивается, причем нелинейно. Таким образом, большим значениям е в системе (1.4) соответствуют малые значения параметра е в системе с запаздывающей обратной связью (1.2).
Заключение
Проведенные исследования показали качественную аналогию поведения линейного осциллятора с нелинейной запаздывающей обратной связью и кольца из соответствующих осцилляторов с локальным однонаправленным нелинейным взаимодействием. В кольце были обнаружены химероподобные режимы, представляющие собой чередование пространственных кластеров, соответствующих почти равновесному состоянию элементов, с кластерами пространственной некогерентности, в которых элементы совершают несинхронные хаотические колебания. Данный тип поведения аналогичен виртуальной химере, обнаруженной в осцилляторе с запаздыванием [15, 16], где роль пространственной координаты играет интервал запаздывания. В отличие от «классических» химер в ансамблях с нелокальной связью, обнаруженная в (1.4) химера не является стационарной. Границы кластеров вращаются вдоль кольца так, что их ширина остается практически постоянной. Динамические (вращающиеся) химеры, возникающие на основе бегущих волн в кольце осцилляторов, наблюдались и при нелокальном взаимодействии (например, в [12]). В целом такие структуры еще мало изучены. Их исследование и сопоставление со «стационарными» химерами, на наш взгляд, само по себе важно и интересно. Однако главный результат проведенных исследований состоит в том, что была показана возможность возникновения химероподобных структур в распределенной системе с локальным характером взаимодействия.
С учетом результатов, полученных в [14], можно сделать вывод, что нелокальный характер взаимодействия элементов ансамбля не является необходимым условием для возникновения химерных структур. Такие структуры возникают при определенных типах ло-
кального взаимодействия (инерционного или нелинейного). Возникают вопросы: насколько типичны химеры для систем с локальными связями и возможны ли они в случае более простых типов локального взаимодействия (например при диффузионной связи)? Для ответа на эти вопросы необходимы дальнейшие исследования.
Авторы выражают благодарность Владимиру Викторовичу Семёнову за ценные советы в процессе проведения исследований и представления результатов.
Список литературы
[1] Кузнецов С. П. Сложная динамика генераторов с запаздывающей обратной связью // Изв. вузов. Радиофизика, 1982, т. 25, №12, с. 1410-1428.
[2] Abrams D. M., Strogatz S. H. Chimera states for coupled oscillators // Phys. Rev. Lett., 2003, vol. 93, no. 17, 174102, 4 pp.
[3] Abrams M., Mirollo R., Strogatz S.H., Wiley D.A. Solvable model for chimera states of coupled oscillators // Phys. Rev. Lett., 2008, vol. 101, no. 8, 084103, 4 pp.
[4] Arecchi F.T., Giacomelli G., Lapucci A., Meucci R. Two-dimensional representation of a delayed dynamical system // Phys. Rev. A, 1992, vol.45, no. 7, R4225-R4228.
[5] Arecchi F. T., Meucci R., Allaria E., Di Garbo A., Tsimring L. S. Delayed self-synchronization in homoclinic chaos // Phys. Rev. E, 2002, vol.65, no.4, 046237, 4pp.
[6] Baldi P., Atiya A.F. How delays affect neural dynamics and learning // IEEE Trans. Neural Netw., 1994, vol.5, no. 4, pp. 612-621.
[7] Bocharov G. A, Rihan F. A. Numerical modelling in biosciences using delay differential equations // J. Comput. Appl. Math., 2000, vol. 125, no. 1, pp. 183-199.
[8] Dahlem M. A., Hiller G., Panchuk A., Scholl E. Dynamics of delay-coupled excitable neural systems // Internat. J. Bifur. Chaos Appl. Sci. Engrg., 2009, vol.19, no. 2, pp. 745-753.
[9] Ford N. J., Baker C.T.H. Qualitative behaviour and stability of solutions of discretised nonlinear Volterra integral equations of convolution type //J. Comput. Appl. Math., 1996, vol.66, no. 1, pp. 213-225.
[10] Giacomelli G., Politi A. Relationship between delayed and spatially extended dynamical systems // Phys. Rev. Lett., 1996, vol.76, no. 15, pp. 2686-2689.
[11] Giacomelli G., Marino F., Zaks M. A., Yanchuk S. Coarsening in a bistable system with long-delayed feedback // Europhys. Lett., 2012, vol.99, no. 5, 58005, 5 pp.
[12] Jaros P., Maistrenko Yu., Kapitaniak T. Chimera states on the route from coherence to rotating waves // Phys. Rev. E, 2015, vol.91, no. 2, 022907, 5 pp.
[13] Kuramoto Y., Battogtokh D. Coexistence of coherence and incoherence in nonlocally coupled phase oscillators // Nonlinear Phenom. Complex Syst., 2002, vol. 5, no. 4, pp. 380-385.
[14] Laing C.R. Chimera in networks with purely local coupling // Phys. Rev. E, 2015, vol.92, no. 5, 050904, 5 pp.
[15] Larger L., Penkovsky B., Maistrenko Yu. Virtual chimera states for delayed-feedback systems // Phys. Rev. Lett., 2013, vol. 111, no. 5, 054103, 5 pp.
[16] Larger L., Penkovsky B., Maistrenko Yu. Laser chimeras as a paradigm for multistable patterns in complex systems // Nat. Commun., 2015, vol. 6, 7752.
[17] Mao X., Yuan C., Zou J. Stochastic differential delay equations of population dynamics //J. Comput. Appl. Math., 2005, vol. 304, no. 1, pp. 296-320.
[18] Mounier H., Rudolph J. Time delay systems // Control systems, robotics, and automation: Vol. 14 / H.Unbehauen (ed.). Oxford: EOLSS, 2014. 9 pp.
[19] Omelchenko I., Maistrenko Yu., Hovel Ph., Schöll E. Loss of coherence in dynamical networks: Spatial chaos and chimera states // Phys. Rev. Lett., 2011, vol. 106, no. 23, 234102, 4 pp.
[20] Omelchenko I., Riemenschneider B., Hovel Ph., Schöll E. Transition from spatial coherence to incoherence in coupled chaotic systems // Phys. Rev. E, 2012, vol. 85, no. 2, 026212, 9 pp.
[21] Omelchenko I., Omel'chenko O. E., Hovel Ph., Scholl E. When nonlocal coupling between oscillators becomes stronger: patched synchrony or multichimera states // Phys. Rev. Lett., 2013, vol.110, no. 22, 224101, 5 pp.
[22] Omel'chenko O.E. Coherence-incoherence patterns in a ring of non-locally coupled phase oscillators // Nonlinearity, 2013, vol.26, no. 9, pp. 2469-2498.
[23] Panaggio M. J., Abrams D. M. Chimera states: Coexistence of coherence and incoherence in networks of coupled oscillators // Nonlinearity, 2014, vol. 28, no. 3, R67-R87.
[24] Pyragas K. Continuous control of chaos by self-controlling feedback // Phys. Lett. A, 1992, vol. 170, no. 6, pp. 421-428.
[25] Pyragas K. Delayed feedback control of chaos // Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci., 2006, vol.364, no. 1846, pp. 2309-2334.
[26] Scholl E., Hovel Ph., Flunkert V., Dahlem M.A. Time-delayed feedback control: From simple models to lasers and neural systems // Complex time-delay systems: Theory and applications / F.M.Atay (ed.). Berlin: Springer, 2010. P. 85-150.
[27] Scholl E., Hiller G., Hövel Ph., Dahlem M.A. Time-delayed feedback in neurosystems // Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci., 2009, vol.367, no. 1891, pp. 1079-1096.
[28] Schmidt L.T., Schönleber K., Krischer K., Garcia-Morales V. Coexistence of synchrony and incoherence in oscillatory media under nonlinear global coupling // Chaos, 2014, vol.24, no. 1, 013102, 7pp.
[29] Sethia G. C., Senan A., Johnston G.L. Amplitude-mediated chimera states // Phys. Rev. E, 2013, vol.88, no. 4, 042917, 5 pp.
[30] Süli E., Mayers D.F. An introduction to numerical analysis. Cambridge: Cambridge Univ. Press, 2003. 444 pp.
[31] Zakharova A., Kapeller M., Scholl E. Chimera death: Symmetry breaking in dynamical networks // Phys. Rev. Lett., 2014, vol. 112, no. 15, 154101, 5pp.
Chimera regimes in a ring of elements with local unidirectional interaction
Igor A. Shepelev1, Tatiana E. Vadivasova2
1,2Saratov State University Astrahanskaia 83, Saratov, 410012, Russia 1igor_sar@li.ru, 2vadivasovate@yandex.ru
Complex spatial structures, called chimeras, are the subject of considerable recent interest. They consist of stationary areas with coherent and incoherent behavior of neighboring elements. A number of problems related to similar structures have not been solved yet. One of these problems concerns the element interaction in ensembles, when stable chimera structures can be observed. Until quite recently it was assumed that one of the most important conditions for the existence of chimeras is the nonlocal character of interaction. However, this assumption is not exactly correct. Chimeras can be realized for special types of local coupling. So, the chimera examples were obtained in ensembles with inertial local coupling. The additional variable is introduced for a coupling specification. It is given by a linear differential equation. Also, the so-called virtual chimeras exist in oscillators with delayed feedback. This allows one to assume that chimera states can be obtained in a ring of local coupling oscillators with unidirectional interaction, which is inertialess, but has a nonlinear character. This assumption is based on
a qualitative similarity between the behaviors of an oscillator with delay feedback and a ring of the same oscillators with local unidirectional coupling.
The basis of this work is the system with delay feedback, which demonstrates the existence of a virtual chimera. The distributed analog is investigated. It is an oscillator ring with unidirectional nonlinear local coupling.
The existence of chimera structures in the ring were found in the special area of parameter changing via computing simulation. This chimera moves in a ring with constant velocity and is similar to the chimera in the system with delay feedback. The area of chimera existence of parameter variations was studied. Regime diagrams were plotted on the plane of control parameters. The scenario of chimera destruction for the coupling increase was shown.
MSC 2010: 34F05
Keywords: oscillator with delayed feedback, distributed system, spatial structure, chimera, dynamical chaos, local coupling
Received March 03, 2016, accepted April 21, 2016
Citation: Rus. J. Nonlin. Dyn., 2016, vol. 12, no. 2, pp. 197-209 (Russian)
