Skip to main content

METHODS article

Front. Phys., 02 June 2023
Sec. Physical Acoustics and Ultrasonics
This article is part of the Research Topic Ultrasound Micromanipulations and Ocean Acoustics: From Human Cells to Marine Structures View all 11 articles

Acoustic scattering of a pair of rigid spheroids based on the T-matrix method

  • 1School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan, Hubei, China
  • 2School of Naval Architecture, Ocean and Energy Power Engineering, Wuhan University of Technology, Wuhan, China
  • 3Hubei Key Laboratory of Naval Architecture & Ocean Engineering Hydrodynamics (HUST), Wuhan, Hubei, China
  • 4Collaborative Innovation Centre for Advanced Ship and Deep-Sea Exploration (CISSE), Shanghai, China

In this study, the T-matrix method combined with the addition theorems of spherical basis functions is applied to semi-analytically compute the underwater far-field acoustic scattering of a pair of rigid spheroids with arbitrary incident angles. The involvement of the addition theorems renders the multiple scattering fields of each spheroid to be translated into an identical origin. The accuracy and convergence property of the proposed method are verified and validated. The interference of specular reflection wave and Franz wave can be spotted from the oscillations of the form function. Furthermore, the propagation paths of specular reflection and Franz waves are quantitatively analyzed in the time domain with conclusions that the Franz waves reach the observation point subsequent to specular reflection waves and the time interval between these two wave series is equal to the time cost of the Franz waves traveling along the sphere surfaces. Finally, the effects of separation distances, aspect ratios (the ratio of the polar radius to equatorial radius), non-dimensional frequencies, and incidence angles of the plane wave on the far-field acoustic scattering of a pair of rigid spheroids are studied by the T-matrix method.

1 Introduction

The study of underwater target acoustic scattering gains attentions from researchers, and the relative research studies are widely applied in engineering practices such as underwater target detection, positioning, imaging, and underwater communication. The mechanism of multiple-target acoustic scattering is more complex than that of single-target scattering due to the existence of multiple scattering. In this study, a pair of rigid spheroids is chosen as the target to investigate the multiple acoustic scattering characteristics.

In the past decades, a series of numerical and analytical methods are proposed to solve the underwater acoustic scattering problem. Numerical methods such as the finite element method (FEM) and boundary element method (BEM) can solve acoustic scattering problems under complex conditions [13]. Recently, various extended methods, like the smooth finite element method and meshfree method [47], are proposed to solve the underwater acoustic scattering issues. However, the computational efficiency of those numerical methods decreases as the frequency increases because of the requirement of very dense meshes. Compared to the numerical methods, the analytical method can provide precise solutions with a faster convergence speed. In addition, the physical understanding of the acoustic scattering wave can be explained by the analytical solutions. Rayleigh first derived the Bessel–Legendre series (mathematical) solution for the acoustic scattering of a sphere by the variable separation method. However, his works are only competent to cases with small non-dimensional frequencies ka [8]. Based on Rayleigh’s work, various extended analytical methods, like the partial wave series (PWS) method and Sommerfeld–Watson transformation (SWT), are proposed to obtain the mathematical expressions of acoustic scattering of the single simple target (such as the sphere and circular cylinder) [911]. Záviška proposed the multipole method to solve the two-dimensional multiple scattering by a series of circular cylinders [12]. Eyges extended the multipole method to the multiple scattering of two spheres [13]. The key to the multipole method is the addition theorems of the spherical basis functions which transform the multiple scattering fields of different targets to a field with a uniform origin [14, 15]. Gabrielli derived the analytical solutions to the acoustic scattering problems of two spheres using the symmetry properties, which are widely used in quantum physics [15, 16]. However, there are some limitations to analytical methods in solving the acoustic scattering of non-spherical targets in the spherical coordinate system. Consequently, the semi-analytical method is applied to study the far-field acoustic scattering of a pair of rigid spheroids.

The T-matrix method, a typical semi-analytical method, is first proposed by Waterman for electromagnetic scattering problems [17] and later extended to the acoustic scattering field [18]. The T-matrix method is defined as the semi-analytical method derived from the Helmholtz integral equation and null-field theory, and the infinite matrix needs to be truncated. The crux of the T-matrix method is to expand all the field quantities by a set of orthogonal basis functions and solve the unknown expansion coefficients. In addition, the T-matrix method is suitable for the acoustic scattering problem with arbitrary incidence and scattering angles. Peterson derived the T-matrix expression for multi-target scattering based on the addition theorems of spherical basis functions and calculates the numerical result of a pair of identical spheres under the plane wave incidence [19]. Most of the published literature works focus on the acoustic scattering of a pair of spheres and spheroids with small aspect ratios (i.e., the ratio of the polar radius to equatorial radius is less than 2). However, it is important to study the acoustic scattering mechanism of a pair of rigid oblate spheroids and prolate spheroids, which are extensively used in hydrodynamics and underwater engineering.

In this work, the addition theorems are embedded in the T-matrix method to investigate the underwater far-field acoustic scattering characteristics of a pair of rigid spheroids with different aspect ratios ensonified by plane waves at different angles. The propagation paths of the returning backscattering waves from a pair of rigid spheres are analyzed by using the geometric and numerical method in the time domain. The structure of this work is as follows: In Section 2, acoustic scatterings of a single rigid spheroid are considered by the traditional T-matrix method; in addition, the addition theorems are embedded in the T-matrix method to investigate the scatterings of a pair of rigid spheroids. In Section 3, some numerical experiments are carried out to verify the accuracy and convergence of the T-matrix method for the acoustic scattering of a pair of rigid spheroids. Furthermore, the effects of the separation distance between spheroids, aspect ratios, non-dimensional frequencies, and incidence angles of the plane wave on the acoustic scattering of a pair of rigid spheroids are investigated, while conclusions are provided in Section 4.

