基于格子boltzmann方法的适用于声学运动边界的比较研究.pdfVIP

第十五届全国计算流体力学会议论文集 比较研究 陈力1，余洋1，侯国祥1 (1，华中科技大学船舶与海洋工程学院，湖北省武汉市洪山区珞瑜路1037号，邮编430074) 摘要：本文在基于格子Boltzmann方法的基础上，对现有的流体与固体的运动边界做了详细的比较分析，并 对其声学适应性做出了评价。文中比较了格子Boltzmann方法中常见的及新兴的曲面流固运动边界，包括 格线反弹格式、插值反弹格式、统一的插值反弹格式、显式和隐式直接力浸入格式、锐化浸入格式和虚拟 流浸入格式。针对以上各种流固运动边界在格点发生变化时产生的新流体点对流体压力产生的虚假噪声， 本文提出了几种改进方案。通过模拟声波的圆柱反射和拖动圆柱绕流的问题，发现统一的插值反弹格式和 隐式直接力浸入格式比较适合做运动声学的模拟，同时证明了本文所提出的改进方案的有效性。 关键词：格子Boltzmann方法；声学模拟；运动边界；反弹边界；浸入式边界 以上的反弹格式都具有二阶精度，由于其简易性及 0引言 高效性，在基于LBM的运动边界问题中被广泛应 用。 BoltzmannMethod． 对于复杂边界问题，浸入式边界方法 格子Boltzmann方法(Lattice LBM)是一种可以高效的处理复杂流动的问题的介(Immersed．boundary 势，LBM被广泛应用之后，IBM也被成功引入 观方法11]。由于LBM可以看作对不可压 LBM学J，并且成功的被用于解决复杂边界的问题， Navier-Stockes方程在微可压条件下的还原，因此的 密度的低可压的性质可以适用于模拟声学现象，经 在固体变形问题和流固耦合问题中很有优势。但是 过一些学者的测试证实LBM可以很好的捕捉非线 Peng等[91发现IBM无法满足无滑移条件，且在计算 性声学中的压力波动【2】。此外，还有天然并行性和 精度上不及插值反弹格式。随后Shu等【101用速度修 处理复杂边晃的简单性等优势。 正的方法使其强制满足了无滑移条件。不久，Luo 对于平直边界，半步长反弹格式(HalfWa)， 等fll】用多重直接力模型修正了这个缺陷，同时 Bounce．Back)是使用广泛、高效且处理简单的格 式，LaddD]将其推广到曲面问题上，提出了格线反 类似的，Tiwari[13】等从有限差分法中引入了虚拟流 弹格式(LinkBounce．Back)。但是该格式用折线近似 Fluid 浸入格式(GhostIBM)，自然的满足了无滑移 物理边界，降低了精度。为了提高精度，Filippova、 条件，同时在流固边界上可以满足Dirichlet或者 Mei和Bouzidi等M】先后提出了类似的插值反弹格 Neumann边界条件。根据插值格式的不同，浸入格 式(InterpolationBounce．Back)。该格式具有无滑移、 适应于曲面边界等特点，但是由于在处理过程中， 式可分为两大类，一类是散布式插值格式(Diffused interface 根据节点与曲面的距离不同，采用硒种不同的插值 scheme)，比如显式和隐式直接力模型，一