Scattering of Stationary and Non-Stationary Sound by Ideal and Elastic Scatterers, Placed near Interface of Media, in Underwater Sound Channel and in Plane Waveguide

Abstract

In the review we obtained a strict solution of problem of the sound diffraction by an elastic spheroidal shell, located near the interface of a liquid medium with an elastic solid medium. Calculation of the scattered sound field for ideal bodies (a spheroid and an elliptical cylinder) is performed. These bodies are placed on the interface between the liquid and the ideal medium; it is shown, that the main role is played not by interaction of scatterers (real and imaginary), but by interference of their scattered fields. The spectrum of the scattered impulse signal for the body in an underwater sound channel is calculated. It is shown, that at large distances the dominante role is played by the spectral characteristic of the channel itself. Based on the method of imaginary sources and imaginary scatterers, the solution of the current study is to solve the diffraction problem of sound pulse signals at ideal (soft) prolate spheroid, which is put in the plane waveguide with the hard elastic bottom. In the work, it is proved that with such a formulation of problems eliminated, there exists possibility of using the method of normal waves because pulses are bundies of energy and can therefore only be distributed to the group velocity which is inherent in just the method of imaginary sources. Calculations made in the article showed that imaginary sources with small numbers exert the effect of total internal reflection, as the result of the reflection coefficient V by the hard elastic bottom which is complex and the real part of V is close to 1.0 which corresponds to V absolutely hard bottom. Found sequences of reflected pulses for the elastic hard bottom and the absolutely hard bottom floor confirmed this approach. In the final part of the article, on the basis of the received results, a solution (the method integral equations) is given, which is a much more complex problem of the diffraction at the elastic non-analytical scatterer, put in the plane waveguide with the hard elastic bottom.

Share and Cite:

Kleshchev, A. (2017) Scattering of Stationary and Non-Stationary Sound by Ideal and Elastic Scatterers, Placed near Interface of Media, in Underwater Sound Channel and in Plane Waveguide. Open Access Library Journal, 4, 1-27. doi: 10.4236/oalib.1104141.

1. Introduction

In the series of problems on the study of the influence of media interfaces on characteristics of sound scattering by bodies, following variants are usually investigated:

a) An interaction of a scatterer with a single interface between media;

b) The scatterer in the field of interfracting modes of the underwater sound channel;

c) The finding of the total scattered field of the system of real and imaginary sources and scatterers of the plane waveguide.

2. Elastic Scatterer near Interface of Media

The interaction of the scatterer with the interface between media is considered in the example of the problem of scattering of sound by an elastic spheroidal body located at the interface between a liquid an elastic medium [1] . Let the elastic gas-filled prolate spheroidal shell be placed near the boundary a liquid― an elastic medium (Figure 1). The axis of rotation of the shell is parallel to the plane boundary. We introduce two systems of spheroidal coordinates ξs, ηs, φs (s = 1, 2), the first of which (s = 1) we associate with the scatterer, the second (s = 2) with the interface plane. The beginning of the Cartesian coordinate system O2 foci of the second spheroidal coordinate system are defined as projections respectively O1 and foci of the first coordinate system on the plane of the boundary Z2Y2, so that inter-focus distance 2h0 is common for both coordinate systems. The interface plane is two coordinate half-planes ( φ 2 = π / 2 and φ 2 = π / 2 ) coordinate system. In order to relate this solution to the solution of the diffraction problem on the elastic spheroidal shell, we simplify the formulation of the problem and assume that the wave vector k of the incident wave is in the plane X1Z1 (and correspondly X2Z2), φ o s = 0 (see Figure 1).

Now along with the potential Φ 1 ( 1 ) of the wave, scattered by the shell, will appear the potential Φ 1 ( 2 ) frothe elastic half-space:

Φ 1 ( 2 ) = 2 m = 0 n m K m , n S ¯ m , n ( C 1 , η 2 ) R m , n ( 1 ) ( C 1 , ξ 2 ) cos m φ 2 . (1)

Figure 1. Elastic spheroidal shell near the liquid-elastic isotropic medium interface.

The potential Φ 1 ( 2 ) decomposes by radial functions of the first kind, this is due to the fact, that foci of the second coordinate system lie in the plane of the interface, physically this means that the interacting of scatterers are distorted by fields of two plane waves: an incident and reflected by an interface. Expansions of potentials of the shell are accompanied by expansions of potentials of an elastic half-space:

U 2 ( 2 ) = 2 m = 1 n m M m , n S ¯ m , n ( C t 2 , η 2 ) R m , n ( 1 ) ( C t 2 , ξ 2 ) sin m φ 2 ; (2)

W 2 ( 2 ) = 2 m = 0 n m N m , n S ¯ m , n ( C t 2 , η 2 ) R m , n ( 1 ) ( C t 2 , ξ 2 ) cos m φ 2 , (3)

Φ 2 ( 2 ) = 2 m = 0 n m L m , n S ¯ m , n ( C l 2 , η 2 ) R m , n ( 1 ) ( C l 2 , ξ 2 ) cos m φ 2 ; (4)

where C l 2 and C t 2 ―are wave dimensions of longitudinal transverse waves respectively.

The potential of the incident wave Φ0 in two coordinate systems has a form [2] [3] :

Φ 0 = 2 m = 0 n m i n ε m S ¯ m , n ( C 1 , η 0 ) S ¯ m , n ( C 1 , η s ) R m , n ( 1 ) ( C 1 , ξ s ) cos m φ s , ( s = 1 , 2 ) . (5)

The potential of a diffracted field Φ Σ = Φ 0 + Φ 1 ( 1 ) + Φ 2 ( 2 ) obeys simultaneously to boundary conditions on a surface of a shell and on a planar interface between a liquid and an elastic medium. To boundary conditions on a surface of a shell are supplemented by conditions on an inter-face between a liquid-elastic medium:

λ 0 k 2 ( Φ 0 + Φ 1 ( 1 ) + Φ 1 ( 2 ) ) = λ 2 k l 2 2 Φ 2 ( 2 ) + 2 μ 2 u φ φ ( 2 ) | φ = π / 2 ; π / 2 ; (6)

