[논문 리뷰] Convergence of a Finite Volume Scheme for a System of Interacting Species with Cross-Diffusion

José A. Carrillo, Francis Filbet외 1명|arXiv (Cornell University)|2018. 04. 12.

Mathematical Biology Tumor Growth참고 문헌 46인용 수 25

한 줄 요약

이 논문은 교차 확산과 비국소 상호작용을 갖는 두 종의 상호작용 시스템에 대해 양수성 보존 유한체적 스킴을 제안한다. 양수성에 기반한 이산 에너지 추정과 컴팩트니스를 통해 저자들은 수치 스킴이 약한 해로 수렴함을 증명하며, 수치 결과는 일阶 공간 정확도와 특이하거나 강한 확산 영역에서도 강인한 행동을 확인한다.

ABSTRACT

In this work we present the convergence of a positivity preserving semi-discrete finite volume scheme for a coupled system of two non-local partial differential equations with cross-diffusion. The key to proving the convergence result is to establish positivity in order to obtain a discrete energy estimate to obtain compactness. We numerically observe the convergence to reference solutions with a first order accuracy in space. Moreover we recover segregated stationary states in spite of the regularising effect of the self-diffusion. However, if the self-diffusion or the cross-diffusion is strong enough, mixing occurs while both densities remain continuous.

연구 동기 및 목표

• 교차 확산을 갖는 두 종의 상호작용 시스템에 대해 양수성과 질량 보존을 보존하는 수치 스킴을 개발한다.
• 비국소 상호작용과 비선형 확산이 존재하는 상황에서도 수치 스킴이 약한 해로 수렴함을 확립한다.
• 강한 자기확산 또는 교차확산 영역, 특히 특이 잠재력이 존재하는 경우의 스킴 행동을 분석한다.
• 수렴 속도를 수치적으로 검증하고, 분리된 상태와 혼합된 상태를 포함한 정적 상태를 재현한다.
• 메타안정 클러스터링과 에너지 감쇠 역학과 같은 복잡한 패턴 형성 시뮬레이션에서 스킴의 강인성을 보여준다.

제안 방법

• 무료경계 조건을 갖는 일차원 영역에서 국소 보존을 보장하기 위해, 반연속적인 유한체적 스킴을 구성한다.
• 밀도와 상호작용 잠재력의 컨볼루션을 포함하는 유량 표현을 통해 비선형 자기확산 및 교차확산 항을 포함한다.
• 이산 해의 양수성을 강제함으로써 이산 에너지 추정을 도출하며, 이는 컴팩트니스 추론에 필수적이다.
• 이산 에너지 함수는 상호작용, 내부 및 확산 항을 포함하며, 연속 에너지의 구조를 모방한다.
• 이산 에너지 추정과 약한 수렴 추론을 통해 컴팩트니스를 확립하고, 극한이 약한 해로 식별됨을 보인다.
• 장거리 인력과 단거리 반발력을 갖는 잠재력(가우시안 및 이차형)을 사용하여 수치 시뮬레이션을 수행한다.

실험 결과

연구 질문

• RQ1 비국소 상호작용이 존재하는 두 종의 교차확산 시스템에 대해, 양수성과 질량 보존을 보존하는 유한체적 스킴을 구성할 수 있는가?
• RQ2 스킴 수준에서 도출된 이산 에너지 추정이 컴팩트니스를 보장하고, 결과적으로 약한 해로의 수렴을 보장하는가?
• RQ3 스킴의 공간 수렴 속도는 얼마이며, 강한 확산 또는 특이 잠재력 조건에서의 성능은 어떠한가?
• RQ4 서로 다른 인력과 반발력이 존재하는 상황에서 스킴이 상호분리, 혼합, 메타안정 클러스터링을 정확히 포착할 수 있는가?
• RQ5 에너지 감쇠는 시간이 지남에 따라 어떻게 변화하며, 이는 지역적 및 비국소 효과 간의 상호작용을 어떻게 드러내는가?

주요 결과

• 수치 실험 결과, 유한체적 스킴은 오차 노름을 통해 일계 공간 수렴을 달성한다.
• 스킴은 해의 양수성을 보존하며, 이는 이산 에너지 추정을 도출하고 컴팩트니스를 증명하는 데 필수적이다.
• 이산 에너지 추정은 스킴 수준에서 확립되어 약한 극한의 식별과 약한 해로의 수렴을 가능하게 한다.
• 수치 시뮬레이션 결과, 스킴은 자기확산이 정규화 작용을 하는 경우에도 분리된 정적 상태의 형성을 성공적으로 포착한다.
• 강한 자기확산 또는 교차확산 조건에서는 혼합이 발생하지만, 두 밀도는 연속성을 유지하며, 분리에서 균질화로의 전이가 나타난다.
• 에너지 감쇠는 인력-인력 상호작용의 장기적 근처에서 지수적으로 감소하는 것으로 관측되나, 시스템에 대한 분석적 증명은 제공되지 않는다.