多孔介质内自然对流的不可压多松弛格子Boltzmann模型研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123568 多孔介质内自然对流的不可压多松弛格 子Boltzmann 模型研究 * 刘清，何雅玲，李庆，何超 （西安交通大学能源与动力工程学院热流科学与工程教育部重点实验室，陕西西安，710049） （Tel: 029 Email: yalinghe@mail.xjtu.edu.cn ） 摘要：建立了模拟多孔介质内自然对流换热的不可压多松弛格子 Boltzmann 模型。本文首先构造了一 个不可压单松弛格子Boltzmann 模型，该模型采用一个基于压力的分布函数求解速度场，而采用内能 分布函数求解温度场。在单松弛模型的基础上，构造了相应的不可压多松弛格子 Boltzmann 模型。利 用Chapman-Enskog 展开方法，可以准确地由该模型推导出多孔介质内不可压缩流动与换热的宏观控 Ra Da 制方程，在理论上消除了已有模型存在的压缩性误差。利用该模型对不同 数、 数和孔隙率下的 多孔介质封闭方腔内自然对流问题进行了数值模拟，模拟结果与相关文献结果吻合良好。 关键词：格子Boltzmann 模型；多松弛模型；不可压缩；多孔介质；自然对流 0 引言 多孔介质内的流动与换热是许多领域都涉及到的基础问题，在原油开采、电子元器 件冷却、核废料的储存、太阳能集热器、地热能源系统等工程实际中有着广泛的应用， 因而对其开展研究具有重要的现实意义。近几十年来，多孔介质内的自然对流换热问题 在理论和应用研究方面都取得了迅速发展，国内外许多学者运用传统数值方法，如有限 差分法、有限容积法、有限元法等对其进行了大量的模拟研究[1,2] 。格子Boltzmann 方法 作为一种有效的介观数值模拟方法，由于其边界条件易于实施，并行性好，在模拟具有 复杂几何结构的问题时具有较大优势，因而在这个方向得到了成功应用。 采用格子Boltzmann 方法对多孔介质内的流动与换热问题进行模拟主要在表征体元 (REV)尺度上进行[2] 。REV 尺度的格子Boltzmann 方法，无需介质的细节结构，计算效 率高，可用于大规模工程计算。基于Nithiarasu 等人[3]发展的多孔介质流动模型，Guo 等人[4]提出了用于求解多孔介质内流动的等温格子Boltzmann 模型，并且成功地将其推 [5] [6] 广到用于求解多孔介质内热流动的热格子Boltzmann 模型 。随后，Seta 等人 采用双 分布函数格子Boltzmann 方法，研究了多孔介质封闭方腔内的自然对流，验证了格子 Boltzmann 方法模拟多孔介质内热流动的有效性，他们发现在相同的网格尺寸下，格子 Boltzmann 方法的计算效率比传统的有限差分法要高。 但是已有的多孔介质格子Boltzmann 模型存在一些缺陷。一方面，流体在多孔介质 内的流动一般是不可压缩的，而已有的格子Boltzmann 模型所导出的宏观控制方程，从 形式上来说是可压缩的，因此模拟不可压缩流体在多孔介质内的流动时，必然会存在一 定的压缩性误差[5] [7] ，这种压缩性误差有可能使数值解偏离实际问题的真实解 ；另一方 面，已有的多孔介质格子Boltzmann 模型都是采用单松弛碰撞算子，在某些情况下，比 如低粘度下模型的稳定性会降低。为了克服已有模型存在的上述缺陷，本文将在双分布 函数模型的基础上采用多松弛格子Boltzmann 方法[8]构造一个用于多孔介质内热流动的 不可压多松弛格子Boltzmann 模型，该模型能够准确地导出不可压缩流体在多孔介质内