Elastic guided wave propagation in a periodic array of multi-layered piezoelectric plates with finite cross-sections

In this study, we present a model study of guided wave dispersion and resonance behavior of an array of piezoelectric plates with arbitrary cross-sections. The objective of this work is to analyze the influence of the geometry of an element of a 1D-array ultrasound transducer on generating multi-resonance frequency so as to increase the frequency bandwidth of the transducers. A semi-analytical finite-element (SAFE) method is used to model guided wave propagation in multi-layered 1D-array ultrasound transducers. Each element of the array is composed of LiNbO₃ piezoelectric material with rectangular or subdiced cross-section. Four-node bilinear finite-elements have been used to discretize the cross-section of the transducer. Dispersion curves showing the dependence of phase and group velocities on the frequency, and mode shapes of propagating modes were obtained for different geometry consurations. A parametric analysis was carried out to determine the effect of the aspect ratio, subdicing, inversion layer and matching layers on the vibrational behavior of 1D-array ultrasound transducers. It was found that the geometry with subdiced cross-section causes more vibration modes compared with the rectangular section. Modal analysis showed that the additional modes correspond to lateral modes of the piezoelectric subdiced section. In addition, some modes have strong normal displacements, which may influence the bandwidth and the pressure field in front of the transducer. In addition, the dispersion curves reveal strong coupling between waveguide modes due to the anisotropy of the piezoelectric crystal. The effect of the matching layers was to cluster extensional and flexural modes within a certain frequency range. Finally, inversion layer is found to have a minor effect on the dispersion curves. This analysis may provide a means to analyze and understand the dynamic response of 1D-array ultrasound transducers.     

Modeling of three-dimensional pressure fields in sonochemical reactors with an inhomogeneous density distribution of cavitation bubbles. Comparison of theoretical and experimental results

During the last 50 years extensive experimental investigation has been carried out on the chemical effects of ultrasound, but limited work has been reported on modeling. This paper presents a new model in which a numerical calculation of the three-dimensional linear sound pressure field distribution in a commonly used sonoreactor containing three transducers is carried out. In this model the inhomogeneous three-dimensional time-dependent wave equation was solved using the finite difference approach. The modeled results are then compared with the experimentally measured values, and the agreement, in general, is found to be good. Further, our modeling studies have an advantage, since they clearly describe the continuous sound pressure field structure, unlike previously reported results in which some information is missing due to limited intermittent measured points.     

Ultrasound fields in attenuating media

For medical ultrasonic imaging and for nondestructive testing, the attenuation of pressure waves and the resulting shift in wave velocity are important features in commonly used transmission media such as biological tissue. An algorithm for the numerical evaluation of pressure field distributions generated by ultrasonic transducers is presented. The attenuation and dispersion of the sound transmission medium are taken into consideration. The sound fields are computed numerically for continuous wave as well as pulse excitation. The transducer has plane or gently curved geometry and is embedded in a plane rigid baffle. The numerically determined pressure fields are presented as 3D plots, as gray-scale images for a fixed time stamp (like a snapshot), or as isobars regarding the maximum values over time for each local point in the area under investigation. The algorithm described here can be utilized as a tool for design of ultrasound transducers, especially array antennas.     

高频聚焦超声声场和温度场的仿真研究

为探究临床常用的7 MHz高频聚焦超声在多层生物组织中的声传播以及毫秒级时间内的生物传热规律问题,基于Westervelt方程和Pennes传热方程,使用有限元方法建立高频聚焦超声辐照多层组织的非线性热黏性声传播及传热模型。首先分析了线性模型和非线性模型之间的差异,然后在非线性模型下探究换能器的参数对声场和温度场的影响。仿真结果显示：在7 MHz频率下,当换能器输出声功率超过5 W时,声波传播的非线性效应不可忽视(p <0.05);当声功率从5 W增大到15 W时,非线性模型与线性模型预测的温度偏差从20%增加到34.703%;高频聚焦超声波的非线性行为比低频更加显著,基频能量向高次谐波转移的程度增大,声功率为10 W和15 W时4次谐波与基波之比分别达到7.33%和12.12%;高频换能器参数的改变对组织中声场和温度场分布的影响较大,换能器焦距从12 mm减小到11.2 mm,焦点处最高温度增加了77%。结果表明,7 MHz聚焦超声的非线性声传播需要考虑到4次谐波的影响。该文提出的多层组织非线性仿真模型可为高频聚焦超声换能器参数优化及制定安全、有效的术前治疗方案提供理论参考。     

