힐 방정식

1. 개요[편집]

힐 방정식 Hill's equations (of relative motion)
다른 이름	힐-클로히시-윌트셔 방정식 (HCW), CW 방정식
유래	힐(1878) 달 이론 · 클로히시-윌트셔(1960) 랑데부
좌표계	LVLH — x 반경 · y 진행 · z 궤도면 법선
가정	기준 궤도가 원 · 상대거리 ≪ 궤도반경 · 이체 중력
해	닫힌 형태 (6×6 상태천이행렬)
대표 해	2×1 타원 (풋볼 궤도) + along-track 표류
동명이인	플로케 이론의 y″ + p(t)y = 0 — 다른 방정식

힐 방정식은 천체역학에서 원궤도를 도는 기준 천체에 고정된 회전좌표계에서, 그 근처를 도는 다른 물체의 상대운동을 기술하는 선형 상미분방정식계다. 힐-클로히시-윌트셔(Hill–Clohessy–Wiltshire) 방정식, 줄여서 HCW 방정식 또는 CW 방정식이라고도 한다. 우주 랑데부, 도킹, 근접운용, 위성 편대비행 설계가 전부 이 세 줄 위에서 이루어진다.

$$\ddot x - 3n^2 x - 2n\dot y = 0, \qquad \ddot y + 2n\dot x = 0, \qquad \ddot z + n^2 z = 0$$

시작하기 전에 이름부터 정리해야 한다. “Hill’s equation”이라는 말은 이 분야에서 서로 다른 두 대상을 가리킨다.

• 이 문서가 다루는 것 — 위의 상대운동 방정식계. 조지 윌리엄 힐이 1878년 Researches in the Lunar Theory 에서 세운 힐 문제(제한 삼체 문제의 국소 극한)를 선형화한 것이다.
• 플로케 이론의 힐 방정식 — 주기계수 2계 스칼라 방정식 $y'' + p(t)\,y = 0$, $p(t+T)=p(t)$. 힐이 1877년 달의 근지점 운동을 다루며 무한행렬식으로 공략한 그 방정식이고, $p(t)=a-2q\cos2t$ 인 가장 단순한 비자명 경우가 마티외 방정식이다. 안정성 선도와 무한행렬식 이야기는 마티외 방정식과 플로케 이론에 있으니 그쪽으로 가면 된다.

둘은 같은 사람이 만들었지만 다른 방정식이다. 이름이 겹치는 것은 힐의 달 이론이 워낙 여러 갈래로 뻗은 탓이고, 문헌에서 “Hill equation”을 만나면 문맥부터 확인해야 한다. 한쪽은 파라메트릭 여기와 불안정 혀 이야기, 다른 쪽은 도킹 이야기다.

기준 궤도 자체의 이체 문제와 “따라잡으려면 감속해야 한다”는 유명한 반직관은 궤도 역학에 이미 정리돼 있으므로, 이 문서는 방정식의 유도, 닫힌해와 상태천이행렬, 상대궤도의 기하, 그리고 가정이 깨질 때의 보정을 다룬다.

2. 힐 문제 — 제한 삼체의 국소 극한[편집]

HCW의 출신을 알면 그 한계도 같이 보인다. 제한 삼체 문제에서 질량비 $\mu = m_2/(m_1+m_2)$ 를 0으로 보내되, 작은 천체 근방을 $\mu^{1/3}$ 로 확대해서 보면 극한 방정식이 남는다. 회전좌표계에서 시간 단위를 $n=1$ 로 잡고 $\rho=\lVert(\xi,\eta,\zeta)\rVert$ 라 하면

$$\ddot\xi - 2\dot\eta = 3\xi - \frac{\xi}{\rho^3},\qquad \ddot\eta + 2\dot\xi = -\frac{\eta}{\rho^3},\qquad \ddot\zeta = -\zeta - \frac{\zeta}{\rho^3}$$

이것이 힐 문제다. 질량비가 사라졌다는 점이 핵심이다 — 지구-달이든 태양-목성이든 토성-타이탄이든, 근방의 동역학은 하나의 매개변수 없는 방정식으로 통일된다. 야코비 상수에 대응하는 보존량이 살아남고, 그 영속도면이 닫히는 거리가 힐 구의 반지름

$$r_H \;=\; a\left(\frac{m}{3M}\right)^{1/3}$$

