空气中强声波非线性传播特性数值模拟研究.pdfVIP

第28卷第5期 声学技术 V01．28．NO．5Pt．2 2009年10月 TechnicalAcoustics Oct．，2009 空气中强声波非线性传播特性的数值模拟研究 张振福1，赵云1，曾新吾2 (1．国防科学技术大学理学院，湖南，长沙。410073，2．国防科学技术大学光电学院。湖南，长沙。410073) onthecharacteristicsofnon—linear of Study propagationhigh ● ●J ·● ● ● l ■ l ● _ - USt|CWaVeln numericalsimulation intensityoca．rIa by 2 ZHANG Xin．wu Zhen．ful，ZHAOYunl，ZENG (1．College 410073，China, ofScience,NationalUniv．ofDefenseTechnology,Changsha Science ofOptoelectric 2．College andEngineering,NationalUniv．ofDefenseTechnology,Changsha410073，China) 1己I皇 1 √l口 虿02p=毒瓦翻+晌 强声波在空气中传播，会产生非线性效应[1】，并 且这种效应会随着传播距离的增加而积累。 矿02p=专万锄删母， 由于强声波在介质中的传播问题的复杂性，通常都 是无法得到解析解，需要用近似的手段来求解，这 鲁：去瓦掘删辞， 就使得数值模拟和实验成为研究非线性声波传播的 有效手段。根据声扰动的特点，本文采用时域有限 可c32p=专c瞄t衰磅M影， 差分法(FDTD)对空气介质中强声波的非线性传播 1 (2∥0—5∥j1+4p孑2一p孑3)+o(彭) 进行了研究，并进行了相应的外场实验。 砰 2Westervelt方程及差分格式 得到 2．1Westervelt方程 fk兀×乙丢《)一乃 ∥=等毛+等气+竺Poco 本文采用的是由流体力学方程直接导出的二阶 可以看出，该时刻的物理量依赖于前面时刻的物理 近似非线线性波动方程——westeⅣelt方程[2】：量，这显然是符合物理规律的。 3模拟结果与分析 V2产虿1可a2p一扣·％+去等=。 前两项分别代表衍射项和线性传播项，第三项由环 现实生活中常见的声源多为活塞式声源，所谓 境介质密度的非均匀性引起，对于均匀介质该项为 活塞式声源，就是一种平面的振子。例如扬声器纸 零，波传播的非线性效果由方程的最后～项控制。