Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 12 May 2022
Sec. Physical Acoustics and Ultrasonics
This article is part of the Research Topic New Trends in Acoustofluidics: Modeling, Experiments, and Applications View all 6 articles

Acoustic Radiation Force and Torque Acting on Asymmetric Objects in Acoustic Bessel Beam of Zeroth Order Within Rayleigh Scattering Limit

  • Centre for Audio, Acoustics and Vibration, University of Technology Sydney, Ultimo, NSW, Australia

Acoustic momentum exchange between objects and the surrounding fluid can be quantified in terms of acoustic radiation force and torque, and depends on several factors including the objects’ geometries. For a one-dimensional plane wave type, the induced torque on the objects with arbitrary shape becomes a function of both, direct polarization and Willis coupling, as a result of shape asymmetry, and has only in-plane components. Here, we investigate, in the Rayleigh scattering limit, the momentum transfer to objects in the non-planar pressure field of an acoustic Bessel beam with axisymmetric wave front. This type of beam is selected since it can be practically realized by an array of transducers that are cylindrically arranged and tilted at the cone angle β which is a proportionality index of the momentum distribution in the transverse and axial propagation directions. The analytical expressions of the radiation force and torque are derived for both symmetric and asymmetric objects. We show the dependence of radiation force and torque on the characteristic parameters β and radial distance from the beam axis. By comparing against the case of a plane travelling plane wave, zero β angle, we demonstrated that the non-planar wavefront of a zeroth order Bessel beam causes an additional radial force and axial torque. We also show that, due to Willis coupling, an asymmetric object experiences greater torques in the θ direction, by minimum of one order of magnitude compared to a plane travelling wave. Further, the components of the partial torques owing to direct polarization and Willis coupling act in the same direction, except for a certain range of cone angle β. Our findings show that a non-planar wavefront, which is quantified by β in the case of a zeroth-order Bessel beam, can be used to control the magnitude and direction of the acoustic radiation force and torque acting on arbitrarily shaped objects, implying that the wavefront should be adjusted according to the object’s shape to impart acoustic momentum in all directions and achieve a desired acoustophoretic response.

1 Introduction

Acoustic manipulation of objects depends on their scattering response, when subjected to an external wave field [15]. Acoustic radiation force and radiation torque originated from the radiated momentum of the scattered waves, with both being proportional to the products of the scattered and of the incident fields [1, 3, 57]. The object-related factors determining these radiation force and torque include the material properties, dimensions, shape asymmetries, internal structure, and absorption capacity [1, 810]. Factors related to the surrounding fluid cover the viscosity, its density and compressibility, boundary conditions and the character of the wave front. Among these factors, the theoretical relation between shape asymmetry and wave front of the incident acoustic field is yet to be investigated.

Acoustic radiation force and torque applied to objects with axisymmetric geometries such as spheres, spheroids, Cassini ovals (red blood cell shape model) and others have been extensively investigated for travelling and standing plane waves, which are 1D propagating waves [1, 2, 5, 10, 1215]. Analytical expressions of the force and torque were derived for such geometries using the partial-wave expansion method [79, 1622]. The same analytical approach was used to study the radiation force and torque acting on axisymmetric objects due to acoustic Bessel beams of zeroth and first orders [11, 2331]. Acoustic Bessel beams are of the 2D type of wave propagation, with an axisymmetric wave front that is characterised by an axis and a cone angle β, which indicates how much of wave propagation is along and normal to the beam axis, as shown in Figure 1A. Compared to the travelling plane waves, spherical objects of certain sizes can experience a pull radiation force in the opposite direction of the wave propagation, when positioned on the beam axis [30, 32]. For off-axis positions, the acoustic radiation force acts in the axial and transverse directions of the beam and there is an additional torque applied to the objects [25, 28]. These effects are examples of demonstrating the potential of controlling the acoustic manipulation of objects by changing the incident wave front.

FIGURE 1
www.frontiersin.org

FIGURE 1. Acoustic Bessel beam of zeroth order in cylindrical coordinate system with beam axis aligned with the z-axis, (A) the wave front variation over a half wavelength in the z-direction, showing a travelling wave propagation, (B) the circular arrangement of sources at infinite distance emitting plane waves at angle β with respect to the z-axis to generate a Bessel beam of zeroth order [11]. The cylindrical coordinates (R, θ, z) of the position vector of an object with arbitrary shape located off-axis and the local Cartesian coordinate system are shown. The wave vector k and its components kR and kz and their relation with the cone angle β are shown in panel (A). The dotted lines in panel (A) shows the pressure nodes (P.N.) of the Bessel beam, implying a standing wave behaviour in the radial R-direction.

Asymmetry in the shape of an object changes the radiation force and torque applied to them by inducing the Willis coupling between the monopole-dipole scattering response and the incident velocity and pressure fields, respectively [3335]. For objects within the Rayleigh scattering limit, the shape asymmetry manifests as non-zero coefficients of the polarizability tensor of up to dipole accuracy [5]. Analytical expressions of the radiation force and torque due to the shape asymmetry show the contribution from Willis coupling effects, compared to direct polarizability coefficients that are the diagonal sub-tensors of the polarizability tensor [5]. Since the acoustic polarization tensor is independent of the excitation type, the radiation force and torque of an object with an arbitrary shape in an acoustic Bessel beam of zeroth order can be investigated directly using the generic force and torque equations in Ref. [5].

