ジュール・トムソン効果 シミュレーター 戻る
熱力学シミュレーター

ジュール・トムソン効果 シミュレーター — 絞り膨張による温度変化

van der Waals 実在気体の等エンタルピー絞り膨張から、ジュール・トムソン係数 μ_JT = (2a/(RT) − b)/c_p、温度変化 ΔT ≈ μ_JT·ΔP、反転温度 T_inv = 2a/(Rb) を実時間計算します。温度 T、圧力降下 ΔP、vdW 定数 a, b を操作すると、絞り模式図の色温度と μ_JT(T) 曲線が連動して更新され、T < T_inv の冷却領域と T > T_inv の昇温領域の境界を直感的に理解できます。

パラメータ設定
温度 T
K
圧力降下 ΔP
atm
vdW 定数 a
J·m³/mol²
vdW 定数 b ×10⁻⁵
m³/mol

既定値は窒素 N₂ 相当の van der Waals 定数(a=0.138 J·m³/mol²、b=3.87×10⁻⁵ m³/mol)、c_p=29.1 J/(mol·K) 固定。R=8.314 J/(mol·K)。1 atm=101 325 Pa。

計算結果
ジュール・トムソン係数 μ_JT
温度変化 ΔT
反転温度 T_inv
領域(T と T_inv の比較)
絞り膨張の模式図と温度差の可視化

左=高圧側 (P_1, T_1)/中央=絞り(多孔栓・キャピラリ・膨張弁)/右=低圧側 (P_2, T_2)。背景色は温度で決まり、μ_JT>0 なら右側が青く(冷却)、μ_JT<0 なら赤く(昇温)表示されます。

μ_JT(T) 曲線と反転温度

横軸=温度 T (K, 50–1200)/縦軸=μ_JT (K/atm)。青曲線=(2a/(RT)−b)/c_p。水平赤破線=μ_JT=0(反転温度 T_inv)。黄丸=現在の動作点。μ_JT>0 領域=冷却、μ_JT<0 領域=昇温。

理論・主要公式

実在気体の等エンタルピー絞り膨張でジュール・トムソン係数 μ_JT は次のように定義されます:

$$\mu_{JT} = \left(\frac{\partial T}{\partial P}\right)_{\!H}$$

van der Waals 状態方程式 $(P+a/V_m^2)(V_m-b)=RT$ を低密度展開した近似で、c_p をモル等圧比熱とすると:

$$\mu_{JT} \approx \frac{1}{c_{p,m}}\!\left(\frac{2a}{R\,T} - b\right)$$

微小圧力降下に対する温度変化と反転温度は次式:

$$\Delta T \approx \mu_{JT}\cdot \Delta P, \quad T_{\mathrm{inv}} = \frac{2a}{R\,b}$$

$a$ は分子間引力の強さ [J·m³/mol²]、$b$ は分子の排除体積 [m³/mol]、$R=8.314$ J/(mol·K)。$T0$(膨張で冷却)、$T>T_{\mathrm{inv}}$ で $\mu_{JT}<0$(膨張で昇温)。本ツールは N₂ 相当の vdW 定数で $T_{\mathrm{inv}}\approx 858$ K を再現し、Linde 液化サイクルの予冷条件を物理的に説明します。

ジュール・トムソン効果 シミュレーターとは

