[논문 리뷰] Numerical solution of fractional order diffusion problems with Neumann boundary conditions

Béla J. Szekeres, Ferenc Izsák|arXiv (Cornell University)|2014. 11. 06.

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

한 줄 요약

이 논문은 이동된 Grünwald–Letnikov 근사법을 사용하여 일차원 분수형 확산 방정식의 노이만 경계 조건을 해결하기 위한 유한 차분 방법을 제시한다. 초기 자료를 동질적 노이만 조건을 만족시키도록 연장하는 수학적 프레임워크를 수립하고, 그로 인한 문제의 잘 정의됨을 증명하며, 분수형 도함수 연산자의 폐쇄형 행렬 표현을 가진 암시적 오일러 스킴의 수렴성을 보여준다.

ABSTRACT

A finite difference numerical method is investigated for fractional order diffusion problems in one space dimension. For this, a mathematical model is developed to incorporate homogeneous Dirichlet and Neumann type boundary conditions. The models are based on an appropriate extension of the initial values. The well-posedness of the obtained initial value problems is proved and it is pointed out that the extensions are compatible with the above boundary conditions. Accordingly, a finite difference scheme is constructed for the Neumann problem using the shifted Grünwald--Letnikov approximation of the fractional order derivatives, which is based on infinite many basis points. The corresponding matrix is expressed in a closed form and the convergence of an appropriate implicit Euler scheme is proved.

연구 동기 및 목표

• 분수형 확산 문제에서 동질적 노이만 및 디리클레 경계 조건을 수학적으로 일관되게 모델링하기 위한 방법 개발.
• 분수형 도함수의 비국소성으로 인해 계산 도메인 외부의 값을 필요로 하는 문제 도전 과제 해결.
• 이동된 Grünwald–Letnikov 공식을 사용하여 노이만 문제에 대한 안정적이고 수렴하는 유한 차분 스킴 구축.
• 확장된 초기 자료와 경계 조건으로 형성된 초기값 문제의 잘 정의됨을 증명.
• 노이만 조건이 적용된 분수형 확산 방정식의 암시적 오일러 이산화에 대한 완전한 오차 분석 제공.

제안 방법

• 주기성과 대칭성을 이용해 초기 자료를 전체 실수선으로 연장함으로써 노이만 경계 조건을 만족시키는 수학적 모델을 구성.
• 리만–리우빌 적분과 도함수를 이용해 실수선 상에서 리에즈 분수형 도함수를 정의함으로써 비국소 연산자의 일관된 다루기 가능하게 함.
• 분수형 도함수를 이산화하기 위해 이동된 Grünwald–Letnikov 근사법을 적용하고, 이의 계수는 이항 전개로부터 유도.
• 결과로 얻어진 이산 연산자는 폐쇄형 행렬 형태로 표현되어 효율적인 수치적 구현 가능.
• 시간에 대해 암시적 오일러 스킴을 적용하고, 에너지 추정과 스펙트럼 분석을 통해 수렴성을 엄밀히 증명.
• 바르나우스 안정성 분석과 삼각함수 항등식을 사용한 증폭 인자 분석을 통해 안정성과 수렴성을 확립.

실험 결과

연구 질문

• RQ1 비국소 도함수를 가진 분수형 확산 방정식에 동질적 노이만 경계 조건을 일관되게 통합할 수 있는 방법은 무엇인가?
• RQ2 분수형 도함수의 맥락에서 노이만 경계 조건과 호환성이 보장되는 초기 자료의 연장 방식은 무엇인가?
• RQ3 이동된 Grünwald–Letnikov 공식에 기반한 유한 차분 스킴이 노이만 문제에 대해 수렴성을 증명할 수 있는가?
• RQ4 이산 분수형 도함수 연산자의 스펙트럼 구조는 무엇이며, 안정성에 어떤 영향을 미치는가?
• RQ5 노이만 조건이 적용된 분수형 확산 방정식에 대해 암시적 오일러 스킴의 수렴 차수는 어떻게 되는가?

주요 결과

• 실수선으로의 초기 자료 연장 방법은 동질적 노이만 경계 조건과 호환성을 보장하고 잘 정의된 초기값 문제의 정의를 가능하게 한다.
• 이산 분수형 도함수 연산자는 폐쇄형 행렬 표현으로 표현되어 있어 효율적인 수치적 해법을 가능하게 한다.
• 노이만 문제에 적용된 암시적 오일러 스킴은 시간에 대해 일관된 1차 수렴성을 가진다.
• 바르나우스 안정성 분 析을 통해 주어진 수식 구조 하에서 스킴이 일관된 안정성을 확보함을 확인.
• 이산 연산자의 고유값이 실수이면서 음이 아니라는 것이 입증되어 안정성 및 수렴성 결과를 지지.
• 스펙트럼 분석과 에너지 추정을 통해 공간 및 시간 모두에서 최적의 수렴 속도를 달성함.