熱力学的相転移のモデル
熱力学的平衡 において 、安定性の必要条件は、圧力 が 体積またはモル体積 とともに増加しないことである。 これは数学的には と表現され 、ここで は 温度である。 [1]
この基本的な安定性要件、および他の 共役変数ペア に対する同様の要件は、一次相 転移 の解析モデルでは満たされない。最も有名なケースは 、 ファンデルワールス方程式である。 [2] [3]
p
{\displaystyle p}
ヴ
=
五
/
いいえ
{\displaystyle v=V/N}
∂
ヴ
p
|
T
<
0
{\displaystyle \partial _{v}p|_{T}}
T
{\displaystyle T}
p
=
R
T
/
(
ヴ
−
b
)
−
1つの
/
ヴ
2
{\displaystyle p=RT/(vb)-a/v^{2}}
ここで、 は次元定数です。この違反は欠陥ではなく、むしろ液体と蒸気を区別する特性において観察される不連続性の原因であり、一次相転移を定義します。
1つの
、
b
、
R
{\displaystyle a,b,R}
図 1: 曲線は、相変化を含む流体の --平面 における等温線で一定です 。曲線のさまざまなセグメントについては、本文で説明します。には と同様に、 から への ジャンプ不連続点があることに注意してください 。
T
{\displaystyle T}
p
{\displaystyle p}
ヴ
{\displaystyle v}
ヴ
{\displaystyle v}
ヴ
メートル
私
ん
{\displaystyle v_{\rm {min}}}
ヴ
メートル
1つの
x
{\displaystyle v_{\rm {max}}}
p
{\displaystyle p}
図 1 は、連続的に微分可能な黒の実線、黒の点線、灰色の破線で描かれた等温線を示しています 。図 1 の点 C の右側の曲線の減少部分は気体を表し、点 E の左側の曲線の減少部分は液体を表します。これらの 2 つの部分は、安定性基準に違反する正の傾きを持つ曲線上の極小値と極大値の間の領域によって分離されています。この数学的基準は、エプスタイン [4] が次のように説明した物理的条件を表現しています 。「曲線の点線部分 [図 1 の破線部分] には、物理的実体がないことは明らかです。実際、曲線のこの部分に対応する状態の流体が、上部がピストンで形成された熱伝導垂直シリンダーに含まれていると想像してみましょう。ピストンはシリンダー内で上下にスライドでき、ガスの圧力と正確に釣り合う負荷をかけます。ピストンから少し重量を取り除くと、平衡がなくなり、ピストンは上向きに動き始めます。しかし、ピストンが動くと、ガスの体積が増加し、それとともに圧力も増加します。ピストンにかかる結果的な力は大きくなり、ピストンは上向きのままになります。したがって、ピストンは動き続け、ガスは等温線の最大値で表される状態に達するまで膨張します。逆に、バランスの取れたピストンに少しでも負荷を加えると、ガスは等温線の最小値に対応する状態に収縮します。」
ヴ
>
b
{\displaystyle v>b}
この状況は、滑らかな表面の上で正確にバランスをとっている物体が、わずかな乱れによって平衡位置から外れ、局所的最小値に達するまで続く状況に似ています。このような状態は説明されているように動的に不安定であり、したがって観察されません。ギャップは、液体から気体への実際の相変化の前兆です。 最大の液体状態と最小の気体状態を区切る 点 E と C (ただし、 ) は スピノーダル 点と呼ばれます。その軌跡は、均質な安定状態が存在できない領域を区切るスピノーダル曲線を形成します。
ヴ
メートル
私
ん
≤
ヴ
≤
ヴ
メートル
1つの
x
{\displaystyle v_{\rm {最小}}\leq v\leq v_{\rm {最大}}}
(
p
メートル
私
ん
、
ヴ
メートル
私
ん
)
{\displaystyle (p_{\rm {min}},v_{\rm {min}})}
(
p
メートル
1つの
x
、
ヴ
メートル
1つの
x
)
{\displaystyle (p_{\rm {max}},v_{\rm {max}})}
∂
ヴ
p
|
T
=
0
{\displaystyle \partial _{v}p|_{T}=0}
実験によれば、一定量の液体が入った容器の容積が加熱され、一定の温度、ある圧力で膨張すると 、蒸気(図 1 の点 と点 でドットで表示)が核形成し、流体は均質ではなくなり、沸騰する液体と凝縮する蒸気の不均質な混合物になる。重力により、沸騰する(飽和)液体 は、同じ飽和温度と圧力で共存する 密度の低い凝縮する(飽和)蒸気から分離される 。加熱が続くと、蒸気の量が 増加し、液体 の量は 減少する。その間、圧力 と温度 は 一定のままで、体積は 増加する。この状況では、混合物のモル体積は、その成分の加重平均である。
p
s
(
T
)
{\displaystyle p_{s}(T)}
f
{\displaystyle f}
g
{\displaystyle g}
v
f
=
V
f
/
N
f
{\displaystyle v_{f}=V_{f}/N_{f}}
v
g
=
V
g
/
N
g
>
v
f
,
{\displaystyle v_{g}=V_{g}/N_{g}>v_{f},}
N
g
{\displaystyle N_{g}}
N
f
=
N
−
N
g
{\displaystyle N_{f}=N-N_{g}}
p
s
{\displaystyle p_{s}}
T
{\displaystyle T}
V
=
V
f
+
V
g
{\displaystyle V=V_{f}+V_{g}}
v
=
V
/
N
=
(
V
f
/
N
f
)
(
N
−
N
g
)
/
N
+
(
V
g
/
N
g
)
(
N
g
/
N
)
=
v
f
(
1
−
x
)
+
v
g
x
{\displaystyle v=V/N=(V_{f}/N_{f})(N-N_{g})/N+(V_{g}/N_{g})(N_{g}/N)=v_{f}(1-x)+v_{g}x}
ここで 、蒸気のモル分率 は 連続的に増加します。しかし、物質自体のモル容積は、与えられた において、液体状態では最大の安定値しか持たず、蒸気状態では最小の安定値しか持ちません 。繰り返しますが、混合物のモル容積は から まで連続的に変化しますが ( 図1の破線で示す)、基礎となる流体はこの特性だけでなく他の特性にも不連続性を持っています。この混合物の状態方程式は、 てこ則 と呼ばれます。 [5] [6] [7]
x
=
N
g
/
N
{\displaystyle x=N_{g}/N}
0
≤
x
≤
1
{\displaystyle 0\leq x\leq 1}
p
(
T
)
{\displaystyle p(T)}
v
f
{\displaystyle v_{f}}
v
g
{\displaystyle v_{g}}
図 1 の曲線の点線部分は 準安定状態 です。長年、このような状態は学術的な興味の対象でした。Callen [8] は 、「1 気圧で 0°C 以下に冷却された水。この状態でビーカーの水に軽く触れると、突然、系が劇的に結晶化します」を例に挙げました。しかし、沸騰熱伝達の研究により、準安定状態はこのプロセスの不可欠な部分として日常的に発生することが明らかになりました。この状態では、加熱面の温度が飽和温度よりも高く、多くの場合大幅に高いため、隣接する液体は過熱されている必要があります。 [9] さらに、非常に高い熱流束で動作するデバイスの出現により、準安定状態、およびそれに関連する熱力学的特性、特に過熱液体状態への関心が高まっています。 [10] さらに、これらの状態がファンデルワールス方程式、および一般的な 3 次方程式によって予測されるという事実は、相転移を説明する上でのその有効性の説得力のある証拠です。ゾンマーフェルトは次のように述べている。 [11]
ファンデルワールス理論が、少なくとも定性的には、 AA` または BB` [ここでは図 1 の BC と FE] の分岐に沿った不安定な [ここでは準安定と呼ばれる] 状態の存在を予測できることは非常に注目に値します。
等面積ルール
における不連続性 と、 物質の 内部エネルギー 、エントロピー などの他の特性は、一次相転移と呼ばれます。 [12] [13] 実験的に観測された唯一の圧力 を特定するには、 別の熱力学的条件が必要です。図 1 から、範囲 のどの圧力でもそれが起こる可能性があることがわかります 。このような条件は、1875 年 2 月 18 日にマクスウェルが英国化学会で行った講義で、巧みな熱力学的議論の中で初めて表明されました [14] (文字 BCDEF を含む図 1 は、彼が説明した曲線です)。
v
{\displaystyle v}
u
{\displaystyle u}
s
{\displaystyle s}
p
s
(
T
)
{\displaystyle p_{s}(T)}
p
m
i
n
≤
p
≤
p
m
a
x
{\displaystyle p_{\rm {min}}\leq p\leq p_{\rm {max}}}
これは、Nature 誌の James Clerk-Maxwell の論文に掲載された図の表現です。
曲線の Cから E までの部分は 、本質的に不安定な点を表しており、したがって実現することはできません。
ここで、媒体が常に均質な状態で 仮想曲線 B C D E Fに沿って Bから F へ通過し 、液体と蒸気の混合物の形で直線経路 F B に沿って戻ると仮定します。温度は最初から最後まで一定であるため、熱は仕事に変換されません。ここで、仕事に変換された熱は、面積 F D E と B C D の差で表されます。したがって、特定の温度での蒸気の最大圧力を決定する条件は、線 B F が 曲線の上下から等しい面積を切り取ることです。
マクスウェルの議論を理解する最も簡単な方法は、彼が示唆したサイクルを温度-モルエントロピー平面上で考えることです。熱力学の入門書 [15]
[16] では必ず、そのような平面上では、任意の曲線の下の面積が物質へのモルあたりの 熱伝達 であり、左から右に行くほど正味で、右から左に行くほど負であるという事実が示されています。さらに、循環プロセスでは、物質への正味の熱伝達は、サイクルの閉じた曲線で囲まれた面積です。彼が考えたサイクルは、同じ温度での2つの灰色の破線の等温線で構成されており、1つはBからF(CDとEを経由)に進み、もう1つはFからBに直接戻るため、2つの線は同一であり、単に逆方向に通過しているだけです。囲まれた面積はゼロなので、 となります 。
q
=
0
{\displaystyle q=0}
さらに、同じテキストでは、これらの曲線の下の面積を圧力-モル体積上にプロットした場合(図 1 参照)、その面積は物質によって 行われた仕事 であり、左から右に向かうにつれて正になり、右から左に向かうにつれて負になる、と説明しています。同様に、サイクル内で行われた正味の仕事は、閉じた曲線で囲まれた面積です。熱力学の第一法則により、サイクルの特殊なケースでは 、マクスウェルが想定したサイクルに対して となるため 、囲まれた面積は I+II=0 であるため(図 1 参照)、I が正で II が負であるため、遷移圧力は 2 つの面積が等しくなるようにする必要があります。
w
=
q
{\displaystyle w=q}
w
=
q
=
0
{\displaystyle w=q=0}
各プロセスで行われた仕事に関して数式で表すと、
∫
v
g
v
f
p
d
v
+
∫
v
f
v
g
p
s
d
v
=
−
∫
v
f
v
g
p
d
v
+
p
s
(
v
g
−
v
f
)
=
0
for
T
=
constant
{\displaystyle \int _{v_{g}}^{v_{f}}\,p\,dv+\int _{v_{f}}^{v_{g}}\,p_{s}\,dv=-\int _{v_{f}}^{v_{g}}\,p\,dv+p_{s}(v_{g}-v_{f})=0\quad {\mbox{for}}\quad T={\mbox{constant}}}
この式は、各状態について書かれた状態方程式と合わせて 、
f
{\displaystyle f}
g
{\displaystyle g}
p
s
=
p
(
v
f
,
T
)
p
s
=
p
(
v
g
,
T
)
{\displaystyle p_{s}=p(v_{f},T)\qquad p_{s}=p(v_{g},T)}
は 4 つの変数 に対する 3 つの方程式であるため 、そのうちの 1 つ、たとえば が与えられれば 、他の 3 つが決まります。言い換えると、 、 およびには 、 相転移が発生する唯一の値があります。
p
s
,
T
,
v
f
,
v
g
{\displaystyle p_{s},T,v_{f},v_{g}}
T
{\displaystyle T}
p
s
(
T
)
{\displaystyle p_{s}(T)}
v
f
(
T
)
{\displaystyle v_{f}(T)}
v
g
(
T
)
{\displaystyle v_{g}(T)}
ギブス基準
講義の最後に、マクスウェルはファンデルワールスの研究を「非常に独創的な論文」と称えて賛辞を述べた後、次のように締めくくった。
しかし、私は、熱力学のこの分野における、米国イェール大学のウィラード・ギブス 教授による最も重要なアメリカの貢献について言及しないわけにはいきません。 彼は、モデルによって物質のさまざまな状態の関係を表す、驚くほど単純で完全に満足のいく方法を与えてくれました。このモデルによって、私や他の人たちの努力が長らく阻んできた問題が、すぐに解決されるかもしれません。
この発言は先見の明があったことが証明された。というのは、1876年から1878年にかけて ギブスは 熱力学に関する決定的な著作 [17]を出版し、その中で 不均質物質 の熱力学的平衡 には、 機械的平衡 (各成分に同じ圧力)と 熱平衡(各成分に同じ温度)に加えて、 物質的平衡 (各成分に同じ化学ポテンシャル) も必要であることを示したからである。1つの物質と2つの相の本例では 、およびに加えて、 物質的平衡には (1つの物質の特殊なケースでは、その化学ポテンシャルはモルギブス関数であり、 )が必要である。 [ 18] この条件は、次のような単純な物理的議論によって導き出すことができる。1 モルを蒸発させるのに必要なエネルギーは 、定温での第二法則から であり 、定圧での第一法則から であり 、これら2つを等しくして整理すると であるため、結果が得られる 。物質平衡の条件は、有名なギブス の相律 、につながります。 ここで は、 物質の数、 相の数、および 状態を特定するために必要な独立した強意変数の数です。 [19] [20] ここで議論されている1つの物質と2つの相の場合、これは 、実験的に観測された数 を与えます。
p
f
=
p
g
=
p
s
{\displaystyle p_{f}=p_{g}=p_{s}}
T
f
=
T
g
=
T
{\displaystyle T_{f}=T_{g}=T}
g
f
=
g
g
{\displaystyle g_{f}=g_{g}}
μ
=
g
≡
G
/
N
{\displaystyle \mu =g\equiv G/N}
g
=
u
+
p
v
−
T
s
{\displaystyle g=u+pv-Ts}
q
v
a
p
=
T
(
s
g
−
s
f
)
{\displaystyle q_{\rm {vap}}=T(s_{g}-s_{f})}
q
v
a
p
=
h
g
−
h
f
{\displaystyle q_{\rm {vap}}=h_{g}-h_{f}}
h
=
u
+
p
v
{\displaystyle h=u+pv}
D
=
n
−
r
+
2
{\displaystyle D=n-r+2}
n
{\displaystyle n}
r
{\displaystyle r}
D
{\displaystyle D}
D
=
1
{\displaystyle D=1}
は熱力学的ポテンシャル関数 であり 、 その微分は [21]
g
(
p
,
T
)
{\displaystyle g(p,T)}
d
g
=
∂
p
g
|
T
d
p
+
∂
T
g
|
p
d
T
=
v
d
p
−
s
d
T
{\displaystyle dg=\partial _{p}g|_{T}\,dp+\partial _{T}g|_{p}\,dT=v\,dp-s\,dT}
これを一定温度で積分すると、
g
(
p
,
T
)
=
g
A
(
T
)
+
∫
p
A
p
v
(
p
¯
,
T
)
d
p
¯
{\displaystyle g(p,T)=g_{A}(T)+\int _{p_{A}}^{p}\,v({\bar {p}},T)\,d{\bar {p}}}
ここで は 積分定数であるが、この定数は等温線ごとに異なるため の関数として表される 。 [22] を評価するに は を反転して を得る必要がある 。しかし、この反転は一意ではないのが相転移現象の性質である。例えば、 について書かれたファンデルワールス方程式は 、
g
A
{\displaystyle g_{A}}
T
{\displaystyle T}
g
{\displaystyle g}
p
=
p
(
v
,
T
)
{\displaystyle p=p(v,T)}
v
=
v
(
p
,
T
)
{\displaystyle v=v(p,T)}
v
{\displaystyle v}
p
v
3
−
(
p
b
+
R
T
)
v
2
+
a
v
−
a
b
=
0
,
{\displaystyle pv^{3}-(pb+RT)v^{2}+av-ab=0,}
図 2: 図 1 に示したのと同じ等温線上のギブス関数。ここでも文字は図 1 と同じ点を表します。
1 またはこの場合は 3 つの実根を持つ 3 次方程式です。したがって、図 2 に示すように、安定状態 (黒の実線で表示)、準安定状態 (黒の点線で表示)、不安定状態 (灰色の破線で表示) の 3 つの曲線があります。
実際、この図は3 次方程式 を解いて積分することによって得られたもので はなく、まず ファンデルワールス方程式に対して解析的に簡単に行える と の 定義から得たもので 、それを をパラメータとして で パラメトリック にプロットしたものである 。安定状態のみを考慮すると、 相転移点では は不連続な 偏微分 、 、 とともに連続である。 エーレンフェスト 分類 では、一次相転移は の一次偏微分が不連続になることを指し 、二次相転移は二次偏微分が不連続になることを指す。 [23]
g
(
v
,
T
)
{\displaystyle g(v,T)}
u
(
v
,
T
)
{\displaystyle u(v,T)}
s
(
v
,
T
)
{\displaystyle s(v,T)}
p
(
v
,
T
)
{\displaystyle p(v,T)}
v
{\displaystyle v}
g
(
p
,
T
)
{\displaystyle g(p,T)}
∂
p
g
|
T
=
v
{\displaystyle \partial _{p}g|_{T}=v}
∂
T
g
|
p
=
−
s
{\displaystyle \partial _{T}g|_{p}=-s}
g
{\displaystyle g}
ギブス基準とマクスウェル基準の関係
先に与えられた飽和液体と蒸気状態の間の
積分式を評価し、この相変化プロセスに物質平衡のギブス基準を適用すると、次のように書くことができる。
g
(
p
,
T
)
{\displaystyle g(p,T)}
g
g
−
g
f
=
∫
p
s
p
m
i
n
v
l
d
p
+
∫
p
m
i
n
p
m
a
x
v
u
d
p
+
∫
p
m
a
x
p
s
v
v
d
p
=
0
{\displaystyle g_{g}-g_{f}=\int _{p_{s}}^{p_{\rm {min}}}\,v_{l}\,dp+\int _{p_{\rm {min}}}^{p_{\rm {max}}}\,v_{u}\,dp+\int _{p_{\rm {max}}}^{p_{s}}\,v_{v}\,dp=0}
ここで、積分は、それぞれ液体 、不安定、 および蒸気、 状態に対応する 3 次方程式の 3 つの実根を使用して 3 つの部分に分割されています。これらの積分は、図 1 を 紙面上で反時計回りに回転させ、次に軸を 中心に回転させて、添付のグラフに示すように曲線の左側の縦座標に が表示される ようにすると、最もよく視覚化さ れます。このビューでは、関数は明らかに多値です。これが 、と の間の動作を記述するために 3 つの実関数が必要な理由です 。次に、中央の積分を 2 つに分割します。
v
l
{\displaystyle v_{l}}
v
u
{\displaystyle v_{u}}
v
v
{\displaystyle v_{v}}
90
∘
{\displaystyle 90^{\circ }}
180
∘
{\displaystyle 180^{\circ }}
v
{\displaystyle v}
v
{\displaystyle v}
v
(
p
,
T
)
{\displaystyle v(p,T)}
p
m
i
n
{\displaystyle p_{\rm {min}}}
p
m
a
x
{\displaystyle p_{\rm {max}}}
g
g
−
g
f
=
∫
p
s
p
m
i
n
v
l
d
p
+
∫
p
m
i
n
p
s
v
u
d
p
+
∫
p
s
p
m
a
x
v
u
d
p
+
∫
p
m
a
x
p
s
v
v
d
p
=
0
{\displaystyle g_{g}-g_{f}=\int _{p_{s}}^{p_{\rm {min}}}\,v_{l}\,dp+\int _{p_{\rm {min}}}^{p_{s}}\,v_{u}\,dp+\int _{p_{s}}^{p_{\rm {max}}}\,v_{u}\,dp+\int _{p_{\rm {max}}}^{p_{s}}\,v_{v}\,dp=0}
ここで最初の 2 つの積分は面積 I で、次の 2 つは面積 II の負数です。2 つの面積を足すと 0 になるため、このギブス基準によれば大きさは等しくなります。これもマクスウェルの等面積則、マクスウェル構成であり、解析的にも示せます 。
d
(
p
v
)
=
p
d
v
+
v
d
p
{\displaystyle d(pv)=pdv+vdp}
d
g
=
v
d
p
−
s
d
T
=
−
p
d
v
+
d
(
p
v
)
−
s
d
T
.
{\displaystyle dg=v\,dp-s\,dT=-p\,dv+d(pv)-s\,dT.}
これをギブス条件で
状態 から 温度一定で積分すると、
f
{\displaystyle f}
g
{\displaystyle g}
g
g
−
g
f
=
−
∫
v
f
v
g
p
(
v
,
T
)
d
v
+
p
s
(
v
g
−
v
f
)
=
0
{\displaystyle g_{g}-g_{f}=-\int _{v_{f}}^{v_{g}}\,p(v,T)\,dv+p_{s}(v_{g}-v_{f})=0}
これはマクスウェルの結果である。この等面積則はヘルムホルツの自由エネルギーを利用して導くこともできる。 [24] いずれにせよ、マクスウェルの構成は物質の平衡のギブス条件から導かれる。しかし、 より基本的ではあるが、幾何学的に理解される等面積則よりも抽象的である。
g
f
=
g
g
{\displaystyle g_{f}=g_{g}}
共通接線の構築
図 3: 図 1 と 2 に示したのと同じ等温線のヘルムホルツ関数と、B と F の両方での接線。ここでも文字は、図 1 と 2 と同じ点を表します。
共存点を決定する別の方法は、ヘルムホルツ電位最小原理に基づいており、これは、 熱貯蔵庫 、およびと 透熱 接触しているシステムでは、すなわち平衡時にヘルムホルツ電位が最小になるというものである。 [25] と同様に 、モルヘルムホルツ関数 も電位関数であり、その微分は次のように なる。 [26]
T
=
T
R
{\displaystyle T=T_{R}}
D
F
=
0
{\displaystyle DF=0}
D
2
F
>
0
{\displaystyle D^{2}F>0}
g
(
p
,
T
)
{\displaystyle g(p,T)}
f
(
v
,
T
)
{\displaystyle f(v,T)}
d
f
=
∂
v
f
|
T
d
v
+
∂
T
f
|
v
d
T
=
−
p
d
v
−
s
d
T
,
{\displaystyle df=\partial _{v}f|_{T}\,dv+\partial _{T}f|_{v}\,dT=-p\,dv-s\,dT,}
この最小原理は安定条件を導く 。 [27] この条件は、系のどの安定状態においても関数が 厳密に 凸で あること、すなわちその近傍において曲線がその接線上かその上にあることを要求する。 [28] さらに、これらの状態では、圧力の前述の安定条件も必然的に満たされる。
∂
2
f
/
∂
v
2
|
T
=
−
∂
p
/
∂
v
|
T
>
0
{\displaystyle \partial ^{2}f/\partial v^{2}|_{T}=-\partial p/\partial v|_{T}>0}
f
{\displaystyle f}
図 1 および 2 と同じ vdW 方程式の亜臨界等温線に対するこの関数のプロットが、図 3 に示されています。この図には、B および F で関数の曲線と二重の (共通の) 接線を持つ (破線/実線) 直線が含まれています。この直線は、
定数 の で 、 と書くことができます 。最後の等式は、関係 、 [29] とから得られます 。 [18] これはすべて、直線上のすべての点、特に点 B と F が同じ値を持つことを意味し、これにより、物質の平衡 と温度と圧力の等式 に関するギブス条件が生成されます。 [30] したがって、この構成は、ギブス条件とマクスウェル構成の両方に相当します。
f
{\displaystyle f}
f
=
f
0
+
∂
v
f
T
v
=
f
0
−
p
v
{\displaystyle f=f_{0}+\partial _{v}f_{T}v=f_{0}-pv}
p
{\displaystyle p}
f
0
=
f
+
p
v
=
g
{\displaystyle f_{0}=f+pv=g}
f
=
u
−
T
s
{\displaystyle f=u-Ts}
g
=
u
−
T
s
+
p
v
{\displaystyle g=u-Ts+pv}
g
,
p
,
T
{\displaystyle g,p,T}
g
f
=
g
g
{\displaystyle g_{f}=g_{g}}
この構成は、ギブスが以前に定義した に基づいており、 [31] [32] ファンデルワールスによって最初に使用されました(彼はこれを二重接線と共通接線の両方と呼びました)。 [33] これは、の等温線が組成変数 とともに共通接面を持つことができる表面を形成するような2成分流体 混合物 を含むように簡単に拡張できるためです 。その後、混合物の相変化の問題を扱う一般的な方法になりました。 [34] [35] [36]
f
(
v
,
T
)
{\displaystyle f(v,T)}
f
(
v
,
x
,
T
)
{\displaystyle f(v,x,T)}
x
=
N
1
/
(
N
1
+
N
2
)
{\displaystyle x=N_{1}/(N_{1}+N_{2})}
ファンデルワールス方程式への応用
ファンデルワールス方程式を飽和液体、 および蒸気に適用すると 、状態は
v
f
<
v
m
i
n
{\displaystyle v_{f}<v_{\rm {min}}}
v
g
>
v
m
a
x
{\displaystyle v_{g}>v_{\rm {max}}}
p
s
=
R
T
s
v
f
−
b
−
a
v
f
2
p
s
=
R
T
s
v
g
−
b
−
a
v
g
2
{\displaystyle p_{s}={\frac {RT_{s}}{v_{f}-b}}-{\frac {a}{v_{f}^{2}}}\qquad p_{s}={\frac {RT_{s}}{v_{g}-b}}-{\frac {a}{v_{g}^{2}}}}
これら2つの方程式は4つの変数を指定しており、について 解くことができます 。その結果は次のようになります。
p
s
,
T
s
{\displaystyle p_{s},T_{s}}
v
f
,
v
g
{\displaystyle v_{f},v_{g}}
p
s
=
p
∗
v
∗
2
[
v
f
v
g
−
v
∗
(
v
f
+
v
g
)
]
v
f
2
v
g
2
T
s
=
T
∗
v
∗
(
v
f
+
v
g
)
(
v
f
−
v
∗
)
(
v
g
−
v
∗
)
v
f
2
v
g
2
{\displaystyle p_{s}=p^{*}{\frac {v^{*2}[v_{f}v_{g}-v^{*}(v_{f}+v_{g})]}{v_{f}^{2}v_{g}^{2}}}\qquad T_{s}=T^{*}{\frac {v^{*}(v_{f}+v_{g})(v_{f}-v^{*})(v_{g}-v^{*})}{v_{f}^{2}v_{g}^{2}}}}
ここで 、、、 は 定数で定義される特性圧力、モル容積、温度である( であることに注意 )。マクスウェル構造をファンデルワールス方程式に適用すると、次の式が得られる。
p
∗
=
a
/
b
2
{\displaystyle p^{*}=a/b^{2}}
v
∗
=
b
{\displaystyle v^{*}=b}
T
∗
=
a
/
(
R
b
)
{\displaystyle T^{*}=a/(Rb)}
p
∗
v
∗
/
T
∗
=
R
{\displaystyle p^{*}v^{*}/T^{*}=R}
−
T
s
ln
(
v
g
−
v
∗
v
f
−
v
∗
)
+
T
∗
v
∗
(
v
g
−
v
f
)
v
f
v
g
+
p
s
(
v
g
−
v
f
)
/
R
=
0
{\displaystyle -T_{s}\ln \left({\frac {v_{g}-v^{*}}{v_{f}-v^{*}}}\right)+T^{*}{\frac {v^{*}(v_{g}-v_{f})}{v_{f}v_{g}}}+p_{s}(v_{g}-v_{f})/R=0}
これら3つの方程式は数値的に解くことができる。これは、 またはのいずれかの値を与えて行われ 、表形式で結果が提示されている。 [37] [38] しかし、これらの方程式は解析的なパラメトリック解も許し、Lenknerによれば [39] 、 Gibbsによって得られた。Lenkner自身がこの解を得るための単純で洗練された方法を考案した。方程式から と を消去し、 が から に 変化する につれて と 0 の 間で変化する 引き伸ばされた無次元密度 で表す 。これにより、
T
s
{\displaystyle T_{s}}
p
s
{\displaystyle p_{s}}
p
s
{\displaystyle p_{s}}
T
s
{\displaystyle T_{s}}
ϱ
=
v
∗
/
(
v
−
v
∗
)
{\displaystyle \varrho =v^{*}/(v-v^{*})}
∞
{\displaystyle \infty }
v
{\displaystyle v}
v
∗
{\displaystyle v^{*}}
∞
{\displaystyle \infty }
ln
(
ϱ
f
ϱ
g
)
=
(
ϱ
f
−
ϱ
g
)
(
ϱ
g
+
ϱ
f
+
2
)
ϱ
f
+
ϱ
g
+
2
ϱ
f
ϱ
g
{\displaystyle \ln \left({\frac {\varrho _{f}}{\varrho _{g}}}\right)={\frac {(\varrho _{f}-\varrho _{g})(\varrho _{g}+\varrho _{f}+2)}{\varrho _{f}+\varrho _{g}+2\varrho _{f}\varrho _{g}}}}
この方程式は超越的であるが、方程式の左辺を と書くことで、単純な解析的パラメトリック解が得られる 。
(
s
g
−
s
f
)
/
R
{\displaystyle (s_{g}-s_{f})/R}
ln
(
ϱ
f
ϱ
g
)
=
s
g
−
s
f
R
=
δ
=
2
y
{\displaystyle \ln \left({\frac {\varrho _{f}}{\varrho _{g}}}\right)={\frac {s_{g}-s_{f}}{R}}=\delta =2y}
そして、 右辺から を 消去すると、の線形方程式 が得られ、その解は
ϱ
f
=
e
δ
ϱ
g
{\displaystyle \varrho _{f}=e^{\delta }\varrho _{g}}
ϱ
f
{\displaystyle \varrho _{f}}
ϱ
g
{\displaystyle \varrho _{g}}
ϱ
g
(
y
)
=
f
(
y
)
e
−
y
and
ϱ
f
(
y
)
=
f
(
y
)
e
y
where
f
(
y
)
=
y
cosh
y
−
sinh
y
sinh
y
cosh
y
−
y
{\displaystyle \varrho _{g}(y)=f(y)e^{-y}\quad {\mbox{and}}\quad \varrho _{f}(y)=f(y)e^{y}\quad {\mbox{where}}\quad f(y)={\frac {y\cosh y-\sinh y}{\sinh y\cosh y-y}}}
したがって、この相転移過程における他のすべての変数を指定する基本変数は である 。飽和問題に対するこの解は、すべての変数を包含するように簡単に拡張できる。
(
s
g
−
s
f
)
/
(
2
R
)
=
y
{\displaystyle (s_{g}-s_{f})/(2R)=y}
T
s
(
y
)
=
T
∗
2
f
(
y
)
[
cosh
y
+
f
(
y
)
]
g
(
y
)
2
p
s
(
y
)
=
p
∗
f
(
y
)
2
[
1
−
f
(
y
)
2
]
g
(
y
)
2
{\displaystyle T_{s}(y)=T^{*}{\frac {2f(y)[\cosh y+f(y)]}{g(y)^{2}}}\qquad p_{s}(y)=p^{*}{\frac {f(y)^{2}[1-f(y)^{2}]}{g(y)^{2}}}}
v
f
(
y
)
=
v
∗
1
+
f
(
y
)
e
y
f
(
y
)
e
y
v
g
(
y
)
=
v
∗
1
+
f
(
y
)
e
−
y
f
(
y
)
e
−
y
{\displaystyle v_{f}(y)=v^{*}{\frac {1+f(y)e^{y}}{f(y)e^{y}}}\quad \quad \quad \quad \quad \quad v_{g}(y)=v^{*}{\frac {1+f(y)e^{-y}}{f(y)e^{-y}}}}
どこ 。
g
(
y
)
=
[
1
+
ϱ
g
(
y
)
]
[
1
+
ϱ
f
(
y
)
]
=
1
+
2
f
(
y
)
cosh
y
+
f
(
y
)
2
{\displaystyle g(y)=[1+\varrho _{g}(y)][1+\varrho _{f}(y)]=1+2f(y)\cosh y+f(y)^{2}}
図 4: -平面 上の解析的な vdW スピノーダル曲線 (青い破線) と共存曲線 (黒い破線) 。
p
{\displaystyle p}
T
{\displaystyle T}
図5: 図1の等温線、 臨界等温線、 超臨界等温線、スピノーダル曲線および共存曲線。
T
r
=
7
/
8
{\displaystyle T_{r}=7/8}
T
r
=
1
{\displaystyle T_{r}=1}
飽和曲線上の他のすべての特性の不連続性の値もこの解から導かれます。 [40]
これらの関数は、ファンデルワールス流体の飽和液体と飽和蒸気状態の軌跡である共存曲線を定義します。図4では、この曲線は青でプロットされ、黒のスピノーダル曲線は次のように計算されます。
T
=
T
∗
(
v
~
−
1
)
2
/
v
~
3
p
=
p
∗
(
v
~
−
2
)
/
v
~
3
{\displaystyle T=T^{*}({\tilde {v}}-1)^{2}/{\tilde {v}}^{3}\quad \quad p=p^{*}({\tilde {v}}-2)/{\tilde {v}}^{3}}
ここで、 はパラメータである。これらのプロットを作成する際に使用される変数は、縮小された(次元のない)変数、 、 、および であり、 下付き文字の量は臨界点の値である。これらは、 、および 臨界点において [41] で定義され、 は測定可能な量である。関係 、 、 は、 解の星の量を、 図で使用される量に変換するために使用される。曲線は、前に参照した数値結果と完全に一致する。スピノーダル曲線の内側の領域では、各点に 2 つの状態、1 つは安定で 1 つは準安定であり、青い曲線の右側の過熱液体または左側の過冷却蒸気のいずれかであるが、スピノーダル曲線の外側では、各点に 1 つの安定状態がある。図 5 では、(黒の点線)スピノーダル曲線の下の領域には均質な安定状態が含まれないが、(赤の点線)共存曲線とスピノーダル曲線の間には、各点に 1 つの準安定状態があり、共存曲線の外側では、各点に 1 つの安定状態がある。 2 つの青い円と 2 つの緑の円は、それぞれの等温線上で飽和した液体と蒸気の状態を示しています。共存曲線の下には、てこ則を満たす不均一状態も至る所で観測されていますが、それらはファンデルワールス方程式の均一状態ではないため、各亜臨界等温線上の飽和点を結ぶ水平線で示されるそれらの存在は表示されません。また、この図の横軸は線形ではなく対数になっています。これは、 液体領域と不安定領域を過度に圧縮することなく、蒸気領域を広く表示するためです 。ただし、このデバイスは領域を歪めるため、図 1 の 2 つの領域 I と II はここでは等しく表示されません。
v
~
=
v
/
v
∗
{\displaystyle {\tilde {v}}=v/v^{*}}
p
r
=
p
/
p
c
{\displaystyle p_{r}=p/p_{c}}
v
r
=
v
/
v
c
{\displaystyle v_{r}=v/v_{c}}
T
r
=
T
/
T
c
{\displaystyle T_{r}=T/T_{c}}
c
{\displaystyle c}
∂
p
/
∂
v
|
T
=
0
{\displaystyle \partial p/\partial v|_{T}=0}
∂
2
p
/
∂
v
2
|
T
=
0
{\displaystyle \partial ^{2}p/\partial _{v^{2}}|_{T}=0}
p
∗
/
p
c
=
27
{\displaystyle p^{*}/p_{c}=27}
v
∗
/
v
c
=
1
/
3
{\displaystyle v^{*}/v_{c}=1/3}
T
∗
/
T
c
=
27
/
8
{\displaystyle T^{*}/T_{c}=27/8}
c
{\displaystyle c}
v
r
{\displaystyle v_{r}}
v
r
{\displaystyle v_{r}}
パラメータ範囲 にわたって 、 は から単調に減少し 、の 極限でのように 0 に近づきます 。したがって 、 の極限で および 、 およびです 。 および の挙動は 、 方程式から得られます。これらの特性は両方とも およびから単調に減少し 、の 極限で のように 0 に近づきます 。これらのことから 、ファンデルワールス飽和蒸気はこの極限では理想気体であることに注意してください。ゾンマーフェルトの言葉を言い換えると、ファンデルワールスによる理論によって、 飽和蒸気が理想気体のように振る舞うとき、実在する気体の飽和蒸気がまさにこのように振る舞うことが予測できるのは注目すべきことです。
0
≤
y
<
∞
{\displaystyle 0\leq y<\infty }
f
(
y
)
{\displaystyle f(y)}
f
(
0
)
=
1
{\displaystyle f(0)=1}
f
(
y
)
∼
2
(
y
−
1
)
e
−
y
{\displaystyle f(y)\sim 2(y-1)e^{-y}}
y
→
∞
{\displaystyle y\rightarrow \infty }
ϱ
f
(
0
)
=
ϱ
g
(
0
)
=
1
/
2
{\displaystyle \varrho _{f}(0)=\varrho _{g}(0)=1/2}
y
→
∞
{\displaystyle y\rightarrow \infty }
ϱ
f
(
y
)
∼
2
(
y
−
1
)
{\displaystyle \varrho _{f}(y)\sim 2(y-1)}
ϱ
g
(
y
)
∼
2
(
y
−
1
)
e
−
2
y
{\displaystyle \varrho _{g}(y)\sim 2(y-1)e^{-2y}}
p
s
(
y
)
{\displaystyle p_{s}(y)}
T
s
(
y
)
{\displaystyle T_{s}(y)}
p
s
(
0
)
=
p
∗
/
27
{\displaystyle p_{s}(0)=p^{*}/27}
T
s
(
0
)
=
8
T
∗
/
27
{\displaystyle T_{s}(0)=8T^{*}/27}
p
s
∼
p
∗
ϱ
g
/
ϱ
f
{\displaystyle p_{s}\sim p^{*}\varrho _{g}/\varrho _{f}}
T
s
∼
T
∗
/
ϱ
f
{\displaystyle T_{s}\sim T^{*}/\varrho _{f}}
y
→
∞
{\displaystyle y\rightarrow \infty }
p
s
∼
(
p
∗
/
T
∗
)
T
s
ϱ
g
∼
(
p
∗
v
∗
/
T
∗
)
T
s
ρ
g
=
ρ
g
R
T
s
{\displaystyle p_{s}\sim (p^{*}/T^{*})T_{s}\varrho _{g}\sim (p^{*}v^{*}/T^{*})T_{s}\rho _{g}=\rho _{g}RT_{s}}
s
g
−
s
f
≫
2
R
{\displaystyle s_{g}-s_{f}\gg 2R}
さらに、 液体のスピノーダル点は負圧で発生し、 この点を説明するために図 4 に等温線が示されています。これは、これらの液体の準安定状態の一部が張力を受けており、温度が低いほど引張応力が大きくなることを意味します。これは直感に反するように思えますが、状況によっては液体が張力を支えることができることが知られています。Tien と Lienhard [42] は これに注目し、次のように書いています。
T
r
<
27
/
32
=
0.84375
{\displaystyle T_{r}<27/32=0.84375}
T
r
=
0.8
{\displaystyle T_{r}=0.8}
ファンデルワールス方程式は、低温では液体が非常に大きな 張力 を受けると予測しているが、この事実から、この方程式を軽視する研究者もいる。近年、この方程式が完全に正しいことを示す測定が行われている。 [43] 溶解したガスのないきれいな液体は、より大きな張力を受ける可能性がある 。
p
c
{\displaystyle p_{c}}
これはファンデルワールス理論のもう一つの興味深い特徴です。
注記
^ カレン、pp 131-135、(1960)
^ ファンデルワールス、p. 174、(1988)
^ エプスタイン、9ページ
^ エプスタイン、10ページ
^ カレン、pp 146-163
^ グッドスタイン、pp 443-452
^ コンデプディ、ディリップ、プリゴジン、イリヤ(2014-12-31)。現代の熱力学:熱機関から散逸構造まで。ジョン・ワイリー&サンズ 。ISBN 9781118371817 。
^ カレン、163ページ
^ リーンハルト IV およびリーンハルト V、pp 467-469
^ リーンハルト、169-187
^ ゾンマーフェルト、66ページ
^ カレン、146ページ
^ グッドスタイン、p 443
^ マクスウェル、JC、pp. 357-359
^ ヴァン・ワイレンとソンタグ、p. 211
^ モランとシャピロ、HN、p. 251
^ ギブス、pp. 62-65
^ ab Callen、pp. 99より
^ ギブス、96-100ページ
^ カレン、pp 163-167
^ カレン、98-100ページ
^ カレン、150ページ
^ カレン、172-173ページ
^ ゾンマーフェルト、A.、67-68 ページ、(1956)
^ カレン、105ページ
^ カレン、120ページ
^ カレン、135ページ
^ レクトリス、429ページ
^ カレン、p.98
^ ファンデルワールス、246ページ
^ デ・ブール、p.8
^ ギブス、pp.87-90
^ ファンデルワールス、pp. 246-247
^ デボア、pp.8-16
^ ウェールズ、デイビッド、ウェールズ (2003)。エネルギーランドスケープ:クラスター、生体分子、ガラスへの応用。ケンブリッジ大学出版局。p. 444。ISBN 9780521814157 。
^ 「一次相転移とスピノーダル分解のダイナミクス」 www.mhkoepf.de . 2019年11月12日 閲覧 。
^ シャムスンダルとリーンハルト、876-880
^ バルフェットとユーバンク、pp. 168-175
^ レックナー、161ページ
^ ジョンストン、16-18ページ
^ ゾンマーフェルト、A.、pp 56-57
^ ティエンとリーンハルト、p.254
^ Temperley, HNV, 静水圧張力下における水の挙動: III , Proc. Phys. Soc. (ロンドン), 59 , p 199, (1947)
参考文献
Barrufet, MA; Eubank, PT (1989). 「三次状態方程式による純粋流体の一般化された飽和特性」. 化学工学教育 . 23 (3): 168–175.
Callen, HB (1960). 熱力学 . NY: John Wiley and Sons.
DeBoer, J. (1974). 「当時のファンデルワールスと現在の復活の開会演説」. Physica . 73 : 1–27. doi :10.1016/0031-8914(74)90223-7.
エプスタイン、PS (1937)。 『熱力学の教科書』 。ニューヨーク:ジョン・ワイリー・アンド・サンズ。
ギブス、JW (1928) [1876, 1878]。「不均質物質の平衡について」。ロングリー、WR、ヴァン・ネーム、RG (編)。J . ウィラード・ギブス著作集 第 1 巻 熱力学 (1948 版)。ニューヘブン: イェール大学出版局。pp. 55–353。
グッドスタイン、DL (1985)[1975]。 物質の状態 。ニューヨーク:ドーバー。
ジョンストン、 DC ( 2014)。 ファンデルワールス流体の熱力学の進歩 。arXiv : 1402.1205。doi : 10.1088 / 978-1-627-05532-1。ISBN 978-1-627-05532-1 。
Lekner, J. (1982). 「ファンデルワールス液体-蒸気共存曲線のパラメトリック解」 Am. J. Phys . 50 (2): 161–163. doi :10.1119/1.12877.
Lienhard IV, JH ; Lienhard V (2019). 熱伝達教科書 (第 5 版). ミネオラ、ニューヨーク州: Dover Publications Inc.
Lienhard, JH (1986). 「過熱液体の特性と挙動」. Lat. Am. J. Heat and Mass Transfer . 10 : 169–187.
マクスウェル 、JC ( 1875)。「物体の分子構成の動的証拠について」。 ネイチャー 。11 (279): 357–359。doi :10.1038/011357a0。
Moran, MJ; Shapiro, HN (2000). 工学熱力学の基礎 (第 4 版). NY: McGraw-Hill.
Rectorys, K. (1969) 「実変数関数の微分積分」 Rectorys, K. (編) 『 応用数学の調査』 マサチューセッツ州ケンブリッジ: MIT 出版。pp. 397–439。
シャムスンダー、N. JH リーンハルト (1983)。 「ファンデルワールス流体の飽和特性と準安定特性」。 カナダ J Chem Eng . 61 (6): 876–880。 土井 :10.1002/cjce.5450610617。
Sommerfeld, A. (1956). Bopp, F.; Meixner, J. (編). 熱力学と統計力学 - 理論物理学講義第5巻 . 翻訳: Kestin, J. NY: Academic Press.
Tien, CL ; Lienhard, JH (1979). 統計熱力学改訂版 . NY: Hemisphere Publishing.
Van Wylen, ジョージア州 ;レバノン州ソンタグ (1973 年)。 古典熱力学の基礎 (第 2 版)。ニューヨーク:ジョン・ワイリー・アンド・サンズ。
ファンデルワールス、JD (1984)。 ローリンソン、JS (編)。 気体と液体の状態の連続性について、JS ローリンソンによる編集と序文付き 。ニューヨーク: ドーバーフェニックスエディション。