zhuzilin's Blog

about

声学基础

date: 2023-09-24
tags: 物理  

线性声学波动方程

声学主要考虑的是机械波在物质间传播相关的性质,而在物质上,一般会先关注气体。在这个前提下,在考虑声波的波动方程时,主要有以下的假设:

  • 质量守恒定律:

    (ρu)+ρt=0\nabla\cdot(\rho\textbf{u}) + \frac{\partial \rho}{\partial t}=0
  • 动量守恒定律:

    p+ρdudt=0\nabla p+\rho\frac{\text{d}\textbf{u}}{\text{d}t}=0

    其中:

    dudt=ut+uξxξxt+uξyξyt+uξzξztdξdt=d(ξx,ξy,ξz)dt=u(ξx,ξy,ξz,t)\begin{aligned} \frac{\text{d}\textbf{u}}{\text{d}t}&=\frac{\partial{\textbf{u}}}{\partial t}+ \frac{\partial{\textbf{u}}}{\partial \xi_x}\frac{\partial{\xi_x}}{\partial t}+ \frac{\partial{\textbf{u}}}{\partial \xi_y}\frac{\partial{\xi_y}}{\partial t}+ \frac{\partial{\textbf{u}}}{\partial \xi_z}\frac{\partial{\xi_z}}{\partial t}\\ \frac{\text{d}\xi}{\text{d}t}&=\frac{\text{d}(\xi_x,\xi_y,\xi_z)}{\text{d}t}=\textbf{u}(\xi_x,\xi_y,\xi_z,t) \end{aligned}
  • 能量守恒定律,这是声学推导中较为特殊的一部分,它假设声音的传播是绝热过程,即没有热量交换,所有外界做功均转变为气体的内能转化。由绝热过程,有:

    p=Ksργp=K_s\rho^{\gamma}

    其中,γ=cp/cv\gamma=c_p/c_v

根据 (4)(4) 有:

pt=pρρt=γpρρt\frac{\partial p}{\partial t}=\frac{\partial p}{\partial \rho}\frac{\partial \rho}{\partial t}=\gamma\frac{p}{\rho}\frac{\partial \rho}{\partial t}

令,p=p0+p,ρ=ρ0+ρp=p_0 + p', \rho=\rho_0 + \rho',其中 p0,ρ0p_0,\rho_0 为 static pressure 和 equilibrium density。并设:

c2=γp0ρ0c^2=\gamma\frac{p_0}{\rho_0}

那么假如我们忽略掉所有的二阶项,例如 ρu\rho'\textbf{u},则 (1)(1) 变为:

(ρ0+ρ)u+ρt=0ρ0u+1c2pt=0\begin{aligned} \nabla\cdot (\rho_0 + \rho')\textbf{u}+\frac{\partial \rho'}{\partial t}&=0\\ \rho_0\nabla\cdot\textbf{u}+\frac{1}{c^2}\frac{\partial p'}{\partial t}&=0 \end{aligned}

(2)(2) 变为:

p+ρ0dudt=0\nabla p'+\rho_0\frac{\text{d}\textbf{u}}{\text{d}t}=0

(8)(8)\nabla,有:

2p1c22pt2=0\nabla^2p-\frac{1}{c^2}\frac{\partial^2 p}{\partial t^2}=0

即得到了线性声学方程(linear acoustic equation)。并由理想气体方程 p=ρRTp=\rho RT ,得到声速为:

c=γp0/ρ0=γRTc=\sqrt{\gamma p_0/\rho_0}=\sqrt{\gamma RT}

即声速仅由温度决定。

上述的 2 个假设(绝热过程和忽略二阶项)合称 acoustic approximation。

声压和 db

声压的测量方式为:

Lp=10log(ps2)avgpref2L_p=10\log{\frac{(p_s^2)_{\text{avg}}}{p_{\text{ref}}^2}}

这里空气中 prefp_{\text{ref}} 一般取 2e-5 Pa。其单位为 db。