ポリトロープ変化シミュレーター 戻る
熱力学

ポリトロープ変化シミュレーター

気体が P·Vⁿ=一定 の法則に沿って状態を変える「ポリトロープ変化」を可視化するツールです。ポリトロープ指数 n をスライダーで動かすと、定圧・等温・断熱・定容の4つの理想過程が1つのモデルで連続的に切り替わり、最終圧力・温度比・仕事・熱量・内部エネルギー変化がリアルタイムで分かります。

パラメータ設定
ポリトロープ指数 n
0=定圧、1=等温、γ=断熱、∞=定容を連続的に切り替える
初期圧力 P₁
kPa
状態1の圧力。自然吸気で約100kPa
初期比体積 v₁
m³/kg
状態1の単位質量あたりの体積
体積比 V₂/V₁
1未満で圧縮、1で変化なし。最終体積を決める
比熱比 γ
c_p/c_v。空気で約1.40、単原子ガスで約1.67
計算結果
最終圧力 P₂ (kPa)
温度比 T₂/T₁
仕事 W (kJ/kg)
熱量 Q (kJ/kg)
内部エネルギー変化 ΔU (kJ/kg)
変化の種類
P-V線図 — 過程アニメーション

太い曲線が現在の指数 n のポリトロープ変化(状態1→2)。薄い曲線は同じ始点を通る等温(n=1)と断熱(n=γ)の参考カーブです。網掛けの面積が気体のする仕事を表します。

P-V曲線(圧力 vs 体積)
仕事 vs ポリトロープ指数 n
理論・主要公式

$$P\,V^{n}=\text{const},\qquad \frac{T_2}{T_1}=\left(\frac{V_1}{V_2}\right)^{n-1}$$

ポリトロープ変化の基本則と温度比。n はポリトロープ指数、V は比体積、T は絶対温度。圧縮(V₂<V₁)では温度が上がる。

$$W=\frac{P_1V_1-P_2V_2}{n-1}\quad(n\ne1)$$

気体が外部にする仕事 W(n≠1)。n=1 の等温変化では 0/0 を避けて W=P₁V₁·ln(V₁/V₂) を用いる。

$$\Delta U=c_v(T_2-T_1),\qquad Q=\Delta U+W$$

内部エネルギー変化と熱力学第一法則による熱量。c_v=R/(γ−1)、R=287 J/(kg·K)。n=0,1,γ,∞ はそれぞれ定圧・等温・断熱・定容変化を再現する。

ポリトロープ変化とは

🙋
熱力学の授業で「定圧変化」「等温変化」「断熱変化」「定容変化」って4つ別々に習ったんですけど、「ポリトロープ変化」っていうのはまた別の5つ目の過程なんですか?
🎓
いい質問だね。実はその逆で、ポリトロープ変化はその4つを全部まとめて1つにした「親玉」みたいな存在なんだ。式は P·Vⁿ=一定 というシンプルなもので、n という1つの数(ポリトロープ指数)を変えるだけで、4つの過程を連続的に行き来できる。左のスライダーで n を 0 から 3 まで動かしてごらん。n=0 で定圧、n=1 で等温、n=γ で断熱、n をうんと大きくすると定容に近づく。4つの「点」が1本の「線」でつながるイメージだよ。
🙋
えっ、たった1つの数でそんなに切り替わるんですか。じゃあ n=1 のときだけ仕事の式が違うのはなぜですか?P-V線図の下の網掛けが面積だってことは分かるんですけど。
🎓
するどい。一般の式 W=(P₁V₁−P₂V₂)/(n−1) を見てみると、n=1 を入れると分母が n−1=0 になってしまう。だからそのまま代入すると 0÷0 で計算できないんだ。でも数学的にちゃんと極限をとると、n=1 のときは W=P₁V₁·ln(V₁/V₂) という対数の式に化ける。これが等温変化の仕事だ。このツールでは n が 1 にごく近いとき自動でこの対数の式に切り替えているから、スライダーを n=1 ぴったりにしても変な数字(NaN)が出ない。
🙋
なるほど。でも実際のエンジンや圧縮機って、ちょうど n=1 とか n=γ になることはあるんですか?
🎓
それがまさにポリトロープ変化が役に立つ理由なんだ。実機はほぼ絶対に「ちょうど」にはならない。例えば空気圧縮機を考えてみよう。完全な等温(n=1)にするには圧縮しながら無限に冷やし続ける必要があるけど、そんなことは不可能だ。逆に完全な断熱(n=γ)にするには熱が壁から1ジュールも逃げてはいけないけど、現実には必ず多少は逃げる。だから実際の圧縮機は n が 1 と γ のあいだ、たとえば 1.25 とか 1.3 くらいの「ポリトロープ」になる。理想のどれにも当てはまらない現実を、ちょうど1と γ のあいだの n で表せるんだよ。
🙋
じゃあ実機の n はどうやって知るんですか?設計するときに勝手に決まるものなんですか?
🎓
測って決めるんだ。実機で圧力と体積を測定して、その P-V のデータを P·Vⁿ=一定 の式に当てはめる。両辺の対数をとると logP と logV が直線になるから、その傾きが n になる。n さえ分かれば、あとはこのツールがやっているように仕事も熱量も温度上昇も全部計算で出せる。圧縮機やエンジンの性能を「ポリトロープ効率」で評価するのは、現場ではごく普通のやり方だよ。1つの指数で現実をまるごと記述できる、これがポリトロープ変化の強みだ。
🙋
最後にひとつ。圧縮したのに「仕事 W」がマイナスで出るのはバグじゃないですよね?
🎓
バグじゃないよ、ちゃんと物理的に正しい。このツールの W は「気体が外部にする仕事」という符号の取り方をしている。圧縮のときは外から気体を押し込むわけだから、気体は仕事を「される」側だ。だから気体がする仕事 W はマイナスになる。デフォルト設定だと W≈−123 kJ/kg、つまり1kgの空気をこの割合で圧縮するのに約123kJを外から与えた、という意味だ。膨張させれば符号は逆になってプラスになる。符号がどっち向きの仕事を指しているかは、熱力学でいつも最初に確認すべき大事なポイントだよ。

