Kelly James F, McGough Robert J
Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, USA.
J Acoust Soc Am. 2008 Apr;123(4):2107-16. doi: 10.1121/1.2885737.
The causal impulse response of the velocity potential for the Stokes wave equation is derived for calculations of transient velocity potential fields generated by circular pistons in viscous media. The causal Green's function is numerically verified using the material impulse response function approach. The causal, lossy impulse response for a baffled circular piston is then calculated within the near field and the far field regions using expressions previously derived for the fast near field method. Transient velocity potential fields in viscous media are computed with the causal, lossy impulse response and compared to results obtained with the lossless impulse response. The numerical error in the computed velocity potential field is quantitatively analyzed for a range of viscous relaxation times and piston radii. Results show that the largest errors are generated in locations near the piston face and for large relaxation times, and errors are relatively small otherwise. Unlike previous frequency-domain methods that require numerical inverse Fourier transforms for the evaluation of the lossy impulse response, the present approach calculates the lossy impulse response directly in the time domain. The results indicate that this causal impulse response is ideal for time-domain calculations that simultaneously account for diffraction and quadratic frequency-dependent attenuation in viscous media.
为了计算粘性介质中圆形活塞产生的瞬态速度势场,推导了斯托克斯波动方程速度势的因果脉冲响应。使用材料脉冲响应函数方法对因果格林函数进行了数值验证。然后,利用先前为快速近场方法推导的表达式,计算了在近场和远场区域内有障板圆形活塞的因果、有损脉冲响应。用因果、有损脉冲响应计算粘性介质中的瞬态速度势场,并与无损脉冲响应得到的结果进行比较。针对一系列粘性弛豫时间和活塞半径,对计算得到的速度势场中的数值误差进行了定量分析。结果表明,最大误差出现在活塞表面附近且弛豫时间较大时,其他情况下误差相对较小。与以往需要对有损脉冲响应进行数值傅里叶逆变换的频域方法不同,本方法直接在时域中计算有损脉冲响应。结果表明,这种因果脉冲响应对于同时考虑粘性介质中的衍射和二次频率相关衰减的时域计算是理想的。