[논문 리뷰] Regularity of the Solution to 1-D Fractional Order Diffusion Equations

Vincent J. Ervin, Norbert Heuer외 1명|arXiv (Cornell University)|2016. 07. 30.

Fractional Differential Equations Solutions참고 문헌 10인용 수 4

한 줄 요약

이 논문은 1차원에서 분수형 확산 연산자의 새로운 해석을 제안하며, 이는 리만-리오빌리 또는 카푸토가 아니며, 일阶 도함수와 분수형 적분의 병합을 통해 정의됨을 주장한다. 이 연산자의 핵을 규명하고, 잭비 다항식이 의사 고유함수임을 증명하며, 가중 잭비 함수를 사용한 스펙트럼 방법을 개발하여 $ (N+2)^{-\beta} $ 비례하는 스펙트럼 수렴 속도를 달성한다. 여기서 $ \beta = \alpha $이며, $ \alpha \in (1,2) $ 이다. 이 방법은 가중 $ L^2 $ 공간에서 최적의 수렴성과 역행성 보장한다.

ABSTRACT

In this article we investigate the solution of the steady-state fractional diffusion equation on a bounded domain in $ eal^{1}$. From an analysis of the underlying model problem, we postulate that the fractional diffusion operator in the modeling equations is neither the Riemann-Liouville nor the Caputo fractional differential operators. We then find a closed form expression for the kernel of the fractional diffusion operator which, in most cases, determines the regularity of the solution. Next we establish that the Jacobi polynomials are pseudo eigenfunctions for the fractional diffusion operator. A spectral type approximation method for the solution of the steady-state fractional diffusion equation is then proposed and studied.

연구 동기 및 목표

• 표준 리만-리오빌리 및 카푸토 정의에 반박하며, 물리적으로 일관된 방식으로 분수형 확산 연산자를 재해석하는 것.
• 해결의 정(regularity)을 규정하는 분수형 확산 연산자의 핵을 규명하는 것.
• 최적의 수렴 속도를 보장하는 가중 잭비 함수를 사용한 스펙트럼 방법을 개발하는 것.
• 제안된 스펙트럼 근사가 가중 $ L^2 $ 공간에서 역행성과 사전 오차 추정치를 확립하는 것.

제안 방법

• 경계에서의 특이성을 피하기 위해, 분수도함수를 $ \mathbf{D}^\alpha u(x) = D \, \mathbf{D}^{-(2-\alpha)} \, D u(x) $ 로 새로운 정의한다.
• 조건을 만족시키기 위해 선택된 $ \beta $ 를 사용해 연산자 $ \mathcal{L}_r^\alpha $ 의 핵을 $ x^{\beta}(1-x)^{\alpha - \beta} $ 로 식별한다.
• 새로운 정의 하에서 잭비 다항식이 분수형 연산자의 의사 고유함수임을 보여준다.
• 가중치 $ \omega(x) = x^\beta(1-x)^{\alpha - \beta} $ 와 계수 $ c_j $ 가 내적을 통해 계산되는 스펙트럼 근사 $ u_N(x) = \omega(x) \sum_{j=0}^N c_j \mathcal{G}_j(x) $ 를 구성한다.
• $ L^2 $ 및 가중 $ L^2 $ 노름에서 사전 오차 추정치를 유도하며, 수렴 속도 $ \|u - u_N\| \leq C (N+2)^{-\alpha} \|f\|_{\omega^*} $ 를 보여준다.
• 스펙트럼 전개를 통해 $ L^2(0,1) $ 에서 $ \mathcal{L}_r^\alpha $ 의 역행성을 확립하며, 해는 $ L^2_{\omega^{-1}}(0,1) $ 에 속한다.

실험 결과

연구 질문

• RQ1 1차원 정적 확산의 맥락에서 분수도함수의 올바른 물리적 해석은 무엇이며, 리만-리오빌리 및 카푸토 정의와 어떻게 다를까?
• RQ2 분수도함수 정의의 선택이 양방향 분수형 확산 방정식의 해의 정(regularity)에 어떤 영향을 미치는가?
• RQ3 가중 잭비 함수를 기반으로 한 스펙트럼 방법이 분수형 확산 문제에서 최적의 수렴 속도를 달성할 수 있는가?
• RQ4 분수형 확산 연산자의 핵의 구조는 무엇이며, 이것이 해의 정(regularity)을 어떻게 결정하는가?
• RQ5 $ \mathcal{L}_r^\alpha $ 가 $ L^2(0,1) $ 에서 역행성이 있는가? 해가 속하는 함수 공간은 무엇인가?

주요 결과

• 분수도함수는 $ \mathbf{D}^\alpha u(x) = D \, \mathbf{D}^{-(2-\alpha)} \, D u(x) $ 로 재정의되며, 이는 경계에서의 특이성을 피하고 물리적 확산을 더 잘 모델링한다.
• $ \mathcal{L}_r^\alpha $ 의 핵은 $ x^\beta(1-x)^{\alpha - \beta} $ 이며, 그 구조는 해의 정(regularity)을 결정하며, 이는 매개변수 $ r $ 에 따라 달라진다.
• 새로운 정의 하에서 잭비 다항식은 분수형 연산자의 의사 고유함수이며, 이는 스펙트럼 근사 가능성을 보장한다.
• 제안된 스펙트럼 방법은 수렴 속도 $ \|u - u_N\| \leq C (N+2)^{-\alpha} \|f\|_{\omega^*} $ 를 달성하며, 이는 이론적 기대와 일치한다.
• $ \mathcal{L}_r^\alpha $ 는 $ L^2(0,1) $ 에서 역행성이 있으며, 해는 $ L^2_{\omega^{-1}}(0,1) $ 에 속하여 문제의 잘 정의됨을 보장한다.
• 수치 실험은 이론적 수렴 속도를 확인하며, 경계에서 갈라지는 갈라짐(Gibbs-type) 진동 현상이 나타나며, 이는 특이성 구조와 일치한다.