感谢网友yiqieke整理此文

背景：

钙钛矿混合导体氧化物同时具有氧离子导电性和电子导电性, 因此已被广泛地用作固体氧化物燃料电池的阴极和透氧膜。这类材料的性能主要由氧在其表面的交换速率以及氧在其体内的扩散速率决定。如果可以提高这两个速率, 固体氧化物燃料电池和透氧膜的操作温度将会降低, 而在较低温度下工作的燃料电池其材料的老化程度将会有所降低, 并且寿命会更长。因此, 如何提高这两个速率就成为了很多燃料电池和透氧膜材料研究者所感兴趣的问题。为此, 对不同的材料我们需要测量其氧表面交换系数和体相扩散系数。

目前测定氧表面交换系数和体相扩散系数的方法主要有同位素交换法和电导弛豫法（ECR,Electrical Conductivity Relaxation）。由于电导弛豫法所需仪器相对简单, 故应用广泛。

原理：

由于钙钛矿混合导体氧化物中存在多价态的过渡金属离子（Fe、Co等）, 在一定温度下突然改变氧化物所处环境的氧分压会造成过渡金属离子的变价, 同时氧化物中的氧（离子）浓度也发生变化。若在初始时刻（t=0t=0）时氧的平衡浓度是C0C0,当外界氧分压发生突变后, 氧在t时刻的浓度为CtCt。C0C0到CtCt的变化是通过氧扩散进出氧化物来实现的, 这种扩散会引起体系电导率sigmasigma的变化。假定体系电导率和氧浓度之间存在线性关系, 则

C0−CinftyoverCt−Cinfty=sigma0−sigmainftyoversigmat−sigmainftyC0−CinftyoverCt−Cinfty=sigma0−sigmainftyoversigmat−sigmainfty

其中, sigma0sigma0为初值时刻（t=0t=0）电导率, sigmatsigmat为氧分压突变后t时刻的电导率, sigmainftysigmainfty为氧浓度重新达到平衡稳定后电导率。

为求得氧的表面交换系数和体相扩散系数, 需要对扩散过程进行求解。

对扩散过程, Fick第一定律指出扩散通量与浓度梯度成正比, 其比例系数为扩散系数：

mathbfJ=−DnablaCmathbfJ=−DnablaC

由质量守恒定律,

partialCoverpartialt+nablacdotmathbfJ=0partialCoverpartialt+nablacdotmathbfJ=0

进而得到Fick第二定律

partialCoverpartialt−nablacdotDnablaC=0partialCoverpartialt−nablacdotDnablaC=0

若扩散系数D为常数, 与浓度无关, 则

partialCoverpartialt=Dnabla2CpartialCoverpartialt=Dnabla2C

对于最简单的一维情况

partialCoverpartialt=Dpartial2Coverpartial2xpartialCoverpartialt=Dpartial2Coverpartial2x

当xin[−a,a]xin[−a,a]时, 其边界条件为

begin{split}J(a)&=-D{partial C over partial x}|_{x=a}=K[C(a)-C(infty)] \\ J(-a)&=-D{partial C over partial x}|_{x=-a}=-K[C(-a)-C(infty)] end{split}begin{split}J(a)&=-D{partial C over partial x}|_{x=a}=K[C(a)-C(infty)] \\ J(-a)&=-D{partial C over partial x}|_{x=-a}=-K[C(-a)-C(infty)] end{split}

其中K为交换系数。设初始浓度C(x,0)=C0C(x,0)=C0,利用本征函数展开, 可将方程的解可写为

C(x,t)−C0overC(infty)−C0=1−sumlimitsinftyi=12Lalphacos(alphaix/a)over(alpha2i+Lalpha2+Lalpha)cosalphaiexp(−tovertaui)C(x,t)−C0overC(infty)−C0=1−sumlimitsi=1infty2Lalphacos(alphaix/a)over(alphai2+Lalpha2+Lalpha)cosalphaiexp(−tovertaui)

时间常数taui=a2overDalpha2itaui=a2overDalphai2, 变量alphaalpha满足方程alphaitan(alphai)=Lalpha=aKoverDalphaitan(alphai)=Lalpha=aKoverD

实际测量时得到的是浓度的平均值, 所以需要求出整个区间上浓度的平均值,

begin{split}{bar C(x,t)-C_0 over C(infty)-C_0} &= 1-{1 over2a}intlimits_{-a}^{a}sumlimits_{i=1}^{infty} {2L_alpha cos(alpha_i x/a)over (alpha_i^2+L_alpha^2+L_alpha)cosalpha_i}exp(-{t over tau_i})dx \\&=1-{1 over 2a}sumlimits_{i=1}^{infty} {2L_alpha over (alpha_i^2+L_alpha^2+L_alpha)cosalpha_i}exp(-{tover tau_i})\\ & quad quad intlimits_{-a}^{a} cos(alpha_i x/a) dx \\ &=1-sumlimits_{i=1}^{infty}{2L_alpha^2 over alpha_i^2 (alpha_i^2+L_alpha^2+L_alpha)}exp(-{t overtau_i}) end{split}begin{split}{bar C(x,t)-C_0 over C(infty)-C_0} &= 1-{1 over2a}intlimits_{-a}^{a}sumlimits_{i=1}^{infty} {2L_alpha cos(alpha_i x/a)over (alpha_i^2+L_alpha^2+L_alpha)cosalpha_i}exp(-{t over tau_i})dx \\&=1-{1 over 2a}sumlimits_{i=1}^{infty} {2L_alpha over (alpha_i^2+L_alpha^2+L_alpha)cosalpha_i}exp(-{tover tau_i})\\ & quad quad intlimits_{-a}^{a} cos(alpha_i x/a) dx \\ &=1-sumlimits_{i=1}^{infty}{2L_alpha^2 over alpha_i^2 (alpha_i^2+L_alpha^2+L_alpha)}exp(-{t overtau_i}) end{split}

