[논문 리뷰] A simple solver for the fractional Laplacian in multiple dimensions

Victor Minden, Lexing Ying|arXiv (Cornell University)|2018. 02. 11.

Fractional Differential Equations Solutions참고 문헌 56인용 수 9

한 줄 요약

이 논문은 정규 격자에서 단순한 특이성 제거와 트라프레즈로이드 법칙을 사용하여 다차원에서 분수라플라시안의 간단하고 두 번째 차수 정확도를 갖는 이산화를 제시한다. 이를 통해 FFT를 통한 빠른 적용이 가능해진다. 표준 라플라시안을 기반으로 한 전처리된 켈리브 솔버를 통해 분수 확산 문제를 효율적으로 해결할 수 있으며, α가 2에 가까워질 때조차도 최적의 수렴성과 뛰어난 성능을 보였다.

ABSTRACT

We present a simple discretization scheme for the hypersingular integral representation of the fractional Laplace operator and solver for the corresponding fractional Laplacian problem. Through singularity subtraction, we obtain a regularized integrand that is amenable to the trapezoidal rule with equispaced nodes, assuming a high degree of regularity in the underlying function (i.e., $u\in C^6(R^d)$). The resulting quadrature scheme gives a discrete operator on a regular grid that is translation-invariant and thus can be applied quickly with the fast Fourier transform. For discretizations of problems related to space-fractional diffusion on bounded domains, we observe that the underlying linear system can be efficiently solved via preconditioned Krylov methods with a preconditioner based on the finite-difference (non-fractional) Laplacian. We show numerical results illustrating the error of our simple scheme as well the efficiency of our preconditioning approach, both for the elliptic (steady-state) fractional diffusion problem and the time-dependent problem.

연구 동기 및 목표

• 다차원에서 분수라플라시안에 대한 단순하고 정확하며 효율적인 이산화 방법을 개발하기.
• 디리클레 경계 조건이 있는 유한 영역에서 공간 분수 확산 문제를 빠르게 해결할 수 있도록 하기.
• α → 2일 때 이산 분수라플라시안 시스템의 불안정성 문제를 효과적으로 완화하기 위해 효과적인 전처리 전략을 도입하기.
• 특히 영역 경계 근처에서 해의 정규성에 제한이 있는 문제에 대해 단순한 구적법의 적용 범위를 확장하기.

제안 방법

• 분수라플라시안의 초특이적 적분 커널을 정규화하기 위해 테일러 전개를 통한 특이성 제거 기법을 사용한다.
• 정규화된 피적분함수는 등간격의 직각좌표 격자에서 트라프레즈로이드 법칙을 사용하여 이산화하여 이동 불변성을 갖는 이산 연산자를 얻는다.
• 이산 연산자는 빠른 푸리에 변환(FFT)을 통해 효율적으로 적용되며, 복잡도가 O(N log N)이 된다.
• 해결된 선형 시스템을 풀기 위해 표준 유한차분 라플라시안을 기반으로 한 전처리기를 켈리브 부분공간 방법과 함께 사용한다.
• 전처리기는 이산 분수라플라시안을 표준 라플라시안과 저랭크 또는 희소 보정 행렬의 합으로 분해하는 데 활용된다.
• 영역의 여집합에 해당하는 격점들을 제거하여 L자형 영역과 같은 차폐된 직각좌표 격자로의 확장을 실현한다.

실험 결과

연구 질문

• RQ1 특이성 제거를 통해 다차원 분수라플라시안에 대해 단순하고 고차수의 구적법을 개발할 수 있는가?
• RQ2 결과로 얻어진 이산 연산자가 이동 불변성을 유지하여 FFT 기반 적용이 가능한가?
• RQ3 α가 2에 가까워질 때 분수라플라시안의 이산 시스템이 악조건이 되는 문제를 표준 라플라시안 기반 전처리기를 통해 효과적으로 완화할 수 있는가?
• RQ4 특히 영역 경계 근처에서 해의 정규성이 제한된 경우, 이 방법의 수렴 속도는 어떠한가?
• RQ5 전처리된 켈리브 솔버의 성능은 다양한 차원과 분수 차수 α에 따라 어떻게 스케일링되는가?

주요 결과

• 매우 매끄러운 함수(예: u ∈ C⁶(Rᵈ))에 대해서는 격자 간격 h에 대해 두 번째 차수의 수렴을 달성한다.
• 더러운 해(일반적으로 유한 영역에서 관찰되는 해)에 대해서는 ℓ² 노름에서 첫 번째 차수의 수렴을 보인다.
• 유한차분 라플라시안 기반 전처리기는 특히 큰 α 값일 경우 켈리브 반복 횟수를 크게 줄였다.
• 두 차원에서는 α = 1.25에서 1.75 사이에서 전처리된 CG 방법이 12~29회의 반복을 필요로 하였으며, 실행 시간은 O(N log N)으로 스케일링되었다.
• 세 차원에서는 메서드가 여전히 강건하였으며, α = 1.75일 때 반복 횟수는 30~170으로 증가했지만, 전처리 효과는 유지되었다.
• 오차 추정에 대한 리처드슨 보간율은 약 1.0~1.5였으며, 이는 테스트된 격자에서 최적의 수렴을 나타낸다.