电导弛豫法测材料的氧表面交换系数和体相扩散系数Word文件下载.docx

1、对不同的材料我们需要测量其氧表面交换系数和体相扩散系数。目前测定氧表面交换系数和体相扩散系数的方法主要有同位素交换法和电导弛豫法（ECR,Electrical Conductivity Relaxation）。由于电导弛豫法所需仪器相对简单,故应用广泛。原理：由于钙钛矿混合导体氧化物中存在多价态的过渡金属离子（Fe、Co等）,在一定温度下突然改变氧化物所处环境的氧分压会造成过渡金属离子的变价,同时氧化物中的氧（离子）浓度也发生变化。若在初始时刻（t=0t=0）时氧的平衡浓度是C0C0,当外界氧分压发生突变后,氧在t时刻的浓度为CtCt。C0C0到CtCt的变化是通过氧扩散进出氧化物来实现的,这

2、种扩散会引起体系电导率sigmasigma的变化。假定体系电导率和氧浓度之间存在线性关系,则C0CinftyoverCtCinfty=sigma0sigmainftyoversigmatsigmainftyC0CinftyoverCtCinfty=sigma0sigmainftyoversigmatsigmainfty其中,sigma0sigma0为初值时刻（t=0t=0）电导率,sigmatsigmat为氧分压突变后t时刻的电导率,sigmainftysigmainfty为氧浓度重新达到平衡稳定后电导率。为求得氧的表面交换系数和体相扩散系数,需要对扩散过程进行求解。对扩散过程, Fick第一

3、定律指出扩散通量与浓度梯度成正比,其比例系数为扩散系数：mathbfJ=DnablaCmathbfJ=DnablaC由质量守恒定律,partialCoverpartialt+nablacdotmathbfJ=0partialCoverpartialt+nablacdotmathbfJ=0进而得到Fick第二定律partialCoverpartialtnablacdotDnablaC=0partialCoverpartialtnablacdotDnablaC=0若扩散系数D为常数,与浓度无关,partialCoverpartialt=Dnabla2CpartialCoverpartialt=Dn

4、abla2C对于最简单的一维情况partialCoverpartialt=Dpartial2Coverpartial2xpartialCoverpartialt=Dpartial2Coverpartial2x当xina,axina,a时,其边界条件为beginsplitJ（a）&=-DpartialCoverpartialx|_x=a=KC（a）-C（infty）J（-a）&x|_x=-a=-KC（-a）-C（infty）endsplitbeginsplitJ（a）&endsplit其中K为交换系数。设初始浓度C（x,0）=C0C（x,0）=C0,利用本征函数展开,可将方程的解可写为C（x,t

5、）C0overC（infty）C0=1sumlimitsinftyi=12Lalphacos（alphaix/a）over（alpha2i+Lalpha2+Lalpha）cosalphaiexp（tovertaui）C（x,t）C0overC（infty）C0=1sumlimitsi=1infty2Lalphacos（alphaix/a）over（alphai2+Lalpha2+Lalpha）cosalphaiexp（tovertaui）时间常数taui=a2overDalpha2itaui=a2overDalphai2,变量alphaalpha满足方程alphaitan（alphai）=La

6、lpha=aKoverDalphaitan（alphai）=Lalpha=aKoverD实际测量时得到的是浓度的平均值,所以需要求出整个区间上浓度的平均值,beginsplitbarC（x,t）-C_0C（infty）-C_0&=1-1over2aintlimits_-aasumlimits_i=1infty2L_alphacos（alpha_ix/a）over（alpha_i2+L_alpha2+L_alpha）cosalpha_iexp（-ttau_i）dx&=1-12asumlimits_i=1infty（alpha_i2+L_alpha2+L_alpha）cosalpha_iexp（-

7、tovertau_i）quadintlimits_-aax/a）dx=1-sumlimits_i=1infty2L_alpha2alpha_i2（alpha_i2+L_alpha2+L_alpha）exp（-tovertau_i）endsplitbeginsplitbar上式有两种极限情形：1.KllD/a,Lalphato0KllD/a,Lalphato0,此时tauitaui增长很快,只取第一项,alpha1tanalpha1simeqalpha21=Lalphaalpha1tanalpha1simeqalpha12=Lalpha,故1-2L_alpha2overalpha_12（alph

8、a_12+L_alpha2+L_alpha）exp（-ttau_1）=1-exp（-Ktovera）2.KggD/a,LalphatoinftyKggD/a,Lalphatoinfty,此时alphai=（i+1over2）pialphai=（i+1over2）pi,barC（x,t）C0overC（infty）C0=18overpi2sumlimitsinftyi=01over（2i+1）2exp（2i+1）2pi2Dtover4a2）barC（x,t）C0overC（infty）C0=18overpi2sumlimitsi=0infty1over（2i+1）2exp（2i+1）2pi2Dt

9、over4a2）按这两种极限情况进行数据处理都无法得到氧交换系数K,还需借助其他的关系式.上面的方程可以推广到三维,设三维区间长宽高分别为2a,2b, 2c,barC（x,t）C0overC（infty）C0=1sumlimitsinftyi=1sumlimitsinftyj=1sumlimitsinftyk=1Aijkexp（tovertauijk）barC（x,t）C0overC（infty）C0=1sumlimitsi=1inftysumlimitsj=1inftysumlimitsk=1inftyAijkexp（tovertauijk）beginsplitA_ijk&=8L_alpha

10、2L_beta2L_gamma2beta_j2gamma_k2（alpha_i2+L_alpha2+L_alpha）（beta_j2+L_beta2+L_beta）（gamma_i2+L_gamma2+L_gamma）1tau_ijkovertau_i+1tau_jtau_kendsplitbeginsplitA_ijk&taui=a2overDalpha2i,tauj=b2overDbeta2j,tauk=c2overDgamma2ktaui=a2overDalphai2,tauj=b2overDbetaj2,tauk=c2overDgammak2alphaitan（alphai）=Lalpha=aKoverD,betajtan（betaj）=Lbeta=bKoverD,gammaktan（gammak）=Lgamma=cKoverDalphaitan（alphai）=Lalpha=aKoverD,betajtan（betaj）=Lbeta=bKoverD,gammaktan（gammak）=Lgamma=cKoverD实验：测量条状测试样品的尺寸。用四探针电导法测定长方体样品的电导率sigma0sigma0,然后瞬间改变样品所在气氛的