動画解説はこちら|YouTube

実在気体|ファン・デル・ワールス式・圧縮係数・対臨界状態・Kayの方法

※本記事の参考文献には、Amazonアソシエイトのリンク(アフィリエイト広告)を含みます。

困っている人
困っている人

高圧の容器に何kg入るのか、理想気体で計算したら現物と合わなかった。

圧縮係数の図を見せられても、どの線を読めばいいのか分からない。

そんな悩みを解決します。

この記事の内容

10 MPa で15%ずれる ── 高圧ガスで理想気体が使えない理由

圧縮係数 Z を読む手順は対臨界に直すだけ

vdW式の a は分子間力、b は分子自身の体積

混合気体は Kayの方法(モル分率で加重平均)

高圧ガス設備の充填量計算、貯槽の在庫管理、安全弁のサイジング。実務で効いてくる補正です。

現場で、この補正はどこで効くか

理想気体で計算すると危険側に外れる場面があります。

場面何が起きるか
高圧容器の充填量理想気体の計算より実際は多く入る(Z<1のとき)
貯槽の在庫管理圧力から量を出すと誤差が出る
安全弁のサイジング吹出量の計算に密度が必要
圧縮機の動力計算体積効率の評価にZが要る
低温での挙動臨界点に近づくと理想気体はまったく使えない

Z < 1 なら「理想気体の計算より多く入る」。 過充填のリスクに直結します。

常温・常圧なら理想気体で十分です。 補正が要るのは高圧・低温のときだけ ── その線引きを知ることが実務では重要です。

理想気体の基本は理想気体の状態方程式と気体定数を参照してください。

なぜ理想気体からずれるのか

理想気体の2つの仮定

① 分子自身の体積はゼロ

② 分子間に力が働かない

どちらも現実には成り立ちません。 そして2つの効果は逆向きに働きます。

効果何が起きるか圧力への影響
分子間力(引力)壁にぶつかる分子が引き戻される実測圧力が下がる → Z < 1
分子自身の体積動ける空間が狭くなる圧力が上がる → Z > 1

だから Z のグラフは下がってから上がる

条件支配的な効果Z
高温・低圧どちらも小さいZ ≒ 1
低温・中圧分子間力Z < 1(最も縮む)
高圧分子自身の体積Z > 1(縮まない)

圧力を上げていくと、Zは一度下がって最小値を通り、また上がって1を超えます。

ファン・デル・ワールスの式

(p + a n²/V²)(V − nb) = nRT

1 mol あたりなら (p + a/Vm²)(Vm − b) = RT

何の補正かどちら向きに効くか
a n²/V²分子間力測定圧力に足す(本来の圧力は高い)
nb分子自身の体積(排除体積)体積から引く(動ける空間は狭い)

a は分子間力の強さ、b は1 mol の分子が占める実効体積です。

物質a[Pa·m⁶/mol²]b[m³/mol]特徴
He0.00352.4×10⁻⁵a が極小(分子間力がほとんどない)
N₂0.1413.9×10⁻⁵──
CO₂0.3644.3×10⁻⁵a が大きい(分極しやすい)
NH₃0.4233.7×10⁻⁵a が大きい(水素結合)
H₂O0.5543.1×10⁻⁵a が最大(強い水素結合)

極性が強いほど a が大きい。 アンモニアや水蒸気は理想気体からのずれが大きい ── 実務ではこれが効きます。

ヘリウムは a がほぼゼロ。だからヘリウムは非常に理想気体に近く、極低温まで液化しません。

圧縮係数 Z ── ずれを1つの数で表す

定義

Z = pV / (nRT)

pV = ZnRT

理想気体なら Z = 1

対応状態の原理 ── これが図表を使える理由

ガスの種類が違っても、臨界点からの相対的な位置が同じなら、Z がほぼ同じになります。

対臨界状態(換算状態)

Tr = T / Tc(対臨界温度)

