火箭推進原理(進階):噴嘴、壅塞流與排氣速度的天花板
從拉瓦爾噴嘴的超音速氣流,到電推進如何鬆綁火箭方程的指數暴政
排氣速度的天花板,藏在一個喇叭口裡
讀過入門篇後,你已經知道推力 $F=\dot{m}\,v_e$、比衝 $I_{sp}=v_e/g_0$,也知道火箭方程裡那個令人窒息的指數,會因為排氣速度 $v_e$ 不夠大而把燃料需求逼到天上。於是一個更尖銳的問題浮現了:$v_e$ 究竟是怎麼被「決定」的? 我們能不能單純把燃燒室燒得更熱、壓力更高,就無限制地提升排氣速度?
答案是不行——而且阻擋我們的,不是工程上的偷工減料,而是一段熱力學與氣體動力學寫死的物理。火箭引擎真正把「高溫高壓的呆滯氣體」轉換成「高速有序氣流」的關鍵元件,是那個你在每張火箭照片裡都看過、卻很少有人細想的喇叭狀噴嘴——收斂—發散噴嘴(converging–diverging nozzle),又稱拉瓦爾噴嘴(de Laval nozzle)。這篇文章要帶你鑽進燃燒室與噴嘴內部,理解 $v_e$ 的上限從哪裡來,再走到化學火箭之外,看電推進如何用完全不同的物理重新分配這場賽局。