( h φ ( 2 ) / h ξ ( 2 ) ) ( / ξ ) ( u φ ( 2 ) / h φ ( 2 ) ) + ( h ξ ( 2 ) / h φ ( 2 ) ) ( / φ ) ( u ξ ( 2 ) / h ξ ( 2 ) ) | φ = π / 2 ; π / 2 = 0 ; (7)

( h φ ( 2 ) / h ξ ( 2 ) ) ( / ξ ) ( u φ ( 2 ) / h φ ( 2 ) ) + ( h ξ ( 2 ) / h φ ( 2 ) ) ( / φ ) ( u ξ ( 2 ) / h ξ ( 2 ) ) | φ = π / 2 ; π / 2 = 0 ; (8)

( h φ ) 1 ( / φ ) ( Φ 0 + Φ 1 ( 1 ) + Φ 1 ( 2 ) ) = ( h φ ( 2 ) ) 1 ( Φ 2 / φ ) + ( h ξ ( 2 ) h η ( 2 ) ) 1 [ ( / ξ ) ( h η ( 2 ) ψ η ( 2 ) ) ( / η ) ( h ξ ( 2 ) ψ ξ ( 2 ) ) ] | φ = π / 2 ; π / 2 , (9)

where λ2 and μ2―Lame coefficients of an elastic half-space; k l ( 2 ) ―a wavt number of a longitudinal wave of an elastic half-space.

In our formulation of a problem because of a panty of a solution respect a plane XZ boundary conditions for φ = π / 2 completely repeat conditions for φ = + π / 2 , that do not provide any additional information in this case. When substituting potential expansions into boundary conditions for a shell and a planar interface when an expansion coefficients is used an addition theorem for wave spheroidal functions [3] ,

R p , q ( 1 ) , ( 3 ) ( C j , ξ j ) S ¯ p , q ( C j , η j ) exp ( i p φ j ) = n = 0 m = n n R m , n ( 1 ) ( C s , ξ s ) × Q ¯ m , n , p , q ( 1 ) , ( 3 ) ( C j , C s ; l ; θ j s ) S ¯ m , n ( C s , η s ) exp ( i p φ s ) , (10)

where

Q ¯ m , n , p , q ( 1 ) , ( 3 ) = 2 ( i ) n q r = 0 , 1 ' t = 0 , 1 ' d r p q ( C j ) d t m n ( C s ) σ = | r + p t m | r + p + t + m ( i ) σ × b ˜ σ ( r + p , p , t + m , m ) Z σ ( 1 ) , ( 3 ) ( k l ) P ˜ σ p m ( θ j s ) exp [ i ( p m ) φ j s ] ;

θjs―polar angle of a point Os―a beginning of a s-th coordinate system in an i-th system (Figure 2); l―a distance between Oj and Os; d r p q ( C j ) and d t m n ( C s ) ―сoefficients of expansions f functions S ¯ p , q ( C j , η j ) and S ¯ m , n ( C s , η s ) by functions P ˜ r p ( η j ) and P ˜ t m ( η s ) , which up to a constant factor coincide with normalized adjoint Legendre functions; Z σ ( 1 ) = j σ ( k l ) ―spherical Bessel functions; Z σ ( 3 ) = h σ ( 1 ) ( k l ) ―spherical Hanrtl functions of a first kind; coefficients b ˜ σ ( r + p , p , t + m , m ) are obtained from coefficients b σ ( r + p , p , t + m , m ) [3] taking into account a relationship P ˜ t m ( η s ) and associated Legendre functions; a prime for signs Σ means that a summation is carried out on an ever basis r and t, if, respectively, g - p and n - m is even and r and t is odd, if g - p and n - m is odd.

A strict solution can be obtained for another orientation of a spheroidal shell with respect to a plane bourder, namely, under a condition of perpendicularity of a rotation axis of a shell to a plane interface between media (Figure 3). We will consider this orientation in more detail, by replacing a spheroidal shell with an ideal prolate soft spheroid and an elastic half-space by an ideal medium (hard or soft) [4] . We will map a scatterer and a source mirror-wise with respect to a boundary and reduce a problem to diffraction of fields of two sources (real and imaginary) by two spheroidal scatterers (real and imaginary). Potentials Φs (S = 1, 2) of waves scattered by spheroids are chosen in a form of expansions (taking into account an axial symmetry) [4] :

Figure 2. Systems of spheroidal coordinates associated interface with two interacting scatterers.

Figure 3. Ideal spheroid at an face.

Φ s ( ξ s , η s ) = n = 0 B 0 , n s S ¯ 0 , n ( C s , η s ) R 0 , n ( 3 ) ( C s , ξ s ) . (11)

Since spheroids (real and imaginary) are ideally soft, then on their surfaces ( ξ 0 = ξ 01 = ξ 02 ) is satisfied a homogeneous Dirichlet:

Φ 0 + s = 1 2 Φ s = 0 | ξ = ξ 0 ; s = 1 , 2 . (12)

A potential of a falling plane wave is given by a decomposition:

Φ 0 ( ξ s , η s ) = 2 n = 0 i n S ¯ 0 , n ( C s , η s ) R 0 , n ( 1 ) ( C s , ξ s ) S ¯ 0 , n ( C s , 1 ) , s = 1 , 2. (13)

Unknown coefficients B 0 , n s of expansions (11) are sought from an infinite system of equations―boundary conditions (12) [4] :

B 0 , n s + q = 0 B 0 , q t R 0 , n ( 1 ) ( C s , ξ 0 s ) [ R 0 , n ( 3 ) ( C s , ξ 0 s ) ] 1 Q ¯ 0 n 0 q ( 3 ) ( C t , C s ; l ; θ t s ) = 2 i n S ¯ 0 n ( C s , 1 ) R 0 , n ( 1 ) ( C s , ξ 0 s ) [ R 0 , n ( 3 ) ( C s , ξ 0 s ) ] 1 , s = 1 , 2 ; t = 1 , 2 ; s t , (14)

where l―a distance between centers of coordinate systems O1 and O2 (Figure 3), in our case:

θ 12 = 0 , θ 21 = π

Q ¯ 0 n 0 q ( 3 ) ( C t , C s ; l ; θ t s ) = 2 i n + q r = 0 , 1 ' j = 0 , 1 ' d r 0 q ( C t ) d j 0 n ( C s ) σ = | r j | r + j i σ b σ ( r , 0 , j , 0 ) h 0 ( 1 ) ( k l ) P σ ( cos θ t s ) ;

b σ ( r , 0 , j , 0 ) = ( r j 00 | σ O ) 2 .

For a regularization of a system (14) with respect to unknown coefficients B 0 , n s we introduce new unknown X 0 , n s from a ratio [4] :

B 0 , n s = R 0 , n ( 1 ) ( C s , ξ 0 s ) X 0 , n s . (15)

As a result, an infinite system (14) for unknown B 0 , n s is reduced to another infinite system of relatively new unknown X 0 , n s [4] :

X 0 , n s + q = 0 X 0 , q t R 0 , q ( 1 ) ( C t , ξ 0 t ) [ R 0 , n ( 3 ) ( C s , ξ 0 s ) ] 1 Q ¯ 0 n 0 q ( 3 ) ( C t , C s ; l ; θ t s ) = 2 i n S ¯ 0 , n ( C s , 1 ) [ R 0 , n ( 3 ) ( C s , ξ 0 s ) ] 1 . (16)

Further we find by a trunkaction method a solution of a regular system (16). Initially we calculated angular scattering functions D s ( θ s ) of two interacting spheroids distorting a field of a plane monochromatic wave. A Figure 4 shows

Figure 4. Modules of angular characteristics of single and interacting spheroids.

moduluses of angular characteristics | D s ( θ s ) | of tnteracting spheroids (a curve 1 refers to a first spheroid, a curve 2―to a second spheroid). A curve 3 depicts in another scale | D ( θ ) | of a single soft spheroid in an infinite medium. A scale had to be changed so that a curve 3 did not merge with curves 1 and 2. A curve 4 characterizes a modulus | D Σ ( θ 1 ) | of a total angular characteristic in coordinates of a first spheroid (Figure 1):

| D Σ ( θ 1 ) | = | D 1 ( θ 1 ) + D 2 ( θ 1 ) exp ( i k l cos θ 1 ) | .

Calculations were carried out at: C 1 = C 2 = 10.0 , ξ 01 = ξ 02 = ξ 0 = 1.005 , a semi-focus distance h 01 = h 02 = 1 м , l = 8 h 0 . An analysis of curves prestnted in Figure 4 shows that when selected parameters (C1, C2, l) an interaction of scatterers turned out to be smail, because of this curves 1, 2, 3 are so close to each other. A main role is played by interference effects (especially in a shadow region), so a curve 4 stands out sharply (again in a shadow region) agatnst a background of other curves. In a second stage (based on a calculation of a scattered field of two spheroids) were calculated angular characteristics D Σ ( θ ) of a soft spheroid ( ξ 0 = 1.005 ; C = 10 ), located at a distance l 0 = 4 h 0 = 4 m from an interface between liquid and an ideal medium. Results of calculations | D Σ ( θ ) | are shown in Figure 5. A curve 1 corresponds to an interface between liquid and soft media. A certain contribution to D Σ ( θ ) is made by an imaginary source and associated with it interference fringes of both spheroids, physically intense scattering in an illuminated region anses because of a reflection of waves scattered by spheroids by a plane interface between media. We must not forget that I this case we are not talking about a distortion of a plane traveling wave (as it was an unlimited medium), about a distortion of a standing wave field real and imaginary sources.

A strict solution also has a problem of scattering of sound by a spheroidal half-body placed on an interface between a liquid and an ideal medium. With a mirror image of a scatterer and a source with respect to a plane boundary, we obtain a spheroidal scatterer located in a field of two sources (real and imaginary).

Figure 5. Modules of angular characteristics of soft spheroid at an interface of media.

A phase of a wave potential from an imaginary source on a plane interface coincides with a phase of a potential of an incident wave in this plane in a liquid borders on an ideally rigid medium and phases of these waves differ at a boundary by 180˚ if a liquid borders on an ideally soft medium. In Figure 6 shows modules of angular characteristics | D ( θ ) | (in .different scales) of spheroidal half-bodies located on a boundary of a liquid with an ideal medium. A curve 1 corresponds to half a hollow steel oblate spheroidal shell placed on a boundary of a liquid an ideally soft medium. Outher radial coordinate of a shell ξ 0 = 0.1005 , internal ξ 1 = 0.07669 , a wave size C = 7.1 , an irradiation angle θ 0 = 0 ˚ (an axisymmetric problem). Two other curves (see Figure 6) relate to modulus of angular characteristics of a soft oblate hemispheroid bordering on soft (a curve 2) and hard (a curve 3) media at a same irradiation angle θ 0 = 0 ˚ and a wave size of a body C = 10.0 , a radial coordinate of a soft oblate spheroid ξ 0 = 0.1005 .

In Figure 7 show relative cross sections for backscattering σ0 of oblate ideal hemispheroids located on interface between media, where n under D ( θ ; φ ) we mean a total angular characteristic (from real and iimaginary sources). Curves of Figure 7 correspond to axially radiation ( θ 0 = 0 ˚ ) of an oblate hemispheroid (soft or hard) with a radial coordinate ξ 0 = 0.1005 . Curves 1 and 3 refer to a hard hemispheroid placed at a boundary with a hard 1 and soft 3 media (a liquid we have designated by a number I, soft―by a number II, hard―a number III). Curves 3 and 4 correspond to a soft hemispheroid located at a boundary with soft 2 and hard 4 media. Curves 5 and 6 give an idea of a variation σ0 of hard 5 and soft 6 oblate spheroids ( ξ 0 = 0.1005 ) in an infinite liquid medium. It is easy to see that if a hemispheroid and an ideal semibounded medium and same (hard

Figure 6. Modules of angular characteristics of spheroidal semi-bodies located an interface between media.

Figure 7. Relative cross sections for back-scattering of oblate spheroids.

or soft, curves 1 and 2) that is, a half body represents only a violation of a shape of a boundary hides this unevenness and σ0 grows slowly and at smail wave dimensions close to zero and vise versa, for different materials of a hemispheroid and a semibouded ideal medium cross-sections σ0 are much larger (curves 3 and 4) in an entre investigated range of wave dimensions. While acurves 5 and 6 tend asymptotically to a value of geometric acoustics, remaining curves increase indefinitely. Mathematical and physical explanations for this phenomenon were given in comments to Figure 5.

