치올콥스키 로켓 방정식: 운동량 보존으로 비추력과 질량비를 계산하는 방법

로켓은 뒤로 가스를 내보내며 앞으로 속도를 얻습니다. 자동차처럼 지면을 밀거나 비행기처럼 대기의 공기를 빨아들일 필요가 없습니다. 연료와 산화제를 싣고 온 로켓은 자기 질량의 일부를 고속으로 분사하고, 로켓과 배기가스 전체의 운동량이 보존되도록 반대 방향으로 가속합니다. 그래서 진공에서도 작동합니다.

이 과정을 계산한 식이 치올콥스키 로켓 방정식입니다. 처음 질량을 m₀, 연소가 끝난 뒤 질량을 m_f, 유효 배기속도를 vₑ라고 하면 이상적인 속도 변화량은 다음과 같습니다.

Δv = vₑ ln(m₀/m_f) = Isp·g₀·ln(m₀/m_f)

짧은 식이지만 로켓 설계의 가장 까다로운 현실이 담겨 있습니다. 속도 변화량은 배기속도에는 비례하지만 질량비에는 자연로그로만 증가합니다. 더 큰 Δv를 얻으려고 추진제를 실으면 그 추진제 자체도 함께 가속해야 하므로 필요한 처음 질량이 지수적으로 커집니다.

치올콥스키 로켓 방정식의 기호부터 정리하기

m₀에는 추진제, 탱크, 엔진, 구조물과 탑재물이 모두 들어갑니다. m_f에는 해당 연소에서 소비한 추진제만 빠지고 구조물과 탑재물, 남은 추진제가 포함됩니다. 흔히 m_f를 건조질량이라고 간단히 부르지만 연소가 일부 추진제만 쓰고 끝난다면 엄밀한 m_f는 완전한 건조질량과 다를 수 있습니다.

왜 F=ma만 그대로 쓰면 부족할까

보통 역학 문제에서는 물체의 질량이 일정하다고 가정합니다. 로켓은 추진제를 계속 밖으로 내보내므로 질량이 시간에 따라 줄어듭니다. 로켓만 떼어 놓고 `F=ma`를 기계적으로 적용하면 빠져나간 질량이 가져가는 운동량을 놓치기 쉽습니다.

해결 방법은 아주 짧은 시간 동안 로켓과 방출된 배기가스를 하나의 계로 보고 방출 전후의 전체 운동량을 비교하는 것입니다. 외력이 없다고 가정하면 이 계의 전체 운동량은 같습니다. 그 결과 질량이 변하는 로켓의 미소 속도 변화식을 얻을 수 있습니다.

운동량 보존으로 미소 변화식을 유도하기

어느 순간 로켓의 질량을 m, 관성계에서 본 속도를 v라고 하겠습니다. 짧은 시간 뒤 로켓은 양의 질량 dm만큼을 뒤로 분사합니다. 남은 로켓의 질량은 `m−dm`, 속도는 `v+dv`가 됩니다. 배기가스가 로켓에 대해 뒤쪽으로 vₑ의 속도로 나간다면 관성계에서 배기가스의 속도는 1차 근사로 `v−vₑ`입니다.

방출 전 운동량 = mv

방출 후 운동량 = (m−dm)(v+dv) + dm(v−vₑ)

외력이 없다면 두 운동량은 같습니다. 오른쪽을 전개하면 `mv + m·dv − dm·dv − vₑ·dm`이 됩니다. dm과 dv를 곱한 항은 두 미소량의 곱이므로 1차 근사에서 제외합니다. 양쪽의 mv도 지우면 다음 관계만 남습니다.

m·dv = vₑ·dm

여기서는 dm을 밖으로 나간 양의 질량으로 정의했습니다. 로켓 자체의 질량 변화 dM은 음수이므로 `dm=−dM`입니다. 로켓 질량을 M으로 표기하면 같은 식을 다음처럼 쓸 수 있습니다.

M·dv = −vₑ·dM

dM이 음수이므로 오른쪽은 양수가 되고 로켓의 속도는 증가합니다. 부호가 헷갈리는 이유는 어떤 교재는 방출된 질량을 양수로 두고, 다른 교재는 로켓의 질량 감소량을 음수 미분으로 쓰기 때문입니다. 정의를 먼저 확인하면 결과는 같습니다.

적분하면 자연로그가 나타난다

