二维球状蓄热体蓄放热性能数值模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093017 二维球状蓄热体蓄放热性能数值模拟研究 1 1 2 2 崔海亭, 吕水燕，邢玉明，吴斌 （1. 河北科技大学机械与电子工程学院, 石家庄 050018） （2. 北京航空航天大学航空科学与工程学院, 北京 100083) （Tel：0311Email:cuiht@126.com） 摘要：本研究采用二维导热微分方程描述球状蓄热体区域，根据数值计算理论及相变储能理论对二维 球状蓄热体的蓄放热过程进行了模拟计算。本研究计算中采用分段线性拟合的方法模拟相变材料复杂 的比热—温度的关系，配合热容法求解相变导热问题。根据研究计算，热容法具有编程方便、便于求 解变物性问题的突出优点。本研究得出了二维球状蓄热体蓄放热性能曲线，为蓄热槽的设计、优化、 评价等提供了参考依据。 关键词：相变材料；二维球状蓄热模型；数值模拟；热容法；蓄放热特性 0. 引言 太阳能具有储量丰富、洁净、可再生等优点，但太阳能资源自身的间歇性与建筑物 的热量需求不相匹配，因此热量储存成为任何太阳能直接供热系统中必不可少的组成部 分。蓄热技术是提高能源利用效率和保护环境的重要技术，采用蓄热技术可以调节能量 供需，利用相变潜热蓄热的物质常称为相变材料(Phase Change Material，PCM) 。由于物 质相变时温度基本不变，而且相变潜热很大，故蓄热密度很高，蓄热设备的尺寸和重量 都比较小，加上物质相变时吸、放热效率很高，故整个蓄热系统具有很高的效率，因而 这种蓄热方式近年来在国际上受到极大的关注，成为目前研究热点[1-7]。因此发展高效 相变蓄热储能技术将带来显著的经济效益，环境效益和社会效益。 1. 模型描述 1.1 物理模型 本研究所采用的模型为直径100mm 的不锈钢金属球体内封装相变温度为58℃的相 变材料（PCM 为固液相变）。金属球体放置在直径为600mm 高度为 1200mm 的圆柱形 桶体内，每层均匀放置19 个，共11 层。 为简化计算做如下假设： 1) 蓄热体之间的相互影响不考虑； 2) 忽略PCM 发生相变时的体积变化，认为密度恒定； 3) 球状PCM 内部的温度分布只与半径及纵截面上的角度有关； 4) 在PCM 内部不考虑对流换热； 1.2 数学模型 由于本研究所用的相变蓄热球体直径较大，故采用二维极坐标下导热微分控制方 资助项目：国家自然科学基金（No），河北省自然科学基金资助项目（No.E2008000700 ） 程： ∂T 1 ∂ ∂T 1 ∂ k ∂T ρc (rk )=+ ( ) ∂t r ∂r ∂r r ∂ r ∂ θ θ 求解相变导热问题的数值解法，常用的有固定步长法，变时间步长法，等温面移动 法、自变量变换法，这些都属于温度法的 范畴，需要分区求解，此外还有焓法和热 容法。 δ n ) 本研究采用热容法是把相变导热问题 r x δ 看成是比热c 随温度T 改变的