In this paper, we investigate analytically the effects of a Bessel beam of zeroth order, as an axisymmetric case of non-planar travelling wave, on the acoustic radiation force and torque exerted on a object with arbitrary shape. This type of incident wave has a pressure anti-node on the axis, which can be considered as an approximate model for focused acoustic beams. It is assumed that there is no phase variation across the wavefront; however, the same analytical approach can be used to obtain the radiation force and torque for a vortex-type non-planar beam [36, 37]. The dependence of the radiation force and torque on the location of an object with respect to the beam axis are studied to provide insight about the acoustophoretic response of objects located off-axis. The objects are considered to be much smaller than the wavelength within the Rayleigh limit, and their polarizability tensor up to monopole-dipole approximation is given. Comparing the radiation force and torque induced by Bessel beam against those of the plane travelling wave provides insights about acoustic manipulation of objects with arbitrary shapes using non-planar engineered beams with potential implications for various fields of science and engineering, including the ultrasonic manipulation of biological structures.

2 Theory

The acoustic wave equations derived from the first-order approximation of the Navier-Stokes equations and in the zero viscosity limit are expressed [1, 2], as follows,

ttp=cf22p,p=ρftv,p=cf2ρ(1)

where p, ρ and v denote the acoustic pressure, density and velocity fields, respectively, cf and ρf are the speed of sound and the mean density of the fluid, respectively, denotes the spatial gradient operator, and t denotes the differentiation with respect to time t. The acoustic fields are time harmonic, e.g., p(x, t) = p(x)ejωt with x being the position vector, ω = 2πf and f denoting the wave frequency. We assume that the object size is within the Rayleigh limit, which is expressed as ka ≪ 1 with k = ω/cf being the wavenumber and a being the characteristic length of the object.1 Then, the scattered pressure field of such small object can be approximated by monopole-dipole partial fields [5, 34, 38], as follows,

psω2ΩMG+ω2ΩDG,G=ejkr4πrejωt,(2)

where M and D are the volumetric monopole and dipole moments, Ω is the volume of the object, and G is the free-space Green’s function of the acoustic wave equation, Eq. 1. By employing the acoustic polarizability, the scattering moments M and D can be expressed in terms of the incident pressure and velocity fields [5], as follows,

MD=α̂pivi,α̂=1Ωα̂ppα̂pvTα̂vpα̂vv,α̂pv=αpvxαpvyαpvz,α̂vp=αvpxαvpyαvpz,α̂vv=αvvxxαvvxyαvvxzαvvyxαvvyyαvvyzαvvzxαvvzyαvvzz,(3)

where α̂ denotes the polarizability tensor, α̂pv and α̂vp are the Willis coupling coefficients, associated with the object’s shape asymmetry, while α̂pp and α̂vv are the direct polarization coefficients. The polarizability tensor is independent of the incident fields, neglecting non-linear scattering response or fluid properties, and describes the object’s scattering response as a function of its geometry, material properties, absorption and dissipation mechanisms within the fluid. A Cartesian description of these sub-tensors with respect to the coordinate system attached to the object, denoted by the hat ( ˆ ), is provided in Eq. 3. Considering the tensor properties under translation and rotation [22], the description of α in cylindrical coordinates, shown in Figure 1B, becomes

α=α̂ppα̂pvTRTRα̂vpRα̂vvRT,R=exeReyeR0exeθeyeθ0001,(4)

where R is the rotation tensor from the local Cartesian coordinates to the global cylindrical coordinates, and e denotes the unit basis vector in the direction indicated by its subscript. For the case of a sphere, the Willis coupling coefficients are zero and α̂vv=α̂vvI, with I denoting the identity tensor. For axisymmetric objects with a plane symmetry in the axis direction, e.g., prolate and oblate spheroids, Cassini ovals etc. [10], the Willis coupling coefficients become zero too, αvvxx=αvvxxαvvzz, assuming the axis of symmetry is in the z-direction, and αvvxy=αvvyx=αvvxz=αvvzx=αvvyz=αvvzy=0. In this study, we consider an object with arbitrary shape such that all the polarizability coefficients in Eq. 3 are non-zero, to provide a general formulation of acoustic radiation force and torque due to Bessel beams. Although we assume that the acoustic polarizability tensor is given, it can be obtained for objects with arbitrary shapes computationally using Boundary Element method [5, 39, 40], or Finite Element method [33, 41]. Here, we consider the zeroth order Bessel beam to derive the analytical expressions; however, the same formalism can be applied to Bessel beams of higher orders, or other non-planar beams regardless of the phase variation across the wavefront.

2.1 Acoustic Bessel Beam of Zeroth Order

The pressure field of an acoustic Bessel beam of zeroth order is expressed [11], as follows,

pi=AJ0ejkzzejωt,J0=J0kRR,R=x2+y2(5)

where J0 denotes the regular Bessel function of zeroth order, and A is the magnitude of the wave. According to the representation of a Bessel beam using Durnin rings [42, 43],

J0kRR=12π02πejkRxcosφ+kRycosφdφ,(6)

this type of incident beam can be generated from the superposition of plane travelling waves coming towards a given axis at an incidence angle of β, as illustrated in Figure 1B. The Bessel beam is of the travelling type along the beam axis, i.e., z-direction, while fixed pressure and velocity nodes are (standing wave profile) along the R-direction. This special acoustic beam is axisymmetric, which implies the same acoustic pressure for points that are at the same radius from the axis and θp = 0. The pressure nodes are indicated in Figure 1A, and a zone in the vicinity of the beam axis is shaded as a region of interest for investigating the radiation force and torque exerted off-axis objects. The velocity field of this Bessel beam and its spatial derivatives are as in the following,

piez=jAkcosβJ0ejkzzejωt,pieR=AksinβJ1ejkzzejωt,viez=AρfcfcosβJ0ejkzzejωt,vieR=jAρfcfsinβJ1ejkzzejωt,vi:ezez=jAkρfcfcos2βJ0ejkzzejωt,vi:eReR=jAkρfcfsin2βJ0J1kRRejkzzejωt,vi:eRez=vi:ezeR=Ak2ρfcfsin2βJ1ejkzzejωt,(7)

