[논문 리뷰] Numerical approximations for fractional elliptic equations via the method of semigroups

Nicole Cusimano, Félix del Teso외 1명|arXiv (Cornell University)|2018. 12. 04.

Differential Equations and Boundary Problems인용 수 7

한 줄 요약

이 논문은 유계 영역에서 스펙트럼 분수 라플라시안 $(-\Delta_{\mathcal{B}})^s$를 포함하는 분수 타원형 방정식을 해결하기 위한 새로운 수치적 방법을 제시한다. 반응-확산 접근법을 활용하여 열 반군을 유한요소로 이산화하고, 특이 적분을 위한 단조적 구적법을 결합함으로써, 매끄러운 데이터에 대해 최대 $O(h^{4})$의 수렴 속도와 $L^2$ 데이터에 대해 $O(h^{2s})$의 수렴 속도를 달성하며, 딜리클레, 뉴먼, 루빈 경계 조건에 대해 높은 정확도와 유연성을 제공한다.

ABSTRACT

We provide a novel approach to the numerical solution of the family of nonlocal elliptic equations $(-\\Delta)^su=f$ in $\\Omega$, subject to some homogeneous boundary conditions $\\mathcal{B}(u)=0$ on $\\partial \\Omega$, where $s\\in(0,1)$, $\\Omega\\subset \\mathbb{R}^n$ is a bounded domain, and $(-\\Delta)^s$ is the spectral fractional Laplacian associated to $\\mathcal{B}$ on $\\partial \\Omega$. We use the solution representation $(-\\Delta)^{-s}f$ together with its singular integral expression given by the method of semigroups. By combining finite element discretizations for the heat semigroup with monotone quadratures for the singular integral we obtain accurate numerical solutions. Roughly speaking, given a datum $f$ in a suitable fractional Sobolev space of order $r\\geq 0$ and the discretization parameter $h>0$, our numerical scheme converges as $O(h^{r+2s})$, providing super quadratic convergence rates up to $O(h^4)$ for sufficiently regular data, or simply $O(h^{2s})$ for merely $f\\in L^2(\\Omega)$. We also extend the proposed framework to the case of nonhomogeneous boundary conditions and support our results with some illustrative numerical tests.

연구 동기 및 목표

• 유계 영역에서 스펙트럼 분수 라플라시안을 갖는 비국소 타원형 문제를 정확하고 효율적으로 해결하기 위한 수치적 스킴을 개발하는 것.
• 동차 및 비동차 경계 조건을 갖는 $(-\Delta_{\mathcal{B}})^s u = f$ 문제를 해결하는 데 도전하는 것.
• 역 분수 라플라시안의 반군 표현을 활용하여 안정적이고 고차수의 이산화를 가능하게 하는 것.
• 데이터 $f$의 정(regularity)에 따라 변화하는 조건 하에서 제안된 방법의 엄밀한 수렴 속도를 확립하는 것.

제안 방법

• 해는 반군 공식 $(-\Delta_{\mathcal{B}})^{-s}f(x) = \frac{1}{\Gamma(s)}\int_0^\infty e^{t\Delta_{\mathcal{B}}}f(x) \frac{dt}{t^{1-s}}$을 통해 표현되며, 이는 분수 역 라플라시안과 열 반군을 연결한다.
• 유한요소 방법을 사용하여 열 반군 $e^{t\Delta_{\mathcal{B}}}f$를 이산화함으로써, 각 시간 $t$에서의 해의 공간적 근사가 가능해진다.
• 시간에 대한 특이 적분을 근사하기 위해 단조적 구적법을 사용하여 안정성과 고차수 수렴을 보장한다.
• 해를 동차 성분과 보정 성분으로 분해하여, 비동차 경계 조건으로의 확장을 수행한다.
• 유한요소 오차 추정과 구적법 오차 경계를 결합한 수렴 분석을 통해 전반적인 수렴 속도를 도출한다.
• 디리클레, 뉴먼, 루빈 경계 조건에 적용 가능하며, 더 일반적인 연산자로의 확장도 가능하다.

실험 결과

연구 질문

• RQ1 반군 방법을 활용하여 유계 영역에서 다양한 경계 조건을 갖는 스펙트럼 분수 라플라시안에 대해 고차수 수치 스킴을 유도할 수 있는가?
• RQ2 반군 공식의 시간 적분에 대해 유한요소 공간 이산화와 고차수 구적법을 조합할 경우 도달할 수 있는 수렴 속도는 무엇인가?
• RQ3 저정규성 공간인 $L^2(\Omega)$에 속하는 데이터에 대해 이 방법은 어떻게 작동하며, 수렴성이 유지되는가?
• RQ4 비동차 경계 조건으로의 확장을 유지하면서 최적의 수렴 속도를 유지할 수 있는가?
• RQ5 비볼록 다각형 영역에서도 이 방법은 안정적인가? 특이점 근처의 메esh 세분화를 통해 최적 수렴 속도를 복구할 수 있는가?

주요 결과

• 제안된 방법은 데이터 $f \in \mathbb{H}_{\mathcal{B}}^r(\Omega)$에 대해 $O(h^{r+2s})$의 수렴 속도를 달성하며, 충분히 매끄러운 데이터에 대해 최대 $O(h^4)$의 수렴 속도를 보인다.
• $f \in L^2(\Omega)$인 경우, 방법은 $O(h^{2s})$의 수렴 속도를 보이며, 주어진 정규성 조건에서 최적이다.
• 열 반군의 유한요소 근사와 단조적 구적법의 조합은 안정성과 고차수 정확도를 보장한다.
• 해를 조화 보정과 동차 문제로 분해함으로써 비동차 경계 조건으로의 확장을 수행한다.
• 이론적 수렴 결과는 다양한 경계 조건과 데이터 정규성에서 예측된 수렴 속도를 실험적으로 확인하는 수치 시험을 통해 지지된다.
• 이 프레임워크는 다른 2차 타원형 연산자와 경계 조건으로도 일반화 가능하며, 구적법 및 유한요소 구성 요소의 최소한의 수정으로 적용 가능하다.