[논문 리뷰] ON TWO SCHEMES FOR FRACTIONAL DIFFUSION AND DIFFUSION-WAVE EQUATIONS

Bangti Jin, Raytcho Lazarov외 1명|arXiv (Cornell University)|2014. 04. 15.

Fractional Differential Equations Solutions참고 문헌 39인용 수 13

한 줄 요약

이 논문은 시간 분수차분 연속형 미분방정식과 시간 분수차분 연속형 파동방정식을 해결하기 위해 갈레르킨 유한요소법과 후진 오일러 및 2차 후진 차분법을 조합한 완전 이산화된 두 가지 스킴—시간에 대해 1차 및 2차 정밀도를 갖는다—을 제안한다. 이 스킴들은 부드럽고 부드럽지 않은 초기 자료 모두에서 최적의 오차 추정치를 달성하며, 수치 실험을 통해 뛰어난 정확성과 기존 방법들과의 경쟁력을 입증한다.

ABSTRACT

We consider the initial/boundary value problem for the fractional diu- sion and diusion-wave equations involving a Caputo fractional derivative in time. We develop two simple fully discrete schemes based on the Galerkin nite element method and the implicit backward Euler method/second-order backward dierence method, and establish error estimates optimal with respect to the regularity of the initial data. These two schemes are rst and second-order accurate in time for both smooth and nonsmooth initial data. Extensive numerical experiments for one and two-dimension problems conrm the convergence analysis. A detailed comparison with several existing time stepping schemes is also performed. The numerical re- sults indicate that the proposed fully discrete schemes are accurate and robust for nonsmooth data, and competitive with existing schemes.

연구 동기 및 목표

• 시간에 대한 캡투로 분수차분 도함수를 포함하는 시간 분수차분 연속형 확산 및 확산-파동 방정식에 대해 완전 이산화된 수치 스킴을 개발한다.
• 초기 자료의 정규성에 따라 최적의 오차 추정치를 확보하며, 부드럽지 않은 경우도 포함한다.
• 초기 자료의 부드러움 여부에 관계없이 시간에 대해 1차 및 2차 정밀도를 유지하는 스킴을 설계한다.
• 기존의 시간 스텝 방법들과의 상세한 비교를 통해 경쟁력과 강건성을 평가한다.
• 1차원 및 2차원에서의 광범위한 수치 실험을 통해 이론적 수렴 속도를 검증한다.

제안 방법

• 변수 계수와 복잡한 기하구조를 다룰 수 있도록 공간 이산화에 갈레르킨 유한요소법을 적용한다.
• 시간에 대해 1차 정밀도를 확보하기 위해 암시적 후진 오일러 방법을 적용한다.
• 향상된 시간 정밀도와 안정성을 확보하기 위해 2차 후진 차분법을 사용한다.
• 공간 반이산화와 시간 스텝 방법을 조합하여 완전 이산화 스킴을 구성한다.
• 초기 자료의 정규성에 따라 최적의 오차 추정치를 유도한다.
• 이론적 수렴 속도를 검증하기 위해 1차원 및 2차원 문제에 대해 스킴을 구현하고 시험한다.

실험 결과

연구 질문

• RQ1 시간 분수차분 연속형 확산 방정식에서 부드럽고 부드럽지 않은 초기 자료 모두에 대해 최적의 수렴 속도를 유지하는 완전 이산화 스킴을 개발할 수 있는가?
• RQ2 기존의 시간 스텝 방법들과 비교해 볼 때, 제안된 스킴은 분수차분 확산 및 확산-파동 문제에서 정확성과 강건성 측면에서 어떻게 비교되는가?
• RQ3 초기 자료가 정규성이 없을 경우, 특히 캡투로 분수차분 도함수의 맥락에서 스킴의 수렴 거동은 어떠한가?
• RQ4 갈레르킨 유한요소법과 후진 시간 이산화의 조합이 다양한 공간 차원에서 안정적이고 정확한 해를 도출할 수 있는가?
• RQ5 수치 결과가 스킴에 대해 유도된 이론적 오차 추정치를 어느 정도 확인하는가?

주요 결과

• 제안된 스킴은 부드럽고 부드럽지 않은 초기 자료 모두에 대해 시간에 대해 1차 및 2차 정밀도를 달성하며, 이론적 기대와 일치한다.
• 최적의 오차 추정치가 확립되어 있으며, 초기 자료의 정규성에 따라 최적의 수렴 속도를 보여준다.
• 1차원 및 2차원에서의 수치 실험은 부드럽고 부드럽지 않은 해 모두에 대해 예측된 수렴 속도를 확인한다.
• 스킴은 부드럽지 않은 초기 자료에 대해 강건하여, 다른 스킴이 성능이 저하될 수 있는 상황에서도 정확성을 유지한다.
• 세부적인 비교를 통해 제안된 스킴이 정확성과 안정성 측면에서 기존의 시간 스텝 방법들과 경쟁 가능하다.
• 구현은 안정적이고 효율적이며, 다차원 문제에 대한 실용적 응용을 지원한다.