where : operator denotes inner product of two second-order tensors, =ReR+1Rθeθ+zez, and the dyadic products of the unit basis vectors are denoted by ezeR, ezez, and eReR. From Eq. 7, it can be seen that all the derivative fields are also axisymmetric, i.e., vi: eReθ = vi: eθez = vi: eθeθ = vieθ = pieθ = 0. These field expressions of the Bessel beam of the zeroth order in Eqs 5, 7 are required to formulate the dependence of the radiation force and torque on cylindrical coordinates (R, θ, z) and cone angle β.

2.2 Acoustic Radiation Force and Torque

Acoustic radiation force and torque are exerted on an object within an incident field as a result of radiation stresses ⟨σ⟩ that can mathematically be expressed as follows,

σ=12κfp212ρfv2Iρfvv,(8)

where κf=1/ρfcf2 is the mean fluid compressibility, ⟨⟩ operator denotes a time-averaging over one wave period [2]. The first term on the right hand side of Eq. 8 is called radiation pressure that is generated by the radiated momentum of the scattered pressure field, and the second term represents the Reynolds stresses indicating the contribution from the object’s surface oscillations to the radiation momentum [1, 2, 5, 8, 9, 44]. Acoustic radiation force and torque are obtained from radiation stresses ⟨σ⟩, as follows,

F=Γσ·ndΓ,T=Γx×σ·ndΓ,(9)

where Γ denotes a surface enclosing the object and n is the outward normal vector of Γ surface. By using the far-field approach and integrating the stresses on any fictitious surface enclosing the object [2, 8, 45], the radiation force and torque exerted on an object with arbitrary shape and size in the Rayleigh scattering limit ka ≪ 1 is expressed [5], as follows,

F=Fd+Fc,Fd=αpp2ρfpi2+jωαvvvivir=0,Fc=1ρfαpvpivivipir=0.T=Td+Tc,Td=jωαvvvi×vir=0,Tc=jωpiαvp×vir=0,(10)

where d and c subscripts indicate the force and torque associated with direct polarization and Willis coupling, respectively. By substituting Eqs 5, 7 into Eq. 10 and using the reciprocity properties of the polarizability tensor [5, 22, 35, 38, 46], the analytical expressions of the force and torque are obtained, as in the following,

Fdez=2Eikρf1κfImαppcosβJ02+kcfReαvvzzcos3βJ02+kcfReαvvRRsin2βcosβJ12,FdeR=2Eikρf1κfReαppsinβJ0J1+kcfReαvvRzsin2βcosβJ02+J12J0J1kRR+kcfImαvvzzcos2βsinβJ0J1kcfImαvvRRsin3βJ0J1J12kRR,(11)
Fcez=2Eikρf2cfReαpvRsinβcosβJ0J1,FceR=2EikρfcfImαpvRsin2βJ02+J12J0J1kRR,(12)
Tdez=2EikρfcfReαvvθzsinβcosβJ0J1+cfImαvvθRsin2βJ12,Tdeθ=2EikρfcfReαvvzzsinβcosβJ0J1cfImαvvzRsin2βJ12+cfImαvvRzcos2βJ02+cfReαvvRRsinβcosβJ0J1,TdeR=2EikρfcfImαvvθzcos2βJ02cfReαvvθRsinβcosβJ0J1,(13)
Tcez=2Eikρf1κfReαvpθsinβJ0J1,Tceθ=2Eikρf1κfReαvpzsinβJ0J1+1κfImαvpRcosβJ02,TceR=2Eikρf1κfImαvpθcosβJ02,(14)

where Ei=A2/4ρfcf2 denotes the energy density of the incident Bessel beam. Both partial forces Fd and Fc have components only in z- and R-directions, which implies objects being pushed or pulled in axial and radial directions in an axisymmetric acoustic field. For the partial torques Td and Tc, in addition to the component in the θ-direction, there exist two more components the z- and R-plane, which means a full 3D torque is induced by an axisymmetric acoustic field. This is in contrast to the case of a standing, plane wave, where only in-plane torque components are generated [5, 22]. Compared to a plane wave, this torque property of the Bessel beam shows that a 3D rotational manipulation an object with arbitrary shape can be realized even using 2D propagation type, axisymmetric beams. The real and imaginary parts of all polarizability coefficients, from the cylindrical description of the polarizability tensor in Eq. 4, are required to obtain the three components of the radiation torque; while only the z- and R-dependent coefficients are required to calculate both partial forces. This indicates the significance of accounting for the shape of the object to find the radiation force and torque in a Bessel beam of zeroth order.

3 Results

To further investigate the changes of Bessel radiation force and torque, we examine the expressions in Eqs 1114 in four limit cases of (I) β = 0, i.e., plane travelling wave, (II) β tends to 0 corresponding to radially weak Bessel beam, (III) R = 0 indicating on-axis location of the object’s centroid, and (IV) kRR ≪ 1, which is the vicinity of the axis, shown by the yellow-shaded strip in Figure 1A. For general off-axis locations, the force and torque should be obtained directly from Eqs 1114.

3.1 Case I: β = 0 - Travelling Plane Wave

This case corresponds to zero radial variation of the wavefront kR = 0, which means the Bessel beam degenerates to a travelling plane wave. By substituting sin β = 0, cos β = 1, J0 = 1, J1 = 0, and J1/kRR ≈ 1/2, the partial radiation forces and torques in Eqs 1114 become

