[논문 리뷰] High-order numerical method for two-dimensional Riesz space fractional advection-dispersion equation

A. Borhanifar, Maria Alessandra Ragusa외 1명|arXiv (Cornell University)|2020. 06. 07.

Fractional Differential Equations Solutions참고 문헌 27인용 수 4

한 줄 요약

이 논문은 두 차원 Riesz 공간 분수 적분 운반-분산 방정식을 해결하기 위해 분수 중심 차분과 교대 방향 음성법(ADI)을 기반으로 한 네째 계수 정확도의 유한 차분 방법을 제안한다. 이 방법은 네째 계수 공간 이산화와 둘레 차수 크랭크-니콜슨 시간 스킴을 결합하여 공간에서 네째 계수 수렴성과 시간에서 둘레 계수 수렴성을 달성하며, 높은 정확도와 안정성을 갖춘 수치 예제를 통해 검증되었다.

ABSTRACT

In this paper, by combining of fractional centered difference approach with alternating direction implicit method, we introduce a mixed difference method for solving two-dimensional Riesz space fractional advection-dispersion equation. The proposed method is a fourth order centered difference operator in spatial directions and second order Crank-Nicolson method in temporal direction. By reviewing the consistency and stability of the method, the convergence of the proposed method is achieved. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed technique.

연구 동기 및 목표

• 두 차원 Riesz 공간 분수 운반-분산 방정식(RSFADE)을 해결하기 위한 고차수 수치 방법을 개발하기 위해, 비국소성과 분수 도함수의 낮은 정규성으로 인해 도전적인 문제를 다루기 위함.
• 분수 중심 차분과 교대 방향 음성법(ADI) 접근법을 조합하여 2D RSFADE에 대한 효율적이고 정확한 수치 스킴의 부족 문제를 해결하기 위함.
• 이론적 분석을 통해 안정성과 수렴성을 보장하면서 공간에서 네째 계수 정확도와 시간에서 둘레 계수 정확도를 달성하기 위함.
• 기존 정확한 해가 알려진 수치 실험을 통해 방법의 성능과 수렴 차수를 종합적으로 검증하기 위함.

제안 방법

• 방정식의 x 및 y 방향에서 Riesz 공간 분수 도함수에 대해 네째 계수 분수 중심 차분 근사법을 적용한다.
• 시간 이산화에 대해 크랭크-니콜슨 방법을 적용하여 둘레 계수 정확도와 무조건적 안정성을 확보한다.
• 이중 차원 문제를 한 차원 부분 문제로 나누어 해결하기 위해 교대 방향 음성(ADI) 기법을 사용하여 계산 효율성을 향상시킨다.
• 고차수 정확도를 유지하면서도 결과 선형 시스템을 효율적으로 해결할 수 있도록 수정된 크랭크-니콜슨 ADI 체계를 구성한다.
• von Neumann 안정성 분석과 절단 오차 추정을 통해 안정성과 일致성 분석을 수행하고, Lax 등가 정리에 의해 수렴성 증명을 이끌어낸다.
• 균일한 공간 및 시간 격자로 구현되며, 문제 크기에 따라 반복적 또는 직접적 방법으로 선형 시스템을 해석한다.

실험 결과

연구 질문

• RQ1 두 차원 Riesz 공간 분수 운반-분산 방정식에 대해 고차수 수치 스킴을 개발할 수 있는가? 이는 공간 정확도와 시간 정확도를 모두 유지할 수 있는가?
• RQ2 분수 중심 차분과 ADI-Crank-Nicolson 스킴의 조합이 2D RSFADE에 대해 안정적이고 수렴 가능한 방법을 제공하는가?
• RQ3 제안된 방법의 실제 수렴 속도는 공간 및 시간 방향에서 얼마이며, 이는 이론적 순서와 일치하는가?
• RQ4 기존 정확한 해가 알려진 2D RSFADE를 해결할 때, 이 방법은 얼마나 효율적이고 정확한가? 특히 격자 세분화가 다양할 경우 어떻게 되는가?

주요 결과

• 공간에서 네째 계수 수렴성이 수치 실험을 통해 확인되었으며, 공간 스텝 크기를 감소시킬수록 수렴 속도가 약 3.9~3.96 수준을 보였다.
• 시간에서 둘레 계수 수렴성이 관측되었으며, 시간 스텝 크기를 세분화할 경우 추정된 수렴 속도는 1.92~1.96 수준이었고, 이는 이론적 순서와 일치하였다.
• 최대 절대 오차는 격자 세분화에 따라 크게 감소하였다. 예를 들어, $ h_x = h_y = 0.00625\theta $ 일 때, 오차는 6.2번 예제에서 $ 8.13 \times 10^{-9} $ 수준으로 감소하였다.
• 변동하는 분수 차수와 확산 계수를 포함한 다양한 파rameter 세트에서도 안정성과 정확도를 유지하여 강건성을 입증하였다.
• 복잡한 소스 항과 비연속적인 초기 조건이 존재하는 경우에도 낮은 오차 수준과 일관된 수렴 행동을 보여, 계산 효율성이 확인되었다.
• 이론적 분석을 통해 일치성과 무조건적 안정성이 확인되었으며, Lax 등가 정리에 의해 수치 스킴의 수렴성이 보장되었다.