[논문 리뷰] A PDE approach to fractional diffusion in general domains: a priori error analysis

Ricardo H. Nochetto, Enrique Otárola외 1명|arXiv (Cornell University)|2013. 02. 04.

Numerical methods in engineering참고 문헌 29인용 수 9

한 줄 요약

이 논문은 일반적인 유계 영역에서 타원형 연산자의 분수거듭제곱을 포함하는 분수 확산 문제를 해결하기 위한 PDE 기반 접근법을 제안한다. 무한 길이 실린더로 확장된 문제와 가중 Sobolev 프레임워크를 사용하여, 유한 요소 이산화에 대한 사전 오차 추정을 유도한다. 이는 최적의 정(regularity)을 보이지만, 준균일 메쉬에서는 최적의 수렴 차수를 넘지 못하고, 이sovotropic 메쉬에서는 준최적의 수렴 차수를 보인다.

ABSTRACT

The purpose of this work is the study of solution techniques for problems involving fractional powers of symmetric coercive elliptic operators in a bounded domain with Dirichlet boundary conditions. These operators can be realized as the Dirichlet to Neumann map for a degenerate/singular elliptic problem posed on a semi-infinite cylinder, which we analyze in the framework of weighted Sobolev spaces. Motivated by the rapid decay of the solution of this problem, we propose a truncation that is suitable for numerical approximation. We discretize this truncation using first degree tensor product finite elements. We derive a priori error estimates in weighted Sobolev spaces. The estimates exhibit optimal regularity but suboptimal order for quasi-uniform meshes. For anisotropic meshes, instead, they are quasi-optimal in both order and regularity. We present numerical experiments to illustrate the method's performance.

연구 동기 및 목표

• 비국소 연산자(예: 분수 라플라시안)에 의해 지배되는 분수 확산 문제를 해결하기 위한 강력한 수치적 방법을 개발하는 것.
• 분수 연산자의 비국소성 문제를 해결하기 위해 문제를 확장된 영역에서 국소적이며 탈구조적/특이적인 타원형 문제로 변환하는 것.
• 확장된 문제의 유한 요소 이산화의 수렴 성질을 분석하고, 특히 가중 Sobolev 공간에서의 오차 추정에 초점을 맞추는 것.
• 준균일 메쉬와 이sovotropic 메쉬의 성능을 비교하여, 확장된 문제의 해를 구할 때 최적의 수렴 속도를 달성하는 데에 어떤 영향을 미치는지 분석하는 것.
• Caffarelli-Silvestre 확장 프레임워크를 통해 라플라시안을 초월하여 일반적인 두 번째 차수, 대칭적, 균일 타원형 연산자로 이 접근법을 일반화하는 것.

제안 방법

• Caffarelli-Silvestre 확장을 사용하여 반무한 실린더 $\mathcal{C} = \Omega \times (0, \infty)$ 상에서 분수 확산 문제를 Dirichlet-to-Neumann 맵 문제로 공식화한다.
• 가중 Sobolev 공간에서 가중치 $y^\alpha$, $\alpha = 1 - 2s$를 포함하는 탈구조적/특이적인 타원형 PDE로 확장된 문제를 모델링한다.
• 해의 $y$-방향에서의 지수적 감쇠 특성을 활용하여, 유한 높이 $\mathcal{C}_\mathcal{Y}$에서 무한 실린더를 잘라내어 수치적 근사를 가능하게 한다.
• 확장된 영역의 형태가 규칙적인 삼각분할 위에서 첫 번째 차수의 텐서 곱 유한 요소를 사용하여 문제를 이산화한다.
• 확장된 문제에 대한 사전 오차 추정을 가중 $H^1$-노름에서 유도하며, 오차를 메쉬 크기와 요소 수에 따라 연결한다.
• 준균일 및 이sovotropic 메쉬의 정밀도 개선에 따른 수렴 행동을 분석하여, 이sovotropic 메쉬에서 더 향상된 수렴 속도를 보임을 보여준다.

실험 결과

연구 질문

• RQ1 일반적인 유계 영역에서 타원형 연산자의 분수거듭제곱을 국소 PDE 공식화를 통해 어떻게 수치적으로 근사할 수 있는가?
• RQ2 준균일 메쉬와 이sovotropic 메쉬를 사용할 때, 확장된 문제의 유한 요소 근사의 수렴 행동은 어떻게 되는가?
• RQ3 반무한 실린더의 잘라내기 작업이 수치 해의 정확성에 어떤 영향을 미치는가?
• RQ4 이 방법을 통해 원래의 분수 확산 문제의 해에서 최적의 정규성과 준최적의 수렴 차수를 달성할 수 있는가?
• RQ5 오차 추정은 기저 타원형 연산자의 스펙트럼 성질과 가중치 $y^\alpha$에 어떻게 의존하는가?

주요 결과

• 확장된 문제에 대해 최적의 정규성을 갖는 오차 추정을 달성하였으며, 해는 가중 $H^1$-공간 $\widetilde{H}^1_L(\mathcal{C}, y^\alpha)$에 속한다.
• 준균일 메쉬의 경우 수렴 차수는 최적이 아니며, $|\log(\#\mathscr{T}_\mathcal{Y})|^s (\#\mathscr{T}_\mathcal{Y})^{-1/(n+1)}$의 비율로 스케일링된다.
• 메쉬의 비대칭성 $\gamma > 3/(1 - \alpha)$ 조건을 만족하는 이sovotropic 메쉬의 경우, 오차 추정은 차수와 정규성 면에서 모두 준최적이며, 최상의 가능한 수렴 속도를 충족한다.
• 잘라내기 오차는 잘라낸 높이 $\mathcal{Y}$에 따라 지수적으로 감쇠되며, $\|\nabla(\mathfrak{u} - v)\|_{L^2(\mathcal{C}, y^\alpha)} \lesssim e^{-\sqrt{\lambda_1}\mathcal{Y}/4}\|f\|_{\mathbb{H}^s(\Omega)'}$의 형태로 표현된다.
• 유한 요소 해 $V_{\mathscr{T}_\mathcal{Y}}$는 가중 에너지 노름에서 진짜 해 $\mathfrak{u}$로 수렴하며, 오차는 메쉬에 의존하는 항에 의해 유계화되어 있다.
• 이 방법은 라플라시안을 초월하여 일반적인 두 번째 차수, 대칭적, 균일 타원형 연산자 $\mathcal{L}$로도 일반화 가능하며, 일반화된 Caffarelli-Silvestre 확장을 통해 가능하다.