压力波是什么波【仿真百科】声学简介（下）

有关声音波动性的最早推测可以追溯到亚里士多德时代，当时的古希腊哲学家已经开始了这方面的研究。他们的灵感来源于小波在水面上的传播，这些波如何与障碍物相互作用，以及振动物体如何发出声音。马兰·梅森（Marin Mersenne，1588-1648）首次提出了振弦的正确（已出版）数学描述（1636 年出版的《Harmonicorum Libri》），由于这一突出贡献，人们通常称他为声学之父。这比伽利略·伽利莱（Galileo Galilei，1564-1642）对力学的开创性研究早了两年。在梅森和伽利略时代，其他物理学家并不相信声音的波动性，有些人认为它是一种物质流（见参考资料 1、21 和 22）。

艾萨克·牛顿（Isaac Newton，1643-1727）在他的著作《自然哲学的数学原理》（Principia）（1687 年出版的第二卷）中对声音在流体中的传播提出了第一个正式的数学描述。他将声音描述为脉冲传播，其中相邻的流体质点相互作用并相互推动。牛顿基于等温假设预测了空气中的声速约为 298 m/s（979 ft/s），这一声速值和理论后来于 1816 年被皮埃尔-西蒙·拉普拉斯（Pierre-Simon Laplace，1749-1827）修正为我们现在所知的 348 m/s（1142 ft/s），这是他通过  因子（其中  表示绝热指数；请参见“流体中的线性波”一节）修正牛顿的预测得到的值。

十八世纪，随着连续介质力学和流体动力学的发展，莱昂哈德·欧拉（Leonhard Euler，1707-1783）、约瑟夫·路易斯·拉格朗日（Joseph-Louis Lagrange，1736-1813）和让·勒朗·达朗贝尔（Jean le Rond d'Alembert，1717-1783）的研究取得了进展，他们为声音的传播提出了完整的理论。大部分现代声音传播理论都建立在这些科学家及其同时代人的工作成果之上。十九世纪，拉普拉斯、克劳德-路易·纳维（Claude-Louis Navier，1785-1836）、西莫恩·德尼·泊松（Siméon Denis Poisson，1781-1840）、奥古斯特·孔特（August Kundt，1839-1894）、赫尔曼·冯·亥姆霍兹（Hermann von Helmholtz，1821-1894）和乔治·斯托克斯（George Stokes，1819-1903）的研究成果推动这一理论继续向前发展。瑞利勋爵（Lord Rayleigh，1842-1919）于 1877 年发表的专著《声学理论》（The Theory of Sound）被认为是这一漫长研究工作（对散射的详细描述等）的完结篇，这一著作常常被认为标志着经典声学到现代声学的转变（参考资料 21-22）。

为了推导所有波动现象的控制方程，您需要从最一般形式的守恒方程开始，即质量守恒、动量守恒和能量守恒。为了封闭方程组，需要为这些方程补充本构关系和热力学状态方程。

根据波的传播介质以及波与波之间的相互作用，存在许多不同类型的波。在下面的章节中，我们将推导流体（液体和气体）中波的控制方程，包括损耗模型和各种假设的详细信息。在此基础上，我们将给出固体中弹性波的控制方程，以及弹性波和压力波在多孔材料中进行组合传播的控制方程。

其中，自变量为时间 （SI 单位：s）和空间坐标 （SI 单位：m）。因变量为密度 （SI 单位：kg/m3）；速度场 （SI 单位：m/s）和熵 （SI 单位：J/kg/K），即单位质量的熵。此外， 为温度（SI 单位：K）， 为黏性耗散函数（W/m3）， 为总应力张量（SI 单位：N/m2）， 为黏性应力张量（SI 单位：N/m2）， 为局部热通量（SI 单位：W/m2）， 表示可能的质量源项（SI 单位：kg/m3）， 表示可能的热源（SI 单位：W/m3）， 是可能的体积力源项（SI 单位：N/m3）。这里的守恒量是密度 、动量  和总熵 。

