[논문 리뷰] Energy and implicit discretization of the Fokker-Planck and Keller-Segel type equations

Luís Almeida, Federica Bubba외 2명|arXiv (Cornell University)|2018. 03. 28.

Mathematical Biology Tumor Growth참고 문헌 17인용 수 5

한 줄 요약

이 논문은 포화-플랑크 및 켈러-세겔 방정식에 대해 에너지 소산, 질량 보존, 양성, 그리고 정 steady 상태 구조를 이산 수준에서 유지하는 두 가지 새로운 암묵적 유한체적 스킴을 제안한다. 기울기 유량 및 색프터-기멜 이산화 전략을 활용하고 상류화를 적용함으로써, 시간 단계 제약 없이 잘 정의되고 안정된 해법을 확보하여 패턴 형성의 정확한 시뮬레이션과 물리적 불안정성과 수치적 오류 사이의 구별을 가능하게 한다.

ABSTRACT

The parabolic-elliptic Keller-Segel equation with sensitivity saturation, because of its pattern formation ability, is a challenge for numerical simulations. We provide two finite-volume schemes whose goals are to preserve, at the discrete level, the fundamental properties of the solutions, namely energy dissipation, steady states, positivity and conservation of total mass. These requirements happen to be critical when it comes to distinguishing between discrete steady states, Turing unstable transient states, numerical artifacts or approximate steady states as obtained by a simple upwind approach. These schemes are obtained either by following closely the gradient flow structure or by a proper exponential rewriting inspired by the Scharfetter-Gummel discretization. An interesting feature is that upwind is also necessary for all the expected properties to be preserved at the semi-discrete level. These schemes are extended to the fully discrete level and this leads us to tune precisely the terms according to explicit or implicit discretizations. Using some appropriate monotony properties (reminiscent of the maximum principle), we prove well-posedness for the scheme as well as all the other requirements. Numerical implementations and simulations illustrate the respective advantages of the three methods we compare.

연구 동기 및 목표

• 포화-플랑크 및 켈러-세겔 방정식에 대해 에너지 소산, 질량 보존, 양성, 그리고 정 steady 상태를 유지하는 수치 스킴을 개발한다.
• 강한 비선형성과 대류 항으로 인한 수치적 과제, 특히 물리적 터닝 불안정성과 수치적 오류를 구별하는 데 도전한다.
• CFL 제약 없이 모든 주요 이산 성질을 유지하는 암묵적 시간 이산화를 구성한다.
• 강한 비선형성을 갖는 해의 정성적 행동을 정확히 포착하는 데 있어 기울기 유량 기반 스킴과 색프터-기멜 기반 스킴의 성능을 비교한다.
• 하향-상향 해석 이론에 기반한 존재성 및 유일성, 그리고 수렴성에 대한 잘 정의된 성질을 증명한다.

제안 방법

• 스킴은 두 가지 다른 대칭화 전략에서 유도된다: 기울기 유량 구조와 이동-확산 방정식에 대한 색프터-기멜 방법.
• 양성, 질량 보존, 에너지 소산을 보존하기 위해 두 스킴에 모두 상류화를 통합한다.
• 완전히 이산화된 스킴은 암묵적 시간 스텝을 사용하여 시간 단계 제약 없이 무조건 안정성을 확보한다.
• 최대 원리와 유사한 단조성 성질을 활용하여 하향-상향 해석 이론을 통해 해의 존재성과 유일성을 증명한다.
• 이산 에너지 소산은 연속적인 에너지 부등식을 반영하는 리아프노프 함수를 구성함으로써 엄밀하게 증명된다.
• 수치적 검증을 통해 표준 상류화 방법 대비 정 steady 상태 및 일시적 패턴을 더 정확하게 해결하는 것으로 나타났다.

실험 결과

연구 질문

• RQ1 암묵적 유한체적 스킴은 포화-플랑크 및 켈러-세겔 방정식에서 에너지 소산, 질량 보존, 양성, 그리고 정 steady 상태 구조를 유지할 수 있는가?
• RQ2 기울기 유량 기반과 색프터-기멜 기반 이산화 방법은 강한 비선형성을 갖는 해의 정성적 행동을 어떻게 비교하여 유지하는가?
• RQ3 상류화만으로도 반디셀 수준에서 핵심 물리적 성질을 유지할 수 있는가?
• RQ4 암묵적 시간 이산화로 인해 CFL 조건을 제거할 수 있으며, 모든 이산 보존 및 안정성 성질을 유지할 수 있는가?
• RQ5 이러한 시스템에서 잡음 진동이나 잘못된 정 steady 상태와 같은 수치적 오류는 어떻게 물리적 패턴 형성과 구별할 수 있는가?

주요 결과

• 제안된 스킴은 반디셀 및 완전히 이산화된 수준에서 에너지 소산을 유지하며, 에너지 기능이 시간에 따라 단조 감소한다.
• 두 스킴 모두 해의 엄격한 양성과 총 질량 보존을 유지하며, 매끄럽지 않거나 피크가 강한 초기 자료에 대해서도 성립한다.
• 스킴은 무조건 안정적이고 잘 정의되어 있으며, 하향-상향 해석 이론을 통해 해의 존재성과 유일성이 증명되었다.
• 색프터-기멜 기반 스킴은 일반적인 민감도 함수 φ(u)에 대해 기울기 유량 접근법의 강력한 대안을 제공한다.
• 수치적 시뮬레이션은 스킴이 터닝 유형의 불안정성과 피크형 정 steady 상태를 정확히 포착하며, 표준 상류화 스킴에서 흔히 발생하는 잡음 진동을 피하는 것으로 나타났다.
• 암묵적 시간 이산화로 인해 안정성이나 정확도를 잃지 않고 큰 시간 단계를 사용할 수 있어, 패턴 형성의 장기 시뮬레이션에 적합하다.