유효 배기속도 vₑ가 연소 중 일정하다고 가정하고 미소 변화식을 정리하면 `dv=−vₑdM/M`입니다. 연소 시작 질량 m₀에서 종료 질량 m_f까지 적분합니다.

∫(v₀→v_f) dv = −vₑ∫(m₀→m_f) dM/M

v_f−v₀ = −vₑ[ln M]_(m₀→m_f) = vₑ ln(m₀/m_f)

따라서 `Δv=v_f−v₀`라고 두면 치올콥스키 로켓 방정식이 완성됩니다. 자연로그가 들어가는 까닭은 연소하는 동안 분모의 로켓 질량 M이 계속 변하고 `1/M`을 적분하기 때문입니다. 같은 양의 추진제를 분사해도 로켓이 가벼워진 뒤에 분사한 추진제가 더 큰 속도 변화를 만듭니다.

델타브이는 실제 속도와 무엇이 다를까

Δv는 단순히 속도계에 표시되는 현재 속도가 아닙니다. 우주선이 엔진 작동으로 만들어 낼 수 있는 속도 변화의 능력 또는 임무에서 필요한 속도 변화의 합을 뜻합니다. 방향을 바꾸는 감속도 Δv를 소비합니다. 예를 들어 같은 크기로 가속했다가 나중에 감속하면 최종 속도가 처음과 같아질 수 있지만 추진제는 두 번의 기동에서 모두 사용됩니다.

이상식의 Δv는 외력이 없고 배기속도가 일정할 때 로켓이 얻는 최대 속도 변화입니다. 지상 발사에서는 중력에 맞서 엔진을 켜 두는 동안 생기는 중력 손실, 대기 저항, 조향 손실이 더해집니다. 따라서 목표 궤도의 속도와 발사체가 준비해야 하는 총 Δv 예산은 같지 않습니다.

유효 배기속도는 단순한 가스 속도가 아니다

로켓 노즐의 추력은 배기가스가 가져가는 운동량뿐 아니라 노즐 출구 압력과 주변 압력의 차이에서도 나옵니다. 질량유량을 ṁ, 실제 노즐 출구 속도를 Vₑ, 출구 면적을 Aₑ, 출구 압력을 pₑ, 주변 압력을 p₀라고 하면 추력은 다음처럼 나타낼 수 있습니다.

F = ṁVₑ + Aₑ(pₑ−p₀)

이 식을 `F=ṁvₑ` 형태로 묶어 주는 값이 유효 배기속도입니다. 즉 `vₑ=F/ṁ`이며 압력 추력까지 속도와 같은 단위로 환산한 성능값입니다. 진공과 해수면에서 같은 엔진의 성능 표기가 달라질 수 있는 이유도 주변 압력 p₀가 달라지기 때문입니다.

비추력 Isp가 초 단위인 이유

비추력은 추력을 추진제의 중량유량으로 나눈 값입니다. 표준중력가속도 `g₀=9.80665 m/s²`를 사용하면 다음 관계가 성립합니다.

Isp = F/(ṁg₀) = vₑ/g₀

추력 F의 단위는 N, 질량유량 ṁ의 단위는 kg/s입니다. `F/(ṁg₀)`의 단위를 정리하면 초가 남습니다. 그래서 비추력은 보통 s로 표시합니다. 하지만 450초 비추력의 엔진이 반드시 450초 동안 연소한다는 뜻은 아닙니다. 같은 추진제 중량유량으로 얼마나 큰 추력을 만드는지를 비교하기 위한 성능 지표입니다.

비추력 300초라면 유효 배기속도는 `300×9.80665≈2,942 m/s`, 450초라면 약 `4,413 m/s`입니다. 같은 질량비에서는 Isp가 높은 엔진이 더 큰 Δv를 냅니다. 다만 엔진의 추력, 추진제 밀도, 탱크 질량, 전력 요구량과 임무시간은 Isp 하나만으로 결정되지 않습니다.

질량비가 2, 4, 8일 때 델타브이 계산

유효 배기속도를 3,000 m/s로 놓고 질량비만 바꿔 보겠습니다. 질량비가 2라는 것은 연소 전 질량이 연소 후 질량의 두 배라는 뜻입니다.

질량비를 2에서 4로, 다시 8과 16으로 두 배씩 늘릴 때마다 Δv는 같은 양인 `vₑln2`만큼 늘어납니다. 반면 초기 질량에서 추진제가 차지하는 비율은 50%, 75%, 87.5%, 93.75%로 빠르게 100%에 가까워집니다. 구조물과 탑재물을 위한 몫이 점점 줄어드는 로켓 방정식의 횡포가 여기서 나타납니다.