Fdez=2Eikρf1κfImαpp+kcfReαvvzz,FdeR=0,(15)
Fcez=0,FceR=0,(16)
Tdez=0Tdeθ=2EikρfcfImαvvRz,TdeR=2EikρfcfImαvvθz,(17)
Tcez=0,Tceθ=2Eikρf1κfImαvpRTceR=2Eikρf1κfImαvpθ.(18)
Equations 1518 are identical to expressions in Eq. [26] of Ref. [5], verifying our formulation of the radiation force and torque under a Bessel beam at the limit case of β = 0. It is also found that Willis coupling effects show zero contribution to the force and only generates torque in the in-plane directions of a travelling plane wave. The forces and torques generated by a travelling plane wave are constant and independent of the location of the objects, meaning that the radiation force and torque fields are uniform and only depend on the polarizability response of the object.

3.2 Case II: β → 0 - Radially Weak Bessel Beam

This case corresponds to the long wavelength limit in the radial R direction, and kR = k sin β. This case is of interest since it represents relatively weaker propagation of acoustic energy in the R direction, i.e., the standing wave profile and fixed pressure and velocity nodes. By using the asymptotic approximation of the Bessel function in this limit, i.e., J0 ≈ 1 and J1kRR/2, sin ββ, cos β ≈ 1 and β2β, Eqs 1114 can be approximated, as follows,

Fdez=2Eikρf1κfImαpp+kcfReαvvzz,FdeR=2Eikρfk2κfReαppβ2R+k2cf2Imαvvzzβ2R+kcf2ReαvvRzβ2,(19)
Fcez=2EikρfkcfReαpvRβ2R,FceR=2Eikρfcf2ImαpvRβ2,(20)
Tdez=2Eikρfkcf2Reαvvθzβ2R,Tdeθ=2Eikρfkcf2Reαvvzz+αvvRRβ2R+cfImαvvRz,TdeR=2EikρfcfImαvvθzkcf2ReαvvθRβ2R,(21)
Tcez=2Eikρfk2κfReαvpθβ2R,Tceθ=2Eikρfk2κfReαvpzβ2R+ImαvpR,TceR=2Eikρf1κfImαvpθ.(22)

For a given cone angle β, we can see from Eqs 1922 that the radiation force and torque spatially depend only on their radial coordinate R, changing linearly, while these fields are independent of the axial coordinate z, due to the travelling wave nature of propagation in this direction. Compared to the limit case of β = 0, it was found that the additional terms proportional to β2 appear in the force and torque expressions TceR, as result of the weak propagation in the R-direction, except for the radial component of the Willis coupling torque TceR that remains unchanged. These analytical expressions reveal the sensitivity of the radiation force and torque, both direct and Willis coupling parts, to the change in the wave front that is the weak Bessel-type undulations. This provides β angle as another degree of freedom for manipulation of objects with arbitrary shape, using this type of radially weak Bessel beams.

3.3 Case III: On-Axis Object R = 0

There is a velocity node, corresponding to maximum pressure, on the axis of the zeroth-order Bessel beam. This is of interest since it is a relatively good example of focused acoustic beams that can be produced using a meta-lens [47] or a phased transducer array [48]. When the centroid of an object is located on the beam axis, the expressions of acoustic radiation force and torque in Eqs 1114 are further simplified, as follows,

Fdez=2Eikρf1κfImαppcosβ+kcfReαvvzzcos3β,FdeR=2Eikρfkcf2ReαvvRzsin2βcosβ,(23)
Fcez=0FceR=2Eikρfcf2ImαpvRsin2β,(24)
Tdez=0Tdeθ=2EikρfcfImαvvRzcos2β,TdeR=2EikρfcfImαvvθzcos2β,(25)
Tcez=0Tceθ=2Eikρf1κfImαvpRcosβ,TceR=2Eikρf1κfImαvpθcosβ.(26)

From Eq. 23, we found that an object on the beam axis is subjected to only an axial force, i.e., in the z-direction as a result of the direct polarization.2 Since the range of values of the cone angle is 0 ≤ β < π/2, cos β > 0 and the direction of the axial force is determined from the polarization coefficients, i.e., Im(αpp) and Re(αvvzz), implying that some objects can experience pull-in effects from a Bessel beam of zeroth order at the ka ≪ 1 limit. The axial components of the Willis-coupling partial force Fc and both direct and Willis coupling torques Td and Tc are zero. Equations 25, 26 show that the dependence of the direct and Willis coupling-induced torques on cone angle β are cos2β and cos β, respectively, implying that partial torque fields Td and Tc are independent of each other and the shape asymmetry needs to be accounted for separately. Similarly, the partial force fields Fd and Fc are independent, as can be seen from Eqs 23, 24. Finally, the most significant contribution of the Willis coupling is the radial component of the radiation force, coming from FceR, meaning that objects with shape asymmetry are pushed away from the beam axis in the direction that Im(αpvR) is non-zero and maximum, if the radiation torque is accounted for. This depends on the orientation of the objects and its shape asymmetry with respect to the beam axis. For axisymmetric shapes with one degree of asymmetry along the axis of revolution, e.g., shapes of a Helmholtz resonator, unhealthy red blood cells [49], or a meta-atom with controlled Willis coupling [22], this radial force is zero if the axis of symmetry of the object is aligned with the beam axis. It is inferred that, in general, objects tends to be off-axis in such acoustic beams; hence, the axisymmetric wavefront provides the same level of control compared to a 2D arbitrary wavefront.

3.4 Case IV: Object in the Vicinity of the Beam Axis kRR ≪ 1

We consider this special off-axis case such that kRR ≪ 1, as shown by the yellow-shaded strip in Figure 1A. This condition applies to long wavelength limit ka ≪ 1, and we assume the radial distance is close to zero (vicinity of beam axis), while 0 ≤ β < π/2. By using the asymptotic approximation of Bessel functions J0 and J1, the radiation force and torque in Eqs 1114 are approximated, as follows,