よくある質問

ポリトロープ変化とは、気体の圧力 P と体積 V が P·Vⁿ = 一定 という関係を保ちながら状態を変える過程のことです。ここで n はポリトロープ指数と呼ばれる1つの数で、この n を変えるだけで熱力学で別々に学ぶ4つの理想過程を連続的に再現できます。n=0 は定圧変化、n=1 は等温変化、n=γ(比熱比)は可逆断熱変化、n→∞ は定容変化に対応します。実機の圧縮や膨張は等温でも断熱でもない中間の挙動をするため、ポリトロープ変化が現実を表すのに最も使われます。
n が1でない一般の場合、気体が外部にする仕事は W = (P₁V₁ − P₂V₂)/(n − 1) で計算します(圧力はPa、比体積は m³/kg)。n がちょうど1の等温変化のときはこの式が0/0になるため、特別な式 W = P₁V₁·ln(V₁/V₂) を使います。圧縮(V₂ < V₁)では気体は外部から仕事をされるため、気体がする仕事 W は負の値になります。本ツールは n がほぼ1のとき自動で等温の式に切り替えるため、NaN を出さずに連続した結果が得られます。
実測した圧力と体積のデータを P·Vⁿ = 一定 の式に当てはめ、両辺の対数をとって直線の傾きから n を求めます。空気圧縮機なら n は1(完全な等温、無限の冷却が必要)と γ≈1.4(完全な断熱、熱損失ゼロ)の間に入り、実際には n=1.2〜1.35 程度になることが多いです。冷却がよく効くほど n は1に近づき、断熱に近いほど γ に近づきます。この n が分かれば、仕事・熱量・温度上昇がすべて計算で求まります。
n=0 では P·V⁰=P=一定 となり圧力が変わらない定圧変化(等圧変化)です。n=1 では P·V=一定 で温度が変わらない等温変化(ボイルの法則)です。n=γ では可逆断熱変化(等エントロピー変化)で、境界を熱が出入りしません。n→∞ では P^(1/n)·V=一定 から体積が変わらない定容変化(等容変化)になります。このように1つの指数 n を連続的に動かすだけで、4つの基本過程をすべて1つの枠組みで扱えるのがポリトロープ変化の最大の利点です。

実世界での応用

空気圧縮機・エアコンプレッサーの性能評価:工場のエアコンプレッサーは、空気を吸い込んで高圧にする典型的なポリトロープ機械です。理想を言えば等温圧縮(n=1)が最も少ない仕事で済みますが、それには無限の冷却が必要で不可能です。実際の圧縮機は壁面冷却やインタークーラーの効き具合に応じて n=1.2〜1.35 程度で運転されます。測定した P-V 線図から n を求め、本ツールのように仕事を計算することで、どれだけ理想からかけ離れているか、冷却を強化すれば何キロワット節約できるかを定量的に評価できます。

内燃機関の圧縮行程・膨張行程:ガソリンエンジンやディーゼルエンジンのシリンダー内で起こる圧縮と膨張も、純粋な断熱変化ではありません。燃焼前の圧縮行程では混合気がシリンダー壁へ熱を逃がし、燃焼後の膨張行程でも熱損失があります。そのためエンジン解析では各行程をポリトロープ変化として扱い、圧縮行程で n≈1.3、膨張行程で n≈1.3〜1.35 といった実測値を使ってサイクル計算を行います。1つの指数で熱損失込みの現実を表現できるのが利点です。

