3단계: 궤도 역학, 다단 로켓, 재사용 착륙 제어

이론적 기초 — 우주에서 "던진다"는 것의 의미

7살 아이에게 우주선이 왜 떨어지지 않냐고 물으면 아마 "엔진이 계속 켜져 있어서요"라고 답할 것이다. 틀렸다. 국제우주정거장(ISS)은 엔진을 끈 채로 지구 주위를 돌고 있다. 그렇다면 왜 떨어지지 않는 걸까? 정확히 말하면, ISS는 계속 떨어지고 있다. 다만 지구가 충분히 둥글기 때문에, ISS가 떨어지는 속도만큼 지구 표면도 휘어져 "도망가는" 것이다. 이 기묘한 자유낙하 상태가 바로 **궤도(Orbit)**의 본질이다. 뉴턴이 사과를 보고 깨달은 것도 이것이었다 — "달도 사과처럼 떨어지고 있다. 단지 옆으로 너무 빠르게 던져진 것뿐이다."

이 직관을 수식으로 단단히 고정시켜야 한다. 뉴턴의 만유인력 법칙은 F = GMm/r² 로 쓰이는데, 여기서 G는 만유인력 상수(6.674 × 10⁻¹¹ N·m²/kg²), M은 중심 천체(예: 지구)의 질량, m은 위성의 질량, r은 중심 간 거리다. 이 힘이 원심력과 균형을 이룰 때 원형 궤도가 만들어진다. 원형 궤도에서 GMm/r² = mv²/r 이므로, 궤도 속도는 v = √(GM/r)이다. 지구 저궤도(고도 약 400km)에서는 약 7.67 km/s — 총알보다 약 10배 빠르다. 이 수식이 보여주는 가장 중요한 사실은 궤도 속도는 위성의 질량에 무관하다는 점이다. 1g짜리 나사못이든 100톤짜리 우주정거장이든 같은 고도에서는 같은 속도로 궤도를 돈다.

[노트 기록] 원형 궤도 속도 유도: F_중력 = F_원심 → GMm/r² = mv²/r → v_c = √(GM/r). 지구: GM = μ = 3.986 × 10¹⁴ m³/s². 이 값(μ)을 **표준 중력 매개변수(Standard Gravitational Parameter)**라 부르며, GM을 별도로 나누지 않고 하나로 묶어 계산 오차를 줄인다.

케플러 법칙은 이 궤도의 기하학적 패턴을 정리한 것이다. 케플러 제1법칙은 행성의 궤도가 타원이며 태양이 한 초점에 위치한다고 말한다. 완벽한 원형 궤도는 이 타원의 특수한 경우(이심률 e = 0)다. 케플러 제3법칙은 훨씬 중요한데, 궤도 주기의 제곱이 반장축(semi-major axis, a)의 세제곱에 비례한다: T² ∝ a³, 정확히는 T = 2π√(a³/μ). 이 법칙의 의미를 음미해보라 — 궤도가 클수록(a가 클수록) 주기가 길고, 따라서 궤도 속도가 느리다. 즉 높은 궤도는 느리고, 낮은 궤도는 빠르다. 이 역설처럼 보이는 사실이 호만 전이를 이해하는 핵심 열쇠다.

에너지 관점에서 보자. 궤도를 도는 물체의 총 역학적 에너지는 운동에너지와 위치에너지의 합인데, ε = v²/2 - μ/r로 쓰인다. 원형 궤도에 v_c = √(μ/r)을 대입하면 ε = -μ/(2a) 가 나온다 (원형 궤도에서 a = r). 이 에너지는 음수다. 음수라는 것은 이 물체가 중력에 속박되어 있다는 뜻이고, ε = 0이 되는 순간(v = √(2μ/r)) 이 물체는 탈출 속도에 도달해 무한히 멀어진다. 높은 궤도는 에너지가 덜 음수이므로(중력 퍼텐셜 에너지가 덜 낮으므로) 에너지적으로 더 높은 상태다. 스스로 질문해보자 — 낮은 궤도에서 높은 궤도로 가려면 에너지를 더해야 하는가, 빼야 하는가? 그리고 로켓은 어떻게 에너지를 더하는가? 이 질문의 답이 지금부터 배울 모든 내용의 출발점이다.

본 내용 1부 — 호만 전이: 연료 없이 우주를 건너는 법

1925년, 독일의 공학자 발터 호만(Walter Hohmann)은 두 원형 궤도 사이를 이동하는 가장 에너지 효율적인 방법을 순수 수학으로 증명했다. 그는 로켓을 직접 만들거나 실험한 적이 없었다. 오직 궤도 역학의 원리만으로, 스페이스X가 수십 년 뒤에 사용할 방법을 예언한 것이다. 그 방법이 바로 **호만 전이(Hohmann Transfer)**다.

