치올콥스키 로켓 방정식

물리 유체역학 마지막 수정: 2026-08-09 04:26:51

1. 개요[편집]

치올콥스키 로켓 방정식 Tsiolkovsky rocket equation
형태	$\Delta v = I_{sp} g_0 \ln(m_0/m_f)$
다른 이름	이상 로켓 방정식, the rocket equation
전제	외력 없음, 배기속도 일정, 1차원
발표	치올콥스키 1903 (유도는 1897)
선행 연구	윌리엄 무어 1813, 고더드·오베르트 독립 재발견
실무 지위	모든 임무 설계의 0번째 계산

로켓 공학의 절반은 이 한 줄이고, 나머지 절반은 이 한 줄을 어떻게든 피해 보려는 시도다.

치올콥스키 로켓 방정식은 외력이 없는 공간에서 일정한 유효배기속도 $c$ 로 추진제를 뿜는 로켓이 얻을 수 있는 속도 변화량 $\Delta v$ 가, 초기 질량과 최종 질량의 비의 자연로그에 비례한다는 관계식이다.

\Delta v = c \ln\frac{m_0}{m_f} = I_{sp}\, g_0 \ln\frac{m_0}{m_f}

$m_0$ 는 추진제를 다 채운 초기 질량, $m_f$ 는 다 태우고 남은 최종 질량, $c = I_{sp} g_0$ 는 유효배기속도다(비추력).

이 식이 잔인한 이유는 오른쪽에 추력이 없다는 데 있다. 엔진이 100 kN을 내든 1 kN을 내든, 같은 추진제를 같은 배기속도로 다 태우면 얻는 $\Delta v$ 가 같다. 추력은 그 $\Delta v$ 를 얼마나 빨리 얻느냐만 정하고, 얼마나 얻느냐는 오직 질량비와 배기속도가 정한다.¹ 로켓 설계 회의가 늘 “그램 단위로 깎기”로 끝나는 이유가 이것이다.

2. 유도 — 순간 정지 좌표계 한 장[편집]

시각 $t$ 에 질량 $m$ , 속도 $v$ 인 로켓이 $dt$ 동안 질량 $-dm > 0$ 의 추진제를 로켓 기준 속도 $c$ 로 뒤로 뿜는다. 외력이 없으면 계 전체의 운동량이 보존되므로, 로켓과 같이 움직이는 관성계에서

m\,dv + c\,dm = 0 \quad \Longrightarrow \quad dv = -c\,\frac{dm}{m}

$m_0$ 에서 $m_f$ 까지 적분하면 곧바로 위 식이 나온다. 유도에 들어간 가정은 셋뿐이다 — 외력 없음(중력·항력 무시), $c$ 일정, 1차원. 세 가정은 전부 실제 발사에서 깨지며, 깨지는 만큼이 뒤에 나올 “손실”이다.

$c$ 가 변하면 로그가 적분으로 바뀔 뿐 구조는 그대로다.

\Delta v = \int_{m_f}^{m_0} c(m)\,\frac{dm}{m}

해면에서 진공까지 올라가며 $c$ 가 변하는 1단 엔진의 실제 $\Delta v$ 는 이 적분으로 계산해야 하고, 그래서 궤적 코드는 매 스텝 고도에 맞는 $I_{sp}$ 를 다시 읽는다.

3. 지수의 폭정[편집]

식을 뒤집으면 필요한 질량비가 나온다.

\frac{m_0}{m_f} = \exp\!\left(\frac{\Delta v}{c}\right)

$\Delta v$ 가 선형으로 늘 때 질량비는 지수로 늘어난다. $c \approx 3400$ m/s(케로신/액체산소급)로 몇 개만 계산해 보면 감이 온다.

요구 $\Delta v$	질량비 $m_0/m_f$	추진제 질량 분율
3 km/s	2.4	58%
6 km/s	5.8	83%
9 km/s	14.1	93%
12 km/s	34.1	97%

여기에 구조계수 $\varepsilon$ (추진제를 담는 탱크·엔진·구조가 단 총질량에서 차지하는 비율)을 얹으면 벽이 보인다. 페이로드를 0으로 밀어붙인 단일 단이 낼 수 있는 최대 $\Delta v$ 는

\Delta v_{\max} = -c \ln \varepsilon

이다. $\varepsilon = 0.08$ , $c = 3400$ m/s면 8.6 km/s. 지구 저궤도 진입에 손실 포함 9.39.5 km/s가 필요하다는 점을 생각하면 페이로드 0으로도 못 간다는 뜻이다. 액체수소( $c \approx 4400$ m/s)로 바꾸면 넘어서긴 하지만 저밀도 극저온 탱크 때문에 $\varepsilon$ 이 나빠져, 남는 페이로드 분율이 12% 언저리로 쪼그라든다. 단일단 궤도진입(SSTO)이 “물리적으로 불가능”은 아니면서 “공학적으로 늘 지는” 이유가 이 산수에 다 들어 있다.

