1. 引言
由于人类频繁的工业生产活动造成温室气体排放和气候反常等问题,近年来国际上通过签订《巴黎协定》 [1] 来寻求合作,以期解决这些问题。截至2019年4月,已有197个国家签署了该协定,其中绝大多数国家倡导开发可再生能源,降低碳排放,从而保护环境。《巴黎协定》实现的一个关键标志是可再生能源生产和技术应用的大规模商业化,并减少传统化石燃料的碳排放 [2]。可再生能源技术可以解决能源需求不断增长的当务之急。随着可再生能源技术得发展,分布式发电是其中一个具有挑战性的前沿技术,众多科学家对发展分布式发电有着浓厚的兴趣 [2] [3] [4] [5]。温差发电(TEPG)由于拥有清洁可再生,无噪声,发电简单成为科学家们研究的热门 [6] [7]。人们熟知的光伏发电技术是收集光子能量并将其转化为电能,而温差发电技术由于物理中的塞贝克效应(Seebeck Effect) [8] 可直接将热能转化为电能。
根据塞贝克效应 [8],热电[TE]效应是通过将温差转换为电压来产生电能,反之亦然。决定热电效率的因素包括温差和热电优值ZT,其中ZT是热电材料特性的最重要参数之一,具有高ZT值的热电系统则具有更高的输出功率。另外,在一些特殊的航天应用中,热电材料具有重量轻、体积小、无噪音、可靠性好等重要性质。此时,这些性质比生产成本和工作效率更为重要。
温差发电是一个很有前景的技术,当材料的ZT值超过4时,温差发电比传统的发电机或光伏电池更有优势。根据文献 [9],ZT值的提高将促进许多技术应用,温差发电技术就是其中之一。当ZT = 4时,与太阳能薄膜发电相比,热电发电效率更高,可达30%。
热电模块由大量热电材料集成。热电材料通常由半导体制成,半导体能将热能直接转化为电能。每对热电偶每度温差可产生几百微伏的电压。所以,当环境存在温差时,热电设备会产生电压,外接负载即可产生电流。刘大为等 [10] 通过将大量的热电模块串并联产生输出功率为1 kW的热电设备。热电能量转换是一个非常有趣的领域,已经吸引了大量关于各种热能转换为电能等课题的研究 [11] [12] [13] [14]。
2. 热电模块原理说明
现如今,温差发电技术还未广泛应用与实际生活,主要是因为热电转换效率不高,而影响热电材料效率最重要的一个因素是热电优值ZT值。对于单一的热电材料,其与热电转换有关的性能可以用材料的热电优值ZT来描述,如下式所示:
(1)
式中:ZT是材料的热电优值,s是塞贝克系数,T是绝对温度,
是电阻率,k是导热系数。
热电模块一般使用N型和P型半导体等不同类型热电材料组成热电偶,一对热电偶具有组合的ZT值。其综合热电优值ZNPT描述如下:
(2)
式中:
和
分别为P型和N型的塞贝克系数。
和
分别是P型和N型的电阻率。
和
分别为P型和N型的导热系数。T是绝对温度。
输出电压ΔV随温差ΔT的增加而增加,并基本呈线性变化。对于N型和P型材料的组合,ΔV的关系式如下所示:
(3)
常数SNP称为热电偶对的塞贝克系数,它等于SP − SN,此常数与热电材料属性有关 [6] [7]。
根据公式(4a)和(4b)可知,热电效率η与ZT值有关。其转换效率η是输出功率P与系统的热输入Qh之比。
(4a)
(4b)
式中:Th、Tc、Tm分别为热电模块热侧温度、冷侧温度和冷热侧平均温度。ZTm是它在平均温度下的热电优值。
由上述可见,提高ZT值是提高热电转换效率的关键因素,而ZT值受塞贝克系数(S)、导热系数和电阻率 [15] [16] [17] [18] [19] 等参数的影响。温差发电技术要在大规模发电领域取得商业上的成功,一个重要的因素是效率需达到一定的级别,例如,
,η可达到20%。热电模块商业化后,其应用意义重大,包括应用于分布式发电、为电子设备提供电能,以及对可能被浪费的热能进行回收利用。为了使热电模块应用于商业化具有可行性,作者为提高其效率,从而获得较大的输出功率,进行了很多研究 [10] [15] [16] [20]。
3. 试验
3.1. 试验原理
图1简要地展示典型的温差发电原理及装置。典型的热电偶由两种不同的半导体材料组成。当热电材料两侧存在温差ΔT,如图1(a)所示,温差可使热电材料电荷迁移并产生电动势。
典型的热电模块(TEM)通常是用大量热电偶组成,如图1(b)所示。此外,图1(c)为温差发电模型示意图。在液化天然气(LNG)加热下,ΔT达到200℃以上时,总输出功率可达百瓦级。

