[논문 리뷰] A Lagrangian scheme for the incompressible Euler equation using optimal transport

Thomas Gallouët, Quentin Mérigot|HAL (Le Centre pour la Communication Scientifique Directe)|2016. 05. 02.

Navier-Stokes equation solutions참고 문헌 15인용 수 6

한 줄 요약

이 논문은 측度를 보존하는 미분형상의 공간에서 아르노르의 기하학적 해석과 브라이엔의 외재적 접근(최적 운반 이론를 활용)을 조합하여 압축성 없는 유동 방정식에 대한 라그랑주 수치 해법을 제안한다. 이 방법은 반산형 최적 운반 해법을 사용하여 비압축성 조건을 페널티로 구현하며, 정규 해에 수렴함을 보이고 2차원 테스트 케이스(벨트라미 흐름, 켈빈-헬름홀츠 불안정성, 레일리-태일러 불안정성 포함)에서 뛰어난 성능을 보인다.

ABSTRACT

We approximate the regular solutions of the incompressible Euler equation by the solution of ODEs on finite-dimensional spaces. Our approach combines Arnold's interpretation of the solution of Euler's equation for incompressible and inviscid fluids as geodesics in the space of measure-preserving diffeomorphisms, and an extrinsic approximation of the equations of geodesics due to Brenier. Using recently developed semi-discrete optimal transport solvers, this approach yields numerical scheme able to handle problems of realistic size in 2D. Our purpose in this article is to establish the convergence of these scheme towards regular solutions of the incompressible Euler equation, and to provide numerical experiments on a few simple testcases in 2D.

연구 동기 및 목표

• 최적 운반을 통한 비압축성 조건 유지 기반의 라그랑주 수치 해법을 개발한다.
• 해법이 압축성 없는 유동 방정식의 정규 해로 수렴하는 이론적 수렴성을 확립한다.
• 실제 입자 수를 고려한 표준 2차원 문제에 대해 해법을 수치적으로 검증한다.
• 다상 유동에 있어 중력과 같은 외부 힘을 포함시키기 위해 방법을 확장한다.
• 켈빈-헬름홀츠 및 레일리-태일러 불안정성과 같은 복잡한 유체 불안정성을 해법이 어떻게 잘 포착하는지 보여준다.

제안 방법

• 압축성 없는 유동 방정식을 측도를 보존하는 미분형상의 공간에서 기하학적 길로 간주한다.
• 비압축성 조건을 측도를 보존하는 매핑 집합으로부터의 제곱 거리에 의해 페널티로 강제하는 형식을 사용한다.
• 브라이엔의 극분해 정리에 기반해 최적 운반을 통해 측도를 보존하는 매핑 집합으로의 거리를 계산한다.
• 반산형 최적 운반 해법을 통해 라그랑주 세포와 그 중심을 계산하여 입자 위치와 속도를 갱신한다.
• 시간 적분은 심플렉틱 오일러 스킴을 사용하며, 속도 갱신에 라그랑주 세포 이동의 중심을 포함한다.
• 외부 힘을 포함시키기 위해 속도 갱신 식에 힘 항을 추가한다.

실험 결과

연구 질문

• RQ1 최적 운반 기반의 라그랑주 해법이 압축성 없는 유동 방정식의 정규 해로 수렴할 수 있는가?
• RQ2 표준 2차원 유체역학 테스트 케이스(예: 벨트라미 흐름, 불안정성)에서 이 해법의 성능은 어떠한가?
• RQ3 중력과 같은 외부 힘이 도입되었을 때도 해법이 안정적이고 정확한가?
• RQ4 대규모 입자 수(예: 300,000개 입자)를 처리하면서도 기하학적·물리적 일관성을 유지할 수 있는가?
• RQ5 입자 분할 방식(정규 분할 대비 최적화된 분할)의 선택이 시뮬레이션의 정확도와 역학에 어떤 영향을 미치는가?

주요 결과

• 적절한 가정 하에 수치 해법이 압축성 없는 유동 방정식의 정규 해로 수렴함을 보이며, 이론적 수렴성을 입증한다.
• N = 900개 입자를 사용한 단위 정사각형에서의 벨트라미 흐름 실험에서 안정적이고 구조적인 라그랑주 세포를 보이며 정상적인 해를 잘 포착한다.
• N = 300,000개 입자를 사용한 켈빈-헬름홀츠 불안정성 실험에서 최적화된 분할 방식이 대칭성을 깨뜨리고 물리적 기대에 부합하는 소용돌이를 형성한다.
• 중력이 작용하는 N = 50,000개 입자로 구성된 레일리-태일러 불안정성 실험에서 경계면 불안정성이 발생함을 확인하여 다상 유동에서의 해법의 강건성을 입증한다.
• 레일리-태일러 케이스의 계산은 단일 랩탑 코어에서 6시간 이내에 완료되어 계산 효율성을 입증한다.