4. 다단화와 최적 Δv 배분[편집]

지수의 벽을 뚫는 유일한 실용 수단이 다단화다. 빈 탱크를 버리면 그 순간 $m_f$ 가 줄어들고, 다음 단은 새 로그를 처음부터 다시 쓴다. 총 $\Delta v$ 는 단별 합이 된다.

\Delta v_{\text{total}} = \sum_i c_i \ln\frac{1+\lambda_i}{\varepsilon_i + \lambda_i}

여기서 $\lambda_i$ 는 단 $i$ 의 페이로드비다. 그럼 총 $\Delta v$ 를 단들에 어떻게 나눠 줘야 최적인가. “1단이 힘드니까 1단에 많이”도 아니고 “공평하게 반반”도 아니다. 전체 발사질량을 고정하고 총 $\Delta v$ 를 최대화하는 제약 최적화 문제이므로 라그랑주 승수법을 쓴다. 승수를 $p$ 라 하면 최적 단 질량비 $n_i = m_{0i}/m_{fi}$ 가 다음 조건을 만족한다.

n_i = \frac{c_i\,p - 1}{c_i\,\varepsilon_i\,p}

$p$ 는 총 $\Delta v$ 제약을 만족하도록 1차원 방정식을 풀어 정한다(뉴턴-랩슨법 몇 번이면 끝난다). 여기서 곧바로 나오는 유용한 따름정리 — 모든 단의 $c_i$ 와 $\varepsilon_i$ 가 같다면 모든 $n_i$ 가 같아지고, 결국 $\Delta v$ 를 단들에 균등 분배하는 것이 최적이다. 실제 발사체는 단마다 엔진과 구조가 다르므로(상단이 대개 $c$ 가 높고 $\varepsilon$ 은 나쁘다) 이 균등 분배에서 벗어나며, 상단 쪽에 $\Delta v$ 를 더 주는 방향으로 기운다.

단 수를 늘리면 이론적으로는 계속 이득이지만 수확이 급격히 체감하고, 분리 기구·배선·화공품이 단마다 새로 붙는다. 실물이 2~3단에서 멈추는 것은 방정식이 아니라 신뢰성과 원가가 정한 경계다.²

5. 이상 방정식이 빼먹은 것들[편집]

발사체 엔지니어가 말하는 “요구 $\Delta v$ “는 궤도속도가 아니다. 이상 방정식이 무시한 항들이 전부 청구서로 돌아오기 때문이다. 지구 저궤도의 경우 궤도속도는 약 7.8 km/s인데 실제로 예산 잡는 값은 9.3~9.5 km/s다. 차액의 정체는 이렇다.

\Delta v_{\text{ideal}} = \Delta v_{\text{orbit}} + \underbrace{\int g\sin\gamma\,dt}_{\text{중력}} + \underbrace{\int \frac{D}{m}dt}_{\text{항력}} + \underbrace{\int \frac{T}{m}(1-\cos\alpha)\,dt}_{\text{조향}} + \underbrace{\int \frac{p_a A_e}{m}dt}_{\text{배압}}

중력 손실 — 가장 크다(전형적으로 1.2~1.5 km/s). 수직으로 서 있는 동안 추력의 일부가 그냥 중력과 비긴다. 이 항이 “추력은 $\Delta v$ 에 안 들어간다”는 앞의 문장을 슬쩍 뒤집는 지점이다. 추력 대 중량비가 낮으면 수직 구간에 오래 머물러 중력 손실이 불어난다. 이상 방정식에는 추력이 없지만 현실에는 있다.
항력 손실 — 의외로 작다(0.1~0.15 km/s). 대기가 두꺼운 구간에서는 아직 느리기 때문이다. 그래서 최대 동압(max-Q) 구간에서 스로틀을 죄는 것은 항력 $\Delta v$ 보다 구조 하중 때문이다.
조향 손실 — 추력 방향과 속도 방향이 어긋난 각 $\alpha$ 의 코사인 손실(0.05~0.1 km/s). 중력 선회를 쓰는 이유가 이것을 0에 가깝게 두기 위해서다.
배압 손실 — 해면에서 노즐 출구를 대기가 밀어 $I_{sp}$ 를 깎는 항(0.2~0.4 km/s). 이 항을 고도 전 구간에서 줄이려는 시도가 에어로스파이크 노즐이다.

반대로 되돌려 받는 항도 있다. 적도 근처에서 동쪽으로 쏘면 지구 자전 속도(적도 465 m/s)를 공짜로 얻는다. 발사장 위도가 임무 설계의 첫 줄에 오는 이유다.

6. 손실을 실제로 구하는 법 — 궤적 적분[편집]

위 적분들은 궤적을 모르면 못 구하고, 궤적은 손실을 모르면 못 구한다. 닫힌 해가 없으니 수치적으로 푼다. 표준은 중력 선회(gravity turn) 상승 궤적 ODE계다.

