[논문 리뷰] Two families of novel second-order fractional numerical formulas and their applications to fractional differential equations

Baoli Yin, Yang Liu외 2명|arXiv (Cornell University)|2019. 06. 04.

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

한 줄 요약

이 논문은 리만-리우빌 분수적 적분 연산자에 대한 새로운 생성 함수를 구성하여, 두 가지 새로운 2차 분수 수치 방법—분수 BT-θ 및 BN-θ 방법—을 제안한다. 이 방법들은 초기 특이성에 대한 보정 항을 통해 최적의 정확도를 확보하며, 2차 수렴을 달성한다. 안정성 분석을 통해 A(ϑ)-안정성이 확인되었으며, 수치 실험을 통해 검증되었다.

ABSTRACT

In this article, we introduce two families of novel fractional $θ$-methods by constructing some new generating functions to discretize the Riemann-Liouville fractional calculus operator $\mathit{I}^α$ with a second order convergence rate. A new fractional BT-$θ$ method connects the fractional BDF2 (when $θ=0$) with fractional trapezoidal rule (when $θ=1/2$), and another novel fractional BN-$θ$ method joins the fractional BDF2 (when $θ=0$) with the second order fractional Newton-Gregory formula (when $θ=1/2$). To deal with the initial singularity, correction terms are added to achieve an optimal convergence order. In addition, stability regions of different $θ$-methods when applied to the Abel equations of the second kind are depicted, which demonstrate the fact that the fractional $θ$-methods are A($\vartheta$)-stable. Finally, numerical experiments are implemented to verify our theoretical result on the convergence analysis.

연구 동기 및 목표

• 리만-리우빌 분수적 적분 연산자에 대해 향상된 정확도와 안정성을 갖춘 새로운 2차 수치 공식을 개발한다.
• θ에 의해 매개변수화된 연속적인 방법의 가닥을 제안하여 기존의 분수적 방법들—FBDF2, FTR, GNGF2—사이의 격차를 메운다.
• 분수적 미분 방정식에서의 초기 특이성 문제를 효과적인 보정 항을 통해 해결하여 최적의 수렴 차수를 유지한다.
• 제안된 θ-방법의 안정성 특성을 분석하며, 특히 야벨형 제2종 방정식에 적용했을 때의 A(ϑ)-안정성을 다룬다.

제안 방법

• 새로운 생성 함수 ω(ξ) = p₁(ξ)^α p₂(ξ)^β 형태를 기반으로 한 분수 BT-θ 및 BN-θ 방법이라는 두 가지 새로운 분수 수치 방법의 가닥을 제안한다.
• BT-θ 방법의 경우 FBDF2 (θ=0)와 FTR (θ=1/2) 사이를 연결하는 생성 함수를 구성하며, BN-θ 방법의 경우 FBDF2 (θ=0)와 GNGF2 (θ=1/2) 사이를 연결한다.
• 생성 함수의 도함수 계수로부터 유도된 재귀 공식을 통해 컨볼루션 가중치 ωₖ를 계산하는 공식을 유도하여 O(N) 복잡도로 효율적인 계산을 가능하게 한다.
• 초기 특이성을 다루기 위해 수치적 방법에 보정 항을 적용하여, 약한 특이성을 갖는 해가 존재하는 상황에서도 2차 수렴을 복원한다.
• 야벨 방정식 제2종에 적용한 θ-방법의 안정성 영역을 유도하고 도식화하여 A(ϑ)-안정성을 확인한다.
• 이론적 수렴 및 안정성 결과를 분수적 미분 방정식에 대한 수치 실험을 통해 검증한다. 이 방정식들은 알려진 해를 갖는다.

실험 결과

연구 질문

• RQ1 FBDF2, FTR, GNGF2와 같은 기존 방법들을 부드럽게 연결하는 통합된 2차 분수 수치 방법의 가닥을 구성할 수 있는가?
• RQ2 리만-리우빌 분수적 적분 연산자에 대해 2차 수렴을 보장하기 위해 생성 함수를 어떻게 설계할 수 있는가?
• RQ3 초기 특이성이 있는 분수적 미분 방정식을 풀 때 최적의 수렴 차수를 복원하기 위해 필요한 보정 항은 무엇인가?
• RQ4 제안된 분수 θ-방법은 A(ϑ)-안정성이 있는가? 그리고 다양한 θ 값에 대해 안정성 영역은 어떻게 행동하는가?
• RQ5 복잡한 적분을 직접 평가하지 않고도 새로운 분수 방법의 컨볼루션 가중치를 계산하기 위한 효율적인 재귀 알고리즘을 유도할 수 있는가?

주요 결과

• 제안된 분수 BT-θ 및 BN-θ 방법은 수렴 분석을 통해 리만-리우빌 분수적 적분 연산자에 대해 2차 수렴을 달성한다.
• 보정 항이 성공적으로 방법에 통합되어, 초기 특이성이 존재하는 상황에서도 최적의 수렴 차수를 유지한다.
• 야벨 방정식 제2종에 적용한 θ-방법의 안정성 영역은 A(ϑ)-안정성을 보여주며, 굳은 분수 문제에 대한 강건성을 확인한다.
• 컨볼루션 가중치 ωₖ를 계산하기 위한 재귀 공식이 도출되었으며, O(N) 계산 복잡도를 갖는 것으로 확인되어 효율적인 구현이 가능하다.
• 수치 실험을 통해 이론적 수렴 속도와 안정성 특성이 확인되었으며, 제안된 방법의 효과성을 검증한다.
• 생성 함수 접근법을 통해 다항식 기반 함수와 조정 가능한 매개변수(θ)를 조합함으로써 새로운 분수적 방법을 체계적으로 구성할 수 있으며, 이는 탄력성과 조정 가능성의 향상을 가능하게 한다.