Abstract
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 [1–3]. Recently, various extended methods, like the smooth finite element method and meshfree method [4–7], 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 [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) [9–11]. 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 (distance from the center point to the endpoint of the spheroid located in the rotation axis ) and equatorial radius is positioned in an ideal fluid medium with density and sound velocity . The angle between the incident plane wave P and axis is . The surface shape function of the spheroid is only related to the polar angle , such that
FIGURE 1
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 of all velocity potentials is omitted. The total wave field is the sum of the incident field and scattering field :
All of the aforementioned velocity potentials satisfy the Helmholtz equation:where 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 aswhere is the velocity potential on the outer surface of the spheroid. denotes the normal velocity on the outer surface of the spheroid. is the free-field Green’s function, which can be expressed as
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 aswhere is the associated Legendre function of order and rank . The azimuthal parity index is either even or odd. is a spherical Hankel function of the first kind. The normalization coefficients can be expressed aswith . The value range of the order index is , and the rank index takes the integer values .
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 , is obtained by substituting the first-kind spherical Bessel function for into Eq. 6. Therefore, the incident and scattered fields can be expanded aswhere is the triple summation symbol, denotes the expanded coefficient of the incident field, and denotes the unknown expanded coefficient of the scattering field. For a plane harmonic wave, the expanded coefficient is
Furthermore, the free-field Green’s function can be expanded using the spherical basis function aswhere and . is the modulus of the point on the surface of the spheroid. is the modulus of the field point . The expanded expression of the unknown surface field iswhere , 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 satisfies the Neumann boundary condition:
Substituting Eqs 6–13 into Eq. 4 yieldswhere
The detailed expression of is obtained by substituting Eq. 6 into Eq. 16:where and . 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 aswhere the transition matrix , showing the influence of the geometry of the spheroid, the incident angle, and the frequency of the plane wave [20, 21], is expressed as
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:
FIGURE 2
The incident field and the free-field Green’s function need to be expanded at origins and to obtain the Q-matrix of each of the spheroids, which are analogous to Eq. 16. When the field point is located inside ,
where . According to Eq. 11, Green’s functions and can be written as
From Eqs 8, 21, where , the incident fields can be expressed as
The addition theorems of the spherical basis functions are used in Eqs 23, 24. The translation properties are as follows [22, 23]:where the coordinates of the vector in the spherical coordinate system are . is the regular form of by replacing the spherical Hankel function with spherical Bessel functions. The matrices are given aswherewhere is the Wigner 3-j symbol [24]; its expression is given as follows:
In this study, the expression of the matrices for the special case of the offset along the z-axis is
Substituting Eqs 22–25 into the second formula of Eq. 20 yields
The expansion coefficients of the unknown surface fields of two spheroids are and , respectively. The boundary of these two spheroids satisfies the Neumann boundary condition. In this way, Eq. 30 can be expressed aswhere and are given by Eq. 17, with the integrals taken over and . Expressing the aforementioned equation in terms of matrices and vectors yieldswhere and the superscript “t” denotes the transposition of the matrix. In the same way, when the field point is located inside ,
Considering the field point located outside the sphere (the smallest external sphere containing both spheroids and with as the center of the sphere), according to Eq. 11, the free-field Green’s functions and can be written as
Substituting Eqs 34, 35, 9 into the first formula of Eq. 20 yields
From Eqs 32, 33, the surface field coefficients and can be solved into expressions in terms of the incident field coefficient . Afterward, the transition matrix can be obtained by substituting the solved and into Eq. 36, such thatwhere I is the identity matrix and is the T-matrix formulation of spheroids , which can be calculated using Eq. 19.
Since the far-field scattering characteristics are majorly considered in the current study, the form function can be defined by [25, 26]where . 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 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 , separation distances , non-dimensional frequencies (), and incidence angles .
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 () 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 . Figure 3 shows the results of versus under the end-on incidence case with the non-dimensional frequency at separation distances and . 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 can be given by [27]
FIGURE 3
The truncation factor in this study is expressed as . When the dimensionless frequency , the values of the truncation factor are taken as 20 and 30 for separation distances of and , respectively. As shown in Figure 3, the form function modulus begins to converge, where is greater than 14 and 26, respectively, and the relative convergence error could reach . Therefore, the expression of the truncation factor 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 of a pair of rigid spheroids with for the separation distance versus with under three incident cases ( and ). The results show that all the curves have good convergence and the solutions are stable and accurate. Similarly, the relative convergence error could reach .
FIGURE 4
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 () for versus the separation distance . 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
FIGURE 5
Panel (b) of Figure 5 displays the far-field backscattering cross section versus the incidence angle of a pair of rigid spheroids for and . 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 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 and . 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 versus the dimensionless frequency of a pair of rigid spheres with two values of the separation distance ( and ) under oblique () incidence and end-on () incidence, respectively. As shown in Panels (a) and (b), it can be seen that the form function (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 tends to be stable as increases for the separation distance of (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
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. and represent specular reflection waves, which propagate backwards immediately when the incident wave reaches the front of the first sphere and the second sphere, respectively. and 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 , , , and at the observation point as , , , and , respectively. The time intervals under end-on incidence can be predicted by the following equations:where denotes the velocity of the specular reflection wave. Likewise, the time intervals under the oblique incidence can be predicted as
FIGURE 7
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 () of pulse waves for the separation distance 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 (), 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 () of the pulse wave for the separation distance 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
TABLE 1
| Time interval | Predicted by the geometric model (ms) | Obtained from Figure 7 (ms) |
|---|---|---|
| t2–t1 | 3.41 | 3.29 |
| t3–t1 | 5.30 | 5.37 |
| t4–t3 | 3.41 | 3.30 |
Time intervals regarding the peak-to-peak intervals in Panel (a) of Figure 7.
TABLE 2
| Time interval | Predicted by the geometric model (ms) | Obtained from Figure 7 (ms) |
|---|---|---|
| t2–t1 | 3.41 | 3.32 |
| t3–t1 | 3.75 | 3.74 |
| t4–t3 | 3.41 | 3.33 |
Time intervals regarding the peak-to-peak intervals in Panel (b) of Figure 7.
TABLE 3
| Time interval | Predicted by the geometric model (ms) | Obtained from Figure 7 (ms) |
|---|---|---|
| t2–t1 = t4–t3 | 3.41 | 3.54 |
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 () of the plane waves for the dimensionless frequency and separation distances and , respectively, are shown in Figures 9A–F. It can be found that the scattering form function (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 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
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 (solid lines) of a pair of rigid prolate () spheroids at three incidence cases () versus the non-dimensional frequency () for the separation distance and , respectively, are shown in Figures 10A–F. The dashed lines represent twice the value of the backscattering form function 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 between the Franz wave and specular reflection wave, which is related to the separation distance , is an integral multiple of the wavelength . The peak-to-peak interval in the backscattering form function curves can be expressed as . Obviously, 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
Figure 11 displays the backscattering form function of a pair of rigid spheroids with different aspect ratios (, , , and ) versus () for the separation distance at three incidence cases (). As shown in Figure 11A, the values of the backscattering form function 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 () 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 (, , , and ) at three incidence cases () of the plane waves for the non-dimensional frequency and separation distances and , 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 () 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 decreases as the aspect ratio increases.
FIGURE 11
FIGURE 12
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 , separation distance , dimensionless frequency (), and incident angle .
Statements
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.
SeybertAFRengarajanTK. 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. 10.1121/1.2024508
2.
HuntJTKnittelMRBarachD. Finite element approach to acoustic radiation from elastic structures. The J Acoust Soc America (1974) 55(2):269–80. 10.1121/1.1914498
3.
SeybertAFSoenarkoBRizzoFJShippyDJ. An advanced computational method for radiation and scattering of acoustic waves in three dimensions. J Acoust Soc America (1985) 77(2):362–8. 10.1121/1.391908
4.
ChaiYBLiWGongZXLiTY. Hybrid smoothed finite element method for two-dimensional underwater acoustic scattering problems. Ocean Eng (2016) 116:129–41. 10.1016/j.oceaneng.2016.02.034
5.
LiWChaiYBLeiMLiTY. Numerical investigation of the edge-based gradient smoothing technique for exterior Helmholtz equation in two dimensions. Comput Structures (2017) 182:149–64. 10.1016/j.compstruc.2016.12.004
6.
YouXGuiQZhangQChaiYLiW. Meshfree simulations of acoustic problems by a radial point interpolation method. Ocean Eng (2020) 218:108202. 10.1016/j.oceaneng.2020.108202
7.
GuiQZhangGChaiYLiW. A finite element method with cover functions for underwater acoustic propagation problems. Ocean Eng (2022) 243:110174. 10.1016/j.oceaneng.2021.110174
8.
RayleighJWSB. The theory of sound. London: Macmillan &Company (1896).
9.
FaranJJJr. Sound scattering by solid cylinders and spheres. The J Acoust Soc America (1951) 23(4):405–18. 10.1121/1.1906780
10.
MorsePMCIngardKU. Theoretical acoustics. New Jersey: Princeton university press (1986).
11.
ZáviškaF. Über die Beugung elektromagnetischer Wellen an parallelen, unendlich langen Kreiszylindern. Annalen der physik (1913) 4(40):1023–56. 10.1002/andp.19133450511
12.
WilliamsKLMarstonPL. Backscattering from an elastic sphere: Sommerfeld–Watson transformation and experimental confirmation. J Acoust Soc America (1985) 78(3):1093–102. 10.1121/1.393028
13.
EygesL. Some nonseparable boundary value problems and the many-body problem. Ann Phys (1957) 2(2):101–28. 10.1016/0003-4916(57)90037-4
14.
SackRA. Three-dimensional addition theorem for arbitrary functions involving expansions in spherical harmonics. J Math Phys (1964) 5:252–9. 10.1063/1.1704115
15.
GabrielliPMercier-FinidoriM. Acoustic scattering by two spheres: Multiple scattering and symmetry considerations. J Sound Vibration (2001) 241(3):423–39. 10.1006/jsvi.2000.3309
16.
GaspardPRiceSA. Exact quantization of the scattering from a classically chaotic repellor. J Chem Phys (1989) 90(4):2255–62. 10.1063/1.456019
17.
WatermanPC. Matrix formulation of electromagnetic scattering. Proceeding of the IEEE (1965) 53(8):805–12. 10.1109/PROC.1965.4058
18.
WatermanPC. New formulation of acoustic scattering. J Acoust Soc America (1969) 45(6):1417–29. 10.1121/1.1911619
19.
PetersonBStrömS. Matrix formulation of acoustic scattering from an arbitrary number of scatterers. J Acoust Soc America (1975) 56:771–80. 10.1121/1.1903325
20.
MishchenkoMI. Electromagnetic scattering by particles and particle groups: An introduction. New York: Cambridge University press (2014).
21.
MishchenkoMIHovenierJWTracisLD. Light scattering by nonspherical particles: Theory, measurements, and applications. Academic Press (1986).
22.
DanosMMaximonLC. Multipole matrix elements of the translation operator. J Math Phys (1965) 6(5):766–78. 10.1063/1.1704333
23.
GonisAButlerWH. Multiple scattering in solids. New York: Springer (2000).
24.
EdmondsAR. Angular momentum in quantum mechanics. Princeton, NJ: Princeton University Press (1974).
25.
NeubauerWGVogtRHDragonetteLR. Acoustic reflection from elastic spheres. I. Steady-state signals. J Acoust Soc America (1974) 55:1123–9. 10.1121/1.1914676
26.
DragonetteLRVogtRHFlaxLNeubauerWG. Acoustic reflection from elastic spheres and rigid spheres and spheroids. II. Transient analysis. J Acoust Soc America (1974) 55:1130–7. 10.1121/1.1914677
27.
LiWChaiYBGongZXMarstonPL. 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. 10.1016/j.jqsrt.2017.06.002
28.
GaunaurdGCHuangHStriforsHC. Acoustic scattering by a pair of spheres. J Acoust Soc America (1995) 98:495–507. 10.1121/1.414447
29.
PillaiTAKVaradanVVVaradanVK. Sound scattering by rigid and elastic infinite elliptical cylinders in water. J Acoust Soc America (1982) 72:1032–7. 10.1121/1.388234
30.
EastlandGCMarstonPL. 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. 10.1121/1.4870240
31.
GongZXLiWMitriFGChaiYBZhaoY. Arbitrary scattering of an acoustical Bessel beam by a rigid spheroid with large aspect-ratio. J Sound Vibration (2016) 383:233–47. 10.1016/j.jsv.2016.08.003
32.
GongZXLiWChaiYBZhaoYMitriFG. T-matrix method for acoustical Bessel beam scattering from a rigid finite cylinder with spheroidal endcaps. Ocean Eng (2017) 129:507–19. 10.1016/j.oceaneng.2016.10.043
33.
LiWLiuGR. 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. 10.1121/1.2040027
34.
TangWL. Highlight model of echoes from sonar targets. Chin J Acoust (1994) 13:131–40. 10.15949/j.cnki.0217-9776.1994.02.004
Summary
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
Volume
11 - 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
Updates
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.