필요한 질량비를 거꾸로 구하기

목표 Δv와 엔진의 비추력을 알고 있다면 방정식을 질량비에 대해 풀 수 있습니다. 양변을 vₑ로 나누고 지수함수를 취하면 다음 식이 됩니다.

m₀/m_f = exp(Δv/vₑ) = exp[Δv/(Isp·g₀)]

예를 들어 Isp가 300초이고 이상적인 Δv가 3,000 m/s 필요하다면 유효 배기속도는 약 2,942 m/s입니다. 질량비는 `exp(3000/2942)≈2.77`입니다. 따라서 최종 질량은 초기 질량의 약 36.1%, 소비한 추진제는 약 63.9%입니다. 이 계산에는 탱크와 엔진을 버리는 단 분리나 외력 손실이 들어 있지 않습니다.

추진제 질량분율 공식

연소에서 소비하는 추진제 질량을 m_p라고 하면 `m_p=m₀−m_f`입니다. 질량비 식을 이용하면 필요한 추진제의 초기 질량 대비 비율을 바로 구할 수 있습니다.

m_p/m₀ = 1−m_f/m₀ = 1−exp[−Δv/(Isp·g₀)]

이 비율은 순수하게 방정식에서 요구하는 소비 추진제 몫입니다. 나머지가 전부 탑재물이 되는 것은 아닙니다. 남은 질량에는 탱크, 엔진, 배관, 단열재, 항법 장치와 구조물이 먼저 들어가야 합니다. 실제 탑재물 비율이 추진제 이외의 비율보다 훨씬 작을 수 있는 이유입니다.

속도 요구량이 조금 늘어도 질량이 크게 늘어나는 이유

로켓 방정식을 질량비로 바꾸면 자연로그의 반대편에 지수함수가 나타납니다. Δv가 같은 양만큼 늘어날 때 필요한 질량비는 일정한 배수로 커집니다. 예를 들어 vₑ만큼 Δv가 증가하면 질량비에는 e≈2.718이 곱해집니다.

종이를 접을 때 두께가 단계마다 두 배가 되는 지수 증가와 수학적 구조가 이어집니다. 종이 문제에서는 횟수가 늘 때 두께가 지수적으로 커지고, 로켓에서는 요구 Δv가 늘 때 필요한 질량비가 지수적으로 커집니다. 두 경우 모두 초반의 작은 차이가 반복되면서 매우 큰 값이 됩니다.

자연로그 때문에 생기는 수익 체감

질량비를 계속 높여도 Δv는 로그로 증가합니다. 질량비 2에서 4로 갈 때 얻는 추가 Δv와 8에서 16으로 갈 때 얻는 추가 Δv는 모두 `vₑln2`입니다. 그러나 뒤쪽 구간에서는 이미 매우 큰 추진제 비율에 다시 추진제를 더해야 합니다.

복리의 두 배 시간을 로그로 근사하는 72의 법칙도 지수와 로그가 서로 역함수라는 점을 이용합니다. 복리에서는 성장률로 시간을 찾고, 로켓 방정식에서는 요구 속도와 배기속도로 질량비를 찾습니다. 적용 분야는 다르지만 로그가 지수적 변화를 거꾸로 푸는 도구라는 점은 같습니다.

다단 로켓이 질량 문제를 줄이는 방법

한 단의 탱크와 엔진은 추진제를 다 쓴 뒤에도 질량으로 남아 다음 가속을 방해합니다. 다단 로켓은 역할을 마친 구조물을 분리해 이후 연소에서 가속해야 할 m_f를 줄입니다. 각 단의 이상적 Δv를 계산해 더하면 전체 Δv를 얻습니다.

Δv_total = Σ[vₑ,i ln(m₀,i/m_f,i)]

각 단의 m₀,i와 m_f,i에는 그 순간 위에 실린 모든 상단과 탑재물이 포함됩니다. 분리한 구조물은 다음 단의 초기 질량에서는 제외됩니다. 단을 무한히 늘리면 항상 좋은 것은 아닙니다. 단마다 엔진, 분리장치와 제어계가 추가되고 신뢰성과 운용 복잡도도 달라집니다.

중력과 궤도 속도는 로켓 방정식 밖에서 계산한다