Figure 1. The operating principle diagram of a TEPG: (a) schematic diagram of a TE device composed of P- and N-type TE legs; (b) a typical TEM integrates a large number of TE devices; (c) the physical diagram of TEPG with an output power of more than 100-Watt
图1. 温差发电设备的工作原理图:(a) PN型热电部件示意图;(b) 热电部件集成的热电模块示意图;(c) 输出功率可达100 W以上的温差发电实物图
温差发电系统可根据实际应用设计成多种多样。在设计温差发电系统时,研究人员应谨慎选择热管理方法,以防热电模块过热损坏。同时热电材料需形成冷却表面和加热表面之间的导热,从而促使载流子迁移,产生电动势。加热热源包括红外热源、液化天然气和炊事用火、加热器、太阳热能、地热能等。上面的选择为各种温差发电模型奠定了应用基础。
3.2. 热电模块原位测试仪
热电模块热电转换效率测试设备是测量热电模块特性的重要装置。金安君等发明了测量精度极高的单热电模块的热电特性测试仪器 [17] [21],具体说明如下。热电模块原位测试仪(ICSTEM)的结构如图2所示。该仪器可以测量以下变量:在发电状态下,热电模块冷端的输出功率P、电流I、电压V、内阻R和热流QC。因此,研究人员得到了以下结果:1) 热电模块热电效率;2) I-V曲线;3) P-V曲线;4) 力响应因子。

Figure 2. Structural diagram of In-situ characterization station of TEM (ICSTEM)
图2. 原位热电测试仪器(ICSTEM)的结构图
热电模块的效率定义为输出功率与设备热侧输入热量的比值。该仪器通过机械结构设计基本上可以确保沿着垂直于模块的z轴方向引导热流。外力由扭矩扳手通过螺栓设定。弹簧和固定螺钉都可以在z方向上分离,以保证热板和冷板之间的精准平行度(误差不超过5微米)。每块板的面积为16 cm × 16 cm。通过冷却剂工质(水或液氮)控制,冷板温度设置在−80℃~25℃,而热侧温度通过加热装置可达到450℃以上。该台设备具有宽泛的温差测量范围。
综上所述,ICSTEM仪器经过精确设计制造,使热量沿z轴 [22] 流动,热电转换效率可以实时方便地现场测量。理论上,冷却能力QC计算方法如下:
(5)
式中:SPN是热电模块塞贝克系数,Th是其热边的温度,Tc是冷边的温度,I是电流,R是电阻,k是热阻。由公式(5)可知当热电模块在冷侧和热侧温度相同时,可达到最大冷却能力。
ICSTEM具有很好的热管理系统。热电模块应平放,并与冷板附近热流传感器的两个边缘相对平行。隔热层是由金属薄层和玻璃纤维片的复合材料制成的多层交替隔热结构。它的多层交替的隔热效果比空气要好,隔热热导率小于0.01 W/M∙K。保温层从12层到20层不等,工作性能比石棉材料好很多,原因如下。一方面,金属箔具有很高的红外反射能力,根据计算,多层膜可以将红外辐射热泄漏的影响降低到可以忽略的程度。另一方面,纤维玻璃导热系数较低,中间有一定的空气,可以显著降低热量泄漏的影响。此时,由z方向热损失引起的误差估计在1%以下。
3.3. 余热回收温差发电系统设计
宁波大学及华能清洁能源研究院温差发电课题组进行了余热回收温差发电的设计、设备开发和测试。采用计算机有限元分析仿真算法及辅助设计出温差发电设备,其方法如图3所示。首先,设计了车辆的排气和冷却系统 [23],以建立温差发电的边界条件。其次,利用传热模型和有限元分析,在Solidworks®模型中建立了 [24] 仿真算法,以实现系统优化。通过优化参数,研究人员利用商业热电模块构建了温差发电系统,并输出直流电流。该温差发电设备最初是在实验室使用柴油发电机产生的废气热流源进行测试,成功提供了能满足要求的电压和输出功率。温差发电设备先后安装在一辆小型皮卡载货汽车 [23] 上,应用于1000公里以上的道路运行测试。

