thermodynamics-workflow について
ポンプシステム向けの熱力学解析ワークフローです。エネルギー・エントロピー収支、熱伝達、効率、サイクル計算を扱い、ポンプと熱力学サイクルを解析する技術者を対象とします。
主な機能
- 01コントロールボリュームと状態点の定義
- 02流体の熱力学的性質の評価
- 03エネルギー・エントロピー収支の計算
- 04効率、損失、エクセルギーの分析
- 05関連するサイクル計算
活用シーン
- 01ポンプ軸動力と流体の温度上昇の計算
- 02ポンプ効率と不可逆損失の評価
- 03ランキンサイクルの給水ポンプ解析
インストール
コマンドをコピーしてターミナルで実行します。
npx skills add https:/インストール後、Claude Code で/thermodynamics-workflowを入力して読み込まれたか確認します。
同梱ファイルequations-reference.mdexamples.md
詳細
熱力学解析ワークフロー
概要
このスキルでは、エネルギー収支、エントロピー解析、熱伝達、サイクル計算など、ポンプシステムの熱力学解析を体系的に行う方法を紹介します。この手法は、ポンプシステムや流体機械に特化して適用される古典熱力学の原理に基づいています。
熱力学解析ワークフロー
1. システムの定義と境界
コントロールボリュームを定義する：
熱力学解析の最初のステップは、システムの境界を明確に定義することです。
ポンプシステムの場合：
- コントロールボリューム： 吸込フランジから吐出フランジまでのポンプケーシング
- 境界： 流体やエネルギーが通過する物理的な面
- インターフェース：
- 入口：吸込フランジ（ステーション 1）
- 出口：吐出フランジ（ステーション 2）
- シャフト：機械エネルギー入力
- 壁面：熱伝達面（該当する場合）
システムの分類：
- 開放系（コントロールボリューム）： ほとんどのポンプ解析
- 定常流： 時間とともに特性が一定
- 過渡流： 起動、停止、サージ状態
記録する仮定：
- 定常運転
- 入口／出口で物性値が一様
- 運動エネルギー変化は無視できる（該当する場合）
- 断熱壁（または熱伝達を指定）
- 化学反応なし
- 非圧縮性流れ（液体の場合）または圧縮性流れ（気体の場合）
座標系：
- 位置エネルギーの基準面を定義（z = 0）
- 流れの方向を指定（正の流れ）
- 計測断面を特定
2. 状態点の特定
熱力学的状態の定義：
純物質の場合、互いに独立した2つの示強性状態量によって状態が定義されます。
状態点1（入口／吸込）：
特定して測定します：
Properties at inlet:
- Pressure: P₁ (Pa, bar, psi)
- Temperature: T₁ (K, °C)
- Velocity: c₁ (m/s)
- Elevation: z₁ (m)
- Phase: Liquid, vapor, two-phase
状態点2（出口／吐出）：
特定して測定します：
Properties at outlet:
- Pressure: P₂ (Pa, bar, psi)
- Temperature: T₂ (K, °C)
- Velocity: c₂ (m/s)
- Elevation: z₂ (m)
- Phase: Liquid, vapor, two-phase
追加の状態点：
複雑なシステムの場合：
- 段間条件（多段ポンプ）
- インペラー出口（ボリュート損失の前）
- 中間熱交換器
- 再循環ループ
物性線図：
状態点を次の線図にプロットします：
- P-h線図： 圧力と比エンタルピーの関係
- T-s線図： 温度と比エントロピーの関係
- P-v線図： 圧力と比体積の関係
- モリエル線図： h-s線図
3. 物性値の評価
流体物性値の参照元：
純物質の場合：
- 水／蒸気： IAPWS-IF97定式化（CoolProp経由）
- 冷媒： NIST REFPROPデータベース
- 炭化水素： Peng-RobinsonまたはSoave-Redlich-Kwong EOS
- 極低温流体： NIST相関式
非圧縮性液体の場合：
Specific volume: v ≈ constant
Specific heat: cp ≈ constant
Internal energy: du = cp·dT
Enthalpy: dh = cp·dT + v·dP
理想気体の場合：
Equation of state: P·v = R·T
Specific heats: cp, cv = f(T)
Internal energy: du = cv·dT
Enthalpy: dh = cp·dT
Entropy: ds = cp·dT/T - R·dP/P
物性値の評価方法：
方法1：CoolPropを使用（推奨）
from CoolProp.CoolProp import PropsSI
# Water at 25°C, 2 bar
h = PropsSI('H', 'T', 298.15, 'P', 200000, 'Water') # J/kg
s = PropsSI('S', 'T', 298.15, 'P', 200000, 'Water') # J/kg-K
rho = PropsSI('D', 'T', 298.15, 'P', 200000, 'Water') # kg/m³
方法2：物性表
- 蒸気表（過熱、飽和、圧縮液）
- 表の値の間を補間
- 線形または二次補間
方法3：経験相関式
- 密度：ρ(T, P)相関式
- 粘度：Andrade式またはVogel式
- 比熱：多項式近似
評価する主要物性値：
各状態点で評価します：
- 比エンタルピー： h（J/kgまたはkJ/kg）
- 比エントロピー： s（J/kg-KまたはkJ/kg-K）
- 比体積： v（m³/kg）または密度ρ（kg/m³）
- 比内部エネルギー： u（J/kgまたはkJ/kg）
- 乾き度： x（二相混合物の場合）
- 粘度： μ（Pa·s）— 損失計算用
- 熱伝導率： k（W/m-K）— 熱伝達用
4. 第一定律の解析（エネルギー収支）
開放系の一般エネルギー方程式：
定常流の制御体積では：
Q̇ - Ẇ = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)]
ここで：
- Q̇ = 熱移動率（W）
- Ẇ = 仕事移動率（W）
- ṁ = 質量流量（kg/s）
- h = 比エンタルピー（J/kg）
- c = 速度（m/s）
- g = 9.81 m/s²
- z = 高さ（m）
ポンプシステムの場合：
慣例上、流体には仕事が加えられる（Ẇ < 0）ため：
Ẇshaft = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)] + Q̇loss
簡略化した形式：
非圧縮性液体（密度一定）：
Ẇshaft/ṁ = (P₂ - P₁)/ρ + (c₂² - c₁²)/2 + g(z₂ - z₁) + q_loss
これを揚程で表すと：
H = (P₂ - P₁)/(ρg) + (c₂² - c₁²)/(2g) + (z₂ - z₁) + h_loss
断熱ポンプ（熱損失なし）：
Ẇshaft = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)]
運動エネルギーと位置エネルギーの変化が無視できる場合：
Ẇshaft = ṁ(h₂ - h₁) + Q̇loss
エネルギー収支の構成要素：
-
エンタルピー変化：
- 圧力上昇：Δh_pressure = ∫v·dP
- 温度上昇：Δh_thermal = cp·ΔT
-
運動エネルギーの変化：
- ポンプでは通常小さい（全体の 5% 未満）
- 高速流の用途では重要です
-
位置エネルギーの変化：
- 立形ポンプでは顕著です
- 水平設置では無視できます
-
熱移動：
- 環境への熱損失：Q̇loss
- 温度測定から求められる場合があります
- 通常、軸動力の 1～5% です
ポンプで送られる流体の温度上昇：
効率の低下による温度上昇：
ΔT = (1 - η)·Ẇshaft / (ṁ·cp)
ここで、η はポンプ効率です。
非圧縮性液体の場合：
ΔT = (1 - η)·g·H / cp
例： 水ポンプ、H = 100 m、η = 0.80、cp = 4186 J/kg-K
ΔT = (1 - 0.80)·9.81·100 / 4186 = 0.47°C
5. 第二法則の解析（エントロピーと効率）
エントロピー収支：
定常流の制御体積では：
ṁ(s₂ - s₁) = Q̇/T_boundary + Ṡgen
ここで：
- s = 比エントロピー（J/kg-K）
- Ṡgen = エントロピー生成率（W/K）
- T_boundary = 境界温度（K）
エントロピー生成：
プロセスの不可逆性を表します：
Ṡgen = ṁ(s₂ - s₁) - Q̇/T_boundary
断熱ポンプの場合：
Ṡgen = ṁ(s₂ - s₁)
等エントロピー効率：
実際のプロセスと理想的な等エントロピープロセスを比較します：
η_isentropic = (h₂s - h₁) / (h₂ - h₁)
ここで：
- h₂s = P₂ における、s₂s = s₁ の比エンタルピー（等エントロピー）
- h₂ = 出口における実際の比エンタルピー
非圧縮性液体の場合：
等エントロピー仕事は次のとおりです：
w_isentropic = ∫v·dP = v·(P₂ - P₁) = (P₂ - P₁)/ρ
エクセルギー解析：
エクセルギー（有効エネルギー）は、最大有効仕事を表します：
ψ = (h - h₀) - T₀(s - s₀) + c²/2 + gz
添字 0 は死状態（周囲条件）を示します。
エクセルギー破壊：
Ẋdestroyed = T₀·Ṡgen
第二法則効率：
η_II = Ẋuseful / Ẋinput = 1 - Ẋdestroyed / Ẋinput
これは、ポンプが利用可能なエネルギーをどれだけ効果的に使うかを表します。
エントロピー生成による損失の内訳：
さまざまな損失機構がエントロピー生成に寄与します：
- 摩擦損失： インペラーおよびボリュート内の流体摩擦
- 衝撃損失： 流れの剥離と入射
- 漏れ損失： すきまを通る再循環
- 円板摩擦： インペラー表面における粘性せん断
- 熱移動： 流体との不可逆な熱移動
6. サイクル計算
ポンプに関連する熱力学サイクル：
ポンプは通常、単一の構成要素ですが、より大きな熱力学サイクルの中で動作します。
ランキンサイクル（蒸気発電所）
構成要素：
- ボイラー（加熱）
- タービン（動力出力）
- 凝縮器（放熱）
- 給水ポンプ（圧力上昇）
ランキンサイクルにおけるポンプの解析：
状態点：
- 状態 3：凝縮器出口（飽和液または過冷却液）
- 状態 4：ボイラー入口（高圧の圧縮液）
ポンプ仕事（等エントロピー）：
w_pump,s = v₃·(P₄ - P₃) = h₄s - h₃
実際のポンプ仕事：
w_pump = w_pump,s / η_pump
サイクル効率への影響：
η_cycle = (w_turbine - w_pump) / q_in
ランキンサイクルでは、ポンプ仕事は通常、タービン仕事の 1～3% です。
ブレイトンサイクル（ガスタービン）
液体ポンプには直接適用されませんが、次の場合に関連します：
- ガスタービンの燃料ポンプ
- 液体ロケット推進剤ポンプ
- 高圧流体システム
ランキンサイクルとの主な違い：
- 作動流体は全過程で気体のままです
- ポンプ輸送ではなく圧縮を行います
- 温度比がより高くなります
冷凍サイクル
関連する用途：
- ヒートポンプの冷媒ポンプ
- チラーシステム
- 極低温ポンプ
構成要素：
- 蒸発器
- コンプレッサー（液体システムの場合はポンプ）
- 凝縮器
- 膨張弁
成績係数（COP）：
COP = Q_useful / W_input
ポンプ仕事は W_input に寄与します。
ポンプへの応用
ポンプで送液される流体の温度上昇
物理的なメカニズム：
ポンプの効率損失により、機械エネルギーが熱エネルギーに変換され、温度が上昇します。
計算方法：
エネルギーバランスから：
ΔT = (h₂ - h₁) / cp - v·(P₂ - P₁) / cp
液体の場合の簡易式：
ΔT = (1 - η)·v·(P₂ - P₁) / cp
または揚程で表すと：
ΔT = (1 - η)·g·H / cp
温度上昇に影響する要因：
- ポンプ効率： η が低い → ΔT が大きい
- 圧力上昇： ΔP が大きい → ΔT が大きい
- 比熱： cp が低い → ΔT が大きい
- 運転点： 設計点から外れる → η が低下 → ΔT が大きい
実用上の影響：
- 低温流体： 周囲に結露が生じる場合があります
- 高温流体： 沸点に近づき、キャビテーションが発生する場合があります
- 粘性流体： 大幅な発熱が生じます（η が低く、摩擦が大きい）
- 揮発性流体： 気化のリスクがあります
- ポリマー： 温度上昇により劣化する場合があります
計算例：
水ポンプ：
- H = 200 m、η = 0.75、cp = 4186 J/kg-K
- ΔT = 0.25·9.81·200/4186 = 1.17°C
作動油：
- ΔP = 200 bar、η = 0.70、ρ = 850 kg/m³、cp = 2000 J/kg-K
- ΔT = 0.30·200e5/(850·2000) = 3.53°C
高粘度流体：
- H = 50 m、η = 0.50（粘度の影響で低い）、cp = 2500 J/kg-K
- ΔT = 0.50·9.81·50/2500 = 0.98°C
効率と損失
総合効率：
η = η_h · η_vol · η_mech
各部の効率：
-
水力効率（η_h）：
- 流体摩擦損失と衝撃損失を考慮します
- 標準的な値：0.85-0.95
-
容積効率（η_vol）：
- 内部漏れを考慮します
- 標準的な値：0.96-0.99
-
機械効率（η_mech）：
- ベアリングとシールの摩擦を考慮します
- 標準的な値：0.95-0.98
熱力学的な損失分析：
各損失メカニズムによってエントロピーが増加します：
摩擦損失：
Ṡ_friction = ṁ·g·h_friction / T_avg
漏れ損失：
Ṡ_leakage = ṁ_leak·(s_discharge - s_suction)
熱伝達損失：
Ṡ_heat = Q̇_loss / T_boundary - Q̇_loss / T_fluid
総エントロピー生成量：
Ṡ_total = Ṡ_friction + Ṡ_leakage + Ṡ_heat + Ṡ_mechanical
エクセルギー効率：
η_exergy = 1 - T₀·Ṡ_total / Ẇshaft
これにより、熱力学的性能をより包括的に把握できます。
熱伝達に関する考慮事項
ポンプにおける熱伝達の形態：
-
熱伝導：
- ポンプケーシングの壁を通じた熱伝導
- シャフトからベアリングへの熱伝導
- インペラー材質を通じた熱伝導
-
対流：
- ポンプで送液される流体による強制対流
- 周囲の空気への自然対流
- 内部再循環ループ
-
放射：
- 高温でない限り、通常は無視できます
- 高温オイルポンプ（>150°C）では無視できない場合があります
熱損失の推定：
ケーシング壁から：
Q̇_loss = h·A·(T_fluid - T_ambient)
ここで：
- h = 総括熱伝達係数（W/m²-K）
- A = 表面積（m²）
- 標準的な h = 5-20 W/m²-K（自然対流、断熱あり）
エネルギーバランスから：
Q̇_loss = Ẇshaft - ṁ·(h₂ - h₁) - ṁ·(c₂² - c₁²)/2 - ṁ·g·(z₂ - z₁)
性能への影響：
- 出口温度を低下させます
- 流体の特性（粘度、密度）に影響します
- 温度勾配や変形を引き起こす場合があります
- 極低温ポンプや高温ポンプでは重要です
熱管理：
冷却要件：
- 外部冷却ジャケット
- メカニカルシール用のフラッシュプラン
- 再循環ループ内の熱交換器
断熱：
- 熱損失を防ぎます（高温流体）
- 結露を防ぎます（低温流体）
- 安全性を高めます（高温表面）
多段ポンプの熱力学
段階化の理由：
大きな揚程要件を複数の羽根車に分割します。
H_total = H_stage1 + H_stage2 + ... + H_stage_n
熱力学的解析：
段間：
State 1 → Impeller 1 → State 2 → Diffuser 1 → State 3 →
Impeller 2 → State 4 → ... → Final discharge
各段の温度上昇：
ΔT_stage = (1 - η_stage)·g·H_stage / cp
全温度上昇：
ΔT_total = Σ ΔT_stage = (1 - η_overall)·g·H_total / cp
段間冷却：
高揚程または高温の用途では、次のようにします。
Q̇_cooling = ṁ·cp·(T_beforecooling - T_aftercooling)
利点：
- 最終吐出温度を下げます
- 流体特性を改善します（粘度を下げます）
- 気化を防ぎます
- ポンプの寿命を延ばします
圧力複合効果：
段を経て圧力が上昇すると：
- 密度が変化する場合があります（圧縮性流体）
- 漏れ流量が増加します（ΔP の増大）
- シール荷重が増加します
- 構造要件が変化します
エネルギー分配：
各段の寄与は均等ではありません：
- 第1段：NPSH 要件が最も高い
- 中間段：効率が最大
- 最終段：圧力荷重が最も高い
ポンプに関係する一般的なサイクル
ランキンサイクル（詳細）
蒸気発電プラントの構成：
Boiler → Turbine → Condenser → Feedwater Pump → Boiler
↓ ↓ ↓ ↓
High P Power Low P High P
High T Output Low T Moderate T
状態点の解析：
状態 1： タービン排気が復水器に流入
- P₁ = 0.05-0.1 bar（真空条件）
- T₁ = P₁ における飽和温度
- x₁ = 0.85-0.95（湿り蒸気）
状態 2： 復水器出口／ポンプ入口
- P₂ = P₁（圧力損失は最小）
- T₂ = T_sat(P₂) - サブクーリング（通常 2-5°C）
- x₂ = 0（飽和液体またはサブクール液体）
状態 3： ポンプ出口／ボイラー入口
- P₃ = 100-300 bar（高圧）
- T₃ ≈ T₂ + ΔT_pump（温度上昇はわずか）
- 圧縮液
状態 4： ボイラー出口／タービン入口
- P₄ = P₃ - ΔP_boiler
- T₄ = 500-600°C（過熱蒸気）
- 過熱蒸気
ポンプ仕事の計算：
等エントロピー仕事：
w_pump,s = h₃s - h₂ = v₂·(P₃ - P₂)
40°C の水の場合（v ≈ 0.001008 m³/kg）：
w_pump,s = 0.001008·(250e5 - 0.08e5) ≈ 2.52 kJ/kg
実際の仕事：
w_pump = w_pump,s / η_pump = 2.52 / 0.85 ≈ 2.96 kJ/kg
温度上昇：
ΔT = (h₃ - h₂) / cp ≈ 0.44 / 4.18 ≈ 0.11°C (negligible)
サイクル効率：
η_Rankine = (w_turbine - w_pump) / q_in
一般的な値：
- 単純サイクル：30-40%
- 再熱サイクル：38-45%
- 再生サイクル：40-50%
ポンプ仕事は小さい（タービン仕事の約 1-2%）ものの、運転に不可欠です。
冷凍サイクル
液体ポンプを用いた蒸気圧縮サイクル：
一部の産業用冷凍システムでは、圧縮機の代わりに、または圧縮機に加えて液体ポンプを使用します。
状態点：
- 蒸発器出口： 低圧、低温、飽和蒸気
- 圧縮機出口： 高圧、高温、過熱蒸気
- 凝縮器出口： 高圧、中程度の温度、飽和液体
- ポンプ出口： 非常に高圧、中程度の温度、圧縮液
- 膨張弁出口： 低圧、低温、二相混合物
ポンプの用途：
- 液過供給システム： 液体冷媒を循環させます
- CO₂ 超臨界システム： 液体 CO₂ をポンプで送ります
- カスケードシステム： 圧力レベル間でポンプ輸送します
COP の計算：
COP = Q_evaporator / (W_compressor + W_pump)
ポンプ仕事は通常、圧縮機仕事に比べて小さいです。
閉ループ熱伝達サイクル
用途：
- 太陽熱システム
- 地熱ヒートポンプ
- 産業廃熱回収
- 地域冷暖房
熱力学的解析：
暖房モード：
Q̇_delivered = ṁ·cp·(T_supply - T_return)
W_pump = ṁ·g·H / η_pump
COP_heating = Q̇_delivered / W_pump
冷房モード：
Q̇_removed = ṁ·cp·(T_return - T_supply)
COP_cooling = Q̇_removed / W_pump
最適化：
- 流量を最小限に抑えます（ポンプ動力を低減します）
- ΔT を最大化します（所定の負荷に必要な流量を減らします）
- 配管摩擦とポンプサイズのバランスを取ります
まとめ
ポンプシステムの熱力学的解析では、次の項目を体系的に適用する必要があります。
- 制御体積の定義： 境界と仮定を明確にします
- 状態点の特定： 物性を完全に評価します
- 物性の評価： 信頼できる情報源から正確な流体物性を取得します
- 第一法則： あらゆる形態のエネルギーを考慮したエネルギーバランス
- 第二法則： エントロピー生成と効率の解析
- サイクル解析： より大きなシステムにおけるポンプの役割を理解します
主な出力：
- 必要動力
- 温度上昇
- 効率の内訳
- 損失メカニズム
- 最適化の機会
参考資料
詳細な方程式については equations-reference.md を、計算例については examples.md を参照してください。