Fdez=2Eikρf1κfImαppcosβ+kcfReαvvzzcos3β+k3cf4ReαvvRRsin4βcosβR2,FdeR=2Eikρfk2κfReαppsin2βR+k2cf2Imαvvzzcos2βsin2βR+kcf2ReαvvRzsin2βcosβ1+k22sin2βR2kcf4ImαvvRRsin4βR,(27)
Fcez=2EikρfkcfReαpvRsin2βcosβR,FceR=2EikρfcfImαpvRsin2β12+k2sin2β4R2,Tdez=2Eikρfkcf2Reαvvθzsin2βcosβR+k2cf4ImαvvθRsin4βR2,(28)
Tdeθ=2Eikρfkcf2Reαvvzzsin2βcosβRk2cf4ImαvvzRsin4βR2+cfImαvvRzcos2β+kcf2ReαvvRRsin2βcosβR,TdeR=2EikρfcfImαvvθzcos2βkcf2ReαvvθRsin2βcosβR,(29)
Tcez=2Eikρfk2κfReαvpθsin2βR,Tceθ=2Eikρfk2κfReαvpzsin2βR+ImαvpRcosβ,TceR=2Eikρf1κfImαvpθcosβ.(30)

From Eq. 30, it was found that the TceR is constant and independent of the R-coordinate. Other components of the partial radiation force and torques change with the radial coordinate R, e.g., proportional to R and R2. There are also some constant terms in Fdez, FceR, Tdeθ, TdeR, Tceθ, and TceR, which provide the base value of these force and torque fields in the kRR ≪ 1 region irrespective of the object’s radial location.

3.5 Numerical Analysis

To investigate the influence of non-planar wavefront of the zeroth order Bessel beam on the acoustic radiation force and torque, we consider an object with asymmetric shape such that major asymmetry occurs along the z-direction for θ = 0 and the polarizability coefficients become

αppΩs=7.2×103+0.0j,αvp=jωρfαpvαpvαpp=8.4×104+7.2×104j3.4×10104.8×105j0.0,,7.4×102j,αvvαpp=1.0,×,1042.39j3.6×1092.89,×,105j0.0,+,1.3×103j1.0,×,1042.39j5.0,×,1082.4×104jSym.1.0,×,1042.43j.(31)

where Ωs denote the volume of the object, and these values satisfy the acoustic energy balance, i.e., they are less than the maximum admissible polarizability [38]. Although these are just examples of the values of polarizability coefficients, they correspond to a sphere of size a with a blind circular hole which leads to the asymmetry along the hole axis [5], i.e., z-direction in this study. The non-zero coefficients with relatively smaller values correspond to the numerical discretization of this object for the calculation of monopole and dipole moments using Boundary Element Method, as outlined in Ref. 5. Nevertheless, these small values indicate the level of asymmetry in x- and y-directions, compared to the intended one in the z-direction.

The ratio of partial forces and torques in cases II to IV against the force and torque of a plane travelling wave, i.e., case I, is considered for comparison. The magnitude of the force and torque for case I are expressed, as follows,

Q=Fd+FcFd+Fc,Z=Td+TcTd+Tc.(32)

First, we investigate the case of weak Bessel beam, i.e., case II and Eqs 1922, by choosing β = 5°, corresponding to kR/k ≈ 0.1. The three components of the partial forces and torques are presented in Figure 2, for relatively large range of 0 < kRR < π.

FIGURE 2
www.frontiersin.org

FIGURE 2. Partial forces, panels (A),(B), and torques, panels (C),(D), due to direct and Willis coupling acoustic polarization for a weak Bessel beam of β = 5° and kR/K ≈ 0.1, with respect to the radial distance from the beam axis at R = 0 for the asymmetric object given by its polarizability coefficients in Eq. 31. The components of the forces and torques in the cylindrical coordinate system (eR, eθ, ez) are shown by different line types. For the forces in panel (A),(B), the θ-component is zero. The magnitude of these forces are shown in the logarithmic scale to indicate the difference in orders of magnitude.

The direct partial force Fd shows a much larger radial component than the axial z-component, which is almost the same as the axial force component of a plane travelling wave, as shown in Figure 2A. Similar differences between components are observed in the Willis-coupling partial force Fc, as shown in Figure 2B; however, the contribution of Fc to the total force is negligible compared to Fd. It is noted that the only non-zero force component for the reference case of a plane travelling wave is the axial component of the direct partial force Fdez. In all cases I to IV, the force component in the θ-direction is always zero due to the axisymmetric wavefront of the Bessel beam. Comparing the partial torques in Figures 2C,D, we observed that the component in the θ-direction is the largest and the Willis coupling contribution is much larger than the direct partial torque. The R-component of the partial torques is of the same of order of magnitude, but the z- and θ-components of Willis coupling torque Tc are much larger than those of Td. The force and torque components changes almost linearly with respect to the radial R-coordinate, except for R → 0 corresponding to region in the vicinity of the beam axis, i.e., R = 0. These results indicate that the effects of asymmetry of the object’s shape manifest as a relatively large torque component in the θ-direction, for Bessel beam of small β angle representing a weak non-planar travelling wave.

Next, the case of objects at the vicinity of the beam axis, i.e., case IV, is shown in Figure 3. The same object with polarizability tensor given in Eq. 31 is placed one radius a away from the beam axis, i.e., R/a = 1 corresponding to kR ≈ 0.03. The force and torque results are shown for cone angle β in the range of 0–90°. As shown in Figure 3A, the direct partial force Fd has a larger component in the radial R-direction. The component in the z-direction undergoes a sign reversal at β ≈ 78°, from opposite (pull) to same (push) as the wave propagation direction along the z-axis. The Willis coupling partial force Fc, shown in Figure 3B, is smaller than the Fd by several orders of magnitude and can be neglected. Comparing the partial torques Td and Tc, shown in Figures 3C,D, respectively, it is observed that the asymmetry in shape results in relatively larger torque components in the θ-direction, and the magnitude increases as the cone angle β approaches π/2, corresponding to the limit of zero axial propagation. Finally, the significant difference between the magnitudes of the force and torque components correspond to the values of the polarizability coefficients for this study, which comes from a discretized geometry of a sphere with a blind circular hole. Nevertheless, they can be considered as an indication of the force and torque acting upon an object with asymmetric shape from several directions. Finally, our results show that the cone angle β can be used as a design parameter for axial manipulation using zeroth order Bessel beam.