2 T-matrix method

In this section, a rigid spheroid and a spheroid pair are considered. The acoustic scatterings of such models under plane wave incidence at an arbitrary angle are investigated using the T-matrix method.

2.1 For a rigid spheroid

As shown in Figure 1, a rigid rotation spheroid with the polar radius b (distance from the center point to the endpoint of the spheroid located in the rotation axis Oz) and equatorial radius a is positioned in an ideal fluid medium with density ρ and sound velocity c. The angle between the incident plane wave P and axis Oz is θi. The surface shape function Sθ of the spheroid is only related to the polar angle θ, such that

Sθ=cos2θ/b2+sin2θ/a21/2.(1)

FIGURE 1
www.frontiersin.org

FIGURE 1. Geometry model for the acoustic scattering of a rigid (A) prolate and (B) oblate spheroid.

The entire wave field can be constructed by the scalar velocity potential because the wave field exists only in the ideal fluid medium. For convenience, in what follows, the monochromatic time factor expiωt of all velocity potentials is omitted. The total wave field Φr is the sum of the incident field Φir and scattering field Φsr:

Φr=Φir+Φsr.(2)

All of the aforementioned velocity potentials satisfy the Helmholtz equation:

2+k2ΦrΦirΦsr=0,(3)

where k=ω/c is the wave number. Waterman discussed the aforementioned three field functions in detail [18]. The outer surface S of the ellipsoid is smooth, so that the Poincare–Huygens principle and Gauss theorem can be applied. The Helmholtz integral representation equation can be expressed as

Φir+SΦ+ngr,rn+Φgr,rdS=ΦrroutsideS0rinsideS,(4)

where Φ+ is the velocity potential on the outer surface of the spheroid. n+Φ denotes the normal velocity on the outer surface of the spheroid. gr,r is the free-field Green’s function, which can be expressed as

gr,r=expikrr/krr.(5)

The crucial point is to expand the whole field quantities with a set of orthogonal basis functions and solve the corresponding unknown coefficients. The scalar spherical basis function is expressed as

Φnmσr=hn1krYnmσθ,φ=hnkrξnmPnmcosθcosmφ,σ=esinmφ,σ=o,(6)

where Pnmcosθ is the associated Legendre function of order n and rank m. The azimuthal parity index σ is either even or odd. hn1kr is a spherical Hankel function of the first kind. The normalization coefficients ξnm can be expressed as

ξnm=εm2n+1nm!1/2×4πn+m!1/2,(7)

with εm=1,m=02,m0. The value range of the order index is n=0,1,,, and the rank index m takes the integer values m=0,1,,n.

The incident and scattered field can be expanded into the form of the weighted sum of the scalar basis function with the expanded coefficients. The regular spherical basis function, denoted by ReΦnmσ, is obtained by substituting the first-kind spherical Bessel function jnkr for hn1kr into Eq. 6. Therefore, the incident and scattered fields can be expanded as

Φir=nmσanmσReΦnmσr,(8)
Φsr=nmσfnmσΦnmσr,(9)

where nmσ=n=0m=0nσ is the triple summation symbol, anmσ denotes the expanded coefficient of the incident field, and fnmσ denotes the unknown expanded coefficient of the scattering field. For a plane harmonic wave, the expanded coefficient is

anmσ=4πξnminPnmcosθicosmφi,σ=esinmφi,σ=o.(10)

Furthermore, the free-field Green’s function gr,r can be expanded using the spherical basis function as

gr,r=nmσik Φnmσr>ReΦnmσr<,(11)

where r>=r,r>rr,r<r and r<=r,r<rr,r>r. r is the modulus of the point r on the surface S of the spheroid. r is the modulus of the field point r. The expanded expression of the unknown surface field Φ+ is

Φ+r=nmσαnmσReΦnmσr,(12)

where αnmσ, which is omitted in the process of deriving the T-matrix, is the expanded coefficient of the unknown surface field.

For a rigid spheroid, the boundary of the spheroid at r=Sθ satisfies the Neumann boundary condition:

n+Φr|r=Sθ=0.(13)

Substituting Eqs 613 into Eq. 4 yields

inmσQnmσ,nmσαnmσ=anmσ,(14)
inmσReQnmσ,nmσαnmσ=fnmσ,(15)

where

Qnmσ,nmσ=kSReΦnmσrnΦnmσrdS.(16)

The detailed expression of Qnmσ,nmσ is obtained by substituting Eq. 6 into Eq. 16:

Qnmσ,nmσ=0πξnmjnkrPnmcosθξnm×hnkrrPnmcosθrθr2hnkr×Pnmcosθθ×r2sinθdθ02πcosmφsinmφcosmφsinmφdφ,(17)

where rθ=rθ/θ and rθ=Sθ. In the expression of the Q-matrix, the symmetry of the scatterer leads to certain elements in the Q-matrix becoming zero. In this study, the rotational symmetry of the spheroid makes the Q-matrix block diagonal of the submatrices, and thus, the computation can be simplified.

From Eqs 14, 15, the relationship between scattering and incident expanded coefficients can be expressed as

fnmσ=Tanmσ,(18)

where the transition matrix T, showing the influence of the geometry of the spheroid, the incident angle, and the frequency of the plane wave [20, 21], is expressed as

T=ReQQ1.(19)

2.2 For a pair of rigid spheroids

In this part, the formula of the T-matrix method for a pair of rigid spheroids immersed in the idea fluid is derived. The geometry of the configuration to be considered is shown in Figure 2. The formula of the T-matrix method for a pair of rigid spheroids is derived exactly like Eq. 4:

