复杂多孔介质腔体内自然对流换热的数值研讨.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：063091 复杂多子L介质腔体内自然对流 换热的数值研究 。曾敏王刚黄自鹏陈秋炀王秋旺 西安交通大学动力工程多相流国家重点实验室西安710049 Tel：029 E·mail：wangqw@mail．xjtu．edu．cn 摘要：本文通过数值模拟研究了由两个弯曲表面及两个水平直表面组成的复杂多孔介质腔体内的自然 对流换热问题。两个弯曲表面的形线分布采用了余弦函数，温度分别保持恒定，上下表面绝热。采用 有限容积方法及Brinkman扩展达西模型，在曲线坐标系下求解多孔介质内的流动换热问题，同时分析 了渗透率、孔隙率等参数对复杂腔体内流动和换热的影响。 关键词：多孔介质复杂腔体自然对流换热有限容积法 1引言 在化工、环境、机械、航空和核能等领域广泛存在着多孔介质内流动和传热问题， 这一问题已引起越来越多研究者的注意，在公开的文献中不断地有关于多孔介质内流动 和换热问题基础研究的报道【J】。在实际应用中，多孔介质内的流动很多均是由于热驱动 引起的，例如物料干燥、化学堆积反应床、热能储存、电子设备冷却和太阳能利用等。 迄今为止，对于含多孔介质规则方腔内自然对流问题研究较多瞳’41，而关于含多孔 介质复杂腔体内自然对流问题的研究却较少降¨，这些研究结果均表明，波纹表面形状对 于多孔介质腔体内自然对流换热有较大的影响。但是在这些文献中，均没有考虑多孑L介 质渗透率和孔隙率的影响，本文在曲线坐标系下，并考虑流体非达西效应，研究了Ra 数、渗透率、孔隙率及相关几何参数对多孔介质复杂腔体内自然对流换热的影响。 2物理问题及数值方法 如图l所示，腔体内充满各向同性、均质的饱和多孔介质，方腔上、下两壁面保持 水平，宽度为w且绝热，腔体的高度为H。左右壁面为弯曲壁面，满足余弦函数分布， 右壁面温度为Th，左壁面温度为T。(Th)。腔体内固体骨架与不可压缩牛顿流体的温度 相等，即采用局部热平衡假设；流体和固体骨架的热物性为常量，同时引入Boussinesq 假设来考虑自然对流效应；采用Brinkrnan扩展达西模型描述多孔介质腔体内的流动。 对多孔介质腔体内牛顿流体的二维层流流动换热问题，其无量纲控制方程如下： V．矿=o (1) ● jD ID． 2u一二二一u (2) —0(矿．vu)=一二，二一+Prv s‘、 OX D口 2矿一旦矿+RaPr0 (3) —l-f矿．vr)=一兰+prv 占‘’ 舀Y Da ．。89 。8e 8二e a二0 ，I、 u瓦坩万2萨+一Oy2 ¨’ 上述方程中涉及的无量纲量定义如下： (U，矿)=(11,V)W／a。，P=pW2／po口：，Pr=}，／口。， (x，】，)=(x，y)／W， 基金项目：国家自然科学基金资助项目 947 Da=x／w 2，Ra=妒(瓦一瓦)∥3／坳。，口=(T—To)／(瓦一瓦) 其中ro=(rh+7-,．)／2，盯卅是多孔介质的有效热扩散系数，其定义式为：