[논문 리뷰] Unconditional Stability Of A Two-Step Fourth-Order Modified Explicit Euler/Crank-Nicolson Approach For Solving Time-Variable Fractional Mobile-Immobile Advection-Dispersion Equation

Eric Ngondiep|arXiv (Cornell University)|2022. 05. 06.

Fractional Differential Equations Solutions인용 수 5

한 줄 요약

이 논문은 시간에 따라 변화하는 분수계수 이동-비이동 확산 방정식에 대해 변수계수 캄프로 분수도수를 갖는 두 단계의 4차 수정형 명시적 오일러/크랭크-니콜슨 스킴을 제안한다. 이 방법은 조건부 안정성 없이 안정성과 수렴 차수 $O(k + h^4)$를 확보하여 기존의 단일 단계 스킴에 비해 뛰어난 효율성과 정확성을 보여준다.

ABSTRACT

This paper considers a two-step fourth-order modified explicit Euler/Crank-Nicolson numerical method for solving the time-variable fractional mobile-immobile advection-dispersion model subjects to suitable initial and boundary conditions. Both stability and error estimates of the new approach are deeply analyzed in the $L^{\infty}(0,T;L^{2})$-norm. The theoretical studies show that the proposed technique is unconditionally stable with convergence of order $O(k+h^{4})$, where $h$ and $k$ are space step and time step, respectively. This result indicate that the two-step fourth-order formulation is more efficient than a broad range of numerical schemes widely studied in the literature for the considered problem. Numerical experiments are performed to verify the unconditional stability and convergence rate of the developed algorithm.

연구 동기 및 목표

• 시간에 따라 변화하는 분수계수 이동-비이동 확산 방정식에 대해 고차수이고 조건부 안정성이 없는 수치적 방법을 개발하기 위해.
• 일반적으로 시간 정확도가 낮거나 조건부 안정성으로 인해 제한되는 기존의 단일 단계 스킴의 한계를 극복하기 위해.
• 계산 효율성을 유지하면서 공간적으로 4차 정확도, 시간적으로 1차 수렴을 달성하기 위해.
• 제안된 스킴에 대해 $L^{ty}(0,T;L^{2})$-노름에서 엄밀한 안정성 및 오차 분석을 제공하기 위해.
• 정확한 해가 있는 수치 실험과 수렴 속도 분석을 통해 이론적 결과를 검증하기 위해.

제안 방법

• 두 단계의 4차 수식은 수정형 명시적 오일러와 크랭크-니콜슨 스킴을 조합하여 정확도와 안정성을 균형 있게 확보한다.
• 시간 도함수는 2차 정확도를 갖는 두 단계 스킴으로 근사하고, 공간 도함수는 4차 고밀도 차분으로 이산화한다.
• 변수계수 캄프로 분수도수는 특이 핵을 포함하는 정의를 기반으로 한 적분 공식을 사용해 근사한다.
• 에너지 방법을 통해 $L^{ty}(0,T;L^{2})$-노름에서 안정성을 분석하여 시간 및 공간 스텝 크기와 무관하게 조건부 안정성이 없음을 입증한다.
• 스킴은 선형 연립방정식으로 표현되며, 효율적인 반복 해법으로 풀 수 있고, 향후 조건부 개선을 위한 예비 조건화 기법 적용 가능.
• 분수도수 항에 대해 반암시적 처리를 통해 안정성을 확보하고 더 큰 시간 스텝 크기를 허용한다.

실험 결과

연구 질문

• RQ1 두 단계의 4차 수치 스킴이 시간에 따라 변화하는 분수계수 이동-비이동 확산 방정식에 대해 조건부 안정성이 없이 안정성을 확보할 수 있는가?
• RQ2 제안된 방법의 수렴 차수는 시간과 공간에서 어떻게 되며, 기존의 단일 단계 스킴보다 뛰어난가?
• RQ3 분수도수가 시간과 공간에 따라 변화할 경우, 이 방법의 안정성과 정확도는 어떻게 되는가?
• RQ4 $L^{ty}(0,T;L^{2})$-노름에서 이론적 오차 한계가 정확한 해가 있는 수치 실험을 통해 검증될 수 있는가?
• RQ5 이 유형의 문제에 대해 널리 사용되는 기존 스킴에 비해 제안된 스킴이 계산적으로 더 효율적인가?

주요 결과

• 제안된 스킴은 $L^{ty}(0,T;L^{2})$-노름에서 조건부 안정성이 없으며, 이는 모든 시간 및 공간 스텝 크기에서 안정성이 유지됨을 의미한다.
• 시간에 대해 1차 수렴($O(k)$)과 공간에 대해 4차 수렴($O(h^4)$)를 달성하여 총 수렴 차수는 $O(k + h^4)$가 된다.
• 수치 실험을 통해 이론적 수렴 속도가 확인되었으며, 여러 테스트 케이스에서 공간 수렴 차수는 약 4.0, 시간 수렴 차수는 3.8~4.1로 관측되었다.
• 예제 1에서 $h = 1/16$일 때 최대 오차는 $1.1488 \times 10^{-6}$이며, 수렴 속도는 4.1029로 4차 공간 정확도가 확인되었다.
• 예제 2에서 $\alpha = 0.25$일 때 $h$가 감소함에 따라 오차는 $1.0918 \times 10^{-2}$에서 $2.5995 \times 10^{-6}$로 감소하며, 수렴 속도는 4.1328로 고차 정확도가 확인되었다.
• 이론적 분석과 수치 결과 모두에서 기존의 광범위하게 사용되는 스킴에 비해 더 효율적이고 빠른 성능을 보였다.