高斯波束对界面附近离轴球形粒子的轴向声辐射力

从声波的散射理论出发，利用级数展开法得到高斯波束的波束因子，推导其对阻抗边界下离轴球形粒子声辐射力.针对刚性球与液体球两种球形粒子进行数值模拟，与自由空间的情况进行比较.讨论边界反射系数、粒子与边界距离、束腰半径以及离轴角度与距离等对声辐射力的影响.仿真结果表明：边界反射系数的增大会引起声辐射力的增加，但不改变峰值的位置；在合适的频率处，可以产生负向声辐射力；声辐射力随粒子与边界距离呈周期性变化；束腰半径的影响主要体现在中高频；随着粒子偏离传播轴的距离和角度增大，声辐射力明显衰减.该研究为利用高斯波束实现对粒子的操纵提供理论基础.     

零阶Bessel驻波场中任意粒子声辐射力和力矩的Born近似

声辐射力和声辐射力矩的计算是实现粒子精准操控的重要基础.基于经典声散射理论的偏波级数展开法较难直接用于复杂模型的研究,而纯数值的方法则不利于进行系统的参数化分析.基于Born近似的基本原理,推导了低频情况下零阶Bessel驻波场中心任意粒子的声辐射力和力矩表达式.在此基础上,以球形粒子、椭球形粒子和柱形粒子为例进行详细地计算,并考虑声参数的非均匀性对声辐射力和力矩的影响.仿真结果表明,在低频范围内Born近似具有很高的精度,随着频率的增加和粒子与流体的阻抗匹配变差,Born近似的精度逐渐下降.对于倾斜放置于零阶Bessel驻波场中的椭球形粒子和柱形粒子,非对称性会导致其受到声辐射力矩的作用.在粒子尺寸远小于波长的情况下,声辐射力特性与粒子的具体形状几乎无关,但声辐射力矩不然.最后,引入周围流体的黏滞效应并对声辐射力的表达式进行了修正.该研究预期可以为生物医学、材料科学等领域利用驻波场声镊子实现微小粒子的精准操控提供一定的理论指导.     

声空化气泡成长及破裂研究

本文对液体内的声空化气泡的成长与破裂过程进行数值计算,得到各种情况下气泡壁的运动情况.通过对不同初始半径、不同频率下声空化气泡运动的计算,得到空化气泡半径小于共振半径,可以增强空化效果,而单一的增强声场的频率并不一定能加强声空化效果,为增强空化效果提供理论依据.研究各种信号作用下声空化气泡成长情况,明确方波信号激励下的...     

不同光聚焦状态的光击穿声辐射特性

强激光通过光学系统聚焦于水下时,击穿水介质辐射强声波信号.光学聚焦特性的不同,使得激光击穿区域形状的差异,导致辐射的声信号在强度、频谱特征上具有较大的差异性.为研究不同光聚焦状态对激光击穿形成的空泡辐射声波的影响,推导了不同聚焦状态下激光空泡壁的运动方程及声辐射模型,构建了激光声实验测量系统,实验研究了不同聚焦状态下光击穿导致的空泡脉动参数、辐射声信号特征的差异性.结论:激光声信号强度与激光能量成线性关系,在激光能量高于一定参数情况下,声信号强度变化不大;在低激光功率下,扩束聚焦与非扩束聚焦,辐射的声信号特性差异不大;在高激光功率条件,扩束聚焦击穿形成的空泡半径大于非扩束聚焦击穿,辐射的声信号强度高于非扩束聚焦;为提高激光击穿辐射卢信号效率,应对激光束进行扩束聚焦.     

Bessel驻波对椭球形液滴的声辐射转矩

在实际的声操控中,由于声辐射力、表面张力和重力的共同作用,液滴往往呈现出椭球的形状,在螺旋声场中会受到力矩的作用而发生转动。从声波的散射理论出发,根据部分波展开法求解得到了椭球形液滴在Bessel驻波场中的声散射系数,并给出了其受到的声辐射转矩的解析式。在此基础上,对椭球形不可压缩液滴和椭球形可压缩液滴分别进行数值计算。仿真结果表明,不可压缩液滴的声辐射转矩与声束半锥角的关系更密切,而可压缩液滴则更依赖于特定的频率;提升Bessel驻波场的阶数有利于增强声辐射转矩的峰值,但在中低频处较难对可压缩液滴产生明显的力矩。该研究结果预期对利用螺旋声场进行液滴的操控具有理论指导作用。      