Φir+S1+S2Φ+ngr,rn+Φgr,rdS=ΦrroutsideS1andS20rinsideS1orS2.(20)

FIGURE 2
www.frontiersin.org

FIGURE 2. Geometry model for the acoustic scattering of a pair of rigid spheroids.

The incident field Φir and the free-field Green’s function gr,r need to be expanded at origins O1 and O2 to obtain the Q-matrix of each of the spheroids, which are analogous to Eq. 16. When the field point r is located inside S1,

r=r1d,r=r1,rr1=r1r1,rr2=r1r2+2d,(21)

where d=O1O=OO2. According to Eq. 11, Green’s functions gr,r1 and gr,r2 can be written as

gr,r1=nmσik Φnmσr1ReΦnmσr1,(22)
gr,r2=nmσik Φnmσr2+2dReΦnmσr1.(23)

From Eqs 8, 21, where r=r1d, the incident fields can be expressed as

Φir=nmσanmσReΦnmσr1d.(24)

The addition theorems of the spherical basis functions are used in Eqs 23, 24. The translation properties are as follows [22, 23]:

ReΦnmσr+d=nmσR^nmσ,nmσdReΦnmσr,Φnmσr+d=nmσR^nmσ,nmσdΦnmσr,r>d,nmσRnmσ,nmσdReΦnmσrr<d,(25)

where the coordinates of the vector d in the spherical coordinate system are d,η,ψ. R^ is the regular form of R by replacing the spherical Hankel function with spherical Bessel functions. The matrices R^nmσ,nmσd are given as

R^nmσ,nmσd=1m2εmεm1/2[1mCnm,nmd,ηcosmmψ+1σCnm,nmd,ηcosm+mψ],σ=σ,R^nmσ,nmσd=1m2εmεm1/2[1σ+mCnm,nmd,ηsinmmψ+Cnm,nmd,ηsinm+mψ],σσ,(26)

where

Cnm,nmd,η=λ=nnn+n1m+n+n+n+λ/22λ+12n+12n+1λmm!λ+mm!1/2×nnλ000nnλmmmmjλkdPλmmcosη,(27)

where is the Wigner 3-j symbol [24]; its expression is given as follows:

j1j2j3m1m2m3=1j1j2m3δm3,m1+m2j1+j2j3!j2+j3j1!j1+j3j2!j1+j2+j3+1!×Πi=1,2,3ji+mi!jimiv1vv!j1+j2j3v!×j2+m2v!j3j1m2+v!j3j2+m1+v!×j1m1v!1,δm3,m1+m2=1,δm3,m1+m2=0,m3m1+m2.(28)

In this study, the expression of the matrices R^nmσ,nmσ±d for the special case of the offset along the z-axis is

R^nmσ,nmσ±d=λ=nnn+n1m+n+n+n+λ/22λ+12n+12n+11/2×nnλ000nnλmm0±1λjλkd×δσσδmm1,m>01,m=0,σ=e1,m=0,σ=o.(29)

Substituting Eqs 2225 into the second formula of Eq. 20 yields

