[논문 리뷰] A Fourier Pseudospectral Method for the "Good" Boussinesq Equation with Second Order Temporal Accuracy

Kelong Cheng, Wenqiang Feng외 2명|arXiv (Cornell University)|2014. 01. 24.

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

한 줄 요약

이 논문은 '좋은' Boussinesq 방정식에 대해 두 번째 차수 정확도를 갖는 완전 이산화된 푸리에 가짜스펙트럴 방법을 제시한다. 여기서 시간 도함수 $u_t$ 를 근사하기 위해 보조 변수 $\psi$ 를 도입함으로써 무조건적 안정성과 향상된 수렴 속도를 달성한다. 이는 이전 연구에서 존재하던 심각한 시간 스텝 제한 $\Delta t \leq Ch^2$ 를 극복하며, $u$ 에 대해 $\ell^{\infty}(0,T^{*};H^{2})$ 오차 추정과 $\psi$ 에 대해 $\ell^{\infty}(0,T^{*};\ell^{2})$ 오차 추정을 달성한다.

ABSTRACT

In this paper, we discuss the nonlinear stability and convergence of a fully discrete Fourier pseudospectral method coupled with a specially designed second order time-stepping for the numerical solution of the "good" Boussinesq equation. Our analysis improves the existing results presented in earlier literature in two ways. First, an $l_\infty(0, T^*; H2)$ convergence for the solution and $l_\infty(0, T^*; l_2)$ convergence for the time-derivative of the solution are obtained in this paper, instead of the $l_\infty(0, T^*; l_2)$ convergence for the solution and the $l_\infty(0, T^*; H^{-2})$ convergence for the time-derivative, given in [17]. In addition, the stability and convergence of this method is shown to be unconditional for the time step in terms of the spatial grid size, compared with a severe restriction time step restriction $Δt \leq Ch^2$ reported in [17].

연구 동기 및 목표

• 두 번째 차수 시간 정확도를 갖는 '좋은' Boussinesq 방정식에 대한 완전 이산화된 푸리에 가짜스펙트럴 스킴을 개발하는 것.
• 이전 연구 [17]에서 존재하던 심각한 시간 스텝 제한 $\Delta t \leq Ch^2$ 를 극복하여, 미세 격자에서의 효율성을 향상시키는 것.
• 보조 변수 $\psi = u_t$ 를 도입함으로써 $\ell^{\infty}(0,T^{*};H^{2})$ 오차 추정을 증명함으로써 $u$ 에 대해 향상된 수렴 속도를 확립하는 것.
• 보조 변수 $\psi$ 를 도입하고 새로운 앨리아싱 오차 제어 기법을 사용함으로써 무조건적 안정성을 달성하는 것.
• 기준 솔로itary 웨이브 해를 통해 수치적으로 두 번째 차수 시간 정확도와 스펙트럴 공간 정확도를 검증하는 것.

제안 방법

• 시간 도함수 $u_t$ 를 근사하기 위해 보조 변수 $\psi$ 를 도입하여, 이차 편미분 방정식을 시간에 대해 일阶계 시스템으로 재구성하는 것.
• 중심 차분을 사용한 $\psi$ 와 중점 근사법을 사용한 $u$ 를 포함하는 두 번째 차수의 Strang 분할 유사 시간 이산화를 적용하여, 두 번째 차수 시간 정확도를 확보하는 것.
• 주기적 경계 조건을 갖는 푸리에 가짜스펙트럴 공간 이산화를 사용하여 공간에서 스펙트럴 정확도를 달성하는 것.
• 비선형 항을 다루고 $H^2$ 노름에서 고차수 수렴을 가능하게 하기 위해 새로운 앨리아싱 오차 제어 추정(보조정리 2.1)을 활용하는 것.
• 에너지 유형 추정을 $L^2$ 및 $H^2$ 노름을 사용하여 도출함으로써, 시간 스텝 제약 없이 비선형 안정성과 수렴성을 증명하는 것.
• 이산 $H^2$ 노름과 푸리에 스펙트럴 노름 사이의 등가성을 활용하여 오차 분석에서 고차수 도함수를 제어하는 것.

실험 결과

연구 질문

• RQ1 완전 이산화된 푸리에 가짜스펙트럴 방법이 '좋은' Boussinesq 방정식에 대해 심각한 시간 스텝 제약 없이 두 번째 차수 시간 정확도를 달성할 수 있는가?
• RQ2 수렴 분석을 개선하여 $u$ 에 대해 $\ell^{\infty}(0,T^{*};H^{2})$ 오차 추정과 $\psi = u_t$ 에 대해 $\ell^{\infty}(0,T^{*};\ell^{2})$ 오차 추정을 달성할 수 있는가?
• RQ3 새로운 공식화와 분석을 통해 이전 연구에서 존재하던 $\Delta t \leq Ch^2$ 시간 스텝 제약을 제거할 수 있는가?
• RQ4 보조 변수 $\psi$ 의 도입이 가짜스펙트럴 스킴의 안정성과 수렴성에 어떤 영향을 미치는가?
• RQ5 수치적으로 확인한 바, 제안된 스킴의 시간 및 공간 정확도의 실제 차수는 얼마인가?

주요 결과

• 제안된 방법은 해 $u$ 에 대해 $\ell^{\infty}(0,T^{*};H^{2})$ 수렴과 시간 도함수 $\psi = u_t$ 에 대해 $\ell^{\infty}(0,T^{*};\ell^{2})$ 수렴을 달성하며, 이는 [17]에서 보고된 $\ell^{\infty}(0,T^{*};\ell^{2})$ 와 $\ell^{\infty}(0,T^{*};H^{-2})$ 수렴보다 향상된 결과이다.
• 시간 스텝 $\Delta t$ 와 공간 격자 크기 $h$ 간의 관계에 대해 무조건적 안정성을 확보하여, 이전 연구에서 존재하던 제한적인 $\Delta t \leq Ch^2$ 조건이 필요 없어졌다.
• 수치 실험 결과, $u$ 와 $\psi$ 에 대해 모두 시간에 대해 두 번째 차수 수렴이 확인되었으며, 오차는 $\mathcal{O}(\Delta t^2)$ 의 스케일링을 보였다. 이는 오차 곡선이 $C N_K^{-2}$ 를 잘 따르는 것으로 확인되었다.
• 공간 수렴은 스펙트럴 정확도를 보이며, 그리드 포인트 수 $N$ 이 증가함에 따라 오차가 $N$ 에 대해 어떤 다항식보다도 더 빠르게 감소한다.
• 앨리아싱 오차 제어 기법(보조정리 2.1)은 고차수 $H^2$ 수렴 추정을 가능하게 하고, 음의 노름 추정이 필요 없도록 하는 데 핵심적인 역할을 한다.
• 중간 변수 $\psi$ 의 사용은 안정성 조건을 향상시키는 것 외에도 일致성과 오차 분석을 단순화하여, 더 강력한 비선형 안정성 증명을 가능하게 하였다.