飽和蒸気圧(Wagner型相関式):
$$ \ln\frac{P_{sat}}{P_c} = \frac{a_1\tau + a_2\tau^{1.5} + a_3\tau^{3} + a_4\tau^{3.5} + a_5\tau^{4} + a_6\tau^{7.5}}{1-\tau},\quad \tau = 1-\frac{T}{T_c} $$湿り蒸気の比エンタルピー・比エントロピー:
$$ h = h_f + x\,h_{fg},\qquad s = s_f + x\,\frac{h_{fg}}{T_{sat}} $$過熱蒸気(簡易近似):
$$ h \approx h_g(T_{sat}) + c_p\,(T - T_{sat}) $$蒸気特性計算機とは
物理モデルと主要な数式
水蒸気の状態を記述する基礎となるのは、状態方程式と熱力学の関係式です。特に、飽和状態(気液平衡)の蒸気圧と温度の関係はAntoine式などの実験式で表されます。
$$ \log_{10}P_{sat}= A - \frac{B}{C + T}$$ここで、$P_{sat}$は飽和蒸気圧[bar]、$T$は温度[°C]、$A, B, C$は物質固有の定数(Antoine定数)です。この式は、シミュレーターで飽和ドームの境界線を計算する基礎となっています。
湿り蒸気領域(飽和ドーム内)の状態を決定する重要なパラメータが乾き度 $x$ です。これを用いて、比エンタルピー $h$ や比エントロピー $s$ などの特性値は、飽和液の値と飽和蒸気の値から線形補間で求められます。
$$ h = h_f + x \cdot h_{fg}$$$h_f$は飽和液の比エンタルピー、$h_{fg}$は蒸発潜熱(飽和蒸気と飽和液のエンタルピー差)です。$x=0$なら全て液体、$x=1$なら全て乾き飽和蒸気となります。シミュレーターでは、ドーム内の点を選ぶとこの乾き度が自動計算されます。
よくある質問
実世界での応用
火力・原子力発電プラント:ランキンサイクルと呼ばれる基本サイクルの熱効率解析に不可欠です。ボイラーで発生させた過熱蒸気をタービンで膨張させる過程をh-s図上で追跡し、取り出せる仕事量や復水器で捨てる熱量を計算します。
化学プラントのプロセス設計:蒸気を反応熱の供給源や蒸留塔の熱源として利用する際、必要な蒸気の温度・圧力・流量を決定します。湿り蒸気の乾き度は、伝熱効率に直接影響するため、厳密な管理が必要です。
船舶・機関車の蒸気機関:往復動式の蒸気機関では、シリンダ内での蒸気の膨張過程を理解し、出力や燃料消費率を予測するために使用されます。歴史的にもモリエ線図は蒸気機関の効率改善に貢献しました。
食品・製薬工業の滅菌工程:高圧蒸気(オートクレーブ)を用いた滅菌では、温度と圧力の関係が殺菌効果を左右します。飽和蒸気表や線図は、確実に滅菌条件を達成するための基準として用いられています。
よくある誤解と注意点
まず、スライダーで温度と圧力を独立に動かせるからといって、どんな組み合わせも実在する状態とは限らないということを押さえよう。例えば、1気圧(約1.013 bar)で120℃の過熱蒸気を選ぶことはできるけど、逆に10 barで80℃の「過熱蒸気」を選ぼうとすると、実はその圧力での飽和温度(約180℃)より低いから、存在しない組み合わせなんだ。ツール上では点が表示されても、それは計算上の外挿で、実際には全てが液体(圧縮水)の状態だよ。実務で使う時は、まずその圧力での飽和温度を確認して、過熱蒸気ならそれ以上、湿り蒸気ならそれ以下、と意識するのがコツだね。
次に、「乾き度x=0.9は蒸気が90%だからほぼ乾いてるでしょ?」という軽い考えは危険。伝熱設計ではこれが大きな落とし穴になる。乾き度0.9の蒸気は、重量比で確かに蒸気が9割だけど、体積比で見ると液滴が占める割合はもっと小さく、一見「乾いている」ように感じる。しかし、配管内を流れる時、これらの液滴が壁面に衝突してウォーターハンマーを引き起こしたり、タービン動翼を侵食したりする原因になる。だから、タービン入口などでは乾き度は1.0(過熱蒸気)が求められるんだ。「乾き度」は単なる割合じゃなく、機器の健全性に直結する品質指標と捉えよう。
最後に、h-s図は平衡状態を前提とした「地図」だという根本的な制約を忘れないで。実際のタービン内での膨張は、ごく短時間で起こるから、蒸気が平衡状態を保つとは限らない(非平衡過程)。特に、湿り蒸気領域での急激な膨張では、蒸気が過飽和状態(メタステブル)になることもある。その場合、図上で等エントロピー変化(断熱可逆)として描かれる線と、実際の膨張線はズレる。このツールで理想的なサイクル効率を算定した後は、必ずタービン効率や配管損失などの現実的な係数を乗せて、実際の性能を見積もる癖をつけよう。
使い方ガイド
- 圧力スライダー(sl-P)を0.001~22MPa(0.01~220 bar)の範囲で調整し、蒸気の作動圧力を設定する
- 温度スライダー(sl-T)を飽和温度以上に設定して過熱度を入力、または飽和状態で固定する
- 湿り蒸気の場合は乾き度スライダー(sl-x)を0~1.0の範囲で調整し、液相混在率を指定する
- リアルタイム計算により比エンタルピーh、比エントロピーs、比容積vが蒸気表にフィッティングした教育用の簡易相関式から自動算出される
- モリエ線図上に現在の状態点がプロット表示され、熱力学サイクルの可視化が可能
具体的な計算例
火力発電所タービン入口:圧力P=8.5MPa、温度T=400℃の過熱蒸気ではh≈3164kJ/kg、s≈6.22kJ/(kg·K)、v≈0.0365m³/kgとなります(本ツールの簡易相関式による値)。一方、圧力P=0.5MPa、乾き度x=0.85の湿り蒸気(凝縮水混在)ではh≈2427kJ/kg、s≈6.01kJ/(kg·K)、v≈0.334m³/kgと算出され、蒸気タービンの効率低下が可視化されます。
準拠規格・前提条件
準拠/参考: 飽和蒸気圧は IAPWS(IAPWS-95/IF97)の Wagner型飽和圧力相関(\(T_c=647.096\,\mathrm{K},\,P_c=22.064\,\mathrm{MPa}\)、規定の6係数)。湿り蒸気 \(h = h_f + x\,h_{fg},\; s = s_f + x\,s_{fg}\)。Clausius–Clapeyron \(dP/dT = h_{fg}/(T(v_g-v_f))\)。
モデルの前提: 飽和圧力曲線は IAPWS 準拠で高精度(100°C→1.014 bar、200°C→15.55 bar 等、誤差<0.1%)。一方 \(h_f,h_{fg},s_f,v_f,v_g\) は蒸気表への低次多項式フィット(数%級の誤差、例:\(h_f\)(200°C)=883 vs 表852 kJ/kg)。過熱蒸気は理想気体+簡易補正。
適用範囲・限界: 教育用の概略計算。臨界点近傍・極低温では誤差増大。表示の状態は外挿でも点が出るため、まずその圧力の飽和温度を確認すること。高精度設計には IAPWS-IF97 準拠ソフトを使用のこと。
実務での注意点
- 飽和圧力以上の温度入力で過熱蒸気領域へ移行し、飽和蒸気線から外れた性質が計算される(例:0.1MPaで飽和温度99.6℃、これ以上で過熱)
- 乾き度x=1.0は飽和蒸気、x=0は飽和液を示し、発電プラント復水器出口はx≈0.88~0.92に設定して蒸気品質を検証
- 比容積vが大きいほどタービン末段翼の応力増加につながるため、圧力低下時のv値監視が重要
- 冷凍サイクルやヒートポンプ設計ではR410A等の代替冷媒物性も同様計算可能、臨界点付近では精度確認が必須