이며, 힐 단위로는 공선 평형점이 $\xi=\pm3^{-1/3}$ 에 온다. 위성 포획, 불규칙 위성의 안정 영역, 원시행성의 물질 공급 반경 판정이 전부 이 방정식의 결과물이다.

이제 $1/\rho^3$ 항을 지우면 — 즉 따라다니는 물체가 서로 끌어당길 만큼 무겁지 않다고 하면 — 정확히 HCW가 남는다. 인공위성 두 대의 상대운동에서는 상호 중력이 완전히 무시할 만하므로 이 생략이 정당하다. 반대로 소행성 주변에서 두 조각의 상대운동을 볼 때는 지우면 안 된다.

3. 방정식의 세 항이 어디서 오는가[편집]

LVLH(Local-Vertical/Local-Horizontal) 좌표계를 기준 위성(chief)에 붙인다. $\hat{\mathbf x}$ 는 지구 중심에서 바깥으로 향하는 반경 방향, $\hat{\mathbf y}$ 는 진행 방향, $\hat{\mathbf z}=\hat{\mathbf x}\times\hat{\mathbf y}$ 는 궤도면 법선(각운동량 방향)이다.¹ 이 좌표계는 각속도 $n=\sqrt{\mu/a^3}$ 로 회전한다.

회전좌표계의 겉보기힘과 중력의 1차 전개를 합치면 계수가 그대로 설명된다.

• 반경 방향 $3n^2x$. 원심력이 $n^2 x$, 중력 $-\mu/r^2$ 를 $r=a$ 에서 전개한 조석 항이 $+2n^2x$. 둘을 더해 $3n^2x$ 다. 부호가 양이라는 것이 중요하다 — 반경 방향은 불안정하다.
• 면외 방향 $-n^2z$. 중력 조석 항 $-\mu z/a^3 = -n^2z$ 만 있고 원심력 성분이 없다. 그래서 면외 운동은 기준 궤도와 정확히 같은 주기의 단순 조화진동이다.
• 진행 방향은 조석 항이 0. 1차 전개에서 진행 방향 중력 기울기가 사라지므로 코리올리 항 $\pm2n\dot{(\cdot)}$ 만 남는다. $y$ 방향에 복원력이 전혀 없다는 사실이 아래에서 볼 표류의 근원이다.

계를 1계로 낮추면 $\dot{\mathbf s}=A\mathbf s$ 꼴의 6차원 선형 시불변계가 되고, $A$ 의 고유값은 $\pm in$ 이 두 쌍, 그리고 0이 두 겹이다. 0 고유값의 조르당 블록(대각화 불가)이 바로 진행 방향의 세속 표류이며, 이 계가 지수 이분적이지 않은 이유이기도 하다.

4. 닫힌해와 상태천이행렬[편집]

선형 시불변계이므로 행렬 지수함수로 정확히 풀린다. 초기 상태 $(x_0,y_0,z_0,\dot x_0,\dot y_0,\dot z_0)$ 에서

$$x(t) = (4-3\cos nt)\,x_0 + \frac{\sin nt}{n}\,\dot x_0 + \frac{2(1-\cos nt)}{n}\,\dot y_0$$ $$y(t) = 6(\sin nt - nt)\,x_0 + y_0 - \frac{2(1-\cos nt)}{n}\,\dot x_0 + \frac{4\sin nt - 3nt}{n}\,\dot y_0$$ $$z(t) = z_0\cos nt + \frac{\dot z_0}{n}\sin nt$$

속도는 이를 미분하면 나온다. 이 여섯 줄을 모으면 6×6 상태천이행렬 $\Phi(t)$ 이고, 면외 성분이 완전히 분리되어 $4\times4$ 와 $2\times2$ 로 블록 대각이다.

$\Phi$ 를 위치·속도 블록으로 쪼개 쓰면 랑데부 유도가 곧바로 나온다.

$$\begin{pmatrix}\mathbf r(t)\\ \dot{\mathbf r}(t)\end{pmatrix} =\begin{pmatrix}\Phi_{rr} & \Phi_{rv}\\ \Phi_{vr} & \Phi_{vv}\end{pmatrix} \begin{pmatrix}\mathbf r_0\\ \dot{\mathbf r}_0\end{pmatrix} \;\Longrightarrow\; \dot{\mathbf r}_0^{\,+} = \Phi_{rv}^{-1}\bigl(\mathbf r_f - \Phi_{rr}\mathbf r_0\bigr)$$

