多孔介质对封闭腔体内传热传质的影响研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093246 多孔介质对封闭腔体内传热传质的影响 1,2 2 2 刘芳 ，陈宝明 ，王丽 1 2 （山东大学，山东 济南 250061； 山东建筑大学，山东 济南 250101） Tel：05318179815，E-mail：fliu@sdjzu.edu.cn 摘要：本文对部分填充多孔介质复合腔体内的自然对流传热传质现象进行研究，考虑复合腔体内多孔 介质与流体交界面上的应力滑移，利用达西扩展方程和N-S 方程及能量和组分传递方程分别对多孔介 质区域和纯流体区域建立数学模型。采用弱约束耦合多孔介质和流体交界面处的应力滑移，利用有限 元法进行数值求解，数值模拟结果表明多孔介质特性：厚度、渗透率、孔隙率对自然对流复合腔体内 传热传质存在不可忽略的影响。 关键词：多孔介质；自然对流；复合腔体；传热传质 0 前 言 多孔介质的流体流动和传热传质过程广泛存在于自然界和工程实际中，其中部分充 填多孔介质的复合腔体问题具有广泛的应用背景。如：干燥过程、太阳能集热器、热管、 核能利用、燃料电池及原油提取和地热工程等。由于流体在不同介质内流动，机理复杂， 因此对部分填充多孔介质腔体内中间交界面处的速度梯度的连续性存在不同的观点。 Beavers 和 Joseph (1967) [1]通过实验研究发现在多孔介质和流体交界面上存在速度滑 移，并提出了速度滑移表达式，但其中的滑移系数需要根据实验结果来确定。尽管文献 [1]的工作具有开创性的意义，但由于未能明确地阐述界面处的动量传递机理，因此当时 速度滑移模型并未得到足够重视。随着对复合区域研究的不断深入，许多研究者对曾经 普遍采用的界面连续性条件提出了质疑，开始对非连续性条件进行分析，文献[1]提出的 速度滑移条件再次引起了研究者的关注，其实验结果也被作为数值模拟的验证。但通过 理论分析认为界面处的速度、温度等物理量应为连续性变化， Ochoa-Tapia 和Whitaker [2,3] （1995a） 利用修正Brinkman 方程开创性地提出了应力滑移模型，并得到了许多研究 者的应用。D.B.Das 等 (2002)[4]利用有限容积法对三维系统流动进行了分析；C. Deng, D.M. Martinez(2005)[5]对底部为多孔介质的平板通道内的粘性流动进行了分析；Hua Tan, Krishna M. Pillai(2009)[6]利用了有限元法对有无应力滑移的斜向流进行了分析；Peng Yu 等 (2009)[7]对微生物反应器中的氧气输运过程的流动及传质进行了分析；Manoranjan Mishra, A. Ramachandra Rao(2005)[8]对流体在四周为多孔介质的二维通道内的流动进行 了研究，以上对于层流问题的研究，研究者均考虑了界面效应的影响，并应用了文献[14] 提出的应力滑移条件。巴西学者M.J.S. de Lemos[9-14] 的科研团队多年来对多孔介质中的 紊流流动展开研究。在已有的研究中，多是对于通道中的强迫对流流动的研究，但涉及 到封闭腔体内自然对流流动的研究较少。 . 因此，本文以污染物气体从多孔介质内部散发至室内空间为背景，采用两区域模型 法建立数学模型，考虑到中间交界面处剪切应力的不连续性，利用有限元法耦合弱约束 通讯作者：陈宝明，男，教授，博士生导师，主要研究方向：复杂体系中的传热传质。 基金项目：国家自然基金（No） 方程进行数值求解，对部分填充的多孔介质对整个腔体内的传热传质影响进行研究。 1 数学模型的建立 将多孔壁面－室内空间物理模型简化为如图 1 所示的二维空间，其高宽比为 1。左 侧为多孔介质区域，右侧为纯流体区域，在水平方向上建立温度、浓度梯度和湿度梯度。 假设所有外壁面均为不可渗透表面，无速度滑移，水平表面均为绝热表面。左侧多孔介 质为各向同