리오믹 비홀로닉 변형으로 보는 Euler 방정식

초록

1913년 A.D. 빌리모비치가 제시한 시간 의존(리오믹) 선형·동차 비홀로닉 제약이 이상적임을 이용해, Euler 방정식에 비홀로닉 제약을 가한 시스템을 분석한다. 제약 f = ω₁ − g(t) ω₂ = 0을 적용하면 에너지 보존을 유지하면서도 운동 방정식이 비자율 형태로 변하고, 이를 적분가능한 1차 미분식으로 환원한다. 특수 경우(상수 g, 대칭 관성 텐서 등)에는 해가 타원함수 형태로 명시적으로 구해지며, 일반 경우에도 수치적·기하학적 구조(불변 측도, 해밀토니안 구조)와 연계된다.

상세 분석

본 논문은 빌리모비치가 1913년에 제시한 “시간 의존 선형·동차 제약은 이상적이다”는 명제를 현대적인 비홀로닉 역학 체계에 적용한다는 점에서 학문적 의미가 크다. 먼저 일반적인 라그랑주계에 비홀로닉 제약 f_j(q, ẋ, t)=0을 도입하고, Chetaev 조건을 통해 반작용력의 일은 영이 됨을 재확인한다. 특히 제약이 ẋ에 선형·동차일 경우, 반작용력 λ ∂f/∂ẋ가 에너지 변화식에 기여하지 않아 전체 기계적 에너지 H가 보존된다. 이 성질을 Euler 방정식 I · ω̇ = I ω × ω에 적용하면, 제약 f = ω₁ − g(t) ω₂ = 0을 통해 비자율 비홀로닉 시스템을 만든다.

제약을 미분하여 λ를 제거하고, 에너지 보존식 T=E와 제약식 f=0을 동시에 풀어 ω₁, ω₂를 ω₃와 시간 함수 g(t), ġ(t)로 표현한다. 결과적으로 ω₃에 대한 1차 비자율 미분식 (11), (12)가 도출되며, 이는 “분리 가능한 변수” 형태이므로 적분은 이론적으로 가능하다. 식 (12)는 제곱근을 포함한 형태이지만, 관성 텐서가 대각화되거나 오프다이아고날 성분이 특정 조건을 만족하면 √항이 사라져 보다 간단한 형태(13)·(14)로 축소된다.

특히 g(t)=α(상수)인 경우, 시스템은 두 개의 독립적 적분 T와 α=ω₁/ω₂를 보유하고, 차원 감소를 통해 ω₃에 대한 단일 ODE로 환원된다. 이때 불변 측도 μ=ρ dω₁dω₂dω₃ (ρ=1/(ω₁ω₂))와 해밀토니안 구조가 존재함을 보이며, 포아송 구조 P와 해밀토니안 T가 X=P dT 형태로 표현된다. 이는 비홀로닉 시스템이지만 특정 대칭 하에서는 사실상 해밀토니안 흐름임을 의미한다.

수치 실험에서는 관성 텐서의 비대각 성분(I₁₂, I₁₃, I₂₃)의 유무가 해의 정성적 형태에 큰 영향을 미침을 확인한다. 대각 텐서(오프다이아고날 성분 0)에서는 ω₃(t)가 주기적·준주기적 패턴을 보이며, g(t)=cos t와 같은 주기적 제약은 궤적을 복합적인 리프 형태로 만든다. 반면 오프다이아고날 성분이 존재하면 궤적이 비대칭적으로 변형되고, 일부 파라미터 구간에서는 해가 실수 영역을 벗어나거나 A>0 조건이 깨져 수치적 불안정이 발생한다.

결론적으로, 빌리모비치의 리오믹 비홀로닉 제약은 에너지 보존을 유지하면서도 시스템을 적분가능하게 만들 수 있는 강력한 도구이며, 관성 텐서와 제약 함수 g(t)의 선택에 따라 해밀토니안 구조, 불변 측도, 그리고 해의 정성적 거동이 크게 달라진다. 향후 연구는 비동질(비이상적) 제약 f=ω₁−g(t)ω₂=a와 같은 경우에 대한 에너지 비보존 현상과 부분적 적분 가능성, 그리고 이러한 시스템을 물리적으로 구현할 수 있는 메커니즘을 탐구하는 방향으로 진행될 수 있다.