시각 $0$ 에 속도를 $\dot{\mathbf r}_0^{\,+}$ 로 바꾸는 임펄스 하나, 시각 $t$ 에 상대속도를 0으로 죽이는 임펄스 하나 — 이것이 2임펄스 랑데부의 전부이며, 총 $\Delta v$ 는 두 임펄스 크기의 합이다. 전이시간 $t$ 를 훑으며 $\Delta v$ 를 그리면 최적 전이시간이 나온다.

함정이 하나 있다. $\Phi_{rv}$ 의 면내 $2\times2$ 부분의 행렬식은

$$\det = \frac{\sin nt\,(4\sin nt - 3nt) - 4(1-\cos nt)^2}{n^2}$$

인데, $nt = 2\pi k$ 에서 정확히 0이 된다. 즉 전이시간이 궤도주기의 정수배면 목표점을 맞출 수 없다 — 정확히 한 바퀴 뒤에는 어떤 초기속도를 줘도 면내 상대위치가 제자리로 돌아온다. 랑데부 코드가 이 근처에서 $\Delta v$ 가 발산하며 죽는 것은 버그가 아니라 물리다. 실무에서는 전이시간을 주기의 정수배에서 몇 % 떨어뜨리거나 다임펄스로 쪼개 피해 간다.

5. 상대궤도의 모양 — 2×1 타원과 표류[편집]

$y(t)$ 에 시간에 비례하는 항이 두 개 들어 있다. 모아 쓰면 세속 표류율은

$$\dot y_{\rm sec} = -\bigl(6n x_0 + 3\dot y_0\bigr) = -3\bigl(\dot y_0 + 2n x_0\bigr)$$

이다. 표류가 없을 조건은 따라서 $\dot y_0 = -2n x_0$ 하나뿐이다. 이 조건을 넣으면 위 해가 놀랍도록 깔끔해진다.

$$x(t) = A\cos(nt-\phi), \qquad y(t) - y_c = -2A\sin(nt-\phi), \qquad z(t) = B\cos(nt-\psi)$$

면내 운동이 반경 방향 반지름 $A$, 진행 방향 반지름 $2A$ 인 타원을 그린다. 장축이 진행 방향이고 축비가 정확히 2:1이라 2×1 타원 또는 생김새 때문에 풋볼 궤도라 부른다. 회전 방향은 궤도운동과 반대다 — 바깥쪽(높은 곳)에 있을 때 뒤로 처지고 안쪽에 있을 때 앞으로 나가는, 궤도 역학의 그 반직관이 타원 한 바퀴로 눈에 보이는 형태다.

표류의 정체는 장반경 차이다. 위 표류율을 궤도 에너지로 바꿔 쓰면 정확히

$$\dot y_{\rm sec} = -\frac{3}{2}\,n\,\delta a$$

가 나온다($\delta a$ 는 두 위성의 궤도 장반경 차이). 케플러 제3법칙 $T\propto a^{3/2}$ 에서 $\delta T/T = \tfrac32\,\delta a/a$ 이므로, 장반경이 다르면 주기가 다르고, 주기가 다르면 매 바퀴마다 진행 방향으로 밀린다. 표류 없는 상대궤도 = 장반경이 정확히 같은 궤도. 이심률이나 궤도면이 아무리 달라도 장반경만 같으면 두 위성은 영원히 서로를 맴돈다. 편대비행 설계의 제1계명이 이것이고, 반대로 GRACE처럼 일부러 앞뒤로 늘어서고 싶으면 장반경을 살짝 어긋내면 된다.

면외 진동 주기가 면내와 정확히 같다는 점도 설계에 그대로 쓰인다. $B=2A$ 로 잡고 위상을 $\psi = \phi \pm \pi/2$ 로 맞추면 진행-법선 평면에서 상대궤도가 반지름 $2A$ 의 원이 된다(투영 원궤도, PCO). 여기에 위상을 다르게 주면 두 위성이 절대 같은 위치를 지나지 않는 나선 편대(TanDEM-X가 쓰는 helix)가 나오고, 이 “지나가지 않음”이 곧 수동적 충돌 안전성이 된다. 안전 타원(safety ellipse)이라는 이름으로 근접운용에서 표준으로 쓰인다.

6. 랑데부와 근접운용[편집]

클로히시와 윌트셔가 1960년에 이 방정식을 다시 세운 것은 순전히 말단 랑데부 유도를 위해서였고, 5년 뒤 제미니 6A호와 7호가 사상 최초의 궤도상 랑데부에 성공했다. 지금도 접근 유도의 기본 골격은 같다.