上式有两种极限情形：

1. KllD/a,Lalphato0KllD/a,Lalphato0, 此时tauitaui增长很快, 只取第一项, alpha1tanalpha1simeqalpha21=Lalphaalpha1tanalpha1simeqalpha12=Lalpha, 故

begin{split}{bar C(x,t)-C_0 over C(infty)-C_0} & = 1-{2L_alpha^2 overalpha_1^2(alpha_1^2+L_alpha^2+L_alpha)}exp(-{t over tau_1}) \\ &=1-exp(-{Ktover a}) end{split}begin{split}{bar C(x,t)-C_0 over C(infty)-C_0} & = 1-{2L_alpha^2 overalpha_1^2(alpha_1^2+L_alpha^2+L_alpha)}exp(-{t over tau_1}) \\ &=1-exp(-{Ktover a}) end{split}

2. KggD/a,LalphatoinftyKggD/a,Lalphatoinfty, 此时alphai=(i+1over2)pialphai=(i+1over2)pi,

barC(x,t)−C0overC(infty)−C0=1−8overpi2sumlimitsinftyi=01over(2i+1)2exp(−(2i+1)2pi2Dtover4a2)barC(x,t)−C0overC(infty)−C0=1−8overpi2sumlimitsi=0infty1over(2i+1)2exp(−(2i+1)2pi2Dtover4a2)

按这两种极限情况进行数据处理都无法得到氧交换系数K, 还需借助其他的关系式.

上面的方程可以推广到三维, 设三维区间长宽高分别为2a,2b, 2c, 则

barC(x,t)−C0overC(infty)−C0=1−sumlimitsinftyi=1sumlimitsinftyj=1sumlimitsinftyk=1Aijkexp(−tovertauijk)barC(x,t)−C0overC(infty)−C0=1−sumlimitsi=1inftysumlimitsj=1inftysumlimitsk=1inftyAijkexp(−tovertauijk)

begin{split}A_{ijk}&={8 L_alpha^2 L_beta^2 L_gamma^2 over alpha_i^2 beta_j^2gamma_k^2 (alpha_i^2+L_alpha^2+L_alpha) (beta_j^2+L_beta^2+L_beta)(gamma_i^2+L_gamma^2+L_gamma) } \\ {1 over tau_{ijk}} &= {1 overtau_i }+{1 over tau_j }+{1 over tau_k } end{split}begin{split}A_{ijk}&={8 L_alpha^2 L_beta^2 L_gamma^2 over alpha_i^2 beta_j^2gamma_k^2 (alpha_i^2+L_alpha^2+L_alpha) (beta_j^2+L_beta^2+L_beta)(gamma_i^2+L_gamma^2+L_gamma) } \\ {1 over tau_{ijk}} &= {1 overtau_i }+{1 over tau_j }+{1 over tau_k } end{split}

taui=a2overDalpha2i,tauj=b2overDbeta2j,tauk=c2overDgamma2ktaui=a2overDalphai2,tauj=b2overDbetaj2,tauk=c2overDgammak2

alphaitan(alphai)=Lalpha=aKoverD,betajtan(betaj)=Lbeta=bKoverD,gammaktan(gammak)=Lgamma=cKoverDalphaitan(alphai)=Lalpha=aKoverD,betajtan(betaj)=Lbeta=bKoverD,gammaktan(gammak)=Lgamma=cKoverD

实验：

测量条状测试样品的尺寸。

用四探针电导法测定长方体样品的电导率sigma0sigma0,然后瞬间改变样品所在气氛的氧分压（比如从0.21 atm变到0.01atm）, 并记录固定时间间隔（如10s）处的电导率值sigmatsigmat。当重新达到平衡后, 记录其电导率值sigmainftysigmainfty。

实验结果与数据处理：

利用所测sigma0,sigmainfty,(t,sigmat)sigma0,sigmainfty,(t,sigmat)计算S(t)=sigma0−sigmainftyoversigmat−sigmainftyS(t)=sigma0−sigmainftyoversigmat−sigmainfty, 对S(t)S(t)进行最小二乘拟合, 便可求出氧表面交换系数K和体相扩散系数D。

这是一个典型的非线性拟合问题, 利用MatLab编写代码时须注意alpha,beta,gammaalpha,beta,gamma需要在每一步拟合时进行求解, 这可借助于fsolve函数实现。

下面是某次实验所得数据及拟合结果。

参考代码：

参考：

王严东, 吕喆, 魏波. 无机材料学报,25(6):635-0, 2010 C.-R. Song, H.-I. Yoo. Solid State Ionics,124:2-299, 1999 

R.A. Cox-Galhotra, S. McIntosh. Solid State Ionics,181:1429-1436, 2010 

M.W. den Otter, A study of OXYGEN TRANSPORT in mixed conducting oxides usingisotopic exchange and conductivity relaxation下载本文