![a177dd839250903254e69ea54cb54adc.png](https://i-blog.csdnimg.cn/blog_migrate/d0d0b54d659ccc89f19e1e349f379603.png)
在大量实验研究的基础上,本文同时采用数值模拟的方法对基于微热管阵列的太阳能空气集热器进行仿真模拟研究,分析计算域内部的流动和传热特性,通过对模型和边界条件的优化和改变,计算不同数值模型,获得与实际最吻合的模拟效果,为空气集热器集热器的性能及结构优化工作提供一定的依据。
![a177dd839250903254e69ea54cb54adc.png](https://i-blog.csdnimg.cn/blog_migrate/d0d0b54d659ccc89f19e1e349f379603.png)
以微热管阵列真空管式太阳能集热器单元为研究对象,建立集热器的 CFD模型。模拟过程中考虑辐射换热作用,吸热膜、玻璃盖板之间的对流和辐射换热作用,风道内空气与冷凝段翅片之间的对流换热以及微热管阵列蒸发段与保温层之间的热传导等传热过程。如图1所示,建立以真空管的简化底平面的圆心为坐标原点的笛卡尔坐标系,建立集热单元三维,稳态计算模型。在数据模拟过程中,取微热管阵列的表观导热系数为105W/m·K。
该模型基本假定如下:
1) 流体不可压缩
2) 认为真空管管壁很薄,忽略热阻
3) 换热与流动充分发展。
4) 各种材料的热物性见表1 。
![2509ce152a94d480dec00e60a3c8d3fd.png](https://i-blog.csdnimg.cn/blog_migrate/5fdfe33c8e04e83386a4a329cc0b7778.jpeg)
图1微热管阵列-真空管式太阳能空气集热器集热单元物理模型
![570db57351631087101d8c2adb2f7bcb.png](https://i-blog.csdnimg.cn/blog_migrate/d0bc3248e1af896300b331f655034198.jpeg)
表1 材料的热物理性质
集热器的整个换热过程遵守质量守恒,能量守恒及动量守恒。集热器内的换热包括:真空管内部空气与微热管蒸发段的流动与辐射换热,风道内空气与冷凝段翅片的强迫对流流动。真空管封闭空腔内空气的自然对流符合Boussinesq假设,计算域内流体流动及传热满足以下基本控制方程,为了便于描述,集热单元内空气流动与换热的三维控制方程采用守恒方程的统一形式,控制方程中各变量φ对应的Γ为扩散系数和Sφ为源项。
![1e898f6899879e471884df6ac5132161.png](https://i-blog.csdnimg.cn/blog_migrate/9f6997f0e5580b031a76fdca71ad48a8.jpeg)
计算模型中的边界条件设置见下表。
![e6fc7fe8cef3e2078df3a5d661448e10.png](https://i-blog.csdnimg.cn/blog_migrate/03f8c007717417b01ee059a516602fe8.jpeg)
![c1737baf07bd94b760958752a7081620.png](https://i-blog.csdnimg.cn/blog_migrate/97def0162ae87d136272fb3c7d117006.jpeg)
式中,ux—x轴向速度, v —径向速度,w —轴向速度;kin —入口处流体的湍流
动能, εin —入口湍流动能的耗散,CD湍流动能常数。
本数值模拟采用Gambit前处理软件进行模型的建立和网格划分。该太阳能空气集热器,真空管长度为1800 mm,内径为47mm,可见在以真空管方向为参考的轴向和径向尺寸差别较大,且风道属于更大尺寸的实体模型,因此计算模型的网格划分采用分体,分区域的方法进行划分(疏密结合)。为减少网格数量,提高计算速度,将风道沿x方向分割为3部分分别进行划分,两端部分采用尺寸为20mm的大网格,中间部分则为尺寸为3、1.5和1mm的小网格。网格部分的划分效果图见下图。选用模型网格数量为121万的划分模型。
![22ceb74fcfe4d9644e6124d466dde573.png](https://i-blog.csdnimg.cn/blog_migrate/ccd587e78299f536961288ebed321c60.jpeg)
模型求解及验证:
如下图所示,横坐标表示太阳辐照度为750W/m2时,不同空气质量流速下的4个工况点的平均集热效率值。1,2,3,4工况对应的空气流速值分别为0.6,0.8,1.2,1.4 m/s。模拟结果风道出口温度整体水平要高于实测数据,这是因为单个集热单元与空气的换热条件要优于15个单元总体的换热条件。15个单元与空气的换热受到前面单元的影响,且随着换热的深入,空气温度逐渐升高,空气与翅片间的温差减小,使得总体换热程度低于15个集热单元的代数和,集热效率同时也高于实验实测值。图5-4为瞬时效率拟合曲线,实验与数值计算误差在7.8%。因此认为该模型精度在工程允许范围内,可用于该类型集热器的数值模拟研究。
![62f47b0ed1aff76ec91d4285385a18c6.png](https://i-blog.csdnimg.cn/blog_migrate/ed63409605b72f5466c3aa01c0da6484.jpeg)
![1476b70b447236b6505fb3b54367d84c.png](https://i-blog.csdnimg.cn/blog_migrate/152c3bd5ee38ca84bfa50c70d2bc1d18.jpeg)
![976e04e49fcecb2d8aa4f1721ba7c77a.png](https://i-blog.csdnimg.cn/blog_migrate/e5aa4ff0893375fb9b9be2a7de6502a9.jpeg)
![d5588a8b03f6cb6772995cd0a2cddc10.png](https://i-blog.csdnimg.cn/blog_migrate/9dd787737ea61c8d64565a5f05c6e292.jpeg)