高温熔融盐为热流体的单罐蓄热性能研究

高温熔融盐为热流体的单罐蓄热性能研究

随着能源危机、温室效应等环境问题日益严重，新能源和可支配能源的开发已成为世界共同主题。由于其对环境没有污染、温室气体和能源再生的优点，太阳能热发电逐渐成为一个潜力巨大的可支配能源技术。简单的太阳能热发电系统存在太阳能发电周期和用电需求周期不匹配,太阳辐射能流密度的突然变化导致的系统热应力的剧烈变化等问题.为了提高发电效率、减少发电成本、提高太阳能热电系统的稳定性和连续性,则需要对太阳能热发电系统增加蓄热装置,以使系统在没有太阳辐射能量的时候能继续满足发电需要.单罐蓄热系统是指作为蓄热介质的冷流体和热流体都储存在一个单罐中,在蓄热或放热过程中,冷流体和热流体会相互接触,在接触区域形成一个温度斜温层.斜温层以上流体保持高温,斜温层以下的流体保持低温,随着换热过程的进行,斜温层会上下移动,最后流出蓄热罐.为了缩短斜温层的距离,防止冷热流体对流混合,增加蓄热量,一般会在罐内填充石英岩或石英砂等材料来增加斜温层的效应.单罐蓄热系统的好处是投资费用比双罐蓄热系统节省了约35%,但是由于冷热流体的导热和对流作用,真正实现温度分层有一定困难.Pacheco等对熔盐作为蓄热介质的单罐槽式太阳能发电站进行了分析研究,对斜温层系统进行了理论模拟和实验分析,结果表明熔盐斜温层单罐蓄热系统是一个可行的蓄热方法.YangZhen等也对熔盐斜温层的热性能和放热效率进行了分析.Abdoly等对水作为蓄热介质的斜温层蓄热系统进行了理论和实验研究.虽然很多学者对蓄热罐填充介质的材料和蓄热效率等进行了研究,但如何较全面地评价蓄热系统性能的研究却较少.本研究对熔盐作为传热流体的斜温层单罐蓄热系统进行了计算分析,从不同角度对蓄热系统的蓄热性能进行了定义和研究,对主要影响蓄热系统性能的填充介质的物性进行了比较分析,为进一步研究系统性能提供了理论计算基础。1填充材料的选择和控制方程高温熔融盐单罐蓄热系统的数值计算模型如图1所示.蓄热罐的高为2m,直径为1m,熔盐和填充介质的初始温度为573K,高温熔盐进口温度为773K.熔盐进口速度为0.001m/s.蓄热初始时刻,蓄热罐中充满了573K的低温熔盐和多孔填充介质,蓄热过程开始后773K的高温熔盐从顶部进入蓄热罐,经过与多孔填充介质的换热后,推动低温熔盐流出蓄热罐.计算过程中考虑了熔盐变物性参数的影响,熔盐的物性参数如下所示:比热为1510J/(kg·K),导热系数为0.571W/(m·K),熔盐的密度与温度的关系为ρ=1938.0-0.732(T-200.0)(1)式中:ρ为密度,T为温度.熔盐的黏度与温度的关系为μ=exp[-4.343-2.0143(lnT-5.011)](2)式中:μ为黏度.蓄热罐中的填充介质为陶瓷、铝硅合金和砂石,主要是考虑3种材料的导热系数、密度和热容对蓄热性能的影响规律,从而选择合适的填充材料.填充材料形状都为球形,保持相同的孔隙率.填充介质的主要物性参数如表1所示,计算控制方程如下.连续性方程:∂(ερ)∂t+Δ⋅(ρu)=0(3)∂(ερ)∂t+Δ⋅(ρu)=0(3)式中:ε为孔隙率;t为时间;u为流体速度.由于流动属于层流范围,动量控制方程采用Darcy定律方程:ρ∂u∂t+ρ∇(u⋅uε)=−∇(εp)−εμKu(4)ρ∂u∂t+ρ∇(u⋅uε)=-∇(εp)-εμΚu(4)式中:p为压力,其中渗透率K可按照下式进行计算:K=ε3d2s15(1−ε)2(5)Κ=ε3ds215(1-ε)2(5)式中:ds为填充颗粒直径.能量方程:(ρcp)m∂T∂t+(ρcp)fu⋅∇T=∇(ρcp)m∂Τ∂t+(ρcp)fu⋅∇Τ=∇·[ue84aSymbollApm∇T](6)式中:(ρcp)m=(1-ε)(ρcp)s+ε(ρcp)f;ue84aSymbollApm=(1-ε)·ue84aSymbollAps+εue84aSymbollApf;cp为比热;ue84aSymbollAp为导热系数;下标s表示填充材料;f表示传热流体;m表示平均值.文中采用Fluent6.2流体力学软件进行数值模拟计算,Gambit网格划分为结构化网格,边界类型定义进口边界为VelocityInlet,出口边界为Outflow,其它边界为Wall,管内流体为Fluid和多孔介质区域.求解采用分离式求解器,层流模型、压力-速度耦合为Simple算法,压力、动量和能量采用二阶迎风格式.采用非稳态迭代求解,总共蓄热时间设定为40min.2评估蓄热性能的方法2.1初始温度tini在单罐蓄热系统中,一般蓄热进口温度和流量保持恒定值,出口温度随着时间而变化.蓄热过程中,当出口温度与熔盐初始温度相同时,说明蓄热罐的蓄热能力最大,当出口温度与进口温度相同时,说明蓄热罐已经没有蓄热能力了.因此文中以熔盐进出口温度差与进口温度和初始温度差的比值的大小来定义蓄热罐的蓄热能力.η1(t)=Tin−Tout(t)Tin−Tini(7)η1(t)=Τin-Τout(t)Τin-Τini(7)式中:η1为蓄热罐蓄热效率;Tini为初始温度;Tin、Tout为流体进出口温度.上式表明在初始时刻,出口温度和初始温度相同,其蓄热能力为最大值1.随着蓄热时间的增加,当出口温度与进口温度相同时,蓄热罐没有继续增加蓄热量的能力,其蓄热能力为0.Zurigat等定义斜温层蓄热器的蓄热效率如下:ηZG=∫td0m˙cp|Tout−Tin|dtmstorecp|Tin−Tini|(8)ηΖG=∫0tdm˙cp|Τout-Τin|dtmstorecp|Τin-Τini|(8)式中:m˙m˙为质量流量;mstore为罐体总质量.从上面蓄热效率的定义可知随蓄热时间增加,其蓄热效率从0逐渐增加到1,反映了瞬时蓄热量与最大蓄热量随时间增加而变化的规律.但其不能反映出蓄热能力的变化情况,当出口温度与进口温度相同时,说明蓄热罐已经达到最大蓄热能力,不能使蓄热量继续增加,但上式定义的蓄热效率却达到最大值1.2.2对蓄热罐体熵控制ds单罐斜温层蓄热过程中,由于高温流体、低温流体和多孔介质材料之间的导热和对流换热过程会使系统的熵增加,随着蓄放热过程时间的增加,其总的熵产值也随着增加.假设蓄热罐与周围环境没有热损失,其系统的熵的变化满足以下方程:ΔSstore=ΔSin-ΔSout+ΔSg(9)式中:ΔSstore、ΔSin、ΔSout、ΔSg、ΔSflow分别表示罐体,流体进、出口,系统和流体熵产.那么ΔSg=ΔSstore-ΔSflow.在蓄热时间为t的蓄热过程中,对熵流,ΔSflow(t)=∫t0m˙(t){sin(t)−sout(t)}dt(10)ΔSflow(t)=∫0tm˙(t){sin(t)-sout(t)}dt(10)ΔSflow(t)=∫t0cp,fuAρflnToutTindt(11)ΔSflow(t)=∫0tcp,fuAρflnΤoutΤindt(11)式中:sin、sout分别表示流体进、出口单位流量的熵产;A为罐体截面积.对蓄热罐熵变,ΔSstore(t)=∫mstore00mstore[s(t)-s(0)]dm(12)ΔSstore(t)=∫l0m˙store{cp,fε+cp,s(1−ε)}ΔSstore(t)=∫0lm˙store{cp,fε+cp,s(1-ε)}·lnT(t)Tinidl(13)lnΤ(t)Τinidl(13)式中:m为质量,l为长度;s为单位质量熵产.随着蓄热过程的进行,系统总熵产逐渐增加,当出口温度与进口温度相同时,蓄热系统没有对流和导热的换热过程,其熵产达到最大值.2.3蓄热有效性的确定在单罐斜温层蓄热系统中,由于冷热流体的导热和对流作用所形成的斜温层会降低系统的有效蓄热量.因此斜温层的存在会影响蓄热罐的蓄热有效性,这里定义蓄热有效性为η2(t)=1−δ(t)L(14)η2(t)=1-δ(t)L(14)式中:δ为斜温层的厚度,L为蓄热罐的总长度.由式(14)可知,在初始时刻,斜温层的厚度为0,其蓄热有效性为1,随着时间的增加,虽然蓄热罐进出口的温度没有变化,但斜温层厚度却逐渐增加,蓄热有效性逐渐下降.因为在蓄热的最后阶段,斜温层有一部分已经流出蓄热罐,因此该公式要求蓄热的时间内能保证斜温层是完整的.3结果和讨论3.1蓄热量t0cp蓄热过程中,熔盐进出蓄热罐的质量流量保持定值,由于熔盐的流动传递给蓄热系统的热量表示为:ΔQflow(t)=∫t0m˙(t){hin(t)−hout(t)}dt(15)ΔQflow(t)=∫0tm˙(t){hin(t)-hout(t)}dt(15)ΔQflow(t)=∫t0cp,fuAρf(Tin-Tout)dt(16)式中:ΔQ为蓄热量;h为焓.因此熔盐放出的热量可以通过计算得到.对蓄热罐体中的多孔介质和流体进行分析,其温度由初始温度升高到最终温度,那么蓄热罐吸热的热量为ΔQstore(t)=[cp,fεALρf+cp,s(1-ε)·ALρp,s]{Taverage(t)-Tinitial}(17)式中:Taverage、Tinitial分别为平均和初始温度.因此蓄热罐吸收的热量可以通过数值模拟得到.根据热量平衡原理,放出的热量应该等于吸收的热量ΔQflow(t)=ΔQstore(t)(18)图2为当多孔介质为陶瓷时,ΔQflow和ΔQstore随着时间变化的规律图.由图可知,ΔQflow和ΔQstore基本保持一致,说明模拟计算过程中蓄热罐的温度分布规律基本是正确的,也验证了模型的有效性.3.2多孔填充材料对蓄热过程的影响图3为当多孔材料为陶瓷的时候,不同蓄热时刻沿着高度方向的温度分布.可以看出在不同的蓄热时刻,高度方向上温度都由773K过渡到573K,而且随着蓄热时间的增加,斜温层逐渐向下移动,最后移出蓄热罐,蓄热罐被高温流体充满,完成蓄热过程.同时可以看出斜温层厚度逐渐增加,也就是温度梯度逐渐减少.图4为不同填充介质的时候,其出口截面的温度随着时间的变化关系.由图可知在蓄热过程的大部分时间,出口截面的温度都为573K.当多孔材料为砂石的时候,其出口截面温度首先上升,进行到1600s时,温度已经达到了773K,完成了整个蓄热过程.当多孔材料为陶瓷时候,温度分布规律与砂石类似,但完成蓄热过程要增加约400s的时间.当多孔材料为铝硅合金时,蓄热进行到大约1000s左右,出口温度由573K逐渐增加,但是到了2400s的时候,其温度仍然没有接近773K,说明其完成蓄热过程需要更长的蓄热时间.3.3多孔介质中蓄热量的影响蓄热罐所能储存的热量是评价蓄热罐蓄热能力的一个重要指标.由图5可知,蓄热开始后的一段时间内,不同多孔材料的蓄热量都相同.当多孔介质为砂石时,首先达到最大蓄热量,但其值最小,这主要是因为砂石的单位体积热容较小的原因.当材料为陶瓷时,其最大蓄热量要大于砂石的蓄热量.当多孔介质为铝硅合金时,最后阶段的蓄热量虽然没有保持定值,但其蓄热量要大于多孔介质为陶瓷时候的蓄热量,这主要是因为铝硅合金的单位体积热容较大.因此,单位体积热容是决定最大蓄热量的主要因素,但同时要考虑其所需要的时间.3.4口温度的大小图6为不同填充介质时的蓄热能力随着时间的变化关系.由于蓄热能力和出口温度有关,因此蓄热能力从一方面反映了出口温度的大小。从图中可以看出,由于导热系数和体积热容的不同造成温度梯度的不同,从而造成蓄热能力的不同.在初始阶段蓄热能力都为1,导热系数大的铝硅合金的蓄热能力首先开始下降,而且其下降斜率较小.导热系数较小的陶瓷和砂石下降斜率相同,但热容较小的砂石首先由1变为0,完成蓄热过程.因此应选择导热系数较小的材料作为填充材料.3.5不同填充介质对总熵产的影响图7为不同蓄热介质在不同蓄热时刻的系统总熵产图.由图可知,熵产随着时间的增加而增加,这主要是时间的增加使不可逆换热过程持续进行.当填充介质为陶瓷和砂石时候的熵产要明显小于填充介质为铝硅合金时所产生的熵产,这主要是因为铝硅合金的导热系数较大,从而使导热传热传播更深,熵产增大,不可逆性增大.陶瓷的熵产要高于砂石的熵产,这主要是由于陶瓷的比热容较大使传热持续的时间更长.可见蓄热量大的介质,其系统的熵产也较大.3.6不同蓄热介质对热有效性的影响图8为基于斜温层厚度定义的蓄热有效性随时间的变化关系.由图可知,随着时间增加,蓄热有效性明显降低,这主要是由于随着蓄热过程的进行,斜温层厚度逐渐增加的缘故.硅铝合金的蓄热有效性明显低于其他介质,说明其斜温层厚度非常大.当蓄热介质为砂石的时候,其蓄热有效性要略低于陶瓷的蓄热有效性,说明在蓄热过程中比热容大的介质的斜温层厚度较小.4不同导热系数时系统的蓄热行为文中从