치올콥스키 식은 엔진과 질량비가 제공하는 이상적 Δv를 알려 줍니다. 어느 고도의 원궤도에 필요한 속도나 행성 탈출 속도는 중력장과 궤도역학으로 따로 구합니다. 뉴턴의 만유인력과 케플러 법칙에서 거리 역제곱 중력을 유도하는 과정은 목표 속도와 궤도 조건을 정하는 쪽에 해당합니다.

임무 설계에서는 궤도 기동마다 필요한 Δv를 계산한 뒤 합계와 여유분을 정합니다. 그 Δv 예산을 로켓 방정식에 넣어 추진제와 질량비를 구합니다. 중력 방정식과 로켓 방정식은 서로 경쟁하는 설명이 아니라 임무의 요구량과 추진기관의 공급 능력을 각각 계산합니다.

이상적 로켓 방정식이 가정하는 조건

• 분석 구간에서 중력과 공기저항 같은 외력의 영향을 무시합니다.
• 유효 배기속도 vₑ가 일정하다고 가정합니다.
• 추진제 분사는 로켓 진행 방향과 정확히 반대라고 가정합니다.
• 로켓은 점질량처럼 다루며 회전과 자세제어 손실을 별도로 계산합니다.
• 질량은 연속적으로 방출되고 상대론적 효과가 필요하지 않은 속도라고 봅니다.
• 한 번의 연소 구간에서는 단 분리처럼 불연속적인 질량 변화를 따로 처리합니다.

우주 공간에서 짧게 하는 기동은 외력 손실이 작아 이상식에 가깝게 계산할 수 있습니다. 반면 지상 발사는 대기와 중력을 통과하고 추력 방향도 계속 바뀌므로 실제 궤적 해석이 필요합니다. 이상식은 쓸모없는 단순화가 아니라 엔진 성능과 질량 구조가 만드는 근본적인 한계를 분리해 보여 주는 기준식입니다.

비추력이 높으면 항상 좋은 엔진일까

같은 질량비라면 비추력이 높을수록 Δv가 커집니다. 그러나 로켓 엔진을 고를 때는 추력도 중요합니다. 이온 추진기는 매우 높은 비추력을 낼 수 있지만 추력이 작아 오랜 시간 작동하는 우주 임무에 적합합니다. 지상에서 무거운 발사체를 들어 올리려면 중량보다 큰 추력이 필요하므로 비추력만 높다고 발사용 1단 엔진이 되는 것은 아닙니다.

추진제 밀도도 탱크 크기와 구조 질량에 영향을 줍니다. 보관 온도, 독성, 점화 방식, 재시동 횟수와 비용도 함께 판단합니다. 비추력은 매우 중요한 비교 지표지만 엔진과 발사체 전체의 우열을 혼자 결정하는 점수는 아닙니다.

추력과 델타브이는 같은 값이 아니다

추력 F는 순간적인 힘이고 단위는 N입니다. Δv는 연소 전체에서 축적할 수 있는 속도 변화량이며 단위는 m/s입니다. 같은 비추력과 질량비라면 추력이 다른 두 엔진도 이상적인 총 Δv는 같을 수 있지만 연소시간은 달라집니다.

추력이 큰 엔진은 추진제를 더 빠르게 분사해 짧은 시간에 Δv를 만듭니다. 추력이 작은 엔진은 같은 추진제와 유효 배기속도라도 더 오래 작동해야 합니다. 중력장에서는 연소시간이 길수록 중력 손실이 커질 수 있으므로 실제 발사에서는 총 Δv뿐 아니라 추력대중량비와 연소시간이 중요합니다.

1903년 치올콥스키가 제시한 우주비행 계산

러시아의 교사이자 우주비행 이론가 콘스탄틴 치올콥스키는 1903년 《로켓 장치에 의한 우주 공간의 탐구》로 번역되는 논문에서 액체 추진 로켓과 우주비행을 수학적으로 다뤘습니다. NASA의 역사 자료는 이 연구에서 배기속도와 질량비가 로켓 속도를 결정하는 관계가 제시됐다고 설명합니다.

치올콥스키가 처음으로 로켓을 만든 사람은 아닙니다. 로켓은 훨씬 이전부터 존재했습니다. 그의 중요한 공헌은 가변질량 운동을 정량화해 우주비행에 필요한 성능을 계산할 수 있는 형태로 제시하고, 액체 추진제와 다단 로켓 같은 개념의 가능성을 체계적으로 논한 데 있습니다.

로켓 방정식에서 자주 생기는 계산 오류

