[논문 리뷰] A Second-Order Crank-Nicolson Method for Subdiffusion

Bangti Jin, Buyang Li외 1명|arXiv (Cornell University)|2016. 07. 23.

Fractional Differential Equations Solutions참고 문헌 20인용 수 3

한 줄 요약

이 논문은 순서 $\alpha \in (0,1)$인 카푸토 분수계 도함수를 갖는 부분확산 방정식에 대해 두 번째 차수 정확도를 갖는 크랭크-니콜슨 유형의 시간 이산화 방법을 제안한다. 후방 오일러의 컨볼루션 쿼드러처와 $\theta = \alpha/2$인 $\theta$-방법을 조합하고 초기 조정을 적용함으로써, 공간에서 갈레르킨 유한요소 방법과 결합할 경우 부드럽고 부드럽지 않은 데이터 모두에서 시간 방향으로 두 번째 차수 수렴을 달성한다.

ABSTRACT

In this work, we analyze a Crank-Nicolson type time stepping scheme for the subdiffusion equation, which involves a Caputo fractional derivative of order $\alpha\in (0,1)$ in time. It combines the backward Euler convolution quadrature with a $ heta$-type method, with the parameter $ heta$ dependent on the fractional order $\alpha$ by $ heta=\alpha/2$, and it naturally generalizes the classical Crank-Nicolson method. We develop essential initial corrections for the fractional Crank-Nicolson scheme, and together with the Galerkin finite element method in space, obtain a fully discrete scheme. A complete error analysis of the fully discrete scheme is provided, and a second-order accuracy in time is established for both smooth and nonsmooth problem data. Extensive numerical experiments are provided to illustrate its accuracy, efficiency and robustness, and further, the comparative study indicates that the scheme is competitive with existing schemes.

연구 동기 및 목표

• 순서 $\alpha \in (0,1)$인 카푸토 분수도함수를 포함하는 부분확산 방정식에 대해 시간 기반으로 두 번째 차수 정확도를 갖는 이산화 방법을 개발하는 것.
• 기본적인 크랭크-니콜슨 방법을 분수계 부분확산으로 일반화하기 위해 $\theta = \alpha/2$인 $\theta$-의존 매개변수를 도입하는 것.
• 부드럽지 않은 초기 자료에서 수렴 차수의 손실 문제를 해결하기 위해 필수적인 초기 조정을 도입하는 것.
• 제안된 시간 스킴을 공간에서 갈레르킨 유한요소 방법과 결합하여 완전히 이산화된 스킴을 구성하는 것.
• 완전한 오차 분석을 제공하고 광범위한 수치 실험을 통해 스킴의 강건성과 효율성을 입증하는 것.

제안 방법

• 시간 이산화 방법은 $\theta = \alpha/2$인 $\theta$-방법과 후방 오일러의 컨볼루션 쿼드러처를 기반으로 하는 크랭크-니콜슨 유형의 스킴을 사용하며, 이는 고전적 방법을 분수계 수준으로 일반화함을 보장한다.
• 초기 조정은 부드럽지 않은 문제 자료에서 두 번째 차수 수렴을 복원하기 위해 유도되고 적용되며, 이에 의해 수렴 속도가 떨어지는 것을 방지한다.
• 공간 이산화 방법은 갈레르킨 유한요소 방법을 사용하여 최적의 안정성과 수렴 성질을 갖는 완전히 이산화된 스킴을 이끈다.
• 오차 분석은 $L^2$ 및 에너지 노름에서 수행되며, 부드럽고 부드럽지 않은 초기 자료 모두에서 시간 방향으로 두 번째 차수 수렴을 입증한다.
• 광범위한 수치 실험을 통해 스킴이 검증되었으며, 기존 방법과의 비교를 통해 정확성, 효율성, 강건성을 입증한다.

실험 결과

연구 질문

• RQ1 카푸토 분수도함수를 갖는 부분확산에 대해 시간 기반으로 두 번째 차수 정확도를 달성하는 크랭크-니콜슨 유형의 스킴을 구성할 수 있는가?
• RQ2 부드럽지 않은 초기 자료에서 분수계 부분확산 스킴의 두 번째 차수 수렴을 복원하기 위해 어떻게 초기 조정을 설계할 수 있는가?
• RQ3 제안된 스킴이 공간에서 갈레르킨 유한요소 방법과 조합되었을 때도 여전히 두 번째 차수 수렴을 유지하는가?
• RQ4 기존의 부분확산 수치 스킴과 비교해 볼 때 제안된 스킴의 정확성과 효율성은 어떠한가?
• RQ5 분수계 수준의 부분확산 환경에 적합하게 고전적 크랭크-니콜슨 방법을 일반화하기 위해 분수계 수준 $\alpha$에 따라 $\theta$ 매개변수를 조정할 수 있는가?

주요 결과

• 제안된 스킴은 부드럽고 부드럽지 않은 초기 자료 모두에서 시간 방향으로 두 번째 차수 수렴을 달성하며, 일반적으로 발생하는 분수계 부분확산에서의 수렴 차수 감소 문제를 해결한다.
• 초기 조정은 특히 부드럽지 않은 자료에서 두 번째 차수 정확도를 유지하는 데 필수적이며, 분석 과정에서 엄밀하게 유도되었다.
• 시간 스킴과 갈레르킨 유한요소 방법을 조합한 완전히 이산화된 스킴은 조건부가 아닌 안정성과 최적의 수렴 차수를 갖는다.
• 수치 실험은 이론적 오차 추정치를 확인하고 다양한 시험 케이스에서 스킴의 강건성을 입증한다.
• 비교 연구 결과 스킴은 정확성과 계산 효율성 측면에서 기존 방법과 경쟁 가능함을 보였다.
• $\theta = \alpha/2$의 매개변수 선택은 고전적 크랭크-니콜슨 방법을 분수계 부분확산 환경으로 효과적으로 일반화한다.