剪切波对HIFU经颅聚焦形成温度场影响的数值仿真研究

在高强度聚焦超声经颅治疗时,既有纵波又有剪切波,为了保障该治疗方法的安全有效性,有必要分析剪切波对HIFU治疗温度场的影响。该文基于人体头颅CT数据和曲率半径为150 mm的256阵元的半球相控换能器建立三维高强度聚焦超声经颅声波传播模型,利用时域有限差分法结合Westervelt声波非线性传播方程、动量方程、质量守恒方程和Pennes生物热传导方程数值仿真其形成温度场,研究在相同输入功率、不同聚焦角度条件下对应阵元数进行激励时,剪切波对换能器形成温度场的影响。结果表明,随换能器聚焦角度减小,在几何焦点处形成的焦域面积逐渐增大,考虑剪切波形成的温度场达到65?C所需时间逐渐延长,焦点前移程度越大;在相同聚焦角度条件下,考虑剪切波的温度场达到65?C所需时间更短,旁瓣更少,在颅骨处的温度更高,对焦点前移几乎没有影响;随换能器聚焦角度减小,考虑剪切波的模型形成的焦域面积变化范围更大;幂指数函数形式对不同聚焦角度下焦域面积大小的拟合优度高,可预测不同聚焦角度换能器形成的焦域面积。     

Sound field prediction of ultrasonic lithotripsy in water with spheroidal beam equations

With converged shock wave,extracorporeal shock wave lithotripsy(ESWL)has become a preferable way to crush human calculi because of its advantages of efficiency and non-intrusion.Nonlinear spheroidal beam equations(SBE)are employed to illustrate the acoustic wave propagation for transducers with a wide aperture angle.To predict the acoustic field distribution precisely,boundary conditions are obtained for the SBE model of the monochromatic wave when the source is located on the focus of an ESWL transducer.Numerical results of the monochromatic wave propagation in water are analyzed and the influences of half-angle,fundamental frequency,and initial pressure are investigated.According to our results,with optimization of these factors,the pressure focal gain of ESWL can be enhanced and the effectiveness of treatment can be improved.     

Exact wave field simulation for finite-volume scattering problems

An exact boundary condition is presented for scattering problems involving spatially limited perturbations of arbitrary magnitude to a background model in generally inhomogeneous acoustic media. The boundary condition decouples the wave propagation on a perturbed domain while maintaining all interactions with the background model, thus eliminating the need to regenerate the wave field response on the full model. The method, which is explicit, relies on a Kirchhoff-type integral extrapolation to update the boundary condition at every time step of the simulation. The Green's functions required for extrapolation through the background model are computed efficiently using wave field interferometry.     

Finite element simulation of the generation and detection by air-coupled transducers of guided waves in viscoelastic and anisotropic materials

The measured characteristics (efficiency and sensitivity) of two air-coupled transducers allow for the prediction of the absolute values of the pressure of the bulk waves generated in air and for the measurement of the pressure of the field radiated in air by guided waves propagating in a structure. With finite element software, the pressure field generated by an air-coupled transducer is simulated by introducing a right-hand side member in the Helmholtz equation, which is used for computing the propagation from the transducer to a plate. The simulated source is rotated in order to impose an angle of incidence with respect to the normal of the plate and generate the corresponding guided mode. Inside the plate, the propagation is simulated with the dynamic equations of equilibrium and a complex stiffness tensor to take into account the viscoelastic anisotropy of the material. For modeling the three-dimensional fields of the guided modes propagating in a two-dimensional non-symmetry plane, a 2.5 dimensional model is introduced. The model computes the value of the pressure field radiated in air by the plates for any guided modes and can predict the detectability of the system for a known defect in a structure. A test bed incorporating two air-coupled transducers is used to generate and receive various guided modes. Two plates made of Perspex and carbon-epoxy composite are tested. The pressure measured by the receiver at various positions is compared to the results of the model to validate it.     

Modeling the propagation of nonlinear three-dimensional acoustic beams in inhomogeneous media

