[논문 리뷰] Symplecticity-preserving continuous-stage Runge-Kutta-Nyström methods

Wensheng Tang, Jing-Jing Zhang|arXiv (Cornell University)|2016. 02. 04.

Numerical methods for differential equations참고 문헌 10인용 수 5

한 줄 요약

이 논문은 분리형 해밀토니안 시스템에서 유래하는 두 번째 차수 상미분방정식을 해결하기 위해 심플렉틱성 보존 연속단계 룬게쿠타-니스트롬(csRKN) 방법을 제안한다. 레지온드르 다항식 전개를 활용하고 충분한 심플렉틱성 조건을 유도함으로써, 고차수, 대각선 음이 아닌 및 양이 아닌 심플렉틱 적분기법을 구축한다. 이는 기존 방법이 복잡한 비선형 순서 조건을 풀어야 하는 데 비해, 해밀토니안 역학의 기하적 구조를 더 효율적으로 유지한다.

ABSTRACT

We develop continuous-stage Runge-Kutta-Nyström (csRKN) methods for solving second order ordinary differential equations (ODEs) in this paper. The second order ODEs are commonly encountered in various fields and some of them can be reduced to the first order ODEs with the form of separable Hamiltonian systems. The symplecticity-preserving numerical algorithm is of interest for solving such special systems. We present a sufficient condition for a csRKN method to be symplecticity-preserving, and by using Legendre polynomial expansion we show a simple way to construct such symplectic RKN type method.

연구 동기 및 목표

• 분리형 해밀토니안 시스템에서 유래하는 두 번째 차수 상미분방정식의 심플렉틱 구조를 보존하는 연속단계 RKN 방법을 개발한다.
• csRKN 방법에서 심플렉틱성의 충분 조건을 설정하여 기하학적 보존 적분기법의 체계적 구성 가능성을 확보한다.
• 기존 방법에서 흔히 발생하는 복잡한 비선형 순서 조건을 풀 필요 없이 고차수 심플렉틱 적분기법의 구성 방식을 단순화한다.
• 동일한 연속단계 공식에서 다양한 적분 규칙을 사용함으로써 여러 심플렉틱 적분기법을 유도할 수 있음을 보여준다.
• 자유 매개변수를 가진 4차의 명시적 및 대각선 음이 아닌 심플렉틱 적분기법을 제공하며, Störmer-Verlet과 같은 기존 방법도 포함한다.

제안 방법

• 내부 단계의 연속단계 공식과 부처 타블로 함수에 대한 의존성 분석을 통해 csRKN 방법에서 심플렉틱성의 충분 조건을 도출한다.
• 계수 함수 $ A_{\tau,\sigma} $, $ \hat{A}_{\tau,\sigma} $, $ B_\tau $, 및 $ \hat{B}_\tau $ 를 레지온드르 다항식 전개로 매개변수화하여, 심플렉틱 방법의 해석적 구성이 가능하도록 한다.
• 확장된 계수에 대해 심플렉틱성 및 순서 조건을 적용함으로써, 자유 매개변수를 결정하기 위한 대수방정식계를 도출한다.
• 가우스 적분 규칙을 사용하여 연속단계 방법을 고전적 RKN 스킴으로 이산화함으로써, 명시적 및 대각선 음이 아닌 적분기법을 도출한다.
• 계수에 특정 대수적 제약 조건을 도입함으로써, 이중 매개변수 체계의 대각선 음이 아닌 심플렉틱 적분기법을 구성한다.
• 계수 행렬 $ \bar{A}_{ij} $ 에 대해 추가로 상삼각형 조건(하삼각형 요소가 0)을 도입함으로써, csRKN 방법에서 명시적 심플렉틱 적분기법을 도출할 수 있음을 보여준다.

실험 결과

연구 질문

• RQ1 분리형 해밀토니안 시스템에서 유도되는 연속단계 RKN 방법을 통해 심플렉틱 구조를 보존할 수 있는가?
• RQ2 csRKN 방법이 심플렉틱임을 보장하는 충분 조건은 무엇이며, 이를 계수 함수의 형태로 어떻게 표현할 수 있는가?
• RQ3 레지온드르 다항식 전개를 어떻게 활용하여 자유 매개변수를 가진 체계적인 csRKN 심플렉틱 방법을 구성할 수 있는가?
• RQ4 동일한 csRKN 공식에서 다양한 적분 규칙을 사용하여 다른 고전적 RKN 스킴을 도출할 수 있는가?
• RQ5 적분 규칙의 순서와 결과로 도출되는 이산 RKN 방법의 순서 사이의 관계는 무엇인가?

주요 결과

• 심플렉틱성의 충분 조건이 도출되어, 복잡한 비선형 순서 조건 전체를 풀지 않고도 심플렉틱 적분기법을 구성할 수 있다.
• 자유 매개변수를 가진 4차 대각선 음이 아닌 심플렉틱 적분기법의 이중 매개변수 체계가 구성되었으며, 표 3 및 식 (5.6)–(5.8)에 명시된 형태로 제공된다.
• 계수 행렬에 추가 제약 조건을 도입함으로써 명시적 심플렉틱 적분기법이 도출되었으며, 이로써 새로운 3개의 명시적 스킴이 유도되었고, 그 중 하나는 잘 알려진 Störmer-Verlet 방법이다.
• 특수한 값 $ \alpha = \frac{1}{8}, \beta = \frac{\sqrt{3}}{24}, \gamma = -\frac{1}{8} $ 에서 Störmer-Verlet 스킴이 복원되며, 이는 기존 결과와의 일致성을 확인한다.
• 고차수 적분기법(예: 3단계, 4차)은 고차수 적분 규칙으로 확장하고 결과로 생기는 대수방정식계를 푸름으로써 구성할 수 있다.
• 기존 방법이 순서 조건에서 유도되는 복잡한 비선형 대수계를 풀어야 하는 데 비해, 이 접근 방식은 고차수 심플렉틱 RKN 방법의 구성 과정을 단순화한다.