Let us turn to passive characteristics (scattering indicatrix) of hemispheroids located on an interface between a liquid and an ideal medium. We imagine for this that on a fairly large distance from a boundary along a line LM a combined system (source-receiver) moves, while we are interested in a reflected signal to a point finding a combined system. A movement of a system is so slow that a Doppler effect can be ignored. Two orientations of half-bodies are possible which admit are depicted in Figure 8 and Figure 9. In Figure 10 (in different scales) presented modules | D Σ 0 ( θ ) | for a prolate soft hemispheroid on a boundary of a liquid with an ideally hard medium (a curve 1, its scale is shown to a left of a vertical axis, ξ 0 = 1.005 , С = 10 ) and for an oblate hard hemispheroid on a boundary of a liquid with an ideally soft medium (a curve 2, its scale is shown to a right of a vertical axis, ξ 0 = 0.1005 ; С = 10 ). A orientation of hemibodies corresponds to Figure 8, directions of rays―by a oblate coordinate system. In Figure 11 shows modules | D Σ 0 ( θ ) | for a prolate hard hemispheroid on a boundary of a liquid with a ideally soft medium (a curve 1, ξ 0 = 1.005 , С = 10 ) and for a oblate hard hemispheroid on a boundary of a liquid with a

Figure 8. The multual orientation of a combined system and hemispheroid (1-st variant).

Figure 9. The mutual orientation of a combined system and hemispheroid (2-nd variant).

ideally sofy medium (a curve 2, ξ 0 = 0.1005 ; С = 10 ). A orientation of hemibodies corresponds to Figure 9, directions of rays―by a prolate coordinate system. An analysis presented scattering indicatrix shows an advantage of a second type orientation (see Figure 9), since in this case appears an intense non-mirror component, associated with sound reflection of an interface. At orientation of 1-st type (see Figure 8), intensive back reflection at θ = 0 ˚ will be masked by reflection from a planar interface.

We will pass from stationary (harmonic) irradiation to nonstationary radiation in a form of sound pulses with a rectangular envelope and monochromatic

Figure 10. Modules of scattering indicatrics of soft hemispheroids.

Figure 11. Modules of scattering indicatrics of hard hemispheroids.

filling. As before, we will consider interfaces of media of three types: a liquid―an elastic bottom; a liquid―an ideally soft medium (Dirichlet condition); a liquid―an ideally hard medium (Neumann condition).

If a scatterer (oblate hard spheroid) is placed at an interface between media (a liquid―an elastic bottom) at an observation point, come first an impulse of a mirror reflection from a scatterer.

Figure 12, a shows a mirror reflection pulse r 1 Ψ S ( t ) for a hard oblate spheroid when it is irradiated at an angle θ 1 = 30 ; a normalized modulus of a spectrum of a pulse r 1 Ψ S ( t ) is seen from Figure 12(b). After a while an observation point P will reach an pulse Ψ Σ ( t ) reflected from an elastic bottom and diffracted on a spheroid.

Figure 13, a shows a pulse Ψ Σ ( t ) , a normalized modulus spectrum | S Σ ( ν ) | of a pulse Ψ Σ ( t ) , presented on a Figure 13(b).

We orient a hard prolate hemispheroid in such a way that it’s a major semiaxis will be in a plane of an interface between media. We calculate mirror-reflected pulses Ψ Σ ( t ) at an angle of incidence θ 1 = 60 for two variants: 1) a hard prolate hemispheroid on a boundary with a hard medium; 2) a hard prolate

(a)(b)

Figure 12. A mirror reflection pulse r 1 Ψ S ( t ) (а); a normalized modulus of a spectrrum | S S ( ν ) | of a pulse r 1 Ψ S ( t ) (b).

(a)(b)

Figure 13. A diffracted pulse Ψ Σ ( t ) (а); a normalized modulus of a spectrum . | S Σ ( ν ) | of a pulse Ψ Σ ( t ) (b).

hemispheroid on a boundary with a soft medium, In Figure 14 is presented a pulse Ψ Σ ( t ) (a) and a normalized modulus of a spectrum | S Σ ( ν ) | for a first variant and in Figure 15 show same characteristics for a second variant

We see in Figure 16 a diffracted pulse Ψ Σ ( t ) and a modulus of its spectrum | S Σ ( ν ) | in a direction of a mirror component for a hard oblate hemispheroid, located on a border of a liquid―an ideally soft mediumin such a way that its major axis lies in a plane of boundary. An irradiation angle remains a same. θ 1 = 30 .

In Figure 17 shows a diffracted pulse Ψ Σ ( t ) in a direction of a mirror component and a modulus of its spektrum | S Σ ( ν ) | for a soft prolate hemispheroid placed on a border of a liquid―an ideally hard medium. Dimensions, orientation of a hemispheroid, a ratio of semi-axis in a same as in Figure 16. We note that for oblate and prolate hemispheroids located on a interface between media a diffracted pulse Ψ Σ ( t ) in a direction of a mirror component is obtained as a result of an interference of two reflected signals (from a foor of a hemispheroid and a boundary itself), herefore it is designated as a diffracted pulse/ Ψ Σ ( t ) .

(a)(b)

Figure 14. A mirror of reflection pulse Ψ Σ ( t ) for a hard prolate hemispheroid on a boundary with a hard medium (a); a normalized modulus of aspectrum | S Σ ( ν ) | (b).

(a)(b)

Figure 15. A mirror reflection pulse Ψ Σ ( t ) for a hard prolate hemispheroid on a boundary with a spft medium (a); a normalized modulus of a spectrum | S Σ ( ν ) | (b).