FIGURE 3
www.frontiersin.org

FIGURE 3. Partial forces, panels (A),(B), and torques, panels (C),(D), due to direct and Willis coupling acoustic polarization in the vicinity of the beam axis at R/a = 1 with a being the nominal radius of the asymmetric object given by its polarizability coefficients in Eq. 31. The results are shown for the entire range of 0 < β < 90° to demonstrate the effects of wave front deviating further from a plane (β = 0) on an asymmetric object of small size ka ≈ 0.03. The components of the forces and torques in the cylindrical coordinate system (eR, eθ, ez) are shown by different line types. For the forces in panel (A),(B), the θ-component is zero. The magnitude of these forces are shown in the logarithmic scale to indicate the difference of orders of magnitude.

4 Conclusion

Acoustic radiation force and torque on objects with arbitrary shape were obtained by applying an acoustic Bessel beam of zeroth order, which has a pressure anti-node on the beam axis and fixed pressure nodes along the radial direction, as an estimate model of engineered non-planar beams using acoustic focusing techniques. To investigate the effects of this non-planar wavefront, analytical expressions of the radiation force and torque were derived by using the far-field approach, employing the polarizability concept, by remaining under the assumption of having a lossless fluid for sub-wavelength objects in the Rayleigh scattering limit. These expressions were verified by comparing against those for plane travelling waves in the limit case of zero cone angle β, which is a measure of energy propagation in the axial and radial directions. For cases of weak radial field, i.e., β tends to zero, it was found that additional terms proportional to β2 emerges, compare to the case of a plane travelling wave, which implies the contribution from the weak radial field. We also showed that an asymmetric object located initially on the beam axis is pushed away from the axis due to the contribution of the Willis coupling partial force. Furthermore, an object located on the axis, where pressure is maximum, can experience a pull-in effect for certain values of direct polarizability coefficients and cone angle β. Finally, for off-axis cases in the vicinity of the beam axis, we found that the additional terms are proportional to R and R2, compared to the case of plane travelling wave, showing the dependence on the radial position in a Bessel beam. These findings are of interest for applications of acoustic tweezing or levitation of thin elastic structures and biological samples with asymmetric geometries, using non-planar acoustic beams that can be produced by novel ultrasound meta-materials, in Space engineering and Bio-material engineering.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author Contributions

SS contributed to conception and design of the study, conducted the theoretical analysis, organized the database, and wrote the first draft of the manuscript. SO contributed to the conception, provided resources and administered the project. Both read and revised the manuscript, and approved the submitted version.

Funding

This research has been financially supported by Australian Research Council Discovery Projects DP200101708 and DP200100358 and by the UTS Faculty of Engineering and IT as well as UTS Techlab.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors acknowledge preliminary discussions with David Powell and Yan Kei Chiang.

Footnotes

1the dimension estimated by, e.g., dividing the object’s volume through its surface area for 3D objects.

2a combination of Im(αpp) cos β and Re(αvvzz)cos3β.

References

1. Doinikov AA. Acoustic Radiation Pressure on a Compressible Sphere in a Viscous Fluid. J. Fluid Mech. (1994) 267:1–22. doi:10.1017/s0022112094001096

CrossRef Full Text | Google Scholar

2. Bruus H. Acoustofluidics 7: The Acoustic Radiation Force on Small Particles. Lab Chip (2012) 12:1014–21. doi:10.1039/c2lc21068a

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Silva GT, Bruus H. Acoustic Interaction Forces between Small Particles in an Ideal Fluid. Phys Rev E Stat Nonlin Soft Matter Phys (2014) 90:063007. doi:10.1103/PhysRevE.90.063007

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lima EB, Leão-Neto JP, Marques AS, Silva GC, Lopes JH, Silva GT. Nonlinear Interaction of Acoustic Waves with a Spheroidal Particle: Radiation Force and Torque Effects. Phys Rev Appl (2020) 13:064048. doi:10.1103/physrevapplied.13.064048

CrossRef Full Text | Google Scholar

5. Sepehrirahnama S, Oberst S, Chiang YK, Powell D. Acoustic Radiation Force and Radiation Torque beyond Particles: Effects of Nonspherical Shape and Willis Coupling. Phys. Rev. E (2021) 104:065003. doi:10.1103/PhysRevE.104.065003

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Settnes M, Bruus H. Forces Acting on a Small Particle in an Acoustical Field in a Viscous Fluid. Phys Rev E Stat Nonlin Soft Matter Phys (2012) 85:016327. doi:10.1103/PhysRevE.85.016327

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sepehrirahnama S, Lim K-M, Chau FS. Numerical Analysis of the Acoustic Radiation Force and Acoustic Streaming Around a Sphere in an Acoustic Standing Wave. Phys Procedia (2015) 70:80–4. doi:10.1016/j.phpro.2015.08.047

CrossRef Full Text | Google Scholar

8. Doinikov AA. Acoustic Radiation Pressure on a Rigid Sphere in a Viscous Fluid. Proc R Soc Lond. Ser A Math Phys Sci (1994) 447:447–66.

Google Scholar