原则上，我们将考虑质量守恒情况，因此，一般情况下 。不过，我们可以使用质量源项的声学扰动（见下文）来表示我们不想详细描述的复杂过程。

编写守恒方程和选择因变量的方法有很多，以上只是其中一种。您可以参见参考资料 1-10 或任何关于流体力学和连续物理学的教科书，获取更多信息。

热力学关系

请勿将恒压热容  和等容热容  分别与等熵声速变量  和等温声速变量  进行混淆。我们假设  且 ，建立密度和熵的微分式

（7）

本构关系

本体黏度项模拟压缩和膨胀黏度效应，这两个效应实际上描述机械压力与热力学压力之间的差异，它们并不总是处于平衡状态。

此外，一般来说，所有材料属性可能都与温度和压力相关。这意味着，我们应该将材料属性看作与空间相关的量。

流体控制方程

这些方程中保留了一个非零质量源项 ，从而使动量方程和能量方程右边分别存在  和  非标准项。在几乎所有的应用中，由于质量守恒得到 ，因此都不包含这些项。

微扰理论、线性化和阶次

声学研究的是微扰的传输和传播。这些微扰可以比背景条件小多个数量级；例如，大气压约为 100,000 Pa，而正常语音信号的振幅则为 。通常情况下，求解上述完整的控制方程是不切实际的。这些方程本质上是非线性的，不仅如此，通过对比以下各项可知，它们在空间和时间上也是多尺度的：

• 空间坐标  和波长 

• 时间  和波动周期 

• 速度场  和声速 

在数值应用中，同时解析不同的物理尺度和时间尺度需要非常高的数值精度。

在大部分实际情况下，我们可以将声学问题化为线性问题，利用微扰理论，在脱离背景属性的情况下对声学进行简化分析。根据微扰展开的阶次、所模拟物理机制方程中保留项的阶数以及其他近似方法，对控制方程进行重新表述将得到不同的声学方程。这些方程可以描述从一般线性声学和亥姆霍兹方程到有关冲击的高级非线性声学模型和方程的一切过程。有关微扰理论的详细信息，请参见参考资料 1-5 和 11。

在微扰理论中，因变量 （压力、温度、速度或密度）可以展开为

其中  是一个小参数。

在大多数声学应用中，因变量都将展开为一阶，从而得到

其中去掉了 。

通常，背景场（零阶）只依赖于空间，或者与声学时间尺度相比，最多只能在时间上缓慢变化。在声流分析等特定应用中，将微扰展开到二阶以分离时间尺度（例如，见参考资料 26）。在许多教科书中，声压扰动用  表示。此处的  用来明确说明这是一阶微扰。

在线性声学中，微扰展开到一阶，控制方程中只保留一阶项（如下所示）。使用泰勒级数将状态方程中与零阶解相关的因变量展开为一阶。例如，在等熵情况（）下，可以得到

在大多数非线性声学（忽略边界层效应时）中，我们将因变量扰动展开到一阶，但控制方程中最多保留到二阶项（请参见下面的非线性声学一节）。在这种情况下，我们使用泰勒级数将状态方程展开到二阶（例如，见参考资料 25）。微扰方法不仅仅依赖于数学方案，而且在很大程度上与需要由控制方程捕获的物理效应有关。

流体中的线性波

在频域中，时间相关性为谐波形式，因变量和源可以用它们的傅里叶分量来表示。例如，压力可以表示为

通过使用一阶泰勒展开（或者通常使用上面给出的密度微分），将密度的扰动分解为压力和温度项。

我们需要理解的是，计算是在线性化点  上进行的。

这两个热力学量就是前面介绍的等温压缩系数  和（等压）热膨胀系数 ，它们与声速（行波速度）有重要的关系

（9）

其中，我们已经明确指出，声速  为等熵声速，是通常给出的表称值，称为声速。等温声速由  给出。

