二維相圖上的物理審判：徒手繪製焦耳-湯姆遜反轉曲線
在上一篇博文《絕熱節流膨脹的熱力學帳本》中，我們利用麥克斯威關係式，優雅地推導出了焦耳-湯姆遜係數的古典熱力學終點公式：
\[\mu_{\text{JT}} = \left(\frac{\partial T}{\partial P}\right)_H = \frac{1}{C_p}\left[T\left(\frac{\partial V}{\partial T}\right)_P - V\right]\]
我們也揭開了理想氣體 \(\mu_{\text{JT}} = 0\) 的宿命。然而，真實氣體的世界從來不是非黑即白的。同一種氣體，在某個溫度壓力下外洩會自發結霜，在另一種極端高壓或高溫下外洩，卻可能反常地飆溫引爆。
在大自然這本複雜的帳本裡，究竟是誰在決定引力與斥力的統治疆界？今天，我們將暫時合上微觀分子動力學的視角，純粹用古典熱力學的偏微分重錘與凡德瓦狀態方程（vdW EOS），在 \(P-T\) 二維平面上，徒手勾勒出那條劃分氣體生與死、冷與熱的臨界邊界——焦耳-湯姆遜反轉曲線（Joule-Thomson Inversion Curve）。
一、 幾何邊界的誕生：何謂反轉狀態？
所謂的反轉狀態（Inversion State），指的是真實氣體在經歷等焓節流膨脹時，其巨觀溫度既不升高也不降低的物理分水嶺。也就是說，此時的溫度隨壓力的微擾變動率為零：
\[\mu_{\text{JT}} = 0 \implies T\left(\frac{\partial V}{\partial T}\right)_P - V = 0\]
如果在以壓力（\(P\)）為橫軸、溫度（\(T\)）為縱軸的二維坐標系中，將所有滿足 \(\mu_{\text{JT}} = 0\) 的狀態點全部連起來，就會勾勒出一條宛如拋物線般的封閉曲線。
這張熱力學幾何地圖，將真實氣體劃分為兩個截然不同的口袋：
- 曲線內部（陰影包絡區）：\(\mu_{\text{JT}} > 0\)。此區間內壓降必然導致降溫（\(dP < 0 \implies dT < 0\)）。這是「引力佔優勢區」，也是工業上低溫液化空氣的黃金戰場。
- 曲線外部（荒漠區）：\(\mu_{\text{JT}} < 0\)。此區間內壓降反而導致升溫（\(dP < 0 \implies dT > 0\)）。這是「硬核斥力佔優勢區」，也就是高壓氫氣外洩自燃的熱力學禁區。
二、 不跳步的硬核代數：用凡德瓦方程定量求解
在化奧決賽或冬令營的考卷上，最經典的魔王題莫過於：「嘗試用凡德瓦方程，推導出反轉曲線的 \(P-T\) 幾何方程。」許多教科書在這裡都會跳步，今天我們就一步步把這筆代數帳算得清清楚楚。
凡德瓦狀態方程形式如下：
\[\left(P + \frac{a}{V_m^2}\right)(V_m - b) = RT \quad \text{--- (式1)}\]
我們的目標是求出 \(\left(\frac{\partial V_m}{\partial T}\right)_P\)。由於這是關於 \(V_m\) 的三次方程，直接解出 \(V_m(T,P)\) 再求導極其愚蠢。我們必須祭出隱函數微分法。
保持壓力 \(P\) 固定（\(dP = 0\)），對（式1）兩邊同時進行關於溫度 \(T\) 的全微分：
\[\left(-\frac{2a}{V_m^3} dV_m\right)(V_m - b) + \left(P + \frac{a}{V_m^2}\right) dV_m = R dT\]
同除以 \(dT\) 並限制等壓條件，提取出 \(\left(\frac{\partial V_m}{\partial T}\right)_P\)：
\[\left(\frac{\partial V_m}{\partial T}\right)_P \left[ \left(P + \frac{a}{V_m^2}\right) - \frac{2a(V_m - b)}{V_m^3} \right] = R \quad \text{--- (式2)}\]
此時，我們帶入反轉曲線的臨界邊界條件：\(\left(\frac{\partial V_m}{\partial T}\right)_P = \frac{V_m}{T}\)。將其帶入（式2）：
\[\frac{V_m}{T} \left[ \left(P + \frac{a}{V_m^2}\right) - \frac{2a(V_m - b)}{V_m^3} \right] = R \implies \left(P + \frac{a}{V_m^2}\right) - \frac{2a(V_m - b)}{V_m^3} = \frac{RT}{V_m} \quad \text{--- (式3)}\]
為了消去棘手的溫度 \(T\)，我們利用原本的凡德瓦方程（式1），將 \(\frac{RT}{V_m}\) 項目展開：
\[\frac{RT}{V_m} = \frac{\left(P + \frac{a}{V_m^2}\right)(V_m - b)}{V_m} = \left(P + \frac{a}{V_m^2}\right) - \frac{b}{V_m}\left(P + \frac{a}{V_m^2}\right) \quad \text{--- (式4)}\]
將（式4）代回（式3）的等號右側：
\[\left(P + \frac{a}{V_m^2}\right) - \frac{2a(V_m - b)}{V_m^3} = \left(P + \frac{a}{V_m^2}\right) - \frac{b}{V_m}\left(P + \frac{a}{V_m^2}\right)\]
兩邊同時消去共通項 \(\left(P + \frac{a}{V_m^2}\right)\)，並移除負號：
\[\frac{2a(V_m - b)}{V_m^3} = \frac{b}{V_m}\left(P + \frac{a}{V_m^2}\right)\]
兩邊同乘以 \(V_m^3\) 展開並移項，我們便得到了反轉曲線上壓力 \(P\) 與摩爾體積 \(V_m\) 的幾何關係式：
\[2a(V_m - b) = b V_m^2 \left(P + \frac{a}{V_m^2}\right) = b P V_m^2 + a b\]
\[ b P V_m^2 = 2a V_m - 3a b \implies { \color{red} P = \frac{2a}{V_m b} - \frac{3a}{V_m^2}} \quad \text{--- (終點公式一)}\]
為了在 \(P-T\) 相圖上作圖，我們需要將（終點公式一）再代回原凡德瓦方程（式1）中消去 \(P\)，即可解鎖溫度的參數方程：
\[RT = \left[ \left(\frac{2a}{V_m b} - \frac{3a}{V_m^2}\right) + \frac{a}{V_m^2} \right] (V_m - b) = \left( \frac{2a}{V_m b} - \frac{2a}{V_m^2} \right) (V_m - b)\]
\[RT = \frac{2a}{V_m^2 b} (V_m - b)(V_m - b) = \frac{2a}{b} \left( \frac{V_m - b}{V_m} \right)^2\]
\[{\color{red} T = \frac{2a}{Rb} \left(1 - \frac{b}{V_m}\right)^2} \quad \text{--- (終點公式二)}\]
透過（終點公式一）與（終點公式二），我們成功用微觀特徵參數（\(a\) 與 \(b\)），聯立鎖定了這條二維幾何拋物線的全部軌跡。
三、 拋物線的靈魂：三個決定性的幾何極限點
這條拋物線在相圖上的特徵位置，直接決定了氣體的物理極限。我們可以透過微積分對這組參數方程進行邊界審判：
1. 最高反轉溫度（Maximum Inversion Temperature, \(T_{\text{inv, max}}\)）
當氣體趨近於極度稀薄的理想狀態時，摩爾體積大到邊界（\(V_m \to \infty\)）。此時由（終點公式一）可知壓力 \(P \to 0\)。將其帶入（終點公式二）：
\[T_{\text{inv, max}} = \lim_{V_m \to \infty} \frac{2a}{Rb} \left(1 - \frac{b}{V_m}\right)^2 = {\color{red}\frac{2a}{Rb}}\]
這就是真實氣體的冷卻天花板。如果一隻氣體的環境溫度高於這個 \(T_{\text{inv, max}}\)，它體內的熱動能將徹底撕裂引力井。不論你如何施加高壓，它在洩漏膨脹時百分之百只會自發飆溫，絕不可能降溫。
2. 最低反轉溫度（Minimum Inversion Temperature, \(T_{\text{inv, min}}\)）
拋物線跌回橫軸 \(P = 0\) 的另一個低溫交點。令（終點公式一）中的 \(P = 0\)：
\[\frac{2a}{V_m b} - \frac{3a}{V_m^2} = 0 \implies V_m = \frac{3}{2}b\]
將其帶回（終點公式二）：
\[T_{\text{inv, min}} = \frac{2a}{Rb} \left(1 - \frac{b}{\frac{3}{2}b}\right)^2 = \frac{2a}{Rb} \left(1 - \frac{2}{3}\right)^2 = { \color{red} \frac{2a}{9Rb}}\]
3. 最大轉折壓力（Maximum Inversion Pressure, \(P_{\text{inv, max}}\)）
這條拋物線在極限右側的頂點。我們對（終點公式一）關於 \(V_m\) 求導並令其為零（\(\frac{dP}{dV_m} = 0\)）：
\[\frac{d}{dV_m}\left(\frac{2a}{V_m b} - \frac{3a}{V_m^2}\right) = -\frac{2a}{V_m^2 b} + \frac{6a}{V_m^3} = 0 \implies V_m = 3b\]
將 \(V_m = 3b\) 帶回壓力公式，得到氣體可以展現冷卻效應的絕對壓力上限：
\[P_{\text{inv, max}} = \frac{2a}{(3b)b} - \frac{3a}{(3b)^2} = \frac{2a}{3b^2} - \frac{3a}{9b^2} = {\color{red}\frac{a}{9b^2}}\]
四、 國手級數據對帳：\(\text{H}_2\) 與 \(\text{N}_2\) 的命運判決
現在，我們把真實世界的凡德瓦常數帶入這三項幾何指標，看看大自然是如何在相圖上狠狠操縱這兩種氣體命命運的。
| 氣體分子 | 凡德瓦 \(a\) (\(\text{L}^2\cdot\text{bar}\cdot\text{mol}^{-2}\)) | 凡德瓦 \(b\) (\(\text{L/mol}\)) | 最高反轉溫度 \(T_{\text{inv, max}}\) | 最大轉折壓力 \(P_{\text{inv, max}}\) |
|---|---|---|---|---|
| 氮氣 (\(\text{N}_2\)) | \(1.370\) | \(0.0387\) | \(851\text{ K } (578^\circ\text{C})\) | \(101.4\text{ bar}\) |
| 氫氣 (\(\text{H}_2\)) | \(0.248\) | \(0.0267\) | \(223\text{ K } (-50^\circ\text{C})\) | \(38.6\text{ bar}\) |
⚖️ 能量會計師的相圖解讀
我們將生活的常溫環境（\(25^\circ\text{C} = 298.15\text{ K}\)）化作這張 \(P-T\) 地圖上的一條水平橫線：
- 氮氣（\(\text{N}_2\)）的城堡守護： 氮氣的最高反轉溫度高達 \(578^\circ\text{C}\)，最低反轉溫度為 \(-178^\circ\text{C}\)。當我們在常溫（\(298\text{ K}\)）下對氮氣進行日常充填與外洩時，這條 \(298\text{ K}\) 的生活橫線，安穩地橫切過氮氣拋物線內部的陰影包絡區。 只要壓力不要極端超越其最大轉折壓力，氮氣外洩的幾何軌跡就是「由右向左」滑動，穩穩地落入降溫區。這就是空壓管線漏氣必然結霜的幾何真相。
- 氫氣（\(\text{H}_2\)）的荒漠宿命： 氫氣因為分子極小、極難極化，凡德瓦 \(a\) 弱到可以忽略，直接導致它的最高反轉溫度封頂在寒冷的 \(-50^\circ\text{C}\)（\(223\text{ K}\)）！ 看清楚了嗎？我們生活的常溫（\(298\text{ K}\)），在幾何相圖上遠遠懸浮在整座氫氣反轉拋物線城堡的天花板之上。
常溫下的氫氣，其初始狀態點從一出發就落在反轉曲線之外的「硬核斥力荒漠」中。不論你把它加壓到 \(350\text{ bar}\) 還是乘用車標準的 \(700\text{ bar}\)，一旦發生洩漏進行等焓膨脹，它在二維相圖上的幾何軌跡，注定只能在拋物線外部的荒漠裡一路狂奔——壓力越跌，溫度越飆，直至點燃自己。
📌 課後思維黑盒子：走向統計力學的伏筆
透過本篇硬核的隱函數代數推導，我們已經不需要任何微觀猜想，純粹用巨觀狀態方程便在 \(P-T\) 相圖上宣告了氫能工程面臨的熱力學宿命。
但作為一名優秀的化奧選手，你的審判思維不該停留在唯象的凡德瓦參數上。你必須進一步追問：凡德瓦方程中那兩個拍腦袋想出來的修正項 \(a\) 與 \(b\)，其數學本質究竟是怎麼從分子間的幾何交互作用勢能中沉澱出來的？第二維里係數與邁爾 \(f\) 函數，又是如何透過統計力學的配分函數，完美收斂到 \(T_{\text{inv}} = \frac{2a}{Rb}\) 這個終點公式的？
下一篇，我們將正式告別巨觀偏微分，解鎖統計力學的終極會計帳本！
Comments
Post a Comment