Исследование диффузионного хаоса
Сидоров С.В.
На примере решения уравнения Курамото-Цузуки приведены результаты численных исследований диффузионного хаоса в окрестности термодинамической ветви системы уравнений "реакция-диффузия". Рассмотрены сценарии перехода к хаосу при решении второй краевой задачи. Показано, что переход к хаосу в фазовом пространстве уравнения Курамото-Цузуки осуществляется через каскады бифуркаций удвоения периода и через субгармонические каскады бифуркаций рождения устойчивых двумерных инвариантных торов, как по внутренней, так и по внешней частоте. Приведено сравнение сценариев перехода к хаосу в фазовом пространстве и в пространстве коэффициентов Фурье.
Ключевые слова: диффузионный хаос; исследование; реакция; моделирование.
Введение
Для моделирования процессов и явлений, описывающих эволюцию распределенных динамических систем широко используются параболические уравнения и системы уравнений вида
u Rm, (1)
где D матрица коэффициентов диффузии, оператор Лапласа, F(u,) векторная в общем случае нелинейная по u функция, скалярный параметр, m размерность системы.
Большое место в многочисленных моделях физики, химии, биологии, а также при моделировании экономических, экологических и социальных задач занимают так называемые уравнения «реакция-диффузия» системы из двух нелинейных дифференциальных уравнений в частных производных
, (2)
диффузионный хаос уравнение
где d1 и d2 коэффициенты диффузии, f() и g() нелинейные относительно переменных X и Y скалярные функции.
Как правило, система уравнений вида (2) имеет, по крайней мере, одно стационарное пространственно однородное и асимптотически устойчивое решение (X0, Y0) при значении параметра < 0, где < 0 значение параметра, соответствующее равновесному состоянию системы. Это решение называют термодинамической ветвью. В окрестности термодинамической ветви при значениях > 0 поведение решений определяется спектром линеаризованной на решении (X0, Y0) краевой задачи в окрестности точки бифуркации 0. В том случае, когда при = 0 только одно простое собственное значение проходит через ноль, возникает пространственно неоднородное стационарное решение - стационарные диссипативные структуры. Эта бифуркация открыта А. Тьюрингом [1] и носит его имя. Если же при = 0 два комплексно сопряженных собственных значения оператора линеаризации пересекают мнимую ось, то имеет место бифуркация рождения пространственно однородных периодических колебаний - бифуркация Андронова-Хопфа [2, 3]. При значительном удалении значений параметра от 0 в системе (2) возникают более сложные, в том числе пространственно неоднородные и непериодические решения, которые получили общее название “диффузионный хаос”. К настоящему времени известно, что вид этих решений зависит от значений коэффициентов диффузии, от соотношения между ними, от граничных условий, а также от формы и размеров пространственной области, в которой ищется решение системы. Однако, механизм образования хаотических решений, теория перехода к диффузионному хаосу и, следовательно, структура решений в системах вида (2) практически не изучены.
Значительный шаг в изучении решений таких систем был сделан в работе [4], где удалось показать, что любое решение системы уравнений реакция-диффузия, возникающее в окрестности термодинамической ветви в результате ее бифуркации при м > м0, удовлетворяет уравнению
(3)
где W W(x,t) = u(x,t) + iv(x,t) комплекснозначная функция, описывающая отклонение от термодинамической ветви, c1 и c2 действительные постоянные, значения которых определяются коэффициентами d1 и d2, функциями f(u, v, м), g(u, v, м) и их производными, вычисленными на термодинамической ветви.
Уравнение (3), называемое уравнением Курамото-Цузуки или зависящим от времени уравнением Гинзбурга-Ландау, играет важную роль в изучении и понимании процессов, происходящих в нелинейных диссипативных средах диффузионного типа, в неравновесных системах, в плоских течениях и т.д. Уже в ранних исследованиях решений уравнения (3) группой А. Самарского и С. Курдюмова [57] особое внимание было обращено на условия появления сложных пространственно неоднородных и непериодических решений и на изучение их свойств. Однако в указанных работах рассматривалась более простая, хотя и обладающая хаотической динамикой, система из трех нелинейных обыкновенных дифференциальных уравнений, полученная в результате редукции уравнения (3) и представления его в так называемом «маломодовом» приближении
(4)
В последнее время появились веские основания для сомнений в правомерности переноса результатов, справедливых для конечномерных систем маломодовых приближений, на бесконечномерные системы. В частности, в работах [8-10] показано, что трехмерные нелинейные диссипативные системы обыкновенных дифференциальных уравнений имеют один общий сценарий перехода к хаосу. Началом любого сценария является каскад бифуркаций удвоения периода циклов - гармонический каскад Фейгенбаума [11], который сходится к хаотическому аттрактору Фейгенбаума. Более сложная структура сингулярного (хаотического) аттрактора создается субгармоническим каскадом бифуркаций рождения устойчивых циклов, период которых определяется согласно порядку Шарковского [12]. Дальнейшее усложнение хаотических аттракторов, рождающихся в точках накопления значений бифуркационного параметра, идет через гомоклинический каскад бифуркаций рождения устойчивых циклов, сходящихся к гомоклиническому контуру - петле сепаратрисы особой точки типа седло-фокуса [10]. Вместе с тем, в системах большей размерности в сценарии перехода к хаосу участвуют двумерные торы и, вполне вероятно, весь субгармонический каскад бифуркаций двумерных торов по одной из частот подобно тому, как это имеет место в пятимерной системе уравнений Лоренца [13].
Каким образом происходит переход к хаосу в уравнении Курамото-Цузуки и какова структура сингулярных (хаотических) аттракторов уравнения (3), пока неизвестно. В настоящей работе приведены результаты численных исследований, которые частично дают ответы на поставленные вопросы.
Порядок Шарковского - это упорядочение в множестве натуральных чисел, имеющее вид
.
В [10] показано, если диссипативная система нелинейных обыкновенных дифференциальных уравнений имеет цикл периода n, то она имеет также цикл любого периода n', такой что в смысле порядка Шарковского.
Переход к хаосу в пространстве коэффициентов Фурье для уравнения Курамото-Цузуки.
При изучении решений уравнения (3) наибольший интерес представляет вторая краевая задача
(5)
Для дальнейшего исследования задачи (5) надо пояснить, каким образом получена система (4) и что по существу представляют решения этой системы. Идея получения уравнений (4) основана на использовании Галеркинских маломодовых аппроксимаций, а именно, на разложении решения W(x,t) = u(x,t) + iv(x,t) уравнения (3) в ряд Фурье. Для второй краевой задачи (5) это разложение имеет вид
.
Так как коэффициенты Фурье достаточно быстро убывают с ростом их номера, то в маломодовом приближении учтены только первые гармоники, т.е.
. (6)
Подстановка (6) в (5) и отбрасывание затем всех членов с более высокими гармониками, содержащих множители , приводит к замкнутой системе обыкновенных дифференциальных уравнений относительно переменных a0(t), a1(t), b0(t), b1(t).
Система (4) получается путем следующей замены переменных:
i 2(t) = ai2(t) + bi 2(t), i = 0, 1,
= 0 2, = 12, = 0 1 ,
с учетом того, что ai = i cos i, bi = i sin i. Таким образом, упрощенная система (4) по существу отражает не решения задачи (5), а решения некоторой другой конечномерной системы уравнений в пространстве коэффициентов Фурье. В этом отношении система (4) аналогична хорошо известной системе Лоренца из трех нелинейных обыкновенных дифференциальных уравнений [14]. В связи с этим здесь уместно отметить, что переход к хаосу в пятимерной системе Лоренца [13] осуществляется через каскад бифуркаций не предельных циклов, а двумерных инвариантных торов по внешней частоте. Причем, сценарий перехода к хаосу в пятимерной системе более простой, чем в трехмерной системе он не содержит гомоклинических каскадов бифуркаций.
Ранее [10] при исследовании механизмов образования хаотических решений для уравнения Курамото-Цуцзуки в маломодовом приближении (4) было показано, что сценарий перехода к хаосу в этой системе не отличается от сценария, общего для всех трехмерных нелинейных диссипативных систем обыкновенных дифференциальных уравнений. Однако существование общего сценария образования хаотических аттракторов не означает, что последние должны иметь одинаковую структуру. Системы дифференциальных уравнений отличаются топологией векторного поля, количеством и типом неподвижных точек, имеют различные области диссипативности в фазовом пространстве и много других отличий. Поэтому не все из указанных выше каскадов бифуркаций удвоения периода циклов, субгармонических и гомоклинических каскадов бифуркаций могут принимать участие в образовании хаотического аттрактора, если последний имеет место. Помимо этого, каскады бифуркаций, порождающих хаотические аттракторы, могут быть как полными, так и неполными, могут прерываться из-за влияния соседних неподвижных точек или других близко расположенных в фазовом пространстве решений. Все это создает большое разнообразие структур хаотических аттракторов.
В частности, возможны ситуации, когда неполным является каскад бифуркаций удвоения периода предельного цикла, и хаотический аттрактор в этом случае отсутствует. Не повторяя результатов исследования системы (4), изложенных в монографии [10], приведем пример неполного каскада бифуркаций удвоения периода предельного цикла в упрощенной системе (4) при фиксированном значении параметров c1 = 2,5, k = 1 и изменении параметра c2 на множестве отрицательных значений.
В области значений c2 [2,87; 0) система (4) имеет неподвижную точку, которая при c2 2,88 теряет устойчивость, и в результате бифуркации Андронова-Хопфа рождается устойчивый предельный цикл (рис. 1, а).
Рис. 1. Проекции на плоскость (0, 1) пространства коэффициентов Фурье решений задачи (4) при значении с1 = 2,5: а) предельный цикл при c2 = 3; б) цикл периода 2 при c2 = 4,16; в) цикл периода 8 при c2 = 4,46.
Далее при c2 4,116 в системе (4) рождается цикл удвоенного периода, при c2 4,37 цикл периода 4, при c2 4,45 цикл периода 8. Однако, как показывает численный эксперимент, с дальнейшим уменьшением значений параметра c2 идет обратный процесс, т.е. цикл периода 8 при значении c2 4,8 становится неустойчивым, и снова появляется цикл периода 4. Затем при c2 5,0 рождается цикл периода 2, который при c2 6,08 также теряет устойчивость, и в системе (4) появляется предельный цикл, а при значении c2 28,3 устойчивая неподвижная точка.
Совсем иначе выглядят каскады бифуркаций решений задачи (5) на отрезке l = (k = 1) в бесконечномерном пространстве коэффициентов Фурье в тех же переменных 0 и 1 и при том же фиксированном значение параметра c1 = 2,5. При значениях параметра c2 (1,851; 0) система имеет неподвижную точку (0, 0, 0, …,0), которой соответствует пространственно однородное периодическое решение задачи (5). При этом траектория в фазовом пространстве имеет форму окружности, что свидетельствует об одинаковой амплитуде колебаний переменных u(t) и v(t), определяемой значением координаты 0. При значении c2 1,851 пространственно однородное периодическое решение теряет устойчивость, и рождается другое устойчивое, но пространственно неоднородное периодическое решение. В пространстве коэффициентов Фурье этой бифуркации соответствует потеря устойчивости точки (0, 0, 0, …,0) и появление другой неподвижной устойчивой точки с отличными от нуля значениями переменных i.
Эта точка остается устойчивой до значения параметра c2 2803, соответствующего бифуркации рождения предельного цикла (рисунок 2а). При величине c2 3,472 происходит бифуркация удвоения периода данного цикла (рисунок 2б), и начинается каскад бифуркаций Фейгенбаума. Так при значении c2 = 3,628 наблюдается цикл периода 4, при c2 = 3,635 цикл периода 8, при c2 = 3,645 цикл периода 16 и т.д. Каскад бифуркаций удвоения периода данного цикла, завершается появлением аттрактора Фейгенбаума при значении параметра c2 3,66 (рисунок 2в).
Рис. 2. Проекции на плоскость (0, 1) пространства коэффициентов Фурье решений задачи (5) при значении с1 = 2,5: а) предельный цикл при значении c2 = 3; б) цикл периода 2 при c2 = 3,5; в) аттрактор Фейгенбаума при c2 = 3,66.
С уменьшением значений параметра c2 решения в пространстве коэффициентов Фурье имеют значительно более сложный хаотический характер. Чтобы понять структуру образовавшихся хаотических аттракторов, рассмотрим каскады бифуркаций решений задачи (5) при значениях параметра c2 > 14,95. При меньших значениях параметра c2 единственным аттрактором в пространстве коэффициентов Фурье является цикл (рис. 3, а), а при значении c2 14,95 происходит бифуркация удвоения периода этого цикла, которая дает начало каскаду бифуркаций рождения циклов удвоенного периода. Так при значении c2 11,92 рождается цикл периода 4, в области значений c2 [11,31; 11,175] лежит цикл периода 8, при значениях c2 [11,31; 11,175] -цикл периода 16, при c2 [11,174; 11,148] -цикл периода 32, в области c2 [11,147 11,141] - цикл периода 64 и т.д. Этот каскад бифуркаций завершается аттрактором Фейгенбаума при значении параметра c2 11,13 (рис.3, в).
Рис.3. Проекции на плоскость (0, 1) пространства коэффициентов Фурье решений задачи (5) при значении с1 = 2,5: а) предельный цикл при c2 = 20; б) цикл периода 2 при c2 = 12,6; в) аттрактор Фейгенбаума при c2 = 11,13.
За аттрактором Фейгенбаума с ростом значений параметра c2 следует субгармонический каскад бифуркаций рождения циклов с периодами, определяемыми порядком Шарковского, о чем свидетельствуют устойчивые циклы с периодом 6 (рисунок 4а) при c2 [10,845; 10,825] и с периодом 5 (рис.4 б) при c2 [10,226; 10,215], а также каскады бифуркаций Фейгенбаума, порожденных этими циклами. Однако в данной области параметров (c1, c2) субгармонический каскад бифуркаций является неполным, так как отсутствует цикл периода 3. При значении c2 9,914 траектория попадает в область притяжения другого решения, субгармонический каскад бифуркаций, порожденный рассматриваемым циклом, обрывается и в системе (5) появляется хаотическое поведение (рисунок 4в).
Рис. 4. Проекции на плоскость (0, 1) пространства коэффициентов Фурье решений задачи (5) при значении с1 = 2,5: а) цикл периода 6 при c2 = 10,82; б) цикл периода 5 при c2 = 10,21; в) хаотический аттрактор при c2 = 9,914.
Наблюдаемое хаотическое поведение в области c2 [9,913; 7,917] может быть обусловлено слиянием решений, порожденных каскадом бифуркаций цикла, рассмотренного выше и цикла, устойчивого при значениях параметра c2 [7,916; 7,884] (рисунок 5а).
Рис. 5. Проекции на плоскость (0, 1) пространства коэффициентов Фурье решений задачи (5) при значении с1 = 2,5: а) цикл при значении c2 = 7,9; б) цикл периода 2 при значении c2 = 7,87; в) цикл периода 3 при значении c2 = 7,849.