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

Félix del Teso|arXiv (Cornell University)|2013. 01. 18.

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

한 줄 요약

이 논문은 분수 마찰 중간 매질 방정식에 대해 (-Δ)^1/2 확산을 갖는 유한 차분 방법을 제안한다. Caffarelli-Silvestre 확장을 통해 비국소 문제를 한 차원 높은 차원에서 국소 문제로 변환한다. 이 방법은 존재성, 유일성, 진짜 해로의 수렴성을 입증하며, 수치 실험을 통해 2차 수렴 속도와 두꺼운 尾행동을 정확히 포착하는 것을 확인한다.

ABSTRACT

We formulate a numerical method to solve the porous medium type equation with fractional diffusion \[\frac{\partial u}{\partial t}+(-Δ)^{1/2} (u^m)=0.\] The problem is posed in $x\in \mathbb{R}^N$, $m\geq 1$ and with nonnegative initial data. The fractional Laplacian is implemented via the so-called Caffarelli-Silvestre extension. We prove existence and uniqueness of the solution of this method and also the convergence to the theoretical solution of the equation. We run numerical experiments on typical initial data.

연구 동기 및 목표

• 분수 마찰 중간 매질 방정식에 대해 (-Δ)^1/2 확산을 갖는 안정적인 수치적 스킴을 개발한다.
• 분수 확산의 비국소성 문제를 Caffarelli-Silvestre 확장을 통해 고차원에서 국소 문제로 변환함으로써 해결한다.
• 최대 원리(최대 원리)를 사용하여 유한 차분 해의 존재성, 유일성, 수렴성을 입증한다.
• 다양한 초기 자료(디랙 델타 함수 및 컴팩트 지지 함수 포함)에 대한 수치 실험을 통해 방법을 검증한다.
• u^m를 초월하여 일반적인 단조 증가 비선형성 φ(u)에 대해, 로그 확산을 포함한 일반적인 비선형성으로 방법을 확장한다.

제안 방법

• 분수 라플라시안은 Caffarelli-Silvestre 확장을 통해 구현되며, 이는 비국소 문제를 R^N × R^+에서의 국소 경계값 문제로 변환한다.
• 확장된 문제를 유한 도메인 Ω = (-X,X)^N × (0,Y)에서 수치적으로 해결하며, 측면 경계 Γ_h에서 동차 딜리클레 조건을 적용한다.
• 확장된 방정식에 대해 유한 차분 스킴을 적용하며, 시간 이산화에는 반-암시적 방법을, 공간 이산화에는 중심 차분을 사용한다.
• 시간 진동은 양수성을 유지하고 이산 최대 원리를 만족하는 단조 스킴을 통해 처리된다.
• 시간 도함수 조건을 φ와 그 역함수를 사용하여 재구성함으로써 일반 φ(u)로의 확장을 달성한다.
• 이산 최대 원리와 안정성 추정을 통해 수렴성을 입증하며, m=1인 경우 정확한 해와 정교한 격자로 오차 분석을 수행한다.

실험 결과

연구 질문

• RQ1 (-Δ)^1/2 확산을 갖는 분수 마찰 중간 매질 방정식에 대해 수렴성과 안정성을 보장하는 유한 차분 방법을 구성할 수 있는가?
• RQ2 Caffarelli-Silvestre 확장은 비국소 문제의 실용적 수치적 해법을 어떻게 가능하게 하는가?
• RQ3 m의 다양한 값과 초기 자료에 대해 제안된 유한 차분 스킴의 수렴 속도는 어떠한가?
• RQ4 이 방법은 로그 확산을 포함한 일반적인 비선형성 φ(u)로 확장될 수 있는가?
• RQ5 이 방법은 분수 확산을 갖는 마찰 중간 매질 방정식의 특징인 두꺼운 尾행동을 얼마나 정확히 포착하는가?

주요 결과

• 이 유한 차분 방법은 이산 최대 원리를 통해 이산 해의 존재성과 유일성을 입증한다.
• 메쉬 크기 Δx와 시간 간격 Δt가 0으로 갈수록 방법은 진짜 해로 수렴한다.
• 수치 실험 결과에서 2차 수렴 속도를 보였다: m=1일 때, 200노드 격자에서 Δx가 1에서 0.25로 감소할 때 오차는 0.011에서 0.00043으로 감소한다.
• m=2일 때도 유사한 정밀도 개선으로 오차가 0.0277에서 0.00085로 감소하여 수렴을 확인하였으며, 감소 속도는 느리지만 여전히 수렴한다.
• 이 방법은 특히 작은 m에 대해 해의 두꺼운 尾행동을 정확히 포착하며, m=1과 m=10일 때 Barenblatt 프로파일을 재현한다.
• φ(u) = log(u+1)로의 확장도 성공적으로 수행되어 일반적인 단조 증가 비선형성에 대한 강건성을 입증한다.