(a)(b)

Figure 16. A mirror reflection pulse Ψ Σ ( t ) for a hard oblate hemispheroid (a); a normalized modulus of a spectrum | S Σ ( ν ) | of a pulse Ψ Σ ( t ) (b).

(a)(b)

Figure 17. A mirror reflection pulse Ψ Σ ( t ) for a soft prolate hemispheroid (a); a normalized modulus of a spectrum | S Σ ( ν ) | of a pulse Ψ Σ ( t ) (b).

A picture of a reflection of sound by a spheroidal body located at a boundary of a liquid―an elastic bottom, we supplement it with angular scattering characteristics R ( θ , φ ) for a stationary sound signal with a fixed frequency. In Figure 18 shows a modulus | R 1 ( θ , φ ) | of an angulare characteristic upon irradiation of a soft prolate spheroid with a plane wave at an angle θ 1 = 60 and a modulus | R Σ ( θ , φ ) | of an angulare characteristic upon irradiation with a plane wave reflected from a bottom.

(a)(b)

Figure 18. Modules of angular characteristics of a sound scattering by an ideal spheroid in free medium (a) and in a presence of a liquid―an elastic bottom boundary.

3. Ideal Spheroid in An Underwater Sound Channel

A study of an effect of two boundaries on a field of a spheroidal scatterer will begin withan ideal spheroid placed in an underwater sound channel with non-reflecting boundaries [5] and irradiated with a pulsed signal with monochromatic filling [6] . At a depth z0 of such a waveguide we place a point source of a rulsed signal and at a horizontal distance r from it and at a depth z2―a spheroidal scatterer (Figure 19(a)). A velocity profile of sound in a symmetrical under-water sound channel is shown in Figure 19(b). A receiver of a scattered signal is compatible with a source, that is we will consider a combined system and find a spectrum of a scattered signal at this point. Since by a condition a velocity of sound depends only on a vertical coordinate z, in a horizontal plane rays are not bent and in each vertical half-plane, passing through a center of a scatterer a sound field will be independent of adjacent half-planes, and each half-plane there will be complex coefficients of excitation of a sound channel modes. A result, in each of these half-planes one observes an interference pattern of modes independent of neighboring half-planes. We are interested in a field in one vertical half-plane P containing a coincident point of a source-receiver and a center of a scatterer. Since a scattered field in this half-plane does not depend on a behavior of this field in all other half-planes, we shall take it to be a same in all vertical yalf-planes and equal to a field in half-plane P. Approximately (without taking into account an effect of a medium on an angular scattering characteristic) a spectrum of a scattered signal at a location of a source will be equal to [6] :

S 2 ( ω ) = ρ 0 2 ( r 1 ) 1 m = 1 M n = 1 N m P m P n , m D n m ( ω ) exp [ i ( κ m r 1 + κ n m r 1 π / 2 ) ] , (17)

where: P m = p m ( ρ 0 ) 1 φ m ( z 2 ) φ m ( z 0 ) ; pm―a mode excitation coefficient m; φ m ( z 2 ) ―a intrinsic wave-guide function determined by boundary conditions;

Figure 19. A mutual arrangement of a spheroidal scatterer and a source in a sound channel (a) and a profile of a velocity of sound in a sound channel with non-reflecting boundaries (b).

ρ0―density at source and receiver depth; P n m = ( 1 / ρ 2 ) φ n ( z 2 ) φ n ( z 0 ) ; ρ2― density at a depth of a center of a scatterer; D n m ( ω ) ―space-transfer function of body scattering for a m-th mode of a source and a n-th mode of a scatterer; κm and κnm―horizontal components of wave numbers of modes of an incident and scattered waves res-pectively; M―a largest admissible source mode; N―a largest admissible scatterer mode for m-th mode of a source.

A dependence of a velocity of sound on a coordinate z for a symmetric waveguide (see Figure 19(b)) has a form [5] :

c ( z ) = ( p | z | + q ) 1 / 2 .

We will find a spectrum S 2 ( ω ) at a combined point of a source and receiver for a ideal soft scatterer in a forn of a prolate spheroid with a coordinate of an outher surface ξ 0 = 1.005 . We place a source (receiver) on an axis of a symmetric waveguide ( z 2 = z 0 = 0 ), a major axis of a spheroid being perpendicular to an axis Z. An interfocus spheroid distance 2h0 is assumed equal to 9.7 m. A source generates a pulse signal with a duration τ 0 = 0.05 s at a frequency f 0 = 400 Hz ( С = 8.0 ). Space-transfer function D n m ( ω ) is determined by a frequency and angular (in a plane XOZ) characteristics of a spheroid scattering. With a chosen velocity profile of a sound (see Figure 19(b)) a largest angle between a wave vectors in an incident and scattered waves is approximately 16˚. Referring to an amplitude-phase angular characteristic of a scattering D ( θ ; φ ) of a soft spheroid in a plane XOZ ( θ = 90 , η = 0 ) (see Figure 20) we note that

Figure 20. Amplitude-phase characteristic of sound scattering by a soft spheroid.

even with a maximum wave size ( C = 10.0 ) used in calculations D ( 90 ; φ ) = | D ( 90 ; φ ) | exp [ i ψ ( 90 ; φ ) ] (an angular characteristic) is practically non-directional within angles ϕ = 0 ÷ 16 (see Figure 20). Such an approach, in which modes of an underwater sound channel are considered in a form of plane waves irradiating our spheroidal scatterer is also approximate. Curve 1 in Figure 20 refers to б, | D ( 90 ; φ ) | , curve 2 characterizes a phase ψ ( 90 ; φ ) increased for all angles φ by π. At lower wave dimensions characteristics will be even closer to circular. Therefore, taking into account a first assumption about uniform scattering in a horizontal plane, we can assume that each incident wave from a set of admissible modes will scatter a body uniformly in all directions (as non-directional scatterer) with an angular constant excitation coefficient D ( ω ) into which a shace-transfer function D n m ( ω ) . As a result, a spectrum of a scattered signal will be calculated by a formula [6] :

