多孔介质方腔内混合对流格子Boltzmann模拟.pdf

第３４卷 第１期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３４，Ｎｏ．１ ２０１７年１月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｊａｎ． ，２０１７ 文章编号：１００１⁃２４６Ｘ（２０１７）０１⁃００３９⁃０８ 多孔介质方腔内混合对流格子Ｂｏｌｔｚｍａｎｎ模拟 １ ２ １ １，∗ １ 王婷婷 ，　 高　 强 ，　 陈　 建 ，　 徐洪涛 ，　 杨　 茉 （１．上海理工大学能源与动力工程学院，上海　 ２０００９３；２．杭州三花微通道换热器有限公司，杭州　 ３１００００） 摘　 要：采用格子Ｂｏｌｔｚｍａｎｎ方法，对多孔介质方腔内的混合对流现象进行研究．方腔内部中心有一发热圆，径宽比 Ｄ／ Ｌ＝０４，冷流体从方腔左下角入口流进，从方腔左上角出口流出，四周壁面绝热．在普朗特数Ｐｒ＝０７１和格拉晓 ４ 夫数Ｇｒ＝１４×１０ 时，分析理查德森数Ｒｉ和达西数Ｄａ对发热圆表面平均努赛尔数Ｎｕ的影响．结果表明：Ｒｉ数位 －３ －５ －２ 于１０ ～１０范围内，Ｎｕ 随Ｒｉ 的增大而减小． Ｄａ越大，Ｒｉ对Ｎｕ的影响越显著；Ｄａ数位于１０ ～１０ 范围内，强 制对流占主导的情况下（Ｒｉ≤０１），Ｎｕ 随着多孔介质的Ｄａ 的增大而增大． 自然对流占主导的情况下（Ｒｉ＝１０）， Ｎｕ对Ｄａ 的变化不敏感． 关键词：多孔介质；方腔；混合对流；格子Ｂｏｌｔｚｍａｎｎ模拟 中图分类号：ＴＫ１２４ 文献标志码：Ａ ０　 引言 多孔介质中的流动与换热问题广泛存在于各种工程问题中，典型的如农业生产中农作物的干燥过程、工 业锅炉中的蓄热与换热、生物工程中的生物过滤、以及新能源技术中的多孔介质燃烧等，因此，对多孔介质内 流动与换热的研究具有重要的意义．如今，采用数值模拟的手段对多孔介质内部流动与换热进行研究已被广 ［１］ 泛采用，其模拟分析方法主要集中在有限容积、有限差分、有限元等传统宏观方法上，如刘云龙 采用压力 －速度耦合的ＳＩＭＰＬＥ算法对方腔内混合对流换热进行模拟，探讨了混合对流向受迫对流过渡的临界参数． ［２］ Ａｌ⁃Ａｍｉｒｉ 采用有限体积法研究了顶盖驱动的二维多孔介质方腔内，热分层稳定的饱和水的流动特性，分析 ［３－４］ 了不同特性参数对高温壁面上平均努赛尔数Ｎｕ的影响．Ｋｈａｎａｆｅｒ等 对单个二维竖直开口方腔内的混合 ［５］ 对流换热问题进行模拟，分析了不同特性参数对等温线和流线分布及Ｎｕ的影响．孙猛等 研究竖直外界来 流方向的三层开口方腔内的混合对流换热问题，分析了不同特性参数对混合对流特性的影响． 与传统的模拟方法不同，格子Ｂｏｌｔｚｍａｎｎ方法基于分子动理论，具有清晰的物理背景，该方法在宏观上 是离散的，在微观上是连续的，因而被称为介观方法．格子Ｂｏｌｔｚｍａｎｎ方法作为一种新兴的数值模拟方法，因 其计算效率较高、演化过程简单、编程容易、具有良好的并行性、以及边界条件处理简单等优点，被广泛应用 ［６］ 于微尺度流动与换热、多孔介质、晶体生长等一些传统方法难以胜任的领域 ．近年来，部分学者使用格子