0=nmσnmσanmσR^nmσ,nmσdReΦnmσr1+nmσikS1ReΦnmσr1Φ+r1×nΦnmσr1n+Φr1Φnmσr1dS+nmσnmσik×Rnmσ,nmσ2dS2ReΦnmσr1Φ+r2nΦnmσr2n+Φr2ReΦnmσ(r2dS.(30)

The expansion coefficients of the unknown surface fields of two spheroids are αnmσ1 and αnmσ2, respectively. The boundary of these two spheroids satisfies the Neumann boundary condition. In this way, Eq. 30 can be expressed as

0=nmσanmσR^nmσ,nmσdinmσQnmσ,nmσ1αnmσ1inmσnmσRnmσ,nmσ2dReQnmσ,nmσ2αnmσ2,(31)

where Qnmσ,nmσ1 and Qnmσ,nmσ2 are given by Eq. 17, with the integrals taken over S1 and S2. Expressing the aforementioned equation in terms of matrices and vectors yields

R^tda=iQ1α1+iR2dReQ2α2,(32)

where R^td=R^d and the superscript “t” denotes the transposition of the matrix. In the same way, when the field point r is located inside S2,

R^tda=iR2dReQ1α1+iQ2α2.(33)

Considering the field point r located outside the sphere SO (the smallest external sphere containing both spheroids S1 and S2 with O as the center of the sphere), according to Eq. 11, the free-field Green’s functions gr,r1 and gr,r2 can be written as

gr,r1=nmσikΦnmσrReΦnmσr1=nmσnmσikΦnmσrR^nmσ,nmσdReΦnmσr1,(34)
gr,r2=nmσikΦnmσrReΦnmσr2=nmσnmσikΦnmσrR^nmσ,nmσdReΦnmσr2.(35)

Substituting Eqs 34, 35, 9 into the first formula of Eq. 20 yields

f=iR^dReQ1α1iR^dReQ2α2.(36)

From Eqs 32, 33, the surface field coefficients α1 and α2 can be solved into expressions in terms of the incident field coefficient a. Afterward, the transition matrix T1,2 can be obtained by substituting the solved α1 and α2 into Eq. 36, such that

T1,2=R^dT1IR2dT2R2dT11×I+R2dT2R^2dR^d+R^dT2IR2dT1R2dT21×I+R2dT1R^2dR^d,(37)

where I is the identity matrix and T1 T2 is the T-matrix formulation of spheroids S1 S2, which can be calculated using Eq. 19.

Since the far-field scattering characteristics are majorly considered in the current study, the form function f can be defined by [25, 26]

fkr0,θ,φ=limr2rr0ΦsΦ0eikr,(38)

where r0=maxa,b. According to Eq. 38, the factors that affect the form function are the observation position of the scattering field, incident angle, and non-dimensional frequency.

3 Numerical examples and results

In this section, the convergence and accuracy of the T-matrix method for calculating the acoustic scattering of a pair of rigid spheres and spheroids are shown by several numerical experiments. Afterward, the monostatic and bistatic acoustic scattering form function modulus f of a pair of rigid spheroids illuminated by the plane wave are calculated by the T-matrix method considering various models with a range of aspect ratios b/a, separation distances 2d, non-dimensional frequencies kr0 (r0=maxa,b), and incidence angles θi.

3.1 Numerical validation

The convergence of the T-matrix method for calculating the acoustic scattering of a pair of rigid spheres under the plane wave end-on incidence (θi=0) is verified. According to the aforementioned theoretical formulation, the infinite summation in the Q-matrix and R-matrix must be truncated at the appropriate limit Nmax. Figure 3 shows the results of f versus Nmax under the end-on incidence case with the non-dimensional frequency kb=5 at separation distances 2d=2b and 4b. It can be seen from Figure 3 that the T-matrix method converges faster in calculating the acoustic scattering problem of a pair of spheres. The relative convergence error of f can be given by [27]

εNmax=fNmaxfNmaxfNmax.(39)

FIGURE 3
www.frontiersin.org

FIGURE 3. Convergence study of the T-matrix method for a pair of rigid spheres. The backscattering modulus f of a pair of rigid spheres was calculated at end-on incidence (θi=0°) with the dimensionless frequency kb=5 versus Nmax for a pair of rigid spheres with separation distances 2d=2b and 4b, respectively.

The truncation factor Nmax in this study is expressed as Nmax=Intk×2d+r0+5. When the dimensionless frequency kb=5, the values of the truncation factor Nmax are taken as 20 and 30 for separation distances of d=b and 2b, respectively. As shown in Figure 3, the form function modulus begins to converge, where Nmax is greater than 14 and 26, respectively, and the relative convergence error could reach 105106. Therefore, the expression of the truncation factor Nmax is verified to be correct and feasible. Furthermore, the convergence and accuracy of the T-matrix method for calculating the acoustic scattering of a pair of rigid spheroids are studied. Figure 4 shows the backscattering modulus f of a pair of rigid spheroids with b/a=4 for the separation distance 2d=4r0 versus Nmax with kb=10 under three incident cases (θi=0°,45°, and 90°). The results show that all the curves have good convergence and the solutions are stable and accurate. Similarly, the relative convergence error could reach 105106.

FIGURE 4
www.frontiersin.org

FIGURE 4. Convergence study of the T-matrix method for a pair of rigid spheroids with aspect ratio b/a=4 for the separation distance 2d=4r0. The backscattering modulus f of a pair of rigid spheroids was calculated with the dimensionless frequency kb=10 versus Nmax under three incident cases (θi=0°,45°, and 90°).

In the following, some numerical experiment results of the acoustic scattering of a pair of rigid spheres calculated by the present method are compared with Peterson’s works [19] and the analytical results. The comparison results are shown in Figure 5. The circles are values according to the data from Peterson’s works, while the solid lines are calculated by the present method. Panel (a) of Figure 5 displays the backscattering cross section of a pair of rigid spheres at the broadside incidence case (θi=90,φi=0) for kb=2 versus the separation distance 2kd. It can be seen from Panel (a) of the figure that the numerical results are in complete agreement with Peterson’s solutions. The scattering cross section is defined by

σNr,θs,φs;θi,φi=4πr2Φs2Φ021πb2.(40)

FIGURE 5
www.frontiersin.org

FIGURE 5. Accuracy study of the T-matrix method for a pair of rigid spheres. (A) Backscattering cross section of a pair of rigid spheres at broadside incidence (θi=90°) for kb=2 versus separation distance 2kd. (B) Backscattering cross section of a pair of rigid spheres for kb=2 and 2kd=4.5 versus the incidence angle θi. (C) Backscattering form function f of a pair of rigid spheres at broadside incidence (θi=90°) versus kb for the separation distance 2d=2b. (D) Backscattering form function f of a pair of rigid spheres at broadside incidence (θi=90°) versus kb for the separation distance 2d=2b.

Panel (b) of Figure 5 displays the far-field backscattering cross section σN versus the incidence angle θi of a pair of rigid spheroids for kb=2 and 2kd=4.5. It is in complete agreement with the results in Figure 7 of Peterson’s works. It can be assured that the T-matrix method is sufficiently accurate in calculating the acoustic scattering of a pair of rigid spheres at arbitrary angular incidence. As shown in panels (c) and (d), the results of the backscattering form function f of a pair of rigid spheres at the broadside incidence case calculated by the present method are exactly the same as those calculated by the analytical method [28] for separation distances 2d=2b and 4b. It can be seen from the aforementioned comparison that the T-matrix method is effective in computing the acoustic scattering of a pair of rigid spheres and spheroids.

3.2 Far-field acoustic scattering properties of a pair of rigid spheres

In this part, the far-field scattering properties of a pair of rigid spheres under the plane wave at arbitrary incident angles are studied by the T-matrix method. Figure 6 displays the backscattering form function f versus the dimensionless frequency ka of a pair of rigid spheres with two values of the separation distance (2d=2b and 4b) under oblique (θi=45°) incidence and end-on (θi=0°) incidence, respectively. As shown in Panels (a) and (b), it can be seen that the form function f (solid line) of a pair of spheres is below the values (dashed line) of twice the form function modulus of a single rigid sphere under plane wave oblique incidence. Furthermore, the backscattering form function curves oscillate faster as the separation distance increases at the oblique incidence case. This phenomenon is similar to the end-on incidence case (see Panels (c) and (d) of Figure 6), but the number of peaks in the backscattering form function curves under the oblique incidence case at the same separation distance is less than the end-on incidence case. In the end-on incidence case, the amplitude of the backscattering form function f tends to be stable as kb increases for the separation distance of 2d=4b (see Panel (d) of Figure 6). The reason for these oscillations in the backscattering form function curves is due to the interference between the specular reflection wave components and Franz wave components in the backscattering response wave [29, 30]. It can be seen that the oscillation period in the backscattering form function curves (solid line) of a pair of rigid spheres is not a constant compared to that (dashed line) of the single rigid sphere [31, 32]. The propagation mechanisms of the specular reflection and Franz wave in the backscattering response wave by a pair of rigid spheres are discussed in the following paragraph.

FIGURE 6
www.frontiersin.org

FIGURE 6. Backscattering form function f (solid lines) of a pair of rigid spheres at oblique incidence (θi=45°) versus kb for the separation distances (A) 2d=2b and (B) 2d=4b. The backscattering form function f of a pair of rigid spheres at end-on incidence (θi=0°) versus kb for separation distances (C) 2d=2b and (D) 2d=4b. The dashed lines represent twice the value of the backscattering form function f of a single rigid sphere.

The path of backscattering waves by a pair of rigid spheres ensonified by the pulse wave from end-on incidence and oblique incidence, respectively, is shown in Panels (a) and (b) of Figure 7. P1 and P3 represent specular reflection waves, which propagate backwards immediately when the incident wave reaches the front of the first sphere and the second sphere, respectively. P2 and P4 represent the backward Franz waves, which diffracted around the shadowed surface of the obstacle, from the first sphere and second sphere, respectively. We define the moments of the arrival of backscattering response wave components P1, P2, P3, and P4 at the observation point as t1, t2, t3, and t4, respectively. The time intervals under end-on incidence can be predicted by the following equations:

t2t1=CD¯+DEF+FG¯c,t3t1=2×CD¯+DH¯/c,t4t3=HI¯+IJK+KL¯/c,(41)

where c denotes the velocity of the specular reflection wave. Likewise, the time intervals under the oblique incidence can be predicted as

t2t1=CD¯+DEF+FG¯/c,t3t1=2×MN¯/c,t4t3=HI¯+IJK+KL¯c.(42)

FIGURE 7
www.frontiersin.org

FIGURE 7. Propagation path of the backscattering waves by a pair of spheres at (A) end-on incidence and (B) oblique incidence.

The ideal fluid medium around the spheres is water. The velocity of the specular reflection wave is 1,510 m/s, and the scattering field point is located at (−30, 0). The time domain response results obtained by the inverse fast Fourier transform (IFFT) of the frequency domain response results [33] of a pair of rigid spheres at three incidence cases (θi=0°,45°,and90°) of pulse waves for the separation distance 2d=4b=4m are shown in Figures 8A–C. Four obvious peaks can be obtained from Figures 8A, B. The specular reflection waves from the first sphere and second sphere lead to two peaks with the largest values. The two peaks with the smaller values are the Franz waves of the first sphere and second sphere, respectively. The values of Franz waves are smaller than those of specular reflection waves due to radiation damping when the Franz wave propagates in the shadowed area of the sphere. It can be seen that there are only two obvious peaks in Panel (c) of Figure 8 because the two spheres do not block each other in the direction of the pulse wave. In the case of broadside incidence (θi=90°), the specular reflection waves of two spheres reach the observing point at the same time, and subsequently, the Franz waves also reach the observing point at another identical time point. Furthermore, the time intervals corresponding to the peak-to-peak intervals for a pair of rigid spheres with three incidence cases (θi=0°,45°,and 90°) of the pulse wave for the separation distance 2d=4b=4m in Figure 8 are listed in Tables 1, 2, 3. The predicted results based on the geometric model are calculated from Eqs 41, 42. As can be seen from Tables 1, 2, 3, the results calculated by the T-matrix method are consistent with the predicted results from the geometric model.

FIGURE 8
www.frontiersin.org

FIGURE 8. Backscattering response waves of a pair of rigid spheres in the time domain for the separation distance 2d=4b=4m under the pulse wave. (A) End-on incidence (θi=0°), (B) oblique incidence (θi=45°), and (C) broadside incidence (θi=90°).

TABLE 1
www.frontiersin.org

TABLE 1. Time intervals regarding the peak-to-peak intervals in Panel (a) of Figure 7.

TABLE 2
www.frontiersin.org

TABLE 2. Time intervals regarding the peak-to-peak intervals in Panel (b) of Figure 7.

TABLE 3
www.frontiersin.org

TABLE 3. Time intervals regarding the peak-to-peak intervals in Panel (c) of Figure 7.

In order to quantitatively study the scattering of a pair of rigid spheres, the bistatic 2D directivity plots of a pair of rigid spheres under three cases (θi=0°,45°,and 90°) of the plane waves for the dimensionless frequency kb=2 and separation distances 2d=2b and 4b, respectively, are shown in Figures 9A–F. It can be found that the scattering form function f (solid line) of a pair of spheres is mostly below the dashed line, which is twice the value of the form function of the single sphere. The maximum values of the scattering form function f in the 2D far-field directivity plots are in the forward scattering direction under three cases. The 2D far-field directivity plots under the oblique incidence are asymmetric in Figures 9C, D because the obstacles are asymmetric along the wave axis in the oblique incidence case. It is clearly shown that the number of petals increased as the separation distance increases at the same incident angles.

FIGURE 9
www.frontiersin.org

FIGURE 9. Bistatic 2D directivity pattern (solid line) of a pair of rigid spheres at end-on (θi=0°) incidence with kb=2 for separation distances (A) 2d=2b and (B) 2d=4b. The bistatic 2D directivity pattern of a pair of rigid spheres at oblique (θi=45°) incidence with kb=2 for separation distances (C) 2d=2b and (D) 2d=4b. The bistatic 2D directivity pattern of a pair of rigid spheres at broadside (θi=90°) incidence with kb=2 for separation distances (E) 2d=2b and (F) 2d=4b. The dashed lines represent twice the value of the 2D directivity pattern of a single rigid sphere under three cases of incident angles (θi=0°, 45°, and 90°) for kb=2.

3.3 Far-field acoustic scattering properties of a pair of rigid spheroids

In this section, the far-field scattering properties of a pair of rigid spheroids under arbitrary incident angles of the plane wave are studied by the T-matrix method. The results of the backscattering form function f (solid lines) of a pair of rigid prolate (b/a=2) spheroids at three incidence cases (θi=0°,45°and 90°) versus the non-dimensional frequency kr0 (r0=maxa,b) for the separation distance 2d=2r0 and 4r0, respectively, are shown in Figures 10A–F. The dashed lines represent twice the value of the backscattering form function f of a single rigid sphere at broadside incidence. As shown in Figures 10A–D, the backscattering form function curves oscillate faster as the separation distance increases in the end-on incidence and oblique incidence cases. This phenomenon is consistent with the results of a pair of rigid sphere scattering in Section 3.2. Those oscillations in backscattering form function curves are caused by the interference between the specular reflection wave components and Franz wave components in the backscattering response wave. The peaks in the backscattering form function curves occur when the wave path difference Δs between the Franz wave and specular reflection wave, which is related to the separation distance 2d, is an integral multiple of the wavelength λ. The peak-to-peak interval Δkr0 in the backscattering form function curves can be expressed as Δkr0=2πaΔs. Obviously, Δkr0 decreases as the separation distance increases, making the backscattering form function curves oscillate faster. In the broadside incidence case, the solid line and dashed line overlap at low frequencies in Panel (e) of Figure 10, implying that there is no interaction between these two spheroids in this frequency band. As the separation distance grows, the solid line and dashed line still overlap at higher frequencies in Panel (f) of Figure 10. This is because the two spheroids do not block each other in the direction of the incident wave, and the interactions between these two spheroids decrease as the separation distance grows.

FIGURE 10
www.frontiersin.org

FIGURE 10. Backscattering form function f (solid lines) of a pair of rigid prolate spheroids under end-on incidence (θi=0°) versus kr0 (r0=maxa,b) for separation distances (A) 2d=2r0 and (B) 2d=4r0. The backscattering form function f (solid lines) of a pair of rigid prolate spheroids under oblique incidence (θi=45°) versus kb for separation distances (C) 2d=2r0 and (D) 2d=4r0. The backscattering form function f of a pair of rigid prolate spheroids under broadside incidence (θi=90°) versus kb for separation distances (E) 2d=2b and (F) 2d=4b. The dashed lines represent twice the value of the backscattering form function f of a single rigid prolate spheroid under three incident cases (θi=0°, 45°, and 90°).

Figure 11 displays the backscattering form function f of a pair of rigid spheroids with different aspect ratios (b/a=1/2, 1, 2, and 4) versus kr0 (r0=maxa,b) for the separation distance 2d=4r0 at three incidence cases (θi=0°,45°,and 90°). As shown in Figure 11A, the values of the backscattering form function f of a pair of rigid spheroids in the end-on incidence case decrease as the aspect ratios grow. This is because the backscattering response wave of a pair of rigid spheroids under plane wave end-on incidence is mainly composed of the specular reflection wave at the front point of the spheroids, and the intensity of the specular reflection wave increases as the radius of curvature at the front point of the spheroids increases [34]. It can be seen from Figure 11 that the number of peaks in the backscattering form function curves for a pair of rigid prolate (b>a) spheroids decreases with the increase in the incidence angle. In order to quantitatively study the scattering of a pair of rigid spheroids under the plane wave, the bistatic 2D directivity plots for a pair of rigid spheroids with different aspect ratios (b/a=1/2, 1, 2, and 4) at three incidence cases (θi=0°,45°,and 90°) of the plane waves for the non-dimensional frequency kb=2 and separation distances 2d=2r0 and 4r0, respectively, are shown in Figures 12A–F. The zoomed-in plots regarding the results of the 2D directivity diagrams of a pair of rigid prolate spheroids (b/a=4) are also given in these panels. It is clearly seen that the number of petals in the 2D directivity plots is mainly influenced by the separation distance, and the number of petals increased as the separation distance increased. The 2D directivity plots are asymmetric in Figures 12C, D because the scatterers are asymmetric along the wave axis in the oblique incidence case, which is different from the end-on incidence and broadside incidence case. In the end-on incidence case (see Figures 12A, B), the 2D directivity plots of a pair of rigid spheroids change from the forward scattering dominant to the sideward scattering dominant, with the increase in the aspect ratio. The maximum values in the 2D far-field directivity plots of a pair of rigid spheroids decrease as the aspect ratio increases in the end-on incidence and oblique incidence cases. In the broadside incidence case (Figures 12E, F), the 2D directivity plots of a pair of oblate spheroids and prolate spheroids are sideward scattering dominant and backward scattering dominant, respectively. It can be found that the 2D far-field directivity plots of a pair of rigid prolate spheroids with different aspect ratios at the same incident angle and separation distance are very similar, and the values of the scattering form function f decreases as the aspect ratio increases.

FIGURE 11
www.frontiersin.org

FIGURE 11. Backscattering form function f of a pair of rigid spheroids with different aspect ratios (b/a=1/2, 1, 2, and 4) versus kr0 for the separation distance 2d=4r0 under the plane wave. (A) End-on incidence (θi=0°), (B) oblique incidence (θi=45°), and (C) broadside incidence (θi=90°).

FIGURE 12
www.frontiersin.org

FIGURE 12. Bistatic 2D directivity pattern of a pair of rigid spheroids with different aspect ratios (b/a=1/2, 1, 2, and 4) at end-on incidence (θi=0°) with kr0=2 for separation distances. (A) 2d=2r0 and (B) 2d=4r0. The bistatic 2D directivity pattern of a pair of rigid spheroids with different aspect ratios (b/a=1/2, 1, 2, and 4) at oblique incidence (θi=45°) with kr0=2 for separation distances (C) 2d=2r0 and (D) 2d=4r0. The bistatic 2D directivity pattern of a pair of rigid spheroids with different aspect ratios (b/a=1/2, 1, 2, and 4) at broadside incidence (θi=90°) with kr0=2 for separation distances (E) 2d=2r0 and (F) 2d=4r0.

4 Conclusion

In this work, the T-matrix method combined with the addition theorems of spherical basis functions is applied to semi-analytically compute the far-field acoustic scattering of a pair of rigid spheroids under the plane wave at an arbitrary incidence angle. It is verified that the T-matrix method can accurately solve the far-field acoustic scattering problem of a pair of rigid spheres with different separation distances under the plane wave of any angle. In addition, some numerical experiments on a pair of rigid (oblate or prolate) spheroids are carried out by the T-matrix method with the following conclusions:

1) The acoustic scattering by a pair of rigid spheroids is more complicated than that of a single rigid spheroid, and the values of the scattering form function of a pair of rigid spheroids are not equal to twice the far-field scattering form function modulus of a single rigid spheroid.

