[논문 리뷰] Uniform error bounds of a time-splitting spectral method for the long-time dynamics of the nonlinear Klein-Gordon equation with weak nonlinearity

Weizhu Bao, Yue Feng외 1명|arXiv (Cornell University)|2020. 01. 29.

Advanced Mathematical Physics Problems참고 문헌 57인용 수 7

한 줄 요약

이 논문은 약한 비선형성($\varepsilon^p$)을 가진 비선형 클라인-고든 방정식에 시간 분할 스펙트럴 방법을 적용한 경우, 모든 $\varepsilon \in (0,1]$ 및 $O(\varepsilon^{-\beta})$ 시간까지의 장시간 시뮬레이션에서 균일한 오차 한계를 확립한다. 비선형성 강도가 $\varepsilon^p$이고 시뮬레이션 시간이 $O(\varepsilon^{-\beta})$에 이르는 조건에서, 큰 시간 스텝을 사용하더라도 시간에 대해 2차 수렴성을 유지함을 증명한다. 이 방법은 공간에서 스펙트럴 정밀도를 달성하고, 모든 $\varepsilon \in (0,1]$에서 균일 수렴성을 보장하며, 실수 및 복소수 진동 해석 영역 모두에서 수치적으로 검증된다.

ABSTRACT

We establish uniform error bounds of time-splitting Fourier pseudospectral (TSFP) methods for the nonlinear Klein--Gordon equation (NKGE) with weak power-type nonlinearity and $O(1)$ initial data, while the nonlinearity strength is characterized by $\varepsilon^{p}$ with a constant $p \in \mathbb{N}^+$ and a dimensionless parameter $\varepsilon \in (0, 1]$, for the long-time dynamics up to the time at $O(\varepsilon^{-β})$ with $0 \leq β\leq p$. In fact, when $0 < \varepsilon \ll 1$, the problem is equivalent to the long-time dynamics of NKGE with small initial data and $O(1)$ nonlinearity strength, while the amplitude of the initial data (and the solution) is at $O(\varepsilon)$. By reformulating the NKGE into a relativistic nonlinear Schrödinger equation, we adapt the TSFP method to discretize it numerically. By using the method of mathematical induction to bound the numerical solution, we prove uniform error bounds at $O(h^{m}+\varepsilon^{p-β}τ^2)$ of the TSFP method with $h$ mesh size, $τ$ time step and $m\ge2$ depending on the regularity of the solution. The error bounds are uniformly accurate for the long-time simulation up to the time at $O(\varepsilon^{-β})$ and uniformly valid for $\varepsilon\in(0,1]$. Especially, the error bounds are uniformly at the second order rate for the large time step $τ= O(\varepsilon^{-(p-β)/2})$ in the parameter regime $0\leβ

연구 동기 및 목표

• 약한 비선형성($\varepsilon^p$)과 $O(1)$ 초기 조건을 갖는 비선형 클라인-고든 방정식(NKGE)의 장시간 역학을 분석하기 위해.
• 모든 $\varepsilon \in (0,1]$ 및 $O(\varepsilon^{-\beta})$까지의 장시간에 걸쳐 유효한 균일한 오차 한계를 확립하기 위해, 시간 분할 푸리에 가짜스펙트럴(TSFP) 방법을 적용한다. 여기서 $0 \leq \beta \leq p$이다.
• 공간 파장이 $O(1)$, 시간 파장이 $O(\varepsilon^\beta)$이고, 파동 속도가 $O(\varepsilon^{-\beta})$인 매우 진동적인 복소 NKGE로 분석을 확장하기 위해.
• 시간 스텝이 $\tau = O(\varepsilon^{-(p-\beta)/2})$로 스케일링되는 경우에도 TSFP 방법이 시간에 대해 2차 수렴성을 유지함을 보여주기 위해.

제안 방법