燃燒室裡其實「沒有速度」
先破除一個關於火箭的直覺陷阱:很多人以為燃燒室是火箭最「猛」的地方,氣體在裡面已經高速衝刺。事實恰恰相反。在燃燒室深處,氣體是接近滯止(stagnation)狀態的——壓力極高(數十到數百個大氣壓)、溫度極高(攝氏三千度以上),但流速很慢。那裡儲存的是焓(enthalpy),是混亂的熱能,不是有方向的動能。
噴嘴的任務,就是一台能量轉換器:把這團混亂的熱焓,盡可能完整地轉換成單一方向的平移動能。我們可以用穩態、絕熱、無摩擦(等熵)的能量守恆來描述這條氣流。沿著流線,單位質量的總焓守恆:
$$h_0 = h + \frac{1}{2}v^2$$
其中 $h_0$ 是滯止焓(燃燒室內,$v\approx 0$ 時的焓),$h$ 與 $v$ 是噴嘴某截面處的焓與流速。把這條式子反解出流速:
$$v = \sqrt{2\left(h_0 - h\right)}$$
對理想氣體 $h=c_p T$,於是
$$v_e = \sqrt{2\,c_p\left(T_0 - T_e\right)}$$
這條式子已經透露出兩個關鍵訊息。第一,要讓 $v_e$ 大,燃燒室溫度 $T_0$ 要高、出口溫度 $T_e$ 要低——也就是說,氣體在噴嘴裡膨脹得越徹底、被「榨得越冷」,吐出的速度就越快。第二,比熱 $c_p$ 越大越好,而 $c_p$ 與氣體分子量成反比。這就從熱力學第一性原理,解釋了入門篇提過卻沒展開的事實:為什麼氫氧引擎的比衝特別高——因為它的燃燒產物是水蒸氣,平均分子量低,$c_p$ 大。
為什麼非得是「喇叭口」:壅塞流的祕密
如果只是要膨脹氣體,為什麼噴嘴要先收窄、再擴張,做成沙漏的形狀?這是整個火箭推進裡最違反直覺、也最優雅的一段流體力學。
關鍵在於可壓縮氣流的行為,會在流速超過音速時整個翻轉。把一維等熵流的連續方程與動量方程合併,可以得到著名的面積—馬赫數關係:
$$\frac{dA}{A} = \left(M^2 - 1\right)\frac{dv}{v}$$
這條式子裡的 $M$ 是馬赫數(Mach number),即流速與當地音速之比。仔細讀它:
- 當 $M<1$(次音速):$M^2-1<0$,要讓氣體加速($dv>0$),截面積必須縮小($dA<0$)。這跟我們捏住水管讓水噴得更快的直覺一致。
- 當 $M>1$(超音速):$M^2-1>0$,要讓氣體繼續加速,截面積反而必須放大。這完全違反日常直覺。
- 當 $M=1$:$dA=0$,截面積取極值,也就是噴嘴最窄的喉部(throat)。
於是答案清楚了:氣體要從次音速一路加速到超音速,幾何上必須先收斂(在次音速段加速直到喉部達到音速),再發散(在超音速段繼續加速)。喉部正是音速與超音速的分界線,一旦那裡達到 $M=1$,就發生了壅塞(choking)——通過喉部的質量流率被「鎖死」,不管你把下游噴嘴出口的背壓再降低多少,流進來的氣體量都不會再增加。
這個壅塞條件也固定了喉部與燃燒室之間的壓力比。對比熱比為 $\gamma$ 的氣體,喉部壓力 $p^*$ 與滯止壓力 $p_0$ 的關係是:
$$\frac{p^*}{p_0} = \left(\frac{2}{\gamma+1}\right)^{\frac{\gamma}{\gamma-1}} \approx 0.53\quad(\gamma\approx 1.2)$$
也就是說,燃燒室壓力一旦超過外界約兩倍,喉部就壅塞了,引擎進入「設計工作狀態」。這正是為什麼火箭引擎能穩定地以可預測的流率運作。
把引擎拆成兩半:c* 與推力係數
工程師處理火箭引擎效能時,有個極漂亮的做法:把整具引擎的性能,乾淨地拆成「燃燒室做得好不好」與「噴嘴做得好不好」兩個獨立指標。這個拆解值得每個學推進的人記在心裡。
第一個指標是特徵速度(characteristic velocity) $c^*$,它只跟燃燒室裡的熱化學有關——推進劑組合、燃燒溫度、產物分子量——而與噴嘴形狀無關:
$$c^* = \frac{p_0\, A_t}{\dot{m}}$$
其中 $A_t$ 是喉部面積。$c^*$ 衡量「在給定的燃燒室壓力下,喉部能放行多少氣體」,本質上是燃燒效率與推進劑能量品質的度量。理論上
$$c^* = \frac{\sqrt{\gamma R T_0 / \mathcal{M}}}{\gamma}\left(\frac{\gamma+1}{2}\right)^{\frac{\gamma+1}{2(\gamma-1)}}$$
$R$ 是通用氣體常數,$\mathcal{M}$ 是產物平均分子量。你又看到了 $T_0/\mathcal{M}$ 這個熟悉的組合:高溫、低分子量,永遠是推進劑配方的聖杯。
第二個指標是推力係數(thrust coefficient) $C_F$,它只反映噴嘴把壓力轉成推力的本事,與燃燒室絕對壓力高低無關:
$$C_F = \frac{F}{p_0\, A_t}$$
於是整具引擎的推力可以寫成一個無比簡潔的乘積:
$$F = C_F\, c^*\, \dot{m} \cdot \frac{1}{c^*}\cdot c^* = C_F\, p_0\, A_t$$
而排氣速度(有效值)則是兩者之積 $v_e = c^* C_F$。這個拆解的威力在於:地面測試時,工程師量測 $c^*$ 來判斷燃燒是否完全(若實測 $c^*$ 低於理論值,問題出在噴注器或混合比),量測 $C_F$ 來判斷噴嘴是否在設計點工作(若偏低,可能是膨脹不當)。兩個病灶被乾淨地分離診斷,這是火箭工程能快速迭代的根基。
膨脹的拉鋸:過度膨脹與膨脹不足
入門篇提過真空型引擎配大噴嘴,這裡我們把背後的物理講透。噴嘴出口壓力 $p_e$ 取決於面積膨脹比 $\epsilon = A_e/A_t$——出口越大,氣體膨脹越充分,$p_e$ 越低、$v_e$ 越高。理想狀態是 $p_e = p_0$(出口壓力恰好等於環境壓力),此時整道氣流的動量被最有效地利用,稱為最佳膨脹(optimal expansion)。
問題是環境壓力 $p_0$ 會變。火箭從海平面(約 101 kPa)一路爬升到真空(約 0),同一個固定形狀的噴嘴不可能對所有高度都剛好最佳:
- 過度膨脹(over-expanded):在海平面,若噴嘴太大,$p_e<p_0$,氣流出口壓力低於外界,外界大氣會把氣流往內擠壓。嚴重時噴流會在噴嘴內部就發生激震分離(shock-induced separation),氣流脫離壁面、產生側向力,可能損壞噴嘴。
- 膨脹不足(under-expanded):在高空或真空,若噴嘴太小,$p_e>p_0$,氣體離開噴嘴後還想繼續膨脹,那部分壓力能就白白浪費了。
這正是為什麼第一節(海平面工作)的引擎噴嘴相對矮胖、上面節(真空工作)的引擎噴嘴又高又大。也是工程師夢寐以求氣動尖塞噴嘴(aerospike)的原因——它讓外界大氣本身充當噴流的「虛擬外壁」,能隨高度自動調整有效膨脹比,在整個飛行包絡內都接近最佳。
動手算一下:噴嘴把氣流加速到幾倍音速?
我們用面積膨脹比反推出口馬赫數,感受一下噴嘴的威力。等熵流的面積—馬赫數關係寫成顯式:
$$\frac{A_e}{A_t} = \frac{1}{M_e}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M_e^2\right)\right]^{\frac{\gamma+1}{2(\gamma-1)}}$$
取一個典型的真空型上面節引擎,膨脹比 $\epsilon = A_e/A_t = 80$,產物比熱比 $\gamma = 1.20$。代入數值反解(這是超越方程,需數值求根),可得出口馬赫數約
$$M_e \approx 4.7$$
也就是說,氣體離開噴嘴時的速度,是它在噴嘴內當地音速的約 4.7 倍。再看溫度的下降幅度。等熵關係給出:
$$\frac{T_e}{T_0} = \left(1+\frac{\gamma-1}{2}M_e^2\right)^{-1} = \left(1 + 0.1\times 4.7^2\right)^{-1} \approx 0.31$$
若燃燒室溫度 $T_0 = 3500\ \text{K}$,則出口溫度只剩約 $1085\ \text{K}$。將近 七成的熱能被轉換成了平移動能——這就是噴嘴的本事:把無方向的高溫亂流,梳理成一道又冷又快、整齊向後的氣箭。把這個 $T_e$ 帶回 $v_e=\sqrt{2c_p(T_0-T_e)}$,對水蒸氣產物($c_p\approx 2500\ \text{J/(kg·K)}$ 量級)估算,$v_e$ 落在 $3\text{–}4\ \text{km/s}$,與真實氫氧引擎的數據相符。
化學火箭撞上的牆:為什麼 $v_e$ 卡在 4.5 km/s
把前面的物理串起來,我們就能正面回答開篇的問題。排氣速度的上限
$$v_{e,\max} = \sqrt{2\,c_p\,T_0}\quad(\text{完全膨脹,}T_e\to 0)$$
裡的 $T_0$ 並非可以任意調高——它由推進劑的燃燒焓決定,是化學鍵重組能量密度的天花板。氫氧反應每公斤產物釋放的能量是固定的,再怎麼提高燃燒室壓力,也只是讓壅塞流率變大、推力變大,並不會提高每單位質量的能量,因此不會突破 $v_e$ 上限。化學能量密度這道牆,把所有化學火箭的 $v_e$ 鎖死在約 $4.5\ \text{km/s}$ 上下。
這就是火箭方程的「暴政」在更深層次的根源:不是工程不夠努力,而是化學鍵裡能裝的能量就那麼多。要真正鬆綁火箭方程,唯一的出路是換掉能量來源——這就把我們帶向電推進與核推進。
電推進:把能量與推進劑分開
化學火箭有個隱藏的束縛:推進劑同時是能量來源,也是被拋出的工作物質。能量密度上限因此直接卡死排氣速度。電推進(electric propulsion)打破了這個綁定——它用太陽能板或核反應爐外部供能,再用電場把推進劑(通常是氙、氪等惰性氣體)加速到極高速度。能量與工質分開,$v_e$ 就不再被化學鍵束縛。
以離子推進器(ion thruster)為例,先把氙原子電離成正離子,再用電位差 $U$ 加速。一個帶電荷 $q$、質量 $m_i$ 的離子,獲得的動能等於電功,由此解出排氣速度:
$$qU = \frac{1}{2}m_i v_e^2 \;\Rightarrow\; v_e = \sqrt{\frac{2qU}{m_i}}$$
代入氙離子($m_i\approx 2.18\times10^{-25}\ \text{kg}$,單電荷)與加速電壓 $U=1200\ \text{V}$:
$$v_e = \sqrt{\frac{2\times 1.6\times10^{-19}\times 1200}{2.18\times10^{-25}}} \approx 4.2\times10^{4}\ \text{m/s}$$
排氣速度約 $42\ \text{km/s}$,比衝高達 $4300$ 秒——比化學火箭高了一個數量級。火箭方程裡那個指數,瞬間從窒息變得溫和。
推力與功率的鐵則:天下沒有白吃的午餐
那為什麼不全用電推進?因為它有個無法迴避的代價,藏在能量與動量的關係裡。噴流帶走的功率(噴流動能率)是
$$P_{jet} = \frac{1}{2}\dot{m}\,v_e^2 = \frac{1}{2}F\,v_e$$
把它反過來解推力:
$$F = \frac{2P_{jet}}{v_e} = \frac{2\,\eta\,P_{in}}{v_e}$$
$\eta$ 是電轉動能的效率。這條式子是電推進的核心鐵則:在功率 $P_{in}$ 受限時,排氣速度越高,推力反而越小。太陽能板能供應的功率是有限的(典型深空探測器只有數千瓦到數十千瓦),於是離子推進器的推力小得可憐——一具典型離子引擎的推力只有數十到數百毫牛頓,大約等於一張紙壓在手上的重量。它無法把任何東西推離地面,只能在太空中持續、安靜地推上數月甚至數年,靠時間累積出巨大的 $\Delta v$。
這就帶出電推進的真正權衡。給定總任務、固定的供電功率與工作時間,比衝並非越高越好。$v_e$ 太低(接近化學火箭),燃料消耗大;$v_e$ 太高,推力太小、加速太慢,任務時間拖太長,而且供電系統(太陽能板、反應爐)本身的質量會反過來吃掉酬載。存在一個讓酬載最大化的最佳比衝,這就是火箭設計裡著名的「功率受限」最佳化問題——化學火箭是「能量受限」,電推進是「功率受限」,兩者面對的是同一個火箭方程,卻被不同的物理瓶頸卡住。
重點回顧
- 火箭的排氣速度由能量守恆決定:$v_e=\sqrt{2c_p(T_0-T_e)}$。氣體在噴嘴裡膨脹得越冷越徹底、產物分子量越低,$v_e$ 越高。
- 拉瓦爾噴嘴先收斂再發散,是因為氣體要從次音速加速到超音速,幾何上必須在最窄的喉部達到 $M=1$(壅塞),之後在發散段持續加速。
- 引擎效能可乾淨拆成只反映燃燒的特徵速度 $c^*$ 與只反映噴嘴的推力係數 $C_F$,$v_e=c^*C_F$,讓燃燒與噴嘴的問題分別診斷。
- 固定形狀噴嘴無法對所有高度都最佳:海平面易過度膨脹、高空易膨脹不足,這決定了第一節與上面節噴嘴的胖瘦差異。
- 化學火箭的 $v_e$ 被燃燒焓鎖在約 $4.5\ \text{km/s}$;電推進把能量與工質分離,$v_e$ 可達數十 km/s,但受 $F=2\eta P/v_e$ 的功率鐵則限制,推力極小,只能在太空中慢慢累積 $\Delta v$。
深入探討(研究所視角)
凍結流與平衡流:化學動力學如何偷走比衝
前面把噴嘴當成單一氣體的等熵膨脹,但真實燃燒產物是一鍋處於化學平衡的混合物——$\text{H}_2\text{O}$、$\text{OH}$、$\text{H}$、$\text{O}$、$\text{H}_2$ 等物種隨溫度動態解離與復合。氣體在噴嘴裡急速冷卻時,這些反應來不及跟上,產生兩種極限情境:
- 平衡流(equilibrium flow):假設反應無限快,物種組成隨溫度瞬時調整。解離的自由基在降溫時放熱復合,把化學能補回平移動能,給出 $v_e$ 的上限。
- 凍結流(frozen flow):假設反應無限慢,組成在喉部就「凍結」不變。復合放出的能量被鎖在化學鍵裡帶走,給出 $v_e$ 的下限。
真實流動介於兩者之間,取決於膨脹的時間尺度與反應的時間尺度之比(一個 Damköhler 數的問題)。高燃燒室壓力會壓低解離度、使流動更接近平衡流,這是提高燃燒室壓力能小幅提升比衝的一個微妙原因(與單純的壅塞流率無關)。準確預測這個效應需要耦合求解化學動力學速率方程與一維可壓縮流,這正是 NASA CEA 等熱化學程式的核心工作。
邊界層、放熱與真實噴嘴的 99 分
理想等熵分析給出的是滿分答案,真實噴嘴拿不到。主要扣分項有三:邊界層摩擦(壁面附近的黏滯損失,且使有效流通面積略縮)、散度損失(氣流並非完全軸向,發散段使部分動量指向側向,鐘形噴嘴的最佳化即是在最小化此損失與縮短噴嘴長度之間取捨)、以及有限速率化學的凍結損失。工程上用一連串小於 1 的效率因子修正理想值。值得一提的是再生冷卻(regenerative cooling):把低溫燃料先繞過噴嘴與燃燒室壁帶走熱再送進燃燒室,這不只是散熱,被預熱的燃料還把熱能帶回了循環,是少數能同時解決熱防護與微幅提升性能的設計。
超越化學:核熱、磁漿與抵達火星的算術
要同時要大推力與高比衝,化學的能量牆必須被打破。核熱火箭(nuclear thermal rocket, NTR)用核分裂反應爐加熱純氫,繞過了燃燒焓的限制——因為 $v_e\propto\sqrt{T_0/\mathcal{M}}$,純氫的 $\mathcal{M}=2$ 遠低於水蒸氣的 18,即使 $T_0$ 受材料熔點限制(約 $2700\ \text{K}$),$I_{sp}$ 仍可達 $850\text{–}1000$ 秒,約為化學的兩倍,而推力仍維持在可觀的數量級。更前沿的磁漿動力火箭(如 VASIMR)用無線電波把電漿加熱到數百萬度,再以磁噴嘴排出,理論上能在飛行中調節比衝與推力的平衡——起飛段用低 $I_{sp}$ 高推力、巡航段切換到高 $I_{sp}$ 省燃料,這是固定形狀化學引擎做不到的自由度。
把這些放回火箭方程,差異是任務尺度級的。從近地軌道送往火星轉移軌道約需 $\Delta v\approx 3.6\ \text{km/s}$。用化學上面節($v_e\approx 4.4\ \text{km/s}$),質量比 $e^{3.6/4.4}\approx 2.3$;用核熱($v_e\approx 9\ \text{km/s}$),質量比降到 $e^{3.6/9}\approx 1.5$;用離子推進($v_e\approx 30\ \text{km/s}$),質量比僅 $e^{3.6/30}\approx 1.13$——燃料占比從超過一半降到一成出頭。代價是離子推進需要月計、年計的點火時間,且要拖著沉重的供電系統。火箭推進的全部故事,最終都收束回同一條方程式:你要嘛帶更多燃料、要嘛找到把每公斤工質拋得更快的物理。而從牛頓第三定律、到拉瓦爾噴嘴的壅塞流、再到磁噴嘴裡的電漿,人類三百年來做的,無非是在這條方程的指數裡,一吋一吋地與宇宙尺度的距離討價還價。