椭球粒子声辐射力计算及分析

针对水下椭球粒子,以声散射理论为基础,采用分波序列的方法,建立了椭球粒子声辐射力的理论计算模型。进而根据声辐射力计算公式,以刚性椭球粒子和液体椭球粒子为例,计算并分析不同Bessel波束作用下椭球粒子的轴向声辐射力函数特征。数值仿真计算结果表明,对于刚性椭球粒子,扁平椭球粒子相对于细长椭球粒子更有助于激发负声辐射力；对于液体椭球粒子,细长椭球粒子相对于扁平椭球粒子更加容易产生负声辐射力；对于不同介质的椭球粒子,不同的入射波束激发的负声辐射力的效果也存在明显的差异。该结果为复杂的尺寸和介质粒子声操控技术提供了理论的可行性。     

自由空间中球形粒子的负向声辐射力

针对一定声场作用下自由空间中的球形粒子,首先分析了声散射过程中的吸收声功率、散射声功率和损失声功率以及三者之间的关系,并通过计算发现了由于参数选取不当导致的负吸收现象。接着从动量守恒定律出发推导了声辐射力的一般表达式,阐释了声辐射力与声能流之间的关系,并从理论和计算两方面验证了负向声辐射力的存在。当负向声辐射力产生时,声波的背向散射被抑制。在此基础上,进一步研究了粒子的偏心特性和流体的黏度这两种常见因素对负向声辐射力的影响。利用球函数的加法公式推导了偏心球的散射系数和声辐射力公式,结果显示偏心距离、粒子的材料等都会显著改变负向声辐射力的产生条件。在低频近似下,由于流体黏度附加的正向声辐射力是否能完全抵消原来的负向声辐射力将决定最终的声辐射力方向。该结果对利用负向声辐射力制成单行波声学镊子来实现对特定粒子的操控有着理论指导意义。      

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

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

近似声斗篷的多层介质覆盖水中弹性球壳的声散射

对水中弹性球壳的低频覆盖层消声效果和消声机理进行研究。声斗篷(Acoustic cloak)是各向异性的超常材料,根据有效介质近似理论(Effective MediumApproximationTheory),采用各向同性多层介质近似各向异性材料的声斗篷,推导出覆盖多层介质水中弹性球壳散射声场的解析表达式,计算了弹性球壳覆盖多层介质前后的散射形态函数、谐振模态以及声场分布,分析了覆盖前后球壳的声散射特性和声传播机理。研究表明,覆盖层内的声波在多层介质之间的声传播方向发生改变,声场形成弯曲变形,声能流绕过目标,这不仅极大的降低了低频散射强度,而且使到达弹性体表面的声场强度非常小,散射频响中除了0阶子波产生的第一个谐振峰外,无法激发弹性球壳的其它谐振模态,有效的抑制弹性球壳的谐振散射。      

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

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

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

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

包覆微米碳颗粒的光散射特性

利用分层介质球光散射理论公式,计算分析了包覆水和黄铜介质层对碳微粒的微分散射特性的影响,比较了包覆前后散射的极化特性。在计算中为保持收敛性,对Bessel函数采用了递推求解,并对求和项数进行了限定。结果表明,在包覆前后,单粒子的微分散射截面发生了变化,而包覆层的影响与其厚度和光在其中的趋肤深度有关;厚度小于趋肤深度的包覆层对粒子的前后向散射的极化特性没有影响,而在其它方向上却影响较大。     

平面行波场中多个粒子受到的声辐射力*

从粒子在行波场中的声散射出发,研究当声场中存在稀疏分布的多个粒子时粒子受到的声辐射力,并且给出了适用于声场中任意位置的粒子声辐射力计算公式。由于声辐射力为非线性力,当声场中存在多个粒子时,直接计算粒子受到的声辐射力非常复杂。结果表明,当声场中存在多个稀疏分布的粒子时,这一多粒子系统可以视为多个单独的双粒子系统的叠加,只需要分别计算各个双粒子系统的声辐射力就可以通过叠加得到声场中任意粒子的声辐射力。这一结果有助于利用声辐射力对微小粒子进行精细操控。     

双球粒子对任意入射单波束及双波束的散射

基于广义Lorenz-Mie理论,研究了双介质球粒子对任意角度入射波束的散射.将入射波束用球矢量波函数展开,推导了双介质球形粒子对任意入射单波束的散射方程;将该方法进一步推广,研究了双波束任意入射时的情况.以Gauss波束为例,对以上理论进行了数值验证,比较分析了单波束及双波束任意入射时散射特性随入射方向、球心距等参数的变化关系. 关键词： 广义波束因子 双球粒子 Gauss波束     

三弹性圆柱壳体多重耦合声辐射计算方法