• 약한 비선형 역학 분석을 단순화하기 위해 NKGE를 상대론적 비선형 슈뢰딩거 방정식으로 재구성하기 위해.
• 재구성된 방정식을 시간과 공간 모두에 대해 시간 분할 푸리에 가짜스펙트럴(TSFP) 방법을 적용하여 이산화하기 위해.
• 수치 해를 유계화하고 $H^1$-노름에서 균일한 오차 추정을 도출하기 위해 수학적 귀납법을 사용하기 위해.
• 오차 한계를 $O(h^m + \varepsilon^{p-\beta}\tau^2)$로 확립하기 위해, 여기서 $h$는 공간 메쉬 크기, $\tau$는 시간 스텝이며, $m \geq 2$는 해의 정칙성에 따라 달라진다.
• 오차 한계가 모든 $\varepsilon \in (0,1]$ 및 장시간 간격 $T = O(\varepsilon^{-\beta})$에 대해 균일하게 유효함을 증명하여, 약한 비선형 영역에서의 강인성을 확보하기 위해.
• 파동 속도가 $O(\varepsilon^{-\beta})$, 공간 파장이 $O(1)$, 시간 파장이 $O(\varepsilon^\beta)$인 복소수이며 매우 진동적인 NKGE로 분석을 확장하기 위해.

실험 결과

연구 질문

• RQ1 약한 비선형성에서 비선형 클라인-고든 방정식(NKGE)의 장시간 시뮬레이션에 대해 시간 분할 스펙트럴 방법에 대해 균일한 오차 한계를 확립할 수 있는가?
• RQ2 시간 스텝이 $\tau = O(\varepsilon^{-(p-\beta)/2})$로 스케일링되는 경우, $0 \leq \beta < p$ 영역에서 TSFP 방법이 시간에 대해 2차 정확도를 유지하는가?
• RQ3 시간 스텝 크기와 총 스텝 수를 고려할 때, $\beta$와 $p$의 다양한 매개변수 영역에서 오차 행동은 어떻게 변화하는가?
• RQ4 파동 속도가 $O(\varepsilon^{-\beta})$, 시간 진동이 $O(\varepsilon^\beta)$인 복소수이며 매우 진동적인 NKGE에 대해 TSFP 방법과 그 오차 한계를 확장할 수 있는가?
• RQ5 모든 $\varepsilon \in (0,1]$에 대해, 특히 $\varepsilon \to 0^+$ 극한에서도 장시간 시뮬레이션에서 오차 한계가 균일하게 유효한가?

주요 결과

• NKGE에 대해 $O(\varepsilon^{-\beta})$까지의 장시간에 걸쳐 모든 $\varepsilon \in (0,1]$에서 균일한 오차 한계 $O(h^m + \varepsilon^{p-\beta}\tau^2)$를 확보한다.
• $0 \leq \beta < p$일 경우, 시간 스텝이 $\tau = O(\varepsilon^{-(p-\beta)/2})$로 크더라도 시간에 대해 2차 수렴성을 유지함을 확인하여, 약한 비선형 영역에서의 강인성을 입증한다.
• 수치 결과는 실수 및 복소 NKGE 모두에서 모든 테스트된 $\varepsilon$ 값(특히 $\varepsilon \ll 1$)에 대해 시간에 대해 2차 수렴성과 스펙트럴 공간 정밀도를 확인한다.
• $\beta = 1$일 경우, $\varepsilon$에 관계없이 항상 균일한 2차 시간 수렴성을 보이며, $\beta = 2$ 및 $\beta = 3$일 경우 각각 $\tau \lesssim \varepsilon^{3/2}$ 및 $\tau \lesssim \varepsilon^3$일 때에만 2차 수렴성이 관찰된다.
• 계산 비용은 $\beta \leq p/2$일 경우 $O(N\ln N)$으로 스케일링되지만, $\beta > p/2$일 경우 $O(N\ln N / \varepsilon^{(3\beta-p)/2})$로 증가하여 매우 진동적인 영역의 과제를 반영한다.
• 모든 $\varepsilon \in (0,1]$에 대해 오차 한계가 균일하게 유효하여, 약한 비선형성 및 소규모 데이터 극한에서의 신뢰성 보장.