2) The peak-to-peak interval of the backscattering response curve obtained by the IFFT in the time domain is consistent with the geometric prediction results, which makes it possible to estimate the geometrical dimension and separation distance of multiple scatterers from the scattering response wave.

3) The parameters affecting the far-field scattering form function modulus of a pair of rigid spheroids are aspect ratio b/a, separation distance 2d, dimensionless frequency kr0 (r0=maxa,b), and incident angle θi.

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 author.

Author contributions

WL conceptualized this investigation. YY performed the formal analysis and validation. YY performed the data analysis with advice from QG, YZ, YC, and WL. YY wrote the manuscript. All authors contributed to the article and approved the submitted version.

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.

References

1. Seybert AF, Rengarajan TK. The use of CHIEF to obtain unique solutions for acoustic radiation using boundary integral equations. The J Acoust Soc America (1987) 81(5):1299–306. doi:10.1121/1.2024508

CrossRef Full Text | Google Scholar

2. Hunt JT, Knittel MR, Barach D. Finite element approach to acoustic radiation from elastic structures. The J Acoust Soc America (1974) 55(2):269–80. doi:10.1121/1.1914498

CrossRef Full Text | Google Scholar

3. Seybert AF, Soenarko B, Rizzo FJ, Shippy DJ. An advanced computational method for radiation and scattering of acoustic waves in three dimensions. J Acoust Soc America (1985) 77(2):362–8. doi:10.1121/1.391908