S 2 ( ω ) = D ( ω ) ρ 0 2 ( r 1 ) 1 m = 1 M n = 1 M P m P n m exp [ i ( κ m r + κ n m r π / 2 ) ] . (18)

4. Diffraction Pulse Sound Signal on the Soft Prolate Spheroid Located in Plane Waveguide with the Hard Elastic Bottom

We turn to a familiar problem of the diffraction of pulses on spheroial bodies I the plane waveguide [7] [8] [9] [10] , retairing the upper boundary condition Dirichlet, waveguide dimensions and scatterer with respect to boundaries, replasing only ideal hard lower boundary on the elastic isotropic bottom. Physical parameters of the lower medium will correspond to фthe isotropic elastic bottom, but in their values, they will be very close to parameters of transversely-isotropic rock―a large gray siltstone [11] . The longitudinal wave velocity in this material is 4750 m/s, the transverse wave velocity―2811 m/s. When used in this case the method of imaginary sources need to enter the reflection coefficient V for the each source [10] , when displaying sources relative to the upper border sources, as before [7] [11] [12] [13] [14] [15] , will change the sign on the opposite, which corresponds to a change of phase by π.

It is known to [10] , that the imaginary sources method boundary conditions are not fulfilled strictly on any of borders of the waveguide even in the case of ideal boundary conditions of Dirichlet and Neumann. For the better fulfillment of these conditions in diffraction problems [7] [11] [12] [13] [14] [15] [16] were introduced imaginary scatterers by mirroring their relatively waveguide boundaries. Likewise conduct imaginary scatterers and in our problem and compare the sequence of reflected pulses [12] [13] [16] in the case of ideal borders and in presence of a hard elastic bottom in the waveguide. In [10] shows that the method of imaginary sources applicable in the case where the reflection coefficient V will be a function of the angle of the incidence of the wave from a source relative to the normal to the boundary. In our case this angle will be determined by the mutual position of the source (real or imaginary) and the scatterers (real or imaginary), which falls the wave from the source. Since the receiver is combined with a real source Q, the sequence of reflected pulses will be determined by the quantity and amplitudes of reflected signals (from different scatterers) having the same propagation time from sources to scatterers and from scatterers to the point Q. Parameters of the waveguide, the position of the real source Q (combined with the receiver) and the real scatterer remained unchanged compared [12] [13] [16] : L = 1000 m, H = 400 m, the real source Q and real scatterer are located at the depth of 200 m., the scatterer is the ideal soft prolate spheroid with the semi-axes ratio a/b = 10 (a = 0.279 m) and its axis of a rotation is perpendicular to the plane of the figure (see Figure 21). The formula for the reflection coefficient V 0 N , where N―the number of a source, is given in [10] . For the calculatio9nof first five of reflected pulses us needed following reflection coefficients: V 03 in the direction of the first (real) scatterer 01, V 05 in the direction of the second (imaginary) scatterer 02, V 06 in the direction of same the second scatterer.

Figure 21. The mutual disposition of the pulse point-soirces and scatterers in the plane waveguide.

As a result of simple calculations with the help of [7] obtain: V 03 = 0.9989 + i 0.0633 ; V 05 = 0.9989 + i 0.0633 ; V 06 = 0.6238 + i 0.7897 . All three coefficients have turned complex, which mesns the total internal reflection at the boundary liquid-hard elastic bottom, therefore all three modules of reflection coefficients are equal 1.0 and real parts of first two coefficients are close to +1.0, which is typical for the boundary liquid-absolutely hard bottom. The resulting sequence of calculations of first five reflected pulses is shown in Figure 22. We compare them to the sequence in Figure 23 for ideal boundaries [12] [13] [16] : 1st and 4th pulsews of Figure 22 are identical with first and second pulses of Figure 23, as for for 2nd, 3rd and 5th pulses in Figure 22 in the case of ideal boundaries and symmetrical location of real a source and a scatterer relatively of boundaries of the waveguide, they are compensated each reflected pulses, i.e. 2nd, 3rd and 5th pulses (see Figure 22) show the difference in sequences of reflected pulses when replacing an absolutely hard bottom on an elastic hard bottom.

A similar pattern is observed for anisotropic bottom, such as silicon, in wich quasi-longitudinal wave velocity of about 8300 m/s and quasi-transverse wave velocity of about 5700 m/s, with the secjnd quasi-transverse wave do not occur because of the problem statement [9] . Because of the high velocities of quasi-longitudinal and quasi-transverse waves total internal reflection effect at the anisotropic bottom manifest itself even more strongly than the isotrohic bottom.

Figure 22. The normalized series of first five reflected impulses in the waveguide with the hard elastic bottom.

Figure 23. The normalized series of first three reflected impulses with the harmonic filling in the point Q.

5. Diffraction Pulse Sound Signal on Non-Analytical Elastic Scatterer Put in Plane Waveguide with Hard Elastic Bottom

Based on the obtained solution, consider a more general problem of the diffraction of the pulsed sound signal on an elastic scatterer as a finite cylindrical shell, supplemented with two hemispherical shells (Figure 24) and placed in the waveguide with an elastic hard bottom (Figure 5), using the method of integral equations [16] - [21] . We note that a similar problem can be solved with the help of other methods: the boundary element method [20] [22] ; T-matrix method [23] ; the method of a potential [24] ; the finite element method [25] ; the method of Green’s functions [26] .

The first stage will solve the problem of the diffraction of a harmonic wave on a such shell. The density of the material of the shell is ρ1, the Lame’s coefficients―λ and μ. The shell was filled in the internal liquid medium with the density ρ2 and the sound velocity C3 and it was placad in the external liquid medium with the density ρ0 and the sound velocity C0. At the shell falls the plane harmonic wave with pressure pi under the angle Θ0 and with the wave vector k .

As was shown in [16] - [21] , the initial equation is integral equation, having the sense of the generalized Huygen’s principle, for the displacement vector u ( r ) of the elastic shell:

u ( r ) = S { t ( r ) G ( r ; r ) u ( r ) [ n ^ Σ ( r ; r ) ] } d S ( r ) , r V , (19)