A three-dimensional model of the forward propagation of nonlinear sound beams in inhomogeneous media, a generalized Khokhlov-Zabolotskaya-Kuznetsov equation, is described. The Texas time-domain code (which accounts for paraxial diffraction, nonlinearity, thermoviscous absorption, and absorption and dispersion associated with multiple relaxation processes) was extended to solve for the propagation of nonlinear beams for the case where all medium properties vary in space. The code was validated with measurements of the nonlinear acoustic field generated by a phased array transducer operating at 2.5 MHz in water. A nonuniform layer of gel was employed to create an inhomogeneous medium. There was good agreement between the code and measurements in capturing the shift in the pressure distribution of both the fundamental and second harmonic due to the gel layer. The results indicate that the numerical tool described here is appropriate for propagation of nonlinear sound beams through weakly inhomogeneous media.     

格子玻尔兹曼方法对不同张角聚焦声束的建模

高强度聚焦超声(HIFU)是一种新型的无创治疗肿瘤新技术,其中换能器声场数值计算能够为HIFU治疗提供重要的依据。传统非线性KZK和SBE模型广泛应用于换能器声场数值计算,但依然存在某些不足。我们采用一种介观尺度的新型流体力学方法,即格子Boltzmann方法(LBM),基于2维9离散速度(D2Q9)格子构建了轴对称多弛豫参数LBM模型,并通过调节弛豫参数分析其对模型的影响;利用该模型对两个具有不同张角的球面聚焦换能器的声场进行数值模拟,并与KZK和SBE模型的计算结果进行比较。结果表明LBM模型能够很好地描述超声波的激发和传播机制,从流体力学的角度描述聚焦声场的分布,具有清晰的物理意义,且计算过程不受换能器张角的限制,在换能器声场的理论分析和模拟计算及其在HIFU治疗中的应用有着积极的意义。      

Numerical simulations of partial elements excitation for hemispherical high-intensity focused ultrasound phased transducer

The hemispherical phased transducer maximizes the coverage of the skull and the ultrasonic energy per unit area of the skull is minimized, thereby reducing the risk of skull burns, but the transducer has a small focal area adjustment range, increasing the focal length of treatment is an urgent question for this type of transducer. In this paper, a three-dimensional high-intensity focused ultrasound (HIFU) transcranial propagation model is established based on the human head structure. The finite difference time domain (FDTD) is combined with the Westervelt acoustic wave nonlinear propagation equation and Penne's biological heat conduction equation for numerical simulation of the sound pressure field and temperature field. Forming a treatable focal area in a small-opening hemispherical transducer with a small amount of numerical simulation calculation focusing at a set position to determine the minimum partial excitation area ratio of focusing. And then, applying these preliminary results to a large-opening diameter hemispherical transducer and the temperature field formed by it or full excitation is studied. The results show that the focus area with the excitation area ratio of less than 22% moves forward to the transducer side when the excitation sound is formed. When the excitation area ratio is greater than or equal to 23%, it focuses at the set position. In the case of partial incentives, using 23% of the partial array, the adjustable range of the treatable focal area formed in the three-dimensional space is larger than that of the full excitation.     

Model for continuously scanning ultrasound vibrometer sensing displacements of randomly rough vibrating surfaces

An analytic model is developed for the time-dependent ultrasound field reflected off a randomly rough vibrating surface for a continuously scanning ultrasound vibrometer system in bistatic configuration. Kirchhoff's approximation to Green's theorem is applied to model the three-dimensional scattering interaction of the ultrasound wave field with the vibrating rough surface. The model incorporates the beam patterns of both the transmitting and receiving ultrasound transducers and the statistical properties of the rough surface. Two methods are applied to the ultrasound system for estimating displacement and velocity amplitudes of an oscillating surface: incoherent Doppler shift spectra and coherent interferometry. Motion of the vibrometer over the randomly rough surface leads to time-dependent scattering noise that causes a randomization of the received signal spectrum. Simulations with the model indicate that surface displacement and velocity estimation are highly dependent upon the scan velocity and projected wavelength of the ultrasound vibrometer relative to the roughness height standard deviation and correlation length scales of the rough surface. The model is applied to determine limiting scan speeds for ultrasound vibrometer measuring ground displacements arising from acoustic or seismic excitation to be used in acoustic landmine confirmation sensing.     

Wide-Angle Parabolic Approximation for Modeling High-Intensity Fields from Strongly Focused Ultrasound Transducers