为研究多圆柱壳组合结构的声辐射特性,采用模态叠加法建立了3个并排无限长弹性圆柱壳的振动声散射耦合物理模型,充分考虑了三圆柱壳的表面振动与散射声场的耦合,其中散射声场可分解为各圆柱壳刚性散射声场和弹性辐射声场的叠加,数学上将各壳间的声场耦合关系通过柱函数加法公式描述。利用该物理模型,分析了多重散射对稳态声场求解结果的影响,比较了三圆柱壳耦合系统与单个圆柱壳系统的辐射声场指向性、声压级及辐射声功率级的差异及其产生机理,结果表明:结构弹性耦合声辐射不仅在低频对总声场有显著影响,在高频范围也不可忽略;另外,针对本文设定参数的组合圆柱壳,在150 Hz以上频段,两旁圆柱壳对中间圆柱壳在正横方位产生了声辐射遮蔽效应,垂直方位则体现声泄漏作用。本文建立的方法可推广到三维空间任意多壳结构的声振耦合建模。      

Axial radiation force of a bessel beam on a sphere and direction reversal of the force

An expression is derived for the radiation force on a sphere placed on the axis of an ideal acoustic Bessel beam propagating in an inviscid fluid. The expression uses the partial-wave coefficients found in the analysis of the scattering when the sphere is placed in a plane wave traveling in the same external fluid. The Bessel beam is characterized by the cone angle beta of its plane wave components where beta=0 gives the limiting case of an ordinary plane wave. Examples are found for fluid spheres where the radiation force reverses in direction so the force is opposite the direction of the beam propagation. Negative axial forces are found to be correlated with conditions giving reduced backscattering by the beam. This condition may also be helpful in the design of acoustic tweezers for biophysical applications. Other potential applications include the manipulation of objects in microgravity. Islands in the (ka, beta) parameter plane having a negative radiation force are calculated for the case of a hexane drop in water. Here k is the wave number and a is the drop radius. Low frequency approximations to the radiation force are noted for rigid, fluid, and elastic solid spheres in an inviscid fluid.     

Negative axial radiation forces on solid spheres and shells in a Bessel beam

Prior computations predict that fluid spheres illuminated by an acoustic Bessel beam can be subjected to a radiation force directed opposite the direction of beam propagation. The prediction of negative acoustic radiation force is extended to the cases of a solid poly(methylmethacrylate) PMMA sphere in water and an empty aluminum spherical shell in water. Compared with the angular scattering patterns for plane wave illumination, the scattering into the back hemisphere is suppressed when the radiation force is negative. This investigation may be helpful in the development of acoustic tweezers and in the development of methods for manipulating objects during space flight.     

Pulling force of acoustic-vortex beams on centered elastic spheres based on the annular transducer model

To solve the difficulty of generating an ideal Bessel beam, an simplified annular transducer model is proposed to study the axial acoustic radiation force(ARF) and the corresponding negative ARF(pulling force) exerted on centered elastic spheres for acoustic-vortex(AV) beams of arbitrary orders. Based on the theory of acoustic scattering, the axial distributions of the velocity potential and the ARF for AV beams of different orders generated by the annular transducers with different physical sizes are simulated. It is proved that the pulling force can be generated by AV beams of arbitrary orders with multiple axial regions. The pulling force is more likely to exert on the sphere with a smaller k_0a(product of the wave number and the radius) for the AV beam with a bigger topological charge due to the strengthened off-axis acoustic scattering. The pulling force decreases with the increase of the axial distance for the sphere with a bigger k_0a.More pulling force areas with wider axial regions can be formed by AV beams using a bigger-sized annular transducer.The theoretical results demonstrate the feasibility of generating the pulling force along the axes of AV beams using the experimentally applicable circular array of planar transducers, and suggest application potentials for multi-position stable object manipulations in biomedical engineering.     

Acoustic radiation force on a sphere in standing and quasi-standing zero-order Bessel beam tweezers

Starting from the exact acoustic scattering from a sphere immersed in an ideal fluid and centered along the propagation axis of a standing or quasi-standing zero-order Bessel beam, explicit partial-wave representations for the radiation force are derived. A standing or a quasi-standing acoustic field is the result of propagating two equal or unequal amplitude zero-order Bessel beams, respectively, along the same axis but in opposite sense. The Bessel beam is characterized by the half-cone angle β of its plane wave components, such that β = 0 represents a plane wave. It is assumed here that the half-cone angle β for each of the counter-propagating acoustic Bessel beams is equal. Fluid, elastic and viscoelastic spheres immersed in water are treated as examples. Results indicate the capability of manipulating spherical targets based on their mechanical and acoustical properties. This condition provides an impetus for further designing acoustic tweezers operating with standing or quasi-standing Bessel acoustic waves. Potential applications include particle manipulation in micro-fluidic lab-on-chips as well as in reduced gravity environments.     

Axial and transverse acoustic radiation forces on a fluid sphere placed arbitrarily in Bessel beam standing wave tweezers

The axial and transverse radiation forces on a fluid sphere placed arbitrarily in the acoustical field of Bessel beams of standing waves are evaluated. The three-dimensional components of the time-averaged force are expressed in terms of the beam-shape coefficients of the incident field and the scattering coefficients of the fluid sphere using a partial-wave expansion (PWE) method. Examples are chosen for which the standing wave field is composed of either a zero-order (non-vortex) Bessel beam, or a first-order Bessel vortex beam. It is shown here, that both transverse and axial forces can push or pull the fluid sphere to an equilibrium position depending on the chosen size parameter ka   (where k$k$ is the wave-number and a$a$ the sphere’s radius). The corresponding results are of particular importance in biophysical applications for the design of lab-on-chip devices operating with Bessel beams standing wave tweezers. Moreover, potential investigations in acoustic levitation and related applications in particle rotation in a vortex beam may benefit from the results of this study.     

