熱力学は、エネルギーの変換と移動を扱う物理学の一分野です。蒸気機関の効率を探求する中で発展し、現代ではエンジン、冷凍機、発電所の設計から宇宙機の熱制御まで、幅広い工学分野の基盤となっています。
本記事では、熱力学を学ぶうえで最初に押さえるべき基本概念——系、状態量、熱力学的プロセス——を整理します。
本記事の内容
- 系の分類（開放系・閉鎖系・孤立系）
- 状態量と過程量の違い
- 代表的な熱力学的プロセス
- PV線図の意味とPythonによる可視化
系（System）の分類
熱力学では、注目する対象を 系（system）、系の外側を 外界（surroundings）と呼びます。系と外界を隔てるものが 境界（boundary）です。
| 系の種類 | 物質の移動 | エネルギーの移動 | 例 |
|---|---|---|---|
| 開放系（open system） | あり | あり | ジェットエンジン、ボイラー |
| 閉鎖系（closed system） | なし | あり | ピストン・シリンダー |
| 孤立系（isolated system） | なし | なし | 完全断熱容器 |
状態量と過程量
状態量（State Property）
系の状態のみで値が決まる量を 状態量 といいます。経路に依存しません。
主な状態量: * 圧力 $P$ [Pa] * 体積 $V$ [m$^3$] * 温度 $T$ [K] * 内部エネルギー $U$ [J] * エンタルピー $H = U + PV$ [J] * エントロピー $S$ [J/K]
数学的には、状態量の微小変化は 完全微分 で表されます。
$$ dU = \left(\frac{\partial U}{\partial T}\right)_V dT + \left(\frac{\partial U}{\partial V}\right)_T dV $$
過程量（Process Quantity）
経路に依存する量を 過程量 といいます。代表的なものは 仕事 $W$ と 熱 $Q$ です。
$$ \delta W = PdV \quad (\text{準静的過程}) $$
$$ \delta Q = TdS \quad (\text{可逆過程}) $$
$\delta$ は不完全微分を表し、経路依存であることを示しています。
状態量と過程量の違い
| 特徴 | 状態量 | 過程量 |
|---|---|---|
| 経路依存性 | なし | あり |
| 微分の表記 | $dU$, $dS$（完全微分） | $\delta W$, $\delta Q$（不完全微分） |
| 循環積分 | $\oint dU = 0$ | $\oint \delta W \neq 0$ |
| 例 | $P$, $V$, $T$, $U$, $S$, $H$ | $W$, $Q$ |
熱平衡と第0法則
熱力学第0法則: 物体Aと物体Bがそれぞれ物体Cと熱平衡にあるならば、AとBも互いに熱平衡にある。
$$ (A \sim C) \land (B \sim C) \Rightarrow (A \sim B) $$
これは温度という概念の存在を保証する法則であり、温度計の原理を支えています。
代表的な熱力学的プロセス
理想気体 $PV = nRT$ を仮定して、代表的なプロセスを整理します。
等温プロセス（Isothermal: $T = \text{const.}$）
$$ PV = nRT = \text{const.} \quad \Rightarrow \quad P \propto \frac{1}{V} $$
仕事:
$$ W = \int_{V_1}^{V_2} PdV = nRT\ln\frac{V_2}{V_1} $$
等圧プロセス（Isobaric: $P = \text{const.}$）
$$ \frac{V}{T} = \frac{nR}{P} = \text{const.} $$
仕事:
$$ W = P(V_2 – V_1) = nR(T_2 – T_1) $$
等積プロセス（Isochoric: $V = \text{const.}$）
$$ \frac{P}{T} = \frac{nR}{V} = \text{const.} $$
体積が変化しないため仕事はゼロ: $W = 0$
断熱プロセス（Adiabatic: $Q = 0$）
断熱条件と理想気体の内部エネルギー変化の関係から、
$$ PV^\gamma = \text{const.} $$
ここで $\gamma = C_P/C_V$ は比熱比です。
仕事:
$$ W = \frac{P_1V_1 – P_2V_2}{\gamma – 1} = \frac{nR(T_1 – T_2)}{\gamma – 1} $$
Pythonでの実装
4つの熱力学的プロセスをPV線図で可視化します。
import numpy as np
import matplotlib.pyplot as plt
# パラメータ
n = 1.0 # モル数
R = 8.314 # 気体定数 [J/(mol·K)]
T1 = 300 # 初期温度 [K]
P1 = 1e5 # 初期圧力 [Pa]
V1 = n * R * T1 / P1 # 初期体積
gamma = 1.4 # 比熱比（二原子分子）
V = np.linspace(V1 * 0.5, V1 * 3.0, 500)
# 等温プロセス
P_isothermal = n * R * T1 / V
# 等圧プロセス
P_isobaric = np.full_like(V, P1)
# 断熱プロセス
P_adiabatic = P1 * (V1 / V)**gamma
# 可視化
plt.figure(figsize=(10, 7))
plt.plot(V * 1e3, P_isothermal / 1e3, 'b-', linewidth=2, label=f'Isothermal ($T$={T1} K)')
plt.plot(V * 1e3, P_isobaric / 1e3, 'r-', linewidth=2, label=f'Isobaric ($P$={P1/1e3:.0f} kPa)')
plt.plot(V * 1e3, P_adiabatic / 1e3, 'g-', linewidth=2, label=f'Adiabatic ($\\gamma$={gamma})')
# 等積プロセス（垂直線）
plt.axvline(x=V1 * 1e3, color='purple', linewidth=2, linestyle='--',
label=f'Isochoric ($V$={V1*1e3:.1f} L)')
# 初期状態のマーク
plt.plot(V1 * 1e3, P1 / 1e3, 'ko', markersize=10, zorder=5, label='Initial state')
plt.xlabel('Volume [L]')
plt.ylabel('Pressure [kPa]')
plt.title('Thermodynamic Processes on PV Diagram')
plt.legend(fontsize=11)
plt.grid(True)
plt.ylim(0, 250)
plt.tight_layout()
plt.show()
PV線図上で、各プロセスの経路が異なることが確認できます。等温過程と断熱過程はともに双曲線的ですが、断熱過程のほうが急峻に圧力が変化します。
まとめ
本記事では、熱力学の基本概念について解説しました。
- 系 は開放系・閉鎖系・孤立系に分類され、物質とエネルギーの境界を通じた移動で区別する
- 状態量（$P$, $V$, $T$, $U$, $S$）は経路に依存せず、過程量（$W$, $Q$）は経路に依存する
- 第0法則 は温度の概念を保証する
- 等温・等圧・等積・断熱の4つの基本プロセスはPV線図上でそれぞれ異なる経路をたどる
次のステップとして、以下の記事も参考にしてください。