Figure 3. Design, development and test flow chart of the TEPG unit
图3. 高功率温差发电系统设计、开发和测试流程图
基于上述方法,研究人员设计并构建了汽车尾气温差发电设备(ATEH)的应用模型。汽车尾气温差发电设备采用“3 + 2”夹层结构,3个冷却水箱夹着2层烟气管道。图4(a)展示了一个集成单元的结构示意图。热电模块夹在冷却水箱和烟气管道中间,每层包含几十个热电模块,参照图4(b),热电模块放置在3 × 9宫格中,即27个格子,每格布置一个模块。每个级联上都有指定的电路连接法,以满足所需的输出规格。如图4(c)所示,热量从左侧流入,冷却水箱位于其顶部、中间和底部的夹层。热量从左侧入口流入汽车尾气温差发电设备,从右侧流出。

Figure 4. Automobile exhaust TEPG device: (a) dual stack in sandwich structure; (b) a batch of TEMs configuration; (c) energy harvest unit
图4. 汽车尾气温差发电设备图。(a) 双级联夹层结构图;(b) 单层热电模块排列图;(c) 热能收集模块图
4. 结果和讨论
4.1. 数据和分析
热电模块包括现有市售模块 [10] 和自制模块 [16] [25]。如图5(a)和图5(b)所示,热电模块的基本电气特性特征如下:图5(a)为电流(I)与电压(V)曲线,即I-V曲线;图5(b)为最大输出功率(P)与电压(V)曲线,即P-V曲线。这些数据绘制了一个热电模块在三种不同温差情况下的典型图像。图5(a)显示了I-V曲线,显示了低内阻电源的典型电气图像;图5(b)显示P-V关系。这些曲线显示了电源输出功率的典型图像;输出功率随温差的升高而增加。图5(c)中发现了热电模块输出功率与施加的夹紧力的关系。在力小的区域,功率随着力显著增加。这条曲线显示出了总变化的20%,这种图像在不同的温差设置中经常可以观察到。功率曲线通常取决于它表面所受的压力,如图所示,选择4800到12,000 (kg/m2)的压力,曲线较为平整,输出功率较为稳定。本研究选择设定外力压力值为6000 (kg/m2)。功率随压力的变化可能是由于接触孔、接触电阻以及热电模块和界面之间的空隙缺陷造成的。通过设置合适的夹紧力,可以优化热电模块的上述缺陷。优化设置可以减少寄生电损耗,可降低接触电阻,减小接触界面的空隙热效应。热电能量转换效率的测量值取决于工作温差,如图5(d)所示。
4.2. 热量采集方法的研究
课题组构建了汽车尾气温差发电设备(ATEH)的三维模型。他们在多个热电模块单元组合的研究基础上,发现不管是串联还是并联,其总输出功率都能理想的叠加,即等于所有单个温差发电 [10] 单元的输出功率之和。汽车尾气温差发电设备的双管道结构既可以引入热气体流过,也可以引入在温差发电系统中流动的冷却介质。如图4(b)所示,为了满足所需的输出规格,几十个热电模块串并联成所需的电路。图4(a)说明含有夹层的汽车尾气温差发电设备组件。如图4(c)所示,热量从左侧流入,冷却水位于其顶部、中部和底部的三层夹层。

