光声光谱检测装置中光声池的数值计算及优化

利用光声光谱技术进行痕量气体的检测具有独特的优势,光声池是系统装置中最为重要的核心部件,它决定着整机性能的优劣.以一圆柱形共振型光声池为研究对象,基于声学与吸收光谱学的基本理论,建立了光声池声场激发的数学模型;利用数值模拟方法对光声池空腔结构进行了声学模态仿真,获得了前8阶声学模态值以及声压可视化振型;在考虑热黏性声学损耗的作用下,对光声池进行了热-声耦合多物理场仿真计算;将仿真结果与解析计算和实验结果进行对比,明确了利用数值模拟方法来计算光声池有关指标的可靠性与可行性;针对光声池的优化问题,提出了一种将响应面代理模型与遗传算法相结合的优化算法,在将原光声池中的谐振腔两端形貌更改为喇叭口形的情况下,通过优化算法获得了以光声池品质因数Q及池常数C_(cell)为最大值寻优的Pareto最优解集;选取一组解进行考察,结果表明,代理模型预测值与数值模拟值指标最大误差仅为1.3%,优化后的新型光声池Q较之前增长了48.9%, C_(cell)增长了34.4%.研究方法可为光声光谱中光声池的优化设计提供参考借鉴.     

Extension of the angular spectrum method to calculate pressure from a spherically curved acoustic source

The angular spectrum method is an accurate and computationally efficient method for modeling acoustic wave propagation. The use of the typical 2D fast Fourier transform algorithm makes this a fast technique but it requires that the source pressure (or velocity) be specified on a plane. Here the angular spectrum method is extended to calculate pressure from a spherical transducer-as used extensively in applications such as magnetic resonance-guided focused ultrasound surgery-to a plane. The approach, called the Ring-Bessel technique, decomposes the curved source into circular rings of increasing radii, each ring a different distance from the intermediate plane, and calculates the angular spectrum of each ring using a Fourier series. Each angular spectrum is then propagated to the intermediate plane where all the propagated angular spectra are summed to obtain the pressure on the plane; subsequent plane-to-plane propagation can be achieved using the traditional angular spectrum method. Since the Ring-Bessel calculations are carried out in the frequency domain, it reduces calculation times by a factor of approximately 24 compared to the Rayleigh-Sommerfeld method and about 82 compared to the Field II technique, while maintaining accuracies of better than 96% as judged by those methods for cases of both solid and phased-array transducers.     

利用维纳滤波改善声透镜光声成像系统的分辨率

为了克服衍射效应对光声成像系统分辨率的限制,需要采用逆卷积方法进行图像反演.从理论上分析了声透镜成像原理,模拟仿真了声透镜的点扩展函数对声透镜成像系统分辨率的影响和维纳滤波解卷积方法复原光声成像的过程,并利用自搭建的声透镜光声成像系统进行了深入的实验研究,得到了物平面上相距4 mm和3 mm的两个黑胶带点的直接成像光声...     

k-space propagation models for acoustically heterogeneous media: application to biomedical photoacoustics

Biomedical applications of photoacoustics, in particular photoacoustic tomography, require efficient models of photoacoustic propagation that can incorporate realistic properties of soft tissue, such as acoustic inhomogeneities both for purposes of simulation and for use in model-based image reconstruction methods. k-space methods are well suited to modeling high-frequency acoustics applications as they require fewer mesh points per wavelength than conventional finite element and finite difference models, and larger time steps can be taken without a loss of stability or accuracy. They are also straightforward to encode numerically, making them appealing as a general tool. The rationale behind k-space methods and the k-space approach to the numerical modeling of photoacoustic waves in fluids are covered in this paper. Three existing k-space models are applied to photoacoustics and demonstrated with examples: an exact model for homogeneous media, a second-order model that can take into account heterogeneous media, and a first-order model that can incorporate absorbing boundary conditions.     

阶梯复合形光声池声-流特性计算与评估

为了提高光声池的检测性能,提出并分析了一种阶梯复合形光声池.以传统圆柱共振型光声池为基准模型,通过对比解析与模拟计算结果,验证了所采用模拟方法的可靠性,基于模拟方法求解并获得了阶梯腔半径、阶梯腔长度和阶梯腔数量对阶梯复合形光声池声-流特性的影响规律.结果表明:减小阶梯复合形光声池中阶梯腔半径,光声信号相对增强,阶梯腔长...     

Numerical study of friction-induced instability and acoustic radiation – Effect of ramp loading on the squeal propensity for a simplified brake model

