28.02.2013 Views

Introduction to Acoustics

Introduction to Acoustics

Introduction to Acoustics

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

74 Part A Propagation of Sound<br />

Part A 3.11<br />

the scattered wave. This term is required <strong>to</strong> be made up<br />

of waves that propagate out from the sphere, which is<br />

centered at the origin. With some generality, one can use<br />

the principle of superposition and set<br />

�<br />

ˆpsc = ˆPinc<br />

ℓ=0<br />

Aℓ Pℓ(cos θ)h (1)<br />

ℓ (kr) , (3.425)<br />

with the fac<strong>to</strong>rs in the terms in the sum representing the<br />

Legendre polynomials and spherical Hankel functions.<br />

The individual terms in the sum are sometimes referred<br />

<strong>to</strong> as partial waves. The analytical task is <strong>to</strong> determine the<br />

coefficients Aℓ. (The expression above applies <strong>to</strong> scattering<br />

from any axisymmetric body, including spheroids<br />

and penetrable objects, where the properties might vary<br />

with r and θ. The discussion here, however, is limited <strong>to</strong><br />

that of the rigid sphere.)<br />

With the plane wave expansion (3.423) inserted for<br />

the direct wave, the sum of incident and scattered waves<br />

takes the form<br />

ˆp = ˆPinc<br />

∞�<br />

ℓ=0<br />

� (2ℓ + 1)i ℓ jℓ(kr) + Aℓh (1)<br />

ℓ (kr)�<br />

× Pℓ(cos θ) . (3.426)<br />

the boundary condition imposed by the rigidity of the<br />

sphere is that ˆvr = 0 at r = a, or equivalently that<br />

∂ ˆp/∂r = 0atr = a. With the aid of the radial component<br />

of the Euler equation (3.175) and of the linear<br />

independence of the various Legendre polynomials, this<br />

yields<br />

d � ℓ<br />

(2ℓ + 1)i jℓ(kr) + Aℓh<br />

dr<br />

(1)<br />

ℓ (kr)� = 0<br />

at r = a . (3.427)<br />

The desired coefficients are consequently<br />

Aℓ =−(2ℓ + 1)i ℓ<br />

� ddr<br />

jℓ(kr) �<br />

r=a<br />

. (3.428)<br />

� ddr h (1)<br />

ℓ (kr)�<br />

r=a<br />

Far-Field Scattering<br />

In the limit of large kr, the asymp<strong>to</strong>tic expression (3.412)<br />

for the spherical Hankel function can be used, and the<br />

scattered wave takes the asymp<strong>to</strong>tic form<br />

where<br />

ˆpsc → ˆPinc f (θ) eikr<br />

kr<br />

f (θ) =<br />

, (3.429)<br />

∞�<br />

(−i) ℓ+1 Aℓ Pℓ(cos θ) (3.430)<br />

ℓ=0<br />

is a complex dimensionless function of only the angle θ.<br />

This asymp<strong>to</strong>tic form holds for scattering from any<br />

axisymmetric object of bounded extent. In the more general<br />

case when the object is not axisymmetric, one should<br />

regard the function f as being also a function of the<br />

azimuthal angle φ,sothat<br />

f (θ) → f (θ,φ) . (3.431)<br />

The prediction in all cases is that the far-field scattered<br />

wave resembles an outgoing spherical wave, but with an<br />

amplitude that depends on the direction of propagation.<br />

The far-field intensity associated with the scattered<br />

wave is asymp<strong>to</strong>tically entirely in the radial direction,<br />

and in accord with (3.281), its time average is given by<br />

Isc = 1 1<br />

2 ρc | ˆpsc| 2 , (3.432)<br />

or,<br />

Isc = 1<br />

2<br />

| ˆPinc| 2<br />

ρc<br />

| f (θ,φ)|<br />

2 1<br />

. (3.433)<br />

2<br />

k 2 r<br />

(ka) 2<br />

0.5 (ka) 2<br />

ka

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!