改进的半空间频率均方声压法计算结构频带振动声辐射

为计算中高频半空间结构频带振动下声辐射问题,推导了基于能量源的半空间频率均方声压法(Frequency Averaged Quadratic Pressure,FAQP),并提出了改进的半空间FAQP法,克服了半空间FAQP法在失效频率下解不唯一的问题。频率带宽为4 Hz的刚性壁面作用下的脉动球和与自由表面相接触的脉动半球的声辐射数值计算,验证了改进的半空间FAQP法对解决失效频率下解不唯一问题的有效性。同时,刚性壁面作用下的脉动球和与自由表面相接触的脉动半球声辐射数值计算结果均表明,在1/3倍频程下,改进的半空间FAQP法与常规边界元方法相比,具有更高的计算精度,更适用于中高频计算。      

基于主动声学超材料的圆柱声隐身斗篷设计研究

基于多层复合材料结构的二维声隐身斗篷设计思想, 利用主动隔膜声学空腔有效密度可以任意控制这一特性, 设计了主动声学超材料下的无限长圆柱声隐身斗篷. 给出了主动隔膜声学空腔单元的声电元件类比模拟电路图和具体的有效密度控制方法. 进行了主动声学超材料声隐身斗篷的结构建模, 并对平面入射波入射下此圆柱隐身斗篷周围声压分布场进行仿真计算. 结果表明, 平面波在一定频率范围内可以毫无阻碍地透过圆柱斗篷, 似乎不存在这种障碍物, 达到声隐身效果. 同时, 计算了主动声材料斗篷下总散射截面随频率变化曲线, 研究了此斗篷隐身效果随频率的变化特性. 本文从主动控制角度探讨实验实现隐身斗篷的技术问题, 有望给声隐身斗篷实验设计提供一条新的技术途径.     

三维浅海下弹性结构声辐射预报的有限元-抛物方程法

利用多物理场耦合有限元法对结构和流体适应性强、抛物方程声场计算高效准确的特点,提出了三维浅海波导下弹性结构声振特性研究的有限元-抛物方程法.该方法采用多物理场耦合有限元理论建立浅海下结构近场声辐射模型,计算局域波导下结构声振信息,并提取深度方向上复声压值作为抛物方程初始值;然后采用隐式差分法求解抛物方程以步进计算结构辐射声场.重点介绍了该方法对浅海下结构声辐射计算的准确性、高效性以及快速收敛性后,对Pekeris波导中有限长弹性圆柱壳的声振特性进行了分析.研究得出,当圆柱壳靠近海面(海底)时,其耦合频率比自由场下的要高(低),当潜深达到一定范围时,与自由场耦合频率基本趋于一致;在低频远场,结构辐射场与同强度点源声场具有一定的等效性,且等效距离随着频率增加而增加;由于辐射声场受结构振动模态、几何尺寸和简正波模式影响,结构辐射场传播的衰减规律按近场声影响区、球面波衰减区、介于球面波和柱面波衰减区、柱面波衰减区四个扩展区依次进行.     

透明液体中激光致声特性

研究透明液体的特性阻抗、激光声的传播距离和激光焦斑半径对激光声信号的影响.通过实验和数值计算发现,在一定范围内,随着特性阻抗的增加,声波的峰值声压、主频和声能也相应增加;随着激光声传播距离的增加,激光声的主频从急速下降变为相对稳定;在保证光击穿的条件下,选择较大的激光焦斑半径,可提高光声能量转换效率.     

椭球三维液态声子晶体完全声波禁带研究

弹性波在声子晶体中传播时，会产生一定的频率禁带，声和振动在禁带中被禁止传播。本文将椭球散射体引入三维声子晶体中，采用平面波展开法计算了该系统的声波禁带结构。对于不同椭球半径的系统，其声波禁带的位置与大小有很大不同。结果表明：当晶格常数和填充率确定时，禁带受到椭球半径的影响；当确定椭球体的某一个半径和填充率时，第一禁带的最大值出现在另2个半径相等的情况下。     

Acoustic scattering by multiple elliptical cylinders using collocation multipole method