完整的线性方程

将上述方程插入控制方程中，仅保留扰动量中的线性项，即可得到完整的线性纳维-斯托克斯方程。将这些项重排得到

其中  是黏性耗散函数。（11）

回想一下我们前面提到的，假设背景场为稳态，或者与扰动相比，至少在时间上变化缓慢。不仅如此，一般情况下，我们假设 （由于流量总质量守恒）且 （对系统使用零熵参考）。不过，我们保留了扰动质量源项 。我们很少包含这一项，但可以用它来对不希望详细描述的复杂过程进行建模。比如，在声学中，脉动球或热注入作用可以用这种质量源项进行很好的近似。当使用二阶微扰方法时，也会出现类质量源项。

在频域中，我们将所有的时间导数  都替换为与  的乘积。在上面给出的方程中，没有对黏度或导热系数等材料属性进行扰动。在无损情况下，上述方程可以降阶为线性欧拉方程。线性欧拉方程和线性纳维-斯托克斯方程构成了气动声学领域的基础。当存在背景流  时，这些方程仍可以描述熵波和涡波等非声波的传播，它们是在背景流速的对流，对这些方程进行数值求解是一项巨大的挑战。

线性热黏性方程

这些方程明确地描述包含热损耗和黏滞损耗的声波传播，是热黏性声学方程，相关问题也称为热黏性声学。在模拟和研究微型声学器件（如移动设备和麦克风）时，为了正确预测器件的特性和响应，必须包含由热黏性方程描述的物理量。

一般线性标量波动方程和亥姆霍兹方程

其中，由于假设背景场的所有时间导数均为零且 ，因此我们使用零阶动量方程来表示 。

一般情况下，当密度  不是常数时，必须使用此方程。对于材料数据为常数并且无热源的情况，该方程可以退化为普通形式 。