• 접근 축. 진행 방향을 따라 다가가는 V-bar 접근과 아래에서 반경 방향으로 올라가는 R-bar 접근이 있다. R-bar 접근이 선호되는 국면이 있는데, HCW 동역학 자체가 접근을 자연 감속시켜 주기 때문이다. 추력 실패 시에도 그냥 멈춰 서는 쪽이 안전하다.
• 유도. 2임펄스 targeting으로 대략적인 궤적을 잡고, 그 위에 글라이드슬로프(거리에 비례해 접근속도를 줄이는 규칙)를 얹은 뒤, 최종 단계는 모델 예측 제어나 리카티 방정식 기반 LQR로 제약(추력 한계, 접근 원뿔, 플룸 회피)을 지키며 수렴시킨다. HCW가 선형이라는 사실이 여기서 값어치를 한다 — 예측 모델이 정확한 닫힌 형태라 최적 제어 문제가 볼록한 이차계획으로 떨어진다.
• 항법. 상대 위치·속도는 관측이 아니라 추정이고, 상대 GPS·라이다·비전 관측을 칼만 필터로 융합한다. 필터의 전파 단계에서 쓰는 것이 바로 위의 $\Phi(t)$ 다.
• 수동 안전성. 추력이 끊겨도 며칠 동안 충돌하지 않는 상대궤도에 머무는 것이 유인 우주선 접근의 규정이며, 이것을 설계하는 언어가 앞 절의 타원 기하다.

7. 가정이 깨지는 지점과 그 보정[편집]

HCW는 원궤도 + 선형화 + 이체 중력이라는 세 가정 위에 서 있고, 실전에서는 셋 다 깨진다.

1. 이심률. 기준 궤도가 타원이면 각속도가 일정하지 않아 계수가 시간에 의존하게 된다. 표준 처방은 독립변수를 시간에서 진근점이각으로 바꾸는 차우너-헴펠(Tschauner–Hempel, TH) 방정식이다(1965). TH는 여전히 선형이지만 주기계수 계라 플로케 이론의 대상이 되고, 오래도록 닫힌 형태 상태천이행렬이 없어 수치적분에 의존해야 했다. 야마나카와 안케르센(2002)이 특이점 없는 간결한 닫힌 형태 STM을 내놓으면서 이심 궤도 편대비행 설계가 실용 궤도에 올랐다. $e\to0$ 극한에서 HCW로 정확히 환원된다.

2. $J2$ 섭동. 지구의 편평도는 승교점 적경과 근지점 인수를 세속 세차시키는데, 그 세차율이 장반경·이심률·경사각에 의존하므로 두 위성이 조금만 달라도 상대궤도가 서서히 뒤틀린다. 특히 면외 성분의 표류가 커서, HCW의 “면외는 얌전한 조화진동”이라는 그림이 며칠 만에 무너진다. 대응은 두 갈래다.

• $J2$ 불변 상대궤도 — 슈아브와 알프리엔드가 정리한 방식으로, 세 개의 평균 궤도요소 세차율(장반경 드리프트·근지점·승교점)이 서로 같아지도록 요소 차이에 제약을 건다. 완벽한 동시 만족은 일반적으로 불가능해서 어느 항을 희생할지 고르는 설계 문제가 된다.
• $J2$ 를 포함한 상태천이행렬 — 김-알프리엔드 STM처럼 평균 궤도요소 공간에서 세차를 해석적으로 넣은 전파 행렬을 쓴다.

3. 비선형성. 선형화에서 버린 항은 상대거리 $\rho$ 에 대해 $O(\rho/a)$ 규모라, 저궤도($a\approx7000\ \mathrm{km}$)에서 1 km 떨어진 두 위성이면 상대오차가 $10^{-4}$ 수준으로 시작하지만 궤도 수에 비례해 누적된다. 수 km를 넘거나 수십 바퀴를 예측해야 하면 HCW는 초기 추정 이상으로 쓸 수 없다.

4. 항력 차이. 저궤도에서 실제로 편대를 무너뜨리는 1순위는 사실 $J2$ 도 비선형성도 아니고 탄도계수 차이에 의한 차등 항력이다. 장반경이 서로 다른 속도로 줄어들면 앞의 공식이 그대로 표류를 낳는다. 자세를 바꿔 단면적을 조절하는 차등 항력 제어가 추진제 없는 편대 유지 수단으로 실제 운용된다.

