[논문 리뷰] A fully discrete low-regularity integrator for the 1D periodic cubic nonlinear Schrödinger equation

Buyang Li, Yifei Wu|arXiv (Cornell University)|2021. 01. 11.

Numerical methods for differential equations참고 문헌 20인용 수 7

한 줄 요약

이 논문은 $ H^1 $ 초기 데이터를 가진 1차원 주기적 입자 비선형 슈뢰딩거 방정식(NLS)에 대해 완전 이산적이고 완전 비상적(fully explicit)인 저규칙성 적분기법을 제안한다. 비틀린 변수(twisted variables), 주파수 절단(frequency truncation), 그리고 FFT 기반 공간 이산화를 통해, CFL 조건 없이 $ L^2 $-노름 오차 한계 $ O(\tau\sqrt{\ln(1/\tau)} + N^{-1}) $를 달성하며, 로그 인자까지 고려한 일차 수렴을 보여준다.

ABSTRACT

A fully discrete and fully explicit low-regularity integrator is constructed for the one-dimensional periodic cubic nonlinear Schrödinger equation. The method can be implemented by using fast Fourier transform with $O(N\ln N)$ operations at every time level, and is proved to have an $L^2$-norm error bound of $O(τ\sqrt{\ln(1/τ)}+N^{-1})$ for $H^1$ initial data, without requiring any CFL condition, where $τ$ and $N$ denote the temporal stepsize and the degree of freedoms in the spatial discretisation, respectively.

연구 동기 및 목표

• 1D 주기적 입자 비선형 슈뢰딩거 방정식에 대해 부드럽지 않은 초기 데이터를 가진 완전 이산적이고 완전 비상적 수치 방법을 개발하는 것.
• 고도의 규칙성 가정이 필요 없이 $ H^1 $ 초기 데이터에 대해 시간과 공간에서 일차 수렴을 달성하는 것.
• 빠른 푸리에 변환(FFT)을 사용하여 계산 효율성이 높은 방법을 구축하는 것 — 각 시간 단계에서 $ O(N\ln N) $의 연산을 필요로 한다.
• CFL 조건이 필요 없이도 저규칙성 가정 하에 최적의 수렴 속도를 유지하는 것.
• 조화 해석 기법을 통해 주파수 도메인에서 공간 및 시간 이산화 오차를 통합적으로 분석하는 것.

제안 방법

• 이 방법은 지수 적분기법 프레임워크 내에서 일致성 오차를 줄이기 위해 비틀린 변수 변환을 사용한다.
• 지수 적분 내에서 저주파수 및 고주파수 성분을 근사하기 위해 각각 다른 조화 해석 기법을 사용한다.
• 공간 이산화에서는 반복적인 주파수 절단과 FFT를 통합한다: $ (4N+1) $-점 FFT를 사용해 두 개의 $ (2N+1) $-항 푸리에 급수의 곱을 계산한 후, 다시 $ (2N+1) $-항으로 절단한다.
• 이 방법은 삼각 함수 보간 오차를 피하고, 공간 이산화 오차가 오직 주파수 절단에서 기인함을 보장한다.
• 이 스킴은 완전 비상적이며 FFT를 사용해 구현되어 각 시간 단계에서 $ O(N\ln N) $의 연산을 필요로 한다.
• 이산 보건 공간 추정과 수치 해의 $ H^1 $-규칙성에 기반해 주파수 도메인에서 안정성과 일致성 분석을 수행한다.

실험 결과

연구 질문

• RQ1 완전 이산적이고 완전 비상적인 적분기법이 $ H^1 $ 초기 데이터를 가진 1D 주기적 입자 비선형 슈뢰딩거 방정식에 대해 $ L^2 $에서 일차 수렴을 달성할 수 있는가?
• RQ2 CFL 조건이 필요 없이 이러한 방법을 구성할 수 있는가?
• RQ3 조화 해석 기법을 통해 주파수 도메인에서 공간 및 시간 이산화 오차를 함께 분석할 수 있는가?
• RQ4 비선형 항에서 FFT를 효율적으로 사용하면서 보간 오차를 유도하지 않을 수 있는가?
• RQ5 초기 데이터의 규칙성이 최소일 경우에도 최적의 수렴 속도를 유지할 수 있는가?

주요 결과

• 제안된 방법은 $ H^1 $ 초기 데이터에 대해 $ O(\tau\sqrt{\ln(1/\tau)} + N^{-1}) $의 $ L^2 $-노름 오차 한계를 달성하며, CFL 조건이 필요하지 않다.
• 수치 실험 결과 $ H^1 $ 초기 데이터에 대해 시간 수렴 속도가 약 $ O(\tau^{1.02}) $로 나타나 이론적 한계와 일치함을 확인하였다.
• $ H^2 $ 초기 데이터의 경우 $ O(\tau + N^{-2}) $ 오차를 달성하여 최적의 공간 수렴을 보였다.
• 공간 이산화 오차는 오직 주파수 절단에서 기인하며, FFT를 통한 보간 오차는 발생하지 않았다.
• $ \tau $와 $ N $에 대해 일관되게 유계인 $ H^1 $-노름을 가지며, 완전 이산 수치 해의 $ H^1 $-노름이 유계이다.
• 이 방법은 계산적으로 효율적이며, FFT를 사용해 각 시간 단계에서 오직 $ O(N\ln N) $의 연산만 필요로 한다.