[논문 리뷰] Finite difference method for a general fractional porous medium equation

Félix del Teso, Juan Luís Vázquez|arXiv (Cornell University)|2013. 07. 09.

Differential Equations and Numerical Methods참고 문헌 14인용 수 21

한 줄 요약

이 논문은 일반적인 지수 $ m \geq 1 $ 와 $ \sigma \in (0,2) $ 를 가진 분수형 다공성 매체 방정식에 대한 유한차분 스킴을 제안한다. 비국소적 분수형 라플라스 연산자를 다루기 위해 Caffarelli-Silvestre 확장 기법을 사용한다. 이 방법은 수렴 차수 $ O(\Delta t + \Delta x^{2-\sigma} + \frac{\Delta x^a}{\Delta t}) $ 를 달성하며, 이때 $ a $ 는 $ \sigma $ 에 따라 달라진다. 이와 함께 이론적 해에 대한 수치해의 존재성, 유일성, 수렴성을 엄밀히 증명한다.

ABSTRACT

We formulate a numerical method to solve the porous medium type equation with fractional diffusion \[ \frac{\partial u}{\partial t}+(-Δ)^{σ/2} (u^m)=0 \] posed for $x\in \mathbb{R}^N$, $t>0$, with $m\geq 1$, $σ\in (0,2)$, and nonnegative initial data $u(x,0)$. We prove existence and uniqueness of the solution of the numerical method and also the convergence to the theoretical solution of the equation with an order depending on $σ$. We also propose a two points approximation to a $σ$-derivative with order $O(h^{2-σ})$.

연구 동기 및 목표

• 비국소적 확산을 포함한 일반적인 분수형 다공성 매체 방정식에 대해 안정적이고 강력한 수치적 스킴을 개발한다.
• 확장 방법을 사용하여 $ \sigma \neq 1 $ 인 경우 분수형 라플라스 연산자를 구현하는 데 발생하는 과제를 다룬다.
• 이론적 해에 대한 수치해의 존재성, 유일성, 수렴성을 증명한다.
• 수렴 속도가 분수형 차수 $ \sigma $ 에 따라 어떻게 달라지는지를 규명한다.
• 유한한 영역에서의 수치해가 영역이 확장될 때 $ \mathbb{R}^N $ 에서의 해로 수렴함을 보장한다.

제안 방법

• 비국소적 확산 항 $ (-\Delta)^{\sigma/2}(u^m) $ 은 Caffarelli-Silvestre 확장을 통해 재구성되어, 비국소 문제를 한 차원 더 높인 국소 PDE 문제로 변환한다.
• 확장된 문제는 유한한 계산 영역 $ \Omega = [-X,X] \times [0,X] $ 에서 적절한 경계 조건을 갖는 유한차분 스킴을 통해 해결된다.
• 확장 공식에 나타나는 가중 미분 $ \partial / \partial y^\sigma $ 에 대해 두 점 유한차분 근사법을 제안하며, 이는 $ O(h^{2-\sigma}) $ 의 정밀도를 가진다.
• 수치 스킴은 안정성과 수렴성을 증명하였으며, 국소 절단 오차는 $ O(\Delta t + \Delta x^{2-\sigma}) $ 로 유계이다.
• 영역 $ \Omega \to \mathbb{R}^N $ 이 되는 조건 하에서, 경계 오차는 Barenblatt 해 추정치를 통해 제어됨으로써 유한한 영역에서의 수치해가 전체 공간 해로 수렴함을 확립한다.
• 이론적 수렴성은 $ \Delta t \leq C(m,f) \Delta x^\sigma $ 라는 가정 하에 증명되었으며, 이 조건은 안정성과 일致성을 보장한다.

실험 결과

연구 질문

• RQ1 일반적인 $ \sigma \in (0,2) $ 와 $ m \geq 1 $ 에 대해 정확하고 안정적인 분수형 다공성 매체 방정식의 해를 근사하는 데 적합한 유한차분 방법을 구성할 수 있는가?
• RQ2 수치 스킴에서 비국소적 분수형 라플라스 연산자를 효율적이고 정확하게 근사하는 방법은 무엇인가?
• RQ3 제안된 수치 스킴의 수렴 속도는 무엇이며, 이는 분수형 차수 $ \sigma $ 에 따라 어떻게 달라지는가?
• RQ4 유한한 영역에서의 수치해가 전체 공간 $ \mathbb{R}^N $ 에서의 해로 수렴하는 방식은 어떠한가?
• RQ5 이론적 해로의 수치해 수렴을 보장하기 위해 도메인 크기와 시간 간격에 어떤 조건이 필요하는가?

주요 결과

• 제안된 유한차분 스킴은 수치해의 존재성과 유일성을 증명하였다.
• 이 방법은 이론적 해로의 수렴 차수 $ O(\Delta t + \Delta x^{2-\sigma} + \frac{\Delta x^a}{\Delta t}) $ 를 달성하며, 이때 $ a $ 는 $ \sigma $ 에 따라 달라진다. $ \sigma \in (0,1] $ 에서는 $ a=2 $, $ \sigma \in (1,3/2] $ 에서는 $ a=3 $, $ \sigma \in (3/2,2) $ 에서는 $ a=4 $ 이다.
• 두 점 근사법을 통한 $ \sigma $-미분은 국소 절단 오차가 $ O(h^{2-\sigma}) $ 수준이며, 이는 주어진 이산화 방법에 대해 최적이다.
• 경계 오차가 $ \Delta x^a $ 로 감소하도록 보장하는 조건 하에, 유한한 영역에서의 수치해가 $ \mathbb{R}^N $ 에서의 해로 수렴함을 입증하였다. 이 조건은 $ |X| \geq L \cdot \Delta x^{-a/(N+\sigma)} $ 를 만족해야 한다.
• 안정성 조건 $ \Delta t \leq C(m,f) \Delta x^\sigma $ 하에서 수렴 속도가 최적임을 증명하였으며, 이 조건은 스킴의 일致성과 안정성을 보장한다.
• 이론적 수렴성은 수치적으로도 검증되었으며, 오차는 시간, 공간, 경계 절단 오차의 합으로 유계이며, 이 모든 오차는 Barenblatt 해 추정치를 통해 제어된다.