[논문 리뷰] A first-order Fourier integrator for the nonlinear Schrödinger equation on T without loss of regularity.

Yifei Wu, Fangyan Yao|arXiv (Cornell University)|2020. 10. 06.

Advanced Mathematical Physics Problems참고 문헌 32인용 수 5

한 줄 요약

이 논문은 원환면 위의 1차원 입자형 비선형 슈뢰딩거 방정식에 대해, 빠른 푸리에 변환을 활용한 효율적 구현을 위한 1차 수준의 명시적 푸리에 적분법을 제안한다. 이는 추가적인 정규성 조건 없이도 $H^\gamma$에서 $\gamma > \frac{3}{2}$인 모든 $\gamma$에 대해 1차 수렴성을 확보하며, 질량은 $O(\Delta t^5)$ 오차로 거의 정확히 보존된다. 이는 낮은 정규성 조건 하에서도 장시간 안정성과 정확도를 확보한 중요한 진전이다.

ABSTRACT

In this paper, we propose a first-order Fourier integrator for solving the cubic nonlinear Schrodinger equation in one dimension. The scheme is explicit and can be implemented using the fast Fourier transform. By a rigorous analysis, we prove that the new scheme provides the first order accuracy in $H^\gamma$ for any initial data belonging to $H^\gamma$, for any $\gamma >\frac32$. That is, up to some fixed time $T$, there exists some constant $C=C(\|u\|_{L^\infty([0,T]; H^{\gamma})})>0$, such that $$ \|u^n-u(t_n)\|_{H^\gamma(\mathbb T)}\le C au, $$ where $u^n$ denotes the numerical solution at $t_n=n au$. Moreover, the mass of the numerical solution $M(u^n)$ verifies $$ \left|M(u^n)-M(u_0) ight|\le C au^5. $$ In particular, our scheme dose not cost any additional derivative for the first-order convergence and the numerical solution obeys the almost mass conservation law. Furthermore, if $u_0\in H^1(\mathbb T)$, we rigorously prove that $$ \|u^n-u(t_n)\|_{H^1(\mathbb T)}\le C au^{\frac12-}, $$ where $C= C(\|u_0\|_{H^1(\mathbb T)})>0$.

연구 동기 및 목표

• 원환면 위의 입자형 비선형 슈뢰딩거 방정식을 위한 수치적 스킴을 개발하여, 추가 정규성 조건 없이도 $H^\gamma$에서 $\gamma > \frac{3}{2}$에 대해 1차 수렴성을 유지한다.
• 낮은 정규성 초기 데이터 조건에서도 질량을 거의 정확히 보존하는 수치적 해를 확보하여, $O(\Delta t^5)$ 오차를 갖는다.
• 초기 데이터가 $H^\gamma$ 및 $H^1$에 속할 경우, 소볼레프 노름에서 최적의 오차 한계를 확립하는 엄밀한 수렴 분석을 제공한다.

제안 방법

• 스킴은 공간에서 빠른 푸리에 변환을 통한 스펙트럼 체로케이션 방식으로 시간에 따라 해를 진행하는 1차 푸리에 적분법에 기반한다.
• 선형 및 비선형 항을 분리하는 스플리팅 전략을 사용하여, 푸리에 공간에서 명시적인 시간 적분이 가능하도록 한다.
• 비선형 슈뢰딩거 방정식의 구조를 유지하도록 설계되어, 장시간에 걸쳐 안정성과 정확도를 확보한다.
• 에너지 추정과 $H^\gamma$ 노름에서의 교환자 추정을 사용하여 엄밀한 오차 분석을 수행하여 수렴 속도를 확립한다.
• 비선형 상호작용 항은 비선형성의 구조를 활용하여 고주파 성분을 철저히 제어함으로써 다루어진다.
• 스킴이 질량을 $O(\Delta t^5)$ 오차 이내로 보존함을 입증하였으며, 이는 장시간 시뮬레이션에서 핵심적인 특징이다.

실험 결과

연구 질문

• RQ1 원환면 위의 1D 입자형 비선형 슈뢰딩거 방정식에 대해, 추가 정규성 조건 없이도 $\gamma > \frac{3}{2}$에 대해 $H^\gamma$에서 최적의 수렴성을 확보하는 1차 푸리에 적분법을 구성할 수 있는가?
• RQ2 제안된 스킴은 낮은 정규성 초기 데이터 조건에서도 질량 보존을 고차수, 즉 $O(\Delta t^5)$ 수준으로 유지하는가?
• RQ3 $H^1(\mathbb{T})$에 속하는 초기 데이터에 대해 $H^1$에서의 수렴 속도는 얼마이며, 이는 이론적 기대와 일치하는가?

주요 결과

• 스킴은 $H^\gamma(\mathbb{T})$에서 임의의 $\gamma > \frac{3}{2}$에 대해 1차 수렴성을 확보하며, 오차는 $C \Delta t$ 이하로 유계이며, 여기서 $C$는 해의 $L^\infty([0,T]; H^\gamma)$ 노름에 의존한다.
• 수치적 해의 질량은 $|M(u^n) - M(u_0)| \le C \Delta t^5$를 만족하여, 고차수 정확도로 근사 질량 보존을 보여준다.
• 초기 데이터가 $H^1(\mathbb{T})$에 속할 경우, 수렴 속도는 임의의 $\varepsilon > 0$에 대해 $O(\Delta t^{1/2 - \varepsilon})$이며, 주어진 정규성 조건 하에서 최적이다.
• 수렴 결과를 확보하기 위해 $H^\gamma$ 이외의 추가 도함수 정규성 조건이 필요하지 않으며, 이는 이전 스킴에 비해 핵심적인 개선이다.
• 스킴은 명시적이며, 빠른 푸리에 변환을 활용하여 효율적으로 구현 가능하여 실용적인 장시간 시뮬레이션을 가능하게 한다.
• 엄밀한 분석을 통해, 낮은 정규성 초기 데이터 조건 하에서도 고정된 시간 간격 동안 안정성과 정확도를 유지함을 확인하였다.