熱力学の基礎 — 熱力学法則・エントロピー・熱サイクルからCAE熱解析まで
熱力学とは何か — CAEとの繋がり
熱力学は「エネルギーの変換と移動の法則」を扱う学問です。エンジンがなぜ動くか、冷蔵庫がなぜ冷えるか、ジェットエンジンの効率の上限はどこかを支配する基礎理論です。CAEエンジニアにとっては、特に圧縮性流体解析・燃焼解析・冷凍サイクル設計の場面で不可欠な知識です。
熱力学って聞くと教科書的でとっつきにくいんですが、熱解析（伝熱）と何が違うんですか？ 同じじゃないですか？
ざっくり言うと、熱力学は「状態量の変化」を扱い、伝熱学は「その変化がどのような過程で起きるか」を扱うんだ。例えばエンジンで考えると、燃料が燃えて何Jの熱量が仕事に変換できるか上限を決めるのが熱力学。その熱がピストンやシリンダーウォールを通じてどう伝わるかの分布を計算するのが伝熱学。CAEでは両方必要で、特にガスタービンや内燃機関のシミュレーションでは切り離せない。
なるほど。で、「系・境界・外界」ってよく出てきますよね。これって何ですか？
解析対象として注目する領域を系（System）、その周りを外界（Surroundings）、系と外界を分ける仮想的な面が境界（Boundary）だ。エンジンのシリンダー内のガスを系とするなら、ピストンとシリンダーウォールが境界で、大気や冷却水が外界。境界を通じて熱と仕事のやり取りが起きる。系の種類は3つあって、外界と何もやり取りしない「孤立系」、エネルギーだけやり取りする「閉じた系（Closed System）」、エネルギーも物質もやり取りする「開いた系（Open System）」がある。CAEで扱うエンジン・タービンはほぼ開いた系だ。
1.1 熱力学 vs 伝熱学の違い
| 観点 | 熱力学 | 伝熱学 |
|---|---|---|
| 扱うもの | 状態量（温度・圧力・エンタルピー） | 過程（伝導・対流・輻射） |
| 時間 | 始状態→終状態（過程の詳細は問わない） | 時間依存の温度分布 |
| 代表的な問い | このエンジンの理論効率は何%か？ | ヒートシンクの表面温度は何℃か？ |
| CAEでの出番 | 圧縮性CFD・燃焼・冷凍サイクル設計 | 熱応力・電子冷却・鋳造プロセス |
1.2 示強性変数 vs 示量性変数
熱力学には2種類の状態量があります：
- 示強性変数（Intensive）：系の量に依存しない量。温度 $T$、圧力 $p$、密度 $\rho$、比エントロピー $s$。2つの系を合わせても値は変わらない（水1Lと水2Lを混ぜても温度は変わらない）。
- 示量性変数（Extensive）：系の量に比例する量。体積 $V$、内部エネルギー $U$、エントロピー $S$。2つの系を合わせると2倍になる。
CAEでは比（単位質量あたり）で扱うことが多く、比内部エネルギー $u = U/m$、比エンタルピー $h = H/m$ などを使います。
熱力学第0・第1法則
2.1 第0法則と温度の定義
第0法則は「A が B と熱平衡にあり、B が C と熱平衡にあるなら、A と C も熱平衡にある」というシンプルな法則です。一見自明ですが、この法則があるからこそ「温度計」という概念が成立します。温度は「熱平衡状態を識別するための状態量」として定義されます。
2.2 第1法則：エネルギー保存
閉じた系（固定質量）に対する第1法則は：
ここで $U$ は内部エネルギー、$Q$ は系に流入する熱量、$W$ は系が外部にする仕事です。符号の取り方は教科書によって異なるので注意が必要です（$W = \int p\,dV$ で外部への仕事を正とするのが工学的標準）。
エンタルピーってよく出てきますけど、内部エネルギーとどう違うんですか？
エンタルピー $H = U + pV$ は「流れのある系」で便利な量だ。ガスがパイプに入ってくるとき、そのガスは内部エネルギー $U$ を持ち込むだけじゃなく、後ろのガスに押されて $pV$ 分の仕事もしている。だから流体として持ち込む「総エネルギー」は $H$ で表すのが自然なんだ。ガスタービンのノズル設計なら「入口エンタルピーと出口エンタルピーの差＝流れの仕事」という使い方をする。冷凍システムの設計でも $h$ の差で冷凍能力や圧縮仕事を計算する。
2.3 開放系（定常流れ）のエネルギーバランス
流体機械（ポンプ、タービン、コンプレッサー、ノズル）では開放系の第1法則を使います：
ここで $\dot{W}_s$ は軸仕事（タービンで取り出す仕事など）、$h$ は比エンタルピー、$V$ は流速、$z$ は高さです。多くの実務問題では運動エネルギー項と位置エネルギー項は省略できます。
2.4 比熱と $\gamma = c_p/c_v$
定圧比熱 $c_p$ と定積比熱 $c_v$ の比 $\gamma$（比熱比）は圧縮性流体解析で極めて重要です：
- 理想気体では $c_p - c_v = R_s$（比気体定数）
- 空気では $\gamma \approx 1.4$（常温）
- 断熱過程の指数：$pV^\gamma = \text{const}$
- 音速：$c_0 = \sqrt{\gamma R_s T}$
例えば航空エンジンの計算では、高温での $\gamma$ の温度依存性（高温ほど $\gamma$ が小さくなる）を無視すると効率計算に数%の誤差が生じます。実用設計では多項式近似や JANAF テーブルを使って温度依存性を取り込みます。
熱力学第2法則とエントロピー
エントロピーって「乱雑さ」って説明されますが、なんか腑に落ちないんですよね。工学的にはどう使うんですか？
工学的にはエントロピーを「不可逆性の指標」として使うのが一番わかりやすい。摩擦で熱が出る、温度差がある場所で熱が伝わる、気体が膨張して混ざる、こういう現象は全部不可逆だ。そしてこれらが起きるたびに宇宙全体のエントロピーは増える。コンプレッサーの設計なら「等エントロピー効率」が何%か、つまり実際の仕事が理想的な断熱圧縮（エントロピー変化ゼロ）に対して何%無駄にしているかを評価する。エントロピーが増えた分だけ「使えるエネルギー」が失われたことになる。
3.1 エントロピーの定義
これは可逆過程（理想的な過程）で系に流入する微小熱量 $\delta Q_\text{rev}$ を温度 $T$ で割ったものです。不可逆過程では常にエントロピーが生成されます。
3.2 Clausius不等式
等号は可逆サイクル（カルノーサイクルなど）のとき成立します。実際の機械では常に $<0$ です。この式が第2法則の数学的表現のひとつです。
3.3 エントロピー生成（不可逆性の源）
エントロピー生成 $\dot{S}_\text{gen}$ を引き起こす主な要因：
| 不可逆性の源 | 実例 | 対策 |
|---|---|---|
| 摩擦（粘性散逸） | パイプ内の圧力損失 | 表面粗さ低減・流路最適化 |
| 有限温度差での熱移動 | 熱交換器の温度クロス | 伝熱面積増加・向流配置 |
| 非平衡混合 | バルブ絞り・急拡大管 | 段階的圧力回復 |
| 化学反応 | 燃焼・腐食 | 反応経路の制御 |
3.4 コンプレッサーの等エントロピー効率
実務でよく使うコンプレッサーの断熱（等エントロピー）効率：
$h_{2s}$ は等エントロピー過程の出口エンタルピー、$h_2$ は実際の出口エンタルピーです。市販の産業用コンプレッサーでは $\eta_c = 0.75 \sim 0.90$ 程度。残りはすべて摩擦等でエントロピー生成として失われます。
理想気体と状態方程式
4.1 理想気体の状態方程式
ここで $n$ はモル数、$R = 8.314\,\text{J/(mol·K)}$ は普遍気体定数、$R_s = R/M$ は比気体定数（$M$ は分子量）です。空気（平均分子量 $M \approx 28.97$）では $R_s \approx 287\,\text{J/(kg·K)}$ です。
実際のCFDでも理想気体を使うんですか？ 高圧・高温になると誤差が出そうで…
航空エンジンや自動車の排気系くらいの圧力と温度なら理想気体で十分な精度が出る。問題になるのは超臨界状態（液化天然ガス輸送や超臨界CO₂サイクル）や高密度気体（超高圧水素タンク）のとき。そういうケースでは van der Waals 方程式や Peng-Robinson 方程式などの実在気体状態方程式を使う。Ansys Fluent や OpenFOAM には実在気体オプションがある。大半の工業的な気体流れシミュレーションはまず理想気体から始めて、ズレが大きければ実在気体に切り替えるというのが実務の流れだ。
4.2 実在気体（van der Waals方程式）
$a$ は分子間引力（圧力補正）、$b$ は分子の有限体積（体積補正）を表します。
4.3 各種準静的過程
| 過程 | 条件 | $W = \int p\,dV$ | $\Delta U$, $Q$ |
|---|---|---|---|
| 等温（Isothermal） | $T = \text{const}$ | $nRT\ln(V_2/V_1)$ | $\Delta U = 0$、$Q = W$ |
| 等圧（Isobaric） | $p = \text{const}$ | $p(V_2 - V_1)$ | $Q = \Delta H = mc_p\Delta T$ |
| 等積（Isochoric） | $V = \text{const}$ | $0$ | $Q = \Delta U = mc_v\Delta T$ |
| 断熱（Adiabatic） | $Q = 0$ | $-\Delta U$ | $pV^\gamma = \text{const}$ |
| ポリトロープ | $pV^n = \text{const}$ | $\frac{p_1V_1 - p_2V_2}{n-1}$ | $n$ は実測で決定 |
ポリトロープ過程は実際のエンジンサイクルをモデル化するときに便利で、$n = 1$ が等温、$n = \gamma$ が断熱、$1 < n < \gamma$ が熱損失のある現実的な圧縮過程に対応します。
熱力学サイクル
カルノーサイクルって「理想的なサイクル」って言われますが、なんでこれが最高効率なんですか？ 直感的に理解したいです。
第2法則から直接導けるんだが、直感的に言うとこうだ。熱機関は高温熱源から熱を受け取って、一部を仕事に変え、残りを低温熱源に捨てる。エントロピーの変化がゼロに保たれる（可逆過程）が最も効率よく仕事に変換できる。カルノーサイクルは等温過程と断熱過程だけで構成されていて、これが唯一エントロピー生成なく動けるサイクルなんだ。「高温側の温度を上げるか低温側の温度を下げるか」しか効率を上げる方法がないというのも第2法則の重要な示唆だ。ジェットエンジンが高温燃焼（1700℃超）にこだわる理由はここにある。
5.1 カルノーサイクル
$T_H$ は高温熱源（K）、$T_L$ は低温熱源（K）。例えば $T_H = 1200\,\text{K}$、$T_L = 300\,\text{K}$ なら $\eta_\text{Carnot} = 75\%$。実際のエンジンはその半分程度の効率しか達成できません。
5.2 オットーサイクル（ガソリンエンジン）
等積燃焼（急速燃焼を等積として近似）を仮定した理想サイクル：
$r_c = V_1/V_2$ は圧縮比（典型的に $r_c = 8 \sim 12$）。$r_c = 10$、$\gamma = 1.4$ で $\eta \approx 60\%$。実際の市販エンジンは摩擦・熱損失・不完全燃焼のため $25 \sim 35\%$ 程度。
5.3 ディーゼルサイクル
$r_\alpha = V_3/V_2$ は遮断比（Cutoff ratio）です。ディーゼルの圧縮比（$r_c = 14 \sim 22$）はガソリンより高く、理論効率で上回ります。大型トラックや船舶エンジンの熱効率が高い理由です。
5.4 ブレイトンサイクル（ガスタービン）
航空エンジン・発電用ガスタービンの基礎サイクル。等圧燃焼と断熱圧縮・膨張で構成されます：
$r_p = p_2/p_1$ は圧力比、$T_1$ は圧縮前温度、$T_2$ は断熱圧縮後温度です。
航空エンジンの設計で「圧力比を上げると効率が上がる」って聞きますが、じゃあ無限に上げれば？ って思うんですけど。
鋭い。圧力比を上げると理論効率は上がる。でも現実には圧縮機出口温度も上がってしまって、タービン入口温度（TIT）を一定に保とうとすると燃焼での熱入力が減って出力が落ちる。また、圧縮機・タービンの等エントロピー効率低下も大きくなる。実用的な最適圧力比があって、最新の民間航空エンジン（GEnx, Trent XWB）では圧力比が 40〜50 に達している。さらに最近は中間冷却（Intercooling）や回生（Recuperation）を組み合わせた改良ブレイトンサイクルで熱効率 50% 超を狙う研究が盛んだ。
5.5 $p$-$V$ 線図と $T$-$s$ 線図
$p$-$V$ 線図のサイクル囲む面積 $= $ 正味仕事 $W_\text{net}$。$T$-$s$ 線図のサイクルが囲む面積 $= $ 正味仕事 $W_\text{net}$（等しい）。$T$-$s$ 線図は特に不可逆性の大小を直感的に把握するのに便利で、実際の損失（面積の変形）を理想サイクルと比較できます。
エクセルギーと有効エネルギー
6.1 エクセルギーとは
エクセルギー（有効エネルギー）は「環境状態（基準状態 $T_0, p_0$）を出発点として、系から取り出せる最大仕事」を表します：
流れのある系（定常流れ）ではフローエクセルギー：
6.2 エクセルギー損失と不可逆性
Gouy-Stodola の定理：エクセルギー損失 $= T_0 \cdot \dot{S}_\text{gen}$
つまりエントロピー生成がある＝有効仕事が失われたということです。エクセルギー解析を使うと、システムのどこで「質の高いエネルギー」が無駄になっているかを定量的に特定できます。
エクセルギー解析って実務でどこに使うんですか？
一番多いのは熱交換器とコンプレッサー・タービンのシステム評価だ。例えば廃熱回収システムで「どの熱交換器がどれだけエネルギーを無駄にしているか」を第1法則（エネルギー収支）だけで評価すると、低温でもエネルギー量は同じように見える。でもエクセルギーで見ると、低温の廃熱は有効仕事への変換効率が低い。だから「高温廃熱を先に使う」「低温廃熱は低温の用途に回す」という設計の優先順位がエクセルギー解析で明確になる。最近は水素プラントや複合サイクル発電所（GTCCプラント）の最適化設計でエクセルギー解析が必須になっている。
相変化と熱力学
7.1 Clausius-Clapeyron方程式
蒸気圧曲線（気液相境界線）の傾きを与える方程式：
$L$ は潜熱、$v_g, v_l$ はそれぞれ気体・液体の比体積です。水を例にとると、$L_\text{water} \approx 2257\,\text{kJ/kg}$（100°C）。高山では大気圧が低いため沸点も下がります（富士山頂では約87°C）。この方程式は冷媒選定・沸騰冷却設計・相変化材料（PCM）のCAEモデリングで使います。
7.2 臨界点と超臨界流体
水の臨界点は $T_c = 647.1\,\text{K}$、$p_c = 22.1\,\text{MPa}$。超臨界状態（臨界点以上）では液体と気体の区別がなくなり、密度は液体並みで粘性は気体並みの特殊な流体状態になります。超臨界CO₂サイクルは次世代原子力・太陽熱発電で注目されています。
7.3 湿り空気線図（Psychrometric Chart）
湿り空気の状態量（乾球温度・湿球温度・相対湿度・比エンタルピー・絶対湿度）の関係を示す線図。空調システムや除湿プロセスのCAE解析では湿り空気の状態変化をこの線図上で追います。Fluent の Mixture model を使った湿り空気解析では、この熱力学的関係が気液間の質量交換モデルに組み込まれています。
CAE熱流体解析との接続
8.1 圧縮性CFDでの完全エネルギー方程式
圧縮性流れの数値解析では、連続式・運動量方程式に加えてエネルギー方程式を連立させます：
$E = e + |\mathbf{u}|^2/2$ は単位体積あたりの全エネルギー（内部エネルギー + 運動エネルギー）、$\kappa$ は熱伝導率、$\boldsymbol{\tau}$ は粘性応力テンソル、$\dot{q}$ は体積発熱源（燃焼など）。
8.2 熱力学的性質の温度依存性
Sutherland の粘性式（気体）：
空気では $\mu_0 = 1.716 \times 10^{-5}\,\text{Pa·s}$、$T_0 = 273.15\,\text{K}$、$S = 110.4\,\text{K}$。Ansys Fluent・OpenFOAM ではデフォルト設定として使用可能です。
8.3 燃焼解析との連成
燃焼 CFD では化学反応による発熱がエネルギー方程式のソース項 $\dot{q}$ に入ります：
$h_k^0$ は化学種 $k$ の標準生成エンタルピー（熱力学データ）、$\dot{\omega}_k$ は反応速度、$M_k$ は分子量です。JANAF データベースから各化学種の $c_p(T)$、$h^0$、$s^0$ を取得して使います。
Fluent で燃焼解析をやろうとしたら「Species Transport」とか「Eddy Dissipation」とかモデルが多すぎて何を選べばいいかわからなくて…
燃焼モデルの選択は確かに複雑だ。まず「何を見たいか」で決める。バーナー出口の温度分布を大まかに見たいなら Eddy Dissipation Model（EDM）が手軽。でも NOx 生成量や火炎構造の詳細が必要なら Flamelet/Progress Variable（FPV）や反応速度ベースのモデルが必要になる。熱力学的には全モデルに共通して JANAF テーブルまたは NASA 多項式で各化学種の $h_k^0(T)$ と $c_{p,k}(T)$ を供給している。まず EDM でエネルギー収支が合うか確認して、そこから精度を上げていく、というのが実務的なアプローチだ。
関連インタラクティブツール
理論を手を動かして確認しよう
- 熱拡散シミュレーター — 1次元熱拡散方程式のアニメーションシミュレーション（FTCS陽解法）
- フィン効率・温度分布ツール — 矩形フィンの効率と温度プロファイルをリアルタイム計算
- 集中容量法ツール — Bi数確認 → 冷却/加熱曲線のアニメーション
- 多層壁面熱抵抗ツール — 各層の厚み・熱伝導率から温度分布・U値を計算
なった
詳しく
報告