[논문 리뷰] Variational formulation of problems involving fractional order differential operators

Bangti Jin, Raytcho Lazarov외 1명|arXiv (Cornell University)|2013. 07. 17.

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

한 줄 요약

이 논문은 순서 $\alpha \in (1,2)$인 캡투 및 리만-리우빌 분수도 도함수를 포함하는 경계값 문제에 대한 변분 형식을 개발하며, 최적의 수렴 속도를 보장하는 강도, 안정성 및 이동 정리들을 확립한다. 주요 기여는 유한요소 방법에서 최적의 수렴 속도를 보장하는 정밀한 유한요소 오차 분석으로, 이는 오직 우변의 정(regularity)에 의존한다. 수치 실험을 통해 검증되었다.

ABSTRACT

In this work, we consider boundary value problems involving Caputo and Riemann-Liouville fractional derivatives of order $α\in(1,2)$ on the unit interval $(0,1)$. These fractional derivatives lead to non-symmetric boundary value problems, which are investigated from a variational point of view. The variational problem for the Riemann-Liouville case is coercive on the space $H_0^{α/2}(0,1)$ but the solutions are less regular, whereas that for the Caputo case involves different test and trial spaces. The numerical analysis of these problems requires the so-called shift theorems which show that the solutions of the variational problem are more regular. The regularity pickup enables one to establish convergence rates of the finite element approximations. Finally, numerical results are presented to illustrate the error estimates.

연구 동기 및 목표

• 순서 $\alpha \in (1,2)$인 캡투 및 리만-리우빌 도함수를 포함하는 분수계수 경계값 문제에 대한 변분 형식을 개발한다.
• 해의 정(regularity)을 향상시키는 보장이 되는 변분 안정성 및 이동 정리들을 확립한다.
• 유한요소 근사의 수렴 속도를 우변 $f$의 부드러움에만 의존하도록 유도한다.
• 이론적 오차 추정치를 수치 실험을 통해 검증한다.

제안 방법

• 특히 리만-리우빌의 경우 $H_0^{\alpha/2}(0,1)$, 캡투의 경우 별개의 시험 및 기저 공간을 사용하여 분수계수 소볼레프 공간을 활용한 약한 문제를 설정한다.
• 해의 잘 정의됨을 보장하기 위해 $H_0^{\alpha/2}(0,1)$에서 리만-리우빌 변분 형식의 강도를 증명한다.
• 해가 우변 항 $f$의 정(regularity) 비례로 정(regularity)을 획득함을 보여주는 이동 정리들을 확립하며, 이는 최적의 유한요소 수렴 속도를 가능하게 한다.
• 균일 메esh에서 다항식 근사와 함께 갈레르킨 유한요소 방법을 적용한다.
• $L^2$ 및 $H^{\alpha/2}$-노름에서 오차 추정치를 유도하며, 이는 $f$의 정(regularity)에 따라 달라진다.
• 다양한 $\alpha$와 $f$의 정(regularity)에서 이론적 수렴 속도를 검증하기 위해 특이한 우변 항을 가진 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 순서 $\alpha \in (1,2)$인 캡투 및 리만-리우빌 분수도 도함수에 대해 일관된 변분 형식을 어떻게 유도할 수 있는가?
• RQ2 이러한 변분 문제의 해에 대해 성립하는 정(regularity) 성질(이동 정리)은 무엇이며, 이는 우변 항 $f$에 어떻게 의존하는가?
• RQ3 이러한 분수계수 문제의 갈레르킨 유한요소 근사에 대해 엄밀하게 확립할 수 있는 수렴 속도는 무엇인가?
• RQ4 $L^2$-노름 오차 추정치가 실질 수렴 속도보다 낮은 이유는 무엇이며, 이 격차는 어떻게 설명되는가?
• RQ5 우변 항 $f$가 특이하고 $\alpha \to 1^+$일 때 해는 원점 근처에서 어떻게 행동하는가? 그리고 이를 해결하기 위한 대안은 무엇인가?

주요 결과

• 리만-리우빌 변분 형식은 $H_0^{\alpha/2}(0,1)$에서 강도를 보이며 잘 정의됨을 보장하지만, 해는 형태 $x^{\alpha-1}$의 특이성을 보이며, 이는 균일 메쉬에서 진동을 유도한다.
• 캡투 형식은 별개의 시험 및 기저 공간이 필요하며, 적절한 함수 공간 설정 하에서 변분 형식이 안정함을 보였다.
• 이동 정리들은 해가 우변 항 $f$의 정(regularity)에 비례하여 정(regularity)을 획득함을 보여주며, 이는 최적의 유한요소 수렴 속도를 가능하게 한다.
• 수치 결과는 이론적 예측과 일치하는 $H^{\alpha/2}$-노름 수렴 속도를 확인하였으며, $\alpha \in (1,2)$ 범위에서 관측된 속도는 약 1.0–1.5 수준이었다.
• $L^2$-노름에서는 이론적 추정치보다 약 반계수 정도 더 높은 실질 수렴 속도를 보였으며, 이는 향후 연구가 필요한 격차를 시사한다.
• 특이한 우변 항 $f(x) = x^{-1/4}$의 경우, $\alpha \to 1^+$에 수렴함에 따라 수렴 속도가 느려지며, 특히 캡투의 경우 더욱 두드러지지만, 여전히 이론적 기대와 일치하는 수렴 속도를 보였다.