[논문 리뷰] Sharp $H^1$-norm error estimates of two time-stepping schemes for reaction-subdiffusion problems

Jincheng Ren, Hong-lin Liao외 2명|arXiv (Cornell University)|2018. 11. 20.

Fractional Differential Equations Solutions참고 문헌 9인용 수 21

한 줄 요약

이 논문은 캐프토 유형 분수 미분을 갖는 반응-하위확산 문제에 대해 비균일 시간 스텝 구조를 갖는 두 가지 스킴—L1 및 분수 차분 크랭크-니콜슨—에 대해 낼맞는 H¹-노름 오차 추정을 수립한다. 개선된 이산 그로워렐 부등식을 도입하고, 시간-공간 오차 분리 기법과 전역 일致성 분석을 결합함으로써, 초기 특이성에도 불구하고 최적 수렴 속도를 회복한다. 이로써 최적의 메esh 그레ading 파rameter를 갖는 비균일 메쉬에서 두 스킴에 대해 H¹-노름 오차 한계를 𝒪(τ^{2−α})의 순서로 확보한다.

ABSTRACT

Due to the intrinsically initial singularity of solution and the discrete convolution form in numerical Caputo derivatives, the traditional $H^1$-norm analysis (corresponding to the case for a classical diffusion equation) to the time approximations of a fractional subdiffusion problem always leads to suboptimal error estimates (a loss of time accuracy). To recover the theoretical accuracy in time, we propose an improved discrete Grönwall inequality and apply it to the well-known L1 formula and a fractional Crank-Nicolson scheme. With the help of a time-space error-splitting technique and the global consistency analysis, sharp $H^1$-norm error estimates of the two nonuniform approaches are established for a reaction-subdiffusion problems. Numerical experiments are included to confirm the sharpness of our analysis.

연구 동기 및 목표

• 캐프토 유형 분수 미분의 초기 특이성과 컨볼루션構조로 인해 기존의 H¹-노름 오차 추정이 비최적화되는 분수 하위확산 문제의 문제를 해결하기 위해.
• 기본 분석에서 일반적으로 상실되는 비균일 시간 이산화에서 H¹-노름의 최적 시간 수렴 속도를 회복하기 위해.
• 리아오 등이 제안한 이론적 프레임워크를 확장하여 일반적인 비균일 메쉬에서 L1 및 분수 크랭k-니콜슨 스킴에 대해 날맞는 H¹-노름 오차 추정을 달성하기 위해.
• 다양한 메쉬 그레딩 파라미터와 분수 차수를 갖는 수치 실험을 통해 이론적 추정의 날맞음 정도를 검증하기 위해.

제안 방법

• 비균일 시간 이산화에서 캐프토 도함수의 초기 특이성과 컨볼루션 특성을 다룰 수 있도록 특화된 개선된 이산 그로워렐 부등식을 제안한다.
• 시간-공간 오차 분리 기법을 적용하여 시간 오차와 공간 오차를 분리함으로써 각 성분을 독립적으로 분석할 수 있도록 한다.
• 전역 일치성 분석을 통해 전체 시간 간격 동안 국소 절단 오차를 제어하며, 특히 t=0 근처에서의 제어를 강화한다.
• 초기 특이성을 해결하기 위해 파라미터 γ를 갖는 계단형 시간 메쉬를 사용하며, τk ∝ tk−1+γ−1로 정의된다.
• 개선된 이산 그로워렐 부등식과 안정성 및 일치성 논증을 결합하여 L1 및 분수 크랭크-니콜슨 스킴에 대해 날맞는 H¹-노름 오차 추정을 수립한다.
• 유한 차분 공간 이산화를 사용하고 메쉬 그레딩 파라미터를 다양하게 설정한 수치 실험을 통해 이론적 결과를 검증한다.

실험 결과

연구 질문

• RQ1 캐프토 도함수를 갖는 반응-하위확산 문제의 비균일 시간 이산화에서 초기 특이성과 컨볼루션 특성에도 불구하고 최적의 H¹-노름 오차 추정을 회복할 수 있는가?
• RQ2 H¹-노름 분석이 t=0 근처에서 해의 비연속성으로 인해 최적 수렴 속도를 달성하지 못하는가? 만약 그렇다면 이를 어떻게 수정할 수 있는가?
• RQ3 개선된 이산 그로워렐 부등식이 L1 및 분수 크랭크-니콜슨 스킴에 대해 H¹-노름에서 최적 시간 수렴을 복원하는 데 효과적으로 활용될 수 있는가?
• RQ4 이러한 스킴에서 H¹-노름에서 이론적 수렴 속도 𝒪(τ^{2−α})를 달성하기 위한 최적의 그레딩 파라미터 γ는 무엇인가?
• RQ5 다양한 메쉬 파rameter와 분수 차수를 갖는 수치 실험을 통해 이론적 H¹-노름 오차 추정이 날맞는가?

주요 결과

• 논문은 초기 특이성에도 불구하고 비균일 시간 메쉬에서 L1 및 분수 크랭크-니콜슨 스킴에 대해 날맞는 H¹-노름 오차 추정 𝒪(τ^{2−α})을 수립하여 최적 수렴 속도를 회복한다.
• 개선된 이산 그로워렐 부등식은 H¹-노름 분석에서의 정확도 손실를 효과적으로 완화하여, 기존 방법이 실패하는 영역에서도 최적 오차 한계를 달성한다.
• 수치 실험 결과, 메쉬 그레딩 파라미터 γ ≥ γ_opt일 경우 수렴 속도가 𝒪(τ^{2−α})에 도달함을 확인하여 이론적 추정의 날맞음 정도를 검증한다.
• L1 스킴의 경우 γσ ≥ 2−α일 때 최적 수렴이 달성되며, 분수 크랭크-니콜슨 스킴의 경우 γσ ≥ 2일 때 최적 수렴이 달성되며, 이에 따라 γ_opt가 결정된다.
• 시간-공간 오차 분리 기법은 시간 오차와 공간 오차를 효과적으로 분리하여, 시간 이산화에 의한 H¹-노름 오차 기여도를 정밀하게 제어하고 분석할 수 있게 한다.
• 수치 결과는 비균일 메쉬가 균일 메쉬에 비해 수렴 속도를 크게 향상시키며, γ ≥ γ_opt일 경우 이론적 예측을 확인함을 보여준다.