太阳能热光伏系统性能分析研究终稿

1、太阳能热光伏系统性能分析研究 陈雪 宣益民 韩玉阁(南京理工大学动力工程学院,江苏 南京 210094摘要:本文建立了太阳能热光伏系统 (STPV 能量吸收及热辐射发电过程的数理模型, 数值 模拟了不同聚光器聚光比时辐射器的温度分布, 研究聚光比及辐射器温度对电池伏安输出特 性及系统效率的影响。 随着聚光比的增大, 辐射器表面温度显著增高, 电池的输出功率及系 统效率也成直线增长; 分析了影响系统效率的关键参数如滤波器、 电池温度和辐射器上下端 部材料表面反射率对系统性能与输出电功率的影响。计算结果证明了 STPV 系统的可行性, 形成了高效 STPV 系统的优化设计方法。关键词:太阳能热光伏

2、,辐射器,系统效率1引言太阳能热光伏系统(STPV 是在热光伏系统的基础上发展而成的一种新型光电转化装 置, 其基本原理是通过聚光器会聚太阳光来加热辐射接收器, 使其达到一定的高温, 再由光 电池将其发出的辐射能转化为电能。 STPV 系统采用清洁可再生的太阳能为热源,具有无污 染、 高转化效率、 体积小便携等优点。 目前, 国外对于太阳能热光伏系统的研究已投入广泛 关注, 如美国 EDTEK 公司 1, 俄罗斯物理研究所 2及美国 NASA 研究所 3, 都已通过 STPV 系统实验测试阶段, 产品正逐步推向市场及军事领域。 而国内对于 STPV 系统的研究还处于 空白阶段,尚未引起足够的重

3、视。一个典型的太阳能热光伏系统由两个基本模块组成:太阳能聚光器和光电转换装置。 而 一个光电转换装置主要包括:热辐射器、 光学滤波器、 光电池, 其它子系统是作为提高系统 的转换效率或者为能量的回收利用而适当增加的,如电池散热片、辅助组件等。本文分别以 STPV 系统中的能量吸收装置 (聚光器+辐射接收器 和热辐射发电模块 (辐 射器+滤波器+电池 为研究对象, 研究光电直接转换过程, 构建辐射接收器吸收太阳能及 瞬态传热过程的数学物理模型, 分析计算不同聚光器聚光比时辐射器表面的温度分布, 同时 对 STPV 系统的光电转化过程进行数值模拟,并分析各个参数对系统性能的影响情况,以期 建立 S

4、TPV 系统的优化设计方法。2 系统性能模型图 1为太阳能热光伏系统结构示意图,包括聚光器、辐射接收器,滤波器,电池组,上 图 1 STPV系统结构示意图Fig.1 Scheme of STPV system基金项目:江苏省自然科学基金重点项目 (No. BK2007726.下反射层等,辐射接收器采用耐高温灰体材料,滤波器为 Si/SiO2光子晶体滤波器 4,电池 为低禁带的族 GaSb 电池, 上下反射层用来反射辐射能以保证辐射器发出的辐射能被 电池充分吸收 1。 辐射器和电池组安装在同轴圆柱面上, 接收器空腔上端采光口直径随聚光 器聚光比及所产生光斑的大小变化,辐射接收器壁厚 2mm ,高

5、 3cm 。电池所在柱面直径为 4cm ,单块电池表面为 2cm×1cm 的方形结构。2.1太阳能聚光器在 STPV 系统中采用卡塞格伦型聚光器 1会聚太阳能量,该聚光器由两级反射镜组成, 一级反射镜为旋转抛物面反射镜, 二级反射镜为旋转双曲面反射镜, 抛物面型反射镜球差为零,而双曲面镜又可将射入其一个焦点的光束严格反射至另一个焦点,没有像差 5,因此一台精密的卡塞格伦型聚光器应将垂直入射的太阳光线严格会聚于一点。 由于两级反射镜加工 中不可避免的会产生误差, 这里假设垂直入射的太阳光线经聚光器会聚后会产生一定大小的 光斑,取聚光器聚光比变化范围为 20008000,一级反射镜入口直

6、径 60cm ,则聚光器所产 生的光斑直径在 6mm 13mm 范围内变化。其中聚光比定义为:2 ' /(a a CR = (1式中, a 一级镜入射孔径, ' a 经聚光器会聚后在焦点处形成的光斑直径假定太阳光线为平行光, 地面可接受的最大太阳辐射强度为 850W/m26。 两级反射镜的 反射率均为 0.921,部分太阳能量在反射过程中损失。2. 2 辐射接收器能量平衡方程的建立辐射接收器吸收由采光口进入的太阳能量, 温度升高, 对电池发出辐射能, 以实现光电 转换。在高温状态下,辐射器、电池、各反射面间的换热主要为表面辐射,可认为辐射器和 电池间的辐射换热发生在密闭腔中,

7、而在接收器空腔中, 存在能量的泄漏, 部分射入的太阳 能量以及辐射器发出的辐射能将由入口处泄漏, 从而引起能量的损失。 对模型作出以下假设:经聚光器会聚而来的太阳光束全部进入接收器空腔对其进行加热, 辐射接收器为耐高温 灰体材料, 其吸收率及发射率均为 0.85; 所有表面均为漫射表面; 只考虑辐射器表面的 热辐射, 由于其他表面的温度均远低于辐射器温度, 根据辐射四次方关系, 忽略上下反射面 及电池表面热辐射是合理的。2. 2.1 辐射器内部节点能量平衡方程针对辐射器内部第 (i, j, k个面元。根据能量守恒定律,得出该面元的热平衡方程: 0321=-+dt dT v Cp Q Q Q (

8、2 式中, 、 Cp 分别为辐射器材料的密度和定压比热容。Q 1为面元 (i, j, k接收到周向相邻面元的导热加热: -+-=+-+ ( (1, , 1, , 11, , 1, , 11k r T T z r k r T T z r Q n kj i n k j i n kj i n k j i (3式中, 1, , +n k j i T (n+1时刻面元 (i,j,k的温度, K ;Q 2为面元 (i, j, k接收到的纵向相邻面元的导热加热:z T T r k r z T T r k r Q n kj i n k j i n kj i n k j i -+-=+-+1, , 1, 1,

9、1, , 1, 1, 2 ( ( (4Q 3为面元 (i, j, k接收到的厚度方向上内外节点的导热加热: r T T z rk r r T T z rk r Q n kj i n k j i n kj i n k j i -+-+=+-+1, , 11, , 1, , 11, , 32 (2 ( (5针对辐射器内表面第 (i, j, k个面元(k=1 ,分析其与外界进行的各种热交换。根据能量 守恒定律,得该面元的热平衡方程:06547=-+dt dTv Cp Q Q Q Q (6(1 Q4为面元 (i, j, k接收到的经聚光器会聚而来的太阳热辐射能,采用射线追踪法跟踪由辐 射接收器入口处进

10、入的太阳光束的传播路径,统计面元 (i, j, k接收的能束数,设为 k j i N , , , 则 该面元吸收的总太阳能量为:N N Q Q kj i sun , , 4= (7(2 Q5为面元 (i, j, k接收到的其它面元的热辐射(包括自己在内所有面元对 (i, j, k的辐射贡 献 :-=m l h unit S T k j i m l h F Q m l h , , 4, , , , , , 5 (8式中, k j i m l h F , , , , -发射面元 (h, l, m对面元 (i, j, k的辐射传递系数;unit S 发射面元的面积。辐射传递系数采用蒙特卡罗法计算 7

11、。 对微元面 (h, l, m发出的 N 条辐射能束进行跟踪, 统计微元 (i, j, k吸收的来自发射微元面的能束数 k j i N , , ,即可确定面元 (h, l, m对面元 (i, j, k的辐射传递系数:N N F kj i k j i m l h , , , , , , =- (9(3 Q 6为面元 (i, j, k接收到的相邻面元的热传导。(4 Q 7为面元 (i, j, k向周围环境辐射散热量:unit k j i i S T Q 4, , 7= (102. 2.3 辐射器外表面节点热平衡方程的建立针对辐射器外表面面元 (i, j, k(k=M, M 为厚度方向划分的总节点数

12、 ,分析其与外界进 行的各种热交换。根据能量守恒定律,得该面元的热平衡方程:01098=-+dt dTv Cp Q Q Q (11(1 Q 8为面元 (i, j, k接收到的其它面元的热辐射,主要为经滤波器返回的能量。 -=+-1000, , 1, , 581 exp(, , , , m l h n m l h unit k j i m l h S T c c F Q (12 式中, 2/4810742. 31m m W c =, K m c =410439. 12辐射传递系数仍然采用蒙特卡罗法, 由于滤波器在不同波段具有不同的透射率, 则不同 波段内的辐射传递系数不同,这就要求在每个微元波长

13、范围内(取 m 02. 0进行射线跟踪, 当辐射能束到达滤波器时, 要产生一随机数与滤波器在该单元波长的透射率比较, 若随机数R (滤波器的光谱性能见 4, 能束可以透过滤波器, 被电池接收 (这里假设透过滤波器的能束均可被电池吸收 ,否则被反射,则继续跟踪直到此能束被一表面吸收为止。根据每个面 元接收到其他所有面元 (包括本面元在内 所发射能束的数量, 则可确定各面元在单元波长 区间的辐射传递系数。(2 Q 9为面元 (i, j, k接收到的相邻面元的热传导。(3 Q 10为面元 (i, j, k向周围环境辐射散热量。热平衡方程式中存在辐射项,属于非线性问题,采用如下方法对方程进行线性化处理

14、:4, , 1, , 3, , 41(, , 34n kj i n k j i n k j i n k j i T T T T -+ (13 式中, nk j i T , , 是上一时刻的节点 (i,j,k的温度,将其与常数项一起当作一个系数,则(13式是关于 1, , +n k j i T 的一个线性表达式。这样原来的非线性问题便转化为线性问题 8。计算时采用松弛迭代法, 将上一时刻算得的温度代入常数项中, 算得各辐射能量,重新求解温度场, 不断重 复以上步骤直到前后两次迭代温度场的平均值差 <1.0E-6时,则达到稳定,结束迭代。2.3 光电转化模块本文计算中所采用的 Si/SiO2

15、光子晶体滤波器光谱性 能如图 2所示 4, 由图可见, 该滤波器在 GaSb 电池可转换波段 0.71.73m 具有较高的透过率,使得辐射器发出的该波段的辐射光子可以透过滤波器到达电池, 转化为电能;而在其余波段透射率较低,因此又可将该部分的辐射光子反射回辐射器重新利用,保持辐射器在一定高温。由此可见,该滤波器和 GaSb 电池较为匹配,可对辐射光谱进行较好的控制,减少电池的废热,实现能量的循环利用。 图 2 Si/SiO2滤波器光谱特性1 (exp(0-=c S L T k qUJ J J (14式中, J L 为光伏电池的短路电流密度,其表达式为 10:=g QEd e hce J b L

16、 0 2/m A (15 式中, b e 为辐射器发出能量的光谱分布 , 为滤波器在不同波段内的透射率; QE 为电池内量子效率,即电池吸收的光子中转化为电子空穴对的百分比; hc 分别为普朗克常 数和真空中的光速。 g为电池的阈值波长。 J S 为反向饱和电流密度 11, 可采用经验公式 :2033/ *exp( 1084. 1(m A T k EgT J c c s -= (16式中, T c 为电池温度, E g 为电池禁带宽度, k 0为玻尔兹曼常数。电池最大输出功率表示为:FF J V P L oc =max (17 式中, V oc 开路电压可表示为 9:1ln(0+=S Lcoc

17、 J J e T k V (18FF 填充因子,其经验公式为 12: ln( ln(ln1 ln(11(Js J J Js J FF LL L-= (19电池填充因子综合了串联电阻及旁路电阻等不理想因素的影响。不能转化部分结合少量透过滤波器的禁带宽度范围以外的辐射光子一起转化为废热,被电池背面的散热器带走,以维持电池在一恒定温度,保证其正常工作性能。系统转化效率为电池输出功率和照射在一级反射镜表面的总太阳能量的比值:的总太阳能 照射在一级反射镜表面 输出电能stpv = (20由电池效率,光谱效率以及集热装置效率决定,分别定义为13: 的总太阳能 照射在一级反射镜表面 辐射器发出的净辐射能=c

18、ollector (21辐 射 器 发 出 的 净 辐 射 能s p e c t r 禁 带 宽 度 的 辐 射 能 通 过 滤 波 器 的 大 于 电 池 = (22 禁带宽度的辐射能 通过滤波器的大于电池 输出电能 cell = (23式中, 辐射器发出的净辐射能为热辐射器在一定温度下发出的全波段的辐射能减去经滤波器 返回并被辐射器重新吸收的能量 14 。3、计算结果与分析通过对以上能量平衡方程的迭代求解, 可以得出辐射器的温度分布随聚光器聚光比的变 化情况; 结合光电池求解经验公式及系统效率计算公式, 可得出聚光比、 光谱控制性能、 电 池温度等参数对系统输出功率及效率的影响。3.1聚光

19、器聚光比对辐射器温度的影响 图 3 聚光比为 2000时辐射器表面温度分布 图 4 聚光比对辐射器表面温度影响 Fig.3 Temperature distribution of emitter with the Fig.4 The influence of the CR on the emitter Concentrator ratio as 2000 temperature图 3示出聚光器聚光比为 2000,即太阳光经会聚后所产生的光斑直径为 1.2cm 时,辐 射器表面的温度分布情况,可以看出,辐射器表面最高温度可达到 1805.15K ,最低为 1715.83K , 分布比较均匀, 温

20、差的减小可使得电池表面热流密度更为均匀, 以避免电池局部 温度过高影响其转化能力及寿命。 表面温度由下至上呈一递减趋势, 由于太阳光线经卡塞格 伦型聚光器两次反射后, 在辐射接收器入口处的入射角很小, 因此大部分光线在照射到辐射 器下端区域后即被吸收, 而只有未被吸收的光线经反射后落到中间及上段区域, 因此各段会 存在一定的温差。 此外由图 3还可以看出, 辐射器表面温度沿周向分布还是比较均匀的。 图 4示出聚光器聚光比对辐射器温度分布的影响,由于表面温度沿周向温差很小,这里仅示出 纵向温度分布以作比较。 随着聚光比的增大, 光斑减小, 在相同的反射镜直径下, 进入接收 器空腔的能量相等, 随

21、着接收器采光口的减小, 由开口处泄漏的太阳能量及辐射能减少, 使 得能量得到更好的利用, 因此随着聚光比的增大, 辐射器的温度逐渐增高, 在聚光比为 8000时,表面最高温度可达到 1892K ,比同样条件下聚光比为 2000的辐射器温度高出 82K, 由此 可以看出, 聚光器的会聚能力对辐射器温度有很大影响, 应尽量减小光斑面积, 使得太阳能 量会聚于焦点。3.2 聚光器聚光比及辐射器温度对系统输出特性的影响图 5 示出了聚光器聚光比对电池输出特性及系统效率的影响，电池温度为 300K。由图 4 可以看出，随着聚光比的增大，辐射器表面温度升高，短路电流密度（和纵轴交点）和开 路电压（和横轴交

22、点）都随之增加，电池最佳工作状态点也逐渐向右上方移动。电池的输出 功率对聚光器聚光比也十分敏感，随着聚光比的增加，电池输出功率成直线增大，聚光器聚 光比为 2000 时，电池输出功率密度为 1.8W/cm2，当聚光比升高到 8000 时，电池输出功率 （a） 图 5 聚光比对系统性能的影响 (b Fig.5 The influence of the CR on the performance of the system 密度增大到 2.9W/cm2。由此可见，应提高聚光器的会聚能 力，在材料允许范围内尽量提高辐射器表面温度，以实现 系统效率的升高。 图 6 示出在各聚光比下， 系统效率随辐射器

23、表面温度的 变化。由图可见，在各个聚光比下，系统效率随辐射器表 面温度的变化趋势类似，随着辐射器温度的升高，系统效 率先增大，这是由于辐射器温度升高，其发出的辐射能量 显著增加，且其发出的辐射光谱和电池的可转换波段更为 匹配，使得透过滤波器的能量增加，因此系统效率增 图 6 系统效率随聚光比及辐射器温度的变化 高。但是当温度升高到一定值时，系统效率 Fig.6 System efficiency verse CR and emitter temperature 反而减小，这是因为随着辐射器温度的升高，由采光口射出的辐射能量增高，当损失的能量 超过一定范围时，系统效率降低。即存在一最优温度可使得

24、系统效率达到最大值。随着聚光 比的增大，该最优温度由 2000suns 时的 1750K 升高到 8000suns 时的 1950K，而最高系统效 率也由 21上升到 34。 3.3 光谱控制对系统的影响 图 7 示出采用性能较为优良的 Si/SiO2 光子晶体 滤波器和不采用滤波器时辐射器表面平均温度及电 池中所产生废热的比较。由图可见，有滤波器时的 辐射器温度要比没有滤波器时辐射器表面平均温度 高出近 300K， 这是由于滤波器可将电池无法转换波 段的光子反射回辐射器以重新利用。同时由图 7 可 以看出，滤波器可以大大减小电池中所产生的废热 量，在含有滤波器的 STPV 系统中，当辐射器表

25、面 图 7 滤波器对辐射器温度及系统性能的影响 温度达到 1890K 时，电池中产生的废热不超过 Fig.7 The influence of the filter on waste 5w/cm2，而不含有滤波器的系统中，当辐射器温度为 heat and emitter temperature 1600K 时，电池中产生的废热就高达 16w/cm2，电池中过多的废热将引起电池性能的下降， 增加了对散热器的需求，由此可见滤波器是至关重要的。 3.4 电池温度的影响 6 图 8 示出聚光器聚光比为 8000， 电池温度在 300380K 范围内变化时的电池输出特性以 及电池效率的变化情况。由图 8

26、（a）可以看出，电池温度的变化对短路电流几乎没有影响, 这由公式（16）亦可看出，但是随着电池温度的升高，电池开路电压逐渐减小，由 300K 时 的 0.49V 下降到 380K 时的 0.38V，电池的最优工况点也逐渐向左下方移动。电池的输出功 率对电池温度也十分敏感，当电池温度为 300K 时，电池输出功率密度为 2.9W/cm2，而当电 (a 图 8 电池温度对系统性能的影响 (b Fig 8 The influence of cell temperature on the STPV performance 池温度升高到 380K 时，系统的输出功率降为 2.3W/cm2。由图 8（b）

27、可以看出，电池温度 的升高对电池效率以及系统总效率也有很大的影响，当电池温度从 300K 升高到 400K 时， 电池效率降低了 15 个百分点，从而使得系统效率降低了 10 个百分点。由此可以看出，电池 温度的升高对系统性能有很大的负作用， 因此必须采用良好的散热系统以消除电池中产生的 废热，同时提高滤光器的性能，尽量使得不可转化波段的辐射光子反射回辐射器重新利用。 3.5 上下反射层性能的影响 图 9 示出辐射器和电池间区域上下端反射层的反射率对辐射器温度及系统效率的影响， 聚光器聚光比取为 8000，电池温度为 300K。由 图可见，随着反射层反射率的减小，辐射器表面 温度及系统效率均成

28、线型下降，随着反射材料对 辐射能的吸收率的增大，部分辐射能无法到达电 池表面，使得电池的输出功率减小，而电池无法 吸收的能量又不能完全反射回辐射器重新吸收， 也使得辐射器的温度降低，造成了能量极大的损 失。 可以看出， 反射层的反射率每降低 10,辐射 器表面平均温度降低 10K，而效率随之降低 1.8。 图 9 反射层的反射率对系统性能的影响 因此在太阳能热光系统的实际搭建中，要注意 Fig.9 The performance of system verse the 该项影响参数，尽量选用高反射率材料以减少能量的损失。 reflectivity of reflector 4、结论 本文构建了

29、 STPV 系统能量吸收及热辐射发电过程的数理模型，讨论了聚光器聚光比对 辐射器温度以及系统性能的影响，同时也分析了辐射器表面温度、电池温度、滤波器等对系 统效率的影响，得出结论如下： 聚光器的会聚能力对辐射器的温度有很大影响， 随着聚光器聚光比的增大， 辐射器表面温 度升高。 辐射器温度是影响系统性能的一个关键参数， 随着辐射器表面温度升高， 电池输出功率及 系统效率也随之升高。在每个聚光比下都存在一最优温度，当超过该极限值，辐射能的泄漏 量大大提高，反而使得系统效率降低。 光谱控制装置的滤光能力对系统性能有着很大的影响， 滤波器在提高辐射器温度的同时也 7 减小了电池中产生的废热。 电池温

30、度的升高将导致其输出功率及系统效率的减小， 因此， 保持电池温度稳定也是提高 系统性能的关键所在。 反射层的反射率降低使得辐射器表面温度及系统效率均有所下降， 因此在实际应用中， 应 尽量选择高反射材料，减少其对辐射能的吸收，提高能量的有效利用。 在进行数值模拟过程中，仍然忽略了一些参数的影响，如辐射器外表面和电池间不能完 全密封，且在该区域中不能保证完全真空，存在一定的对流换热，电池计算模型比较理想等 等，因此对系统模型的分析还有待进一步完善。 参考文献 1 Horne E. Hybrid Thermophotovoltaic Power Systems. Consultant Report

31、, 2002. 2 Andreev V M, et al. Solar thermophotovoltaic converters based on tungsten emitters. Transactions of the ASME, 2007, 129: 298-303. 3 Stone K. W, et al. Solar thermophotovoltaic power experiments at McDonnell Douglas. Proceedings of the 29th Intersociety Energy Conversion Engineering Confere

32、nce, 1994, 4:1692-1696. 4 刘广平, 宣益民, 韩玉阁.一维光子晶体在热光伏技术中的应用, 光子学报, 2007, 37(1: 115-119. 5 徐南荣, 卞南华. 红外辐射与制导. 北京: 国防工业出版社，1997. 6 沈辉, 曾祖勤. 太阳能光伏发电技术. 北京: 化学工业出版社, 2004. 7 卞伯绘. 辐射换热的分析与计算. 北京: 清华大学出版社, 1988. 8 宣益民. 计算传热学. 南京: 南京理工大学出版社, 1994. 9 刘恩科, 朱秉生等. 半导体物理学. 西安: 西安交通大学出版社, 1997. 10 Qiu K, Haydena A

33、C S, Mauk M G, et al. Generation of electricity using InGaAsSb and GaSb TPV cells in combustion-driven radiant sources. Solar Energy Materials and Solar Cells, 2006, 68(81: 68-81. 11 Ferguson L G., Frass L M. Theoretical study of Gasb PV cell efficiency as a function of temperature. Solar Energy Mat

34、erals and Solar cells, 1995, 39: 11-18. 12 MacMurray D W. Modeling and performance of microscale thermophotovoltaic energy conversion devices. IEEE Transactions on Energy Conversion, 2002, 17(1:130-141. 13 Zenker M. et al. Efficiency and power density potential of combustion-driven thermophotovoltaic sys