• 상용로그 log₁₀을 넣습니다. 기본식의 ln은 자연로그입니다.
• m₀와 m_f의 순서를 바꿉니다. 정상적인 연소에서는 m₀>m_f이므로 로그값이 양수여야 합니다.
• 비추력 Isp를 그대로 배기속도로 사용합니다. 초 단위 Isp에는 g₀를 곱해야 m/s가 됩니다.
• Isp의 초를 엔진 연소시간으로 해석합니다.
• 최종 질량을 탑재물 질량만으로 둡니다. 구조물과 엔진도 남습니다.
• 추력과 Δv를 같은 성능값으로 봅니다. 하나는 힘이고 다른 하나는 속도 변화량입니다.
• 이상적 Δv를 실제 지상 발사 후 속도와 같다고 봅니다. 중력·항력·조향 손실을 별도로 고려해야 합니다.
• 여러 단의 전체 질량비 하나만 넣어 단 분리 효과를 계산합니다. 각 단의 연소와 분리를 순서대로 계산해야 합니다.

치올콥스키 로켓 방정식 FAQ

로켓은 우주에서 밀어낼 공기가 없는데 어떻게 움직이나요?

로켓은 공기를 미는 것이 아니라 탑재한 추진제를 뒤로 분사합니다. 배기가스가 뒤쪽 운동량을 얻으면 로켓은 전체 운동량 보존에 따라 앞쪽 운동량을 얻습니다. 따라서 진공에서도 추진할 수 있습니다.

질량비가 10이면 추진제가 정확히 90%인가요?

m₀/m_f=10이고 해당 연소에서 m_f가 추진제를 제외한 최종 질량이라면 소비한 추진제는 초기 질량의 90%입니다. 그러나 실제 단계에 미사용 추진제가 남거나 계산 구간이 일부 연소만 포함하면 단순한 건조질량 비율과 다를 수 있습니다.

비추력 450초는 무엇을 뜻하나요?

유효 배기속도로 환산하면 `450×9.80665≈4,413 m/s`라는 뜻입니다. 450초 동안만 작동한다는 뜻은 아닙니다. 정의상 추력을 추진제의 중량유량으로 나누기 때문에 초 단위가 남습니다.

질량비를 두 배로 하면 델타브이도 두 배가 되나요?

아닙니다. Δv는 질량비의 자연로그에 비례합니다. 질량비를 R에서 2R로 늘리면 추가되는 Δv는 언제나 `vₑln2`입니다. 질량비 자체가 두 배가 되어도 전체 Δv는 일반적으로 두 배가 되지 않습니다.

로켓 방정식은 대기 중에서도 사용할 수 있나요?

엔진과 질량비가 주는 이상적 Δv를 계산하는 기준식으로 사용할 수 있습니다. 다만 실제 궤적에는 고도에 따른 유효 배기속도 변화, 중력 손실, 공기저항과 조향 손실을 더해야 하므로 로켓 방정식 하나만으로 실제 상승 속도와 고도를 구할 수는 없습니다.

운동량 보존에서 질량비까지 핵심 정리

로켓과 배기가스 전체에 운동량 보존을 적용하면 `Mdv=−vₑdM`을 얻습니다. 유효 배기속도가 일정할 때 초기 질량 m₀에서 최종 질량 m_f까지 적분하면 `Δv=vₑln(m₀/m_f)`가 됩니다. 비추력은 `Isp=vₑ/g₀`이므로 같은 식을 `Δv=Isp·g₀·ln(m₀/m_f)`로 쓸 수 있습니다.

질량비를 역산하면 `m₀/m_f=exp[Δv/(Isp·g₀)]`입니다. 필요한 Δv가 늘수록 초기 질량은 지수적으로 커지고, 질량비를 높여 얻는 Δv는 로그로 증가합니다. 이 관계 때문에 높은 비추력, 가벼운 구조, 다단화와 정확한 Δv 예산이 모두 중요합니다. 치올콥스키 방정식은 실제 비행의 모든 손실을 계산하지는 않지만 로켓이 피할 수 없는 질량의 대가를 가장 간결하게 보여 줍니다.

참고 자료

• NASA Glenn Research Center, Ideal Rocket Equation
• NASA Glenn Research Center, Specific Impulse
• NASA Glenn Research Center, Mass Ratios
• NASA Science, Basics of Space Flight: Gravity and Mechanics
• NASA, Rockets Educator Guide
• NASA History, Challenge to Apollo: Origins of the Soviet Space Program