-
According to the theories of the solid physics and irreversible thermodynamics, the performance characteristics of a novel high-efficient graphene thermionic power device (TPD) are studied. The temperature of the cathode plate and anode plate are determined by solving the energy balance equation of hot and cold sides of the TPD. The effects of the output voltage and the work function of the cathode on the volt-ampere characteristics of the TPD and the temperature of the two electrodes are analyzed to determine the parametric characteristics of the TPD at the maximum power density and efficiency. The power density and efficiency are compromised, and the parametric optimal designs are given. The influence of the temperature of heat source at high temperature on optimization performance is analyzed. The results obtained here can provide theoretical guidance for developing the thermionic energy conversion devices.
-
Keywords:
- thermionic power device /
- graphene /
- work function /
- performance optimization
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的热电转换效率和功率密度可分别表示为[16,17]
η=P/QH (1) 和
P=QH−QL=VJ, (2) 式中
V=(ΦC−ΦA)/q 是TPD的输出电压,ΦC 和ΦA 分别是阴极和阳极的功函数, q是一个单位正电荷的电量,J=JC−JA 是净电流密度. 阴极表面激发的电流密度JC 可表为[16,17]JC=A∗T3Cexp[−(ΦC−EF)/(kBTC)], (3) 其中,
EF 是石墨烯费米能级,A∗=qk3B/(πℏ3υ2F)= 0.01158 A⋅cm–2⋅K–3是修正Richardson-Dushman (RD)常数,ℏ=h/(2π) 是约化普朗克常数,kB 是玻尔兹曼常数,υF≈108 cm⋅s–1为无质量的费米子的费米速度.阳极的热离子激发电流密度
JA 由传统RD方程决定[20—22], 即:JA=AT2Aexp[−ΦA/(kBTA)], (4) 其中,
A=4πqmk3B/h3=120 A⋅cm–2⋅K–2是一般RD常数, m为静止电子质量.根据图1和热力学第一定律, 两个极板所满足的能量平衡方程可分别表示为[17,20]
QH=[JC(ΦC+3kBTC)−JA(ΦC+2kBTA)]/q+εσ(T4C−T4A) (5) 和
QL=[JC(ΦA+3kBTC)−JA(ΦA+2kBTA)]/q+εσ(T4C−T4A), (6) 其中,
εσ(T4C−T4A) 为两个极板之间的辐射热损失,ε=1/(ε−1C+ε−1A−1) 两个极板之间的等效热发射率,εC 和εA 分别是阴极和阳极表面的热发射率,σ 为斯特潘-玻尔兹曼常数;(ΦC+3kBTC) 是阴极表面逸出的电子所携带的热流率,(ΦC+2kBTA) 是阳极表面逸出的电子到达阴极所携带的热流率,3kBTC 和2kBTA 分别是阴极和阳极中电子越过势垒后的平均热动能;(ΦA+2kBTA) 是阳极表面逸出的每个电子所携带的热流率,(ΦA+3kBTC) 是离开阴极的电子到达阳极所携带的热流率.假设TPD高温端和低温端的热传递服从牛顿冷却定律,
QH 和QL 表示为[22]QH=UH(TH−TC) (7) 和
QL=UL(TA−TE), (8) 其中,
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 给定不同 时, (a)阴极温度和(b)阳极温度随输出电压变化曲线及(c)伏安特性曲线, 其中ΦC ,ε=0.20 ,EF=0.80eV ,UH=UL=0.2W⋅cm−2⋅K−1 和TH=1500K TE=300K Fig. 2. (a) Curves of the cathode temperature, and (b) the anode temperature varying with the output voltage, and (c) the volt-ampere characteristic for given values of , whereΦC ,ε=0.20 ,EF=0.80eV ,UH=UL=0.2W⋅cm−2⋅K−1 , andTH=1500K .TE=300K 图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 , 即:FF=JoptVoptJscVoc=PoptJscVoc, (9) 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<VP 时,Popt 和ηopt 随V的减小而降低; 而当V>Vη 时,Popt 和ηopt 随V的增加而减小. 因此,VP⩽V⩽Vη 和ΦC,P⩽ΦC⩽ΦC,η 为TPD的优化区间. 当电压和功函数处于该优化区间时, 图4(c)性能特征曲线中负斜率部分是功率密度和效率的优化区间, 即:Pη⩽Popt⩽Pmax 和ηP⩽η⩽ηmax , 其中Pη 和ηP 分别是η=ηmax 时的功率密度和P=Pmax 时的效率.表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].图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<TH,η 作为高温热源的优化运行区间, 在该区间内增加高温热源温度TH 可显著提高Pmax , 但过高的温度将会使TPD阴极板的稳定性和使用寿命降低, 也会使TPD的效率下降, 在实际设计时, 可选取温度在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的设计难度.4. 结 论
本文详细研究了TPD的性能特性与参数优化设计, 结果表明通过优化电压和阴极板功函数可显著提高TPD的功率密度和效率, 同时通过折衷考虑功率和效率, 可确定TPD性能参数的优化区间. 另外, 本文也分析了高温热源温度对TPD优化性能特性的影响, 确定了高温热源温度的优化范围. 本文所得结果可为高性能热电功率器件的设计与研制提供理论指导.
[1] 王长宏, 林涛, 曾志环 2014 物理学报 63 197201
Google Scholar
Wang C H, Lin T, Zeng Z H 2014 Acta Phys. Sin. 63 197201
Google Scholar
[2] 周敏, 黄荣进, 李来风, 蒋明波, 吴智雄 2010 物理学报 59 7314
Google Scholar
Zhou M, Huang R J, Li L F, Jiang M B, Wu Z X 2010 Acta Phys. Sin. 59 7314
Google Scholar
[3] Zeng T 2006 Appl. Phys. Lett. 88 153104
Google Scholar
[4] Lee J H, Bargatin I, Melosh N A, Howe R T 2012 Appl. Phys. Lett. 100 173904
Google Scholar
[5] Lee J I, Jeong Y H, No H C, Hannebauer R, Yoo S K 2009 Appl. Phys. Lett. 95 223107
Google Scholar
[6] Smith J R, Bibro G L, Nemanich R J 2007 Phys. Rev. B 76 245327
Google Scholar
[7] Liao T 2019 IEEE Electron Device Lett. 40 115
Google Scholar
[8] Meir S, Stephanos C, Geballe T H, Mannhart J 2013 J. Renew. Sustain. Energy 5 043127
Google Scholar
[9] Olawole O C, De D K 2018 J. Photon. Energy 8 018001
[10] Post A D, King B V, Kisi E H 2017 Appli. Therm. Eng. 117 245
Google Scholar
[11] Xiao L, Wu S Y, Yang S L 2018 Int. J. Energy Res. 42 656
Google Scholar
[12] Liao T, Chen X, Lin B, Chen J 2016 Appl. Phys. Lett. 108 033901
Google Scholar
[13] Schwede J W, Bargatin I, Riley D C, et al. 2010 Nature Mater. 9 762
Google Scholar
[14] Zhang X, Zhang Y, Ye Z, Li W, Liao T, Chen J 2018 IEEE Electron Device Lett. 39 383
[15] Yang Z, Peng W, Li W, Chen X, Chen J 2018 J. Appl. Phys. 124 154501
Google Scholar
[16] Liang S, Ang L K 2015 Phys. Rev. Appl. 3 014002
Google Scholar
[17] Liang S, Liu B, Hu W, Zhou K, Ang L K 2017 Sci. Rep. 7 46211
Google Scholar
[18] Misra S, Upadhyay K M, Mishra S K 2017 J. Appl. Phys. 121 065102
Google Scholar
[19] Mishra S K, Kahaly M U, Misra S 2017 Int. J. Therm. Sci. 121 358
Google Scholar
[20] 廖天军, 杨智敏, 林比宏 2014 中国科学: 物理学 力学 天文学 44 125
Liao T J, Yang Z M, Lin B H 2014 Sci. Sin.: Phys. Mech. Astron. 44 125
[21] Wang Y, Su S, Lin B, Chen J 2013 J. Appl. Phys. 114 053502
Google Scholar
[22] Chen L, Ding Z, Sun F 2010 J. Appl. Phys. 107 104507
Google Scholar
[23] Zhang X, Pan Y, Chen J 2017 IEEE Trans. Electron Device 64 4594
Google Scholar
[24] Kuznetsov V I, Ender A Y, Babanin V I 2018 J. Appl. Phys. 124 044502
Google Scholar
-
图 2 给定不同
ΦC 时, (a)阴极温度和(b)阳极温度随输出电压变化曲线及(c)伏安特性曲线, 其中ε=0.20 ,EF=0.80eV ,UH=UL=0.2W⋅cm−2⋅K−1 ,TH=1500K 和TE=300K Fig. 2. (a) Curves of the cathode temperature, and (b) the anode temperature varying with the output voltage, and (c) the volt-ampere characteristic for given values of
ΦC , whereε=0.20 ,EF=0.80eV ,UH=UL=0.2W⋅cm−2⋅K−1 ,TH=1500K , andTE=300K . -
[1] 王长宏, 林涛, 曾志环 2014 物理学报 63 197201
Google Scholar
Wang C H, Lin T, Zeng Z H 2014 Acta Phys. Sin. 63 197201
Google Scholar
[2] 周敏, 黄荣进, 李来风, 蒋明波, 吴智雄 2010 物理学报 59 7314
Google Scholar
Zhou M, Huang R J, Li L F, Jiang M B, Wu Z X 2010 Acta Phys. Sin. 59 7314
Google Scholar
[3] Zeng T 2006 Appl. Phys. Lett. 88 153104
Google Scholar
[4] Lee J H, Bargatin I, Melosh N A, Howe R T 2012 Appl. Phys. Lett. 100 173904
Google Scholar
[5] Lee J I, Jeong Y H, No H C, Hannebauer R, Yoo S K 2009 Appl. Phys. Lett. 95 223107
Google Scholar
[6] Smith J R, Bibro G L, Nemanich R J 2007 Phys. Rev. B 76 245327
Google Scholar
[7] Liao T 2019 IEEE Electron Device Lett. 40 115
Google Scholar
[8] Meir S, Stephanos C, Geballe T H, Mannhart J 2013 J. Renew. Sustain. Energy 5 043127
Google Scholar
[9] Olawole O C, De D K 2018 J. Photon. Energy 8 018001
[10] Post A D, King B V, Kisi E H 2017 Appli. Therm. Eng. 117 245
Google Scholar
[11] Xiao L, Wu S Y, Yang S L 2018 Int. J. Energy Res. 42 656
Google Scholar
[12] Liao T, Chen X, Lin B, Chen J 2016 Appl. Phys. Lett. 108 033901
Google Scholar
[13] Schwede J W, Bargatin I, Riley D C, et al. 2010 Nature Mater. 9 762
Google Scholar
[14] Zhang X, Zhang Y, Ye Z, Li W, Liao T, Chen J 2018 IEEE Electron Device Lett. 39 383
[15] Yang Z, Peng W, Li W, Chen X, Chen J 2018 J. Appl. Phys. 124 154501
Google Scholar
[16] Liang S, Ang L K 2015 Phys. Rev. Appl. 3 014002
Google Scholar
[17] Liang S, Liu B, Hu W, Zhou K, Ang L K 2017 Sci. Rep. 7 46211
Google Scholar
[18] Misra S, Upadhyay K M, Mishra S K 2017 J. Appl. Phys. 121 065102
Google Scholar
[19] Mishra S K, Kahaly M U, Misra S 2017 Int. J. Therm. Sci. 121 358
Google Scholar
[20] 廖天军, 杨智敏, 林比宏 2014 中国科学: 物理学 力学 天文学 44 125
Liao T J, Yang Z M, Lin B H 2014 Sci. Sin.: Phys. Mech. Astron. 44 125
[21] Wang Y, Su S, Lin B, Chen J 2013 J. Appl. Phys. 114 053502
Google Scholar
[22] Chen L, Ding Z, Sun F 2010 J. Appl. Phys. 107 104507
Google Scholar
[23] Zhang X, Pan Y, Chen J 2017 IEEE Trans. Electron Device 64 4594
Google Scholar
[24] Kuznetsov V I, Ender A Y, Babanin V I 2018 J. Appl. Phys. 124 044502
Google Scholar
计量
- 文章访问数: 8404
- PDF下载量: 64