新型高效热离子功率器件的性能特性研究*
2019-10-09廖天军林比宏王宇珲
廖天军 林比宏 王宇珲
1) (重庆理工大学物理与能源系,重庆 400054)
2) (华侨大学物理系,厦门 361021)
1 引 言
开发热电能源转换器件是解决能源紧缺和环境污染的重要途径[1,2].热离子功率器件(TPD)是利用材料表面热电子发射将热能转化为电能的一种装置,它具有体积小、功率密度大、热电转换效率高等特点,在民用、军事、航空航天等领域具有重要应用,它可利用诸如汽车尾气、工业余热、矿物燃料、太阳能、核能等多种形式的热源驱动,实现热电转换,其中以核燃料为热源的TPD已应用于星际考察等空间技术.TPD主要由阴极和阳极组成,两个电极被真空间隙隔开[2−7].阴极经外部高温热源加热后温度升高,其内部的电子获得能量,一部分电子可以克服金属表面“势垒”的障碍,摆脱金属原子核的束缚,逸出金属表面,再通过真空间隙聚集在阳极,电子通过外部负载回到阴极,构成回路,从而把加热阴极的部分热能转变成负载上消耗的电能[4,5].TPD的热电转换效率的上限由卡诺效率所决定,阴极和阳极的温度相当于进出口温度[8].基于热离子发射理论,近年来有人提出聚光TPD太阳能电池[9−13]、聚光石墨烯TPD太阳能电池的物理模型[14]和聚光多层石墨烯TPD太阳能电池的物理模型[15].
石墨烯具有奇特的光学、电学、力学特性,在材料学、微纳加工和能源存储与转换等方面具有广阔应用前景.石墨烯在高温下工作的TPD的热离子发射性能优于金属材料,因此,石墨烯更适合作为TPD的阴极材料.Liang和Ang[16]提出了新型高效石墨烯TPD的物理模型,给出了器件的工作机理,研究了单层石墨烯的热离子发射现象,导出功率密度和效率的解析表达式,获得了许多有利于热能开发与利用的新结论,为石墨烯TPD的发展奠定了理论基础.Liang等[17]还设计了一类采用范德瓦耳斯异质结构夹在两个石墨烯电极之间的固态TPD,研究发现高温热源的温度可在400—500 K范围内实现高能量转换效率.Mishra等[18]提出了一种单层石墨烯热电子发射(忽略衬底效应)的物理模型,并探讨了其作为TPD的阴极的应用,且推导出关于单层石墨烯热离子激发电流密度的解析式,研究发现热离子激发电流是温度、功函数和费米能级等参数的函数,该公式不同于传统的Richardson-Dushman方程.Mishra等[19]还研究了有限温度下多层石墨烯热电子发射的理论模型,并论证了其作为阴极的热离子能量转换方案的可行性.
尽管国内外学者深入研究了石墨烯TPD的工作机理和性能特性,而对不可逆因素对器件优化运行的影响和参数的优化设计方面的研究还相对较少.本文应用固体物理和不可逆热力学理论,研究TPD的性能特性,推导出TPD的功率密度和效率所满足的表达式,研究TPD的伏安特性以及阴极板温度与输出电压的依赖关系,综合功率密度和效率随输出电压与集电极功函数变化的情况,给出参数的优化设计策略.探讨外部热源的温度对器件优化运行的影响.
2 石墨烯TPD的功率密度和效率
石墨烯TPD的结构示意图见图1,其中石墨烯材料覆盖在金属钨表面构成阴极板,而阳极由金属钨材料构成,TC和TA分别为两个极板的温度,TH和TE分别为高温热源和环境的温度,QH为高温热源与阴极板之间的热传递量,QL为阳极与环境之间的热传递量.石墨烯TPD的工作机理是高温热源通过导热材料与阴极接触,使其获得高温,极板内部的一部分电子获得能量而逸出表面,然后通过真空间隙被阳极收集,产生热离子激发电流密度JC,阳极板吸收来自阴极的热电子和光子,导致其温度升高,产生反向热离子激发电流密度JA.由于JC>JA,净电子流将通过外部负载RL回到阴极,构成电子循环回路,实现热电能量转换.值得注意的是金属钨具有较高功函数(约4.5 eV),应用目前的材料工艺技术,在阳极板表面喷涂负电子亲和势材料铯原子,可将其功函数降至1.5 eV,从而提高TPD输出电压[4].采用以上所述的相同工艺也可降低石墨烯功函数,从而提高阴极板热离子激发电流的能力[13].
图1 石墨烯TPD示意图Fig.1.Schematic diagram of a graphene-based TPD.
在图1的模型中,TPD的热电转换效率和功率密度可分别表示为[16,17]
和
式中V=(ΦC−ΦA)/q是TPD的输出电压,ΦC和ΦA分别是阴极和阳极的功函数,q是一个单位正电荷的电量,J=JC−JA是净电流密度.阴极表面激发的电流密度JC 可表为[16,17]
其中,EF是石墨烯费米能级0.01158 A⋅cm–2⋅K–3是修正Richardson-Dushman(RD)常数,=h/(2π) 是约化普朗克常数,kB是玻尔兹曼常数,υF≈108cm⋅s–1为无质量的费米子的费米速度.
阳极的热离子激发电流密度JA由传统RD方程决定[20−22],即:
根据图1和热力学第一定律,两个极板所满足的能量平衡方程可分别表示为[17,20]
和
假设TPD高温端和低温端的热传递服从牛顿冷却定律,QH和QL表示为[22]
和
其中,UH和UL分别为高低温端的热导系数.给定相关参数如ε,TH,TE等,联立(2)—(8)式可数值求解出两极板的温度TC和TA,将TC和TA代入(1)—(4)式,可得石墨烯TPD的电流密度、功率密度和效率.
3 TPD的性能特性与优化设计
当给定三个不同阴极功函数时,阴极和阳极的温度随输出电压变化的关系曲线及伏安特性见图2,其中,ε=0.20 ,UH=UL=0.2 W⋅cm–2⋅K–1,TH=1500 K和TE=300 K.图2(a)显示阴极板的工作温度随电压的增加几乎保持不变,而阴极功函数越大则阴极工作温度越高,这是由于功函数ΦC越大,对电子的束缚能力越强,只有在较高的温度下,电子才能从极板表面逸出.
图2(b)显示阳极温度随着电压的升高而降低,当ΦC=2.0 eV时,阳极板温度的变化最为明显,从图2(b)还可以看出ΦC越大,阳极温度则越低.结合图2(a)和图2(b),可看出ΦC越大,两个极板的温差TC−TA则越高.从图2(c)可以看出TPD的伏安特性与传统太阳能电池的伏安特性曲线相似,随着输出电压的增加,电流密度先保持不变,然后随着电压的进一步升高而迅速降低.不难发现阴极功函数越高,回路中的电流密度越小,这说明降低极板功函数的重要性.除了上面所述降低极板功函数的方法外,具有负电子亲和势特征的典型材料还有Ⅲ—Ⅴ族化合物半导体材料砷化镓、氮化镓及其多元化合物铝镓砷、铟镓氮也可用于降低极板功函数.ΦC的大小为真空能级EV与费米能级EF之差,即:ΦC=EV−EF,由此可见调节石墨烯的费米能级也可降低ΦC.相关研究表明可通过化学掺杂和外加门电压等方法调节石墨烯的费米能级[16].从图2(c)还可以看出TPD的电压较小而电流密度较大,从而只能驱动较小电压的负载,因此,在实际设计TPD系统时,TPD的电极面积要尽可能小,同时需要将多个TPD串联以提高输出电压和带负载能力.当电路短路时,V=0V ,可得到TPD的短路电流密度Jsc; 当电路开路时,J=0 A⋅cm–2,可得到TPD的开路电压Voc.在伏安特性曲线中可找到一对优化值Vopt和Jopt使功率密度达到其优化值Popt,根据4个重要参数Jsc,Voc,Vopt和Jopt可确定一个表征TPD性能高低的填充因子FF,即:
图2 给定不同 ΦC 时,(a)阴极温度和(b)阳极温度随输出电压变化曲线及(c)伏安特性曲线,其中 ε=0.20 ,EF=0.80eV ,UH=UL=0.2W·cm−2·K−1,TH=1500K 和TE=300KFig.2.(a) Curves of the cathode temperature,and (b) the anode temperature varying with the output voltage,and (c) the voltampere characteristic for given values of ΦC ,where ε=0.20 ,EF=0.80eV ,UH=UL=0.2W·cm−2·K−1 ,TH=1500K ,and TE=300K.
FF的值越接近1则表示TPD的性能越好.
利用(1)—(6)式和数值计算,可得TPD的功率密度和效率随阴极功函数和输出电压变化的等高线图,如图3所示,其中相关参数的取值与图2相同.给定电压V时,根据V=(ΦC−ΦA)/q及(3)和(4)式可知,当ΦC增大时,ΦA随ΦC线性增大,导致热离子激发电流密度JC和JA减小,而两个极板的温差TC−TA随ΦC的增加而增加,从而导致净电流密度J=JC−JA随V的变化存在极值.另外,图2(c)已证明给定ΦC时,存在优化功率密度Popt和相应的优化值Vopt和Jopt,进而可以同时优化阴极功函数和电压得到最大功率密度Pmax和效率ηmax.尽管如此,图3显示在不同的阴极功函数和输出电压取得最大功率密度和效率,在取ηmax时的阴极功函数和输出电压均大于取Pmax时的阴极功函数和电压,这主要是由TPD内外部的不可逆热损失所导致的.
为了进一步确定TPD参数的优化区间,给定阴极功函数,通过优化电压V,可得到优化功率密度Popt和效率ηopt,如图4(a)和图4(b)所示,其中相关参数的取值与图2相同.图4(a)和图4(b)显示优化功率密度和效率随电压的增加先增加后减小,当ΦC=ΦC,P和ΦC=ΦC,η时,功率密度和效率分别达到最大值Pmax和ηmax,并且ΦC,P<ΦC,η.根据ΦC,P和ΦC,η可确定TPD处于Pmax和ηmax 时电压的最佳值为VP和Vη.同时,图4也显示了优化功率密度和效率时功函数的取值均是电压的单调递增函数.当V
表1为本文和文献[23]在优化性能时重要参数取值的对比,其中相关参数取值相等.与文献[23]相比,本文首次同时考虑了石墨烯TPD高温和低温两端与热源的不可逆传热对两个极板温度的影响,对TPD的输出电压和两个极板功函数等参数进行了优化分析,获得了它们的优化值.因此,本文对石墨烯TPD性能的优化方法更加普适.从表1可看出本文所研究石墨烯TPD的性能得到显著提高,可以获得比文献[23]更高的效率,最大效率可达到60%.主要原因是文献[23]给定阴极板功函数ΦC=3.0 eV,对输出电压优化得到Pmax和ηmax,而本文是同时优化电压和阴极板功函数,得到Pmax和ηmax.为了与文献[23]所获得的结果相比较,当本文也给定ΦC=3.0 eV时,Pmax和ηmax分别为2.16 A⋅cm–2和 0.264 ,可见所获得的性能优于文献[23]的Pmax和ηmax,其主要原因是本文根据石墨烯RD方程,修正了阴极和阳极热流表达式,而文献[23]未考虑这方面因素.表1显示本文阳极功函数的优化区间为0.59—0.79 eV,在该区间TPD可获得理论上的最优性能,其最大效率是文献[23]所给出的最大效率的两倍.必须指出的是,在此区间的阳极功函数与阳极功函数最低值1.50 eV具有较大差异,而文献[23]阳极功函数的优化区间为1.75—2.21 eV,比较接近目前实验上可实现的阳极功函数的数值,所以文献[23]的模型在实验上相对更容易实现.然而,文献[23]对石墨烯TPD性能的优化只是局域优化,没有获得装置的最优理论性能.本文的优化方法可应用于文献[23]的模型的优化.随着科学技术的发展,通过新技术还可以使阳极功函数进一步降低,从而在实验上实现本文所提出的最优模型.另外,最新研究表明在TPD内部添加钡和铯蒸汽,通过改变蒸汽压强,可调节极板功函数[24].
图3 (a) TPD的功率密度和(b)效率随输出电压和阴极板功函数变化的三维图Fig.3.Three-dimensional graphs of (a) the power density and (b) the efficiency varying with the output voltage and the work function of the cathode.
图4 (a) TPD的优化功率密度和阴极功函数,(b)优化效率和阴极功函数随电压变化的曲线以及(c)性能特征曲线ηopt-PoptFig.4.Curves of (a) the optimal power density and work function,(b) the optimal efficiency and work function varying with the voltage,and (c) the performance characteristic of TPD.
表1 本文和文献[23]在优化性能时重要参数的取值.Table 1.Values of key parameters at the optimum performances for the present work and the Ref.[23].
图5(a)显示Pmax随TH的升高呈线性增加趋势,而ηmax随TH的升高先增加后减小,存在优化温度TH,η使ηmax达到其最大值ηM.图5(b)和图5(c)显示VP,Vη,ΦC,P和ΦC,η是TH 的单调递增函数,而优化区间Vη-VP和ΦC,η-ΦC,P随TH的增加而基本保持不变.若选择TH 图6显示Pmax和ηmax是EF的单调递增函数,因此,增加EF可提升Pmax和ηmax.需要指出的是,目前石墨烯的EF可调范围为0.50—0.85 eV[16].数值模拟结果还表明VP,Vη,ΦC,P和ΦC,η也是EF 的单调递增函数,这说明增加EF不仅可以提高TPD的输出电压,而且较大的ΦC,P和ΦC,η同时也降低了TPD的设计难度. 图5 (a) Pmax和ηmax ,(b) VP和Vη ,(c) ΦC,P和ΦC,η 随TH的变化Fig.5.Curves of (a) Pmax and ηmax ,(b) VP and Vη ,and(c) ΦC,P and ΦC,η as a function of TH. 图6 Pmax和ηmax 随 EF 变化的曲线,其中TH=1500KFig.6.Curves of Pmax and ηmax as a function of EF ,where TH=1500K. 本文详细研究了TPD的性能特性与参数优化设计,结果表明通过优化电压和阴极板功函数可显著提高TPD的功率密度和效率,同时通过折衷考虑功率和效率,可确定TPD性能参数的优化区间.另外,本文也分析了高温热源温度对TPD优化性能特性的影响,确定了高温热源温度的优化范围.本文所得结果可为高性能热电功率器件的设计与研制提供理论指导.
4 结 论
