基于格子Botlzmann方法的三维多孔介质中多相导热过程的.PDFVIP

中国科学 E 辑 : 技术科学 2009 年 第 39 卷 第 10 期 : 1743 ~ 1750 SCIENCE IN CHINA PRESS 基于格子Botlzmann 方法的三维多孔介质中多相 导热过程的研究 * 赵凯, 李强, 宣益民 南京理工大学动力工程学院, 南京 210094 * E-mail: liqiang@ 收稿日期: 2008-07-01; 接受日期: 2008-11-20 国家自然科学基金项目(批准号:和新世纪优秀人才支持计划资助 摘要 建立了模拟三维多孔介质中多相导热过程的格子Boltzmann模型, 结合三维多孔介 关键词 质重构算法, 从孔隙尺度对CPL/LHP 毛细芯中的瞬态导热过程进行了模拟, 计算了不同条件 多孔介质 下毛细芯内温度分布, 分析了热负荷、孔隙率等因素对导热过程的影响, 并计算了三维多孔 格子Boltzmann 介质的有效导热系数. 孔隙尺度 有效导热系数 多孔介质中能量、动量和质量的传递现象遍及于 介质内的导热问题, 而最近出现的格子Boltzmann 方 工业生产的各个领域, 例如: 石油开采, 食品、药材 法可以自动跟踪界面, 从而可以从孔隙尺度精细地 和农产品对流干燥, 各种工业换热设备中的多孔材 描述多孔介质内的能量传递过程. 因此, 本文建立了 料里的蒸发以及毛细抽吸过程等[1,2], 对其内部能量 描述三维多孔介质中多相导热过程的格子Boltzmann 传递过程的深入研究有着十分重要的意义. 以往, 关 模型, 以毛细泵回路和环路热管(CPL/LHP) 的启动过 于多孔介质内部的能量传递过程的研究存在以下不 程为研究背景, 从孔隙尺度对多孔毛细芯内的瞬态 足[3~6] : 1) 大多是建立在连续性假设之上, 以有限差 导热过程进行了模拟, 分析了不同热负荷、孔隙率等 分、有限体积、有限元等为离散方法, 基于N-S 方程 因素对其达到相变点时间的影响, 同时利用本文提 的宏观尺度的数值模型, 该类模型忽略了多孔介质 出的模型对各种多孔介质的有效导热系数进行了计 内部骨架和孔隙的空间拓扑结构对能量传递过程的 算. 影响[7,8]; 2) 在传统CFD 中对复杂结构的有效导热系 1 基本模型 数选取都是根据经验公式, 实际上, 多孔介质的有效 导热系数不仅与骨架和孔隙的份额有关, 而且与其 1.1 控制方程 空间结构有着密切联系. 仅以份额来决定有效导热 描述无热源的两相(骨架和流体) 导热问题的控 [9,10] 系数显然不够准确 . 制方程组如下: 由于多孔介质内部具有很多随机的骨架与流体 ⎛ ∂T ⎞ 2 (ρ cp ) ⎜ ⎟ k =∇ T , (1)