[논문 리뷰] Positivity-preserving and energy-dissipative finite difference schemes for the Fokker-Planck and Keller-Segel equations

Jingwei Hu, Xiangxiong Zhang|arXiv (Cornell University)|2021. 03. 31.

Mathematical Biology Tumor Growth참고 문헌 22인용 수 4

한 줄 요약

이 논문은 포화-플랑크 및 켈러-세겔 방정식에 대해 공간에서 두 번째 및 네 번째 차수 정확도를 갖는 반-암시적 유한차분 스킴을 제안한다. 이 스킴들은 양수성 보존 및 에너지 소산 성질을 가지며, 네 번째 차수 스킴은 시간 스텝과 메쉬 제약 조건이 약간의 조건을 만족할 경우 이러한 성질을 확보한다. 이는 양수성과 에너지 감소를 동시에 보존하는 최초의 고차수 공간 이산화 방법이다.

ABSTRACT

In this work, we introduce semi-implicit or implicit finite difference schemes for the continuity equation with a gradient flow structure. Examples of such equations include the linear Fokker-Planck equation and the Keller-Segel equations. The two proposed schemes are first order accurate in time, explicitly solvable, and second order and fourth order accurate in space, which are obtained via finite difference implementation of the classical continuous finite element method. The fully discrete schemes are proven positivity-preserving and energy-dissipative: the second order scheme can achieve so unconditionally; the fourth order scheme only requires a mild time step and mesh size constraint. Furthermore, the fourth order scheme is the first high order spatial discretization that can achieve both positivity and energy decay properties, which is suitable for long time simulation and to obtain accurate steady state solutions.

연구 동기 및 목표

• 포화-플랑크 및 켈러-세겔과 같은 경량 흐름 방정식에 대해 장시간 안정성을 확보하는 구조 보존 수치 스킴을 개발한다.
• 이산 해가 비음이 아니고 에너지가 단조 감소하도록 보장하여 연속적인 변분 구조를 모방한다.
• 양수성과 에너지 소산을 유지하는 고차수 공간 이산화 방법(두 번째 및 네 번째 차수)을 구축하여 시간 스텝에 대한 엄격한 제약 없이 적용 가능하게 한다.
• 화학성향 및 확산 모델의 정확한 장시간 시뮬레이션과 신뢰할 수 있는 안정 상태 계산을 가능하게 한다.

제안 방법

• 안정성을 확보하기 위해 반-암시적 1차 시간 이산화를 사용하며, 시스템을 완전히 암시적 또는 반-암시적으로 풀어 해를 도출한다.
• 균일한 메쉬 위에서 선형 및 2차 다항식을 사용한 연속 유한요소 방법의 유한차분 근사로 제2 및 제4차 공간 정확도를 확보한다.
• 스킴은 연속 방정식의 등가 변수 계수 확산 형태에서 유도된다: ∂ₜρ = ∇·(ℳ∇(ρ/ℳ)), 여기서 ℳ = exp(−(𝒱 + 𝒲∗ρ)).
• 양수성은 이산 스킴에 대한 단조성 조건을 통해 강제하며, 이는 이산 최대원리 증명을 통해 입증된다.
• 에너지 소산은 이산 자유 에너지 기능이 시간에 따라 감소함을 보여줌으로써 확립되며, 이는 연속적인 에너지 감소를 반영한다.
• 제4차 스킴은 단조성 확보를 위해 시간 스텝에 하한선 제약만 필요하며, 상한선 제약은 필요로 하지 않는다.

실험 결과

연구 질문

• RQ1 고차수 유한차분 스킴이 포화-플랑크 및 켈러-세겔 방정식에 대해 양수성과 에너지 소산을 동시에 보존할 수 있는가?
• RQ2 단조성과 양수성을 보장하기 위해 필요한 시간 및 공간 이산화 제약 조건은 무엇인가?
• RQ3 제안된 스킴은 장시간 시뮬레이션 및 폭발 근처 영역에서 어떻게 성능을 발휘하는가?
• RQ4 메쉬 제약 조건이 완화되었을 때도 제4차 스킴이 양수성과 에너지 감소를 유지할 수 있는가?

주요 결과

• 제2차 스킴은 시간 스텝 제약 없이 조건부로 양수성 보존 및 에너지 소산을 보장하여 안정성을 확보한다.
• 제4차 스킴은 시간 스텝이 하한선으로 제약되어야만 양수성과 에너지 소산이 유지되며, 상한선 제약은 필요로 하지 않는다.
• 수치 결과는 두 스킴 모두 비임계 초깃질량 조건에서 T ≈ 13.52 시점에 ∥ρⁿ⁺¹ − ρⁿ∥∞ ≤ 10⁻⁸ 수준으로 안정 상태에 수렴함을 보여준다.
• 초임계 초깃질량 조건에서는 제4차 스킴이 제2차 스킴보다 폭발 역학을 더 정확하게 포착하며, 특히 T = 0.8 시점에서 두드러진 성능을 보인다.
• 에너지 변화도는 폭발 이후에도 두 스킴 모두 단조 감소함을 확인하였으며, 에너지 소산 또는 양수성 위반은 발생하지 않았다.
• 제4차 스킴은 양수성과 에너지 감소를 동시에 보존하는 최초의 고차수 공간 이산화 방법으로, 정확한 장시간 시뮬레이션을 가능하게 한다.