ガスタービンと圧縮機・タービン段の設計:ガスタービンの軸流圧縮機やタービンの各段でも、流れと羽根のあいだの摩擦や熱の出入りのため、過程は理想の断熱からずれます。これをポリトロープ効率(多段に分けても一定とみなせる効率)で表すのが業界の標準的なやり方です。本ツールで n を γ から少しずらすと熱量 Q がゼロから外れていく様子が見られ、これが「断熱からのずれ」を意味します。

熱力学教育とサイクル理解の土台:ポリトロープ変化は、定圧・等温・断熱・定容の4過程を「バラバラの暗記」ではなく「1つの指数 n の連続体」として理解するための最良の教材です。オットーサイクルやブレイトンサイクルなどの各行程は、すべてポリトロープ変化の特別な場合として位置づけられます。本ツールで n を連続的に動かし、P-V 線図のカーブの曲がり方や仕事・熱量の符号がどう変わるかを見ることで、サイクル全体の理解が一段深まります。

よくある誤解と注意点

まず多い誤解が、「ポリトロープ指数 n は気体の物性で決まる定数だ」というものです。比熱比 γ は気体の種類(空気、二酸化炭素、ヘリウムなど)で決まる物性値ですが、n は物性値ではありません。n は「その過程で熱がどれだけ出入りするか」という過程の性質を表す数で、同じ空気でも冷却をよく効かせれば n は1に近づき、断熱に近づければ γ に近づきます。n は気体ではなく過程に属する量だ、という区別を最初にはっきりさせてください。本ツールで γ と n を別々のスライダーにしているのはそのためです。

次に、「仕事 W と熱量 Q の符号の取り方を混同する」こと。熱力学第一法則には Q = ΔU + W(W は気体がする仕事)という工学系の流儀と、ΔU = Q − W のように書く流儀、さらに W を「気体がされる仕事」と定義する教科書もあり、符号の約束が混在しています。本ツールは「W = 気体が外部にする仕事」で統一しており、圧縮では W が負、Q = ΔU + W で熱量を求めています。デフォルト設定では ΔU が正(温度上昇で内部エネルギー増加)、W が負(圧縮)、その和の Q が負(正味では放熱)になります。他の資料の数値と比べるときは、必ずどちらの符号約束かを確認してください。

最後に、「ポリトロープ変化なら何でも可逆だと思い込む」こと。P·Vⁿ=一定 という式は、過程の各点で気体が平衡状態にある(準静的である)ことを暗黙に仮定しています。現実の急速な圧縮や、強い乱流・衝撃をともなう流れは厳密には準静的ではなく、摩擦による不可逆性をともないます。実機データに当てはめて得た n は「見かけのポリトロープ指数」であり、その背後には不可逆損失が隠れています。n=γ ちょうどが「可逆断熱(等エントロピー)」を意味するのに対し、実機で測った n が γ に近くても、それは熱損失が小さいだけで完全に可逆とは限らない、という点に注意してください。

使い方ガイド

  1. ポリトロープ指数n(0.0~2.0)を設定します。n=1.0で等温変化、n=1.4で断熱変化(空気)、n=1.67で単原子気体の断熱変化となります
  2. 初期圧力P₁(50~500 kPa)と初期体積V₁(0.1~10 m³/kg)を入力します
  3. 体積比(V₂/V₁ = 0.2~5.0)を指定すると、P·Vⁿ=一定の法則から最終状態が自動計算されます
  4. シミュレーターがP₂、温度比T₂/T₁、仕事W、熱量Q、内部エネルギー変化ΔUをリアルタイムで出力します

具体的な計算例

空気(理想気体、R=287 J/kg·K)を圧縮シリンダー内で処理する場合、初期条件:P₁=100 kPa、T₁=300 K、V₁=1.0 m³/kgとし、n=1.4(断熱圧縮)、体積比V₂/V₁=0.5とします。P·Vⁿ=一定より、P₂=100×(1/0.5)^1.4≈189.3 kPaとなります。温度比はT₂/T₁=(P₂/P₁)^((n-1)/n)=1.189が得られ、T₂=356.7 Kです。このとき圧縮仕事W=-92.8 kJ/kgとなり、断熱変化なのでQ=0、内部エネルギー変化ΔU=92.8 kJ/kgです。

実務での注意点

  1. ターボチャージャーやコンプレッサー設計では、実際の効率を考慮するため計算結果に0.75~0.85のポリトロープ効率を乗じます
  2. ディーゼルエンジンの燃焼室内膨張では、n=1.2~1.3程度の値を採用する場合が多く、完全な断熱変化(n=1.4)より現実的です
  3. 水蒸気やCO₂などの実在気体では、ポリトロープ指数が温度・圧力範囲で変動するため、狭い範囲に限定して使用してください
  4. 体積比が極端に大きい(>3.0)場合、計算精度低下の可能性があるため、段階的に分割計算することを推奨します