熱力学は、エネルギーの変換や移動に関する法則を体系的に扱う学問です。火力発電所でのタービン駆動、ロケットエンジンの推進、エアコンや冷蔵庫の冷凍サイクルなど、工学のあらゆる場面で熱力学の原理が活用されています。
本記事では、熱力学を学び始めるうえで最初に押さえておくべき基本概念を整理します。系の分類、状態量と過程量の違い、温度・圧力・内部エネルギー・エンタルピーの定義、そして理想気体の状態方程式まで、丁寧に解説していきます。
本記事の内容
- 熱力学とは何か
- 系（システム）の分類
- 状態量と過程量の違い
- 温度・圧力の定義
- 内部エネルギーとエンタルピー
- 状態変化の種類
- 理想気体の状態方程式
- Pythonで各状態変化をpV図に可視化
熱力学とは
熱力学（Thermodynamics）は、エネルギーの変換と移動の法則を記述する物理学・工学の一分野です。
大雑把に言うと、「熱」と「仕事」というエネルギーの形態がどのように変換され、どのような制約のもとで変化が起きるかを扱う学問です。歴史的には、蒸気機関の効率向上を目的として18〜19世紀に発展しました。
熱力学が基盤となっている工学分野は非常に多く、代表的なものとして以下が挙げられます。
| 分野 | 熱力学の応用例 |
|---|---|
| 推進工学 | ロケットエンジン、ジェットエンジンの燃焼とノズル膨張 |
| 発電工学 | 火力・原子力発電のランキンサイクル |
| 冷凍・空調 | エアコン、冷蔵庫の冷凍サイクル |
| 化学工学 | 化学反応における平衡と自発性 |
| 材料工学 | 相変態、材料の熱処理 |
これらの工学分野を理解するための第一歩が、本記事で扱う熱力学の基本概念です。
系（システム）の分類
熱力学では、まず注目する対象とその外側を区別するために、系（system） と 外界（surroundings） という概念を導入します。系と外界を合わせたものを 宇宙（universe） と呼びます。
系は、その境界が物質やエネルギーをどのようにやり取りするかによって、3種類に分類されます。
開放系（Open system）
物質とエネルギーの両方が境界を通じて出入りできる系です。ジェットエンジンやタービンのように、作動流体が流入・流出する系がこれに該当します。
閉鎖系（Closed system）
エネルギー（熱や仕事）は出入りできるが、物質の出入りがない 系です。ピストン-シリンダ装置で密封された気体が典型例です。気体がピストンを押して仕事をしたり、外部から加熱されたりしますが、気体自体は系の外に出ません。
孤立系（Isolated system）
物質もエネルギーも境界を通じて出入りしない系です。完全に断熱された容器がこれに近い理想的なモデルです。実際にはエネルギーの出入りをゼロにすることは難しいため、理論的な極限として使われます。
これらの分類は、熱力学の法則を適用する際に「何がどう保存されるか」を考えるための出発点になります。
状態量と過程量の違い
熱力学で扱う物理量は、状態量（state variable） と 過程量（path variable / process variable） に分けられます。この区別は非常に重要です。
状態量
系の現在の状態だけで値が一意に決まる量を 状態量 と言います。過去にどのような経路をたどったかに依存しません。温度 $T$、圧力 $p$、体積 $V$、内部エネルギー $U$、エンタルピー $H$ などがこれに該当します。
イメージとしては、山の「標高」に近いです。山頂の標高は、どのルートで登っても同じ値になりますよね。状態量も同様に、どのような過程（プロセス）を経ても、同じ状態にたどり着けば同じ値を取ります。
数学的には、状態量の微小変化は完全微分（exact differential） で表されます。例えば内部エネルギー $U$ の微小変化は $dU$ と書きます。
過程量
系がある状態から別の状態へ変化する際の経路（過程）に依存する量を 過程量 と言います。代表的なものは 熱 $Q$ と 仕事 $W$ です。
山登りで言えば、「歩いた距離」に相当します。同じ山頂に着くとしても、急な直登ルートと緩やかな迂回ルートでは歩いた距離が異なります。同様に、同じ始状態と終状態の間でも、過程が違えば授受される熱や仕事の量は異なります。
過程量の微小変化は不完全微分として $\delta Q$、$\delta W$ のように記号 $\delta$ を用いて書かれることが多いです。
温度の定義と第0法則
温度は日常的に馴染み深い概念ですが、熱力学において厳密に定義するには熱力学第0法則（Zeroth Law of Thermodynamics） が必要です。
熱力学第0法則は次のように述べられます。
物体Aと物体Bがそれぞれ物体Cと熱平衡にあるならば、物体Aと物体Bも互いに熱平衡にある。
ここで 熱平衡 とは、2つの物体を接触させたときに正味の熱の流れがない状態を指します。
この法則は一見当たり前に思えますが、これがあるからこそ「温度計」という共通の尺度（物体C）を使って、別々の物体の温度を比較できるわけです。つまり、第0法則は温度という状態量が定義できる根拠を与えています。
温度のスケールとしては、以下のものが代表的です。
| スケール | 記号 | 基準 |
|---|---|---|
| セルシウス温度 | $T\,[\mathrm{^\circ C}]$ | 水の凝固点を0、沸点を100 |
| 絶対温度（ケルビン） | $T\,[\mathrm{K}]$ | 分子運動が停止する理論的な下限を0 |
両者の関係は次の通りです。
$$ T\,[\mathrm{K}] = T\,[\mathrm{^\circ C}] + 273.15 $$
熱力学の計算では、通常 絶対温度（ケルビン） を使います。
圧力の定義
圧力 $p$ は、単位面積あたりに作用する力として定義されます。
$$ p = \frac{F}{A} $$
ここで、$F$ は面に垂直に作用する力 $[\mathrm{N}]$、$A$ は面積 $[\mathrm{m^2}]$ です。SI単位はパスカル $[\mathrm{Pa}] = [\mathrm{N/m^2}]$ です。
工学では、絶対圧力（absolute pressure）とゲージ圧力（gauge pressure）の区別が重要です。
$$ p_{\mathrm{abs}} = p_{\mathrm{gauge}} + p_{\mathrm{atm}} $$
ここで $p_{\mathrm{atm}} \approx 101325\,\mathrm{Pa}$ は標準大気圧です。熱力学の式では絶対圧力を使います。
内部エネルギー $U$
内部エネルギー $U$ は、系を構成する分子の運動エネルギーとポテンシャルエネルギーの総和です。これは状態量であり、系の巨視的な運動（系全体が動いている等）は含みません。
直感的に言えば、ある温度・圧力・体積の状態にある気体が「内部に蓄えているエネルギーの総量」です。
理想気体の場合、分子間の相互作用を無視するため、内部エネルギーは温度のみの関数になります。
$$ U = U(T) $$
定積モル比熱 $c_v$（体積一定のもとで温度を1 K上げるのに必要な単位物質量あたりの熱量）を用いると、内部エネルギーの変化は次のように書けます。
$$ \Delta U = n c_v \Delta T $$
ここで $n$ はモル数、$\Delta T$ は温度変化です。
エンタルピー $H = U + pV$
エンタルピー $H$ は、次のように定義される状態量です。
$$ H = U + pV $$
ここで $U$ は内部エネルギー、$p$ は圧力、$V$ は体積です。
エンタルピーがなぜ便利かというと、定圧過程（圧力一定の変化）における系が吸収する熱量を直接表すからです。これを確認してみましょう。
定圧過程における熱力学第1法則は次のようになります。
$$ \delta Q = dU + p\,dV $$
圧力 $p$ が一定であるとき、$d(pV) = p\,dV + V\,dp$ のうち $dp = 0$ なので $d(pV) = p\,dV$ が成り立ちます。これを代入すると、
$$ \begin{align} \delta Q &= dU + p\,dV \\ &= dU + d(pV) \\ &= d(U + pV) \\ &= dH \end{align} $$
となります。つまり、定圧過程では系に加えられた熱がそのままエンタルピーの変化に等しくなるということです。
$$ Q_p = \Delta H $$
定圧モル比熱 $c_p$ を用いると、理想気体のエンタルピー変化は次のようになります。
$$ \Delta H = n c_p \Delta T $$
また、理想気体では定圧モル比熱と定積モル比熱の間にマイヤーの関係式が成り立ちます。
$$ c_p – c_v = R $$
ここで $R = 8.314\,\mathrm{J/(mol \cdot K)}$ は気体定数です。
状態変化の種類
気体の状態変化には、代表的なものとして以下の4種類があります。
定容変化（Isochoric process）
体積 $V$ が一定の状態変化です。気体が膨張も圧縮もしないため、系は仕事をしません（$W = 0$）。加えた熱はすべて内部エネルギーの変化になります。
$$ Q_v = \Delta U = n c_v \Delta T $$
定圧変化（Isobaric process）
圧力 $p$ が一定の状態変化です。先ほど導出したように、加えた熱はエンタルピーの変化に等しくなります。
$$ Q_p = \Delta H = n c_p \Delta T $$
等温変化（Isothermal process）
温度 $T$ が一定の状態変化です。理想気体では内部エネルギーが温度のみの関数なので $\Delta U = 0$ となり、加えた熱はすべて仕事に変換されます。
等温過程で理想気体の状態方程式 $pV = nRT$ から、$p = nRT / V$（$T$ は定数）が成り立つので、pV図上では反比例の曲線（双曲線） になります。
$$ pV = \text{const.} $$
断熱変化（Adiabatic process）
系と外界の間で熱のやり取りがない変化です（$Q = 0$）。断熱過程では、比熱比 $\gamma = c_p / c_v$ を用いて次の関係が成り立ちます。
$$ pV^{\gamma} = \text{const.} $$
この導出を確認しましょう。断熱過程では $\delta Q = 0$ なので、熱力学第1法則より、
$$ 0 = dU + p\,dV $$
理想気体に対して $dU = n c_v\,dT$ ですから、
$$ n c_v\,dT + p\,dV = 0 $$
また、理想気体の状態方程式 $pV = nRT$ の両辺を微分すると、
$$ p\,dV + V\,dp = nR\,dT $$
ここから $dT$ を求めると、
$$ dT = \frac{p\,dV + V\,dp}{nR} $$
これを先ほどの式に代入します。
$$ n c_v \cdot \frac{p\,dV + V\,dp}{nR} + p\,dV = 0 $$
$$ \frac{c_v}{R}(p\,dV + V\,dp) + p\,dV = 0 $$
$$ \frac{c_v}{R} \cdot p\,dV + \frac{c_v}{R} \cdot V\,dp + p\,dV = 0 $$
$p\,dV$ の項をまとめます。
$$ p\,dV \left(\frac{c_v}{R} + 1\right) + \frac{c_v}{R} \cdot V\,dp = 0 $$
マイヤーの関係式 $c_p – c_v = R$ より $c_v / R + 1 = (c_v + R) / R = c_p / R$ なので、
$$ \frac{c_p}{R} \cdot p\,dV + \frac{c_v}{R} \cdot V\,dp = 0 $$
両辺を $pV \cdot c_v / R$ で割ると、
$$ \frac{c_p}{c_v} \cdot \frac{dV}{V} + \frac{dp}{p} = 0 $$
$\gamma = c_p / c_v$ とおけば、
$$ \gamma \frac{dV}{V} + \frac{dp}{p} = 0 $$
両辺を積分すると、
$$ \gamma \ln V + \ln p = \text{const.} $$
$$ \ln(pV^{\gamma}) = \text{const.} $$
$$ pV^{\gamma} = \text{const.} $$
これが断熱過程における圧力と体積の関係式です。$\gamma > 1$ であるため、断熱変化の曲線はpV図上で等温変化の曲線より急勾配になります。
理想気体の状態方程式
ここまでの議論の多くで使ってきた理想気体の状態方程式を改めて整理しておきましょう。
$$ pV = nRT $$
ここで、
| 記号 | 意味 | 単位 |
|---|---|---|
| $p$ | 圧力 | $\mathrm{Pa}$ |
| $V$ | 体積 | $\mathrm{m^3}$ |
| $n$ | モル数 | $\mathrm{mol}$ |
| $R$ | 気体定数 | $8.314\,\mathrm{J/(mol \cdot K)}$ |
| $T$ | 絶対温度 | $\mathrm{K}$ |
この式は、分子間相互作用を無視し、分子自体の体積も無視できるという2つの仮定のもとで成り立つ近似モデルです。低圧・高温の条件では実在気体もこの式によく従います。
PythonでpV図を可視化する
理想気体の4つの状態変化（定容・定圧・等温・断熱）をpV図上に描画してみましょう。初期状態を $p_1 = 200\,\mathrm{kPa}$、$V_1 = 0.5\,\mathrm{m^3}$ とし、1 mol の単原子理想気体（$\gamma = 5/3$）を想定します。
import numpy as np
import matplotlib.pyplot as plt
# --- パラメータ設定 ---
p1 = 200e3 # 初期圧力 [Pa]
V1 = 0.5 # 初期体積 [m^3]
n = 1.0 # モル数 [mol]
R = 8.314 # 気体定数 [J/(mol·K)]
gamma = 5.0 / 3 # 比熱比（単原子理想気体）
# 初期温度を状態方程式から計算
T1 = p1 * V1 / (n * R)
# 体積の範囲
V = np.linspace(0.2, 1.5, 500)
# --- 各状態変化における圧力 ---
# 1. 定容変化（V = V1 で一定、圧力が変化）
# pV図上では垂直線になるため、圧力の範囲を指定して描画
p_isochoric = np.linspace(50e3, 400e3, 500)
V_isochoric = np.full_like(p_isochoric, V1)
# 2. 定圧変化（p = p1 で一定）
p_isobaric = np.full_like(V, p1)
# 3. 等温変化（T = T1 で一定 → pV = nRT1 → p = nRT1 / V）
p_isothermal = n * R * T1 / V
# 4. 断熱変化（pV^γ = p1 * V1^γ → p = p1 * (V1/V)^γ）
C_adiabat = p1 * V1**gamma
p_adiabatic = C_adiabat / V**gamma
# --- 可視化 ---
fig, ax = plt.subplots(figsize=(10, 7))
ax.plot(V_isochoric * 1e3, p_isochoric / 1e3,
color='tab:blue', linewidth=2.5, label='Isochoric (constant $V$)')
ax.plot(V * 1e3, p_isobaric / 1e3,
color='tab:orange', linewidth=2.5, label='Isobaric (constant $p$)')
ax.plot(V * 1e3, p_isothermal / 1e3,
color='tab:green', linewidth=2.5, label='Isothermal (constant $T$)')
ax.plot(V * 1e3, p_adiabatic / 1e3,
color='tab:red', linewidth=2.5, label='Adiabatic ($pV^\\gamma$ = const.)')
# 初期状態を点で表示
ax.plot(V1 * 1e3, p1 / 1e3, 'ko', markersize=8, zorder=5, label='Initial state')
ax.set_xlabel('Volume $V$ [L]', fontsize=14)
ax.set_ylabel('Pressure $p$ [kPa]', fontsize=14)
ax.set_title('$p$-$V$ Diagram: Four Thermodynamic Processes', fontsize=15)
ax.set_xlim(100, 1500)
ax.set_ylim(0, 450)
ax.legend(fontsize=12)
ax.grid(True, alpha=0.3)
plt.tight_layout()
plt.show()
このコードを実行すると、4つの状態変化がpV図上でどのような曲線を描くかが一目で分かります。特に、等温変化（緑）と断熱変化（赤）はどちらも体積が増えると圧力が減少しますが、断熱変化の方が急勾配になっていることが確認できます。これは、断熱変化では膨張に伴い温度が低下する分、さらに圧力が下がるためです。
まとめ
本記事では、熱力学の基本概念について解説しました。
- 系の分類: 開放系（物質・エネルギー出入り可）、閉鎖系（エネルギーのみ出入り可）、孤立系（出入りなし）の3種類がある
- 状態量と過程量: 状態量は経路に依存せず状態だけで決まる量（$T$, $p$, $V$, $U$, $H$ 等）、過程量は経路に依存する量（$Q$, $W$）
- 温度と圧力: 温度は第0法則によって定義が正当化される。圧力は単位面積あたりの力
- 内部エネルギー $U$: 系が内部に蓄えるエネルギー。理想気体では温度のみの関数
- エンタルピー $H = U + pV$: 定圧過程での熱量の授受を直接表す便利な状態量
- 4つの状態変化: 定容・定圧・等温・断熱の各プロセスにおける基本関係式を整理した
- 理想気体の状態方程式: $pV = nRT$
次の記事では、エネルギー保存の法則を熱力学的に定式化した 熱力学第1法則 について、詳しく解説していきます。