平面驻波下小球声悬浮机理的理论分析和格子Boltzmann.doc

第八届全国流体力学学术会议 2014年9月18～21日 甘肃兰州 文章编号： CSTAM2014-B01-0263 标题：平面驻波下小球声悬浮机理的理论分析和格子Boltzmann数值模拟 作者：任虎，陈效鹏 单位：西北工业大学 Copyright ? 2014 版权所有 中国力学学会地址: 北京市北四环西路15号　邮政编码:100190 　Address: No.15 Beisihuanxi Road, Beijing 10 任虎*，陈效鹏*，2） *（西北工业大学，陕西西安 710072） 摘要 利用格子Boltzmann（LBM）在弱可压缩流体动力学模拟方面的特性，本文开展了小球（3D）在平面驻波作用下所致流场的数值模拟，并研究了不同粘性、不同障碍物尺寸对流场及举力的影响。通过与经典基于势流理论的Gorkov模型相对比，LBM方法显示了较高的模拟精度；同时验证了LBM可以被用于粘性介质声传播机制、声致辐射压力的预测。本研究对声悬浮机理的定量分析，对于声悬浮设备的研制，及声悬浮应用技术的发展有促进作用。 关键词 LBM，边界条件，平面驻波，声悬浮 引 言 众所周知，微小物体在声场中会受到时间无关的净力，这是声悬浮力源，一般被称作“声辐射力”。在无粘的条件下，声辐射力主要是由声场中一个时间无关的压力梯度在几何边界上积分得来的。 这个时间无关的压力梯度称作“辐射声压”。最早计算辐射力的是King5。 他用一阶线性波动方程作为基本流动，在二阶方程中使用一阶方程的结果作为驱动，得到了无粘可压的辐射力结果。而Gorcov1得到的结果比King5还要简单。他在考虑一阶波动时忽略了声场中小物体的影响，而只在结果上做一些修正，理由是其尺寸比波长要小得多。在有粘的情况下，问题会变得复杂。由于粘性和边界作用，流体内部还会产生净流动，叫做“声致冲流”，而声致冲流会影响流体和小物体之间的能量和动量交换10。关于声致冲流的分析，其实比声辐射力更早得到关注，最早可追溯到Rayleigh7。他用一个不可压缩模型渐进分析的方法得到了昆特管中的声致冲流结果，然而结论只关注了离边界较远的流场作用。其后Schlichting首先发展了边界层内部冲流的理论分析方法，而Stuart甚至使用了双层的边界层模型2。Westervelt虽然使用了跟Rayleigh相似的方法，但是进一步突破了不可压缩的假设2。 关于声致冲流的总结有很多，比如Said和Philip2，及Riley2，还有最近的Sadhal8。然而，为了简化分析的难度，理论分析部分在本文中只限于无粘假设和没有冲流的情况。 数值模拟中本文呈现的是格子Boltzmann方法的结果。这种方法非常适用于微可压缩有粘流动，因此对于粘性介质声学问题可以得到比较好的结果，从后文的结果对比也可看出。关于格子Boltzmann方法的详细介绍，读者可以参考wolf等的文献4, 6, 9, 11。在数值模拟中，我们单独研究了粘性和微小物体尺寸对于辐射力的影响并与无粘的理论结果进行了对比。 本研究对声悬浮机理的定量分析，对于声悬浮设备的研制，及声悬浮应用技术的发展有促进作用。 1 渐进分析方法 1.1 基本方程和无量纲化 首先，针对本问题的特性，粘性只在较薄的边界层起显著作用，我们用无粘模型来分析整个流场是合理和方便的。无粘的质量和动量方程可以表示为： (1) (2) 接下来，我们根据问题的特点引入无量纲化单位系统。长度单位，假设具有波长的量级；时间单位，假设具有基本波动的周期量级，即；而速度单位，与实际自变量处于同一量级，声速单位具有实际声速量级。考虑状态方程, 所以方程(1)和(2)各项都含有一次的密度因子，因此无量纲化不需要对密度进行处理，而直接无量纲化。方程(1)和(2)可以写成： (3) (4) 其中，是一个小量。而第二方程右端的，即具有“1”的量级。 1.2 摄动法展开 既然是一个小量，那么，以其为摄动参数，可以将自变量按量级展开为： (5a) (5b) (5c) 将方程(5)代入方程(3)和(4)得到各阶方程，其中零阶方程为 (3-I) (4-I) 一阶方程为： (3-II) (4-II) 二阶方程为 (3-III) (4-III) 以上摄动展开法主要参考了Eckart3的过程。 2 辐射力的计算 2.1 时间无关的压力