孔隙率阶梯分布多孔介质内自然对流传热与流动.docVIP

孔隙率阶梯分布多孔介质内自然对流传热与流动研究 杨伟a，赵柄翔a，薛思瀚a 辽宁工程技术大学a.建筑工程学院，辽宁省，阜新市，123000 摘要：以孔隙率阶梯分布多孔介质内传热与流动为研究对象。采用局部热平衡假设，考虑流体密度随温度变化，引入Brinkman-Forchheimer的扩展Darcy模型进行修正，建立水平方向孔隙率阶梯分布多孔介质自然对流数值模型，采用有限体积法进行计算。利用自主研制的实验台对数值模型进行实验验证，分析出孔隙率的不同阶梯分布情况对传热与流动的影响。综合数值计算和实验结论表明：相同温差条件下，大孔隙率侧为高温壁面对流传热达到稳态的时间要小于小孔隙率侧为高温壁面时对流传热达到稳态的时间；腔体内最大温度出现在高温壁面处，靠近高温壁面的孔隙率越大，高温部分体积越大，随着Ra增大，温度等值线由平缓顺时针曲线变为近似垂直于高低温壁面；流体速度场分布主要受孔隙率分布影响，最大流速随温差和Ra增大而增大，且最大流速出现在大孔隙率靠近壁面处, 与是否靠近高、低温壁面无关；高温壁面随着Ra增大而增大，当lgRa4时， Ra的增加对的变化影响较小。 关键词：孔隙率；阶梯分布；多孔介质；自然对流；传热；流动 中图分类号：TK124 文献标志码：A 0 引言 多孔介质内流动和传热作为流体力学和传热领域的重要现象，在自然界及工农业生产中广泛存在。多孔介质集热墙、多孔介质太阳能吸热器、多孔介质保温材料、大规模集成电路设备的放热等都包含多孔介质内流动与传热现象。因此对内部流动传热过程的深入研究有着十分重要的意义。[1]对底部恒温加热水平多孔介质结构进行了数值计算与实验，实验发现有与数值计算结果不同的对流现象，说明其数值模型有待提高。陈凯[2]等通过在方腔内布置固体颗粒模拟多孔介质结构，并采用虚拟区域方法求解多孔介质中的流场和温度场。于立章[3]等对紊流及紊流过渡区范围内骨架发热多孔介质竖直通道内的非达西强制对流换热进行了数值模拟，发现了孔隙有效雷诺数、表面热流密度和冷却剂入口温度的变化对多孔介质流道内流动阻力及换热特性的影响。张志军[4]等提出了在平行平板间填充颗粒的同时，运用板面加开微型纵槽构成复合多孔介质，以改善流体流动与传热综合性能。而对于多孔介质内孔隙大小的不同分布情况，进行对流传热与流动的研究较少。 孔隙率阶梯分布多孔介质在自然界和生活中比较常见，如多层保温墙体、不同类型岩体、多孔介质填充层等，本研究建立了孔隙率阶梯分布的多孔介质模型，对其自然对流下的流动与传热进行了实验研究与数值模拟。详细分析了孔隙率阶梯分布下，孔隙率大小和雷诺数Ra对多孔介质内自然对流的传热与流动影响，及雷诺数Ra对高温壁面平均努谢尔数的影响。本文实验模型对于孔隙率为30%-50%的土壤、岩石、多孔介质墙体及保温材料有很好的代表性，对于多孔介质内温度控制、强化传热、削弱传热提供研究基础。 1数学模型 1.1 数学模型及数值计算 物理模型见图1。左表面温度为高温壁面Th、右表面温度为低温壁面Tc [5] ，上下前后壁面为绝热边界条件，假设局部热平衡；流体密度随温度变化[6]，其他热物性为常量，在模型壁面上速度均采用无滑移边界条件[7]。 图1 物理模型 Fig.1 Physical Model 对充满流体的多孔介质控制方程引用Brinkman-Forchheimer extended Darcy模型[8]进行修正。 连续性方程、动量方程、能量方程为: （1） （2） （3） （4） 上述方程定义为：为密度, kg?m-3；、分别为x、y?s-1；为时间，s；为孔隙率；为有效粘度，pa? s；为多孔介质渗透率，, m2；为惯性系数, ；为重力加速度，m?s-2；为热膨胀系数，k-1,；为温度, k；为低温壁面温度, k；为定压比热容，J?kg?K；为有效导热系数,为导热系数，Wm-1K-1, ，下标s为固体，下标f为流体。 为了揭示固体和流体传热的普遍规律，将方程（1）～（4）转换为无量纲方程。引入无量纲量定义为：无量纲坐标，，无量纲速度，，无量纲压力，无量纲密度无量纲时间，；无量纲准则达西数，瑞利数，无量纲温度,为任一点温度值,k。， ，，。可得到如下的无量纲控制方程： （5） （6） （7） （8） 无量纲边界条件为： 时，，。 时，，，，；，，；，，。 2实验研究 2.1 实验装置和内容 实验装置如图2所示，两侧为高低温水箱、中间为多孔介质方腔体（长0.25m，宽0.2m，高0.13m）；多孔介质方腔体尺寸同图1物理模型。水箱和方腔之间为玻璃薄层，厚度为0.8cm，水箱和方腔周围用保温材料包围