where t ( r ) = n ^ T ( r ) is the stress vector; n ^ n ^ ( r ) = n ( r ) is the single vector of the external along the relation to S normal; T ( r ) is the stress tensor of the isotropic material; G ( r ; r ) is the displacement Green’s tensor; Σ ( r ; r ) is the stress Green’s tensor; if r concerns to the point of the surface S, in the left part of the Equation (19) will stand u ( r ) / 2 .

The displacement vector u ( r ) , the stress tensor T ( r ) , the displacement Green’s tensor (Figure 25).

Figure 24. The elastic shell in the form of the terminal cylinder with the semi-spheres.

Figure 25. The mutual disposition of the impulse point-sources and scatterers in the plane waweguide.

The displacement vector u ( r ) , the stress tensor T ( r ) , the displacement Green’s tensor G ( r ; r ) and the stress Green’s tensor Σ ( r ; r ) were connected between themselves by the following correlations [16] [21] :

T ( r ) = λ I u ( r ) + μ ( u + u ) , (20)

where I = I L + I T ; I L = ( ) / 2 ; I L I T = 0 ; I T = [ ( I ) ] / 2 , I L and I T are the longitudinal and transverse single tensors for the Hamilton’s operator ;

( r ; r ) = λ I G ( r ; r ) + μ [ G ( r ; r ) + G ( r ; r ) ] ; (21)

G ( r ; r ) = ( 1 / 4 π ρ t ω 2 ) { k 2 I g ( k 2 | r r | ) + [ g ( k 1 | r r | ) g ( k 2 | r r | ) ] } , (22)

where k1 and k2 are the wave numbers of the longitudinal and transwerse waves in the material of the shell; g ( k 2 | r r | ) = exp ( i k 2 | r r | ) / 4 π | r r | is the Green’s function.

The second integral equation presents the Kirchhoff integral for the diffracted pressure p Σ ( P 1 ) in the external medium:

C ( P 1 ) p Σ ( P 1 ) = S a { p Σ ( Q ) ( / n ) [ exp ( i k r 0 / r 0 ) ] [ exp ( i k r 0 / r 0 ) ] ρ 0 ω 2 ( u n ) } d S a + 4 π p i ( P 1 ) , (23)

where p Σ ( P 1 ) = p i ( P 1 ) + p s ( P 1 ) ; p s ( P 1 ) is the scattered pressure in the point P1; C(P1) is the numerical coefficient, equal 2π, if P 1 S a and 4π, if P1 out Sa; Sa is the external surface of the shell; Q is the point of the external surface of the shell.

For the pressure p 2 ( M 1 ) in the internal liquid medium in the point M1 is got the third integral equation:

C ( M 1 ) p 2 ( M 1 ) = S b { p 2 ( Q ) ( / n ) [ exp ( i k r 3 ) / r 3 ] [ exp ( i k r 3 ) / r 3 ] ρ 0 ω 2 ( u n ) } d S b , (24)

where Q is the point of the internal surface of the shell;

C ( M 1 ) = { 4 π , if M 1 out S b ; 2 π , if M 1 S b ;

Sb is the internal surface of the shell.

To the integral Equations (19), (23) and (24) are added the boundary conditions on the external (Sa) and internal (Sb) surfaces of the shell:

1) at the both surfaces of the shell the tangent stresses are equally null:

τ i | S a = 0 ; τ i | S b = 0 ; i = 1 , 2 ; (25)

2) the normal stress σ n at the external surface of the shell is equally the diffracted pressure pΣ, but at the internal surface is equally the pressure p2

σ n | S a = p Σ ; σ n | S b = p 2 ; (26)

In the conformity with the conditions (7) and (8) the stress vector t ( r ) in the Equation (1) is equal:

t ( r ) = p Σ n | S a ; t ( r ) = p 2 n | S b ; (27)

3) the continuity of the normal component of the displacement at the both boundaries of the shell:

u n = ( 1 / ρ 0 ω 2 ) ( p Σ / n ) | S a ; u n = ( 1 / ρ 2 ω 2 ) ( p 2 / n ) | S b . } (28)

The substitution of integral Equations (19), (23) and (24) in the boundary conditions gives the system of equations in terms of unknown functions pΣ, p2 and the components of the displacement vector u at the both surfaces of the shell. To obtain numerical solution of this system the integral equations are replaced the quadrature formulas and the grid of the nodal points is chosen at both surfaces of the shell as well as it has be done for the ideal non-analytical scatterers [16] [17] .

For choosing boundary conditions we will have the integrals of the two types: the integrals with the isolated special point and the integrals which are considered of the sence of the principal meaning. The method of the calculation of the second types was described in [16] .

Thus calculated reflection characteristics of the harmonic signal with frequency ν can determine the spectral reflectance function S S ( 2 π ν ) and it can help be applying a Fourier transform we obtain a temporary function of the reflected pulse Ψ S ( t ) [27] :

Ψ S ( t ) = 1 π Re 0 S S ( 2 π ν ) e + i 2 π ν t d ( 2 π ν ) (29)

Similarly using spectral reflectance characteristics of elastic bodies of spheroidal form [16] [28] [29] [30] [31] , we can compute sequences of reflected pulses in the waveguide with hard elastic bottom and for their.

6. Conclusions

In the first part of the review we investigate an interaction of a scatterer and an interface between media; it is shown that a main role in this is played by interference effects. The second part of the review is devoted to a study of a spektrum of a scattered field of an ideal prolate spheroid placed in an underwater sound channel with non-reflecting boundaries. The third part of the review demonstrates the effect of a bottom structure on a series of pulses, reflected from a spheroidal body located in a plane waveguide.

As a result of the research we can draw three conclusions:

1) During studying propagation and diffraction of pulse signals in a plane waveguide, one needs to use the method of imaginary sources as pulses like bundles of energy spread to any directions (including and along the axis of the waveguide) with the group velocity which does not exceed the sound velocity, namely the group velocity based the method of imaginary sources;

2) Replacing the hard elastic bottom on the absolutely hard bottom is acceptable to those sources (real and imaginary) from which waves in the fall to the hard elastic bottom try total internal reflection.

