[논문 리뷰] Three-Point Compact Approximation for the Caputo Fractional Derivative

Yuri Dimitrov|arXiv (Cornell University)|2015. 10. 06.

Fractional Differential Equations Solutions참고 문헌 18인용 수 6

한 줄 요약

이 논문은 공간과 시간에서 세 번째 순서의 정확도를 가지며, 리만 제타 함수 보정을 사용하여 L1 근사치를 수정함으로써 $O(h^{3-\beta})$ 수렴성을 달성하는 캡투 분수도수의 삼점(compact) 유한차분 스킴을 제안한다. 이 방법은 분수차 하위확산 및 이완 방정식에 대해 안정성과 정확도를 향상시키면서도 다양한 $\alpha$ 값에서 단순한 계수 구조를 유지한다.

ABSTRACT

In this paper we derive the fourth-order asymptotic expansions of the trapezoidal approximation for the fractional integral and the $L1$ approximation for the Caputo derivative. We use the expansion of the $L1$ approximation to obtain the three point compact approximation for the Caputo derivative \begin{equation*} \dfrac{1}{Γ(2-α)h^α}\sum_{k=0}^{n} δ_k^{(α)} y_{n-k}=\dfrac{13}{12}y^{(α)}_n-\dfrac{1}{6}y^{(α)}_{n-1}+\dfrac{1}{12}y^{(α)}_{n-2}+O\left(h^{3-α} ight), \end{equation*} with weights $δ_0^{(α)}=1-ζ(α-1),\; δ_n^{(α)}=(n-1)^{1 -α}-n^{1-α},$ $$ δ_1^{(α)}=2^{1-α}-2+2ζ(α-1),\; δ_2^{(α)}=1-2^{2-α}+3^{1-α}-ζ(α-1),$$ $$δ_k^{(α)}=(k-1)^{1-α}-2k^{1-a}+(k+1)^{1-α},\quad (k=3\cdots,n-1),$$ where $y$ is a differentiable function which satisfies $y'(0)=0$. The numerical solutions of the fractional relaxation and the time-fractional subdiffusion equations are discussed.

연구 동기 및 목표

• 표준 L1 근사치의 $O(h^{2-\alpha})$ 정확도를 초월하는 고차수 compact 유한차분 스킴을 개발하는 것.
• 다양한 $\alpha$ 값에서 일관된 계수를 유지하면서도 고차수 수렴성을 확보하는 삼점 compact 근사치를 유도하는 것.
• 스킴을 분수차 이완 및 시간 분수 하위확산 방정식에 적용하여 개선된 수치 성능를 입증하는 것.
• L1 가중치와 제타 함수 보정을 사용하여 compact 근사치를 구성하는 프레임워크를 수립하는 것.

제안 방법

• Trapzoidal 규칙과 Caputo 도수의 L1 근사치에 대한 네 번째 순서 점근 전개를 유도한다.
• 첫 번째 세 개의 가중치를 $\zeta(\alpha-1)$를 사용하여 조정함으로써 L1 근사치를 수정하고, 삼점 스티프를 얻는다.
• Compact 스킴을 구성한다: $\frac{1}{\Gamma(2-\alpha)h^{\alpha}}\sum_{k=0}^{n}\delta_k^{(\alpha)}y_{n-k} = \frac{13}{12}y^{(\alpha)}_n - \frac{1}{6}y^{(\alpha)}_{n-1} + \frac{1}{12}y^{(\alpha)}_{n-2} + O(h^{3-\alpha})$.
• 리만 제타 함수와 멱법칙 차분을 사용하여 계수 $\delta_k^{(\alpha)}$를 유도하며, 모든 $\alpha \in (0,1)$에 대해 일관성을 확보한다.
• 수치적으로 분수차 이완 및 시간 분수 하위확산 방정식을 해결하기 위해 스킴을 적용한다.
• 제조된 해를 사용하여 방법을 검증하며, 이는 $u(x,t) = (1+2x^2+3x^3)t^{3-\alpha}$ 및 $u(x,t) = \sin x E_\alpha(-t^\alpha)$와 같은 알려진 해석적 형태를 포함한다.

실험 결과

연구 질문

• RQ1 L1 가중치를 수정하여 $O(h^{3-\alpha})$ 정확도를 가지는 삼점 compact 유한차분 스킴을 캡투 도수에 대해 구성할 수 있는가?
• RQ2 첫 번째 세 개의 계수에 $\zeta(\alpha-1)$를 포함시키는 것이 L1 근사치의 수렴 순서를 어떻게 향상시키는가?
• RQ3 기본 L1 및 수정된 L1 방법과 비교하여 이 스킴이 분수차 하위확산 및 이완 방정식에서 성능이 어떻게 되는가?
• RQ4 다양한 $\alpha$ 값에서 일관된 계수 구조를 유지하면서도 고차수 정확도를 확보할 수 있는가?

주요 결과

• 제안된 compact 스킴은 $O(h^{3-\alpha})$ 수렴성을 달성하여 표준 L1 방법의 $O(h^{2-\alpha})$ 정확도보다 크게 향상된다.
• 분수차 이완 방정식에 대한 수치 결과는 $\alpha = 0.4$일 때 공간 수렴 순서가 약 2.0, 시간 수렴 순서가 약 2.17임을 보여준다.
• $\alpha = 0.5$인 하위확산 방정식의 경우, 공간 수렴 순서는 약 2.0에 가까워지고 시간 수렴 순서는 약 2.5에 이르며, 높은 정확도를 나타낸다.
• 해의 조건 $y'(0) = 0$이 만족될 경우, 다양한 $\alpha$ 값에서 일관된 계수 구조를 유지한다.
• 모든 시험 케이스에서 표준 L1 및 수정된 L1 방법보다 공간 및 시간 수렴 속도가 뛰어나 성능이 뛰어나다.
• 이 스킴은 L1 가중치와 $\zeta(\alpha-1)$ 보정을 사용하여 $O(h^{3-\alpha})$ 정확도를 달성하는 첫 번째 삼점 compact 유한차분 방법이다.