아이디어는 단순하다. 저궤도(r₁)에서 고궤도(r₂)로 가고 싶다면, 두 궤도를 부드럽게 연결하는 **타원형 전이 궤도(Transfer Ellipse)**를 이용한다. 이 타원의 근지점(periapsis)은 r₁이고 원지점(apoapsis)은 r₂다. 즉 전이 타원의 반장축은 a_transfer = (r₁ + r₂) / 2 이 된다. 로켓은 두 번의 짧은 엔진 점화, 즉 두 번의 ΔV(델타 V, 속도 변화량) 를 가한다. 첫 번째 점화(ΔV₁)는 저궤도에서 전이 타원으로 진입하기 위해, 두 번째 점화(ΔV₂)는 고궤도에 도달했을 때 타원 궤도에서 원형 고궤도로 순환화(circularize)하기 위해서다.

각 ΔV를 계산하기 위해 **비스-비바 방정식(Vis-Viva Equation)**을 사용한다. 이 이름은 라틴어로 "살아있는 힘"을 뜻하며, 17세기 라이프니츠가 에너지를 묘사하기 위해 사용했다. 현대적으로는 v² = μ(2/r - 1/a) 로 표현된다. 이 방정식은 어떤 타원 궤도에서든, 특정 위치 r에서의 속도를 반장축 a만 알면 계산할 수 있게 해준다. 원형 궤도(a = r)에 대입하면 v = √(μ/r)이 회복되는 것을 확인해보라.

[노트 기록] 호만 전이 ΔV 계산 공식:

전이 타원 반장축: a_t = (r₁ + r₂) / 2
전이 타원 근지점 속도(저궤도에서): v_t1 = √(μ(2/r₁ - 1/a_t))
저궤도 원형 속도: v_c1 = √(μ/r₁)
ΔV₁ = v_t1 - v_c1 (가속, 양수)
전이 타원 원지점 속도(고궤도 도착 시): v_t2 = √(μ(2/r₂ - 1/a_t))
고궤도 원형 속도: v_c2 = √(μ/r₂)
ΔV₂ = v_c2 - v_t2 (가속, 양수)
총 ΔV = ΔV₁ + ΔV₂

여기서 잠깐, 앞서 배운 개념을 떠올려보자. 고궤도는 느리다고 했다. 그런데 전이 타원의 원지점 속도(v_t2)는 고궤도 원형 속도(v_c2)보다 느리다. 왜냐하면 전이 타원은 에너지가 고궤도 원형 궤도보다 낮기 때문이다(a_t < r₂). 따라서 ΔV₂도 양수, 즉 추가 가속이 필요하다. 우주에서는 높은 곳으로 갈수록 두 번 모두 가속해야 한다. 이것이 직관과 다른 점이다 — 많은 사람들이 "올라가려면 한 번 가속하고, 거기서 브레이크 걸면 되지 않나?"라고 생각하지만, 실제로는 두 번 다 가속이다. 왜 그런지 에너지 보존으로 스스로 증명해보길 권한다.

전이 시간(Transfer Time), 즉 전이 타원을 절반 돌아 고궤도에 도달하는 데 걸리는 시간도 계산할 수 있다. t_transfer = π√(a_t³/μ) — 전체 주기의 절반이기 때문이다. 지구-화성 호만 전이의 경우, 이 시간은 약 259일이다. 스페이스X가 화성까지 몇 달이 걸린다고 말하는 것이 이 계산에서 나온다.

본 내용 2부 — 다단 로켓: 왜 버리면서 올라가는가

자, 이제 ΔV가 얼마나 필요한지 안다. 그렇다면 로켓은 그 ΔV를 어떻게 만들어내는가? 그리고 왜 팰컨9은 단계적으로 분리되는가? 이 답은 1903년 러시아 수학자 콘스탄틴 치올콥스키(Konstantin Tsiolkovsky)가 제시했다. 그는 우주 여행을 꿈꿨지만 실제로 로켓을 만들지 않았다. 그가 남긴 것은 치올콥스키 로켓 방정식(Tsiolkovsky Rocket Equation) 하나였다.

ΔV = v_e × ln(m₀ / m_f)