Figure 5. Test result diagram: (a) I-V curve; (b) P-V curve; (c) P-F curve; (d) the TE efficiency shows dependence on the temperature differential
图5. 测试结果图。(a) I-V曲线;(b) P-V曲线;(c) P-F曲线;(d) 热电效率与温差关系曲线

Figure 6. Temperature distribution of thermoelectric modules (3 × 9)
图6. 热电模块温度分布图(3 × 9宫格)
以设计模型为例,给出了一个典型的XY平面温度变化,如图6所示。温度变化符合预期,这是由于接近热电模块边缘的热气体的粘度比中心高,而且模块边缘的气体温度较低。沿y轴方向靠近中心的温度应较高。
构建并测试各种模型。通过实验室对温差发电系统(ATEH)的测试,总输出功率达到1 kW (315℃时1080 W)。实验测得最大输出功率和测试范围内温度关系图7所示:
除上述实验室试验外,还成功地制造了样机,并将其安装在一辆小型货车皮卡上进行道路测试。当地一家制造商 [23] 提供了一款皮卡,其发动机型号为萨博V系列,发动机型号为4Q20M。研究人员在这辆车里安装了1个或2个单元,而在这个实验室里安装了4个单元的汽车尾管能量采集器。每一个比较实验,道路测试行驶了1000公里 [23]。测试显示,使用汽车尾气余热利用装置,汽油行驶里程显著提高。
4.3. 多级联的温差发电
多级联温差发电技术使用先进的制造,把多种块体热电材料在合适温度区进行恰当的层叠;从而获得高性能的优化热电材料。这些优化材料可达到比目前市场上的热电设备更高的热电转换效率。原因是现有热电材料工作温度范围较小,其中效率相对较低。在目前的情况下,热电优值ZT普遍不够高。
利用MatlabTM软件 [25] 对P型和N型热电偶的各节段长度分别进行建模分析并优化。将P型或N型热电偶按热流方向等温差分成m + 1段。如果m足够大(例如,大于1000),则每段热电偶材料性质(塞贝克系数、电阻率和导热系数)可视为常数,并根据温度赋值。以P型热电偶为例(N型热电偶也可以用同样的方法模拟计算),该热电偶由两部分组成:第1部分从第1段到x段,使用热电材料Pb0.94Ag0.01La0.05;第二部分,从(x + 1)段到(m + 1)段,使用了热电材料Bi2Te3。
分别为冷侧和热侧温度的第i段长度。假设横截面面积(A)和传热流量(Q0)为单位值,然后根据材料性质得到各段P型热电偶的长度(
)。每个部分都被视为一个热电电池,因此整个P型热电偶由这些部分串联而成。第i节的输出功率可以由文献 [26] 所示方式得到。该系统可以改进如下。如图8所示,研究了热电器件的多级联结构,以提高ZT,达到较大的热电效率。横坐标表示高温段长度与总长度(x)之比,高温段的最佳长度比随热侧温度的升高而增大。当x值很小时(热电偶的大部分是由低温材料制成的),热侧温度越高并不意味着效率越高。这是因为低温材料的热电性能在温度超过其适当的工作范围后迅速下降,总效率被拉低,就像给电池添加一个劣质电池一样。
图8显示了P型和N型热电材料的效率与长度比(x)的关系,其中x为分段热电偶对所处长度比中的节点位置。仿真模拟的依据是上述方程(2),(4a, 4b)。图8(a)为设置三种温度下P型效率随x的函数关系;图8(b)为N型效率与x的函数关系,P型在550℃时热电效率为11.8%,而在550℃时N型在热侧效率为9.7%。注意,变量x是沿着每个双级联热电段的长度比。对多级联热电器件的结构进行了大量的计算机建模仿真模拟和实验测试。上述模拟依据为公式(2-4a, 4b)。