9. Sepehrirhnama S, Lim K-M. Generalized Potential Theory for Close-Range Acoustic Interactions in the Rayleigh Limit. Phys. Rev. E (2020) 102:043307. doi:10.1103/PhysRevE.102.043307

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lima EB, Silva GT. Mean Acoustic Fields Exerted on a Subwavelength Axisymmetric Particle. J Acoust Soc Am (2021) 150:376–84. doi:10.1121/10.0005625

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Marston PL. Axial Radiation Force of a Bessel Beam on a Sphere and Direction Reversal of the Force. J Acoust Soc Am (2006) 120:3518–24. doi:10.1121/1.2361185

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Dual J, Hahn P, Leibacher I, Möller D, Schwarz T, Wang J. Acoustofluidics 19: Ultrasonic Microrobotics in Cavities: Devices and Numerical Simulation. Lab Chip (2012) 12:4010–21. doi:10.1039/c2lc40733g

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Leão-Neto JP, Hoyos M, Aider J-L, Silva GT. Acoustic Radiation Force and Torque on Spheroidal Particles in an Ideal Cylindrical Chamber. J Acoust Soc Am (2021) 149:285–95. doi:10.1121/10.0003046

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Foresti D, Nabavi M, Poulikakos D. On the Acoustic Levitation Stability Behaviour of Spherical and Ellipsoidal Particles. J. Fluid Mech. (2012) 709:581–92. doi:10.1017/jfm.2012.350

CrossRef Full Text | Google Scholar

15. Jerome TS, Ilinskii YA, Zabolotskaya EA, Hamilton MF. Acoustic Radiation Force on a Compressible Spheroid. J Acoust Soc Am (2020) 148:2403–15. doi:10.1121/10.0002277

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hasegawa T, Yosioka K. Acoustic‐Radiation Force on a Solid Elastic Sphere. J Acoust Soc Am (1969) 46:1139–43. doi:10.1121/1.1911832

CrossRef Full Text | Google Scholar

17. Lopes JH, Azarpeyvand M, Silva GT. Acoustic Interaction Forces and Torques Acting on Suspended Spheres in an Ideal Fluid. IEEE Trans. Ultrason, Ferroelect., Freq. Contr. (2016) 63:186–97. doi:10.1109/tuffc.2015.2494693

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Sepehrirahnama S, Lim K-M, Chau FS. Numerical Study of Interparticle Radiation Force Acting on Rigid Spheres in a Standing Wave. J Acoust Soc Am (2015) 137:2614–22. doi:10.1121/1.4916968

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sepehrirahnama S, Chau FS, Lim KM. Numerical Calculation of Acoustic Radiation Forces Acting on a Sphere in a Viscous Fluid. Phys Rev E Stat Nonlin Soft Matter Phys (2015) 92:063309. doi:10.1103/PhysRevE.92.063309

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Gong Z, Baudoin M. Acoustic Radiation Torque on a Particle in a Fluid: An Angular Spectrum Based Compact Expression. J Acoust Soc Am (2020) 148:3131–40. doi:10.1121/10.0002491

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gong Z, Baudoin M. Equivalence between Angular Spectrum-Based and Multipole Expansion-Based Formulas of the Acoustic Radiation Force and Torque. J Acoust Soc Am (2021) 149:3469–82. doi:10.1121/10.0005004

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Sepehrirahnama S, Oberst S, Chiang Y, Powell D. Willis Coupling-Induced Acoustic Radiation Force and Torque Reversal (2021). 01354. arXiv preprint (2021) arXiv:2110. doi:10.48550/arXiv.2110.01354

CrossRef Full Text | Google Scholar

23. Marston PL. Axial Radiation Force of a Bessel Beam on a Sphere, Direction Reversal of the Force, and Solid Sphere Examples. J Acoust Soc Am (2007) 121:3109. doi:10.1121/1.4808510

CrossRef Full Text | Google Scholar

24. Mitri FG. Acoustic Radiation Force of High-Order Bessel Beam Standing Wave Tweezers on a Rigid Sphere. Ultrasonics (2009) 49:794–8. doi:10.1016/j.ultras.2009.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Silva GT. Off-axis Scattering of an Ultrasound Bessel Beam by a Sphere. IEEE Trans. Ultrason, Ferroelect., Freq. Contr. (2011) 58:298–304. doi:10.1109/tuffc.2011.1807

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Silva GT, Lobo TP, Mitri FG. Radiation Torque Produced by an Arbitrary Acoustic Wave. Epl (2012) 97:54003. doi:10.1209/0295-5075/97/54003

CrossRef Full Text | Google Scholar

27. Wijaya FB, Lim K-M. Numerical Calculation of Acoustic Radiation Force and Torque Acting on Rigid Non-spherical Particles. Acta Acustica united Acustica (2015) 101:531–42. doi:10.3813/aaa.918850

CrossRef Full Text | Google Scholar

28. Wijaya FB, Sepehrirahnama S, Lim K-M. Interparticle Force and Torque on Rigid Spheroidal Particles in Acoustophoresis. Wave Motion (2018) 81:28–45. doi:10.1016/j.wavemoti.2018.06.004

CrossRef Full Text | Google Scholar

29. Sepehrirahnama S, Lim KM. Acoustophoretic Agglomeration Patterns of Particulate Phase in a Host Fluid. Microfluid Nanofluidics (2020) 24:1–14. doi:10.1007/s10404-020-02397-5

CrossRef Full Text | Google Scholar

30. Zhang L, Marston PL. Geometrical Interpretation of Negative Radiation Forces of Acoustical Bessel Beams on Spheres. Phys. Rev. E (2011) 84:035601. doi:10.1103/PhysRevE.84.035601

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Fan X-D, Zhang L. Trapping Force of Acoustical Bessel Beams on a Sphere and Stable Tractor Beams. Phys. Rev. Appl (2019) 11:014055. doi:10.1103/PhysRevApplied.11.014055

