基于格子Boltzmann方法纳米流体流动与传热机理的研究.pdfVIP

基于格子 Boltzmann 方法的纳米流体流动与 传热机理研究 李强，宣益民 (南京理工大学动力学院 808 教研室，江苏 南京 210094) 摘 要：纳米流体是胶体悬浮物，它与普通固液两相流相比，具有不同的动力学和热力学性 质。本文建立了模拟纳米流体流动与传热的两相流格子 Boltzmann 模型，数值模拟结果符合 理论分析和实验结果。 关键词：格子 Boltzmann；纳米流体；两相流格子 Boltzmann 模型 1 引 言 纳米技术的发展为热科学和热能工程技术带来了新的挑战与机遇。1995 年，美国 Argonne 国家实验室的 Choi[1]等人提出一个崭新的概念—纳米流体：即以一定的方式和比 例在液体中添加纳米级金属或金属氧化物粒子，形成一类新的传热冷却工质。实验证明[1-3]， 在液体中添加纳米粒子可使流体的粘性系数、导热系数和对流换热系数均有显著增加，表现 出不同于普通两相流体的性质。 由于纳米粒子的小尺寸效应，悬浮的纳米粒子的行为更接近于流体，传统的两相流的方 法不适用于这种微小尺度下的不连续两相流动传热问题，而在格子气自动机模型基础上发展 起来的格子 Boltzmann 方法(简称 LB)则为解决此类微尺度两相流动传热问题提供了先进的 工具，可以用来研究纳米流体流动状态下的流动与传热特性在介观尺度上揭示纳米流体在运 动状态下的传热机理。 本文借鉴 J.M.Buick 和 C.A.Greated[4]在单相流中处理重力的方法及 S.Chen,X.Shan[5]提 出的两相流格子Boltzmann 模型，推导了适用于纳米流体的基于九点正方形模型的两相流格 子 Boltzmann 流动与传热模型，并据此研究了纳米流体的输运参数，从介观层次研究纳米流 体强化传热特性。 2 两相流格子 Boltzmann 流动模型 格子 Boltzmann 方法从离散的角度出发，通过与宏观的联系描述流体的热力学性质，具 有全离散、并行性强、无数值误差等优点，可以用于模拟复杂的流动与传热。格子 Boltzmann 方法在过去的几年里得到了很大的发展，已被成功地用于流体动力学的各个领域。我们简要 介绍一下基于九速四方格子的两相流格子 Boltzmann 模型。 在九速正方形的离散网格上的二维空间流场中，其速度向量的定义如下： ei (0,0), i 0 (i −1)π (i −1)π (cos , sin c i ( 1) ei 1 ) , 1,...4 2 2 (2i 1) (2i 1) − π − π 2 (cos , sin c i ei2 ) , 5,...8 4 4 各相具有速度为e 的粒子的分布密度定义为f σ (x, t) ，平衡分布函数为f iσeq (x, t) 。假 i i -1- 设各相流体之间没有内外力作用，则其格子 Boltzmann 演化方程可表示为[5]： σ σ 1 σ σ