CrossRef Full Text | Google Scholar

4. Chai YB, Li W, Gong ZX, Li TY. Hybrid smoothed finite element method for two-dimensional underwater acoustic scattering problems. Ocean Eng (2016) 116:129–41. doi:10.1016/j.oceaneng.2016.02.034

CrossRef Full Text | Google Scholar

5. Li W, Chai YB, Lei M, Li TY. Numerical investigation of the edge-based gradient smoothing technique for exterior Helmholtz equation in two dimensions. Comput Structures (2017) 182:149–64. doi:10.1016/j.compstruc.2016.12.004

CrossRef Full Text | Google Scholar

6. You X, Gui Q, Zhang Q, Chai Y, Li W. Meshfree simulations of acoustic problems by a radial point interpolation method. Ocean Eng (2020) 218:108202. doi:10.1016/j.oceaneng.2020.108202

CrossRef Full Text | Google Scholar

7. Gui Q, Zhang G, Chai Y, Li W. A finite element method with cover functions for underwater acoustic propagation problems. Ocean Eng (2022) 243:110174. doi:10.1016/j.oceaneng.2021.110174

CrossRef Full Text | Google Scholar

8. Rayleigh JWSB. The theory of sound. London: Macmillan &Company (1896).

Google Scholar

9. Faran JJ. Sound scattering by solid cylinders and spheres. The J Acoust Soc America (1951) 23(4):405–18. doi:10.1121/1.1906780