이 방정식에서 v_e는 배기 속도(exhaust velocity), m₀는 초기 질량(연료 포함), m_f는 최종 질량(연료 소진 후)이다. ln은 자연로그다. 배기 속도는 **비추력(Specific Impulse, Isp)**으로도 표현하는데, v_e = Isp × g₀ (g₀ = 9.80665 m/s²)다. Isp는 단위 무게의 추진제가 1초 동안 만들어내는 추력으로, 단위는 초(second)다. 이것이 로켓 엔진의 효율을 나타내는 가장 표준적인 지표다. 자동차의 연비(km/L)와 비슷한 개념이다.

[노트 기록] 주요 추진제 조합의 Isp (진공 기준):

LOX/LH₂ (액체산소/액체수소): ~450 s — 고효율, 수소의 낮은 밀도가 단점
LOX/RP-1 (케로신): ~335 s — 팰컨9, 사투른V 1단 사용, 취급이 쉬움
N₂O₄/UDMH (사산화질소/비대칭 디메틸히드라진): ~340 s — 하이퍼골릭(자연 발화), 위성 기동에 많이 사용
고체 추진제: ~280 s — 간단하지만 Isp 낮고 중단 불가

이제 핵심 문제로 들어간다. 로켓 방정식의 ln(m₀/mf) 부분을 보라. 질량비(mass ratio)를 높이면 ΔV가 늘어난다. 하지만 연료를 더 넣으면 초기 질량이 늘어나고, 그 연료를 가속시키기 위한 연료가 또 필요하고, 그 연료를 가속시키기 위한 연료가 또 필요하다 — 이것을 **추진제 질량의 지수적 폭증 문제(Tyranny of the Rocket Equation)**라고 부른다. 지구 저궤도(LEO)까지 가는 데 필요한 ΔV는 약 9,400 m/s인데, Isp = 335s 엔진 하나로 이를 달성하려면 질량비가 e^(9400/3286) ≈ 17.3이 나온다. 즉 연료가 전체 질량의 (1 - 1/17.3) ≈ 94%를 차지해야 한다. 실제 1단 로켓이 90%+ 이상 연료인 이유가 이것이다.

다단(Multistage) 설계는 이 문제를 우아하게 해결한다. 연료가 다 탄 빈 탱크와 엔진(이를 **구조 질량(Structural Mass)**이라 한다)을 버림으로써, 이후 단계에서는 훨씬 가벼운 초기 질량으로 다시 높은 질량비를 달성한다. 각 단의 ΔV를 모두 더하면 전체 ΔV가 된다. 수학적으로, n단 로켓의 총 ΔV는:

ΔV_total = Σᵢ vₑᵢ × ln(m₀ᵢ / mfᵢ)

단 수가 늘어날수록 효율이 좋아지지만, 동시에 복잡도와 비용도 늘어난다. 그렇다면 최적 단 수는 몇 개인가? 이론적으로 모든 단의 mass ratio가 동일할 때 효율이 극대화되지만, 실제로는 각 단의 Isp, 구조 질량 비율(ε = m_structural / m_propellant), 그리고 경제성이 모두 다르기 때문에 수치 최적화를 사용한다. 팰컨9가 2단, 사투른V가 3단인 것은 이런 트레이드오프의 결과다.

**구조 질량 비율(ε, epsilon)**은 설계의 현실적 한계를 반영하는 중요한 변수다. 아무리 비우려 해도 탱크, 엔진, 배선 등의 무게가 남기 때문이다. ε = m_구조 / (m_구조 + m_추진제) 로 정의되며, 현대 로켓에서는 보통 0.05~0.15 수준이다. 재사용 로켓은 착륙 다리, 그리드 핀 등 추가 구조물 때문에 ε이 일반적으로 더 크다. 이것이 재사용 로켓이 성능 면에서는 소모성 로켓보다 불리한 이유이기도 하다 — 그럼에도 비용 절감이 이를 상쇄한다.

본 내용 3부 — 재사용 착륙: 중력과의 마지막 협상

팰컨9의 1단 부스터가 발사대로 되돌아오는 장면을 본 적이 있다면, 그것이 얼마나 터무니없이 정밀한 제어인지 실감할 것이다. 수십 미터 높이의 금속 기둥이 수직으로 서서, 정확히 목표 지점 위에서 멈추는 것이다. 이 과정을 제어하는 시스템은 크게 세 분야로 나뉜다: 유도(Guidance), 항법(Navigation), 제어(Control) — 통칭 GNC 시스템.

