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

Оригинал: "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 · собственный разбор · 3 мин
Численный анализ запаздывающего уравнения климата впервые показал, как исчезают квазипериодические колебания при складывании торов.

Контекст

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

Методы

Авторы используют модель ГЗТ — запаздывающее дифференциальное уравнение с отрицательной обратной связью и периодическим воздействием, где переменная h описывает отклонение глубины термоклина в восточной части Тихого океана. Для анализа бифуркаций они применяют программный пакет 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 год.

Ключевые числа

  • задержка обратной связи τ_n: 0.9395 лет (около 11 месяцев)
  • резонанс: 2:7 (два цикла обратной связи за семь лет)
  • амплитуда сезонного воздействия c: 2.976 (в безразмерных единицах)
  • период сезонного воздействия: 1 год
  • длительность дрейфа параметра в примере: 2000 лет (медленный) и 100 лет (быстрый)
Понятия
обратная связь Бьеркнеса модель осциллятора с подзарядкой южная осцилляция эль-ниньо модель Джина–Тиммермана теория запаздывающего осциллятора