有限容积法与格子Boltzmann方法耦合模拟传热流动问题.PDFVIP

2010 年 第 55 卷 第 32 期：3128 ~ 3140 《中国科学》杂志社 论 文 SCIENCE CHINA PRESS 有限容积法与格子Boltzmann 方法耦合模拟传热 流动问题 栾辉宝, 徐辉, 陈黎, 陶文铨* 西安交通大学能源与动力工程学院, 西安 710049 * 联系人, E-mail: wqtao@ 2010-06-09 收稿, 2010-08-10 接受 国家自然科学基金重点项目资助 摘要 自然界和工程领域中的许多物理现象的发生通常涵盖几个数量级的几何空间及时 关键词 间范围, 我们将其统称为多尺度物理现象. 在模拟多尺度问题时, 如果仅采用宏观方法, 有限容积法 则会存在一些不足, 如无法预知微小部分的细节以及引入复杂的经验关联式; 如果仅采 格子Boltzmann 方法 用介观/微观方法, 则需要消耗大量的计算资源. 构造宏观-介观、宏观-微观、宏观-介观- 多尺度 耦合 微观等多种层次上方法的耦合体系, 可以在很大程度上克服这些不足. 构造了宏观有限 多孔介质 容积法(FVM)与介观格子Boltzmann 方法(LBM) 的耦合模型(CFVLBM), 给出了由宏观物 理量重构密度分布函数和温度分布函数的两个重构算子, 解决了LBM 与宏观方法耦合的 关键难题. 选取二维、三维典型传热流动问题对耦合模型进行了考核, 计算结果同基准解 符合得很好. 最后将 CFVLBM 应用于计算多孔介质内的复杂流动问题. 研究表明, 基于 文中重构算子的CFVLBM 可以准确有效地应用于模拟传热流动问题. 世界的本质是多尺度的[1]. 大千世界里, 物质本 类[3,4], 即过程多尺度问题与系统多尺度问题. 湍流 身是由原子构成的, 但是其几何外观往往会很大; 原 传热是典型的过程多尺度问题, 最大的涡可以与通 子的运动过程发生在飞秒(10−15 s)级别, 然而事物演 道尺寸相比拟, 而小尺度的涡比通道尺寸小数个数 变的过程又非常缓慢. 在生物和化学科学里, 分子尺 量级, 在湍流传热中不同尺度的涡相互作用构成了 度上的不同形态产生了生物体尺度上的复杂现象; 极其复杂的流动与传热过程; 质子交换膜燃料电池 在固体破坏中, 不同尺度的微损伤相互作用产生了 中发生的传递过程包括了反应气体从双极板通道(尺 更大尺度上裂纹, 进而导致材料破坏; 在流体力学 度为毫米- 厘米量级)穿过扩散层到达催化层(尺度为 中, 不同时空尺度的涡相互作用构成复杂的流动图 微米量级), 几何尺度经历了 3~4 个数量级的变化, 案[2]. 空间多尺度常常伴随着时间多尺度, 两者共同 在这里宏观尺度上的流动与微观尺度上的扩散、电化 构成了缤纷复杂的客观世界. 古人已经对多尺度现 学反应的耦合过程组成了一个典型的过程多尺度现 象有了很多的描述, 如“一尺之棰, 日取其半, 万世 象. 系统多尺度问题以数据处理中心的电子器件散 不竭”( 《庄子·天下》), “夫自细视大者不尽, 自大视 热问题为例, 从整个机房冷却气流的合理组织到机 细者不明”( 《庄子·秋水》), “天街小雨润如酥, 草色 箱、印刷板、封装以及芯片的冷却, 几何尺度经历了 遥看近却无”(韩愈《早春呈水部张十八员外》). 与热 近十个数量级的变化. 在过程多尺度中, 微观、介观 量传递过程有关的多尺度问题基本上可