基于mcr模拟的高温高压太阳能吸热器传热过程研究.docxVIP

基于mcr模拟的高温高压太阳能吸热器传热过程研究 0 技术难点分析 近年来，该算法在加热能源中得到了迅速发展，大量实验和运营数据充分证明了该算法的技术可行性和商业应用前景。在塔式太阳能热发电系统中,太阳能吸热器是实现塔式太阳能热发电最为关键的核心技术之一。而有压腔式吸热器具有换热效率高的优点,是未来发展的方向,但也存在着一些技术难点,如:耐高温吸热器选材及制造工艺、耐高温石英玻璃窗口制造工艺难点以及吸热器腔体内气体流动与换热及其停留时间与换热面积的确定。前两个方面有赖于材料学科及制造业的发展,而第三个问题需要本领域科研人员在吸热器结构型式和换热机理等方面进行深入研究。 高温高压有压腔式吸热器内聚光传热过程是一个极其复杂的光能聚集、转换以及复杂耦合传热的过程。为揭示其光热转换机理、聚光传热特性以及流动换热性能,本文在西班牙PSA REFOS项目所研究的模块化有压腔式吸热器基础上,建立有压腔式太阳能吸热器轴对称计算模型,提出MCRT方法及其与FVM耦合方法,进而来揭示高温高压腔式太阳能吸热器内复杂耦合换热过程与换热机理,为太阳能吸热器的优化设计奠定一定的基础。 1 模型和方法的计算 1.1 陶瓷多孔介质模型 本文以西班牙PSA REFOS项目所研究的模块化有压腔式吸热器为研究对象,其结构示意简图如图1所示,包括二次聚光器与有压腔式吸热器两部分。为便于计算,将有压式吸热器内陶瓷多孔介质形状近似为抛物线型、将二次聚光器进出口多边形假定为圆形。即假定该吸热器为完全对称的旋转体结构,适合于采用轴对称模型计算方法。简化后的计算模型如图2所示。吸热器主要由椭圆形石英玻璃管窗、抛物线型陶瓷多孔介质、腔体空间、腔体壁面以及绝热管壁等组成。为使表达简洁,图2中将围成吸热器腔体的各壁面进行了编号,Wall 001-Wall 005分别为吸热器腔体左侧壁面、玻璃窗外表面、腔体右侧壁面、出口及延长管道壁面以及玻璃窗内表面。 1.2 有压腔式吸热器内对流换过程模拟控制方程 多孔介质接收器内太阳能热流密度分布采用MCRT方法自编程序进行模拟计算。在此之前,作者已运用MCRT方法研究了槽式太阳能吸热器聚光集热传热特性。计算中同时采用直角坐标系、柱坐标系及球坐标系,且假设光束垂直于二次聚光器进口平面入射、二次聚光器为镜面反射,考虑太阳光不平行度。MCRT模拟计算多孔介质接收器内太阳能热流密度分布程序框图见图3左侧MCRT方框。 获得接收器内热流密度分布之后,可将其与计算流体与传热的有限容积FVM法(FLUENT软件)进行耦合(参见图3 COUPLING与FVM方框),以研究整个有压腔式吸热器内光热转换机理、聚光传热特性以及流动换热性能。这包括图2中所示:①陶瓷多孔介质内的光热转换过程;②陶瓷多孔介质内的对流换热过程;③吸热器腔壁间的热辐射换热过程;④石英玻璃窗外大空间自然对流(或混合对流)换热过程;⑤石英玻璃窗向外热辐射换热过程。 对于各个过程的模拟,实现方法如下:对于①,采用源项法,将由MCRT计算所统计得到的太阳能热流密度分布通过跨接程序赋值给陶瓷多孔介质区域;对于②,通过求解修正的多孔介质区域轴对称坐标系控制方程即FLUENT内部多孔介质模型来模拟多孔介质内的对流换热过程;对于③,采用数值方法即FLUENT内部S2S辐射模型计算多表面系统辐射换热;对于④,先通过数值模拟方法计算出玻璃窗外大空间自然对流换热系数,然后作为本计算模型之边界条件;对于⑤,在玻璃窗外表面采用四次方定律来计算其对外热辐射。 综上所述,有压腔式吸热器内对流换热数值模拟控制方程包括轴对称坐标系下的连续性方程、动量方程、带源项能量方程以及RNG k-ε两方程,多孔介质区域为修正的轴对称坐标系控制方程。各边界条件为:入口采用质量流量入口边界条件,质量流量为1.98 kg.s-1、入口温度为894.15 K;石英玻璃窗外表面设置自然对流换热系数及四次方热辐射边界条件;出口及其延长管道壁面为保温绝热边界条件;出口采用充分发展假设,为保证假设正确性将出口段延长了5倍管道直径长度(参见图2);壁面上采取壁面函数法。 计算中对计算网格进行了网格独立性考核,使得计算结果与计算网格无关。 2 热流密度分布 自编采用蒙特卡罗射线追踪法(MCRT)求热流分布的计算程序,在前述适当假设前提下获得了PSA在1999年9月24日13:04试验工况下多孔介质接收器内太阳能热流分布,如图4所示。 由图4可见,太阳热流密度在多孔介质接收器内沿径向Y方向是对称的,而沿对称轴X方向呈连续性降低趋势,且多孔介质区域内右侧热流密度梯度变化较大。这是因为太阳光在多孔介质内是沿X方向被吸收的,且在前段较薄多孔介质层内,吸收了大量太阳光使得热流密度较高,而越到后段未被前段吸收而通过的光线越少则热流密度越低。由图4还可见,在多孔介质中