마네프 포텐셜과 대기 항력 하에서 위성 궤도 변동 분석

초록

본 논문은 마네프 중력 포텐셜과 속도 제곱에 비례하는 대기 항력을 동시에 고려한 위성의 원형 궤도 변화를 이론적으로 분석한다. 고도에 따라 지수적으로 감소하는 대기 밀도를 가정하고, 항력 파라미터와 마네프 상수의 영향을 받아 반경 감소율을 도출하였다. 수치 예시를 통해 마네프 항력이 뉴턴 항력 대비 궤도 감쇠에 미치는 차이를 정량화하였다.

상세 분석

이 연구는 전통적인 뉴턴 중력 포텐셜 V_N = −GM/r에 추가적인 1/r² 항을 포함하는 마네프 포텐셜 V_M = −GM/r + β/r² (β는 마네프 상수) 를 채택함으로써, 고전역학과 일반상대론 사이의 중간 형태를 모델링한다. 마네프 항은 특히 근접 궤도에서의 섭동을 설명하는 데 유용하며, 기존 연구에서는 별 궤도 전이, 블랙홀 근처 물체의 궤도, 그리고 원자 물리학에서의 유사 포텐셜로 활용돼 왔다.

논문은 위성이 원형 궤도 r₀ 를 유지한다는 가정 하에, 대기 밀도 ρ(h) = ρ₀ exp(−h/H) (여기서 h 는 고도, H 는 스케일 높이) 로 표현되는 지수 대기 모델을 도입한다. 항력은 F_D = ½ C_D A ρ v² 방향은 속도와 반대이며, 여기서 C_D 는 항력계수, A 는 단면적, v 는 궤도 속도이다. 마네프 포텐셜에 의해 속도는 뉴턴 경우보다 약간 증가하므로, 항력 크기도 변한다.

에너지 보존식에 항력에 의한 손실 dE/dt = −F_D v 를 포함시키고, 마네프 포텐셜에 대한 원심력 균형식 v²/r = GM/r² − 2β/r³ 를 이용해 반경 변화율 dr/dt 를 도출한다. 결과적으로

dr/dt = − (C_D A ρ₀ exp(−h/H) v³) / (2 GM − 4β/r)

와 같은 형태를 얻는다. 여기서 분모의 2GM − 4β/r 는 마네프 항이 궤도 에너지에 미치는 보정이다. β > 0이면 유효 중력이 증가해 항력에 의한 감쇠가 다소 억제되고, β < 0이면 반대로 감쇠가 가속된다.

수치 해석에서는 지구 평균값 GM = 3.986×10¹⁴ m³/s², ρ₀ = 1.225 kg/m³, H ≈ 8.5 km, C_D ≈ 2.2, A ≈ 10 m², β 값을 10⁹ ~ 10¹¹ m³/s² 범위에서 변동시켰다. 결과는 다음과 같다. (1) β가 10⁹ m³/s² 수준일 때, 반경 감소율은 뉴턴 경우 대비 약 5 % 감소한다. (2) β가 10¹¹ m³/s²에 가까워지면 감소율이 30 % 이상 억제된다. 이는 마네프 항이 고도 의존적인 항력 효과와 상호작용해 궤도 수명을 연장시킬 수 있음을 시사한다.

또한, 항력 파라미터 C_D·A·ρ₀ 의 변동에 대한 민감도 분석을 수행했으며, 항력 계수가 10 % 증가하면 dr/dt이 거의 선형적으로 10 % 증가한다. 반면 β의 변화는 비선형적이며, 특히 r이 작아질수록(저궤도) β의 영향이 크게 증폭된다. 이는 저궤도 위성 운영 시 마네프 항을 고려한 궤도 유지 전략이 필요함을 강조한다.

결론적으로, 마네프 포텐셜은 단순한 뉴턴 중력 모델에 비해 항력에 의한 궤도 감쇠를 정량적으로 수정할 수 있는 유용한 도구이며, 특히 저궤도에서의 장기 미션 설계에 적용 가능하다.