[논문 리뷰] User guide on Hopf bifurcation and time periodic orbits with pde2path

Hannes Uecker|arXiv (Cornell University)|2019. 08. 02.

Nonlinear Dynamics and Pattern Formation참고 문헌 29인용 수 6

한 줄 요약

이 논문은 유한 요소를 사용하는 PDE를 위한 MATLAB 도구상자인 pde2path에서 히프 분기 분석 및 주기적 궤도 연속을 구현하기 위한 종합적인 사용자 가이드를 제시한다. 주기적 해를 탐지하고 추적하기 위한 수치적 방법을 상세히 기술하며, 대칭을 고려한 연속, 플로케트 다항식을 통한 분기 탐지, O(2) 대칭, 질량 보존 및 단계 조건을 갖는 시스템을 위한 전용 루틴을 포함하여, PDE에서 히프 분기 분 branches 및 2차 분기의 안정적인 계산을 가능하게 한다.

ABSTRACT

We explain the setup for using the pde2path libraries for Hopf bifurcation and continuation of branches of periodic orbits and give implementation details of the associated demo directories. See [Uecker, Comm. in Comp. Phys., 2019] for a description of the basic algorithms and the mathematical background of the examples. Additionally we explain the treatment of Hopf bifurcations in systems with continuous symmetries, including the continuation of traveling waves and rotating waves in O(2) equivariant systems as relative equilibria, the continuation of Hopf bifurcation points via extended systems, and some simple setups for the bifurcation from periodic orbits associated to critical Floquet multipliers going through +-1.

연구 동기 및 목표

• 시간 주기적 해를 갖는 PDE에서 히프 분기를 탐지하고 연속하는 데에 pde2path를 실용적이고 구현 중심으로 사용할 수 있도록 안내하는 것.
• 지속적인 대칭을 갖는 시스템, 예를 들어 O(2) 등장성, 질량 보존, 단계 불변성 등을 다룰 수 있도록 pde2path의 기능을 확장하는 것.
• 히프 분기점과 주기적 궤도에서의 분기 전환을 가능하게 하여, 주기 분할 및 히프-피치포크 분기를 포함한 연속을 가능하게 하는 것.
• 표준 초기값 해법이 실패하는 비산산성 및 잘못된 설정의 PDE(예: 후방 확산을 포함한 최적 제어 문제)에서도 계속 방법이 유용한 바, 이를 지원하는 것.
• 주기적 궤도 계산을 위한 고급 루틴, 예를 들어 호 길이 연속, 메쉬 적응, 후행 다항식 계산을 문서화하고 표준화하는 것.

제안 방법

• PDE를 대규모 ODE 시스템으로 유한 요소 공간 반할산분해하여 pde2path의 내장된 해법을 통해 정적 상태 및 주기적 궤도의 연속을 가능하게 한다.
• 시간에 대한 선형화 방법을 적용하고 암시적 시간 적분 및 호 길이 연속을 사용하여 주기적 궤도를 추적하며, 효율적인 계산을 위해 조화 균형 방법과 TOM(시간 연속 방법)을 활용한다.
• 확장된 시스템을 통해 히프 분기점을 찾기 위해 고유값 조건과 단계 조건을 포함한 증강된 시스템을 풀이한다.
• 플로케트 이론을 적용하여 단조행렬과 그 고유값(플로케트 다항식)을 계산하며, 직접 계산 및 후행 계산을 위한 루틴을 제공한다.
• 대칭 시스템에서 발생하는 단계 모호성을 해결하기 위해 단계 조건(예: 적분 또는 콜로케이션 기반)을 도입하여 상대 평형 상태 및 주기적 궤도의 안정적 연속을 가능하게 한다.
• 예측-수정 기반의 스킴을 사용하여 확장된 시스템에 대해 뉴턴 반복을 수행하며, 확장된 시스템의 야코비안을 활용하고 시간에 대한 메쉬 적응을 통한 정밀도 향상을 도모한다.

실험 결과

연구 질문

• RQ1 유한 요소 방법을 사용하여 복잡한 시공간 역학을 보이는 PDE에서 히프 분기를 신뢰성 있게 탐지하고 연속하는 데에 어떤 수치 전략이 필요한가?
• RQ2 지속적인 대칭을 갖는 시스템(예: 이동 또는 회전 불변성)에서 주기적 궤도를 연속하기 위해 필요한 수치 전략은 무엇인가?
• RQ3 플로케트 다항식 분석을 통해 주기적 궤도에서 발생하는 분기, 예를 들어 주기 분할 및 히프-피치포크 분기를 어떻게 탐지하고 국소화할 수 있는가?
• RQ4 비산산성 또는 잘못된 설정의 PDE(예: 후방 확산을 포함한 최적 제어 문제)를 다루기 위해 연속 프레임워크에 어떤 수정이 필요한가?
• RQ5 대칭 시스템에서 주기적 궤도의 연속을 안정화하기 위해 단계 조건을 어떻게 적절히 설정하고 구현할 수 있는가?

주요 결과

• pde2path 프레임워크는 1D, 2D, 3D PDE에서 히프 분기의 탐지 및 연속을 성공적으로 가능하게 하였으며, 복잡한 기니즈부르크-랜드의 방정식 및 브루셀레이터 시스템을 포함한다.
• O(2) 대칭을 갖는 시스템에서는 이동파와 정적파를 구분할 수 있으며, 상대 평형 상태 및 상대 주기적 궤도를 정확히 계산한다.
• 단계 조건의 사용은 질량 보존 및 이동 불변성 시스템(예: 쿠라모토-시바시크스키 및 FHN 유형 모델)에서 주기적 궤도의 안정적 연속을 가능하게 한다.
• 정규형 계수를 활용한 예측-수정 기반 스킴을 통해 히프 궤도에서의 분기 전환을 실현하였으며, 주기 분할 및 히프-피치포크 분기를 탐지할 수 있었다.
• 후행 플로케트 다항식 계산(floqap) 및 메쉬 적응(hopftref)의 구현은 장시간 주기적 궤도 연속에서 정확성과 견고성을 향상시켰다.
• 이론적으로 잘못된 설정의 PDE, 예를 들어 최적 오염 제어 문제에서도 표준 시간 적분이 실패하지만, 연속 방법은 여전히 유용하므로 본 방법은 효과적이다.