Figure 8. The diagram of P-type and N-type material efficiency and length ratio: (a) the P-type efficiency shows dependence on x; (b) the N-type efficiency depends on x.
图8. P型和N型材料效率与长度比关系图。(a) P型材料效率与长度比关系;(b) N型材料效率与长度比关系图
实验测试采用的热电材料性能列于表1。通过对市场上可用热电材料的选择,所选每级热电材料具有相应的温度,并且热电装置可以增加工作温度区间。热电材料的选择如表1所示。按照双级联优化值构建热电器件的测试效率达到10.6%,远高于单个级联BiTe效率。
如表1所示,双级的叠层热电材料是基于Pb-Te材料和Bi-Te材料组成。研究人员对P型和N型的基于Pb-Te为基础的热电材料,在中等温度范围内进行了优化;对基于Bi-Te为基础的热电材料,在低温范围内进行了优化。P型和N型分段热电偶上的基于Pb-Te材料和Bi-Te材料的双级联结构需通过高温硬焊接连接在一起,否则无法焊接在一起。双级联热电系统明显提高了ZT值。表1列出了其覆盖其高温工作区间和各种相关热电材料的优值系数ZT值。

Table 1. The variation law of ZT value of different materials with temperature. (P1-PbTe: Pb0.94Sr0.04Na0.02Te, N1-PbTe: Pb0.94Ag0.01 La0.05Te, P2-BiTe: Bi0.5Sb1.5Te3, N2-BiTe: Bi2Te3)
表1. 不同材料ZT值随温度的变化规律(P1-PbTe: Pb0.94Sr0.04Na0.02Te, N1-PbTe: Pb0.94Ag0.01 La0.05Te, P2-BiTe: Bi0.5Sb1.5Te3, N2-BiTe: Bi2Te3)
本文所制造的热电模块是基于BiTe和PbTe的热电材料组成。通过实测数据和基于公式(5)的简单计算,得出仿真模拟的效率值。
通过大量的材料调研及实践,得到优化的系统设计方案。他们发现通过对相关制造方法的研究,可以获得比之前的技术 [27] 更高的温差区、更高的能量密度、更好的热电输出效率等性能。
多级联温差发电系统在提高系统输出效率方面具有优势。如果单价的成本继续从现有热电模块的价格下降,理论上,当
时温差发电具有较好的性价比并可能实现大规模的商业应用 [27]。
ZT值与材料的电导率、导热系数、塞贝克系数等性能有关。分段热电偶所采用的热电材料的ZT值列于表1中。基于PbTe的热电材料的高温度工作区间为573 K到873 K,此时P型的热电优值系数范围约为1.5~2.1,N型为0.9~1.4。基于BiTe的热电材料的低温工作区间为室温到573 K,此时P型的热电优值范围约为0.7~0.9,N型为0.8~1.1。提高ZT值来改进热电材料已经引起广泛深入的研究 [15] [20] [26] [28]。
5. 总结
研究团队基于现有的材料而设计的装置,在实验室中研制出就双级联装置而取得10%以上的热电效率;就三级联装置的建模,效率达19%。本文中调研了多种温差发电仪器或装置的设计和模型;应用实验及仿真计算的方法,创新设计了高功率输出的热电系统;研发获得了1 kW以上的输出功率。他们在实验室里成功搭建先进的原位测试仪器,系统地研究热电响应的转换效率、I-V曲线、P-V曲线、及输出功率与外力的响应曲线。研究者还提出了一整套热电装置的设计及架构特性。最后,他们的研究方法可成功地获得高效率、提升热电装置的输出功率。
NOTES
*通讯作者。