[논문 리뷰] Hamilton-Pontryagin Integrators on Lie Groups: Introduction and Structure-Preserving Properties

Nawaf Bou‐Rabee, Jerrold E. Marsden|ArXiv.org|2008. 01. 07.

Numerical methods for differential equations참고 문헌 39인용 수 9

한 줄 요약

이 논문은 런게쿠타-밀스타인-카스(RKMK) 방법을 사용하여 해밀턴-포트리아진 변분 원리를 이산화하여 리 군 위에서 새로운 유형의 변분 분할 룬게쿠타(VPRK) 적분기를 도입한다. 이 방법은 심플렉틱 오일러 및 슈타르머-베르트 적분기들을 리 군으로 일반화하며, 운동량 맵과 심플렉틱 구조를 유지하며, 에너지나 운동량을 보존하는 방법들보다 장기적으로 더 뛰어난 안정성과 정확성을 보인다.

ABSTRACT

In this paper structure-preserving time-integrators for rigid body-type mechanical systems are derived from a discrete Hamilton-Pontryagin variational principle. From this principle one can derive a novel class of variational partitioned Runge-Kutta methods on Lie groups. Included among these integrators are generalizations of symplectic Euler and Störmer-Verlet integrators from flat spaces to Lie groups. Because of their variational design, these integrators preserve a discrete momentum map (in the presence of symmetry) and a symplectic form. In a companion paper, we perform a numerical analysis of these methods and report on numerical experiments on the rigid body and chaotic dynamics of an underwater vehicle. The numerics reveal that these variational integrators possess structure-preserving properties that methods designed to preserve momentum (using the coadjoint action of the Lie group) and energy (for example, by projection) lack.

연구 동기 및 목표

• 기하적 성질이 보장되는 리 군 위의 기계 시스템을 위한 기하적 성질 유지를 보장하는 시간 적분기를 개발한다.
• 비平坦한 구성 공간에서 운동학적 제약 조건을 변분 원리를 사용하여 이산화하는 데 도전한다.
• 평탄한 공간에서의 심플렉틱 및 변분 적분기(예: 슈타르머-베르트, 심플렉틱 오일러)를 평탄한 공간에서 리 군으로 일반화한다.
• 변분 설계를 통해 대칭 하에서 운동량 맵과 심플렉틱 구조를 보존한다.
• 강체 및 수중 차량에서의 수치 분석과 실험을 위한 기초를 마련한다(제2부에서 다룸)

제안 방법

• 왼쪽으로 선형화된 해밀턴-포트리아진 작용 원리를 리 군 위에 설정하여, 왼쪽으로 불변인 라그랑지안에 대해 오일러-포이아르 및 리-포아송 역학을 통합한다.
• 군 값 궤적을 다루기 위해 s단계의 런게쿠타-문테-카스(RKMK) 방법을 사용하여 운동학적 제약 조건을 이산화한다.
• 부처 표본 계수를 사용하여 가중 라그랑지안과 운동량 쌍화된 운동학적 제약 조건을 조합한 이산 작용 합을 구성한다.
• 해결된 이산 오일러-라그랑주 방정식의 해로서 변분 분할 룬게쿠타(VPRK) 적분기를 유도한다.
• 운동량 맵의 진화를 표현하고 보존 성질를 검증하기 위해 코어지조인트 작용과 동반 표현을 사용한다.
• 군 소속을 유지하기 위해 패드에 근사식과 프로젝터(예: SO(n)에 대한 반대칭화)를 구현한다.

실험 결과

연구 질문

• RQ1 해밀턴-포트리아진 변분 원리는 리 군 위에서 기하적 성질을 유지하는 적분기를 구성하기 위해 일반화될 수 있는가?
• RQ2 런게쿠타 방법을 사용하여 리 군 위의 운동학적 제약 조건을 기하적 구조를 유지하면서 어떻게 이산화할 수 있는가?
• RQ3 유도된 VPRK 적분기는 심플렉틱성과 운동량 맵을 보존하는가? 그리고 전통적인 에너지 또는 운동량 보존 방법보다 우월한가?
• RQ4 이 적분기의 정확도 순서는 무엇이며, 이는 기초가 되는 변분 구조와 어떻게 관련이 있는가?
• RQ5 비가역 시스템(예: 턴테이블 위의 강체)에서 이 적분기는 대칭적 적분기와 비교해 어떻게 성능을 보이는가?

주요 결과

• 제안된 VPRK 적분기는 대칭 하에서 심플렉틱 형식과 운동량 맵을 모두 보존하여 장기적 안정성을 보장한다.
• 이 방법들은 변분 구조를 통한 탄성 다발 위의 작용 원리를 통해 심플렉틱 오일러 및 슈타르머-베르트 적분기를 리 군으로 일반화한다.
• 이산 작용 합은 부처 표본 계수와 재구성 방정식의 RKMK 이산화를 사용하여 구성된다.
• 운동량 맵의 진화는 이산화된 리-포아송 괄호의 형태로 도출되며, 동반 작용을 통해 보존 성질가 검증된다.
• 이 적분기는 두 번째 차수 정확도를 달성하며, 이는 Marsden & West(2001)의 변분적 정확도 순서 증명을 통해 입증된다.
• 제2부의 수치 실험 결과, 혼돈 동역학과 비가역 시스템에서 에너지 또는 운동량 보존 방법보다 뛰어난 성능을 보였다.