CrossRef Full Text | Google Scholar

10. Morse PMC, Ingard KU. Theoretical acoustics. New Jersey: Princeton university press (1986).

Google Scholar

11. Záviška F. Über die Beugung elektromagnetischer Wellen an parallelen, unendlich langen Kreiszylindern. Annalen der physik (1913) 4(40):1023–56. doi:10.1002/andp.19133450511

CrossRef Full Text | Google Scholar

12. Williams KL, Marston PL. Backscattering from an elastic sphere: Sommerfeld–Watson transformation and experimental confirmation. J Acoust Soc America (1985) 78(3):1093–102. doi:10.1121/1.393028

CrossRef Full Text | Google Scholar

13. Eyges L. Some nonseparable boundary value problems and the many-body problem. Ann Phys (1957) 2(2):101–28. doi:10.1016/0003-4916(57)90037-4

CrossRef Full Text | Google Scholar

14. Sack RA. Three-dimensional addition theorem for arbitrary functions involving expansions in spherical harmonics. J Math Phys (1964) 5:252–9. doi:10.1063/1.1704115

CrossRef Full Text | Google Scholar

15. Gabrielli P, Mercier-Finidori M. Acoustic scattering by two spheres: Multiple scattering and symmetry considerations. J Sound Vibration (2001) 241(3):423–39. doi:10.1006/jsvi.2000.3309

