[논문 리뷰] Fast Fourier-like Mapped Chebyshev Spectral-Galerkin Methods for PDEs with Integral Fractional Laplacian in Unbounded Domains

Changtao Sheng, Jie Shen외 3명|arXiv (Cornell University)|2019. 08. 27.

Fractional Differential Equations Solutions참고 문헌 51인용 수 5

한 줄 요약

이 논문은 덴포르-타일러 표현과 푸리에 유사 이중수직 매핑 체비셰프 함수(MCFs)를 사용하여 무한도 영역 내 적분 분수 라플라스 연산자를 갖는 PDE를 해결하기 위한 빠른 스펙트럼-갈레르킨 방법을 제안한다. 이 방법은 $ O((N\log_2 N)^d) $의 준최적 복잡도를 달성하고 분수 라플라스 연산자를 완전히 대각화하여 고차 정확도를 확보하며 비선형 분수 PDE(예: 분수 비선형 셰딩거 방정식)의 효율적 해법을 가능하게 한다.

ABSTRACT

In this paper, we propose a fast spectral-Galerkin method for solving PDEs involving integral fractional Laplacian in $\mathbb{R}^d$, which is built upon two essential components: (i) the Dunford-Taylor formulation of the fractional Laplacian; and (ii) Fourier-like bi-orthogonal mapped Chebyshev functions (MCFs) as basis functions. As a result, the fractional Laplacian can be fully diagonalised, and the complexity of solving an elliptic fractional PDE is quasi-optimal, i.e., $O((N\log_2N)^d)$ with $N$ being the number of modes in each spatial direction. Ample numerical tests for various decaying exact solutions show that the convergence of the fast solver perfectly matches the order of theoretical error estimates. With a suitable time-discretization, the fast solver can be directly applied to a large class of nonlinear fractional PDEs. As an example, we solve the fractional nonlinear Schr{ö}dinger equation by using the fourth-order time-splitting method together with the proposed MCF-spectral-Galerkin method.

연구 동기 및 목표

• 무한도 영역 내 적분 분수 라플라스 연산자를 포함하는 편미분방정식을 효율적이고 정확하게 해결하기 위한 수치적 방법을 개발하는 것.
• 다중 차원에서 분수 라플라스 연산자의 비국소성과 특이성으로 인한 계산적 과제를 극복하는 것.
• 감소하는 해를 갖는 문제에서 스펙트럼 정확도를 유지하면서 준최적 복잡도를 달성하는 것.
• 선형 시스템을 완전히 대각화함으로써 비선형 분수 PDE의 효율적 시간 적분을 가능하게 하는 것.

제안 방법

• 이 방법은 분수 라플라스 연산자의 덴포르-타일러 적분 표현을 활용하여 리졸베이트 연산자를 포함하는 경로 적분으로 표현한다.
• 해의 무한도 영역 내 감쇠 성질을 다루기 위해 설계된 이중수직 매핑 체비셰프 함수(MCFs)를 기저 함수로 사용한다.
• MCFs의 이중수직성 덕분에 스펙트럼-갈레르킨 프레임워크 내에서 분수 라플라스 연산자의 약한 형태가 완전히 대각화된다.
• 결과적으로 얻어진 선형 시스템은 대각 구조를 가지며, 이는 $ N $을 차원당 모드 수로 하여 $ O((N\log_2 N)^d) $의 복잡도로 해를 구할 수 있게 한다.
• 이 방법은 4차 시간 분할 스킴을 통해 시간에 의존하는 문제로 확장되어 고차 정확도를 유지한다.
• 해의 알려진 감쇠 정확해를 갖는 문제에 대한 광범위한 수치 실험을 통해 방법의 타당성이 검증되었다.

실험 결과

연구 질문

• RQ1 $ \mathbb{R}^d $ 내에서 적분 분수 라플라스 연산자는 무한도 영역에서 스펙트럼 방법을 통해 효율적이고 정확하게 근사될 수 있는가?
• RQ2 분수 라플라스 연산자의 비국소성과 특이성은 스펙트럼-갈레르킨 프레임워크 내에서 완전히 대각화될 수 있는가?
• RQ3 감쇠하는 해를 갖는 분수 PDE를 위한 빠른 스펙트럼 방법의 도달 가능한 수렴 속도와 계산 복잡도는 무엇인가?
• RQ4 제안된 방법은 고차 정확도를 유지하면서 비선형 분수 PDE로 확장될 수 있는가?

주요 결과

• 제안된 방법은 $ O((N\log_2 N)^d) $의 준최적 복잡도를 달성하여 기존 스펙트럼 방법 대비 계산 비용을 크게 감소시킨다.
• 수치 실험 결과에서 수렴 속도가 이론적 오차 추정치와 일치함을 확인하여 감쇠하는 해에 대해 스펙트럼 정확도를 입증한다.
• 분수 라플라스 연산자가 MCF 기저에서 완전히 대각화되어 각 시간 단계에서 선형 시스템을 효율적으로 해석할 수 있다.
• 4차 시간 분할 스킴을 사용하여 분수 비선형 셰딩거 방정식에 성공적으로 적용되었으며, 고차 정확도를 유지하였다.
• 이론적 분석과 수치 결과는 다차원 공간 문제에 대해 이 방법의 강건성과 효율성을 확인한다.