[논문 리뷰] A variational finite volume scheme for Wasserstein gradient flows

Clément Cancès, Thomas Gallouët외 1명|arXiv (Cornell University)|2019. 07. 18.

Geometric Analysis and Curvature Flows참고 문헌 68인용 수 27

한 줄 요약

이 논문은 JKO 스킴을 통한 시간 이산화와 공간 이산화된 상류 이동도를 갖는 이중점 유량 근사법을 조합하여 워셔슈타인 기울기 흐름을 해결하기 위한 새로운 변분 유한체적 스킴을 제안한다. 이 방법은 연속 문제의 변분 구조를 유지하며, 음이 아닌 값을 보장하고 에너지 감소를 보장하며, 선형 포커-플랑크 방정식에 대해 수렴성이 입증된다. 시간 및 공간에서 일阶 정확도를 달성하고 다양한 에너지와 초기 프로파일에 대해 강건하다.

ABSTRACT

We propose a variational finite volume scheme to approximate the solutions to Wasserstein gradient flows. The time discretization is based on an implicit linearization of the Wasserstein distance expressed thanks to Benamou-Brenier formula, whereas space discretization relies on upstream mobility two-point flux approximation finite volumes. Our scheme is based on a first discretize then optimize approach in order to preserve the variational structure of the continuous model at the discrete level. Our scheme can be applied to a wide range of energies, guarantees non-negativity of the discrete solutions as well as decay of the energy. We show that our scheme admits a unique solution whatever the convex energy involved in the continuous problem, and we prove its convergence in the case of the linear Fokker-Planck equation with positive initial density. Numerical illustrations show that it is first order accurate in both time and space, and robust with respect to both the energy and the initial profile.

연구 동기 및 목표

• PDE(포커-플랑크 방정식 및 다공성 매체 방정식 등)에서 나타나는 워셔슈타인 기울기 흐름을 위한 강건하고 구조를 유지하는 수치적 스킴을 개발하기 위해.
• 먼저 이산화한 후 최적화하는 접근 방식을 통해 연속 모델의 변분 구조를 이산 수준에서 유지하기 위해.
• 일반적인 볼록 에너지에 대해 이산 음이 아닌 성질과 에너지 감소를 보장하기 위해.
• 양의 초기 밀도를 가진 선형 포커-플랑크 방정식에 대해 스킴의 수렴성을 입증하기 위해.
• 다양한 물리적 상황에서의 수치 실험을 통해 강건성과 일阶 정확도를 입증하기 위해.

제안 방법

• 시간 이산화는 JKO 반이산화 기반으로 하며, 벤아모-브레니에 역동적 공식을 통해 워셔슈타인 거리가 포함된 최소화 문제로 문제를 재구성한다.
• 공간 이산화에는 상류 이동도를 갖는 이중점 유량 근사법을 사용한 유한체적 방법을 적용하여 대류 지배 흐름을 처리한다.
• 스킴은 전진 연속 방정식과 후진 해밀턴-자코비 방정식을 포함하는 원시-이중 최적화 문제에서 유도되며, 이는 평균장 게임 시스템으로 이어진다.
• 질량 보존 조건 하에 이산 에너지를 최소화하고, 이중 문제는 이중 변수(해밀턴-자코비 포텐셜)에 대한 고정점 반복을 통해 해결한다.
• 모든 볼록 에너지에 대해 유일한 해가 존재하도록 보장하며, 이산 벤아모-브레니에 공식의 버전을 통해 변분 구조를 유지한다.
• 비구조 메쉬 위에서 조각별 상수 유한체적 수식을 사용하여 스킴을 구현하고, 최적성 조건을 만족시키기 위해 이중 변수를 반복적으로 갱신한다.

실험 결과

연구 질문

• RQ1 이산 수준에서 워셔슈타인 기울기 흐름의 변분 구조를 유지하는 유한체적 스킴을 구성할 수 있는가?
• RQ2 제안된 스킴은 일반적인 볼록 에너지에 대해 이산 음이 아닌 성질과 에너지 감소를 보장하는가?
• RQ3 양의 초기 밀도를 가진 선형 포커-플랑크 방정식에 대해 스킴이 수렴하는가?
• RQ4 스킴의 시간 및 공간 수렴 속도는 얼마이며, 다양한 초기 프로파일과 에너지 함수에 대해 얼마나 강건한가?
• RQ5 이 스킴은 다공성 매체 내 염수 침투를 모델링하는 방정식계로 확장될 수 있는가?

주요 결과

• 모든 볼록 에너지 함수에 대해 스킴이 유일한 해를 가지며, 이는 에너지 유형에 관계없이 수치적 안정성을 보장한다.
• 선형 포커-플랑크 방정식의 경우, 시간 간격과 메쉬 크기가 0으로 갈수록 연속 문제의 약한 해로 수렴한다.
• 수치 결과는 시간 및 공간에서 일阶 수렴을 보이며, 메쉬가 더 미세해지고 시간 간격이 작아질수록 수렴 프로파일이 향상된다.
• 표준 유한체적 스킴보다 에너지 감소 속도가 더 빠르며, 이는 각 단계에서 에너지 감소를 최대화하도록 설계된 것과 일치한다.
• 메쉬가 미세해지고 시간 간격이 줄어들수록 스킴의 에너지 감소가 연속 시스템의 진짜 감소에 수렴한다.
• 이론적 수렴성이 입증되지 않은 시스템에 대해서도 스킴은 다공성 매체 방정식 및 염수 침투와 같은 복잡한 물리 현상을 성공적으로 시뮬레이션한다.