CrossRef Full Text | Google Scholar

32. Marston PL. Negative Axial Radiation Forces on Solid Spheres and Shells in a Bessel Beam. J Acoust Soc Am (2007) 122:3162–5. doi:10.1121/1.2799501

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jordaan J, Punzet S, Melnikov A, Sanches A, Oberst S, Marburg S. Measuring Monopole and Dipole Polarizability of Acoustic Meta-Atoms. Appl Phys Lett (2018) 113:224102. doi:10.1063/1.5052661

CrossRef Full Text | Google Scholar

34. Su X, Norris AN. Retrieval Method for the Bianisotropic Polarizability Tensor of Willis Acoustic Scatterers. Phys Rev B (2018) 98:174305. doi:10.1103/physrevb.98.174305

CrossRef Full Text | Google Scholar

35. Melnikov A, Chiang YK, Quan L, Oberst S, Alù A, Marburg S, et al. Acoustic Meta-Atom with Experimentally Verified Maximum Willis Coupling. Nat Commun (2019) 10:1–7. doi:10.1038/s41467-019-10915-5

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Gong Z, Baudoin M. Three-dimensional Trapping and Assembly of Small Particles with Synchronized Spherical Acoustical Vortices. Phys Rev Appl (2020) 14:064002. doi:10.1103/physrevapplied.14.064002

CrossRef Full Text | Google Scholar

37. Baudoin M, Thomas JL, Sahely RA, Gerbedoen JC, Gong Z, Sivery A, et al. Spatially Selective Manipulation of Cells with Single-Beam Acoustical Tweezers. Nat Commun (2020) 11:1–10. doi:10.1038/s41467-020-18000-y

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Quan L, Ra’di Y, Sounas DL, Alù A. Maximum Willis Coupling in Acoustic Scatterers. Phys Rev Lett (2018) 120:254301. doi:10.1103/physrevlett.120.254301

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Klaseboer E, Sepehrirahnama S, Chan DY. Space-time Domain Solutions of the Wave Equation by a Non-singular Boundary Integral Method and Fourier Transform. J Acoust Soc Am (2017) 142:697–707. doi:10.1121/1.4996860

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Gong Z, Marston PL, Li W. T-matrix Evaluation of Three-Dimensional Acoustic Radiation Forces on Nonspherical Objects in Bessel Beams with Arbitrary Order and Location. Phys Rev E (2019) 99:063004. doi:10.1103/PhysRevE.99.063004

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Chiang YK, Oberst S, Melnikov A, Quan L, Marburg S, Alù A, et al. Reconfigurable Acoustic Metagrating for High-Efficiency Anomalous Reflection. Phys Rev Appl (2020) 13:064067. doi:10.1103/physrevapplied.13.064067

CrossRef Full Text | Google Scholar

42. Jjja D. Exact Solutions for Nondiffracting Beams. I. The Scalar Theory. JOSA A (1987) 4:651–4.

Google Scholar

43. Marston PL. Scattering of a Bessel Beam by a Sphere. J Acoust Soc Am (2007) 121:753–8. doi:10.1121/1.2404931

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Doinikov AA. Acoustic Radiation Interparticle Forces in a Compressible Fluid. J Fluid Mech (2001) 444:1–21. doi:10.1017/s0022112001005055

CrossRef Full Text | Google Scholar

45. Sepehrirahnama S, Chau FS, Lim KM. Effects of Viscosity and Acoustic Streaming on the Interparticle Radiation Force between Rigid Spheres in a Standing Wave. Phys Rev E (2016) 93:023307. doi:10.1103/PhysRevE.93.023307

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Sieck CF, Alù A, Haberman MR. Origins of Willis Coupling and Acoustic Bianisotropy in Acoustic Metamaterials through Source-Driven Homogenization. Phys Rev B (2017) 96:104303. doi:10.1103/physrevb.96.104303

CrossRef Full Text | Google Scholar

47. Chiang YK, Quan L, Peng Y, Sepehrirahnama S, Oberst S, Alù A, et al. Scalable Metagrating for Efficient Ultrasonic Focusing. Phys Rev Appl (2021) 16:064014. doi:10.1103/physrevapplied.16.064014

CrossRef Full Text | Google Scholar

48. Polychronopoulos S, Memoli G. Acoustic Levitation with Optimized Reflective Metamaterials. Sci Rep (2020) 10:1–10. doi:10.1038/s41598-020-60978-4

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Lim G, Wortis M, Mukhopadhyay R. Stomatocyte–discocyte–echinocyte Sequence of the Human Red Blood Cell: Evidence for the Bilayer–Couple Hypothesis from Membrane Mechanics. Proc Natl Acad Sci (2002) 99:16766–9. doi:10.1073/pnas.202617299

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Willis coupling, acoustic levitation, acoustic Bessel beam, acoustic beam forming, biological structures

Citation: Sepehrirahnama S and Oberst S (2022) Acoustic Radiation Force and Torque Acting on Asymmetric Objects in Acoustic Bessel Beam of Zeroth Order Within Rayleigh Scattering Limit. Front. Phys. 10:897648. doi: 10.3389/fphy.2022.897648

Received: 16 March 2022; Accepted: 19 April 2022;
Published: 12 May 2022.

Edited by:

Ashis Sen, Indian Institute of Technology Madras, India

Reviewed by:

Chen Shen, Rowan University, United States
Yan-Feng Wang, Tianjin University, China

Copyright © 2022 Sepehrirahnama and Oberst. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shahrokh Sepehrirahnama, shahrokh.sepehrirahnama@uts.edu.au; Sebastian Oberst, sebastian.oberst@uts.edu.au

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.