[논문 리뷰] Structure preserving schemes for nonlinear Fokker-Planck equations with anisotropic diffusion

Nadia Loy, Mattia Zanella|arXiv (Cornell University)|2019. 05. 08.

Advanced Thermodynamics and Statistical Mechanics인용 수 1

한 줄 요약

이 논문은 비선형 Fokker-Planck 방정식에 대한 고차수, 구조를 유지하는 수치적 스킴을 제안하며, 이는 해의 비음성, 엔트로피 소산, 정상 상태 근사의 정확성을 유지한다. 스킴은 일시적 영역에서 최소 제2차 정밀도를 확보하고 장시간에 걸쳐 고차수 정밀도를 확보하여 집단 행동 및 생명 과학 모델에서의 비정상적 확산의 견고한 시뮬레이션을 가능하게 한다.

ABSTRACT

In this work we propose novel numerical schemes for nonlinear Fokker-Planck-type equations with anisotropic diffusion matrix that preserve fundamental structural properties like non negativity of the solution, entropy dissipation and which guarantees an arbitrarily accurate approximation of the steady state of the problem. All the methods presented are at least second order accurate in the transient regimes and high order for large times. Applications of the schemes to models for collective phenomena and life sciences are considered, in these examples anomalous diffusion is often observed and must be taken into account in realistic models.

연구 동기 및 목표

• 비선형 Fokker-Planck 방정식에 대한 수치적 스킴을 개발하여, 비음성, 엔트로피 소산, 정확한 정상 상태 근사와 같은 핵심 구조적 성질을 유지한다. 이는 비대칭 확산을 수반한다.
• 표준 스킴이 물리적 일관성을 유지하지 못할 수 있는 집단 현상 및 생명 과학 모델에서의 비정상적 확산을 모델링하는 데 도전한다.
• 일시적 영역에서 최소 제2차 정밀도를 확보하고 장시간에 걸쳐 고차수 정밀도를 확보하여 장기 시뮬레이션의 정밀도를 향상시킨다.

제안 방법

• 비음성과 엔트로피 구조를 유지하기 위해, 반음성 시간 이산화와 유한체적 공간 이산화를 조합한 방법을 사용한다.
• 해의 변수에 비선형 변환을 적용하여 비음성을 보장하고, 비대칭 확산 존재 시에 스킴의 안정성을 확보한다.
• 유동에 대해 고차수 재구성 기법을 사용하여 시간에 대해 제2차 정밀도와 장시간에 대해 고차수 정밀도를 확보한다.
• 연속 방정식의 구조에서 유도된 이산 엔트로피 부등식을 통해 엔트로피 안정성을 확보한다.
• 복잡한 비대칭 확산 텐서가 존재하더라도, 스킴은 구조적 설계에 의해 정확한 정상 상태로 수렴할 수 있으며, 임의의 정밀도로 수렴한다.
• 집단 운동 및 생물학적 시스템에서 비정상적 확산을 보이는 모델에 대한 적용을 통해 방법의 타당성을 검증한다.

실험 결과

연구 질문

• RQ1 비선형 Fokker-Planck 방정식에 대한 수치적 스킴은 비대칭 확산이 존재하는 상황에서 해의 비음성을 유지하면서도 고차수 정밀도를 확보할 수 있는가?
• RQ2 이산 스킴에서 어떤 구조 보존 기법이 연속 문제의 특성과 유사하게 엔트로피 소산을 보장하는가?
• RQ3 복잡한 비대칭 확산이 존재하더라도, 스킴이 임의의 정밀도로 정상 상태를 정확히 근사할 수 있는가?
• RQ4 스킴은 일시적 영역과 장시간 영역에서 어느 정도의 고차수 정밀도를 유지하는가?
• RQ5 제안된 방법은 집단 행동 및 생명 과학 모델에서 비정상적 확산을 효과적으로 시뮬레이션할 수 있는가?

주요 결과

• 제안된 스킴은 임의의 비대칭 확산 행렬 조건에서도 해의 비음성을 유지하여 물리적 일관성을 보장한다.
• 스킴은 이산 엔트로피 소산 법칙을 만족하여 연속 시스템의 열역학적 거동를 모방한다.
• 일시적 영역에서는 제2차 정밀도를 확보하고 장시간에 걸쳐 고차수 정밀도를 관찰하여 장기 시뮬레이션의 신뢰성을 향상시킨다.
• 스킴은 시간 간격의 크기에 관계없이 구조적 설계 덕분에 정확한 정상 상태로 임의의 정밀도로 수렴한다.
• 집단 운동 및 생명 과학 모델에 대한 적용 결과, 비정상적 확산 현상을 견고하고 정확하게 캐릭터라이즈하는 데 성공하였다.