This paper presents a numerical study of the influence of loading conditions on the vibrational and acoustic responses of a disc brake system subjected to squeal. A simplified model composed of a circular disc and a pad is proposed. Nonlinear effects of contact and friction over the frictional interface are modelled with a cubic law and a classical Coulomb?s law with a constant friction coefficient. The stability analysis of this system shows the presence of two instabilities with one and two unstable modes that lead to friction-induced nonlinear vibrations and squeal noise. Nonlinear time analysis by temporal integration is conducted for two cases of loadings and initial conditions: a static load near the associated sliding equilibrium and a slow and a fast ramp loading. The analysis of the time responses shows that a sufficiently fast ramp loading can destabilize a stable configuration and generate nonlinear vibrations. Moreover, the fast ramp loading applied for the two unstable cases generates higher amplitudes of velocity than for the static load cases. The frequency analysis shows that the fast ramp loading generates a more complex spectrum than for the static load with the appearance of new resonance peaks. The acoustic responses for these cases are estimated by applying the multi-frequency acoustic calculation method based on the Fourier series decomposition of the velocity and the Boundary Element Method. Squeal noise emissions for the fast ramp loading present lower or higher levels than for the static load due to the different amplitudes of velocities. Moreover, the directivity is more complex for the fast ramp loading due to the appearance of new harmonic components in the velocity spectrum. Finally, the sound pressure convergence study shows that only the first harmonic components are sufficient to well describe the acoustic response.     

Semi-analytical computation of the acoustic field of a segment of a cylindrically concave transducer in lossless and attenuating media

Conventional ultrasound transducers used for medical diagnosis generally consist of linearly aligned rectangular apertures with elements that are focused in one plane. While traditional beamforming is easily accomplished with such transducers, the development of quantitative, physics-based imaging methods, such as tomography, requires an accurate, and computationally efficient, model of the field radiated by the transducer. The field can be expressed in terms of the Helmholtz-Kirchhoff integral; however, its direct numerical evaluation is a computationally intensive task. Here, a fast semianalytical method based on Stepanishen's spatial impulse response formulation [J. Acoust. Soc. Am. 49, 1627-1638 (1971)] is developed to compute the acoustic field of a rectangular element of cylindrically concave transducers in a homogeneous medium. The pressure field, for, lossless and attenuating media, is expressed as a superposition of Bessel functions, which can be evaluated rapidly. In particular, the coefficients of the Bessel series are frequency independent and need only be evaluated once for a given transducer. A speed up of two orders of magnitude is obtained compared to an optimized direct numerical integration. The numerical results are compared with Field II and the Fresnel approximation.     

声纳导流罩吸声障板声学特性的研究

本文研究内容属于舰艇声兼容性技术研究的一部分。应用声学软件SYSNOISE计算了,声纳定向发射时某试验用声纳导流罩及尾部吸声障板的声场特性,同时进行相同工况的水下模型试验。计算结果与试验结果吻合良好,说明应用声学软件对声纳声场进行计算是一种可行的方法,有可能推广至其他的水下模型的声场计算。结果数据反映了吸声障板对声纳定向发射时声场的影响,为舰艇总体声纳导流罩内吸声障板的结构设计、声兼容设计提供依据。     

Fast prediction of aerodynamic noise induced by the flow around a cylinder based on deep neural network

Accurate and fast prediction of aerodynamic noise has always been a research hotspot in fluid mechanics and aeroacoustics. The conventional prediction methods based on numerical simulation often demand huge computational resources, which are difficult to balance between accuracy and efficiency. Here, we present a data-driven deep neural network (DNN) method to realize fast aerodynamic noise prediction while maintaining accuracy. The proposed deep learning method can predict the spatial distributions of aerodynamic noise information under different working conditions. Based on the large eddy simulation turbulence model and the Ffowcs Williams-Hawkings acoustic analogy theory, a dataset composed of 1216 samples is established. With reference to the deep learning method, a DNN framework is proposed to map the relationship between spatial coordinates, inlet velocity and overall sound pressure level. The root-mean-square-errors of prediction are below 0.82 dB in the test dataset, and the directivity of aerodynamic noise predicted by the DNN framework are basically consistent with the numerical simulation. This work paves a novel way for fast prediction of aerodynamic noise with high accuracy and has application potential in acoustic field prediction.     

浅海波导中目标散射的简正波-Kirchhoff近似混合方法