착륙의 핵심 난제는 연료-시간의 트레이드오프다. 일찍 엔진을 켜서 천천히 내려오면 연료가 많이 소모된다. 반대로 늦게까지 기다렸다 켜면 속도가 너무 빠른 상태에서 감속을 시작해야 한다. 이 극한의 경우가 자살 버닝(Suicide Burn) 혹은 **호버슬램(Hoverslam)**이다. 엔진을 최대 추력으로 켜면서 정확히 지면에 닿는 순간 속도가 0이 되도록 타이밍을 맞추는 것이다. 팰컨9 실제 착륙이 이 방식을 사용한다. 왜 이것이 최적인가? 연료를 중력에 저항하면서 "허공에 낭비"하는 시간을 최소화하기 때문이다. 이를 중력 손실(Gravity Loss) 이라 하며, 호버링 상태에서는 추진제의 100%가 중력 손실로 사라진다.

[노트 기록] 자살 버닝 착화점 추정 (1D 단순화): 지면으로부터 높이 h, 현재 속도(하강) v₀일 때, 최대 감속도가 a_max = F_thrust/m - g라면, 착화까지 남은 거리 d = v₀² / (2 × a_max)이다. 이를 **최소 착화 고도(Minimum Ignition Altitude, MIA)**라 한다. 실제로는 대기 저항, 연료 소모에 따른 질량 변화(Isp에 의한 m 감소), 6자유도 자세 변화 등이 모두 포함되므로 수치 적분으로 계산한다.

자세 제어(Attitude Control)는 더 복잡하다. 로켓은 수직으로 착지해야 하므로, 낙하 중 발생하는 기울기를 실시간으로 교정해야 한다. 팰컨9는 세 가지 방법을 조합한다. 첫째, 그리드 핀(Grid Fin) — 로켓 상단에 달린 격자 모양 핀으로, 대기 중에서 공기역학적 토크를 만들어 자세를 조정한다 (1단계에서 배운 양력/항력 개념이 여기서 다시 등장한다). 둘째, 추력 벡터 제어(Thrust Vector Control, TVC) — 엔진 노즐의 방향을 짐벌(Gimbal)로 기울여 추력 방향을 바꾼다. 노즐이 로켓 중심축에서 각도 θ만큼 기울면, 측방향 추력 F_side = F × sin(θ), 축방향 추력 F_axial = F × cos(θ)가 발생한다. 셋째, 질소 가스 냉추진기(Cold Gas Thrusters) — 극히 정밀한 자세 보정에 사용한다.

이 자세 제어의 수학적 심장은 **PID 제어기(Proportional-Integral-Derivative Controller)**다. PID는 공학에서 가장 널리 쓰이는 피드백 제어 알고리즘이다. 목표 자세각(θ_desired)과 현재 자세각(θ_current)의 차이를 **오차(error, e = θ_desired - θ_current)**라 하면, PID 제어 출력은:

u(t) = Kp × e(t) + Ki × ∫e(t)dt + Kd × de/dt

비례 항(P) 은 현재 오차에 즉각 반응한다. 오차가 크면 강하게 교정한다. 적분 항(I) 은 오차가 누적된 "편향(bias)"을 제거한다. 작은 오차가 계속 같은 방향으로 남아 있을 때 이를 서서히 제거한다. 미분 항(D) 은 오차의 변화 속도에 반응한다. 오차가 빠르게 줄어들고 있다면 제어를 약화시켜 오버슈트(overshoot), 즉 목표를 지나쳐 반대쪽으로 진동하는 현상을 막는다. Kp, Ki, Kd 세 이득(gain)을 적절히 조율하는 것을 **튜닝(Tuning)**이라 하며, 잘못 튜닝하면 로켓이 발산적으로 진동하다 분해된다. 이것이 팰컨9 초기 착륙 시험에서 종종 폭발이 발생했던 이유 중 하나다.

프로젝트 — 스스로 계산하고 설계하라

아래 세 문제는 각각 학습목표 ①②③에 대응하며, 순서대로 풀도록 설계되었다. 앞 문제의 계산 결과가 다음 문제의 입력이 되는 경우가 있으므로, 순서를 지켜라. 계산기와 노트는 자유롭게 사용하고, 로그 및 제곱근 계산은 계산기를 활용해도 된다. 단, 각 단계의 물리적 의미를 반드시 한 문장씩 서술하라.

[프로젝트 1] 화성 전이 궤도 설계 (학습목표 ①)

스페이스X의 스타십이 지구 저궤도에서 출발해 화성으로 가는 호만 전이 궤도를 계획하고 있다. 다음 데이터를 사용하라: 태양의 표준 중력 매개변수 μ_sun = 1.327 × 10²⁰ m³/s², 지구 태양 궤도 반경 r₁ = 1.496 × 10¹¹ m, 화성 태양 궤도 반경 r₂ = 2.279 × 10¹¹ m. (지구와 화성 모두 원형 궤도로 가정)