pr = p / pc(対臨界圧力)

だから1枚の図表(一般化圧縮係数線図)で、すべてのガスに対応できます。

主なガスの臨界定数

ガスTc[K]Tc[℃]pc[MPa]
He5.2−2680.23
H₂33.2−2401.30
N₂126.2−1473.39
O₂154.6−1195.04
CH₄190.6−82.64.60
C₂H₆305.432.24.88
CO₂304.231.07.38
C₃H₈369.896.74.25
NH₃405.5132.411.35
H₂O647.1374.022.06

臨界温度が実務の分水嶺です。

臨界温度の意味

臨界温度より上では、どれだけ加圧しても液化しない

→ N₂・O₂・H₂・He は常温で圧縮ガス(液化しない)

→ C₃H₈・NH₃ は常温で液化ガス(加圧すれば液になる)

CO₂ の臨界温度は31℃。夏場は液化しにくく、冬場は液化する ── 実際に季節で挙動が変わります。

Z を読む手順

  1. Tr = T/Tcpr = p/pc を計算する
  2. 図表で Tr の曲線を選ぶ
  3. pr の位置で縦軸の Z を読む
  4. pV = ZnRT に代入する

数値で確かめる

例1:高圧窒素の充填量

内容積 47 L の容器に、15 MPaA・25℃ で窒素を充填した。
理想気体で計算した場合と、Z = 1.02 で補正した場合の質量は?

〔解答〕

Tr = 298.15 / 126.2 = 2.36(臨界温度よりはるか上)

pr = 15 / 3.39 = 4.42

→ 図表から Z ≒ 1.02(高圧なので Z > 1)

理想気体:n = 15×10⁶ × 0.047 / (8.314 × 298.15) = 284.5 mol → 7.97 kg

補正あり:n = 284.5 / 1.02 = 278.9 mol7.81 kg

Z > 1 なので、理想気体の計算より実際は少ししか入りません。 差は2%。

例2:低温の CO₂ ── ここは差が大きい

CO₂ を 5 MPaA・20℃ で扱う。理想気体との差は?

〔解答〕

Tr = 293.15 / 304.2 = 0.96(臨界温度のすぐ下)

pr = 5 / 7.38 = 0.68

→ この領域では Z ≒ 0.3〜0.4

理想気体の2〜3倍の量が入る

臨界点付近では理想気体はまったく使えません。 3倍近く外れます。

この領域では実際には液化していることも多く、状態そのものを確認する必要があります。

例3:Kayの方法(混合気体)

仮臨界定数

Tc,mix = Σ(yᵢ Tc,ᵢ)

pc,mix = Σ(yᵢ pc,ᵢ)

yᵢ はモル分率

メタン70 mol%、エタン30 mol% の混合気体を、250 K・6.0 MPaA で扱う。Z を読む座標は?

〔解答〕

Tc,mix = 0.70×190.6 + 0.30×305.4 = 133.4 + 91.6 = 225.0 K

pc,mix = 0.70×4.60 + 0.30×4.88 = 3.22 + 1.46 = 4.68 MPa

Tr = 250 / 225.0 = 1.11

pr = 6.0 / 4.68 = 1.28

Kayの方法は近似です。 成分の性質が大きく違う(分子サイズや極性の差が大きい)と誤差が出ます。

実務で精度が要るときは、状態方程式(Peng-Robinson、SRKなど)をプロセスシミュレータで解きます。Kayは手計算での概算用です。

実務での使いどころと、失敗例

使いどころ:どこから補正するか

実務的な目安

pr < 0.1(臨界圧力の1割以下)→ 理想気体でよい

Tr > 2(臨界温度の2倍以上)→ 高圧でも Z は1に近い

Tr < 1.2 かつ pr > 0.5必ず補正する

窒素・酸素・空気は常温なら Tr が2以上なので、10 MPa 程度までは誤差数%で済みます。