Acoustic radiation force of a quasi-Gaussian beam imparted to a solid spherical scatterer in a fluid

The radiation force generated upon the scattering of a quasi-Gaussian acoustic beam on a homo-geneous elastic sphere in a fluid is investigated. It is shown that the force depends nonmonotonically on the ratio between the sphere’s diameter and the beam’s waist. For a given beam power, the radiation force has its maximum value when the diameters are roughly egual to each other. This is due to the resonant excitation of shear waves on the sphere’s surface under the impact of acoustic wave in the surrounding fluid.     

Acoustic radiation force of high-order Bessel beam standing wave tweezers on a rigid sphere

Background and objective

Particle manipulation using the acoustic radiation force of Bessel beams is an active field of research. In a previous investigation, [F.G. Mitri, Acoustic radiation force on a sphere in standing and quasi-standing zero-order Bessel beam tweezers, Annals of Physics 323 (2008) 1604–1620] an expression for the radiation force of a zero-order Bessel beam standing wave experienced by a sphere was derived. The present work extends the analysis of the radiation force to the case of a high-order Bessel beam (HOBB) of positive order m having an angular dependence on the phase ?.

Method

The derivation for the general expression of the force is based on the formulation for the total acoustic scattering field of a HOBB by a sphere [F.G. Mitri, Acoustic scattering of a high-order Bessel beam by an elastic sphere, Annals of Physics 323 (2008) 2840–2850; F.G. Mitri, Equivalence of expressions for the acoustic scattering of a progressive high order Bessel beam by an elastic sphere, IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control 56 (2009) 1100–1103] to derive the general expression for the radiation force function Y_Jm,st(ka,β,m), which is the radiation force per unit characteristic energy density and unit cross-sectional surface. The radiation force function is expressed as a generalized partial wave series involving the half-cone angle β of the wave-number components and the order m of the HOBB.

Results

Numerical results for the radiation force function of a first and a second-order Bessel beam standing wave incident upon a rigid sphere immersed in non-viscous water are computed. The rigid sphere calculations for Y_Jm,st(ka,β,m) show that the force is generally directed to a pressure node when m is a positive even integer number (i.e. Y_Jm,st(ka,β,m)>0), whereas the force is generally directed toward a pressure antinode when m is a positive odd integer number (i.e. Y_Jm,st(ka,β,m)<0).

Conclusion

An expression is derived for the radiation force on a rigid sphere placed along the axis of an ideal non-diffracting HOBB of acoustic standing (or stationary) waves propagating in an ideal fluid. The formulation includes results of a previous work done for a zero-order Bessel beam standing wave (m = 0). The proposed theory is of particular interest essentially due to its inherent value as a canonical problem in particle manipulation using the acoustic radiation force of a HOBB standing wave on a sphere. It may also serve as the benchmark for comparison to other solutions obtained by strictly numerical or asymptotic approaches.     

Acoustic radiation force of a Bessel beam on a porous sphere

The possibility of using acoustic Bessel beams to produce an axial pulling force on porous particles is examined in an exact manner. The mathematical model utilizes the appropriate partial-wave expansion method in spherical coordinates, while Biot's model is used to describe the wave motion within the poroelastic medium. Of particular interest here is to examine the feasibility of using Bessel beams for (a) acoustic manipulation of fine porous particles and (b) suppression of particle resonances. To verify the viability of the technique, the radiation force and scattering form-function are calculated for aluminum and silica foams at various porosities. Inspection of the results has shown that acoustic manipulation of low porosity (<0.3) spheres is similar to that of solid elastic spheres, but this behavior significantly changes at higher porosities. Results have also shown a strong correlation between the backscattered form-function and the regions of negative radiation force. It has also been observed that the high-order resonances of the particle can be effectively suppressed by choosing the beam conical angle such that the acoustic contribution from that particular mode vanishes. This investigation may be helpful in the development of acoustic tweezers for manipulation of micro-porous drug delivery carrier and contrast agents.     

Acoustic radiation force on an air bubble and soft fluid spheres in ideal liquids: Example of a high-order Bessel beam of quasi-standing waves

The partial wave series for the scattering of a high-order Bessel beam (HOBB) of acoustic quasi-standing waves by an air bubble and fluid spheres immersed in water and centered on the axis of the beam is applied to the calculation of the acoustic radiation force. A HOBB refers to a type of beam having an axial amplitude null and an azimuthal phase gradient. Radiation force examples obtained through numerical evaluation of the radiation force function are computed for an air bubble, a hexane, a red blood and mercury fluid spheres in water. The examples were selected to illustrate conditions having progressive, standing and quasi-standing waves with appropriate selection of the waves’ amplitude ratio. An especially noteworthy result is the lack of a specific vibrational mode contribution to the radiation force determined by appropriate selection of the HOBB parameters.