格子Boltzmann方法应用于气动声学研究.pdfVIP

第 3O卷 第6期 计 算 物 理 V0l_3O．NO．6 2013年 11月 CHINESE JOURNALOFCOMPUTATIONALPHYSICS NOV．。2013 文章编号：1001—246X(2013)06—0808-07 格子 Boltzmann方法应用于气动声学研究 邓义求， 唐 政 ， 董宇红 (上海大学 ，上海市应用数学和力学研究所 ，上海 200072) 摘 要：应用格子 Boltzmann方法(LBM)对不 同类型的气动声学 问题进行数值研究．通过模拟一维平面声波和二 维点源声波的传播 ，得到沿传播方向的声压脉动，其振荡幅值和衰减趋势与理论值相吻合．其次进行声波衍射和干 涉现象的数值模拟．最后，模拟处于流场中运动声源辐射声场的多普勒效应．模拟结果 说明LBM方法能较好地模 拟低马赫数下的声学问题 ，包括声压脉动的传播 ，声波 的波动特性以及流动与声波 间的相互作用． 关键词 ：气动声学；格子波尔兹曼方法 ；声辐射 ；多普勒效应 中图分类号 ：0422 文献标识码 ：A O 引言 随着工业和科技的高速发展 ，民用和军用交通工具 的速度在不断的提升 ，而高速的行驶也带来 了日趋严 重的噪声问题．这些噪声不仅会困扰驾驶员、乘客和周 围生活的居民，还会引起交通工具结构的疲劳和破坏 ， 同时也会影响军用设备的隐蔽性 ，理解、预测以及控制气动噪声的研究正 日益受到人们的重视． Lighthill 在 1952年开创性地建立 了流动产生声音 的最一般理论 ，奠定了气动噪声 的理论基础 ，之 后人们对气动噪声计算方法的研究始终没有停止．随着计算机效能的显著提高，计算气动声学成为研究 和预测气动噪声问题的重要途径 ，其 目的是用计算的手段预测流动诱导 的声源或流体 中的运动固体所辐 射 的近场和远场噪声．采用 的数值方法必须能够正确计算所产生 的声源信息，同时准确预测声场辐射．根 据 Lighthill的声 比拟理论 ，流场和流动诱导的声场始 自于同一基本方程——Navier．Stokes方程 ，因此理论 上可 以在一个统一的区域 (近场和远场)中，同时计算声源在流场 中的产生和声场 的辐射．这是所谓 的直 接数值模拟方法，它可以直观而清晰地研究流体与声波间的相互作用．但是，由于流动的脉动谱域宽阔， 声能与涡能、声压与流场宏观压力在量级上的巨大差异，要求计算方法对波动频率、幅值、波阵面以及波 传播速度能够准确模拟，因而对网格尺度、时间步长以及数值方法的色散和耗散误差及稳定性都有着很 高的要求． 格子 Bohzmann方法 (1atticeBoltzmannmethod，LBM)是近二十年来 由格子气 自动机 (LGA)发展而来 的一种高效的模拟方法 ，该方法遵循质量、动量和能量守恒等基本物理规律，用简单的介观模型来模 拟流动和波动 ．一方面 ，尽管该方法在时空上是二阶精度 ，但它的理论基础跟 Navier．Stokes方程不 同，最 新的研究表明 ，LBM方法的耗散和色散误差小于从 Navier-Stokes方程 出发 的高阶格式 (如三阶 Runge． Kutta格式)．另一方面，相 比传统的计算流体力学 (CFD)方法，LBM方法有着算法实施简单，物理图象清 晰，边界处理方便 ，压力可直接求解等优势 ，该方法 的另一个显著优点是非常适合于并行计算．这些优点 使其有望较好地解决复杂流动和气动声学模拟所面临的困难．1998年 ，Buick等 首次将 LBM方法应用 于模拟声音的传播 ，以D2Q7模型对声波衰减进行 了数值研究，模拟获得 的波动频率及衰减系数与理论 值吻合 良好．王勇和何雅玲等 在 Buick等 工作的基础上，研究了一维平面声波和二维管道声波的衰 减 ，得到了与理论值吻合的速度分布和压力梯度分布 ，并给 出了声源 的激发声压级．Dong等 。。基于格 收稿 日期 ：2012—12—03；修 回日期 ：2013—04—06 基金项 目：国家 自然科学基金 11272198)资助项 目 作者简介 ：邓义求 (1985一)，男 ，硕士 ，主要从事流动噪声的数值模拟 通讯作者 ：董宇红(1968一)，E-mail：dongyh