[논문 리뷰] Resonances in long time integration of semi linear Hamiltonian PDEs

Erwan Faou, Benoît Grébert|ArXiv.org|2009. 04. 09.

Numerical methods for differential equations참고 문헌 10인용 수 4

한 줄 요약

이 논문은 반선형 해밀턴 PDE에 대한 표준 심플렉틱 적분법이 선형 연산자의 고주파 성분이 잘리지 않은 이상 수치적 공진현상에 시달릴 수 있음을 보여준다. 주파수 컷오프를 도입하고 분할법을 사용함으로써, 날카로운 CFL 유사 조건 하에서 행동량과 소볼레프 노름의 장기적 보존을 증명한다. 이는 모든 시간 간격에서 부당한 공진현상을 피하는 데 성공한다.

ABSTRACT

We consider a class of Hamiltonian PDEs that can be split into a linear unbounded operator and a regular non linear part, and we analyze their numerical discretizations by symplectic methods when the initial value is small in Sobolev norms. The goal of this work is twofold: First we show how standard approximation methods cannot in general avoid resonances issues, and we give numerical examples of pathological behavior for the midpoint rule and implicit-explicit integrators. Such phenomena can be avoided by suitable truncations of the linear unbounded operator combined with classical splitting methods. We then give a sharp bound for the cut-off depending on the time step. Using a new normal form result, we show the long time preservation of the actions for such schemes for all values of the time step, provided the initial continuous system does not exhibit resonant frequencies.

연구 동기 및 목표

• 반선형 해밀턴 PDE의 장기적 시뮬레이션에서 중점 규칙 및 IMEX 스킴과 같은 표준 심플렉틱 적분법이 수치적 공진현상으로 인해 실패할 수 있는 이유를 규명하는 것.
• 시간 간격과 고유주파수 간의 비선형 관계로 인해 수치적 공진현상이 발생하며, 특히 $ h\tilde{\theta} \to 2\tilde{\theta}\tilde{\theta}\tilde{\theta} $ 일 때 $ \tilde{\theta} \neq 0 $ 이면 물리적 보존 성질이 손실된다는 것을 보여주는 것.
• 선형 연산자의 고주파 모드를 잘라내기 전에 분할법을 적용함으로써 이러한 부당한 공진현상을 제거할 수 있음을 입증하는 것.
• 시간 간격 $ h $ 에 따라 달라지는 컷오프 주파수에 대한 날카로운 경계를 유도하여 장기적 안정성과 행동량 보존을 보장하는 것.
• 원래 시스템에 비공진 조건이 만족될 경우, 수정된 스킴이 물리적 보존 법칙을 도입하지 않고도 유지함을 증명하는 것.

제안 방법

• 해밀턴 PDE를 Sobolev 공간 프레임워크 내에서 $ H = H_0 + P $ 로 형식화하며, $ H_0 $ 는 선형 비유계 연산자이고 $ P $ 는 정규 비선형성임을 가정한다.
• 선형 연산자 $ H_0 $ 의 스펙트럼에 주파수 컷오프 $ K $ 를 도입하여 $ |j| > Kh^{-1} $ 일 때 $ \tilde{\theta}_j = 0 $ 으로 설정함으로써 고주파 성분을 제거한다.
• 분할법 형태인 $ \tilde{\theta}_K(hH_0) \times \text{flow of } P $ 를 적용하며, 여기서 $ \tilde{\theta}_K $ 는 선형 흐름의 잘린 형태이다.
• 소수의 나눗셈을 제어하기 위해 정규형 변환을 사용하며, 비공진 조건을 통해 유계성을 확보한다. 이는 $ |e^{ih\tilde{\theta}^h(\boldsymbol{j})} - 1| $ 의 형태를 갖는다.
• homological 방정식 $ (e^{ih\tilde{\theta}^h(\boldsymbol{j})} - 1)\tilde{\theta}_{\boldsymbol{j}} + hZ_{\boldsymbol{j}} = hG_{\boldsymbol{j}} $ 에 대한 추정을 유도하며, 이는 $ \tilde{\theta}^h(\boldsymbol{j}) $ 와 $ \boldsymbol{j} $ 의 색인 구조에 따라 달라진다.
• 정규형의 나머지 항들이 $ \tilde{\theta}^{r+1} $ 의 순서를 가지며, 초기 자료 크기가 $ \tilde{\theta} $ 일 때 $ \tilde{\theta}^{-r} $ 까지의 장기적 안정성을 보장한다.

실험 결과

연구 질문

• RQ1 표준 심플렉틱 적분법, 예를 들어 중점 규칙 또는 IMEX 스킴이 반선형 해밀턴 PDE의 장기적 통합에서 수치적 공진현상을 피할 수 있는가?
• RQ2 시간 간격 $ h $ 가 CFL 조건에 의해 제약을 받지 않을 경우, 심플렉틱 스킴에서 수치적 공진현상은 무엇으로 인해 발생하는가?
• RQ3 선형 연산자 $ H_0 $ 의 고주파 모드를 잘라내는 것은 분할법의 장기적 행동에 어떤 영향을 미치는가?
• RQ4 행동량과 소볼레프 노름의 장기적 보존을 보장하기 위해 시간 간격 $ h $ 에 따라 결정되는 컷오프 주파수 $ K $ 의 날카로운 경계는 무엇인가?
• RQ5 정규형 논증을 사용하여 심플렉틱 적분법의 비선형 위상 관계에서 발생하는 소수의 나눗셈을 제어할 수 있는가?

주요 결과

• 표준 심플렉틱 적분법, 예를 들어 중점 규칙이나 IMEX 스킴에서 $ h\tilde{\theta}^h(\boldsymbol{j}) \to 2\tilde{\theta}\tilde{\theta}\tilde{\theta} $ 일 때 $ \tilde{\theta} \neq 0 $ 이면 원래 시스템이 비공진이더라도 수치적 공진현상이 발생할 수 있다.
• 전체 시스템에 중점 규칙을 적용할 경우, 시간 간격 $ h $ 에 따라 물리적 보존 법칙을 파괴하거나 예상치 못하게 행동량을 유지하는 인위적 공진현상이 발생할 수 있다.
• 선형 연산자 $ H_0 $ 를 주파수 $ |j| > Kh^{-1} $ 에서 잘라내는 것은 부당한 공진현상을 제거하고 장기적 안정성을 복원한다.
• 시간 간격 $ h $ 에 따라 결정되는 날카로운 경계 $ K \thicksim h^{-1} $ 가 도출되었으며, $ K \thicksim h^{-1} $ 일 때 초기 자료 크기가 $ \tilde{\theta} $ 인 경우 행동량과 소볼레프 노름이 $ \tilde{\theta}^{-r} $ 까지 보존된다.
• 정규형의 동차 방정식은 $ |e^{ih\tilde{\theta}^h(\boldsymbol{j})} - 1| \thicksim h\tilde{\theta}^h(\boldsymbol{j}) $ 일 때 해를 갖는다. 이는 컷오프 조건과 비공진 가정 하에서 성립한다.
• 정규형의 나머지는 $ \tilde{\theta}^{r+1} $ 의 순서를 가지며, 변환은 $ \tilde{\theta}_{\boldsymbol{j}} \notin \tilde{\theta}_{\boldsymbol{j}} $ 의 범주에 머물며 장기적 제어를 보장한다.