[논문 리뷰] A Finite Element Method for the Fractional Sturm-Liouville Problem

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

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

한 줄 요약

이 논문은 순서 α ∈ (1,2)인 Caputo 또는 Riemann-Liouville 도함수를 갖는 분수계 스투름-리우빌 문제를 해결하기 위한 새로운 유한요소법(FEM)을 제안한다. 새로운 변분형식을 바탕으로, 고유값 근사에서 2차 수렴성을 달성하고 엄밀한 오차 한계를 제공하여, 부드럽지 않은 잠재력 조건이 존재하더라도 다수의 고유값을 동시에 정확하게 계산할 수 있다.

ABSTRACT

In this work, we propose an efficient finite element method for solving fractional Sturm-Liouville problems involving either the Caputo or Riemann-Liouville derivative of order $α\in(1,2)$ on the unit interval $(0,1)$. It is based on novel variational formulations of the eigenvalue problem. Error estimates are provided for the finite element approximations of the eigenvalues. Numerical results are presented to illustrate the efficiency and accuracy of the method. The results indicate that the method can achieve a second-order convergence for both fractional derivatives, and can provide accurate approximations to multiple eigenvalues simultaneously.

연구 동기 및 목표

• Caputo 또는 Riemann-Liouville 도함수를 포함하는 분수계 스투름-리우빌 문제를 위한 강력하고 정확한 유한요소법을 개발하는 것.
• 다수의 고유값을 동시에 계산할 수 있는 효율적이고 안정적이며 수학적으로 타당한 수치적 방법의 부족을 해결하는 것.
• 일반적인 잠재함수 q(x)를 갖는 문제에 대해 고유값 근사에 대한 엄밀한 오차 추정을 제공하는 것.
• 분수계 스투름-리우빌 문제에서 고유값의 渐近적 행동, 분기 현상, 교차 성질 등의 복잡한 분석적 성질을 수치적으로 탐구할 수 있도록 하는 것.

제안 방법

• Caputo 및 Riemann-Liouville 도함수를 고려한 분수미분 연산자의 새로운 변분형식을 기반으로 하는 방법.
• 적절한 함수 공간을 사용: Riemann-Liouville의 경우 U = V = Ḣ^{α/2}(D)이며, Caputo의 경우 다른 공간을 사용하여 약한 형식의 안정성을 확보한다.
• 약한 형식은 모든 v ∈ V에 대해 a(u,v) = λ(u,v)로 유도되며, 이는 일반화된 고유값 문제로 이어진다.
• 조밀한 메esh 위에서 일차 다항식 근사를 사용하여 유한차원 부분공간 Uh ⊂ U 및 Vh ⊂ V를 구성한다.
• 결과적으로 얻어진 이산 일반화된 고유값 문제는 표준 수치선형대수 기법을 통해 해결된다.
• 정규성 이론과 보간 추정을 사용하여 오차 한계를 도출하였으며, q에 대한 부드러움 조건 하에 Riemann-Liouville의 경우 수렴률 r < α−1, Caputo의 경우 r < min(α+s,2)−1/2를 확보한다.

실험 결과

연구 질문

• RQ1 Caputo 및 Riemann-Liouville 분수계 스투름-리우빌 문제에 대해 둘 다 2차 수렴성을 달성하는 유한요소법을 개발할 수 있는가?
• RQ2 부드럽지 않은 잠재력 조건이 존재할 경우, 고유값 근사에 대한 최적의 수렴률은 무엇인가?
• RQ3 분기 상황에서 고유값과 고유함수는 어떻게 행동하는가? 그리고 수치적 방법이 이러한 복잡한 행동을 신뢰성 있게 포착할 수 있는가?
• RQ4 이 방법은 엄밀한 오차 한계를 보장하면서 다수의 고유값을 동시에 정확하게 근사할 수 있는가?
• RQ5 잠재항이 분수계 스투름-리우빌 문제에서 복소수 공액 고유값 쌍의 안정성과 구조에 어떤 영향을 미치는가?

주요 결과

• 수치 실험을 통해 Caputo 및 Riemann-Liouville 모두에서 유한요소법이 2차 수렴성을 확보함을 확인하였다.
• Riemann-Liouville의 경우, 고유값 근사 오차는 |λ − λh| ≤ Ch^{r}를 만족하며, r < α−1이다.
• Caputo의 경우, q ∈ Ḣ^s(D) ∩ L∞(D) 이며 s ∈ [0,1]이고 α+s > 3/2일 조건 하에 오차 한계는 |λ − λh| ≤ Ch^{r}이며, r < min(α+s,2)−1/2이다.
• 이 방법은 복소수 공액 쌍을 포함한 다수의 고유값을 동시에 성공적으로 계산하였으며, 임계 α 값 근처에서 분기 현상을 포착하였다.
• 수치 결과는 세 번째 고유함수가 두 번째 고유함수와 거의 구분되지 않음을 보여주며, 특정 α 값에서 고유함수의 구조에서 고유성의 상실 가능성을 시사한다.
• q2와 같은 잠재항의 존재는 분기 안정성을 해치며, 복소수 공액 고유값을 분리하고 고유함수의 형태를 변화시켜, 수치적 방법의 민감성을 부각시킨다.