这个方程组构成了推导[压力声学](Summary of Acoustics Pillar#pressure_acoustics)方程的基础，其中的数学问题得到简化，只需求解包含一个因变量（压力）的方程。请注意，类似的公式也可以用速度势或密度来表示。现在，我们可以重新表述标量波动方程。利用对连续性求时间导数以及对动量方程求散度的常用技巧，并结合压力的时间导数，可以得到

（13）

这是标量波动方程的最一般形式，其源项适用于密度和声速均非恒定的情况。在频域中，通过将时间导数替换成与  的乘积，可以重新表述亥姆霍兹方程，从而得到

（14）

其中，波数为 。

在材料数据与深度相关的水声学等领域，这种控制方程的一般形式至关重要。在多物理场环境中，消声器系统中的温度分布  就是这样一个例子，它会影响材料属性， 和 。这种情况下，使用正确的公式同样非常重要。

阻抗

阻抗的概念在机械工程和声学中尤其重要，它是一个频域概念，定义为边界或点上的力与流动变量之比。阻抗边界条件可用于在不显式建模的情况下施加边界属性。因此，该边界条件概括了声硬和声软边界条件，可用于分析许多介于这两个极端之间的情况。与频率相关的阻抗可用于模拟或简化复杂的机械系统或吸收边界的属性。人体耳膜（鼓膜）的阻抗便是其中一个例子，我们可以对这种阻抗进行详细的描述和研究。

在声学中，我们通常使用三种不同的阻抗概念：声阻抗、比阻抗和机械阻抗。声阻抗 （SI 单位：Pa*s/m3）定义为边界上的平均压力  与通过边界的体积流率 （SI 单温：m3/s）之比，

这个量有时也称为面法向阻抗 ，通常用于指定或定义边界条件。

比声阻抗 （SI 单位：Pa*s/m）定义为某个点上的压力  与粒子速度  之比

由每个点的压力与粒子速度的关系，不同的波型可以通过其各自的特性比声阻抗来定义。无损平面行波的特性（比声）阻抗由已知量  给出。最后，机械阻抗  定义为边界上的作用力 （SI 单位：N）与粒子速度之比

强度

声波的传播与能量的传递或能量的流动有关。声强大小 （或简称强度大小）定义为单位时间内通过单位面积的时间平均能量，其中面积的法向为波的传播方向。一般来说，强度是描述能量传递大小和方向的矢量 ，通过瞬时强度的时间平均值计算

其中，瞬时强度  的精确形式取决于求解的波动方程。

 的表达式遵循声能推论（论述声波能量守恒；见参考资料 2 和 9），并具有不同的复杂性。积分时间  取决于声信号的类型。对于噪声来说，这是一段很长的时间，而对于谐波信号，则是信号周期。对于在无损静止流体中平面传播的波，其强度由下式给出

在更广泛的情况下，对于在频域求解的静态流体中的波动问题（求解亥姆霍兹方程），瞬时强度（强度为 ）由下式给出

其中  表示共轭复数。

在参考资料 9 中，Myers 给出了适用于所有线性声学的声能推论的非常普遍的形式。

在分析能量传递和耗散时，了解系统的强度（矢量场分布）非常有用。上图显示通过吸收式消声器系统的强度。对上述消声器的入口和出口等边界上的强度进行积分，可以给出系统的总功率输入/输出。具体来说，知道边界上入射波和反射波的组合，我们可以计算传输损耗，即：发射功率与入射功率之比。对于扬声器系统，对封闭表面上法向强度的积分可以定义辐射功率。

举例来说，计算和预测人体组织或多孔材料中吸收的声能也会涉及强度。简化的半解析表达式往往需要知道强度，比如平面波吸收定律 ，其中  为耗散功率密度（SI 单位：W/m3）， 表示波数的虚部。

边界条件

在求解控制方程时，我们需要使用适当的边界条件在数学上形成封闭。这些条件描述声场如何产生，及其在边界上的特性，以及在无限域中辐射时的特性。

在没有平均背景流的情况下，声场与其周围环境相互作用的边界处会显示声学的多物理性质。在弹性结构的边界处，声场将受到结构加速度的影响，与此同时，结构也将受到压力载荷的影响。这就是典型的声振耦合或声-结构相互作用。

从数学上讲，界面上施加了两个条件：一个是表示位移（或速度）连续性的运动学条件，另一个是表示应力连续性的动力学条件。在扬声器等声源的表面，不仅会产生声音，声音还会通过附加质量影响扬声器，最后再耦合回电源。

声-结构相互作用示例，其中的声音由箱内的轮系总成发出。

正如前面讨论的，阻抗模型可以用作边界条件，以简化的方式对复杂系统的特性进行建模；例如，我们可以指定多孔吸收体表面的阻抗，而无需进行详细建模（参考资料 14 和 15）。在计算声学中，阻抗条件被广泛用于降低模型的复杂性。由于阻抗模型可以近似表征复杂系统的真实特性并提供简化描述，因此，我们可以将这些模型视为工程关系。在求解复杂问题时，使用阻抗条件常常可以给出比较准确的初步近似模型。

声波进入或离开系统时的特性可以用辐射条件或端口条件来描述。在数值仿真中，我们不能总是使用理想化条件，因此这类条件显得尤为重要。在波导等系统的入口处，我们应该允许波进入波导，并允许反射波离开波导。这是辐射或非反射边界条件（NRBC）的一个例子。在管道的入口和出口处，此类条件可以作为端口条件来实现。在开放系统中，波可以向无穷远处自由传播，由于波与源和散射体的距离趋于无穷，因此必须满足索末菲辐射条件。

在数值仿真中，为了避免可能影响求解精度的虚假数值反射，正确实现和表述这一条件非常重要。Bayliss、Gunzburger 和 Turkel（参考资料 23）给出了此类条件的经典示例。

一种可以替代边界条件的方法是使用泡沫材料层来模拟无限开放问题。在这样的域中，出射波在几乎没有反射的情况下被抑制。我们熟知的完美匹配层（PML）就是这样一种吸收层方法。对于实际坐标拉伸与人工数值扩散相结合的时域问题，还存在其他的表述方法。通常，我们把这些表述称为吸收层。

等效流体模型

当声波在多孔材料或狭窄管道等物体中传播时，所涉及的损耗机制可能会非常复杂，我们很难根据基本原理进行描述和建模。对于可以详细描述损耗机制的其他情况，也很难找到问题的解析解或数值解。在此类情况下，使用等效流体模型以均匀化方式描述损耗过程比较有利。这种模型的典型应用是在求解亥姆霍兹方程时，通过定义复值密度  和声速 （这两个参数通常与频率有关）来引入损耗。

我们可以使用许多等效流体模型来描述多孔材料中的损耗；例如 Delany-Bazley-Miki 模型（参考资料 14、15 和 24）和 Johnson-Champoux-Allard 模型（参考资料 14）。其原理是描述等效流-固组合系统（饱和流体和多孔基体）的频率相关有效流体密度  和有效流体的体积模量 。此类模型的有效性局限于特定的参数范围，这些范围表示刚性或柔性多孔基体近似。在数值仿真中，等效流体模型的计算要求低于求解 Biot 方程的完整多孔弹性波模型。

当声波在狭窄区域中传播时，我们需要考虑与声边界层相关的损耗（黏滞损耗和热损耗）。在大多数情况下，这需要求解完整的线性热黏性声学方程。对于窄缝或等截面波导，我们可以使用均质模型，这些模型通常称为低约化频率模型（请参见“热黏性声学”一节）。

此外，均质等效流体模型还可用于求解水声学，其中使用衰减系数  来定义损耗。计算值通常与频率有关，以奈培/米、分贝/公里或分贝/波长（参考资料 12）表示。

固体中的弹性波

当声压波与固体相互作用时，流体的压力波动在固体域上产生流体载荷，并且结构加速度作为通过流-固边界的法向加速度影响流体域。声音在固体中的传播通过固体形状和结构的小幅弹性振动来实现。

固体中线弹性波传播的控制方程是从一般的结构力学控制方程推导出来的。我们在应力和应变小扰动的极限下对方程进行线性化处理和公式化表示。在固体材料中，应力与应变张量之间最普适的线性关系式可以写为

其中， 为柯西应力张量， 为应变张量， 为（四阶）弹性张量。

这是适用于线弹性材料的胡克定律。对于小扰动，应变张量（矢量形式）由下式给出

其中  是位移场（SI 单位：m）。

弹性波动方程由牛顿第二定律（动量守恒）得出，给定为

其中， 表示材料密度（SI 单位：kg/m3）， 和  是可能的源。

其中  表示角频率。

多孔弹性波

多孔弹性波描述压力波和弹性波在多孔材料中的组合传播：压力波在孔隙中的饱和流体中传播，弹性波在多孔基体框架中传播。在多孔基体框架几乎静止（刚性）的极限下，我们可以使用均质等效流体模型来描述声学特性。在柔性极限下也是如此，其中框架跟随流体产生运动。不过，这种描述并不适用于所有频率或材料参数。此外，当多孔材料与振动的固体表面发生接触时，还需要包含弹性波。Biot 理论给出了一个模型，描述了求解位移  和压力  的组合传播。

在 Biot 于 1956 年发表的开创性著作中，他将线弹性经典理论扩展到了充满流体的饱和多孔介质中（见参考资料 16-18）。在他的理论中，体积模量和压缩系数与波频率无关，因此可以作为常数参数处理。多孔基体通过线弹性进行描述。我们可以通过考虑孔隙中饱和流体的黏度和热传导引起的频率相关损耗来引入阻尼。包含热效应的 Biot 方程的公式有时称为 Biot-Allard 理论（见参考资料 14、19 和 20）。控制方程的一种形式为

其中  是流体密度； 是流-固系统的平均密度； 是 Biot-Willis 系数（SI 单位：1），描述流体框架从刚性  到柔性  的耦合特性。在此公式中，Biot-Willis 系数假设为空间常数。

复值密度和波数定义如下

其中， 是迂曲度（SI 单位：1）或结构形式因子（高频极限）， 是孔隙率（SI 单位：1）， 是与流阻率相关的渗透率（SI 单位：m2）， 是基体的排水体积模量（真空时的体积模量）。

黏度函数  和流体压缩系数  这两个频率相关的量分别包含有关黏滞损耗和热损耗的信息。因此，我们可以将控制方程识别为具有对角应力源项（通过压力定义）的弹性波动方程（15） 和力源项（通过位移定义）的亥姆霍兹方程（14）。

为了使用声学控制方程（如上文“控制方程的详细推导”一节中所述）来预测换能器的特性或房间的响应等，我们需要在计算框架中对这些方程进行数值求解。这些方程只存在有限数量的解析解，这些解不适用于实际的几何、复杂边界条件或多物理场应用。求解声学控制方程的技术多种多样，但每种方法都有各自的优缺点。通常来说，不同的应用需要使用不同的方法。

在求解声学问题时，我们要做的第一个选择是，在频域还是时域求解波动方程。当然，这是一个宽泛的问题，原因是有些方法同时依赖于这两种方法。对于所有的公式，我们都需要所用的方法能同时在时域/频域和空间中解析波，同时还需要考虑频率组成、波长、几何细节、可能的声学边界层等各种因素。

由此，我们可以使用以下多种方法中的一种来求解方程：

• 有限元法（FEM）

• 积分法，如边界元法（BEM）

• 间断伽辽金（DG）法

• 高频方法，如射线追踪和波束追踪

• 时域有限差分（FDTD）

• 基于波的方法

• 有限体积法

• 波数积分技术

• 固有模态

• 抛物线方程法

请注意，在气动声学领域，有一个专门的学科叫做“计算气动声学”（CAA），用于求解气流与声学之间的相互作用。

压力声学

在压力声学中，只使用一个因变量来表示声场：压力。在时域求解标量波动方程（13），并在频域求解亥姆霍兹方程（14）。从某种意义上说，这就是经典声学。

对于许多情况下的声学，我们都可以使用这两个方程进行描述，比如，散射和辐射问题、消声器系统、汽车车厢声学，等等。这两个方程的一般形式可用于模拟具有非恒定密度和声速的问题，例如分析水声学或具有背景温度分布的多物理场问题。

我们可以在各种边界条件下求解这两个方程。在频域中，使用等效流体模型（均质）能够以简化的方式在多孔介质和窄波导中添加复杂的损耗效应。阻抗边界条件可以用来模拟与复杂环境的相互作用。

声学和多物理场

声-结构相互作用

有关固体中的弹性波与流体中的压力波相互作用的多物理场分析称为声-结构相互作用或振动声学。这是声学中一种特殊形式的流-固耦合（FSI）。

在大多数实际应用中，声音由振动表面产生或通过结构进行传递，因此声-结构相互作用非常重要。这些应用包括消声器、扬声器、隔音问题和机器声学的设计和优化。

多孔弹性波

在多孔材料中，为了分析多孔基体中的弹性波与饱和孔隙流体中的压力波的耦合传播，需要求解 Biot 方程，其中包含孔隙流体对组合波传播的阻尼作用。

多孔弹性应用的例子包括弹性波在岩石、土壤和沉积物中的传播，以及微粒过滤器、消声器和车厢多孔衬里的声衰减属性建模。

压电声学

压电声学是一种多物理现象，用于解释如何通过压电材料产生和检测声信号。压电换能器是压电声学最重要的应用领域，例如，水下声呐应用中需要使用此类换能器。

管道中的声波

在管道系统中，声波的传播具有主导方向，此时可以简化数学描述。我们可以对管道横截面的影响进行积分，将三维问题简化为一维方程组，从而可以对大型管道系统进行显著简化的数值处理。不仅如此，还可以通过修改流体的有效压缩系数来包含管道弹性壁的影响。

电声学

电声学主要研究如何通过机电声换能器产生和检测声信号，包括扬声器、耳机、麦克风和 MEMS 换能器；近年来，移动设备中音频组件的建模和设计（在很大程度上）都属于这一范畴。

对这些设备的建模体现了实际的多物理场应用，其中，物理场之间的双向耦合非常重要。仿真可以包含电磁学、结构力学和声学的完整详细模型，然后通过有限元法求解模型，也可以将 FEM 与使用电路/SPICE 图的集总表示相结合进行求解。

热黏性声学

在模拟或研究移动设备、麦克风、助听器和穿孔板等微型声学器件时，为了正确预测这些器件的特性和响应，必须包含热黏性声学方程描述的物理效应，其中，与黏性声学和热声边界层相关的损耗尤为重要。

在此类应用中，决定其是否为小器件的因素不是特征几何尺寸，而是这些层的厚度。此外，热黏性声学还可以精确预测电容式麦克风等器件的小缝隙中从绝热到等温特性的过渡。

气动声学

气动声学的研究涉及背景平均流与声场之间的相互作用，包括湍流或周期性流动结构产生的声音，以及流体流动对声场的影响。

应用示例包括消声器设计，其中平均流会显著影响消声器性能；科里奥利流量计；以及喷气发动机噪音的产生和传播。

在解决气动声学问题时，我们可以求解许多控制方程公式，包括声学类比方程、线性公式，以及求解完全非线性流动问题的所有细节。部分示例如下：

• Lighthill 声学类比方程

• 线性纳维-斯托克斯方程

• 线性欧拉方程

• 线性势流方程

• 声学扰动方程（APE）

• 完全可压缩流动方程的直接数值仿真（DNS）

几何声学

几何声学提供求解高频极限下的声学问题的方法，此时所研究的特征波长远小于几何特征。其中包括假设声音沿着射线或管道传播的射线追踪方法，这些方法不仅用于对房间、音乐厅和车厢声学进行建模，还广泛用于水声学、室外噪声传播和声音映射等多个方面。

能量法

另一类用于求解大规模声学问题的方法称为能量法，包括统计能量分析（SEA）方法、能量有限元法（EFEM）以及其他形式的声学扩散方程。这些方法不解析声波，而仅限于声波和振动的能量传递，适用于汽车和房屋等复杂的耦合系统以及纯室内声学。

水声学

水声学主要研究声音在海洋和其他水体中的传播，包括研究与海底（海床）和海面以及与潜艇等水下物体的相互作用。由于声音在海洋中进行长距离传播，因此，衰减、声速和密度随深度的变化等影响变得非常重要。典型的应用领域包括声呐、地震勘探、水下通信和海洋生物学。

非线性声学

当声波的振幅变得足够大时，最常用的线性假设不再有效。例如，行波的波谷和波峰将不再以相同的速度运动，信号也会失真。涡度、非常大的时间尺度和尺寸差异等其他影响因素也会引起非线性现象。为了捕捉这些影响，我们需要使用非线性声学公式，或者可能需要求解完全可压缩流动方程。

在超声造影和扫描等需要高振幅的超声应用中，我们常常会遇到非线性效应。实际上，某些成像技术依靠非线性和产生二次谐波来提高分辨率。我们可以在扬声器端口和穿孔板中观察到因旋涡产生引起的非线性效应，这些位置会产生声音失真。在微流体（声流体）中，由于时间尺度差异非常大，导致非线性贡献累积，从而产生声流和辐射力等声泳效应。

1. D.T. Blackstock, Fundamentals of Physical Acoustics, John Wiley & Sons, 2000.

2. A.D. Pierce, Acoustics, 'An Introduction to Its Physical Principles and Applications,' Acoustical Society of America, 1991.

3. S. Temkin, 'Elements of Acoustics,' Acoustical Society of America, 2001.

4. P.M. Morse and K.U. Ingard, Theoretical Acoustics, McGraw-Hill (Princeton University Press), 1986.

5. H. Bruus, Theoretical Microfluidics, Oxford University Press, 2010.

6. G.K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, 2000.

7. L.D. Landau and E. M. Lifshitz, Fluid Mechanics, Course of Theoretical Physics, Volume 6, Butterworth-Heinemann, 2003.

8. B. Lautrup, Physics of Continuous Matter, Exotic and Every Day Phenomena in the Macroscopic World, 2nd ed., CRC Press, 2011.

9. M.K. Myers, 'Transport of energy by disturbances in arbitrary steady flow,' J. Fluid Mech. 226, pp. 383–400, 1991.

10. F.V. Hunt, 'Notes on the Exact Equations Governing the Propagation of Sound in Fluids,' J. Acoust. Am. 27, pp. 1019–1039, 1955.

11. C.M. Bender and S.A. Orszag, Advanced Mathematical Methods for Scientists and Engineers: Asymptotic Methods and Perturbation Theory, Springer, 1999.

12. F.B. Jensen, W.A. Kuperman, M.B. Porter, and H. Schmidt, Computational Ocean Acoustics, Springer, 2011.

13. H. Kuttruff, Room Acoustics, CRC Press, 2009.

14. J.F. Allard and N. Atalla, Propagation of Sound in Porous Media, Modeling Sound Absorbing Materials, 2nd Edition, John Wiley & Sons, 2009.

15. T.J. Cox and P. D’Antonio, Acoustic Absorbers and Diffusers, Taylor and Francis, 2nd ed., 2009.

16. M.A. Biot, “Theory of Propagation of Elastic Waves in a Fluid-saturated Porous Solid. I. Low-frequency Range,” J. Acoust. Soc. Am., vol. 28, pp. 168–178, 1956.

17. M.A. Biot, “Theory of Propagation of Elastic Waves in a Fluid-saturated Porous Solid. II. Higher Frequency Range,” J. Acoust. Soc. Am., vol. 28, pp. 179–191, 1956.

18. M.A. Biot, “Generalized Theory of Acoustic Propagation in Porous Dissipative Media,” J. Acoust. Soc. Am., vol. 34, pp. 1254–1264, 1962.

19. J.F. Atalla, R. Panneton, and P. Debergue, 'A mixed displacement-pressure formulation for poroelastic materials,' J. Acoust. Soc. Am., vol. 104, pp. 1444–1452, 1998.

20. J.F. Atalla, M.A. Hamdi, and R. Panneton, 'Enhanced weak integral formulation for the mixed (u-p) poroelastic equations,' J. Acoust. Soc. Am. vol. 109, pp. 3065–3068, 2001.

21. J.W.S. Rayleigh, The Theory of Sound: Introduction by R.B. Lindsay, Dover Publications, 1945.

22. R.B. Lindsay, 'The Story of Acoustics,' J. Acoust. Soc. Am., vol. 39, pp. 629–644, 1966.

23. A. Bayliss, M. Gunzburger, and E. Turkel, “Boundary Conditions for the Numerical Solution of Elliptic Equations in Exterior Regions,” SIAM J. Appl. Math., vol. 42, no. 2, pp. 430–451, 1982.

24. Y. Miki, “Acoustical properties of porous materials – modifications of Delany-Bazley models,” J. Acoust. Soc. Jpn (E), vol. 11, no. 1, 1990.

25. M.F. Hamilton and D.T. Blackstock (edited by), Nonlinear Acoustics, Acoustical Society of America, 2008.

26. P. Barkholt Muller and H. Bruus, 'Numerical study of thermoviscous effects in ultrasound-induced acoustic streaming in microchannels,' Phys. Rev. E 90, 043016 1–12, 2014.