A novel numerical algorithm based on the wide-angle parabolic approximation is developed for modeling linear and nonlinear fields generated by axially symmetric ultrasound transducers. An example of a strongly focused single-element transducer is used to compare the results of ultrasound field simulations based on the Westervelt equation, Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation with differently modified boundary condition, and nonlinear wide-angle parabolic equation. It is demonstrated that having a computational speed comparable to modeling the KZK equation, the use of wide-angle parabolic approximation makes it possible to obtain solutions for highly focused ultrasound beams that are closer in accuracy to solutions based on the Westervelt equation.     

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.     

Multiple scattering using the Foldy-Twersky integral equation

The importance of multiple scattering processes for the propagation of ultrasound through dispersions is assessed. Under certain well defined conditions, the Foldy-Twersky integral equation can be used to find a value for the amplitude and complex wavenumber of the multiply scattered plane wave collected by either the receiving or the emitting transducer. Expressions for the amplitude and wavenumber are evaluated in terms of particle number density, wavelength of ultrasound and the forward and backward scattering amplitudes. In the case where the wavelength is much greater than the particle size, explicit expressions for the real and imaginary parts of the wavenumber are provided, in terms of the scattering and absorption cross-sections, as far as terms in the square of the number density of particles.     

The development of a shock-tube based characterization technique for air-coupled ultrasonic probes

The present paper proposes a new characterization technique for air-coupled ultrasound probes. The technique is based on a shock tube to generate a controlled pressure wave to calibrate transducers within their operating frequency range. The aim is to generate a high frequency pressure wave (at least up to 200 kHz) with the low energy levels typical of commonly used air-coupled ultrasound probes. A dedicated shock-tube has been designed and tested to assess calibration performances. The sensor transfer function has been measured by using a pressure transducer as reference.     

Sound propagation and scattering in media with random inhomogeneities of sound speed, density and medium velocity

This review presents both classical and new results of the theory of sound propagation in media with random inhomogeneities of sound speed, density and medium velocity (mainly in the atmosphere and ocean). An equation for a sound wave in a moving inhomogeneous medium is presented, which has a wider range of applicability than those used before. Starting from this equation, the statistical characteristics of the sound field in a moving random medium are calculated using Born-approximation, ray, Rytov and parabolic-equation methods, and the theory of multiple scattering. The results obtained show, in particular, that certain equations previously widely used in the theory of sound propagation in moving random media must now be revised. The theory presented can be used not only to calculate the statistical characteristics of sound waves in the turbulent atmosphere or ocean but also to solve inverse problems and develop new remote-sensing methods. A number of practical problems of sound propagation in moving random media are listed and the further development of this field of acoustics is considered.     

A reflectivity equation for acoustic media with absorption

The wave equation which governs the propagation of dilatational stress (pressure) in viscoacoustic media can be reduced to the lossless acoustic wave equation by the introduction of a complex traveltime coordinate. An interesting property of this equation is that it contains a term which is proportional to the reflectivity function in its single scattering approximation.     

Study of the second harmonics in nonlinear scattering of interacting acoustic waves by a rigid cylinder

The field formed by the second harmonics in the case of the nonlinear scattering of interacting plane acoustic waves by a rigid cylinder is considered. The method of successive approximations is used to obtain the solutions to the inhomogeneous wave equation in the first and second approximations. Asymptotic expressions are derived for the components of the total acoustic pressure of the second harmonics, and the scattering diagrams for these components are presented.     

热粘弹波在变温非均匀合金熔体中的传播

机械波在金属凝固过程中传播的定量计算一直是一个难题,主要原因就是在这个过程中的熔体结构非常复杂.本研究考虑到熔体的变温、非均匀和粘弹性的特点,采用Kelvin粘弹性介质模型,建立了具有粘热损失特性的热粘弹性波动方程,通过隐式有限差分方法对波动方程进行求解,并以ZL203A合金熔体为研究对象,探究了热粘弹波在变温非均匀介质中的传播规律.结果表明:热粘弹波从合金熔体的低温区向高温区传播时,非均匀的温度场对波的传播有较大影响;相反,当波从合金熔体的高温区向低温区传播时,非均匀的温度场对波的传播几乎没有影响.热粘弹波在合金熔体中的衰减系数随频率的增大呈线性增大,而随温度的升高先增大后减小,在熔体的枝晶搭接温度附近达到最大值.