28.02.2013 Views

Introduction to Acoustics

Introduction to Acoustics

Introduction to Acoustics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

is the Wronskian for the Bessel and Neumann functions.<br />

Here the primes denote differentiation with respect <strong>to</strong> the<br />

argument η = kw of the indicated function. One can derive<br />

from the differential equation (3.453) that these two<br />

functions independently satisfy the Wronskian relation<br />

ηW(J0, Y0) = constant , (3.495)<br />

and this is so regardless of the specialized definitions of<br />

the two functions. The constant can be evaluated from<br />

values and derivatives at any given point. In particular,<br />

one can use the asymp<strong>to</strong>tic expressions, and the derived<br />

constant is 2/π. Thus one has<br />

Iw,av = 1 |C|<br />

π<br />

2<br />

� �<br />

1 1 �<br />

= �<br />

d ˆm �2<br />

ρcw ρcw 16π �ω2 �<br />

dℓ � . (3.496)<br />

3.12.3 Radially Oscillating Cylinder<br />

A classic example for the radiation of cylindrical waves<br />

is that where a cylinder has nominal radius a, andan<br />

instantaneous radius<br />

wcyl(t) = a + vo<br />

sin(ωt) . (3.497)<br />

ω<br />

Here vo/ω is the amplitude of the deviations of the actual<br />

radius from the nominal value a. For the linear acoustics<br />

idealization <strong>to</strong> be valid, it is required that this deviation<br />

be substantially less than a,sothat<br />

vo ≪ ωa . (3.498)<br />

The boundary condition on the fluid dynamic equations<br />

should ideally be<br />

vw = vo cos(ωt) at w = wcyl(t) , (3.499)<br />

but (also in keeping with the linear acoustics idealization)<br />

it is replaced by<br />

vw = vo cos(ωt) at w = a . (3.500)<br />

The corresponding boundary condition on the complex<br />

amplitude is<br />

ˆvw = vo at w = a . (3.501)<br />

If one takes the complex amplitude of the acoustic<br />

pressure <strong>to</strong> be of the form<br />

ˆp = AH (1)<br />

0 (kw) , (3.502)<br />

Basic Linear <strong>Acoustics</strong> 3.12 Cylindrical Waves 81<br />

where A is a complex number <strong>to</strong> be determined, then<br />

the radial component of the fluid velocity is<br />

ˆvw = A d<br />

iωρ dw H(1)<br />

0 (kw) . (3.503)<br />

The imposing of the boundary condition determines the<br />

value of A and one obtains<br />

ˆp = iωρvo<br />

H (1)<br />

0 (kw)<br />

� d<br />

dw H (1)<br />

0 (kw)�<br />

w=a<br />

. (3.504)<br />

3.12.4 Transversely Oscillating Cylinder<br />

If the cylinder is rigid and oscillating back and forth<br />

in the x-direction, the analysis is similar <strong>to</strong> that for<br />

the transversely oscillating sphere. The complex amplitude<br />

of the acoustic part of the pressure should have<br />

the general form<br />

ˆp = B d<br />

dx H(1)<br />

d<br />

0 (kw) = B cos φ<br />

dw H(1)<br />

0 (kw) ,<br />

(3.505)<br />

where B is <strong>to</strong> be determined from the boundary condition.<br />

The latter can also be written<br />

ˆp =−kBcos φH (1)<br />

1 (kw) , (3.506)<br />

in accord with the recursion relation (3.471). The latter<br />

version is of the standard form for a fac<strong>to</strong>red solution of<br />

the Helmholtz equation in cylindrical coordinates.<br />

Euler’s equation gives the radial component of the<br />

fluid velocity as<br />

ˆvw = 1<br />

iωρ<br />

B cos φ d2<br />

dw<br />

H(1)<br />

2 0<br />

(kw) . (3.507)<br />

The appropriate approximate boundary condition is that<br />

ˆvw = vc cos φ at w = a . (3.508)<br />

This allows one <strong>to</strong> identify<br />

B =<br />

iωρvc<br />

� d 2<br />

dw 2 H (1)<br />

0 (kw)�<br />

w=a<br />

. (3.509)<br />

Part A 3.12

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

Saved successfully!

Ooh no, something went wrong!