[논문 리뷰] Multiple timescales and the parametrisation method in geometric singular perturbation theory

Ian Lizarraga, Bob Rink외 1명|arXiv (Cornell University)|2020. 09. 22.

Combustion and flame dynamics참고 문헌 39인용 수 7

한 줄 요약

이 논문은 기하학적 미소 섭동 이론에서 좌표에 종속되지 않는 매개변수화 방법을 제안하여, 공액 방정식을 반복적으로 풀어 임의의 정밀도로 느린 다변량과 그의 빠른 섬유 번들( fibre bundles)을 계산한다. 이 방법은 세 개 이상의 시간스케일을 가진 시스템에서 숨겨진 시간스케일을 드러내며, 반응 네트워크 모델에 성공적으로 적용되어, 시간스케일의 수나 안정집합의 구조에 대한 사전 지식이 없이도 상향식(top-down) 접근을 가능하게 한다.

ABSTRACT

We present a novel method for computing slow manifolds and their fast fibre bundles in geometric singular perturbation problems. This coordinate-independent method is inspired by the parametrisation method introduced by Cabr\'e, Fontich and de la Llave. By iteratively solving a so-called conjugacy equation, our method simultaneously computes parametrisations of slow manifolds and fast fibre bundles, as well as the dynamics on these objects, to arbitrarily high degrees of accuracy. We show the power of this top-down method for the study of systems with multiple (i.e., three or more) timescales. In particular, we highlight the emergence of hidden timescales and show how our method can uncover these surprising multiple timescale structures. We also apply our parametrisation method to several reaction network problems.

연구 동기 및 목표

• 다중 시간스케일을 가진 기하학적 미소 섭동 문제에서 느린 다변량과 빠른 섬유 번들을 좌표에 종속되지 않게 계산하기 위한 방법을 개발하는 것.
• Cabré, Fontich, de la Llave의 매개변수화 방법을 세 개 이상의 시간스케일을 가진 섭동 문제로 확장하는 것.
• 표준 형태의 시스템 (1.1)에서 명백하지 않은 숨겨진 시간스케일을 드러내는 것, 특히 작은 매개변수들이 상호의존적인 시스템에서의 경우.
• 시간스케일의 수나 안정집합 상태에 대한 사전 지식이 없이도 중첩된 불변 다변량과 그 역학을 체계적이고 알고리즘적으로 근사화하는 방법을 제공하는 것.
• 반응 네트워크 문제에 이 방법을 적용하여, 특히 정규 초구형성 상실 상황에서 느리고 초느린 다변량을 식별하는 데서의 유용성을 입증하는 것.

제안 방법

• 느린 다변량과 빠른 섬유 번들을 식별하는 문제를, 좌표 차트에서 모델 시스템의 역학과 연결하는 공액 방정식을 푸는 것으로 재구성한다.
• 반복적 방법을 사용하여 공액 방정식을 풀어 작은 매개변수 ε에 대해 임의의 차수까지 느린 다변량과 빠른 섬유 번들의 매개변수화를 계산한다.
• 매 반복 단계에서 해를 해석적으로 구하는 호모로지 방정식을 풀며, 느린 흐름의 거의 정적 성질 덕분에 해에 대한 명시적 공식을 도출할 수 있다.
• 이 방법은 상향식으로 작동한다: 먼저 주된 느린 다변량 Sε₀를 계산한 후, 체계적으로 중첩된 초느린 다변량 Sε₁, Sε₂ 등과 그에 관련된 빠른 섬유를 드러낸다.
• 두 번째 차수의 느린 다변량 보정과 첫 번째 차수의 빠른 섬유 보정을 위해 기호 계산(예: Mathematica)을 활용하며, 고차수 계산의 효율성을 높이기 위해 자동 미분을 사용할 것을 제안한다.
• 비판적 다변량이 정상 초구형성이며, 매끄러운 사상 φ₀를 통해 임bedded될 수 있는 시스템에 적용 가능하며, 이는 축소된 역학의 국소 좌표 표현을 가능하게 한다.

실험 결과

연구 질문

• RQ1 세 개 이상의 시간스케일을 가진 다중 시간스케일 시스템에서 느린 다변량과 그의 빠른 섬유 번들을 ε에 대해 임의의 정밀도로 어떻게 계산할 수 있는가?
• RQ2 숨겨진 시간스케일이 미소 섭동 시스템에서 어떻게 발생하며, 어떻게 이를 시스템적으로 드러낼 수 있는가?
• RQ3 표준 다중 시간스케일 프레임워크가 실패하는 상호의존적인 작은 매개변수를 가진 시스템으로 이 매개변수화 방법을 확장할 수 있는가?
• RQ4 이 방법은 반응 네트워크 모델에서 어떻게 작동하는가, 특히 느리고 초느린 다변량을 식별하는 데서 어떤 성능을 보이는가?
• RQ5 공액 방정식은 기하학적 미소 섭동 이론에서 좌표에 종속되지 않는 불변 다변량의 계산을 어떻게 가능하게 하는가?

주요 결과

• 매개변수화 방법은 ε에 대해 임의의 차수까지 느린 다변량과 빠른 섬유 번들을 성공적으로 계산하며, 다변량의 두 번째 차수 보정과 섬유의 첫 번째 차수 보정을 기호적으로 계산하였다.
• 이 방법은 숨겨진 시간스케일을 드러내며, 예를 들어 임베딩된 Van der Pol 시스템 (5.1)에서는 느린 다변량 Sε₀ 상에서 리듬성 진동이 나타나며 느린 다변량과 초느린 다변량의 역학을 연결한다.
• 반응 네트워크 예시 (4.10)에서는 안정적인 노드가 존재하며 제3의 시간스케일이 없음을 드러내며, 축소된 문제에서 평형점의 곡선이 존재하지 않음을 통해 숨겨진 시간스케일이 존재하지 않음을 규명한다.
• 의존적인 작은 매개변수를 가진 시스템은 Cardin & Teixeira (2019)의 프레임워크를 넘어서는 복잡한 다중 시간스케일 구조를 보일 수 있으며, 이는 더 일반적인 기하학적 정의가 필요함을 시사한다.
• 특히 접합점 근처가 아닌 영역에서 정규 초구형성 상실 상황에서도 이 방법은 강건하게 작동하며, 느린 다변량과 초느린 다변량 및 그 역학을 여전히 근사할 수 있다.
• 저자들은 Gi와 Hi의 고차수 계산에서 기호 도함수의 계산 부담을 줄이기 위해 자동 미분을 사용할 것을 제안한다.