This paper presents the collocation multipole method for the acoustic scattering induced by multiple elliptical cylinders subjected to an incident plane sound wave. To satisfy the Helmholtz equation in the elliptical coordinate system, the scattered acoustic field is formulated in terms of angular and radial Mathieu functions which also satisfy the radiation condition at infinity. The sound-soft or sound-hard boundary condition is satisfied by uniformly collocating points on the boundaries. For the sound-hard or Neumann conditions, the normal derivative of the acoustic pressure is determined by using the appropriate directional derivative without requiring the addition theorem of Mathieu functions. By truncating the multipole expansion, a finite linear algebraic system is derived and the scattered field can then be determined according to the given incident acoustic wave. Once the total field is calculated as the sum of the incident field and the scattered field, the near field acoustic pressure along the scatterers and the far field scattering pattern can be determined. For the acoustic scattering of one elliptical cylinder, the proposed results match well with the analytical solutions. The proposed scattered fields induced by two and three elliptical–cylindrical scatterers are critically compared with those provided by the boundary element method to validate the present method. Finally, the effects of the convexity of an elliptical scatterer, the separation between scatterers and the incident wave number and angle on the acoustic scattering are investigated.     

Acoustic radiation force on a rigid cylinder near rigid corner boundaries exerted by a Gaussian beam field

Acoustic manipulation is one of the well-known technologies of particle control and a top research in acoustic field. Calculation of acoustic radiation force on a particle nearby boundaries is one of the critical tasks, as it approximates realistic applications. Nevertheless, it is quite difficult to solve the problem by theoretical method when the boundary conditions are intricate. In this study, we present a finite element method numerical model for the acoustic radiation force exerting on a rigid cylindrical particle immersed in fluid near a rigid corner. The effects of the boundaries on acoustic radiation force of a rigid cylinder are analyzed with particular emphasis on the non-dimensional frequency and the distance from the center of cylinder to each boundary. The results reveal that these parameters play important roles in acoustic manipulation for particle-nearby complicated rigid boundaries. This study verifies the feasibility of numerical analysis on the issue of acoustic radiation force calculation close to complex boundaries, which may provide a new idea on analyzing the acoustic particle manipulation in confined space.     

水下任意刚性散射体对Bessel波的散射特性分析

本文利用Bessel波的谐波展开式, 采用T矩阵方法的推导思路, 建立了水下任意刚性散射体在Bessel波照射下的声散射场计算公式. 以水下刚性椭球体和两端附连半球的刚性圆柱体为例, 计算了在不同波锥角β 下的反向散射形态函数, 同时, 依据镜反射波和绕行波的干涉物理模型, 给出了预报Bessel波照射下的反向散射形态函数峰峰间隔值的计算模型. 仿真结果表明本文提出的Bessel波照射下反向散射形态函数峰峰间隔值预报方法是准确有效的, 同时也说明, 本文建立的基于T矩阵法计算水下任意刚性散射体在Bessel波束下的声散射场方法是有效的, 这拓展了T矩阵法的应用领域.     

Theoretical and numerical calculations for the time-averaged acoustic force and torque acting on a rigid cylinder of arbitrary size in a low viscosity fluid

In this paper, theoretical calculations as well as numerical simulations are performed for the time-averaged acoustic force and torque on a rigid cylinder of arbitrary size in a fluid with low viscosity, i.e., the acoustic boundary layer is thin compared to the cylinder radius. An exact analytical solution and its approximation are proposed in the form of an infinite series including Bessel functions. These solutions can be evaluated easily by a mathematical software package such as mathematica and matlab. Three types of incident waves, plane traveling wave, plane standing wave, and dual orthogonal standing waves, are investigated in detail. It is found that for a small particle, the viscous effects for an incident standing wave may be neglected but those for an incident traveling wave are notable. A nonzero viscous torque is experienced by the rigid cylinder when subjected to dual orthogonal standing waves with a phase shift even when the cylinder is located at equilibrium positions without imposed acoustic forces. Furthermore, numerical simulations are carried out based on the FVM algorithm to verify the proposed theoretical formulas. The theoretical results and the numerical ones agree with each other very well in all the cases considered.     

Axial acoustic radiation force of progressive cylindrical diverging waves on a rigid and a soft cylinder immersed in an ideal compressible fluid

Background and motivation