\dot v = \frac{T}{m} - \frac{D}{m} - g\sin\gamma, \qquad v\dot\gamma = -\left(g - \frac{v^2}{r}\right)\cos\gamma

\dot h = v\sin\gamma, \qquad \dot m = -\frac{T}{c}

$\gamma$ 가 비행경로각이다. 중력 선회의 요체는 초기에 아주 작은 피치오버 한 번만 주고 그 뒤로는 중력이 알아서 기수를 눕히게 두는 것 — 받음각이 0에 가깝게 유지되므로 조향 손실과 횡방향 공력 하중이 동시에 최소가 된다.

여기서 실무적 함정이 두 개 있다.

초기 조건 민감도. 피치오버 시각과 각도를 1° 바꾸면 궤도 진입 조건이 통째로 달라진다. 사실상 슈팅법으로 푸는 경계값 문제이고, 그래서 상승 궤적 설계는 최적 제어 문제로 정식화하는 것이 정석이다. 평지·일정추력 근사에서 변분법이 주는 유명한 답이 선형 탄젠트 조향 법칙 — 최적 추력 방향각의 탄젠트가 시간의 1차 함수라는 것이다.
강성과 불연속. 단 분리, 페어링 분리, 스로틀 변경마다 우변이 불연속으로 튄다. 고정 스텝 룽게-쿠타법으로 이벤트를 통과하면 오차가 그 자리에서 생기므로, 이벤트 검출을 붙여 스텝을 정확히 끊고 다시 시작해야 한다. 대기 모델·추력 테이블 보간의 불연속도 같은 취급을 받아야 한다.

궤적을 다 적분하고 나면 실제 소모 추진제에서 역산한 “달성 $\Delta v$ “와, 이상 방정식이 그 질량비에서 약속했던 값의 차이가 나온다. 그 차이가 손실의 정의다. 손실은 물리 상수가 아니라 궤적의 함수라는 뜻이며, 그래서 “중력 손실 1.5 km/s”라는 숫자는 언제나 특정 발사체·특정 궤적에 붙는 꼬리표다.

7. 역사와 오해[편집]

치올콥스키(Константин Циолковский)는 1897년에 이 관계를 유도해 1903년 논문 「반작용 장치를 이용한 우주 공간의 탐구」에서 발표했다. 다만 최초는 아니다. 1813년 영국의 윌리엄 무어(William Moore)가 로켓 운동에 관한 논고에서 사실상 같은 결과를 얻었다는 것이 정설이고, 이후 로버트 고더드와 헤르만 오베르트도 독립적으로 재발견했다. 치올콥스키의 진짜 기여는 방정식 자체보다 그 방정식이 무엇을 요구하는지를 끝까지 밀어붙인 것 이다 — 액체 추진제, 다단 로켓(“우주 열차”), 궤도 스테이션까지 이 한 줄의 논리적 귀결로 제시했다.

흔한 오해 몇 가지도 짚어 두자.

“ $g_0$ 가 들어 있으니 중력장 안에서만 성립한다” — 아니다. $g_0 = 9.80665$ m/s²는 비추력을 초 단위로 쓰기 위한 정의 상수일 뿐이고, $c$ 로 쓰면 사라진다.
“이온엔진은 $I_{sp}$ 가 10배니 10배 좋다” — $\Delta v$ 는 그렇지만 그 $\Delta v$ 를 얻는 데 몇 달이 걸린다. 발사에는 추력 대 중량비가 필요하다.
“오베르트 효과는 로켓 방정식을 위반한다” — 아니다. 근점에서 연소하면 같은 $\Delta v$ 가 더 큰 궤도 에너지 변화를 만드는 것이지, $\Delta v$ 자체가 늘어나는 것이 아니다.³

8. 관련 문서[편집]

그래서 “우리 엔진 추력이 경쟁사보다 20% 높다”는 홍보 문구는 $\Delta v$ 에 대해 아무 말도 하지 않는다. 물론 중력 손실을 통해 간접적으로는 영향을 준다 — 다만 그건 이상 방정식 바깥의 이야기이고, 홍보 담당자가 그 구분을 하고 있을 확률은 높지 않다. ↩
단 수 $N$ 을 무한대로 보내면 연속적으로 빈 구조를 버리는 “연속 다단 로켓”의 극한이 나오고, 이론 성능은 계속 좋아진다. 실제로 이걸 근사한 물건이 다발 부스터를 순차 분리하는 구성이다. 물론 분리 기구가 하나 늘 때마다 실패 모드도 하나 는다는 사실은 방정식에 안 적혀 있다. ↩
오베르트 효과의 정체는 운동에너지가 속도의 제곱이라는 것뿐이다. $\Delta(v^2/2) = v\Delta v + \Delta v^2/2$ 이므로 이미 빠를 때( $v$ 가 클 때) 같은 $\Delta v$ 를 쓰면 에너지 이득이 크다. 로켓 방정식이 아니라 고등학교 물리가 하는 일이다. ↩