[논문 리뷰] Numerical solutions to boundary value problem for anomalous diffusion equation with Riesz-Feller fractional operator

Mariusz Ciesielski, Jacek Leszczyński|arXiv (Cornell University)|2006. 07. 05.

Fractional Differential Equations Solutions참고 문헌 8인용 수 28

한 줄 요약

이 논문은 일차원 공간에서 Riesz-Feller 분수 도함수를 포함하는 경계값 문제를 해결하기 위한 분수 유한차분법(FFDM)을 제안한다. 이는 비정상 확산의 수치 시뮬레이션을 가능하게 한다. 이 방법은 고전적 유한차분 스킴을 일반화한 것으로, 주요 결과로는 α < 2일 때 비선형 온도 프로파일이 나타나며, 실험적 나노튜브 열전도 데이터를 α = 0.35 및 θ = -0.055로 정확하게 피팅함을 보여준다.

ABSTRACT

In this paper, we present a numerical solution to an ordinary differential equation of a fractional order in one-dimensional space. The solution to this equation can describe a steady state of the process of anomalous diffusion. The process arises from interactions within complex and non-homogeneous background. We present a numerical method which is based on the finite differences method. We consider a boundary value problem (Dirichlet conditions) for an equation with the Riesz-Feller fractional derivative. In the final part of this paper, same simulation results are shown. We present an example of non-linear temperature profiles in nanotubes which can be approximated by a solution to the fractional differential equation.

연구 동기 및 목표

• 비정상 확산에서 Riesz-Feller 분수 도함수를 포함하는 경계값 문제를 해결하기 위한 수치적 방법을 개발한다.
• 표준 확산이 실패하는 복잡하고 비균질적인 매질(예: 나노튜브)에서 정적인 온도 프로파일을 모델링한다.
• 장거리 공간적 의존성과 무거운 尾部 입자 점프를 포착하는 분수계수 도함수로 고전적 유한차분 방법을 확장한다.
• 나노튜브에서의 실험 데이터와 수치적 해를 비교하여 방법을 검증한다.
• 기울기 파ameter θ와 순서 α가 해의 대칭성과 형태에 미치는 영향을 조사한다.

제안 방법

• 공간적 Riesz-Feller 분수 도함수 연산자를 이산화하기 위해 유한차분법(FDM)을 사용한다.
• α와 θ에 따라 달라지는 계수를 가진 함수 값의 가중합을 이용해 Riesz-Feller 도함수의 수치적 근사를 유도한다.
• 행렬 A와 벡터 B를 사용한 선형계 A·T = B를 구성하며, 행렬 A는 분수 도함수 스텐실을 표현하고, 벡터 B는 딜리클레 경계 조건을 포함한다.
• α = 2 및 θ = 0일 경우, 스킴은 두 번째 도함수에 대한 고전적 중심차분 방법으로 축소된다.
• 이 방법은 국소적이지 않은 행동을 고려한다: 각 점의 도함수는 국소적 이웃 값뿐만 아니라 전체 도메인의 값에 의존한다.
• 선형계를 풀어 해를 계산하며, 분수 도함수의 비국소성으로 인해 경계 조건이 모든 내부 노드에 영향을 미친다.

실험 결과

연구 질문

• RQ1 유한차분 스킴을 사용하여 경계값 문제에서 Riesz-Feller 분수 도함수를 정확하게 근사할 수 있는가?
• RQ2 분수 순서 α와 기울기 파ameter θ를 변화시켰을 때 정적 온도 프로파일의 형태에 어떤 영향을 미치는가?
• RQ3 제안된 분수 유한차분 방법이 표준 확산이 실패하는 나노튜브에서 실험적 온도 프로파일을 정확하게 재현할 수 있는가?
• RQ4 α가 1에 수렴하고 θ가 ±1에 수렴할 경우 해의 행동은 어떻게 변화하며, 이는 어떤 물리적 해석을 갖는가?
• RQ5 분수 도함수 파ameter와 비정상 확산에서 장꼬리 확률 분포의 발생 간의 관계는 무엇인가?

주요 결과

• α = 2 및 θ = 0일 경우 해는 선형이며, 고전적 열방정식과의 일致성을 확인한다.
• α < 2일 경우 해는 비선형 프로파일을 보이며, 이는 비정상 확산 행동을 나타낸다.
• α → 1+ 및 θ → ±1+일 때 해는 일阶 파동방정식의 정적 상태에 수렴한다.
• 기울기 파ameter θ는 해에 비대칭성을 도입하며, θ ∈ (0,1)일 경우 비대칭 프로파일을 생성한다.
• α = 0.35 및 θ = -0.055인 모델이 장과 리(2005)의 실험적 나노튜브 온도 데이터와 가장 잘 일치한다.
• 분수 도함수는 장꼬리 입자 점프를 포착하여, 가우시안 통계로는 기대할 수 없는 희귀하지만 극단적인 사건의 모델링이 가능하다.