5. 그래서 현대의 표준. 오늘날 편대비행 설계는 데카르트 상대 상태 대신 상대 궤도요소(ROE)를 상태변수로 쓰는 쪽이 주류다. 준-비특이 ROE는 $J2$ 세차와 차등 항력이 요소 공간에서 거의 선형·해석적으로 표현되고, 각 요소가 “표류”, “이심률 벡터 분리” 같은 기하적 의미를 그대로 갖는다는 것이 강점이다. PRISMA와 TanDEM-X의 편대 유지가 이 방식으로 운용됐다. HCW는 죽은 것이 아니라 ROE 표현의 원궤도·선형 극한으로 자리를 옮겼다.²

8. 어디서 실물로 쓰이나[편집]

• 유인 우주선 도킹. ISS 방문 우주선(드래곤, 시그너스, 프로그레스, HTV)의 접근 회랑과 중단 규정이 전부 HCW 기하로 정의된다.
• 과학 편대. GRACE 계열의 진행 방향 편대(장반경 차이를 의도적으로 유지), TanDEM-X의 나선 편대(면외 위상으로 충돌 회피), MMS의 사면체 편대(네 위성의 상대 기하 자체가 관측 장비).
• 궤도상 서비싱. 정지궤도 위성에 도킹해 수명을 연장하는 임무들이 접근 최종 단계에서 이 방정식을 쓴다. 정지궤도는 $n$ 이 작아 표류가 느리고 상대 동역학이 온순한 편이라 근접운용에 유리하다.
• 잔해 제거·검사. 비협조 목표(회전하는 잔해)에 접근할 때는 HCW로 상대궤도를 설계하고, 목표의 자세 운동은 별개의 강체 문제로 붙인다.

한 줄로 요약하면 이렇다. HCW는 정확한 모델이 아니라 정확한 언어다. 실제 궤적은 N체 문제 수준의 고정밀 모델과 심플렉틱 적분기로 검증하지만, 설계자가 “여기에 2×1 타원을 하나 놓고, 표류를 죽이려면 장반경을 맞추고, 안전을 위해 면외 위상을 90도 준다”고 말할 때 쓰는 언어가 이 세 줄이다. 60년이 지나도록 자리를 지키는 이유가 정확도가 아니라 말이 되는 좌표계라는 데 있다.³

9. 관련 문서[편집]

• 제한 삼체 문제 · 궤도 역학 · N체 문제 · 힐 구
• 마티외 방정식 · 플로케 이론 — 동명이인 쪽 힐 방정식
• 랑데부 · 편대비행 · 상태천이행렬
• 최적 제어 · 모델 예측 제어 · 리카티 방정식 · 칼만 필터
• 심플렉틱 적분기 · 행렬 지수함수 · 전기추진

10. Footnotes[편집]

1. 좌표 관례가 문헌마다 다르니 남의 코드를 가져다 쓸 때 반드시 확인해야 한다. 이 문서의 “x 반경, y 진행, z 법선” 말고도 “x 진행, y 법선 반대, z 천저” 같은 배열이 항공우주 실무에서 흔히 쓰인다. 부호 하나 틀리면 감속해야 할 때 가속하는 유도 법칙이 나오고, 그 결과는 시뮬레이터에서 목표를 향해 우아하게 멀어지는 그래프로 나타난다. ↩

2. “1960년 논문의 방정식을 아직도 쓴다”는 말은 반은 맞고 반은 틀리다. 개념 설명과 초기 설계는 HCW로 하지만, 실제 운용 코드에서 상태를 전파하는 것은 ROE + 섭동 모델이다. 학부 과제로 HCW를 짜 본 사람이 현업에 와서 “왜 여기엔 CW가 없죠?”라고 묻는 것이 이 바닥의 통과의례. ↩

3. 힐이 1878년 논문을 쓸 때 상상한 응용은 달의 운동 이론이었지, 80여 년 뒤 사람이 탄 캡슐 두 개를 궤도에서 붙이는 일이 아니었다. 순수 천체역학에서 나온 극한 문제가 그대로 유인 우주비행의 설계 도구가 된 사례로 자주 인용된다. 물론 클로히시와 윌트셔는 힐의 결과를 참고한 게 아니라 필요해서 다시 유도했고, 이름이 셋 다 붙은 것은 나중 일이다. ↩