提出了计算浅海波导中复杂目标散射的数值方法:简正波-Kirchhoff近似混合方法。通过把目标散射的Kirchhoff近似方法和简正波声传播模型相结合,可对大尺寸复杂目标在浅海波导中的散射声场进行计算。以浅海波导中刚性球体散射的解析解为标准解验证了本方法,说明简正波-Kirchhoff近似混合方法是有一定计算精度的工程预报方法。数值计算Pekeris浅海波导中球体目标散射声场在深度-频率平面和距离-频率平面上的二维干涉结构及其与自由空间中的差异。进而通过FFT获得目标时域回波随深度的分布图,分析浅海波导中目标姿态、声速剖面对Benchmark潜艇目标散射的影响。      

基于声学扫描振镜的超声/光声双模态成像技术

超声/光声双模态成像技术因其同时兼具超声的高分辨率结构成像和光声的高对比度功能成像优势,极大地推动了光声成像技术的临床应用推广.传统超声/光声双模态成像技术多基于超声成像所用阵列探头同时收集光声信号,系统结构紧凑且无需图像配准,操作便捷.但该类设备使用阵列探头和多通道数据采集,使得其成本较高;且成像结果易受通道一致性差异影响.本文提出了一种基于声学扫描振镜的超声/光声双模态成像技术,该技术采用单个超声换能器结合一维声学扫描振镜进行快速声束扫描,实现超声/光声双模态成像,是一种小型化、低成本的双模态快速成像技术.本文开展了系列仿体和活体成像研究,实验结果表明:系统有效成像范围为15.6 mm,超声和光声成像B扫描速度分别为1.0 s^–1和0.1 s^–1 (光声成像速度主要受制于脉冲激光器重复频率).基于本文所提技术研究,有助于进一步推动超声/光声双模态成像技术的临床转化和普及;也为基于超声信号检测的多模态成像技术提供了一种低成本、小型化和快速声信号检测的参考方案.     

Design, simulation and structural optimization of a longitudinal acoustic resonator for trace gas detection using laser photoacoustic spectroscopy (LPAS)

The design of the acoustic resonator is critical for the optimization of the sensitivity of laser photoacoustic spectroscopy (LPAS) in trace gas detection applications. In this paper, an LC circuit model is used for the simulation of a 1D acoustic resonator. This acoustic resonator is designed for CO photoacoustic spectroscopy. The effects of the structural parameters, quality factor and resonant frequency on the performance of the device are theoretically analyzed. The role of the buffer volume as an acoustic filter is investigated and optimized dimensions of the buffer volume, to achieve minimum noise transmission coefficient, are calculated. The effects of the ambient temperature, variety of pressure and gas flow velocity on the resonant frequency of photoacoustic resonator and PA signal are simulated. The temperature dependence of the microphone sensitivity is also introduced.     

Coupled hydrodynamic-acoustic modeling of sound generated by impacting cylindrical water jets

A coupled hydrodynamic-acoustic model describing acoustic propagation in a fluid containing multiple bubbles is proposed and applied to simulate noise generated by impacting water jets. The total pressure is decomposed into a "hydrodynamic" part and an "acoustic" part and computed using different schemes. The hydrodynamic pressure field is calculated independently using a generalized hydrodynamic model, and the pressure variations serve as sources in the wave equation for the acoustic pressure. A numerical algorithm developed to calculate wave propagation in an irregular region is used to account for the existence of the cavities. Noise generated by the impact of two cylindrical water jets is predicted. The computed near-field pressure is compared with the experimental data.     

板结构的声辐射模态分析及其高分辨率法向振动速度重建

提出了一种基于声辐射模态的速度基向量构建方法,该速度基向量不受网格划分的影响,可用于高分辨率的板结构法向振动速度重建。首先对板表面稀疏网格的声辐射模态进行计算,再以声辐射模态和模态系数构建板法向振动速度分布的基向量,然后由声场测量声压求解基向量系数,最后由该系数和加密网格的速度基向量重建高分辨率的板法向振动速度分布。以简支板声源进行仿真计算,当测量声压信噪比为30 dB时,低频的法向振动速度重建误差最低可达3.7%;以固支板声源在消声室中进行实验验证,131.5 Hz振动频率下的重建误差低于7%。该方法实现了只需要少量声压测量点即可精确重建板声源更高分辨率的法向振动速度分布。      

A new pressure-dominant approximation model for acoustic–structure interaction

A high accuracy approximation modeling approach for the acoustic–structure interaction problem with a shell structure is presented in this paper. The new approximation model aims to accurately reveal the relationship between pressure and velocity in the acoustic field. The main idea of this model is to separate the velocity terms into a combination of velocity and pressure by using a weighting parameter. The modal analysis was performed to find an appropriate weighting parameter for the new model for the spherical case. The stability range of the model is limited during this process. An approximation model was coupled with the equation of motion of a spherical shell to check the performance of the model. Responses of a spherical shell excited by a plane step wave and cosine-type incident pressure from the new model were compared to the exact solution and solutions from former approximation models such as Doubly Asymptotic Approximations (DAAs). The new proposed model can approximate high accuracy responses in both early and late time.     

