[논문 리뷰] Discrete momentum maps for lattice EPDiff

Colin J. Cotter, Darryl D. Holm|arXiv (Cornell University)|2006. 02. 14.

Nonlinear Waves and Solitons참고 문헌 13인용 수 12

한 줄 요약

이 논문은 미분형상의 에일러-포이아르 방정식(EPDiff)에 대한 변분 입자-격자(VPM) 이산화를 개발하여, 라그랑주 입자에 대한 속도장의 왼쪽 및 오른쪽 작용에 대한 이산 운동량 맵을 유도한다. 이 방법은 수치적으로 검증된 이산 킬빈의 정리에 의해 순환을 정확히 유지하며, 운동량 보존을 고려한 비선형 내부파 상호작용을 기하학적으로 일관된 프레임워크로 시뮬레이션할 수 있다.

ABSTRACT

We focus on the spatial discretization produced by the Variational Particle-Mesh (VPM) method for a prototype fluid equation the known as the EPDiff equation}, which is short for Euler-Poincaré equation associated with the diffeomorphism group (of $\mathbb{R}^d$, or of a $d$-dimensional manifold $Ω$). The EPDiff equation admits measure valued solutions, whose dynamics are determined by the momentum maps for the left and right actions of the diffeomorphisms on embedded subspaces of $\mathbb{R}^d$. The discrete VPM analogs of those dynamics are studied here. Our main results are: (i) a variational formulation for the VPM method, expressed in terms of a constrained variational principle principle for the Lagrangian particles, whose velocities are restricted to a distribution $D_{\VPM}$ which is a finite-dimensional subspace of the Lie algebra of vector fields on $Ω$; (ii) a corresponding constrained variational principle on the fixed Eulerian grid which gives a discrete version of the Euler-Poincaré equation; and (iii) discrete versions of the momentum maps for the left and right actions of diffeomorphisms on the space of solutions.

연구 동기 및 목표

• 변분 입자-격자(VPM) 방법을 사용하여 EPDiff 방정식의 기하학적으로 일관된 공간 이산화를 개발한다.
• 반연속 EPDiff에 대한 제약 조건이 붙은 변분 원리를 유도하여 라그랑주 입자 운동과 유레지안 격자 운동을 연결한다.
• 미분형상의 작용에 대한 입자 구성에 대한 이산 운동량 맵을 구축한다.
• VPM 이산화에서 켈빈의 순환 법칙을 유지하는 이산 순환 정리를 수립한다.
• 미래 작업에서 질량 밀도와 같은 운반 물리량을 포함시키기 위해 프레임워크를 확장한다. 이는 분산 효과를 모델링할 수 있도록 한다.

제안 방법

• 영역 위의 벡터장의 유한차원 부분공간으로 입자 속도를 제약 조건을 두어 VPM에 대한 제약 조건이 붙은 변분 원리를 수립한다.
• 변분 원리로부터 유후리안 격자 위의 이산 올레르-포이아르 방정식을 도출하여 연속적인 EPDiff 구조와의 일관성을 확보한다.
• 왼쪽 작용 운동량 맵을 입자 운동량과 형태 함수로부터 표준 VPM 격자 운동량 재구성으로 정의한다.
• 입자 재표기와 관련된 오른쪽 작용 운동량 맵을 도입하여 격자 속도는 유지하면서 운동량을 재할당한다.
• 흐름 지도의 야코비안을 사용하여 순환 고리의 이산 선 요소를 구성하고, 고리 적분의 수치적 근사를 가능하게 한다.
• 이산 변형 행렬을 사용하여 선 요소를 진화시키고 순환 적분을 계산하며, 반올림 오차 이외에는 정확히 보존됨을 보여준다.

실험 결과

연구 질문

• RQ1 VPM 방법은 어떻게 제약 조건이 붙은 변분 원리로 공식화될 수 있으며, 이는 EPDiff의 기하학적 구조를 유지하는가?
• RQ2 미분형상의 입자 구성에 대한 왼쪽 및 오른쪽 작용에 대한 이산 운동량 맵의 이산 형태는 무엇인가?
• RQ3 VPM 프레임워크에서 켈빈의 순환 법칙을 유지하는 이산 순환 정리를 도출할 수 있는가?
• RQ4 이산 운동량 맵은 수치적 방법에서 운동량 보존과 입자 재표기의 물리적 해석과 어떻게 관련되는가?
• RQ5 이산 설정에서 순환 적분을 위한 무한소 선 요소를 진화시키는 데 야코비안의 역할은 무엇인가?

주요 결과

• VPM 방법은 유후리안 격자 위에서 이산 올레르-포이아르 방정식을 만족하여 연속적인 EPDiff 역학과의 일관성을 확보한다.
• 이산 설정에서의 왼쪽 작용 운동량 맵은 표준 VPM 격자 운동량 재구성 공식과 정확히 일치한다.
• 오른쪽 작용 운동량 맵은 격자 속도를 유지하면서 입자 운동량을 재할당하는 이산 재표기 변환과 대응한다.
• 이산 순환 적분은 시뮬레이션 동안 정확히 보존되며, 수치 결과에서 확인된 바에 따르면 변동은 반올림 오차 이외에는 없다.
• 이 방법은 흐름장에서 선 요소의 연신 현상을 성공적으로 포착하며, 시간이 지남에 따라 고리의 이송을 보여주는 그림을 통해 시각화된다.
• 이 프레임워크는 질량 밀도 및 기타 운반 물리량을 포함시키기 위한 VPM 방법의 확장을 위한 기초를 제공하며, 분산 효과의 모델링이 가능하게 한다.