[논문 리뷰] An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data

Bangti Jin, Raytcho Lazarov외 1명|arXiv (Cornell University)|2015. 01. 01.

Fractional Differential Equations Solutions참고 문헌 24인용 수 57

한 줄 요약

이 논문은 순서 $\alpha \in (0,1)$인 캐프토 정수미분의 부분확산 방정식에 대해 L1 스킴을 재검토하며, 해가 시간에 대해 $C^2$ 정규성 없이도 성립하는 $O(\tau)$ 수렴 속도를 확립한다. 이는 부드럽지 않은 초기 자료에 대해서도 유효하다. 분석은 일반적인 섹터형 연산자와 시간-공간 분수확산 방정식으로 확장되며, 수치 실험을 통해 오차 추정치의 정확성과 강건성을 확인한다.

ABSTRACT

The subdiffusion equation with a Caputo fractional derivative of order $α\in(0,1)$ in time arises in a wide variety of practical applications, and it is often adopted to model anomalous subdiffusion processes in heterogeneous media. The L1 scheme is one of the most popular and successful numerical methods for discretizing the Caputo fractional derivative in time. The scheme was analyzed earlier independently by Lin and Xu (2007) and Sun and Wu (2006), and an $O(τ^{2-α})$ convergence rate was established, under the assumption that the solution is twice continuously differentiable in time. However, in view of the smoothing property of the subdiffusion equation, this regularity condition is restrictive, since it does not hold even for the homogeneous problem with a smooth initial data. In this work, we revisit the error analysis of the scheme, and establish an $O(τ)$ convergence rate for both smooth and nonsmooth initial data. The analysis is valid for more general sectorial operators. In particular, the L1 scheme is applied to one-dimensional space-time fractional diffusion equations, which involves also a Riemann-Liouville derivative of order $β\in(3/2,2)$ in space, and error estimates are provided for the fully discrete scheme. Numerical experiments are provided to verify the sharpness of the error estimates, and robustness of the scheme with respect to data regularity.

연구 동기 및 목표

• 기존 L1 스킴 수렴 분석에서 시간에 대한 $C^2$ 정규성 가정에 기인한 한계를 해결하기 위해, 이는 부드럽지 않은 초기 자료에서 실패한다.
• 초기 자료에 대한 최소한의 정규성 조건 하에서 L1 스킴의 엄밀한 $O(\tau)$ 수렴 속도를 확립하기 위해.
• 공간-시간 분수확산 방정식에서 발생하는 연산자들을 포함한 일반적인 섹터형 연산자로 오차 분석을 확장하기 위해.
• 이론적 결과를 수치 실험을 통해 검증하여 오차 추정치의 정확성과 강건성을 입증하기 위해.

제안 방법

• 시간 도함수의 조각별 선형 근사에 기반한 캐프토 분수미분에 대한 L1 스킴 재구성.
• 라플라스 변환 기법을 활용한 생성 함수 접근을 통해 局부 절단 오차를 분석하기 위해.
• 섹터형 연산자의 성질과 분수확산 반그룹의 스무딩 효과를 이용해 오차 한계 유도하기.
• 리만-리우빌 공간 도함수를 가진 공간-시간 분수확산 방정식에 대해 완전 이산 유한요소 또는 유한차분 스킴에 분석 적용하기.
• 시간 및 공간 오차를 제어하기 위해 $L^2$-노름에서 에너지 방법과 안정성 추정치 사용하기.
• 균일 및 비균일 시간 메esh에서의 수치 실험을 통해 이론적 수렴 속도 검증하기.

실험 결과

연구 질문

• RQ1 기존의 $C^2$ 정규성 가정이 실패하는 부드럽지 않은 초기 자료를 가진 부분확산에서 L1 스킴이 일阶 수렴을 달성할 수 있는가?
• RQ2 해가 $t=0$에서 특이성을 보일 경우 L1 스킴의 최적 수렴 속도는 무엇인가?
• RQ3 초기 자료의 정규성에 따라 오차 행동은 어떻게 달라지며, 특히 $t \to 0$ 근처에서 어떻게 변화하는가?
• RQ4 오차 분석은 공간-시간 분수확산 방정식에서 발생하는 연산자들을 포함한 일반적인 섹터형 연산자로 확장 가능한가?
• RQ5 공간-시간 분수확산 모델의 $\alpha$ 및 $\beta$ 값에 관계없이 $O(\tau)$ 수렴 속도가 강건한가?

주요 결과

• 해가 시간에 대해 $C^2$ 정규성이 없더라도, 부드럽고 부드럽지 않은 초기 자료에 대해 L1 스킴은 $L^2$-노름에서 $O(\tau)$ 수렴 속도를 달성한다.
• $\alpha \in (0,1)$에 대해 수렴 속도는 강건하며, 수치 실험을 통해 $t=0.1$, $t=0.01$, $t=0.001$에서 일阶 수렴이 확인된다.
• 공간-시간 분수확산 방정식에서 $\beta \in (3/2,2)$일 경우, 적절한 정규성 조건 하에 완전 이산 스킴은 $O(\tau + h^2)$ 오차 추정치를 달성한다.
• 부드럽지 않은 초기 자료의 경우 오차는 $t \to 0$에서 $t^{\alpha}$와 유사하게 행동하며, 이는 정리 4.1의 이론적 추정치의 정확성을 확인한다.
• 수치 결과는 $\alpha$, $\beta$, 초기 자료 유형에 관계없이 수렴 속도가 약 1.0을 유지함을 보여주며, 이는 강건성을 시사한다.
• 이론적 분석은 일반적인 섹터형 연산자에 대해 유효하며, 이는 L1 스킴의 적용 가능 범위를 표준 타원형 연산자 초과로 확장한다.