화성 중력장 내 준원형 동결궤도에 대한 해석적 탐구

초록

본 논문은 화성의 비구형 중력장을 고려한 준원형 동결궤도(편심도와 근점인자 ω가 평균적으로 일정한 궤도)를 두 가지 해석적 방법—라그랑주 방정식 기반과 리 변환 기반—으로 분석한다. J₉까지의 조석항을 포함한 80×80 중력모델을 이용해 검증했으며, 네 개의 동결궤도 계열(세 개는 작은 편심도, 하나는 임계궤도 근처)을 발견했다. 대부분의 궤도는 ω=270°(J₃와 J₂ 부호 동일)이며, 저고도·저궤도에서는 ω=90°도 가능함을 제시한다. 또한 특수히 편심도가 거의 0에 가까운 궤도를 찾아 정찰 임무에 유리함을 강조하고, 단일주기 모노드리 매트릭스 트레이스를 이용해 안정성을 평가하였다.

상세 분석

이 연구는 화성 궤도 설계에서 핵심적인 ‘동결궤도’를 정밀히 규명하기 위해 두 가지 전통적인 해석 기법을 적용하였다. 첫 번째는 라그랑주 행성 방정식을 이용해 평균화된 섭동식(특히 J₂, J₃ 등 짝·홀수 차수 조석항)의 장기 변화를 도출하고, 편심도와 근점인자 ω의 평균 변화율을 0으로 만드는 제약조건을 얻는다. 여기서 ω=90° 혹은 270°가 자연스럽게 등장하는데, 이는 J₃가 J₂와 같은 부호를 가져 홀수 차수 항이 ω=270°에서 편심도 감소 효과를 주기 때문이다. 식(11)·(12)를 통해 특정 반경(a≈3597 km)과 궤도경사(i)에서 요구되는 편심도 e를 계산하면, i가 0~90° 구간에 네 개의 해가 존재함을 확인한다. 특히 i≈12° 근처에서 e가 거의 0에 수렴하는 ‘극소 편심도’ 궤도가 나타나며, 이는 정밀 정찰에 이상적이다.

두 번째 방법은 리 변환(Lie transforms)을 이용해 해밀토니안 시스템을 정규화하고, 위상 평면(ω‑e)에서 고정점(동결궤도)의 존재와 안정성을 직접 시각화한다. 정규화된 해밀토니안을 2차까지 보존하면서 등고선(컨투어) 분석을 수행하면, 앞서 라그랑주 방식에서 얻은 네 개의 계열이 동일하게 나타난다. 특히 임계궤도(63.43°) 근처에서는 ω=90°가 자연스럽게 고정점으로 나타나며, 이는 기존 지구·달 연구와 일치한다.

수치 검증은 GMM‑2B 모델(80×80 차수·차수)에서 1년간 적분해 수행했으며, 라그랑주 기반으로 도출한 초기조건( a=3597 km, i=50°, e≈0.00746, ω=270°)에 대해 편심도와 ω가 각각 ±0.0025, ±20° 범위 내에서 진동함을 확인했다. 이는 ‘준동결’ 수준으로 충분히 안정적임을 의미한다. 또한 J₉까지 포함한 조석항만으로도 오차가 수십도 이하로 제한되므로, 고차 항을 무시해도 실용적인 설계가 가능함을 입증한다.

안정성 평가는 단일주기 모노드리 매트릭스의 트레이스를 이용했으며, 0°~90° 범위의 모든 궤도에서 트레이스가 |Tr|<2를 만족, 즉 선형적으로 안정함을 보였다. 다만 고궤도·고경사 구간에서는 비선형 효과가 커질 가능성이 있으므로, 향후 전천후 모델(삼체, 대기, 태양복사압)과 결합한 수치 보정이 필요하다.

이 논문의 주요 공헌은 (1) 화성 특유의 J₂·J₃ 부호 관계가 ω=270°를 우세하게 만든다는 물리적 해석, (2) 저고도·저경사에서 ω=90°가 가능한 새로운 동결궤도 계열 제시, (3) 극소 편심도 궤도의 존재를 밝혀 정찰 임무에 활용 가능성을 제시, (4) 라그랑주와 리 변환 두 방법이 정량적으로 일치함을 검증함으로써 설계 초기조건 제공에 대한 신뢰성을 확보한 점이다. 향후 연구는 비구형 중력 외에 대기·태양압·삼체 효과를 포함한 전천후 시뮬레이션과, 실제 미션에 적용 가능한 연료 최소화 최적화 절차를 개발하는 방향으로 진행될 필요가 있다.