在物理学 和热力学 中，状态方程 （英语：Equation of state 物态方程 ，表达了热力学系统中若干个态函数 参量之间的关系。[ 1] 温度 、压强 、体积 和内能 。状态方程在描述流体 、混合流体、固体甚至是研究恒星 内部都十分有用。
物态方程最显著的作用是它可以通过已知条件来预测气体和液体的状态。具有这样的目的的、较为简单的物态方程是理想气体状态方程 ，这个方程在低压强、中等温度的情况下，可以粗略地估算气体的状态。然而，当压强升、温度降低时，这个方程的精确度会降低，并且不能预测气体液化为液体。因此，科学家研究出了一系列针对气体和液体的更为精确的物态方程。目前，没有一个单独的方程可以精确地预测所有物质在所有情况下的状态。
除了上述的针对气体和液体的物态方程，也有可以预测固体体积，甚至可以预测固体从一种晶态变为另一种的转变。对于恒星、中子星 内部，也有专门的模型来描述其物态变化。与之相关的则是理想流体的状态方程 。
玻意耳定律 可能是第一个有关物质状态的方程表达。1662年，这位著名的爱尔兰物理学家、化学家罗伯特·波义耳 利用一端封闭、管内装有气体的J形管进行了一系列实验。管里还加入了水银 封闭，使管内气体的数量保持恒定。然后，波义耳仔细地测量了气体的体积。气体的压强可以通过测量J形管长短两端水银页面的高度来得到。通过这一系列实验，波义耳注意到气体体积与压强成反比关系，这可以表达为以下的数学形式：
p
V
=
c
o
n
s
t
a
n
t
.
{\displaystyle pV=\mathrm {constant} .\,\!}
以上的关系式也被埃德姆·马略特 发现，这个定律也被称为“玻意耳－马略特定律”。不过，马略特直到1676年都没有发表这个发现。
1787年，法国物理学家雅克·查理 （Jacques Charles ）发现氧气 、氮气 、氢气 、二氧化碳 和其他气体，当温度升高80开尔文（热力学温标 ）、其他条件一致时，气体膨胀的体积相等。约瑟夫·路易·盖-吕萨克 也发表了相似的实验结果，表明了体积与温度之间的线性关系：
V
1
T
1
=
V
2
T
2
.
{\displaystyle {\frac {V_{1}}{T_{1}}}={\frac {V_{2}}{T_{2}}}.}
上面这个方程被称为“查理定律 ”。
道尔顿分压定律 指出，
n
{\displaystyle n}
p
total
=
p
1
+
p
2
+
⋯
+
p
n
=
p
total
=
∑
i
=
1
n
p
i
.
{\displaystyle p_{\text{total}}=p_{1}+p_{2}+\cdots +p_{n}=p_{\text{total}}=\sum _{i=1}^{n}p_{i}.}
理想气体定律，又称理想气体状态方程 ，是1834年法国物理学家埃米里·克拉珀龙（法语：Benoît Paul Émile Clapeyron ）综合波义耳和查里的发现得到的一个方程。最初，这个定律被表达为
P
V
m
=
R
(
T
C
+
267
)
{\displaystyle PV_{m}=R(T_{C}+267)}
R
{\displaystyle R}
气体常数 。然而，不久后的研究表明267应该被修正为接近273.2，因为当时定义0摄氏度等于273.15开尔文，于是有以下的式子：
p
V
m
=
R
(
T
C
+
273.15
)
.
{\displaystyle \ pV_{m}=R(T_{C}+273.15).}
1873年，荷兰物理学家范德华 给出了有限体积被连续分子充满情况下的物态方程。[ 2] 雷德利希-邝氏方程 （英语：Redlich–Kwong equation of state ）等对三次方程的修正。
对于一个系统中的给提物质，温度、体积和压强不是相互独立的量。它们被下面的方程约束：
f
(
p
,
V
,
T
)
=
0
.
{\displaystyle {\ f(p,V,T)=0}.}
式中的变量定义在下面给出。一些始终被用到的单位还将被使用，尽管国际单位制 （SI ）现在被推荐使用。特别指出，温度采用的是热力学温标。
p
{\displaystyle \ p}
V
{\displaystyle \ V}
n
{\displaystyle \ n}
V
m
{\displaystyle \ V_{m}}
V
n
{\displaystyle {\frac {V}{n}}}
摩尔体积 摩尔 物质（通常指气体或液体）的体积
T
{\displaystyle \ T}
R
{\displaystyle \ R}
p
c
{\displaystyle \ p_{c}}
V
c
{\displaystyle \ V_{c}}
T
c
{\displaystyle \ T_{c}}
经典理想气体状态方程可以写作：
p
V
=
n
R
T
.
{\displaystyle {\ pV=nRT}.}
也可以表达为以下形式：
p
=
ρ
(
γ
−
1
)
e
{\displaystyle {\ p=\rho (\gamma -1)e}}
这里，
ρ
{\displaystyle \rho }
密度 ，
γ
=
C
p
/
C
v
{\displaystyle \gamma =C_{p}/C_{v}}
绝热指数 ，
e
=
C
v
T
{\displaystyle e=C_{v}T}
C
v
{\displaystyle C_{v}}
C
p
{\displaystyle C_{p}}
范德华方程 写作：
(
p
+
a
V
m
2
)
(
V
m
−
b
)
=
R
T
{\displaystyle {\left(p+{\frac {a}{V_{m}^{2}}}\right)\left(V_{m}-b\right)=RT}}
这里
V
m
{\displaystyle V_{m}}
a
{\displaystyle a}
b
{\displaystyle b}
临界点 时的参数
p
c
,
T
c
{\displaystyle p_{c},T_{c}}
V
c
{\displaystyle V_{c}}
V
c
{\displaystyle V_{c}}
a
=
3
p
c
V
c
2
{\displaystyle a=3p_{c}\,V_{c}^{2}}
b
=
V
c
3
.
{\displaystyle b={\frac {V_{c}}{3}}.}
也可以写作：
a
=
27
(
R
T
c
)
2
64
p
c
{\displaystyle a={\frac {27(R\,T_{c})^{2}}{64p_{c}}}}
b
=
R
T
c
8
p
c
.
{\displaystyle b={\frac {R\,T_{c}}{8p_{c}}}.}
范德华方程于1873年被提出来，是最早的几个能够比理想气体状态方程更为准确的物态方程之一。在这一里程碑意义的方程中，
a
{\displaystyle a}
attraction parameter ），而
b
{\displaystyle b}
repulsion parameter ）或有效摩尔体积（effective molecular volume ）。尽管这个方程确实比理想气体状态方程更加优秀，并且确实能够预测液体的相态 ，其与实验值的较高吻合度却仅仅发生在有液体形成的情况。由于历史原因，范德华方程在教科书和论文中不断被引用，以现在的眼光评判的话，范德华方程仍然是过时的。其他一些稍微复杂一点的物态方程表现出更高的精确度。
范德华方程可以被认为是一个改进的理想气体状态方程，原因如下：
分子 被看做是具有体积的微粒，而不是一个根据其物质决定的点。因此
V
{\displaystyle V}
V
−
b
{\displaystyle V-b}
V
{\displaystyle V}
尽管理想气体分子间不进行相互碰撞，我们考虑分子与其他分子互相吸引的距离在几个分子半径的距离内。这对分子内的物质影响不大，但是却使分子表面的物质在半径方向被向里吸引。因此，利用理想气体状态方程计算式，我们人为地使得外层表面的压强减小，因此我们使用
p
{\displaystyle p}
p
{\displaystyle p}
p
{\displaystyle p}
ρ
2
{\displaystyle \rho ^{2}}
1
V
m
2
{\displaystyle {\frac {1}{V_{m}^{2}}}}
利用这些无量纲的状态参量，例如，Vr =Vm /Vc ，Pr =P/Pc 以及Tr =T/Tc ，范德华方程可以表示为：
(
P
r
+
3
V
r
2
)
(
3
V
r
−
1
)
=
8
T
r
{\displaystyle {\left(P_{r}+{\frac {3}{V_{r}^{2}}}\right)\left(3V_{r}-1\right)=8T_{r}}}
这个方程形式的好处在于，对于给定的Tr 和Pr ，无量纲体积可以直接利用三次方程 解得：
V
r
3
−
(
1
3
+
8
T
r
3
P
r
)
V
r
2
+
3
V
r
P
r
−
1
P
r
=
0
{\displaystyle {V_{r}^{3}-\left({\frac {1}{3}}+{\frac {8T_{r}}{3P_{r}}}\right)V_{r}^{2}+{\frac {3V_{r}}{P_{r}}}-{\frac {1}{P_{r}}}=0}}
对于Pr <1且Tr <1，系统处于蒸发-液态平衡。在上面的情况里，这个三次方程可以导出3个根。最大的根以及最小的根分别为气体和液体的无量纲体积。
[ 编辑 ]
p
=
R
T
V
m
−
b
−
a
T
V
m
(
V
m
+
b
)
{\displaystyle {p={\frac {R\,T}{V_{m}-b}}-{\frac {a}{{\sqrt {T}}\,V_{m}\left(V_{m}+b\right)}}}}
a
=
0.42748
R
2
T
c
2.5
p
c
{\displaystyle a={\frac {0.42748\,R^{2}\,T_{c}^{\,2.5}}{p_{c}}}}
b
=
0.08662
R
T
c
p
c
{\displaystyle b={\frac {0.08662\,R\,T_{c}}{p_{c}}}}
雷德利希-邝氏方程于1949年被提出，是对于当时其他物态方程的一个较大的改进形式。由于它相对简单的数学形式，它仍然引起了很多人的兴趣。虽然比范德华方程更优越，但是因为它在涉及液相的情况里表现得不太好，因此不能用于精确计算气液平衡。然而，它仍然可以在涉及单独的液相的情况里发挥作用。
当压强对临界压强的比值小于温度对临界温度比值的一半时，即：
p
p
c
<
T
2
T
c
.
{\displaystyle {\frac {p}{p_{c}}}<{\frac {T}{2T_{c}}}.}
雷德利希-邝氏方程能够充分胜任计算气相属性。
p
=
R
T
V
m
−
b
−
a
α
V
m
(
V
m
+
b
)
{\displaystyle p={\frac {R\,T}{V_{m}-b}}-{\frac {a\,\alpha }{V_{m}\left(V_{m}+b\right)}}}
a
=
0.42748
R
2
T
c
2
P
c
{\displaystyle a={\frac {0.42748\,R^{2}\,T_{c}^{2}}{P_{c}}}}
b
=
0.08664
R
T
c
P
c
{\displaystyle b={\frac {0.08664\,R\,T_{c}}{P_{c}}}}
α
=
(
1
+
(
0.48508
+
1.55171
ω
−
0.15613
ω
2
)
(
1
−
T
r
0.5
)
)
2
{\displaystyle \alpha =\left(1+\left(0.48508+1.55171\,\omega -0.15613\,\omega ^{2}\right)\left(1-T_{r}^{\,0.5}\right)\right)^{2}}
T
r
=
T
T
c
{\displaystyle T_{r}={\frac {T}{T_{c}}}}
这里，
ω
{\displaystyle \omega }
偏离系数 。索阿维修正的原始形式为：
α
=
(
1
+
(
0.48
+
1.574
ω
−
0.176
ω
2
)
(
1
−
T
r
0.5
)
)
2
{\displaystyle \alpha =\left(1+\left(0.48+1.574\,\omega -0.176\,\omega ^{2}\right)\left(1-T_{r}^{\,0.5}\right)\right)^{2}}
在1972年，索阿维(G. Soave ）[ 3] T )项替换为一个与温度和偏离系数有关的函数α(T,ω)，替换后的方程也被称作“索阿维-雷德利希-邝氏方程”。引入α函数是为了与烃 的蒸发压强相适应，事实证明这个方程在这些物质的相关计算中更加精确。
[ 编辑 ]
p
=
R
T
V
m
−
b
−
a
α
V
m
2
+
2
b
V
m
−
b
2
{\displaystyle p={\frac {R\,T}{V_{m}-b}}-{\frac {a\,\alpha }{V_{m}^{2}+2bV_{m}-b^{2}}}}
a
=
0.45724
R
2
T
c
2
p
c
{\displaystyle a={\frac {0.45724\,R^{2}\,T_{c}^{2}}{p_{c}}}}
b
=
0.07780
R
T
c
p
c
{\displaystyle b={\frac {0.07780\,R\,T_{c}}{p_{c}}}}
α
=
(
1
+
(
0.37464
+
1.54226
ω
−
0.26992
ω
2
)
(
1
−
T
r
0.5
)
)
2
{\displaystyle \alpha =\left(1+\left(0.37464+1.54226\,\omega -0.26992\,\omega ^{2}\right)\left(1-T_{r}^{\,0.5}\right)\right)^{2}}
T
r
=
T
T
c
{\displaystyle T_{r}={\frac {T}{T_{c}}}}
可以写作多项式的形式：
A
=
a
α
p
R
2
T
2
{\displaystyle A={\frac {a\alpha p}{R^{2}\,T^{2}}}}
B
=
b
p
R
T
{\displaystyle B={\frac {bp}{RT}}}
Z
3
−
(
1
−
B
)
Z
2
+
(
A
−
2
B
−
3
B
2
)
Z
−
(
A
B
−
B
2
−
B
3
)
=
0
{\displaystyle Z^{3}-(1-B)\ Z^{2}+(A-2B-3B^{2})\ Z-(AB-B^{2}-B^{3})=0\,\!}
这里，
ω
{\displaystyle \omega }
R
{\displaystyle R}
气体常数 ，
Z
=
P
V
R
T
{\displaystyle Z={\frac {PV}{RT}}}
压缩因子 。
彭-罗宾逊方程于1976年被导出，主要是为了以下的几个目的[ 4]
参数应该可以用临界参量和偏离系数表达；
模型能够在临界状态附近体现较高的精确度，尤其是对于压缩因子和流体密度的计算；
混杂的规定不应该引用比二元相互作用参量更加复杂的参数，而应该独立于温度、压强和物质组分；
方程应该能够适用于所有流体（例如天然气）的性质。 虽然它总的来说在预测许多流体密度方面表现得更好，但是在大多数部分，彭-罗宾逊方程的效果和索阿维的方程相似，尤其是在针对非极化物质的计算方面。
[ 编辑 ] Elliott, Suresh, Donohue方程（简称：ESD方程）于1990年被提出。[ 5]
p
V
m
R
T
=
Z
=
1
+
Z
r
e
p
+
Z
a
t
t
{\displaystyle {\frac {pV_{m}}{RT}}=Z=1+Z^{\rm {rep}}+Z^{\rm {att}}}
这里：
Z
r
e
p
=
4
c
η
1
−
1.9
η
{\displaystyle Z^{\rm {rep}}={\frac {4c\eta }{1-1.9\eta }}}
Z
a
t
t
=
−
z
m
q
η
Y
1
+
k
1
η
Y
{\displaystyle Z^{\rm {att}}=-{\frac {z_{m}q\eta Y}{1+k_{1}\eta Y}}}
且，
c
{\displaystyle c}
c
=
1
{\displaystyle c=1}
对于非球状分子，下面的关系被用到：
c
=
1
+
3.535
ω
+
0.533
ω
2
{\displaystyle c=1+3.535\omega +0.533\omega ^{2}}
ω
{\displaystyle \omega }
减少的分子数量密度
η
{\displaystyle \eta }
η
=
v
∗
n
V
{\displaystyle \eta ={\frac {v^{*}n}{V}}}
这里
v
∗
{\displaystyle v^{*}}
n
{\displaystyle n}
V
{\displaystyle V}
与尺寸相关的特征参数与形状参数
c
{\displaystyle c}
v
∗
=
k
T
c
P
c
Φ
{\displaystyle v^{*}={\frac {kT_{c}}{P_{c}}}\Phi }
这里
Φ
=
0.0312
+
0.087
(
c
−
1
)
+
0.008
(
c
−
1
)
2
1.000
+
2.455
(
c
−
1
)
+
0.732
(
c
−
1
)
2
{\displaystyle \Phi ={\frac {0.0312+0.087(c-1)+0.008(c-1)^{2}}{1.000+2.455(c-1)+0.732(c-1)^{2}}}}
k
{\displaystyle k}
波尔兹曼常数 。应用时，应当注意波尔兹曼常数和气体常数之间的关系，并要观察到分子数量可以用阿伏伽德罗常数 和摩尔质量 表示。简化的数量密度
η
{\displaystyle \eta }
η
=
R
T
c
P
c
Φ
1
V
m
.
{\displaystyle \eta ={\frac {RT_{c}}{P_{c}}}\Phi {\frac {1}{V_{m}}}.}
形状参数
q
{\displaystyle q}
Y
=
exp
(
ϵ
k
T
)
−
k
2
{\displaystyle Y=\exp \left({\frac {\epsilon }{kT}}\right)-k_{2}}
这里
ϵ
{\displaystyle \epsilon }
势阱 的“深度”，又可以由以下式子得到：
ϵ
k
=
1.000
+
0.945
(
c
−
1
)
+
0.134
(
c
−
1
)
2
1.023
+
2.225
(
c
−
1
)
+
0.478
(
c
−
1
)
2
{\displaystyle {\frac {\epsilon }{k}}={\frac {1.000+0.945(c-1)+0.134(c-1)^{2}}{1.023+2.225(c-1)+0.478(c-1)^{2}}}}
z
m
{\displaystyle z_{m}}
k
1
{\displaystyle k_{1}}
k
2
{\displaystyle k_{2}}
k
3
{\displaystyle k_{3}}
z
m
=
9.49
{\displaystyle z_{m}=9.49}
k
1
=
1.7745
{\displaystyle k_{1}=1.7745}
k
2
=
1.0617
{\displaystyle k_{2}=1.0617}
k
3
=
1.90476.
{\displaystyle k_{3}=1.90476.}
关于这个模型的细节，可以参阅J.R. Elliott Jr在1990年的论文。[ 5]
p
(
V
−
b
)
=
R
T
e
−
a
/
R
T
V
{\displaystyle \ p(V-b)=RTe^{-a/RTV}}
这里
a
{\displaystyle a}
b
{\displaystyle b}
可以简化为：
T
c
=
a
4
R
b
,
p
c
=
a
4
b
2
e
2
,
V
c
=
2
b
.
{\displaystyle \ T_{c}={\frac {a}{4Rb}},\ p_{c}={\frac {a}{4b^{2}e^{2}}},\ V_{c}=2b.}
p
V
m
R
T
=
1
+
B
V
m
+
C
V
m
2
+
D
V
m
3
+
…
{\displaystyle {\frac {pV_{m}}{RT}}=1+{\frac {B}{V_{m}}}+{\frac {C}{V_{m}^{2}}}+{\frac {D}{V_{m}^{3}}}+\dots }
B
=
−
V
c
{\displaystyle B=-V_{c}\,}
C
=
V
c
2
3
{\displaystyle C={\frac {V_{c}^{2}}{3}}}
尽管形式上不是最方便的物态方程，维里状态方程（Virial equation of state ）仍然十分重要，因为它可以直接由统计力学 推导出。这个方程也被称作昂内斯 （Heike Kamerlingh Onnes ）方程。如果关于分子内部力的数学假设设定得适当，那么就能得到每一个维里系数的理论表达。在上面的式子中，
B
{\displaystyle B}
C
{\displaystyle C}
B
,
C
,
D
{\displaystyle B,C,D}
通过应用范德瓦尔斯方程中的参数
a
,
b
{\displaystyle a,b}
B
=
0
{\displaystyle B=0}
B
=
b
−
a
R
T
.
{\displaystyle B=b-{\frac {a}{RT}}.}
p
=
ρ
R
T
+
(
B
0
R
T
−
A
0
−
C
0
T
2
+
D
0
T
3
−
E
0
T
4
)
ρ
2
+
(
b
R
T
−
a
−
d
T
)
ρ
3
+
α
(
a
+
d
T
)
ρ
6
+
c
ρ
3
T
2
(
1
+
γ
ρ
2
)
exp
(
−
γ
ρ
2
)
{\displaystyle p=\rho RT+\left(B_{0}RT-A_{0}-{\frac {C_{0}}{T^{2}}}+{\frac {D_{0}}{T^{3}}}-{\frac {E_{0}}{T^{4}}}\right)\rho ^{2}+\left(bRT-a-{\frac {d}{T}}\right)\rho ^{3}+\alpha \left(a+{\frac {d}{T}}\right)\rho ^{6}+{\frac {c\rho ^{3}}{T^{2}}}\left(1+\gamma \rho ^{2}\right)\exp \left(-\gamma \rho ^{2}\right)}
这里
p
{\displaystyle p}
ρ
{\displaystyle \rho }
1973年K.E. Starling出版的一部专著记有15种物质的参数值。[ 6]
亥姆霍兹物态方程（Multiparameter equations of state (MEOS) ）可以用来精确地计算纯净的流体（包括液态和气态）。亥姆霍兹函数为理想气体项和剩余项的和。这两项在简化的温度和密度很明晰，由此：
a
(
T
,
ρ
)
R
T
=
a
o
(
T
,
ρ
)
+
a
r
(
T
,
ρ
)
R
T
=
α
o
(
τ
,
δ
)
+
α
r
(
τ
,
δ
)
{\displaystyle {\frac {a(T,\rho )}{RT}}={\frac {a^{o}(T,\rho )+a^{r}(T,\rho )}{RT}}=\alpha ^{o}(\tau ,\delta )+\alpha ^{r}(\tau ,\delta )}
这里
τ
=
T
r
T
,
δ
=
ρ
ρ
r
{\displaystyle \tau ={\frac {T_{r}}{T}},\delta ={\frac {\rho }{\rho _{r}}}}
极端相对的流体具有以下的物态方程：
p
=
ρ
m
c
s
2
{\displaystyle p=\rho _{m}c_{s}^{2}}
这里
p
{\displaystyle p}
ρ
m
{\displaystyle \rho _{m}}
c
s
{\displaystyle c_{s}}
音速 。
对于理想的玻色气体 ，其物态方程为：
p
V
m
=
R
T
Li
α
+
1
(
z
)
ζ
(
α
)
(
T
T
c
)
α
{\displaystyle pV_{m}=RT~{\frac {{\textrm {Li}}_{\alpha +1}(z)}{\zeta (\alpha )}}\left({\frac {T}{T_{c}}}\right)^{\alpha }}
这里
α
{\displaystyle \alpha }
α
=
3
/
2
{\displaystyle \alpha =3/2}
z
{\displaystyle z}
exp
(
μ
/
k
T
)
{\displaystyle \exp(\mu /kT)}
μ
{\displaystyle \mu }
化学势 ，
L
i
{\displaystyle Li}
多重对数函数
，
ξ
{\displaystyle \xi }
黎曼ζ函数 ，
T
c
{\displaystyle T_{c}}
T c 玻色-爱因斯坦凝聚 开始形成时的临界温度。
^ Perrot, Pierre. A to Z of Thermodynamics. Oxford University Press. 1998. ISBN 0-19-856552-6 ^ van der Waals, J. D. On the Continuity of the Gaseous and Liquid States (doctoral dissertation). Universiteit Leiden. 1873. ^ Soave, G. Equilibrium Constants from a Modified Redlich-Kwong Equation of
State, Chem. Eng. Sci.,1 9 7 2, 27, 1197-1203
^ Peng, D. Y., and Robinson, D. B. A New Two-Constant Equation of State. Industrial and Engineering Chemistry: Fundamentals. 1976, 15 : 59–64. doi:10.1021/i160057a011 ^ 5.0 5.1 J. Richard Jr. Elliott, S. Jayaraman Suresh, Marc D. Donohue. A Simple Equation of State for Nonspherical and Associating Molecules. Ind. Eng. Chem. Res. 1990, 29 (7): 1476–1485. doi:10.1021/ie00103a057 ^ K.E. Starling. Fluid Properties for Light Petroleum Systems. Gulf Publishing Company . 1973.
