纳米流体多相流动的多尺度模拟方法 Multiscale Simulation of Nanofluid Multiphase Flows.pdf

CHINESE COM物PUTATIO理NALPHYSICS 第200鬻9年月1嚣1期 JO计URNAL算OF 。赫。．’猫 訾j镒 【文章编号] 1001．246x(2∞9)06．0849．08 纳米流体多相流动的多尺度模拟方法 周陆军， 宣益民， 李 强 (南京理工大学动力工程学院，江苏南京210094) [摘 要] 针对纳米流体多相流动的微观特征，提出一种基于格子Boltzmann的多尺度耦合方法，在速度和悬浮 纳米粒子分布变化比较剧烈的区域采用细网格多相模型，在其它区域视纳米流体为均匀混合的单相流体，使用 粗网格单相模型．为保证不同尺度区域之间的物理信息(参数)的准确传递，运用质量和动量守恒原理，建立 跨区域的边界耦合模型．几个算例表明，该方法既可以反映纳米流体流动的微观特征，又能提高计算效率，与 单纯使用多相模型相比，节省大量时间． [关键词] 格子Boltzmann；多尺度模拟；纳米流体；多相流 [中图分类号】0359 【文献标识码] A O 引言 宏观的桥梁，研究其在不同尺度中的运用对复杂流动问题的处理十分有意义．LB方法最早使用的均匀网格 对于大部分工程问题已经足够，然而此方法在处理某些复杂流动(如机翼尾流、多相流、台阶流等)问题时，受 到精度和效率的限制．随后国内外很多学者对LB方法进行改进，提出了基于有限差分、有限容积、有限元等 思想的网格处理技术¨’2]，如多重网格、无网格方法，以及常用的空间插值方法b圳和区域分类法伸“州等，其 中空间插值方法的特点是网格在整个计算空间是变化的，迁移得到的分布不在格点上，格点上的分布是对迁 移后的分布插值求得的，其计算精度和效率与插值方法的选取有关；区域分类法是不同区域采用不同的网 格，在同一区域采用统一的计算网格，其计算精度和效率与区域网格划分尺度及区域之间的耦合相关． 上述LB多尺度方法均针对单相流体，很少在多相流动中应用．在模拟多相流动时，由于考虑各相之间 复杂的作用，严重影响了计算效率；同时为反映多相流体的微观信息，网格尺度划分较小，导致多相模型只能 模拟比较小的空间．实际工程问题的计算区域往往比较大，使用多相模型在计算时间上不再合适，如果采用 单相模型就不能反映微观信息．本文针对多相流动问题的特点，根据单相多尺度耦合方法的思想，提出了一 种新的LB多尺度耦合模型．将计算区域分为粗细两种网格尺度，在粗网格区域采用单相模型，细网格区域 采用多相模型．本模型假设多相流体在大网格尺度上 密度均匀分布，可作为单相处理．如图l所示，在网格 尺度100rim时各相密度分布是波动的，要作为两相处 000 理；当网格尺度到1 nm时，可认为每个网格上的密 度是统一的，可作为单相流体计算．此假设与能量最小 多尺度模型类似u¨，区别是，本文是基于LB方法的耦 合．通常在速度、密度等变化剧烈的区域采用细网格多 相模型，在其他区域采用粗网格单相模型．该方法在保 证较高的计算效率与精度的同时，又能刻画多相区域 图1纳米流体在不同网格尺度下的分布 内的微观信息，为解决某些大尺度甚至三维空间的多 of Fig．1Distributionnanofluidonmultiscale 相流动问题开辟了新的思路． [收稿13期】2008—05—19；[修回日期】2008—08—12 [基金项目】国家自然科学基金重点(No资助项目 【作者简介】周陆军(1982一)，男，山东潍坊，博士生，主要从事纳米流体多尺度研究． 万方数据万方数据 850 计 算 物 理 第26卷 1 模型介绍 LB最早使用的是单相模型，然后推广到多相模型．以下对LB模型作基本介绍． 1．1 LB单相模型