🙋
エアコンの膨張弁って冷媒が「シューッ」と通るだけで冷たくなるんですよね?仕事もしてないのに、なんで温度が下がるんですか?
🎓
それがまさにジュール・トムソン効果だ。実在気体を細い絞り(多孔栓・キャピラリ・膨張弁)に通すと、外に仕事をしないのにエンタルピー一定で温度が変わる。式は μ_JT = (∂T/∂P)_H で、van der Waals 近似だと (2a/(RT) − b)/c_p になる。本ツールの既定値(N₂ 相当、T=300 K, a=0.138, b=3.87e-5, c_p=29.1)で計算すると μ_JT ≈ 0.251 K/atm。100 atm の圧力降下なら ΔT ≈ −25.1 K の冷却だ。理想気体なら H=c_p·T で圧力に依存しないからゼロ。実在気体の分子間引力(a 項)が「分子同士を引き離す仕事」で内部エネルギーを消費するのが冷却の正体だ。
🙋
じゃあ温度が高いと冷えなくなるんですか?スライダーを T=900 K まで上げたら ΔT がプラスになりました…えっ、逆に温度が上がってる?!
🎓
そう、それが反転温度 T_inv = 2a/(Rb)を越えた領域だ。本ツールの既定値だと T_inv ≈ 858 K。T < T_inv では a 項(分子間引力)が支配的で μ_JT > 0(冷える)、T > T_inv では b 項(排除体積)が支配的で μ_JT < 0(昇温する)。だから水素(T_inv ≈ 202 K)は常温で絞っても加熱されてしまう。液化するには液体窒素などで先に T_inv 以下まで予冷してから絞らないといけない。1898 年に Dewar が水素を液化したのも、Linde サイクルでこの予冷工夫を組み合わせたからだ。
🙋
理想気体と何が違うんでしたっけ?理想気体は絞ったら温度変わらない、って習った気がします。
🎓
よく覚えてる。理想気体は H = c_p·T で圧力に依存しないから、等エンタルピー絞りなら ΔT=0。逆に言うと「絞ったら温度が変わる」ことこそが、その気体が理想気体から外れている証拠だ。van der Waals 方程式 (P + a/V²)(V−b) = RT の a, b 項が小さい He や H₂ は反転温度が低く(He は 45 K!)、a, b が大きい CO₂ や水蒸気は反転温度が高い(CO₂ で約 1500 K)。室温で「絞れば冷える」と直感できるのは、実は N₂, O₂, CO₂, R32 などごく限られた気体だけだ。
🙋
なるほど〜。あと「ΔP は負の値」というのも変な気がします。スライダーが −500 atm まで動くのは、圧力降下を負で表現してるんですか?
🎓
そうだ。ΔP = P_2 − P_1 と定義すると、膨張(高圧→低圧)では P_2 < P_1 だから ΔP < 0。式 ΔT ≈ μ_JT · ΔP に直接放り込めるよう、本ツールでは負の値を採用している。例えば「圧縮機 100 atm → 大気圧」なら ΔP ≈ −99 atm。スライダーを 0 から −100 まで動かすと、線形に ΔT が変化するのが見えるはずだ。実機の膨張弁は数 atm 〜 数十 atm の降下で運用されるが、Linde サイクルの初段は 200 atm から大気圧近くまで一気に絞り、ΔT ≈ −50 K 級の冷却で液化温度に近づける。
🙋
右下のμ_JT(T) 曲線で、横軸を T_inv が横切るのを見てから、もう一度反転温度の意味が腑に落ちました。a と b のスライダーで気体の種類を切り替えると T_inv 位置が移動するのも面白いです。
🎓
いいね。例えば a = 0.366(CO₂ 相当), b = 4.29e-5 にすると T_inv ≈ 2050 K、ヘリウム相当(a = 0.0035, b = 2.4e-5)なら T_inv ≈ 35 K。本ツールでこれらを試してみると「なぜ He は液化が難しいか」「なぜ CO₂ は室温の絞り膨張だけでドライアイスができるか」が一発で見える。反転温度の概念を持つだけで、極低温工学・冷凍機・天然ガス液化(LNG)の本質が同じ枠組みで理解できるのがジュール・トムソンの偉大なところだ。

よくある質問

ジュール・トムソン効果は、実在気体を細い絞り(throttle)や多孔栓を通して断熱的に膨張させたとき、温度が変化する現象です。理想気体ではエンタルピー一定で温度は変わりませんが、分子間引力(vdW 定数 a)と排除体積(b)を持つ実在気体では μ_JT = (∂T/∂P)_H ≠ 0 になります。van der Waals の近似式では μ_JT ≈ (2a/(RT) − b)/c_p と表され、本ツールの既定値(T=300 K, a=0.138 J·m³/mol², b=3.87×10⁻⁵ m³/mol, c_p=29.1 J/(mol·K))では μ_JT ≈ 0.251 K/atm。100 atm 圧力降下で ΔT ≈ −25.1 K の冷却が得られます。気体液化(Linde サイクル)や家庭用エアコンの冷媒制御弁の基本原理です。
反転温度 T_inv = 2a/(Rb) は μ_JT の符号が反転する境界温度です。T < T_inv では分子間引力の影響が支配的で 2a/(RT) > b、μ_JT > 0 となり膨張で冷却します。T > T_inv では排除体積の影響が支配的で μ_JT < 0 となり、膨張で逆に温度が上がります。本ツール既定値では T_inv ≈ 858 K。窒素 N₂ の実測 T_inv は約 621 K、酸素は 764 K、水素はわずか 202 K と低く、室温(300 K)の水素は膨張で加熱されるため、液化にはまず予冷が必要です。He は T_inv ≈ 45 K で、より厳しい条件が要求されます。
蒸気圧縮冷凍サイクルでは、高圧液冷媒(R32, R410A 等)が膨張弁(キャピラリチューブ)を通って蒸発器に入る瞬間、ジュール・トムソン絞り膨張で温度が急降下します。これにより蒸発器内で大気から熱を吸収し、室内を冷やす仕組みです。冷媒は T < T_inv の領域で使われるよう設計されており、典型的に 10 atm → 3 atm の圧力降下で 15〜30 K の冷却が起きます。本ツールで ΔP スライダーを 0 から −100 atm に動かすと、μ_JT × ΔP に比例した温度降下が線形に表示され、膨張弁設計の感度を体感できます。
断熱膨張(タービン・ピストン)はガスが外部に仕事をして内部エネルギーを失うため、理想気体でも温度降下が起きます(Pv^γ=const)。一方ジュール・トムソン絞り膨張は等エンタルピー(H=const)で、外部に仕事をしない代わりに分子間引力に対抗するための内部エネルギーを消費して冷却します。理想気体では H = c_p·T のみで圧力に依存しないため絞り膨張で温度は変わりません。実在気体で初めて μ_JT ≠ 0 となるのが本質的な違いです。Linde の液化サイクル(1895)はこの効果を活用、Claude サイクル(1902)は断熱膨張機を組み合わせ効率を改善しました。

