一、PEMFC基础之热力学
- 1.内能U、焓H、熵S、吉布斯自由能G、赫姆霍兹自由能F关系图
- 2.可逆电压与温度和压力的关系
- 3.能斯特方程
- 4.燃料电池效率
1.内能U、焓H、熵S、吉布斯自由能G、赫姆霍兹自由能F关系图
2.可逆电压与温度和压力的关系
标准状态可逆电压Er计算:
E
r
=
−
△
G
n
F
E_{r}=-\frac{\bigtriangleup G}{nF}
Er=−nF△G
与温度的关系:
E
T
=
E
r
+
△
S
n
F
(
T
−
298
)
E_{T}=E_{r}+\frac{\triangle S}{nF}\left ( T-298 \right )
ET=Er+nF△S(T−298)
燃料电池中∆S为负值,因此可逆电压随温度升高而下降。
与压力的关系:
E
P
=
E
r
−
△
V
n
F
(
P
−
P
0
)
E_{P}=E_{r}-\frac{\triangle V}{nF}\left ( P-P_{0} \right )
EP=Er−nF△V(P−P0)
根据勒沙特列原理,燃料电池中压力增大∆V为负值,因此可逆电压随着压力增大而增大。
3.能斯特方程
E = E r − R T l n ( P H 2 P O 2 0.5 P H 2 O ) = E r − R T n F l n α H 2 O α H 2 α O 2 0.5 E=E_{r}-RT\frac{}{}ln\left ( \frac{P_{H2}P_{O2}^{0.5}}{P_{H2O}} \right )=E_{r}-\frac{RT}{nF}ln\frac{\alpha _{H2O}}{\alpha _{H2}\alpha _{O2}^{0.5}} E=Er−RTln(PH2OPH2PO20.5)=Er−nFRTlnαH2αO20.5αH2O
用压力计算时:
饱和水蒸气压力(Bar)计算公式:
l
g
P
S
A
T
=
−
2.1794
+
0.02953
T
−
9.1837
×
1
0
−
5
T
2
+
1.4454
×
1
0
−
7
T
3
lgP_{SAT}=-2.1794+0.02953T-9.1837\times 10^{-5}T^{2}+1.4454\times 10^{-7}T^{3}
lgPSAT=−2.1794+0.02953T−9.1837×10−5T2+1.4454×10−7T3
相对湿度RH非100%时,水蒸气压力:
P
H
2
O
=
P
S
A
T
⋅
R
H
P_H2O=P_{SAT}\cdot RH
PH2O=PSAT⋅RH
阳极氢气压力计算公式:
P
H
2
=
0.5
×
P
H
2
,
i
n
e
x
p
(
1.653
i
T
K
1.334
)
−
P
H
2
O
P_H2=0.5\times \frac{P_{H2,in}}{exp\left ( \frac{1.653i}{T_{K}^{1.334}} \right )}-P_{H2O}
PH2=0.5×exp(TK1.3341.653i)PH2,in−PH2O
阴极氧气压力计算公式:
P
O
2
=
P
a
i
r
e
x
p
(
4.192
i
T
K
1.334
)
−
P
H
2
O
P_O2=\frac{P_{air}}{exp\left ( \frac{4.192i}{T_{K}^{1.334}} \right )}-P_{H2O}
PO2=exp(TK1.3344.192i)Pair−PH2O
注意:
对于许多实际的目的,你可以假设在使用PEM燃料电池时,水的活性(a)和相对湿度(RH)大致相等。这是因为水在膜和催化剂层中的平衡蒸汽压通常接近于纯水的平衡蒸汽压。
用活度α计算时:
对于理想气体:根据乌拉尔定律,活度α为气体分压/标准态饱和压力。
1个大气压下,空气中氧气活度大约为0.21;2个大气压下活度约为0.42。
对于非理想气体:路易斯引入了活度系数γ来描述与理想状态的偏离。
对于纯组分物质,例如铂、氢气、液态水等,活度α为1。
注意:
a. 在阴极n=4;阳极n=2。
b.尽管在能斯特方程中温度为变量,但并没有完全解释可逆电压随温度的变化关系,如果T≠T0,需要同时使用上述与温度的关系来修正可逆电压。
4.燃料电池效率
1.理想的可逆效率计算公式:
η
f
u
e
l
c
e
l
l
=
△
G
△
H
\eta _{fuel_cell}=\frac{\bigtriangleup G}{\triangle H}
ηfuelcell=△H△G
标准状态下可逆效率为83%
2.实际的可逆效率计算公式:
η
r
e
a
l
=
η
t
h
e
r
a
m
×
η
v
o
l
t
a
g
e
×
η
f
u
e
l
=
(
△
G
△
H
)
(
V
E
)
(
1
λ
)
\eta _{real}=\eta _{theram}\times \eta _{voltage}\times \eta _{fuel}=\left ( \frac{\bigtriangleup G}{\bigtriangleup H} \right )\left ( \frac{V}{E} \right )\left ( \frac{1}{\lambda } \right )
ηreal=ηtheram×ηvoltage×ηfuel=(△H△G)(EV)(λ1)