(1) 전이 타원의 반장축 a_t를 구하라.

(2) 비스-비바 방정식을 이용해, 지구 위치(r₁)에서 전이 타원에 진입하기 위한 속도 v_t1과 지구의 원형 궤도 속도 v_c1을 각각 구하고, ΔV₁을 계산하라.

(3) 화성 위치(r₂)에서 전이 타원의 속도 v_t2와 화성 원형 궤도 속도 v_c2를 구하고, ΔV₂를 계산하라.

(4) 총 ΔV_total을 구하고, 이것이 지구 궤도에서의 탈출 속도와 어떤 관계가 있는지 설명하라. (힌트: 지구 탈출 속도 v_esc = √(2μ_earth/r_LEO)와 비교하라)

(5) 전이 비행 시간 t를 구하라. 실제 화성 탐사선 발사 주기(약 26개월)가 이 전이 시간과 어떤 관련이 있는지 논리적으로 추론하라.

[프로젝트 2] 스타십 다단 추진 최적화 (학습목표 ②)

가상의 2단 로켓이 총 ΔV = 9,200 m/s를 달성해야 한다. 스타십과 유사하게 1단은 LOX/CH₄(메탄) 엔진 (Isp₁ = 330 s), 2단도 LOX/CH₄ (Isp₂ = 360 s 진공 기준)을 사용한다. 로켓의 최종 페이로드 질량 m_payload = 100 t이고, 각 단의 구조 질량 비율 ε = 0.08이다. (ε = m_structural / (m_structural + m_propellant))

(1) ΔV를 1단에 4,200 m/s, 2단에 5,000 m/s로 배분했을 때, 각 단에서 필요한 질량비 (m₀ᵢ / mfᵢ)를 구하라.

(2) 페이로드에서 거꾸로 계산하여 (2단부터 1단 순서로) 각 단의 연료 질량, 구조 질량, 초기 질량을 구하라.

(3) 전체 발사 질량 대비 페이로드 비율(Payload Fraction)을 계산하라.

(4) 이제 ΔV 배분을 1단 5,000 m/s / 2단 4,200 m/s로 바꾸어 같은 계산을 반복하라. 페이로드 비율이 어떻게 달라지는가?

(5) 두 경우를 비교하여, Isp가 낮은 1단에 더 많은 ΔV를 배분하는 것이 유리한가, 불리한가? 그 이유를 치올콥스키 방정식과 연결지어 설명하라.

(6) (심화) 만약 3단 로켓으로 변경한다면, 총 ΔV 배분의 최적화 방향에 대해 직관적으로 논하라. (계산 불요, 논리만)

[프로젝트 3] 팰컨9 착륙 시뮬레이션 (학습목표 ③)

팰컨9 1단 부스터가 귀환 중이다. 지면에서 고도 h = 3,000 m, 수직 하강 속도 v₀ = 200 m/s (하강 방향 양수), 현재 질량 m = 22,000 kg이다. 머린(Merlin) 엔진 하나의 최대 추력 F = 845 kN, g = 9.81 m/s², 이 고도에서의 공기 저항은 무시한다.

(1) 엔진 점화 없이 자유낙하한다면, 지면에 도달하는 데 걸리는 시간과 충돌 속도를 구하라.

(2) 엔진을 최대 추력으로 점화할 때 로켓이 받는 순 감속도 a_net를 구하라. (추력이 중력보다 커야 착륙 가능함에 유의)

(3) 자살 버닝 방식으로 착륙하기 위한 최소 착화 고도(MIA)를 구하라. 즉, 현재 속도 v₀에서 최대 감속도 a_net로 감속할 때 속도를 0으로 만드는 데 필요한 거리 d를 계산하고, 착화 시점의 고도(h - d 혹은 d, 어느 기준인지 명확히)를 구하라.

(4) 위 MIA보다 1,000 m 높은 곳에서 점화했다면 어떤 일이 벌어지는지 서술하라. 지면 도달 시 남은 속도(개략적으로)는 어느 방향인가?

(5) 착지 직전 로켓이 동쪽으로 2°기울어진 것이 감지되었다. TVC를 이용해 이를 교정하려 한다. 교정 토크를 만들기 위해 노즐을 기울이는 방향(동쪽 vs 서쪽)과 그 이유를 설명하라. 엔진 짐벌 각도가 1°일 때 발생하는 측방향 힘 F_side를 계산하라.

(6) (설계 문제) PID 제어기에서 Kd (미분 이득)이 너무 작으면 어떤 현상이 발생하는가? 반대로 너무 크면? 각각의 경우를 착륙 시나리오에서 구체적으로 묘사하라.