- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
高温相变蓄热材料熔化时间与空气出口温度的关系
在许多地方,太阳能作为一种重要的绿色能源使用广泛,尤其是在含有热热机的热能源发电系统的应用(集热器收集能量加热的工艺,驱动热机自发地发电)。由于太阳能具有间歇性的特点,如何有效地储存太阳能是目前研究的热点之一。太阳能热储存方式主要有三种:显热蓄热、相变蓄热和化学反应热蓄热。显热蓄热结构简单,已在太阳能热发电系统中得到广泛的应用(如在Solar One中,导热油被用作显热蓄热材料;在Solar Two中,熔融盐被用作显热蓄热材料)。由于相变蓄热技术上的复杂性,目前在太阳能热发电系统中没有实际应用,只停留在实验阶段。但是,相对于显热蓄热,相变蓄热具有相变潜热大、相变温区窄、可显著降低蓄热系统的尺寸等优点,是国内外许多学者研究的热点之一。
目前对相变蓄热器的研究主要集中在低温(温度不超过393 K),大多被应用在食品干燥和建筑采暖;而对高温相变蓄热器的研究较少。本文在前人研究的基础上,考虑传热流体的流动阻力、动量损失,建立了高温壳管式相变蓄热器的二维物理和数学模型,基于焓法模型,用有限容积法数值模拟相变材料的固—液界面位置和温度随时间的变化情况,分析流速和入口温度对熔化时间的影响,并对流速和入口温度进行了优化,讨论了相变材料温度、空气出口温度和固—液界面位置的变化规律。
1 相变材料的选择
高温相变材料是熔点温度范围一般在393 K~1 123 K的纯盐、混合盐、碱和金属,其中最有应用前景的是混合熔融盐。方明等人对多种熔融盐混合物应用在相变热存储设备中的效果进行了探讨,得出了混合熔融盐与单组份的熔融盐相比具有很大的优势,因此本文选用的相变材料为Li2CO3(20)—Na2CO3(60)—K2CO3(20) (质量百分比)的混合熔融盐。其热物性如下:密度2 380 Kg/m3;导热系数1.73 W/m·K (固)1.83 W/m·K(液);比热容1 590 J/kg·K(固)1 880 J/kg·K(液);熔解热283 kJ/kg;熔点823 K。
2 物理和数学模型的构建
2.1 材料密度与结构厚度设计
物理模型如图1所示,为单级壳管式结构。
管子的内径为20 mm,管长1.2 m,相变材料厚10 mm,被封装在壳体之间。传热流体(空气)在内径中流过与相变蓄热材料交换热量。
在保证数值模型计算结果准确性的前提下,做如下假设:
(1)忽略管壁的厚度,内壁设为传热壁,忽略管壁的热阻,并且外壁、侧管壁以及底部设为绝热壁;
(2)相变材料具有均质和各向同性。
边界条件和初始条件:入口条件设置为流速入口,管壁和管内空气的初始温度为298K。
2.2 注重流变性的物理计算模型
空气的质量守恒式(1)与动量守恒方程式(2)。
?ρ?t+?(ρu)?x=0(1)?ρ?t+?(ρu)?x=0(1)
?(ρu)?t+?(ρuu)?x=??x[?(p+23ρk)+τ]+Fi(2)?(ρu)?t+?(ρuu)?x=??x[-(p+23ρk)+τ]+Fi(2)
式中:ρ为空气的密度,u为空气流速,p为静压,k为湍流动能,τ为应力张量,其中:
Fi=0.3164Ra0.25L2Riu22g(3)Fi=0.3164Ra0.25L2Riu22g(3)
式(3)中:L为管长,Fi为空气流经管径时的沿程损失,Ra是雷诺数,Ri为圆管的内径。
对相变材料熔化时考虑对流传热,其能量方程为:
(ρcp)p?θ?t=??x(kp?θ?x)+1r??x(r?θ?r)?ρpΔH?f?t(ρcp)p?θ?t=??x(kp?θ?x)+1r??x(r?θ?r)-ρpΔΗ?f?t(4)
θ=T?Tm(5)?????f=0,0f1,f=1,θ0θ=0θ0(6)θ=Τ-Τm(5){f=0,θ00f1,θ=0f=1,θ0(6)
式中:Cp为相变材料的比热容,Kp为相变材料的导热系数,H表示熔解热,T、Tm分别代表相变材料的温度和熔点。f是相变材料中液相所占的比例,计算相变材料有效导热系数的公式如下:
kp=fke+(1?f)ks(7)kekl=CRan(8)kp=fke+(1-f)ks(7)kekl=CRan(8)
C值的大小依据传热流体的入口温度,取值如下:
C=?????0.24,0.18,0.16,θin≤1010θin≤2020θin(9)C={0.24,θin≤100.18,10θin≤200.16,20θin(9)
式(9)中:θin为传热流体入口温度与相变材料熔点温度之差,Ks表示相变材料在固态时的导热系数,Kl为相变材料在液态时的导热系数,Ke是相变材料的有效导热系数,在本文中C取0.16。
计算流体力学软件Fluent,可以直观的仿真出所需要的结果,将 GAMBIT软件划分的网格(将计算区域离散化)导入到Fluent中进行数值运
文档评论(0)