[논문 리뷰] A Note on Symplectic Algorithms

Hao Yang Guo, Y. Q. Li외 1명|ArXiv.org|2001. 04. 06.

Cellular Automata and Applications인용 수 5

한 줄 요약

이 논문은 비가환 미분법과 오일러-라그랑주 코homology를 사용하여 유한차원 해밀턴계에 대한 라그랑주 형식에서 해밀턴형 알고리즘을 제안한다. 삼각형 적분법이 해의 공간뿐 아니라 함수 공간에서도 대칭형을 유지함을 보이며, 유지되는 대칭형 두차형이 닫혀 있고 비퇴화되어 있음을 증명함으로써, 해밀턴형 프레임워크를 초월하여 대칭형 보존을 일반화한다.

ABSTRACT

We present the symplectic algorithm in the Lagrangian formalism for the Hamiltonian systems by virtue of the noncommutative differential calculus with respect to the discrete time and the Euler--Lagrange cohomological concepts. We also show that the trapezoidal integrator is symplectic in certain sense.

연구 동기 및 목표

• 유한차원 해밀턴계에 대해 기존의 해밀턴형 대칭형 보존 기법을 보완하는 라그랑주 형식에서의 대칭형 알고리즘을 수립하기 위해.
• 표준 미분법이 성립하지 않는 비레니츠 법칙을 따르지 않는 시간 도함수로 인해 수치적 적분에서 대칭형을 보존하는 데 도전하는 문제를 해결하기 위해, 해의 공간에서 함수 공간으로 개념을 확장하기 위해.
• 삼각형 적분법이 해의 궤적 위가 아닌 함수 공간에서 일반화된 의미로 대칭형인지를 보이며, 함수 공간에서 유지되는 대칭형 두차형을 유도함으로써 이를 입증하기 위해.
• 이론적 기반을 제공하기 위해 이론적 기반을 제공하기 위해 오일러-라그랑주 코hom로지 개념을 활용하여 이산 시간에서의 대칭형 보존 분석을 체계화하기 위해.

제안 방법

• 비가환 미분법(NCDC)을 이산 시간 격자에 적용하여 레이니츠 법칙을 만족하지 않는 시간 도함수를 다루기 위해.
• 오일러-라그랑주(EL) 일차형과 코hom로지 개념—영 EL 일차형, 코 boundary EL 일차형, EL 조건—을 도입하여 대칭형 보존을 분석하기 위해.
• 기본 일차형 θ = (∂L/∂q̇)dq 에서 유도된 대칭형 두차형 ω = dθ 를 정의하고, EL 조건 dE = 0 하에서 그 시간 도함수가 0이 됨을 보여주기 위해.
• DEL 일차형을 구성하여 삼각형 방법의 유지되는 대칭형 두차형을 도출하고, 그 두차형이 닫혀 있고 비퇴화되어 있음을 보여주기 위해.
• 변수변환 (p̃, q̃) 을 사용하여 유지되는 대칭형 두차형을 dṕ ∧ dq̃ 로 표현함으로써 그 대칭형 성격을 확인하기 위해.
• 분리 가능 및 비분리 가능 해밀턴계 양쪽에 이 방법을 적용하여 각 경우에 대해 유지되는 대칭형 두차형의 명시적 표현을 유도하기 위해.

실험 결과

연구 질문

• RQ1 표준 해밀턴형 프레임워크 외에 라그랑주 형식에서 대칭형 보존을 엄밀하게 정의하고 분석할 수 있는가?
• RQ2 표준 미분법이 비레니츠 법칙을 위반함으로써 실패할 경우, 이산 시간 적분법의 대칭형을 어떻게 분석할 수 있는가?
• RQ3 삼각형 적분법은 해의 궤적 위가 아닌 함수 공간에서 일반화된 의미로 대칭형인가?
• RQ4 오일러-라그랑주 코hom로지가 이산 시간에서 대칭형 적분법을 특성화하는 데 어떤 역할을 하는가?
• RQ5 주어진 수치적 방법에 대해 유지되는 대칭형 두차형을 명시적으로 구성하고, 그것이 닫혀 있고 비퇴화되어 있음을 보일 수 있는가?

주요 결과

• 삼각형 적분법은 해의 궤적 위가 아닌 함수 공간에서도 대칭형을 보존하며, 그 이유는 유지되는 대칭형 두차형 ω_D^{(n+1)} = ω_D^{(n)} 이 닫혀 있고 비퇴화되어 있음을 보장하기 때문이다.
• 분리 가능한 해밀턴함수 H = T(p) + V(q) 에서 유지되는 대칭형 두차형은 ω^{(n+1)} = (dp^{(n+1)})^T (1 + h²/4 T_pp^{(n+1)} V_{qq}^{(n+1)}) dq^{(n+1)} 이다.
• 삼각형 방법의 대칭형은 변환 p̃^{(n)} = p^{(n)} - h/2 H_q^{(n)}, q̃^{(n)} = q^{(n)} + h/2 H_p^{(n)} 하에서 표준 대칭형 dṕ ∧ dq̃ 와 동치이다.
• 삼각형 방법의 유지되는 대칭형 두차형은 명시적으로 ω_D^{(n)} = (dp^{(n)})^T (1 + h²/4 H_pp^{(n)} H_{qq}^{(n)} - h²/4 H_{pq}^{(n)} H_{pq}^{(n)} - H_{pq}^{(n)} H_{pp}^{(n)}) ∧ dq^{(n)} - h²/4 (dq^{(n)})^T H_{qq}^{(n)} H_{pq}^{(n)} ∧ dq^{(n)} 로 주어진다.
• EL 코hom로지 조건 dE = 0 는 대칭형 보존을 위한 필수 및 필요조건이며, 코hom로지와 수치적 대칭형 사이의 연결고리가 된다.
• 이 방법은 해의 공간에서 함수 공간으로 분석을 승격시킴으로써 라그랑주 형식으로 대칭형 보존을 일반화하는데 성공하였으며, 이산 시간 격자 위에서 엄밀한 미분법을 가능하게 한다.