基于时间反转法的相控换能器声场的仿真研究

相控换能器具有焦距可调的优势。本文以82阵元相控换能器建立的3D数值仿真模型为例,基于时间反转法提取阵元的激励信号,利用时域有限差分(FDTD)法对Westervelt声波非线性传播方程进行声场数值仿真,研究不同阵元分布、偏离声轴的距离、设定焦距大小对形成声场的影响,可调控范围及其消除旁瓣方法。研究结果表明,随机分布相控阵可明显降低声场中的旁瓣;随着偏离声轴距离的增加,主瓣声压幅值逐渐减小,旁瓣与主瓣的最大声强比值r逐渐增大,且沿声轴的可调控范围逐渐减小;随声轴方向上设定焦距的增加,主瓣声压幅值先增大后减小,r值先减小后增大;基于时间反转法的高声压旁瓣消除法可在一定程度上扩大相控阵声场的可调控范围。     

Photoacoustic investigations of the optical absorption spectra of porous silicon layers on a silicon backing

This paper presents a series of experimental photoacoustic spectra of porous silicon layers on crystalline silicon and their numerical analysis performed in the proposed two-layer model. The goal of the analysis was to calculate the optical absorption spectra of porous silicon from the photoacoustic spectra of porous silicon layers on a silicon background. The experimental character of the observed absorption band associated with the porous silicon was revealed. This is the first attempt at a theoretical interpretation of the photoacoustic spectra of porous silicon on a silicon backing.     

一种预测扩散声场在非规则形状刚性壁面附近干涉图样的方法

扩散声场会在反射边界附近形成干涉图样,研究方法包括平面波模型、简正模态分析、渐进模态分析等,但仅适用于尺度远大于声波波长的矩形声腔。提出一种预测扩散声场在非规则刚性壁面结构附近形成的干涉图样的数值方法,表明结构附近“受挡”声压的互谱矩阵取决于:(1)假定该结构在自由空间中振动辐射声音时其表面法向振速到表面及场点声压的边界元系数矩阵;(2)假定结构置于自由空间中且表面刚性时,点声源辐射声波入射到结构表面上产生的散射声场的边界元系数矩阵;(3)扩散声场均方声压。仿真表明,该途径预测的干涉图样与理论值完全吻合。该预测方法还可用于混响环境下声源附近直达声压均方值的空间分布估计,为混响环境下设备的声源定位提供帮助。      

Highly directional acoustic receivers

The theoretical directivity of a single combined acoustic receiver, a device that can measure many quantities of an acoustic field at a collocated point, is presented here. The formulation is developed using a Taylor series expansion of acoustic pressure about the origin of a Cartesian coordinate system. For example, the quantities measured by a second-order combined receiver, denoted a dyadic sensor, are acoustic pressure, the three orthogonal components of acoustic particle velocity, and the nine spatial gradients of the velocity vector. The power series expansion, which can be of any order, is cast into an expression that defines the directivity of a single receiving element. It is shown that a single highly directional dyadic sensor can have a directivity index of up to 9.5 dB. However, there is a price to pay with highly directive sensors; these sensors can be significantly more sensitive to nonacoustic noise sources.     

Linear time domain model of the acoustic potential field

A new time domain formulation of the acoustic wave is developed to avoid approximating assumptions of the linearized scalar wave equation that limit its validity to low Mach particle velocity modeling or to a smooth potential field in a stationary medium. The proposed model offers precision of the moving frame while retaining the form of the widely used linearized scalar wave equation although with respect to modified coordinates. It is applicable to field calculations involving transient waves with unlimited particle velocity, propagating in inhomogenous fluids or in those with time varying density. The model is based on the exact flux continuity equation and the equation of motion, both using the moving reference frame. The resulting closed-form free space scalar wave equation employing total derivatives is converted back to the partial differential form by using modified independent variables. The modified variables are related to the common coordinates of space and time following integral expressions involving transient particle velocity representing wave radiated by each point of a stationary source. Consequently, transient field produced by complex surface velocity sources can be calculated following existing surface integrals of the radiation theory although using modified coordinates. The use of the proposed model is presented in a numerical simulation of a transient velocity source vibrating at selected magnitudes, leading to the determination of the propagating pressure and velocity wave at any point.