CrossRef Full Text | Google Scholar

16. Gaspard P, Rice SA. Exact quantization of the scattering from a classically chaotic repellor. J Chem Phys (1989) 90(4):2255–62. doi:10.1063/1.456019

CrossRef Full Text | Google Scholar

17. Waterman PC. Matrix formulation of electromagnetic scattering. Proceeding of the IEEE (1965) 53(8):805–12. doi:10.1109/PROC.1965.4058

CrossRef Full Text | Google Scholar

18. Waterman PC. New formulation of acoustic scattering. J Acoust Soc America (1969) 45(6):1417–29. doi:10.1121/1.1911619

CrossRef Full Text | Google Scholar

19. Peterson B, Ström S. Matrix formulation of acoustic scattering from an arbitrary number of scatterers. J Acoust Soc America (1975) 56:771–80. doi:10.1121/1.1903325

CrossRef Full Text | Google Scholar

20. Mishchenko MI. Electromagnetic scattering by particles and particle groups: An introduction. New York: Cambridge University press (2014).

Google Scholar

21. Mishchenko MI, Hovenier JW, Tracis LD. Light scattering by nonspherical particles: Theory, measurements, and applications. Academic Press (1986).

Google Scholar

22. Danos M, Maximon LC. Multipole matrix elements of the translation operator. J Math Phys (1965) 6(5):766–78. doi:10.1063/1.1704333

CrossRef Full Text | Google Scholar

23. Gonis A, Butler WH. Multiple scattering in solids. New York: Springer (2000).

Google Scholar

24. Edmonds AR. Angular momentum in quantum mechanics. Princeton, NJ: Princeton University Press (1974).

Google Scholar

25. Neubauer WG, Vogt RH, Dragonette LR. Acoustic reflection from elastic spheres. I. Steady-state signals. J Acoust Soc America (1974) 55:1123–9. doi:10.1121/1.1914676

CrossRef Full Text | Google Scholar

26. Dragonette LR, Vogt RH, Flax L, Neubauer WG. Acoustic reflection from elastic spheres and rigid spheres and spheroids. II. Transient analysis. J Acoust Soc America (1974) 55:1130–7. doi:10.1121/1.1914677

CrossRef Full Text | Google Scholar

27. Li W, Chai YB, Gong ZX, Marston PL. Analysis of forward scattering of an acoustical zeroth-order Bessel beam from rigid complicated (aspherical) structures. J Quantitative Spectrosc Radiative Transfer (2017) 200:146–62. doi:10.1016/j.jqsrt.2017.06.002

CrossRef Full Text | Google Scholar

28. Gaunaurd GC, Huang H, Strifors HC. Acoustic scattering by a pair of spheres. J Acoust Soc America (1995) 98:495–507. doi:10.1121/1.414447

CrossRef Full Text | Google Scholar

29. Pillai TAK, Varadan VV, Varadan VK. Sound scattering by rigid and elastic infinite elliptical cylinders in water. J Acoust Soc America (1982) 72:1032–7. doi:10.1121/1.388234

CrossRef Full Text | Google Scholar

30. Eastland GC, Marston PL. Enhanced backscattering in water by partially exposed cylinders at free surfaces associated with an acoustic Franz wave. J Acoust Soc America (2014) 135:2489–92. doi:10.1121/1.4870240

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Gong ZX, Li W, Mitri FG, Chai YB, Zhao Y. Arbitrary scattering of an acoustical Bessel beam by a rigid spheroid with large aspect-ratio. J Sound Vibration (2016) 383:233–47. doi:10.1016/j.jsv.2016.08.003

CrossRef Full Text | Google Scholar

32. Gong ZX, Li W, Chai YB, Zhao Y, Mitri FG. T-matrix method for acoustical Bessel beam scattering from a rigid finite cylinder with spheroidal endcaps. Ocean Eng (2017) 129:507–19. doi:10.1016/j.oceaneng.2016.10.043

CrossRef Full Text | Google Scholar

33. Li W, Liu GR. Estimation of radius and thickness of a thin spherical shell in water using the midfrequency enhancement of a short tone burst response. J Acoust Soc America (2005) 118:2147–53. doi:10.1121/1.2040027

CrossRef Full Text | Google Scholar

34. Tang WL. Highlight model of echoes from sonar targets. Chin J Acoust (1994) 13:131–40. doi:10.15949/j.cnki.0217-9776.1994.02.004

CrossRef Full Text | Google Scholar

Keywords: T-matrix method, addition theorems, a pair of rigid spheroids, acoustic scattering, far-field form function

Citation: Yang Y, Gui Q, Zhang Y, Chai Y and Li W (2023) Acoustic scattering of a pair of rigid spheroids based on the T-matrix method. Front. Phys. 11:1170811. doi: 10.3389/fphy.2023.1170811

Received: 21 February 2023; Accepted: 17 May 2023;
Published: 02 June 2023.

Edited by:

Glauber T. Silva, Federal University of Alagoas, Brazil

Reviewed by:

Vladimir Rabinovich, National Polytechnic Institute (IPN), Mexico
Bernhard Johan Hoenders, University of Groningen, Netherlands

Copyright © 2023 Yang, Gui, Zhang, Chai and Li. 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: Wei Li, hustliw@hust.edu.cn

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.