プロパン・アンモニア・CO₂ は常温で Tr が1前後。ここは要注意です。

失敗例:貯槽の在庫を圧力で管理する

何が起きるか

気相部の圧力から在庫を推定する運用にしていた

外気温が下がって Tr が下がる → Z が小さくなる

同じ圧力でも実際の物質量は増えている

在庫を過小評価してしまう

季節で誤差が変わります。 温度補正なしの圧力管理は、液化ガスでも圧縮ガスでも危険です。

失敗例:安全弁のサイジングで密度を間違える

安全弁の必要吹出面積は、吹出条件での密度に依存します。

理想気体で計算すると

Z < 1 の領域では実際の密度は理想気体より大きい

→ 同じ質量流量でも体積流量は小さい

→ 理想気体で計算すると過大な面積になる(安全側だが不経済)

逆に Z > 1 の領域では過小になります。 こちらは危険側です。

安全弁の設計は安全弁の設計と口径計算、設置基準は安全弁の設置基準を参照してください。

超臨界流体という領域

臨界温度・臨界圧力の両方を超えると、気体でも液体でもない状態になります。

性質超臨界流体は
密度液体に近い(溶解力が大きい)
粘度・拡散性気体に近い(浸透しやすい)
界面存在しない(気液の区別がない)

超臨界CO₂による抽出(コーヒーのカフェイン除去など)はこの性質を使っています。

プラントで意図せず超臨界に入ることもあります。 気液分離を前提にした設計が成り立たなくなるので、運転条件の確認が必要です。

フガシティー ── 平衡計算のための補正

Z が体積の補正なら、フガシティーは圧力(化学ポテンシャル)の補正です。

定義

f = φ p(φ:フガシティー係数)

活量 a = f / f°

理想気体なら φ = 1、f = p

圧縮係数 Zフガシティー係数 φ
補正の対象体積・密度化学ポテンシャル
使う場面充填量・在庫・流量化学平衡・気液平衡
理想気体Z = 1φ = 1

高圧下の平衡計算では、分圧ではなくフガシティーを使わないと合いません。 アンモニア合成のような高圧反応では必須です.

平衡定数を活量で定義しておくと、理想気体でも実在気体でも同じ形の式が使えます。

平衡定数の温度依存性学識・型6|化学平衡

甲種化学ではこう出ます

学識の型2にあたります。計算としては令和2〜6年度の60大問中1回ですが、記述では3回出ています。

年度・問内容
R4 国試 問1Kayの方法 → 対臨界 → 圧縮係数 → 質量計算
R2 国試 問6(1)vdW式を分子間力と排除体積で説明記述
R3 検定 問6(1)vdW式と2つの補正項記述
R4 検定 問6(1)理想気体とは何か記述

記述で「排除体積」「分子間力」という語を書けるかが問われます。

解き方と答案例は学識・型2|実在気体、記述対策全般は学識 問6の記述問題にまとめています。

ガスごとの臨界定数はガス各論①ガス各論②でも扱っています。

まとめ

  • 理想気体の仮定は①分子体積ゼロ ②分子間力なし。両方とも現実には成り立たない
  • 2つの効果は逆向き。だから Z は下がってから上がる
  • vdW式の a は分子間力、b は分子自身の体積(排除体積)
  • 極性が強いほど a が大きい(NH₃・H₂O はずれが大きい)
  • Tr = T/Tc、pr = p/pc に直して図表を引く
  • 臨界温度より上では液化しない。N₂・O₂・H₂ は常温で圧縮ガス
  • Tr < 1.2 かつ pr > 0.5 なら必ず補正
  • 混合気体はKayの方法(モル分率で加重平均、近似)
  • Z は体積、フガシティー φ は平衡の補正

気体の溶解はヘンリーの法則とガスの溶解度、状態変化は熱力学第一法則と4つの状態変化へ。

基本に戻る場合は理想気体の状態方程式と気体定数を参照してください。