声悬浮过程格子Boltzmann方法研究.PDFVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 16 (2014) 164301 声悬浮过程的格子Boltzmann方法研究∗ 解文军 滕鹏飞 (西北工业大学空间应用物理与化学教育部重点实验室, 西安 710072) ( 2014 年1 月17 日收到; 2014 年3 月13 日收到修改稿) 采用轴对称多弛豫时间格子Boltzmann (LB) 方法, 研究了圆柱形封闭谐振腔中圆盘形样品的声悬浮 过程. 模拟结果表明, (001) 模式下谐振腔的共振长度 , 在谐振腔中心引入样品后共振漂移量 , 这与线性声学理论计算结果基本相符. 声悬浮力的LB 模拟过程包含了黏滞性效应和共振漂移效 应, 所获得的模拟结果与理论公式计算值在量值上一致, 而且其在细节上更符合实验现象. 此外, LB 模拟还 揭示出了声悬浮过程中的声压波形畸变、声流和声辐射压等非线性声学效应. 关键词: 声悬浮, 格子Boltzmann 方法, 非线性效应 PACS: 43.25.Uv, 47.11.–j DOI: 10.7498/aps.63.164301 域[20−22] , 在模拟非线性声学现象方面也取得了一 1 引 言 些进展. Buick 等[23] 采用LB 方法, 研究了大振幅 声波的波形畸变. 文献[24, 25] 研究了非线性声传 声悬浮是高声强条件下的一种非线性声学效 播过程中的声流和声辐射压现象. 在声悬浮研究方 应, 同时也是一项先进的材料无容器处理技术[12] . 面, Barrios 和Rechtman[26] 采用LB 方法对单轴式 近十年来, 声悬浮技术得到了越来越广泛的应用, 声悬浮过程进行了初步探索. 但是, 这些研究都是 其应用领域包括无容器材料制备加工[3−5] 、物性 基于二维格子模型, 不能完全反映三维物理过程的 [67] [89] [1011] 测定 、液滴动力学 、分析化学 和生物医 真实情况. 学[1213] 等. 作为一种非线性现象, 声悬浮过程中 近年来, 基于直角坐标系下二维格子的轴对 的声辐射压、声流等基本物理效应也得到了越来越 称三维LB 模型获得了较快的发展[27−29] . 由于利 多的关注[1415] . 在线性声学近似下, 声压和流速作 用了所研究问题的轴对称性, 计算效率得到充分 周期性变化, 两者的时间平均值均为零, 因此不会 提高. 另外, 传统的LB 方法通常采用单弛豫时间 出现稳态的声辐射压和声流. 为了揭示这些非线性 (SRT)碰撞模型[19] , 其缺点是在低黏度情况下会 效应, 必须利用原始的非线性流体动力学方程和状 出现数值不稳定性. 为了提高数值稳定性, 人们发 [16] [30] 态方程 . 但是, 流体动力学方程组的求解十分困 展了多弛豫时间(MRT)碰撞模型 , 并将其推广 难, 特别是在涉及声场与悬浮样品相互作用的情况 到轴对称LB 模拟中[2829] . 下. 数值模拟方法提供了一种可行途径. 采用最新发展的一种轴对称MRT-LB 方法, 格子Boltzmann (LB) 方法是近二十年发展起 本文对圆柱形封闭谐振