Resonance scattering and radiation force calculations for an elastic cylinder using the translational addition theorem for cylindrical wave functions F G Mitri Citation: AIP Advances 5, 097205 (2015); doi: 10.1063/1.4931916 View online: http://dx.doi.org/10.1063/1.4931916 View Table of Contents: http://aip.scitation.org/toc/adv/5/9 Published by the American Institute of Physics AIP ADVANCES 5, 097205 (2015) Resonance scattering and radiation force calculations for an elastic cylinder using the translational addition theorem for cylindrical wave functions F G Mitria Chevron, Area 52 Technology – ETC, Santa Fe, NM 87508, United States (Received 27 July 2015; accepted 16 September 2015; published online 23 September 2015) The standard Resonance Scattering Theory (RST) of plane waves is extended for the case of any two-dimensional (2D) arbitrarily-shaped monochromatic beam incident upon an elastic cylinder with arbitrary location using an exact methodology based on Graf’s translational addition theorem for the cylindrical wave functions The analysis is exact as it does not require numerical integration procedures The formulation is valid for any cylinder of finite size and material that is immersed in a nonviscous fluid Partial-wave series expansions (PWSEs) for the incident, internal and scattered linear pressure fields are derived, and the analysis is further extended to obtain generalized expressions for the on-axis and off-axis acoustic radiation force components The wave-fields are expressed using generalized PWSEs involving the beam-shape coefficients (BSCs) and the scattering coefficients of the cylinder The off-axial BSCs are expressed analytically in terms of an infinite PWSE with emphasis on the translational offset distance d Numerical computations are considered for a zeroth-order quasi-Gaussian beam chosen as an example to illustrate the analysis Acoustic resonance scattering directivity diagrams are calculated by subtracting an appropriate background from the expression of the scattered pressure field In addition, computations for the radiation force exerted on an elastic cylinder centered on the axis of wave propagation of the beam, and shifted off-axially are analyzed and discussed C 2015 Author(s) All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License [http://dx.doi.org/10.1063/1.4931916] I INTRODUCTION The present study concerns the extension of the standard formalism of the acoustic Resonance Scattering Theory (RST) for plane wave incident upon an elastic cylinder,1–3 to the case of 2D acoustical beams of arbitrary character The extension of the classical formalism is necessary, since the fundamental features of the results which are being analyzed here will be present in a large variety of applications involving scattering processes, ranging from the nondestructive evaluation and imaging of elastic cylinders, to biomedical and biophysical applications using finite beams as opposed to plane waves of infinite extent It is of some importance to develop a formalism which accounts for the nature and finite character of the incident field, since the resulting resonance scattering process may be enhanced or suppressed depending on the beam’s parameters, as recently demonstrated for the case of an elastic sphere (or a shell) placed arbitrarily in an acoustical beam of arbitrary wave-front.4,5 Nevertheless, analyzing the arbitrary scattering with a cylinder using the formalism devoted to spheres4,5 as a first approximation, leads to significant inaccuracies in the numerical predictions Therefore, it is of particular importance to develop a method applicable to cylinders in order to numerically predict and compute the arbitrary scattering, and other phenomena such as the resulting acoustic radiation force for various applications a Electronic mail: F.G.Mitri@ieee.org 2158-3226/2015/5(9)/097205/12 5, 097205-1 © Author(s) 2015 097205-2 F G Mitri AIP Advances 5, 097205 (2015) Despite the extensive analyses related to plane wave scattering by an infinitely-long cylinder,2,6–9 earlier works considered the case of beams, based on the angular spectrum decomposition of plane waves, requiring the numerical evaluation of indefinite integrals.10,11 On the other hand, the present analysis is based on the partial-wave series expansion (PWSE) method (known also as normal-mode decomposition in Fourier series) in cylindrical coordinates, and the evaluation of the beam-shape coefficients (BSCs) stemming from Graf’s additional theorem for the cylindrical wave functions, without the need of numerical integration procedures, used previously in the method of the angular spectrum decomposition into plane waves, or in the computation of the acoustic radiation force on a rigid (sound impenetrable) cylinder.12 It is also worth mentioning that the formalisms for the electromagnetic scattering of optical beams by (dielectric or perfectly conducting) cylindrical objects have used the angular spectrum and PWSE methods with numerical integration in several studies.13–18 Nonetheless, the extension to predict and numerically compute the resonance scattering as well as the resulting acoustic radiation force using an exact method such as Graf’s additional theorem, remains to be accomplished, which is developed here in this investigation The analysis is based on the PWSE method in cylindrical coordinates and the evaluation of the on- and off-axis beam-shape coefficients (BSCs) using the translational addition theorem Normalized form functions are derived and acoustic (resonance) scattering directivity diagrams are calculated by subtracting an appropriate background to isolate the pure resonances of an aluminum cylinder placed arbitrarily in the field of an acoustical zeroth-order quasi-Gaussian beam, chosen as an example Numerical predictions for the axial and transverse acoustic radiation force components are also provided Particular emphasis is given on the effect of shifting the cylinder arbitrarily with respect to the incident waves, as well as the focusing properties of the illuminating beam II GENERALIZED THEORY OF RESONANCE SCATTERING FOR A CYLINDER IN A NON-VISCOUS FLUID Consider an acoustical beam propagating in a nonviscous fluid of density ρ and a speed c, and incident upon an infinitely-long cylinder of radius a and density ρc The center of the cylinder coincides with the origin of a cylindrical coordinate system (r, θ, z), and the incident beam is of arbitrary shape (Fig 1) The acoustic field is described by its incident pressure field Pinc that is a solution of the Helmholtz wave equation, ∇2 + k Pinc = 0, (1) where the wave number k = ω/c, ω is the angular frequency, and c is the compressional speed of sound in the surrounding fluid For the case of a monochromatic beam, the most general separation of variables (non-singular) solution of (1) in cylindrical coordinates is (p 430 in Ref 19), Pinc (r, θ, z) = P0 +∞ bn Jn (k r r) ei(k z z+nθ), (2) n=−∞ where P0 is the pressure amplitude in the absence of the cylinder, bn are the BSCs that will be determined subsequently, Jn (·) is the cylindrical Bessel function of the first kind, k r and k z are the radial and axial wave-numbers, respectively, such that k = k r2 + k z2.20 A time-harmonic variation in the form of e−iωt is assumed, but suppressed for convenience from (2) since the space-dependent pressure field is only considered It is important to emphasize that (2) (and the subsequent equations) is only applicable to 2D beams, with a discrete (single-valued) k z For the more general case of beams bounded also in the z-direction, the mathematical expression for the incident wave-field must include an integration over all possible values of k z (i.e., continuous spectrum) (See Eq.(4.46) in Ref 20) Regarding the determination of the BSCs bn , there exist two different methods that could be used for this purpose In the first method, bn can be obtained after integrating both sides of (2) using 097205-3 F G Mitri AIP Advances 5, 097205 (2015) FIG The schematic describing the interaction of a 2D arbitrary-shaped acoustical beam with a cylinder of radius a of arbitrary location immersed in a non-viscous fluid The cylindrical coordinate system (r, θ, z) is referenced to the center of the cylinder, while the system of coordinates (r 0, θ 0, z 0) is referenced to the center of the beam the orthogonality relationship 2π ′ ei(n−n )θ dθ = 2πδ nn′, (3) where δ i j is the Kronecker delta function, so that the BSCs are expressed after some algebraic manipulation as, e−ik z z bn = 2πP0 Jn (k r r) 2π Pinc (r, θ, z) e−inθ dθ (4) Since the BSCs describe the incident beam’s characteristics in the cylindrical coordinate system independently of the presence of the cylinder and its size, it is important to note from Eq (4) that the axial wave number with discrete value k z is solely determined by the form of the incident pressure field in order to cancel out any z-dependence in the BSCs Moreover, the integral must be proportional to Jn (k r r) to eliminate any r-dependence as well From a numerical analysis standpoint, however, the apparent r-dependence in the denominator in Eq (4) may lead to mathematical indeterminacies if k r r is not adequately chosen so as to avoid the zeros of Jn (k r r) A way to circumvent this difficulty requires evaluating the BSCs over a virtual control cylindrical surface of radius Rc that encloses the cylinder of radius a such that Rc ≫ a This approach may only be used for beams that satisfy (1) and propagate in an isotropic fluid medium in which attenuation and diffraction effects in the space between the virtual cylindrical surface and the cylinder are negligible Therefore, the BSCs evaluated at r = Rc are equivalent to those evaluated at r = a Note that evaluation of the BSCs as given by Eq (4) requires numerical integration21 for a given incident pressure field Pinc (r, θ, z) In the second method, which will be used here, evaluation of the BSCs is accomplished using the exact method of Graf’s translational addition theorem for the cylindrical wave functions (Ch in Ref 22) The method requires deriving the off-axial BSCs bn in the coordinates system centered 097205-4 F G Mitri AIP Advances 5, 097205 (2015) on the cylinder from a known expression of the BSCs (denoted by b p ) in the coordinates system centered on the incident beam’s axis At point M (see Fig 1), the incident field, written in the system of coordinates (r 0, θ 0, z0), can be expressed as, Pinc (r 0, θ 0, z0) = P0 +∞ b p Jp (k r r 0) ei(k z z 0+pθ0) (5) p=−∞ Based on Graf’s addition theorem (Ch in Ref 22) and the geometry for the distances and angles as shown in Fig 1, the function Jp (k r r 0) ei pθ0 in (5) can be expanded as a PWSE as, Jp (k r r 0) ei pθ0 = +∞ Jn−p (k r d) Jn (k r r) ei(n−p)(θ−φ)ei pθ (6) n=−∞ Substituting (6) into (5), and equating the result with (2), leads to the expression for the off-axial BSCs bn in terms of the axial BSCs b p as, bn = +∞ b p Jn−p (k r d) e−i(n−p)φ , (7) p=−∞ where d is the distance that separates systems, φ is termed the shift-angle, ) the two coordinates ( −1 y0−yoff 2 x off + yoff , where x off and yoff determine the which is defined as φ = tan x 0−x off , and d = amount of offset in the axial and lateral directions, respectively Eq (7) constitutes the generalized result obtained without any approximations and without numerical integration It can be applicable to any 2D beam for which a closed-form expression for the on-axis BSCs is known, in order to evaluate the off-axis BSCs in the translated system of coordinates Wave propagation in the cylinder’s core material is represented by suitable solutions of the equation of motion of a solid elastic medium, which may be written as23 ∂ 2U c , (8) ∂t where λ c and µc are the Lamé coefficients The parameter Uc represents the vector displacement that is expressed as a sum of the gradient of a scalar potential Φc representing longitudinal waves and the curl of a vector potential Ψc representing shear waves as follows, (λ c + 2àc ) 2Uc + àc ì ( ì Uc ) = ρc Uc = ∇Φc + (∇ × Ψc ) (9) The displacement equations are satisfied if the potentials Φc and Ψc satisfy the Helmholtz equations for the solid medium, ∇2 + k 2L Φc = 0, (10) ∇ + (11) k S2 Ψc = 0, where k L = ω/cL = ω/[(λ c + 2µc ) /ρc ]1/2, and k S = ω/cS = ω/[µc /ρc ]1/2, refer to the longitudinal and shear wave numbers in the cylinder, respectively Generally, the vector potential Ψc (ψr , ψθ , ψ z ) (describing horizontally polarized (SH) and vertically polarized (SV) shear-waves) has three scalar components,24,25 but to ensure the uniqueness of the solution of (8), the vector potential has to satisfy the gauge invariance ∇ · Ψc = 0,26,27 thus, it is a solenoidal vector field This allows defining appropriate potentials related to the longitudinal and shear waves inside the cylinder, which are represented in cylindrical coordinates by, +∞ P0 bn An Ωn (κ Lr) ei(k z z+nθ), ρc ω2 n=−∞ (12) +∞ P0 bn Bn Ωn+1 (κ Sr) ei(k z z+nθ), i ρc ω2 n=−∞ (13) Φc (r, θ, z) = ψr (r, θ, z) = 097205-5 F G Mitri AIP Advances 5, 097205 (2015) ψθ (r, θ, z) = − ψ z (r, θ, z) = +∞ P0 bn Bn Ωn+1 (κ Sr) ei(k z z+nθ), ρc ω2 n=−∞ +∞ P0 bn Dn Ωn (κ Sr) ei(k z z+nθ), i ρc ω2 n=−∞ (14) (15) where κ L and κ S are defined in Table I, and the functions Ωn (·) represent the Bessel functions of the first kind Jn (·), or the modified Bessel functions of the first kind In (·) respectively, according to the values of α = arctan(k z /k r ) Eqs (13) and (14) show that both potentials ψr and ψθ depend on the coefficients Bn , which is anticipated from the property of the gauge invariance Note also that instead of finding the three components of the vector potential Ψc , one may express it in terms of a pair of scalar (Debye) potentials (Ch 13 in Ref 28),23,29 to obtain directly two scalar potentials representing the SH and SV waves, respectively.30 This method has been employed in the generalized formalism developed for spheres.4 It is emphasized that both methods are commensurate with the same result as long as an elastic isotropic cylinder is considered Nevertheless, for an anisotropic cylinder, different equations have been established when decomposing the vector potential in terms of scalar potentials (See Section in Ref 31) Upon the interaction of the incident beam with the cylinder, outgoing cylindrically diverging scattered waves are produced in the surrounding medium, which can be represented by Psc (r, θ, z) = P0 +∞ bnCn Hn(1) (k r r) ei(k z z+nθ), (16) n=−∞ where Hn(1) (·) is the cylindrical Hankel function of the first kind of order n, and Cn are the scattering coefficients that are determined by applying the boundary conditions of continuity of radial stresses and displacements, and the nullity of the tangential stresses since the surrounding fluid is nonviscous Denoting the differential operators ∂/∂ (r, θ, z) by ∂r,θ, z , and using the expressions for the potentials given by Eqs (12)-(15), the components of the displacements (ur , uθ , uz ) and stresses (σr r , σr θ , σr z ) are written, respectively, as, ur = ∂r Φc + ∂θ ψ z − ∂z ψθ , r uθ = ∂θ Φc + ∂z ψr − ∂r ψ z , r uz = ∂z Φc + [∂r (rψθ ) − ∂θ ψr ] , r σr r = λ c ∆ + 2µc ∂r ur , σr θ = µc ∂r uθ + (∂θ ur − uθ ) , r (17) (18) (19) (20) (21) TABLE I The expression for the parameters κ L , κ S , Ωn , δ and γ appearing in the matrix elements of the scattering coefficients C n according to the longitudinal and shear waves coupling angles θ L, S = sin−1(c/c L, S ) κL κS Ω n (κ L a) Ω n (κ S a) δ γ α ≤ θL k − k z2 , L k S2 − k z2 , θL ≤ α < θS 2, k z2 − k L k S2 − k z2 , θS ≤ α 2, k z2 − k L k z2 − k S2 , J n (κ L a) J n (κ S a) 1 I n (κ L a) J n (κ S a) −1 I n (κ L a) I n (κ S a) −1 −1 097205-6 F G Mitri AIP Advances 5, 097205 (2015) σr z = µc (∂z ur + ∂r uz ) , (22) where ∆ in (20) is given by, ∆ = ∂r ur + (∂θ uθ + ur ) + ∂z uz r (23) At the interface fluid-cylinder for r = a, the following boundary conditions are applied, • continuity of the radial displacement, ur | r =a = ∂r (Pinc + Psc) , ρω r =a (24a) • continuity of the radial stress, σr r | r =a = − (Pinc + Psc)| r =a , (24b) • nullity of the tangential and shear stresses, σr θ | r =a = 0, (24c) σr z | r =a = (24d) These boundary conditions lead to four linear equations with four unknowns coefficients An , Bn , Cn and Dn The general solution for the scattering coefficients Cn for the scattered pressure field in the fluid medium, which appear in Eq (16) is given by Cn = M1 M2 det M3 M4 m12 m22 m32 m42 m13 m23 m33 m43 m14 m24 m34 m44 m11 m21 det m31 m41 m12 m22 m32 m42 m13 m23 m33 m43 m14 m24 m34 m44 , (25) where Mi and m j k are the dimensionless elements of the determinants given explicitly in Appendix A of Ref 32 These coefficients are found to equal those obtained previously from the study of the acoustic resonance scattering of plane progressive waves by an elastic infinite cylinder at oblique incidence.32–35 Note also that the interior (resonance) compressional field of the cylinder which depends on An , as well as the interior (resonance) SH and SV fields which depend on Bn and Dn , respectively, can be determined in a straightforward manner by solving the matrix equation m j k {Cn , An , Bn , Dn } = {M1, M2, M3, M4} by Cramer’s rule The corresponding interior displacement potential fields for the compressional and shear waves inside the elastic cylinder were evaluated in the context of the standard RST for cylinders.36 The plots for the resonating potential fields inside the elastic scattering cylinder displayed “pure” resonance peaks with no interference with any specular reflection background A Generalized Total Form Function in the Far-Field In the region far from the cylinder, a dimensionless representation of the acoustic scattered pressure [given by (16)] in the surrounding fluid medium is defined by a generalized (steady-state) total form function as, (2r/a) [Psc (r, θ, z) /P0] e−i(k z z+k r r ) (26) f ∞ (ka, θ) = kr→ ∞ 097205-7 F G Mitri AIP Advances 5, 097205 (2015) limit for the cylindrical Hankel function of the first kind, Hn(1) (k r r) Using thei(kasymptotic ≈ 2/(πk r r) e r r −nπ/2−π/4), (26) is expressed as, f ∞ (ka, θ) = √ +∞ i −n bnCn einθ iπk r a n=−∞ (27) B Generalized Resonance Form Function in the Far-Field The generalized total form function (27) is obtained as a PWSE, based on the linear theory of elastodynamics.23 When the form function is plotted as a function of the dimensionless size parameter ka, variations in the form of resonance peaks or dips appear The standard RST formalism for plane waves1 showed that resonance peaks appear because of a constructive interference of inner elastic (resonance) waves with inelastic (non-resonant) outer surface waves circumnavigating the cylinder’s surface in the surrounding fluid (also known as Franz’ waves7), which cause the specular reflection echo in experimental signals.37 To isolate the pure resonances, the standard RST suggested that appropriate backgrounds be subtracted from the total form function.38 In particular, it has been found that for dense cylinder materials, the background is closely approximated by the modal amplitudes of rigid-body scattering, whereas for very soft materials, it approaches the soft-body scattering amplitude However, for materials with density comparable to the external fluid, (or for example a dense but quite thin cylindrical shell material), an intermediate background38 has been adequately defined and subtracted from the total form function to properly identify the resonances The main purpose of identifying and isolating resonances is to carry out remote classification of a target in underwater applications, or to characterize some of the physical and mechanical properties of the scatterer The generalized resonance form function for a cylinder (or cylindrical shell) is defined by, f ∞res (ka, θ) = √ +∞ r, s, i inθ e , i −n bn Cn − Cn iπk r a n=−∞ (28) where Cnr, s,i are the scattering coefficients for a perfectly rigid, a perfectly soft or an intermediate cylinder material The scattering coefficients for the rigid and soft cylinder are given, respectively, ′ by,7 Cnr = −Jn′ (k r a) /Hn(1) (k r a), where the prime indicates a derivative with respect to the argument, and Cns = −Jn (k r a) /Hn(1) (k r a) For near-field scattering calculations, the process of isolation of the resonances can be readily applied to (16), in a procedure similar to what is performed in Eq (28) III ACOUSTIC RADIATION FORCE COMPONENTS FOR A CYLINDER IN A 2D ARBITRARY-SHAPED BEAM IN A NON-VISCOUS FLUID Now that the expressions for the incident and scattered fields are established, evaluation of the acoustic radiation force components becomes possible In a non-viscous fluid, it proves advantageous to consider the far-field scattering in the evaluation of the force39 that is expressed as 2π ⟨F⟩ = ρk r ℜ {Φi s } dS, (29) kr→ ∞ where, Φi s = Φ∗sc [(i/k r ) ∂r Φinc − Φinc − Φsc] (30) The incident and scattered velocity potentials Φinc and Φsc in (30) can be expressed in terms of the incident and scattered pressure fields, given by (2) and (16), respectively, such that Φ{inc, sc} = P{inc, sc}/(i ρω) The cylinder’s differential surface is dS = nLr dθ where L is the characteristic length of the cylinder The outward normal unit vector is n = cos θe x + sin θe y , where e x and e y are the unit vectors in the Cartesian coordinates system, the symbol ⟨·⟩ denotes time-averaging, the superscript * denotes a complex conjugate 097205-8 F G Mitri AIP Advances 5, 097205 (2015) Taking the asymptotic limits for Eq (2) and Eq (16) in the far-field scattering limit (i.e kr → ∞), and substituting them into Eq (29) using Eq (30) as well as the property of the following angular integral, (δ n, n+1 + δ n, n−1) cos θ ′ , dθ = π ei(n −n)θ (31) (δ ) i − δ sin θ n, n+1 n, n−1 where δij is the Kronecker delta function, and manipulating the result, the on-axis and off-axis components of the acoustic radiation force are expressed, respectively, as, 2π ex F x = ⟨F⟩ · ey Fy Yx = Yy Sc E0, (32) where Sc = 2aL is the cross-sectional surface of the cylinder, E0 = 21 ρk r2|φ0|2 is a characteristic energy density, and the functions Y{x, y} are the on- and off-axis radiation force functions, that normalize the radiation force components Their expressions are given after some algebraic manipulation, respectively, by, +∞ ∗ ∗ ∗ ∗ , (1 ) ℑ b + C b C − b C (33) n n n+1 n+1 n−1 n−1 k r a n=−∞ +∞ ∗ ∗ ∗ ∗ , (1 ) ℜ b + C b C + b C (34) Yy = − n n n+1 n+1 n−1 n−1 k r a n=−∞ where the symbols ℜ{.} and ℑ{.} denote the real and imaginary parts of a complex number, respectively Yx = IV NUMERICAL RESULTS AND DISCUSSIONS A Numerical Validation The analysis is started by validating a numerical code that is developed to compute a normalized form function, f norm,∞ (ka, θ) = +∞ −n i bnCn einθ , n=−∞ (35) based on the aforementioned generalized acoustic scattering theory of an arbitrary-shaped beam incident upon a cylinder This normalized function (35) is equivalent to a simplified expression given previously for the case of plane waves (p 411 in Ref 6), which allows direct comparison and validation of the results with those given in Ref The truncation order Nmax for the index n in the series is chosen such that Cn+Nmax/C0 ∼ 10−8, for n = 1, 2, , which yields a negligible truncation error As an initial test of the arbitrary incident beam theory and the computer program, the case of infinite plane wave incidence40 is assumed, such that b p = i p The computed BSCs for the infinite plane wave case using (7) are used to evaluate the magnitude of the normalized form function (35), and compare it with the result obtained using the standard RST f norm,∞ (ka, θ) is evaluated for an aluminum elastic cylinder (ρc = 2700 kg/m3, cL = 6200 m/s, cS = 3170 m/s) immersed in water (ρ = 1000 kg/m3, c = 1470 m/s) for ka = 5, and the results are displayed in Fig As noticed, perfect agreement is observed An additional test of the validity of the arbitrary incident beam theory and the computer program considered an incident plane wave such that the cylinder is shifted off-axially (along the x-direction or the y-direction) The computed scattering directivity diagram 097205-9 F G Mitri AIP Advances 5, 097205 (2015) FIG The plot for the modulus of the normalized form function given by Eq (35) for an aluminum cylinder immersed in water at ka = in infinite plane waves with k z = The plot is computed using the standard elastic scattering theory (See Fig 12 in Ref 6), while the dotted circles correspond to the computations obtained using the translational addition theorem Perfect agreement is observed The arrow on the left-hand side denotes the direction of the incident waves showed also perfect agreement with Fig 2, such that the shift did not affect the scattering, as expected for the plane wave case B Off-Axial (Resonance) Scattering of a Cylindrical Quasi-Gaussian Beam by an Elastic Aluminum Cylinder The present theory is now illustrated by considering numerical computations for the scattering directivity patterns of a zeroth-order quasi-Gaussian beam with k z = and a dimensionless waist kw0 = 7, chosen as an example, and corresponding to a focused beam (See Fig in Ref 41) The total and resonance form-functions, given respectively by Eq (28) and Eq (29) are numerically evaluated, such that the on-axis BSCs40 are expressed as b p = i p I p (k r x R )/I0(k r x R ) Fig 3(a) displays the polar scattering directivity patterns for an aluminum cylinder centered on the beam’s axis for ka = The directivity pattern is affected by the choice of kw0, such that the scattering lobes around θ = 30◦ and 330◦ observed in Fig are reduced in Fig 3(a) The effect of shifting the cylinder off-axially along the y−direction in arbitrary units is shown in panel (b) for yoff = 0.3, while panel (c) corresponds to yoff = 0.6 The resulting asymmetry in the directivity patterns is obvious when the cylinder is shifted off the beam’s axis Moreover, the patterns in (b) and (c) are asymmetric with respect to one another as expected Panels (d)-(f) show the corresponding resonance directivity patterns, which clearly show a quadrupole resonance vibrational mode for the aluminum cylinder at the selected ka Note that the radiating pattern is affected by the translational shift of the beam from the center of the cylinder The analysis is further extended to calculate the total backscattering (θ = 180◦) form function as given by Eq (27) as well as its resonance counterpart, given by Eq (28), for the aluminum elastic cylinder placed in the field of a zeroth-order quasi-Gaussian beam Magnitude plots are evaluated in the non-dimensional frequency and transverse shift ranges defined, respectively, as < ka ≤ 10, and −5 ≤ yoff ≤ Panels (a) and (b) of Fig display the results for | f ∞ (ka, θ = π)| and f ∞res (ka, θ = π) , respectively, for kw0 = In panels (a), (b), the maxima in the plots are observed at the center of the beam (i.e., yoff = 0) Moreover, panel (b) clearly shows the effects of background subtraction such that the pure resonances appear clearly in the plot 097205-10 F G Mitri AIP Advances 5, 097205 (2015) FIG Panel (a) shows the directivity pattern for the on-axis scattering of a zeroth-order quasi-Gaussian beam by an elastic aluminum cylinder submerged in water at k a = 5, and kw0 = Panel (b) shows the effect of shifting the cylinder along the transverse direction by an arbitrary shift yoff = 0.3 One notices the asymmetry in the scattering as a result of the shift Panel (c) shows the effect of shifting the cylinder further along the (positive) transverse direction by an arbitrary shift yoff = 0.6 As in panels (a)-(c), panels (d)-(e) correspond to the resonance directivity diagrams for the aluminum cylinder, revealing a quadrupole resonance vibrational mode The arrows on the left-hand side of each of the panels denote the direction of the beam incidence FIG Panel (a) shows the magnitude plot for the total backscattering (θ = π) form function as given by Eq (27), whereas panel (b) displays its resonance counterpart as given by Eq (28) for an aluminum elastic cylinder immersed in water and placed in the field of a zeroth-order quasi-Gaussian beam with kw0 = The axis yoff = corresponds to a zero offset; that is, the beam is centered on the cylinder C Axial and transverse acoustic radiation force components The axial and transverse acoustic radiation force functions given by Eqs (33), (34) are evaluated numerically for the aluminum elastic cylinder, illuminated by a zeroth-order quasi-Gaussian beam with a non-dimensional beam waist kw0 = Similarly to Fig 4, the non-dimensional frequency and transverse shift ranges defined, respectively, as < ka ≤ 10, and −5 ≤ yoff ≤ 5, are considered Panels (a) and (b) of Fig display the results for the axial Yx and transverse Yy radiation force functions, respectively As observed in panel (a), the Yx plot is maximal at the center of the beam such that yoff = Due to symmetry considerations in this limit, the transverse radiation force function Yy vanishes, as shown in panel (b) of Fig Note that the plot for Yx is always positive and symmetric with respect to the axis yoff = 0, whereas the plot for Yy is asymmetric therein, and varies between positive and negative values determined by the amount of the shift, ka and kw0 Since the axial radiation force function Yx > 0, there is no axial negative (pulling) force acting upon 097205-11 F G Mitri AIP Advances 5, 097205 (2015) FIG Panel (a) shows the axial radiation force function plot and panel (b) displays the transverse component for an aluminum elastic cylinder immersed in water and placed in the field of a zeroth-order quasi-Gaussian beam with kw0 = Note that the plot for the axial component shown in panel (a) is symmetric with respect to the offset axis yoff = 0, whereas the transverse component shown in panel (b) is asymmetric with respect to it the cylinder, which was also confirmed in Ref 41 Consequently an axial trapping cannot be accomplished with the zeroth-order cylindrical quasi-Gaussian beam This has also been demonstrated for the spherical quasi-Gaussian beam incident upon a rigid or an elastic sphere.42,43 Nonetheless, since the transverse radiation force function Yy alternates between positive and negative values as yoff and ka take non-zero values, it may be feasible to produce a lateral tweezing with the zeroth-order cylindrical quasi-Gaussian beam V CONCLUSION In this work, the standard Resonance Scattering Theory (RST) of plane waves is extended for the case of any two-dimensional (2D) arbitrarily-shaped monochromatic beam incident upon an elastic cylinder with arbitrary location using an exact methodology based on Graf’s translational addition theorem for the cylindrical wave functions The analysis is exact as it does not require numerical integration procedures The method is applicable to any beam that is an exact solution of Helmholtz equation for any (elastic) cylinder of finite radius located arbitrarily in the beam The present analysis provides an important tool to characterize the arbitrary scattering and compute the acoustic radiation force components on an elastic cylinder illuminated by beams of arbitrary wavefront Appropriate PWSEs for the incident and scattered pressure fields are derived, and the off-axis BSCs are obtained stemming from the on-axis BSCs using the translational addition theorem for cylindrical wave functions Numerical computations illustrate the analysis with particular emphasis on the offset of the cylinder with respect to the incident beam The present benchmark analytical solutions could be used to validate those obtained by strictly numerical methods using the FEM, the FDTDM or other numerical tools The method for predicting the arbitrary acoustic scattering and acoustic radiation force finds important applications in various fields in acoustics especially when the physical phenomena under investigation involve the interaction of cylindrical particles with arbitrary-shaped 2D acoustical beams L Flax, G.C Gaunaurd, and H Überall, in Phys Acoustics, edited by W.P Mason (Academic Press, New York, 1981), Vol XV, pp 191-294 G.C Gaunaurd, “Elastic and acoustic resonance wave scattering,” Applied Mechanics Reviews 42, 143-192 (1989) N D Veksler, Resonance Acoustic Spectroscopy (Springer-Verlag, Berlin, 1993), p 282 F G Mitri, “Generalized theory of resonance excitation by sound scattering from an elastic spherical shell in a nonviscous fluid,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control 59, 1781-1790 (2012) F G Mitri, “Generalized Theory of Resonance Scattering (GTRS) Using the Translational Addition Theorem for Spherical Wave Functions,” IEEE Trans Ultrason Ferroelectr Freq Control 61, 1880-1888 (2014) J.J Faran, “Sound scattering by solid cylinders and spheres,” J Acoust Soc Am 23, 405-418 (1951) H Überall, in Phys Acoustics, edited by W.P Mason (Academic Press, New York, 1973), Vol X, pp 1-60 D L Jain and R P Kanwal, “Scattering of acoustic, electromagnetic and elastic SH waves by two-dimensional obstacles,” Annals of Physics 91, 1-39 (1975) G.C Gaunaurd and M F Werby, “Acoustic resonance scattering by submerged elastic shells,” Applied Mechanics Reviews 43, 171-208 (1990) 10 S Zeroug and L B Felsen, “Nonspecular Reflection of 2-Dimensional and 3-Dimensional Acoustic Beams from FluidImmersed Cylindrically Layered Elastic Structures,” J Acoust Soc Am 98, 584-598 (1995) 097205-12 11 F G Mitri AIP Advances 5, 097205 (2015) T J Cloutier, A Safaeinili, D E Chimenti, S Zeroug, and L B Felsen, “Ultrasonic beam reflection from fluid-loaded cylindrical shells,” J Appl Phys 83, 2408-2419 (1998) 12 M Azarpeyvand and M Azarpeyvand, “Acoustic radiation force on a rigid cylinder in a focused Gaussian beam,” J Sound Vib 332, 2338-2349 (2013) 13 N G Alexopoulos and P Pyong, “Scattering of waves with normal amplitude distribution from cylinders,” IEEE Transactions on Antennas and Propagation 20, 216-217 (1972) 14 T Kojima and Y Yanagiuchi, “Scattering of an offset two-dimensional Gaussian beam wave by a cylinder,” J Appl Phys 50, 41-46 (1979) 15 S Kozaki, “Scattering of a Gaussian beam by an inhomogeneous dielectric cylinder,” Journal of the Optical Society of America 72, 1470-1474 (1982) 16 G Gouesbet and G Gréhan, “Interaction Between Shaped Beams and an Infinite Cylinder, including a discussion of Gaussian beams,” Particle & Particle Systems Characterization 11, 299-308 (1994) 17 E Zimmermann, B Krattiger, R Dändliker, and N Souli, “Scattering of an off-axis Gaussian beam by a dielectric cylinder compared with a rigorous electromagnetic approach,” Journal of the Optical Society of America A 12, 398-403 (1995) 18 J A Lock, “Scattering of a diagonally incident focused Gaussian beam by an infinitely long homogeneous circular cylinder,” Journal of the Optical Society of America A 14, 640-652 (1997) 19 E Skudrzyk, The Foundations of Acoustics (Springer-Verlag, Wien, New York, 1971), p 790 20 E G Williams, Fourier Acoustics (Academic Press, London, 1999), pp 115-148 21 F G Mitri, “Interaction of an acoustical 2D-beam with an elastic cylinder with arbitrary location in a non-viscous fluid,” Ultrasonics 62, 244-252 (2015) 22 M Abramowitz and I A Stegun, Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965) 23 A.C Eringen and E.S Suhubi, Elastodynamics: Linear Theory (Academic Press, New York-London, 1975) 24 T.R Meeker and A H Meitzler, in Phys Acoustics, edited by W.P Mason (Academic Press, New York, 1964), Vol I, Pt A 25 X Bao, H Cao, and H Überall, “Resonances and surface waves in the scattering of an obliquely incident acoustic field by an infinite elastic cylinder,” J Acoust Soc Am 87, 106-110 (1990) 26 L Flax, V K Varadan, and V.V Varadan, “Scattering of an obliquely incident acoustic wave by an infinite cylinder,” J Acoust Soc Am 68, 1832-1835 (1980) 27 T Li and M Ueda, “Sound scattering of a plane wave obliquely incident on a cylinder,” J Acoust Soc Am 86, 2363-2368 (1989) 28 P.M Morse and H Feshbach, Methods of theoretical physics (McGraw-Hill Book Co., New York, 1953) 29 C G Gray, “Multipole expansions of electromagnetic fields using Debye potentials,” American Journal of Physics 46, 169-179 (1978) 30 W Huang, Y J Wang, and S I Rokhlin, “Oblique scattering of an elastic wave from a multilayered cylinder in a solid Transfer matrix approach,” J Acoust Soc Am 99, 2742-2754 (1996) 31 F Honarvar, E Enjilela, A N Sinclair, and S A Mirnezami, “Wave propagation in transversely isotropic cylinders,” International Journal of Solids and Structures 44, 5236-5246 (2007) 32 F G Mitri, “Acoustic backscattering enhancements resulting from the interaction of an obliquely incident plane wave with an infinite cylinder,” Ultrasonics 50, 675-682 (2010) 33 F Léon, F Lecroq, D Décultot, and G Maze, “Scattering of an obliquely incident acoustic wave by an infinite hollow cylindrical shell,” J Acoust Soc Am 91, 1388-1397 (1992) 34 J.M Conoir, P Rembert, O Lenoir, and J.L Izbicki, “Relation between surface helical waves and elastic cylinder resonances,” J Acoust Soc Am 93, 1300-1307 (1993) 35 P L Marston, “Approximate meridional leaky ray amplitudes for tilted cylinders: End-backscattering enhancements and comparisons with exact theory for infinite solid cylinders,” J Acoust Soc Am 102, 358-369 (1997) 36 L Flax, P K Raju, J George, and H Überall, “Resonating fields inside elastic scattering objects,” The Journal of the Acoustical Society of America 68, 973-979 (1980) 37 F G Mitri, J F Greenleaf, Z E A Fellah, and M Fatemi, “Investigating the absolute phase information in acoustic wave resonance scattering,” Ultrasonics 48, 209-219 (2008) 38 J.D Murphy, J George, and H Überall, “Isolation of the resonant component in acoustic scattering from fluid-loaded cylindrical shells,” Wave Motion 1, 141-147 (1979) 39 F G Mitri, “Theoretical and experimental determination of the acoustic radiation force acting on an elastic cylinder in a plane progressive wave - far-field derivation approach,” New Journal of Physics 8, 138 (2006) 40 F.G Mitri, “Optical theorem for two-dimensional (2D) scalar monochromatic acoustical beams in cylindrical coordinates,” Ultrasonics 62, 20-26 (2015) 41 F G Mitri, Z E A Fellah, and G T Silva, “Pseudo-Gaussian Cylindrical Acoustical Beam - Axial Scattering and Radiation Force on an Elastic Cylinder,” J Sound Vib 333, 7326–7332 (2014) 42 F G Mitri, “Interaction of an acoustical Quasi-Gaussian beam with a rigid sphere: linear axial scattering, instantaneous force, and time-averaged radiation force,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control 59, 2347-2351 (2012) 43 F G Mitri and Z E A Fellah, “Mechanism of the quasi-zero axial acoustic radiation force experienced by elastic and viscoelastic spheres in the field of a quasi-Gaussian beam and particle tweezing,” Ultrasonics 54, 351-357 (2014)