Previous works investigating the radiation force of diverging spherical progressive waves incident upon spherical particles have demonstrated the direction of reversal of the force when the particle is subjected to a curved wave-front. In this communication, the analysis is extended to the case of diverging cylindrical progressive waves incident upon a rigid or a soft cylinder in a non-viscous fluid with explicit calculations for the radiation force function (which is the radiation force per unit energy density and unit cross-sectional surface) not shown in [F.G. Mitri, Ultrasonics 50 (2010) 620-627].

Method

A closed-form solution presented previously in [F.G. Mitri, Ultrasonics 50 (2010) 620-627] is used to plot the radiation force function with particular emphasis on the difference from the results of incident plane progressive waves versus the size parameter ka (k is the wave number and a is the cylinder’s radius) and the distance of the cylinder from the acoustic source r₀.

Results

Radiation force function calculations for the rigid cylinder unexpectedly reveal that under specific conditions determined by the frequency of the acoustic field, the radius of the cylinder, as well as the distance to the acoustic source, the force becomes attractive (negative force). In addition, the numerical results show that the radiation force on a rigid cylinder does not generally obey the inverse-distance law with respect to the distance from the source.

Conclusion and potential applications

These results suggest that it may be possible, under specific conditions, to pull a cylindrical structure back toward the acoustic source using progressive cylindrical diverging waves. They may also provide a means to predict the radiation force required to manipulate non-destructively a single cylindrical structure. Potential applications include the design of a new generation of acoustic tweezers operating using a single beam of progressive waves (in contrast to the traditional version of acoustical tweezers in which an acoustic standing wave field is produced using two counter-propagating acoustic fields) for investigations in the field of flow cytometry, particle manipulation and entrapment.     

Pseudo-Gaussian cylindrical acoustical beam – Axial scattering and radiation force on an elastic cylinder

Making use of the addition theorem for the cylindrical wave functions and the complex-source-point method in cylindrical coordinates, an exact solution to the Helmholtz equation is derived, which corresponds to a tightly focused (or collimated) cylindrical quasi-Gaussian beam with arbitrary waist. The solution is termed “quasi-Gaussian” to make a distinction from the standard Gaussian beam solution obtained in the paraxial approximation. The advantage of introducing this new solution is the efficient and fast computational modeling of tightly focused or quasi-collimated cylindrical wave-fronts depending on the dimensionless waist parameter kw₀, where k is the wavenumber of the acoustical radiation. Moreover, a closed-form partial-wave series expansion is obtained for the incident field, which has the property that the axial scattering (i.e. along the direction of wave propagation) and the axial acoustic radiation force (which is a time-averaged quantity) on a cylinder, can be calculated without any approximations in the limit of linear acoustical waves in a nonviscous fluid. Examples are found where the extinction in the radiation force function plot is found to be correlated with conditions giving reduction of the backscattering from an elastic cylinder. Those results are useful in beam-forming design, particle manipulation in acoustic tweezers operating with focused cylindrical beams, and the prediction of the scattering and radiation forces on a cylindrical particle or liquid bridges.     

Acoustic elliptical cylindrical cloaks

By making a comparison between the acoustic equations and the 2-dimensional (2D) Maxwell equations, we obtain the material parameter equations (MPE) for acoustic elliptical cylindrical cloaks. Both the theoretical results and the numerical results indicate that an elliptical cylindrical cloak can realize perfect acoustic invisibility when the spatial distributions of mass density and bulk modulus are exactly configured according to the proposed equations. The present work is the meaningful exploration of designing acoustic cloaks that are neither sphere nor circular cylinder in shape, and opens up possibilities for making complex and multiplex acoustic cloaks with simple models such as spheres, circular or elliptic cylinders.     

椭圆柱体对任意入射极化的电磁波的散射

基于电磁场的多尺度变换理论,得到了以导体椭圆柱为例的电磁波任意极化,任意垂直入射到目标上的散射场的解析式.将所得结果应用于计算圆柱目标的散射场,结果与文献完全一致.对椭圆柱体的散射宽度随入射波方位角,电磁波频率以及目标尺寸的变化进行了仿真.结果表明,在垂直于电磁波的方向上观测时,椭圆柱的垂直尺寸对散射有较大的影响,电磁波水平极化时散射最强.所用算法适用于介质椭圆柱等目标的散射特性研究,所得结果为目标尺寸、形状遥感电磁测量等实际应用提供了理论依据.