[논문 리뷰] Uniformly accurate exponential-type integrators for klein-gordon equations with asymptotic convergence to classical splitting schemes in the nonlinear schroedinger limit

Simon Baumstark, Erwan Faou외 1명|arXiv (Cornell University)|2016. 06. 15.

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

한 줄 요약

이 논문은 상대론적($c=1$) 및 고주기 비상대론적($c\gg1$) 모든 영역에서 시간에 대해 균일하게 정확한 지수형 적분법을 케링-고든 방정식에 도입한다. 이는 점근 전개나 $c$에 의존하는 스텝 크기 제약 없이 제1차 및 제2차 수렴성을 유지한다. $c\to\infty$일 때는 비선형 슈뢰딩거 방정식의 고전적 리와 스트랑 분할 방법으로 수렴하며, 최소한의 정규성 조건 하에서 균일 정확성을 달성한다.

ABSTRACT

We introduce efficient and robust exponential-type integrators for Klein-Gordon equations which resolve the solution in the relativistic regime as well as in the highly-oscillatory non-relativistic regime without any step-size restriction, and under the same regularity assumptions on the initial data required for the integration of the corresponding limit system. In contrast to previous works we do not employ any asymptotic/multiscale expansion of the solution. This allows us derive uniform convergent schemes under far weaker regularity assumptions on the exact solution. In particular, the newly derived exponential-type integrators of first-, respectively, second-order converge in the non-relativistic limit to the classical Lie, respectively, Strang splitting in the nonlinear Schr{ö}dinger limit.

연구 동기 및 목표

• 고주기 비상대론적 영역($c\gg1$)에서 스텝 크기 제약 없이 케링-고든 방정식을 수치적으로 해결하는 데 도전한다.
• 기존 방법들이 점근 전개나 다스케일 전개에 의존하고 $c$에 의존하는 안정성 제약을 겪는 한계를 극복한다.
• 시간에 대해 균일하게 수렴하는 적분법을 개발하고, 비선형 슈뢰딩거 극한에서 고전적 분할 방법으로 수렴하도록 한다.
• 이전의 균일 정확한 방법보다 더 약한 정규성 조건 하에서도 최적의 수렴 속도를 달성한다.
• 제1차 및 제2차 방법 모두에서 $c\to\infty$ 극한에서 고전적 분할 방법(Lie 및 Strang)과 일致하도록 보장한다.

제안 방법

• 해의 점근 전개를 피하기 위해 제2차 형식의 케링-고든 방정식에 대한 지수형 적분법을 제안한다.
• 방정식의 구조에 기반한 새로운 분할 전략을 사용하여 고주기 항을 직접 지수 적분법으로 처리한다.
• 제1차 및 제2차 방법(식 28 및 81)을 도출하여, $H^7$ 및 $H^6$ 정규성 조건 하에서 모든 $c\geq1$에 대해 시간에 대해 균일하게 수렴함을 보장한다.
• 고전적 리 및 스트랑 분할 방법이 비선형 슈뢰딩거 극한 시스템(2)에 대해 $c\to\infty$일 때 수렴함을 입증한다.
• 공간 이산화에 푸리에 가짜스펙트럴 방법을 적용하여 $K=2^{10}$ 모드를 사용해 공간에서 스펙트럼 정밀도를 확보한다.
• 지수 적분법의 수정된 형태를 적용하여 $c$에 의존하는 분산 관계를 $\langle\nabla\rangle_c = (I - c^{-2}\Delta)^{1/2}$를 통해 포함시킨다.

실험 결과

연구 질문

• RQ1 케링-고든 방정식에 대해 $c$의 모든 값, 특히 비상대론적 극한에서 균일 수렴성을 유지하는 지수형 적분법을 구성할 수 있는가?
• RQ2 이러한 적분법이 $c\to\infty$일 때 비선형 슈뢰딩거 극한에서 고전적 리 및 스트랑 분할 방법으로 수렴하는가?
• RQ3 다스케일 또는 점근 전개 기반 방법이 요구하는 것보다 더 약한 정규성 조건 하에서도 균일 수렴을 달성할 수 있는가?
• RQ4 지수 적분법 방법과 고전적 분할 방법 사이의 관계는 $c\to\infty$ 영역에서 어떻게 되는가?
• RQ5 수치 실험에서 다양한 영역에서 수렴 차수와 $c$-독립성 측면에서 방법의 성능은 어떠한가?

주요 결과

• 제1차 지수 적분법(28)은 $H^7$ 해에 대해 $\sup_t \|z(t)\|_{H^7} + c^{-2}\|\partial_t z(t)\|_{H^6} \leq 1$ 조건 하에 $c$에 대해 균일하게 $\mathcal{O}(\tau)$ 수렴성을 $H^1$ 노름에서 확보한다.
• 제2차 지수 적분법(81)은 충분히 스무스한 해에 대해 $H^7$ 공간에서 $c$에 대해 균일하게 $\mathcal{O}(\tau^2)$ 수렴성을 확보한다.
• 제2차 방법(81)은 $c\to\infty$일 때 비선형 슈뢰딩거 극한 시스템(2)의 고전적 스트랑 분할 방법으로 수렴하며, 오차는 $\|u_*^n - u_{*,\infty}^n\|_r \leq c^{-1}k_r(M_3,T)$를 만족한다. 조건은 $\tau c \geq 1$이다.
• 수치 실험은 $c = 1, 5, 10, \dots, 10000$ 범위에서 $H^1$ 노름에서 제1차 및 제2차 수렴성을 확인하며, 기울기가 이론적 순서와 일致한다.
• 스텝 크기 제약 없이 $c=10000$일 때에도 방법이 안정적이고 정확하게 유지되며, 고주기 영역에서의 강건성을 입증한다.
• 기존의 가우츠키 유형 방법이 요구하는 $c^4\tau^2$ 전역 오차를 피하며, 이는 $c^2\tau < 1$ 조건이 필요하기에 이전 방법의 핵심 한계를 극복한다.