3) We have adopted the model of image sources and image scatterer is gute acceptable (due to internal reflection), at least, for first five calculated reflected pulses in a plane wave-guide with hard elastic bottom.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Kleshchev, A.A. (1977) Scattering of Sound by Spheroidal Bodies Located at Interface between Media. Soviet Physics. Acoustics, 23, 404-410.
[2] Sommerfeld, A. (1950) Differential Equations in Partial Derivatives of Physics. Foreign Literature, M., 456 p.
[3] Ivanov, E.A. (1968) Diffraction of Electromagnetic Waves on Two Bodies. Science and Technology, Minsk, 583 p.
[4] Kleshchev, A.A. (1979) Scattering of Sound by Spheroidal Body Located at Interface between Media. Soviet Physics. Acoustics, 25, 143-145.
[5] Tolstoy, I. and Kley, K.S. (1969) Acoustic of Ocean. World, M., 304 p.
[6] Kleshchev, A.A. and Klyukin, I.I. (1974) The Spectral Characteristics of the Scattering of the Sound by Body, Placed in the Sound Channel. Soviet Physics. Acoustics, 20, 470-473.
[7] Kleshchev, A.A. (2015) Diffraction of Pulse Sound Signals on Elastic Bodies of Spheroidal Form Put in Plane Waveguide. MIT, Vol. 2, 77-81.
[8] Kleshchev, A.A. (2017) Analytical and Non-Analytical Scatterers on Plane Waveguide with Hard Elastic Bottom, Irradiated by Pulse Sound Signal. American Journal of Modern Physics, 6, 51-56.
[9] Lekhnitskiei, S.G. (1977) Theory of Elasticity of Anisotropic Elastic Body. Science, M., 416 p.
[10] Brechovakikh, L.M. (1957) Waves in Laminated Mediums. Publ. Acad. Scien. SSSR, M., 502 p.
[11] Kleshchev, A.A. (2015) Pulse Sound Signals Diffraction on Elastic Bodies of Аnalytical and Non-Analytical Forms, Put in Plane Waveguide. Zeitschrift fur Naturforschung A, 70, 419-427.
[12] Kleshchev, A.A. and Kuznetsova, E.I. (2011) Scattering of Pulse Sound Signals by the Spheroidal Body, Put in Plane Wavegui de. Coll. Proc. Russ. Acoust. Soc. XXIV Session. GEOS, M., Vol. 1, 743-745.
[13] Kleshchev, A.A. and Kuznetsova, E.I. (2012) Diffraction of Impulse Signals on Spheroidal Body, Put in Plane Waveguide. International Journal of Theoretical and Mathematical Physics, 2, 211-214.
[14] Kleshchev, A.A. (2013) Diffraction of Pulse Sound Signals on Elastic Spheroidal Shell, Put in Plane Waveguide. Advanced Studies in Theorrsetical Physics, 7, 697-705.
[15] Kleshchev, A.A. (2014) Diffraction of Sound Signals at Elastic Shell of Non-Analytical Form Put in Plane Waveguide. Advances in Signal Processing, 2, 46-49.
[16] Kleshchev, A.A. (2012) Hydroacoustic Scatterers. Prima, S.-Pb., 268 p.
[17] Кleshchev, A.A. (1989) Scattering of Sound by Ideal Bodies of Non-Analytical Form. Тr. Lenin. Korablestr. Inst. Generalship. Syst., 95-99.
[18] Кleshchev, А.А. (1993) Method of Integral Equations in Problem of Sound Diffraction on Elastic Shell of Non-Analytical Form. Journal of Technical Acoustics, 2, 65-66.
[19] Кleshchev, А.А. (2013) Мethod of Integral Equations in Problem of Sound Diffraction on Bodies of Non-analytical Form. Мarine Messenger, 2, 94-98.
[20] Seybert, A.F., Wu, T.W. and Wu, X.F. (1988) Radiation and Scattering of Acoustic Waves from Elastic Solids and Shells using the Boundary Element Method. The Journal of the Acoustical Society of America, 84, 1906-1912.
[21] Podstrigach, J.S. and Poddubnjak (1986) Scattering of Sound Beams of Spherical and Cylindrical Form. Naukova Dumka, Кiev, 264 p.
[22] Brebbia, C.A. and Walker, S. (1982) Boundary Element Techniques in Engineering. World, М., 242 p.
[23] Peterson, B. and Strom, S. (1975) Matrix Formulation of Acoustic Scattering from Multilayered Scatterers. Journal of the American Statistical Association, 57, 2-13.
[24] Кupradze, V.D. (1963) Methods of Potential in Theory of Elasticity. Fizmatgiz, М., 472 p.
[25] Dushin, A.Y., Il’menkov, S.L., Kleshchev, A.A. and Postnov, V.A. (1989) Use of Finite El Sement Method to Solution of Problems of Sound Radiating by Elastic Shells. Proc. All-Union Symp. Interaction of Acoustical Waves with Elastic Bodies, Таllinn, 89-91.
[26] Il’menkov, S.L., Kleshchev, A.A. and Klimenkov, A.S. (2014) The Green’s Function Method in the Problem of Sound Diffraction by an Elastic Shell of Noncanonical Shape. Acoustical Physics, 60, 579-586.
[27] Kharcevich, A.A. (1957) Spectrum and Analisis. GITTL, М., 236 p.
[28] Кleshchev, А.А. (1986) Three-Dimensional and Two-Dimensional (Axis-Symmetrical Characteristics of Elastic Spheroidal Scatterers. Akusticheskii Zhurnal, 32, 268-270.
[29] Кleshchev, A.A. and Klyukin, I.I. (1987) Principles of Hydroacoustics. Shipbuilding, L., 224 p.
[30] Кleshchev, A.A. (2012) Debye and Debye-Type Potentials in Diffraction, Radiation and Elastic Wave Propagation Problems. Acoustical Physics, 58, 338-341.
[31] Кleshchev, A.A. (2014) Resonance Scattering of Sound by Spheroidal Elastic Bodies and Shells. Acoustical Physics, 60, 253-261.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.