Как ломается оркестр: пузыри Шенсинера в климатической модели ✎

Оригинал: "Chenciner bubbles and torus break-up in a periodically forced delay differential equation"
· Andrew Keane, Bernd Krauskopf
arXiv:1708.02334 · 2017-08-07 · arXiv non-exclusive · собственный разбор · 2 мин
Численный эксперимент впервые показал, как в уравнении с запаздыванием исчезают квазипериодические колебания, когда два тора складываются и пропадают.

Дирижёр оркестра постепенно теряет контроль над темпом. Сначала вторая скрипка чуть отстаёт, потом виолончели ускоряются, и вот уже каждый музыкант играет в своём ритме, хотя ещё минуту назад все следовали одной мелодии. Примерно так в климатической системе Эль-Ниньо две независимые частоты — собственный ритм обратной связи и навязанный сезонный цикл — могут сосуществовать годами, а затем резко рассогласоваться.

Внутри этой кухни работает модель ГЗТ: простое уравнение с запаздыванием, где скорость изменения глубины термоклина зависит от её значения в прошлом. Когда параметры попадают в резонансный язык 2:7 (два цикла обратной связи за семь лет), в фазовом пространстве возникают торы — поверхности, по которым траектория движется, как нить по бублику, наматывая витки с двумя разными периодами.

Тор — это геометрическое место точек, где система может находиться бесконечно долго, не повторяясь точно, но и не убегая в хаос.

Теоретики давно предсказали, что при изменении параметров два таких тора могут сблизиться и исчезнуть, оставив после себя область сложной динамики — пузырь Шенсинера. Но одно дело нарисовать схему на бумаге, и совсем другое — увидеть, как она разворачивается в уравнении с бесконечномерным фазовым пространством. Авторы работы использовали программный пакет DDE-Biftool, чтобы продолжать периодические решения и вычислять их мультипликаторы Флоке — числа, показывающие, насколько устойчив каждый цикл.

Внутри пузыря для резонанса 2:7 обнаружилась структура, хорошо согласующаяся с теоретической картиной Баесенса и Маккея: кривые нейтральных сёдел, гомоклинических и гетероклинических переходов, а также кривые складывания торов. Но нашлось и отличие: бифуркация тора оказалась субкритической, то есть рождался седловой тор, а не устойчивый. Это значит, что при пересечении критической линии система не плавно переходит на новый режим, а срывается в резкое переключение, как если бы оркестр не постепенно расстроился, а разом перешёл на какофонию.

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

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

Результаты могут помочь в интерпретации климатических временных рядов и в поиске предвестников резких переходов в системах с несколькими частотами, хотя прямого практического применения пока нет. Дирижёр, потерявший контроль, — это не только красивая метафора, но и предупреждение: если мы научимся распознавать момент, когда торы готовы сложиться, то, возможно, успеем расслышать первые фальшивые ноты в климатическом оркестре.

Резонансный язык 2:7 в модели означает, что за 7 лет совершается ровно 2 цикла обратной связи, и именно этот язык оказался самым крупным в рассматриваемой области параметров, что облегчило численный анализ.

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

\dot{h}(t) = -b \tanh[\kappa h(t-\tau_n)] + c \cos(2\pi t)
Скорость изменения глубины термоклина зависит от её значения в прошлом (задержка τ_n) и от сезонного воздействия с периодом 1 год.
Понятия
обратная связь Бьеркнеса модель осциллятора с подзарядкой южная осцилляция эль-ниньо модель Джина–Тиммермана теория запаздывающего осциллятора