実世界での応用

家庭用エアコン・冷蔵庫の膨張弁:蒸気圧縮冷凍サイクルでは凝縮器を出た高圧液冷媒(R32 で 30〜40 atm)が膨張弁またはキャピラリチューブを通って蒸発器(3〜10 atm)に入ります。この絞り膨張で冷媒は気液混合の二相流となり、温度が −5〜−25 °C まで急降下。蒸発器で室内空気から熱を吸い取って気化し、再び圧縮機に戻ります。本ツールで R32(a ≈ 0.50, b ≈ 6.5e-5 SI 換算)を再現するには a, b スライダーを調整してみると、N₂ より μ_JT が大きく、より少ない圧力降下で大きな冷却が得られることが分かります。

天然ガス液化プラント(LNG)と空気液化:LNG 工場では天然ガス(主成分 CH₄, a ≈ 0.228, b ≈ 4.28e-5)を 70〜80 atm に圧縮し、ジュール・トムソン弁で大気圧近くまで一気に絞り −162 °C の液化温度に到達させます。Linde-Hampson プロセス(1895)が原型で、現代の大規模 LNG プラント(カタール RasGas, 米国 Sabine Pass 等)も基本構造は同じ。空気液化では Linde 単段では効率が悪いため、Claude サイクル(1902, 断熱膨張タービン併用)や Heylandt サイクルが採用されています。

クライオスタット・量子コンピュータ用希釈冷凍機:超伝導量子ビットを動作させる希釈冷凍機(dilution refrigerator)の予冷段では、液体窒素(77 K)→ 液体ヘリウム(4.2 K)→ ³He/⁴He 希釈段(数 mK)と多段冷却します。Helium は T_inv が約 45 K と極めて低いため、液体窒素で予冷してから初めてジュール・トムソン弁で液化可能になります。本ツールで a を 0.0035 まで下げ、b を 2.4e-5 にすると He 相当の T_inv ≈ 35 K が再現でき、なぜ He の液化が 1908 年(Kamerlingh Onnes)まで待たねばならなかったかが分かります。

燃料電池車・水素ステーションの充填技術:水素は T_inv ≈ 202 K と低いため、常温(300 K)で圧縮水素タンクを急速充填すると逆に温度が上がる(「充填昇温」)。70 MPa(約 700 atm)の高圧タンクへ短時間で充填すると 80〜90 °C まで温度上昇し、タンク強度低下のリスクがあるため、SAE J2601 規格では充填前のプリクーリング(−40 °C)が義務化されています。本ツールで T=300 K, a=0.0247(H₂ 相当), b=2.66e-5 にすると μ_JT < 0 となり、水素特有の「絞ったら温まる」挙動が確認できます。

よくある誤解と注意点

最も典型的な誤解は、「ガスを絞れば必ず冷える」と考えてしまうことです。実際には温度と気体種に依存し、T > T_inv では逆に加熱されます。水素・ヘリウム・ネオンは室温で「絞っても冷えない」気体の代表例で、これらを液化するには事前に予冷が必要です。Linde が水素液化に成功したのも、彼が事前に液体窒素で水素を 77 K まで冷却しておいた工夫があったからです。本ツールで T=300 K にして a=0.025(H₂ 相当)に下げると、T_inv が下がって μ_JT < 0 になり「冷えない水素」を体験できます。

次に、「ジュール・トムソン効果は仕事をしている」と勘違いするケースです。実際には等エンタルピー過程で、外部仕事はゼロ。冷却の源は内部エネルギー(分子間引力を引き離すための位置エネルギー)です。「仕事を伴う断熱膨張(タービン)」と「仕事を伴わない絞り膨張」を混同しないこと。前者は理想気体でも温度降下が起き効率が高い(Claude サイクル)、後者は実在性に依存しエネルギー効率は低いがコンパクトで信頼性が高い(Linde サイクル)。家庭用エアコンが膨張タービンでなく膨張弁を使うのは、可動部の少ない簡便さが優先されるからです。

最後に、「van der Waals 近似で全領域カバーできる」という過信。実際には低温・高圧域(液化点近傍)では vdW 近似は精度不足で、Redlich-Kwong、Peng-Robinson、Lee-Kesler 等の修正状態方程式が必要です。さらに本ツールの μ_JT 式は「低密度展開の主要項」のみで、完全形は μ_JT = (T(∂V/∂T)_P − V) / c_p で計算します。GERG-2008 や NIST REFPROP 等の精密データベースでは多変数フィットで再現精度を確保。本ツールは「反転温度の存在と符号変化」を理解するための教育ツールで、実機プロセス設計には専用ソフト(HYSYS, Aspen Plus)を併用してください。