[논문 리뷰] Is Symplectic-Energy-Momentum Integration Well-Posed?

Yosi Shibberu|ArXiv.org|2006. 08. 04.

Numerical methods for differential equations참고 문헌 19인용 수 3

한 줄 요약

이 논문은 확장된 위상공간에서 이산시간 해밀턴(DTH) 방정식을 분석함으로써 심플렉틱-에너지-모멘텀(SEM) 적분기의 잘 정의됨을 조사한다. 비선형 해밀토니안 시스템, 특히 비선형 단자에 대해서도, 위상공간의 특정 점에서 임의로 작은 시간 스텝에서도 DTH 방정식의 해가 존재하지 않을 수 있음을 증명하며, 이는 SEM 적분이 심플렉틱성과 에너지 보존 성질을 갖더라도 항상 잘 정의된다는 가정에 도전한다.

ABSTRACT

We provide new existence and uniqueness results for the discrete-time Hamilton (DTH) equations of a symplectic-energy-momentum (SEM) integrator. In particular, we identify points in extended-phase space where the DTH equations of SEM integration have no solution for arbitrarily small time steps. We use the nonlinear pendulum to illustrate the main ideas.

연구 동기 및 목표

• 심플렉틱-에너지-모멘텀(SEM) 적분이 이산시간 해밀턴(DTH) 방정식의 가역성 분석을 통해 잘 정의되어 있는지 조사하기.
• DTH 방정식이 임의로 작은 시간 스텝에서도 해가 존재하지 않는, 확장된 위상공간 내 특정 점을 특정하기.
• 선형 시스템을 초월해 비선형 해밀토니안 시스템으로까지 DTH 역학의 존재성 및 유일성 결과를 확장하기.
• 변분형식으로 기반한 SEM 적분이 이론적으로 유혹적이지만 해를 생성하지 못하는 조건을 명확히 하기.

제안 방법

• 시간을 일반화된 좌표로 간주하고 그에 해당하는 운동량을 포함한 확장된 위상공간에서 변분 원리를 사용하여 이산시간 해밀토니안(DTH) 역학을 수립한다.
• 시간 스텝의 중간점에서 해밀토니안에 대한 제약 조건을 포함한 정적 작용 원리로부터 DTH 방정식을 유도한다.
• 작용과 에너지 제약 조건과 관련된 함수 $ g(s, z_k) $ 의 행동을 분석함으로써 유도된 비선형 대수방정식의 가역성을 분석한다.
• 중간값 정리와 단조성 논증을 사용하여 시간 스텝 매개변수 $ s_k $ 에 대한 해의 존재성과 유일성을 결정한다.
• 테일러 전개에서 유도된 경계와 리프시츠 유형 추정치를 사용하여 해가 존재하지 않을 수 있는 위상공간 영역을 확립한다.
• 비선형 단자를 구체적인 예로 사용하여, 특정 에너지 및 운동량 조건 하에서 해의 존재가 실패하는 것을 설명한다.

실험 결과

연구 질문

• RQ1 SEM 적분의 DTH 방정식이 임의로 작은 시간 스텝에서도 해가 존재하지 않는 조건은 무엇인가?
• RQ2 확장된 위상공간 내에서는 시간 스텝이 무한소일지라도 DTH 방정식이 잘 정의되지 않는 점이 존재하는가?
• RQ3 DTH 수식에서 에너지 및 운동량 제약 조건이 이산 방정식의 가역성에 어떤 영향을 미치는가?
• RQ4 해의 존재 실패는 $ \mathcal{H}_k $, $ \psi_k $, 및 $ \Lambda_k $ 와 같은 위상공간 변수로 정량화할 수 있는가?
• RQ5 해의 존재 실패는 지에의 정리에 기반한 SEM 적분기의 이론적 정당성과 모순되는가?

주요 결과

• 확장된 위상공간의 특정 점에서 비선형 해밀토니안 시스템에 대해서도, 임의로 작은 시간 스텝에서도 SEM 적분의 DTH 방정식 해가 존재하지 않을 수 있다.
• 비선형 단자에 대해 $ \mathcal{H}_k / \psi_k > \frac{2}{3}(\psi_k / \psi_k')^2 $ 라면, $ (0, S_k) $ 내에서 해가 존재하지 않음을 나타내며, 이는 해가 존재하지 않는 임계값을 의미한다.
• 만약 $ S_k \geq 6 $ 이고 $ \mathcal{H}_k / \psi_k > \frac{2}{3}(\psi_k / \psi_k')^2 $ 라면, $ (0, S_k) $ 내에서 해가 존재하지 않으며, 특정 조건 하에 $ (0, 6/5) $ 와 $ (6/5, S_k) $ 에서 각각 두 개의 해가 공존할 수 있다.
• $ \mathcal{H}_k / \psi_k > \frac{1}{16}\Lambda_k^2(2 - S_k) $ 이고 $ S_k < 6/5 $ 라면, $ (-S_k, 0) $ 에서 해가 존재하지 않음을 나타내며, 이는 해가 존재하지 않는 영역을 의미한다.
• 함수 $ g(s, z_k) $ 는 $ s $ 에 대해 엄격하게 단조적이며, 이는 해가 존재할 경우 유일성을 보장한다. 이는 레마 10 (ii) 및 (iii) 에 의해 입증된다.
• $ \mathcal{H}_k / \psi_k' < -\frac{1}{16}\Lambda_k^3 $ 이면, $ [-\Lambda_k, \Lambda_k] $ 내에서 해가 존재하지 않음을 나타내며, 이는 시간 스텝 매개변수에서의 잘 정의되지 않는 영역을 보여준다.