[논문 리뷰] Finite difference schemes for multi-term time-fractional mixed diffusion-wave equations

Zhaopeng Hao, Guang Lin|arXiv (Cornell University)|2016. 07. 24.

Fractional Differential Equations Solutions참고 문헌 12인용 수 7

한 줄 요약

이 논문은 (1,2) 범위의 순서를 가진 Caputo 도함수를 갖는 다항 시간 분수형 혼합 확산-파동 방정식에 대해 1차 정확도를 갖는 L2 유형의 유한차분 스킴을 제안한다. 수정된 L2 근사와 고밀도 차분 스킴을 조합하고 에너지 방법을 적용함으로써, 시간에 대해 1차 정확도, 공간에 대해 4차 정확도를 갖는 조건부 수렴을 증명하였다. 수치 예제와 선형로그 복잡도를 갖는 빠른 해법을 통해 검증되었다.

ABSTRACT

The multi-term time-fractional mixed diffusion-wave equations (TFMDWEs) are considered and the numerical method with its error analysis is presented in this paper. First, a $L2$ approximation is proved with first order accuracy to the Caputo fractional derivative of order $β\in (1,2).$ Then the approximation is applied to solve a one-dimensional TFMDWE and an implicit, compact difference scheme is constructed. Next, a rigorous error analysis of the proposed scheme is carried out by employing the energy method, and it is proved to be convergent with first order accuracy in time and fourth order in space, respectively. In addition, some results for the distributed order and two-dimensional extensions are also reported in this work. Subsequently, a practical fast solver with linearithmic complexity is provided with partial diagonalization technique. Finally, several numerical examples are given to demonstrate the accuracy and efficiency of proposed schemes.

연구 동기 및 목표

• 이상적 확산을 모델링하는 데 중요한데도 불구하고 연구가 부족한 (1,2) 범위의 분수 차수를 갖는 다항 시간 분수형 혼합 확산-파동 방정식(TFMDWEs)에 대해 강력한 수치적 방법을 개발하는 것.
• 분수 도함수의 비국소성과 기억 요구 사항에 의해 발생하는 과제를 해결하기 위해, 엄밀한 오차 분석을 갖는 안정적이고 정확한 유한차분 스킴을 구축하는 것.
• 해법을 분포형 차수 방정식과 2차원 문제로 확장하여 광범위한 적용 가능성을 확보하는 것.
• 테일로프 구조를 가진 선형 시스템을 처리하기 위해 부분 대각화 기반으로 선형로그 복잡도를 갖는 빠른 해법을 설계하는 것.
• 에너지 방법을 사용한 완전한 오차 분석을 제공하여 조건부 수렴과 안정성을 확립하는 것.

제안 방법

• 순서 β ∈ (1,2)인 Caputo 분수 도함수에 대해 수정된 L2 근사를 유도하여 시간에 대해 1차 정확도를 달성한다.
• 수정된 L2 근사와 낮은 차수 도함수를 위한 L1 근사를 조합하여 완전히 암시적인 고밀도 유한차분 스킴을 형성한다.
• 이산 시스템에 에너지 방법을 적용하고, 새로운 이산 분수 Poincaré-Friedrichs-Sobolev 부등식을 사용하여 안정성과 수렴성을 확립한다.
• 선형 시스템의 토플리츠 구조를 다루기 위해 부분 대각화 기반 빠른 해법을 개발하여 O(N log N) 복잡도를 달성한다.
• 스킴을 순서가 (0,2) 범위인 분포형 차수 방정식과 2차원 문제로 확장하여 고차 정밀도를 유지한다.
• 이론적 수렴 속도를 검증하고 1D 및 2D 문제에서의 효율성을 입증하기 위해 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 순서가 (1,2) 범위인 Caputo 분수 도함수에 대해 1차 정확도를 갖는 L2 유형 근사를 엄밀하게 유도할 수 있는가?
• RQ2 다항 시간 분수형 혼합 확산-파동 방정식에 대해 고밀도 유한차분 스킴을 구성하고 조건부 수렴성을 증명할 수 있는가?
• RQ3 스킴의 수렴 속도는 시간과 공간에서 어떻게 되며, 에너지 방법을 사용하여 이를 증명할 수 있는가?
• RQ4 결과로 생기는 선형 시스템에 대해 선형로그 복잡도를 갖는 빠른 해법을 설계할 수 있는가?
• RQ5 스킴은 분포형 차수 방정식과 2차원 문제에서 어떻게 성능을 발휘하는가?

주요 결과

• 수정된 L2 근사가 순서 β ∈ (1,2)인 Caputo 도함수에 대해 1차 정확도를 달성하며, 엄밀한 오차 한계가 증명되었다.
• 에너지 방법을 사용하여 고밀도 유한차분 스킴이 시간에 대해 1차 정확도, 공간에 대해 4차 정확도를 갖는 조건부 수렴성을 증명하였다.
• 이산 분수 Poincaré-Friedrichs-Sobolev 부등식이 확립되었으며, 이는 수렴 증명에서 핵심적인 역할을 하였다.
• 수치 결과는 1D 및 2D 문제에서 이론적 수렴 속도가 시간과 공간 모두에서 확인됨을 보여주었다.
• O(N log N) 복잡도를 갖는 빠른 해법이 성공적으로 구현되어 표준 해법 대비 계산 비용을 크게 감소시켰다.
• 스킴은 순서가 (0,2) 범위인 분포형 차수 방정식에서도 강력한 성능을 발휘하여 광범위한 적용 가능성을 확인하였다.