[논문 리뷰] Long time error analysis of the fourth-order compact finite difference methods for the nonlinear Klein-Gordon equation with weak nonlinearity

Yue Feng|arXiv (Cornell University)|2020. 03. 09.

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

한 줄 요약

이 논문은 약한 비선형성($\varepsilon \in (0,1]$)을 가지는 비선형 클라인-고든 방정식에 대해 제4차(compact) 유한차분(4cFD) 방법을 제안한다. 비선형성은 $\varepsilon^p$로 스케일링된다. 장기 시간 $O(\varepsilon^{-p})$까지 엄밀한 오차 한계를 도출하였으며, 최적의 정확도를 확보하기 위해 $h = O(\varepsilon^{p/4})$ 및 $\tau = O(\varepsilon^{p/2})$가 필요하다는 것을 보여주며, 제2차 방법 대비 뛰어난 공간 해상도를 보이고 이론적 수렴 속도가 수치적으로도 정확히 확인된다.

ABSTRACT

We present the fourth-order compact finite difference (4cFD) discretizations for the long time dynamics of the nonlinear Klein-Gordon equation (NKGE), while the nonlinearity strength is characterized by $\varepsilon^p$ with a constant $p \in \mathbb{N}^+$ and a dimensionless parameter $\varepsilon \in (0, 1]$. Based on analytical results of the life-span of the solution, rigorous error bounds of the 4cFD methods are carried out up to the time at $O(\varepsilon^{-p})$. We pay particular attention to how error bounds depend explicitly on the mesh size $h$ and time step $τ$ as well as the small parameter $\varepsilon \in (0, 1]$, which indicate that, in order to obtain `correct' numerical solutions up to the time at $O(\varepsilon^{-p})$, the $\varepsilon$-scalability (or meshing strategy requirement) of the 4cFD methods should be taken as: $h = O(\varepsilon^{p/4})$ and $τ= O(\varepsilon^{p/2})$. It has better spatial resolution capacity than the classical second order central difference methods. By a rescaling in time, it is equivalent to an oscillatory NKGE whose solution propagates waves with wavelength at $O(1)$ in space and $O(\varepsilon^p)$ in time. It is straightforward to get the error bounds of the oscillatory NKGE in the fixed time. Finally, numerical results are provided to confirm our theoretical analysis.

연구 동기 및 목표

• 약한 비선형성($\varepsilon^p$로 특징지어지는) 비선형 클라인-고든 방정식(NKGE)에 제4차 유한차분(4cFD) 방법을 적용한 장기 오차 한계를 개발하고 분석하는 것.
• 장기 시간 $O(\varepsilon^{-p})$까지 정확도를 유지하기 위해 필요한 메쉬 전략($h$ 및 $\tau$)을 규명하는 것.
• 4cFD 방법의 $\varepsilon$-스케일러빌리티를 확립하여 기존 제2차 중심차분 방법 대비 향상된 공간 해상도를 보여주는 것.
• 시간 스케일링을 통해 비선형 클라인-고든 방정식의 진동형으로 일반화된 오차 분석을 확장하여 고정 시간 오차 한계를 도출하는 것.
• 이론적 오차 추정치의 정확성을 수치 실험으로 검증하고 그 날카운성(sharpness)을 확인하는 것.

제안 방법

• 약한 비선형성 $\varepsilon^p u^{p+1}$을 가진 NKGE에 제4차 정확도를 가진 공간 스텐실과 제2차 시간 이산화 방법을 적용한 4cFD 방법을 적용한다.
• 원래 NKGE를 $O(1)$의 공간 파장과 $O(\varepsilon^p)$의 시간 파장을 가지는 진동형으로 변환하기 위해 시간 스케일링 변환을 적용한다.
• 장기 시간 $O(\varepsilon^{-p})$까지 4cFD 스킴에 대한 엄밀한 오차 한계를 유도하며, 이는 $h$, $\tau$, $\varepsilon$에 명시적으로 의존한다.
• 오차 성장 억제를 위해 에너지 방법과 점근 전개를 활용한 분석을 수행한다.
• 이 방법을 진동형 NKGE로 확장하여 고정 시간 오차 한계를 간편하게 도출할 수 있도록 한다.
• 수치적 검증을 위해 고정밀 지수파 통합 푸리에 가짜스펙트럴 방법을 기준 해답으로 사용하여 공간 및 시간 오차를 계산한다.

실험 결과

연구 질문

• RQ1 약한 비선형성을 가진 비선형 클라인-고든 방정식에 대해 제4차 유한차분 방법의 장기 오차 거동은 어떠한가?
• RQ2 장기 시간 $O(\varepsilon^{-p})$까지 정확도를 유지하기 위해 필요한 메쉬 크기 $h$와 시간 간격 $\tau$는 무엇인가?
• RQ3 이 문제에 대해 4cFD 방법은 기존 제2차 중심차분 방법 대비 공간 해상도에서 어떤가?
• RQ4 시간 스케일링을 통해 도출된 진동형 NKGE에 대해 오차 한계를 확장할 수 있는가?
• RQ5 수치 실험을 통해 이론적 오차 추정치가 날카로운가?

주요 결과

• 고정된 $\varepsilon > 0$에 대해 4cFD 방법은 공간에서 제4차 수렴, 시간에서 제2차 수렴을 보이며, 이는 정확도 설계 순서를 확인한다.
• 매우 진동성이 강한 영역($\varepsilon \ll 1$)에서는 $h \lesssim \varepsilon^{p/4}$ 및 $\tau \lesssim \varepsilon^{p/2}$일 때에만 최적의 수렴이 관찰되며, 이는 이론적 $\varepsilon$-스케일러빌리티를 검증한다.
• 필요한 메쉬 크기 $h = O(\varepsilon^{p/4})$는 기존 제2차 중심차분 방법($h = O(\varepsilon^{p/2})$ 필요)보다 더 뛰어난 공간 해상도를 제공한다.
• 수치 결과로 수렴 속도가 이론적 예측과 일치함을 확인하여 오차 한계가 날카롭다는 것을 입증한다.
• 진동형 NKGE의 형태는 고정 시간 오차 한계를 직접 도출할 수 있게 하여 장기 동역학 분석을 단순화한다.
• 수치 실험에서 $p=1$ 및 $p=2$의 경우, 유도된 메쉬 전략에 따라 4cFD 방법이 $O(\varepsilon^{-p})$까지 장기적으로도 높은 정확도를 유지함을 확인한다.