熱力学 では、 残差特性は、 同じ 密度、温度、組成で考慮された実際の流体特性と理想気体特性 の 差 として 定義され、通常は次のように表され
ます 。
バツ
(
T
、
五
、
ん
)
=
バツ
私
d
(
T
、
五
、
ん
)
+
バツ
r
e
s
(
T
、
五
、
ん
)
{\displaystyle X(T,V,n)=X^{id}(T,V,n)+X^{res}(T,V,n)}
ここで、 は 与えられた温度、体積、モル数における 熱力学的性質 であり、 は 理想気体 の特性値であり 、 は 残余特性である。基準状態は、通常、理想気体の値への寄与に組み込まれ、次のように表される。
バツ
{\displaystyle X}
バツ
私
d
{\displaystyle X^{id}}
バツ
r
e
s
{\displaystyle X^{res}}
バツ
私
d
(
T
、
五
、
ん
)
=
バツ
∘
、
私
d
(
T
、
ん
)
+
Δ
私
d
バツ
(
T
、
五
、
ん
)
{\displaystyle X^{id}(T,V,n)=X^{\circ ,id}(T,n)+\Delta _{id}X(T,V,n)}
ここで、 は 基準状態(通常は 1 bar における純粋な理想気体種)における の値であり、 は この基準状態からの
における理想気体のプロパティの偏差です。
バツ
∘
、
私
d
{\displaystyle X^{\circ,id}}
バツ
{\displaystyle X}
Δ
私
d
バツ
{\displaystyle \Delta _{id}X}
(
T
、
五
、
ん
)
{\displaystyle (T,V,n)}
残留特性を過剰特性 と混同しないでください。過剰特性は 、通常、理想気体システムではない何らかの参照システムからの熱力学的特性の偏差として定義されます。過剰特性と過剰モデル ( 活動係数モデルとも呼ばれる) は通常、 精錬所 、 ポリマーブレンド 、 電解質 などの厳密に液相のシステムに関係しますが 、残留特性は、 気液平衡 が優勢なシステム、またはガスと液体の両方が対象となるシステムをモデル化するために一般的に使用される 状態方程式 と密接に関連しています。一部のアプリケーションでは、活動係数モデルと状態方程式は、 活動係数 と フガシティーを 表すために一般的に使用される記号を指す「 - モデル」(読み: ガンマ-ファイ) と呼ばれるものに統合されます 。
γ
{\displaystyle \gamma}
ϕ
{\displaystyle \phi}
意義
状態方程式 の開発と実装において 、残余特性の概念は有用である。これは、非理想性から生じる流体の挙動と理想気体の特性から生じる挙動を区別することを可能にするからである。例えば、等積 熱容量は 次のように与えられる
。
C
五
=
(
∂
あなた
∂
T
)
五
、
ん
=
T
(
∂
ス
∂
T
)
五
、
ん
=
T
[
(
∂
ス
私
d
∂
T
)
五
、
ん
+
(
∂
ス
r
e
s
∂
T
)
五
、
ん
]
=
C
五
私
d
+
C
五
r
e
s
{\displaystyle C_{V}=\left({\frac {\partial U}{\partial T}}\right)_{V,n}=T\left({\frac {\partial S}{\partial T}}\right)_{V,n}=T\left[\left({\frac {\partial S^{id}}{\partial T}}\right)_{V,n}+\left({\frac {\partial S^{res}}{\partial T}}\right)_{V,n}\right]=C_{V}^{id}+C_{V}^{res}}
ここで、理想気体の熱容量は、 非常に低い圧力で熱容量を測定することによって実験的に測定できます。測定後は通常、 Shomate方程式などの 多項式 を使用して表されます。残留熱容量は次のように表されます
。
C
五
私
d
{\displaystyle C_{V}^{id}}
T
(
∂
ス
r
e
s
∂
T
)
五
、
ん
=
−
T
(
∂
2
あ
r
e
s
∂
T
2
)
五
、
ん
{\displaystyle T\left({\frac {\partial S^{res}}{\partial T}}\right)_{V,n}=-T\left({\frac {\partial ^{2}A^{res}}{\partial T^{2}}}\right)_{V,n}}
、
理想的な寄与は状態方程式とは無関係であるため、熱容量を予測または相関させる際の特定の状態方程式の精度は、残差寄与のみを考慮することによって評価できます。
均衡計算において
流体の相平衡(すなわち、液体-蒸気平衡または液体-液体平衡)において、 フガシティ係数 の概念は極めて重要である。相 、、、 … からなるシステムの平衡条件は、 化学平衡の条件は、
α
{\displaystyle \alpha}
β
{\displaystyle \beta}
γ
{\displaystyle \gamma}
x
私
α
Φ
私
α
=
x
私
β
Φ
私
β
=
x
私
γ
Φ
私
γ
=
。
。
。
{\displaystyle x_{i}^{\alpha }\Phi _{i}^{\alpha }=x_{i}^{\beta }\Phi _{i}^{\beta }=x_{i}^{\gamma }\Phi _{i}^{\gamma }=\quad ...}
すべての種について 、 は 相 における 種のモル分率 、 は 相 における 種のフガシティ係数である 。フガシティ係数は次のように定義される。
私
{\displaystyle i}
x
私
じゅう
{\displaystyle x_{i}^{j}}
私
{\displaystyle i}
じゅう
{\displaystyle j}
Φ
私
じゅう
{\displaystyle \Phi_{i}^{j}}
私
{\displaystyle i}
じゅう
{\displaystyle j}
μ
私
=
μ
私
∘
+
R
T
行
Φ
私
x
私
p
p
∘
{\displaystyle \mu _{i}=\mu _{i}^{\circ }+RT\ln {\frac {\Phi _{i}x_{i}p}{p^{\circ }}}}
残留化学ポテンシャルに直接関係しており、
μ
私
=
μ
私
私
d
+
μ
私
r
e
s
=
μ
私
∘
+
R
T
行
x
私
p
p
∘
+
μ
私
r
e
s
⟹
μ
私
r
e
s
=
R
T
行
Φ
私
{\displaystyle \mu _{i}=\mu _{i}^{id}+\mu _{i}^{res}=\mu _{i}^{\circ }+RT\ln {\frac {x_{i}p}{p^{\circ }}}+\mu _{i}^{res}\implies \mu _{i}^{res}=RT\ln \Phi _{i}}
、
したがって、 であるため、総ヘルムホルツエネルギーではなく、残差 ヘルムホルツエネルギー を正確に記述することが、システムの平衡状態を正確に計算するための鍵である
ことがわかります。
μ
私
r
e
s
=
(
∂
あ
r
e
s
∂
ん
私
)
T
、
五
{\displaystyle \mu _{i}^{res}=\left({\frac {\partial A^{res}}{\partial n_{i}}}\right)_{T,V}}
残余エントロピースケーリング
流体の 残留 エントロピーには特別な意味があります。1976 年に Yasha Rosenfeld は 画期的な論文を発表し、純粋な液体の 輸送係数を 残留エントロピーの関数として表すと、2 つの変数 (温度と圧力、または温度と密度) の関数としてではなく、単変量関数として扱うことができることを示しました。 [1] この発見は残留エントロピー スケーリングの概念につながり、今日に至るまで大量の研究を促し、輸送係数を残留エントロピーの関数としてモデル化するさまざまなアプローチが研究されてきました。 [2] 残留エントロピー スケーリングは、現在でも非常に活発に研究されている分野です。
変数セットへの依存
実状態( )における 任意の実状態変数 は、 を評価する か を 評価するかに依存しないが 、残差特性は一般に変数セットに依存することに注意する必要がある。すなわち、
バツ
{\displaystyle X}
T
、
五
、
p
、
ん
{\displaystyle T,V,p,n}
バツ
(
T
、
p
、
ん
)
{\displaystyle X(T,p,n)}
バツ
(
T
、
五
、
ん
)
{\displaystyle X(T,V,n)}
バツ
r
e
s
(
T
、
p
、
ん
)
≠
バツ
r
e
s
(
T
、
五
、
ん
)
{\displaystyle X^{res}(T,p,n)\neq X^{res}(T,V,n)}
これは、実状態が一般に有効な 理想気体 状態ではないという事実から生じ 、特性の理想部分は変数セットに応じて異なります。たとえば、純粋な流体の 化学ポテンシャル を考えてみましょう。この状態は 理想気体の法則を満たしていませんが、実際の流体の実状態である可能性があります。温度、圧力、モル数の関数として計算された理想気体の化学ポテンシャルは、
(
T
、
五
、
p
、
ん
)
{\displaystyle (T,V,p,n)}
(
T
、
五
、
p
、
ん
)
{\displaystyle (T,V,p,n)}
μ
私
d
(
T
、
p
、
ん
)
=
μ
∘
+
R
T
行
p
p
∘
{\displaystyle \mu^{id}(T,p,n)=\mu^{\circ}+RT\ln {\frac{p}{p^{\circ}}}
、
これを濃度( )の関数として計算すると 、
c
=
ん
/
五
{\displaystyle c=n/V}
μ
私
d
(
T
、
五
、
ん
)
=
μ
∘
+
R
T
行
c
c
∘
{\displaystyle \mu^{id}(T,V,n)=\mu^{\circ}+RT\ln {\frac{c}{c^{\circ}}}
、
そのような
μ
私
d
(
T
、
p
、
ん
)
−
μ
私
d
(
T
、
五
、
ん
)
=
R
T
行
p
p
∘
−
R
T
行
c
c
∘
=
R
T
行
p
五
ん
R
T
=
R
T
行
ず
{\displaystyle \mu^{id}(T,p,n)-\mu^{id}(T,V,n)=RT\ln {\frac {p}{p^{\circ }}}-RT\ln {\frac {c}{c^{\circ }}}=RT\ln {\frac {pV}{nRT}}=RT\ln Z}
、
ここで 、 は 圧縮率を表す。この結果、
p
∘
=
c
∘
R
T
{\displaystyle p^{\circ}=c^{\circ}RT}
ず
{\displaystyle Z}
μ
私
(
T
、
p
、
ん
)
−
μ
私
(
T
、
五
、
ん
)
=
0
⟹
μ
私
r
e
s
(
T
、
五
、
ん
)
−
μ
私
r
e
s
(
T
、
p
、
ん
)
=
R
T
行
ず
{\displaystyle \mu _{i}(T,p,n)-\mu _{i}(T,V,n)=0\implies \mu _{i}^{res}(T,V,n)-\mu _{i}^{res}(T,p,n)=RT\ln Z}
。
実用計算
実際には、最も重要な残差特性は残差 ヘルムホルツエネルギー である。その理由は、他の残差特性は残差ヘルムホルツエネルギーから様々な導関数として計算できるからである( マクスウェル関係を 参照)。
(
∂
あ
∂
五
)
T
、
ん
=
(
∂
あ
私
d
∂
五
)
T
、
ん
+
(
∂
あ
r
e
s
∂
五
)
T
、
ん
⟺
(
∂
あ
r
e
s
∂
五
)
T
、
ん
=
(
∂
あ
∂
五
)
T
、
ん
−
(
∂
あ
私
d
∂
五
)
T
、
ん
=
−
p
(
T
、
五
、
ん
)
−
(
−
p
私
d
(
T
、
五
、
ん
)
)
{\displaystyle \left({\frac {\partial A}{\partial V}}\right)_{T,n}=\left({\frac {\partial A^{id}}{\partial V}}\right)_{T,n}+\left({\frac {\partial A^{res}}{\partial V}}\right)_{T,n}\iff \left({\frac {\partial A^{res}}{\partial V}}\right)_{T,n}=\left({\frac {\partial A}{\partial V}}\right)_{T,n}-\left({\frac {\partial A^{id}}{\partial V}}\right)_{T,n}=-p(T,V,n)-(-p^{id}(T,V,n))}
さらに 、あらゆる流体は無限体積の限界において理想気体に減少するため、
A
r
e
s
(
T
,
V
′
,
n
)
−
A
r
e
s
(
T
,
V
=
∞
,
n
)
=
∫
∞
V
′
(
∂
A
r
e
s
∂
V
)
d
V
=
∫
∞
V
′
p
i
d
(
T
,
V
,
n
)
−
p
(
T
,
V
,
n
)
d
V
{\displaystyle A^{res}(T,V',n)-A^{res}(T,V=\infty ,n)=\int _{\infty }^{V'}\left({\frac {\partial A^{res}}{\partial V}}\right)dV=\int _{\infty }^{V'}p^{id}(T,V,n)-p(T,V,n)dV}
A
(
T
,
V
=
∞
,
n
)
=
A
i
d
(
T
,
V
=
∞
,
n
)
⟺
A
r
e
s
(
T
,
V
=
∞
,
n
)
=
0
{\displaystyle A(T,V=\infty ,n)=A^{id}(T,V=\infty ,n)\iff A^{res}(T,V=\infty ,n)=0}
。
したがって、 ファンデルワールス状態方程式 のような 圧力 が明示的な 状態方程式 については、次のように計算できる。
A
r
e
s
(
T
,
V
,
n
)
=
∫
∞
V
n
R
T
V
′
−
p
(
T
,
V
′
,
n
)
d
V
′
{\displaystyle A^{res}(T,V,n)=\int _{\infty }^{V}{\frac {nRT}{V'}}-p(T,V',n)dV'}
。
しかし、 SAFT などの状態方程式を展開する現代的なアプローチでは、圧力について明示的な方程式を展開するよりも、 の方程式を直接展開することによって状態方程式を展開する方が簡単であることがわかっています 。
A
r
e
s
{\displaystyle A^{res}}
参考文献
^ Rosenfeld, Yaakov (1977-06-01). 「単純なシステムの輸送係数と内部エントロピーの関係」. Physical Review A. 15 ( 6): 2545–2549. doi :10.1103/PhysRevA.15.2545.
^ Bell, Ian H.; Hellmann, Robert; Harvey, Allan H. ( 2020-03-12 ). 「輸送特性の残差エントロピースケーリングのゼロ密度限界」。Journal of Chemical & Engineering Data。65 ( 3): 1038–1050。doi : 10.1021/acs.jced.9b00455。ISSN 0021-9568 。
JM Smith、HCVan Ness、MM Abbot 化学工学熱力学入門 2000、McGraw-Hill 第 6 版 ISBN 0-07-240296-2
ロバート・ペリー、ドン・W・グリーン ペリーの化学技術者ハンドブック 2007 McGraw-Hill 第 8 版 ISBN 0-07-142294-3