电导弛豫法测材料的氧表面交换系数和体相扩散系数及Matlab使用方法

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

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

3、氧浓度之间存在线性关系, 则：其中,0为初值时刻（t=0）电导率, t为氧分压突变后t时刻的电导率, 为氧浓度重新达到平衡稳定后电导率。为求得氧的表面交换系数和体相扩散系数, 需要对扩散过程进行求解。对扩散过程, Fick第一定律指出扩散通量与浓度梯度成正比, 其比例系数为扩散系数：由质量守恒定律：进而得到Fick第二定律：若扩散系数D为常数, 与浓度无关, 则：对于最简单的一维情况：当xa,axa,a时, 其边界条件为：其中K为交换系数。设初始浓度C(x,0)=C0C(x,0)=C0,利用本征函数展开, 

4、;可将方程的解可写为：时间常数i=a2D2ii=a2Di2, 变量满足方程itan(i)=L=aKDitan(i)=L=aKD实际测量时得到的是浓度的平均值, 所以需要求出整个区间上浓度的平均值上式有两种极限情形：1. KD/a,L0,KD/a,L0, 此时i增长很快, 只取第一项, 1tan121=L1tan112=L, 故2. KD/a,L,KD/a,L, 此i=(i+12)i=(i+12)按这两种极限情况进行数据处理都无法得到氧交换系数K, 还需借助其他的关系式.上面的方程可以推广到三维, 

5、;设三维区间长宽高分别为2a,2b, 2c, 则：实验：测量条状测试样品的尺寸。用四探针电导法测定长方体样品的电导率0,然后瞬间改变样品所在气氛的氧分压（比如从0.21 atm变到0.01atm）, 并记录固定时间间隔（如10s）处的电导率值t。当重新达到平衡后, 记录其电导率值。实验结果与数据处理：1、打开origin软件编辑1，2，3作为X轴，电阻作为Y轴。复制数据到Y轴。2、X：去掉前面的60个点，后每隔多少点取一个平均值则【colA-60】*平均值个数。Y：【1/colB】之后减去第一行的数据（取前小数点后两位），之后再除以最后一行的数据（同样取小数点后两位）。3、复制Y轴数据，打开excel中去平均值文件。改变算法后再黏贴数据回origin，之后复制XY轴的数据到txt文档D并保存。4、打开Matlab软件，在其中打开fit2文件夹，找到cond-calc.m（第二个）改掉其中的长宽高，之后关掉此窗口，如图1。再打开程序命令.txt复制下图蓝色部分，并关掉窗口，如图2。图1图25、在右边workspace中找到绿色向下箭头标志，打开之前第三步保存的数据txt文档D，如图1。会弹出图2窗口，点next后点finish如图3。图1图2图36、在command window中黏贴第四步图