11.07.2014 Views

Flute acoustics: measurement, modelling and design - School of ...

Flute acoustics: measurement, modelling and design - School of ...

Flute acoustics: measurement, modelling and design - School of ...

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.

2.2. ACOUSTICS OF WOODWIND INSTRUMENTS 15<br />

infinitesimal width dx <strong>and</strong> area S is E(x)dx where the energy function E(x) is given by<br />

E(x) = 1 2<br />

2.2.3 Wave propagation in pipes <strong>and</strong> horns<br />

( p 2 S<br />

ρc 2 + ρU 2 )<br />

. (2.2)<br />

S<br />

The acoustic pressure wave inside wind instruments with rigid walls can be determined generally<br />

by solving the wave equation<br />

∇ 2 p = 1 ∂ 2 p<br />

c 2 ∂t 2 (2.3)<br />

subject to the boundary condition that n·∇p = 0 at the instrument walls, where n is a unit vector<br />

normal to the wall. Simple solutions to this equation exist for several geometries, including<br />

cylindrical tubes <strong>and</strong> conical horns (Fletcher & Rossing 1998). Approximate solutions can be<br />

found for other geometries, such as the exponential horn. The acoustic flow velocity u is related<br />

to the pressure according to Newton’s law<br />

u =<br />

i ∇p (2.4)<br />

ρω<br />

where ω = 2πf is the angular frequency.<br />

The solution to (2.3) in cylindrical polar coordinates for a cylindrical pipe <strong>of</strong> radius a is proportional<br />

to exp[i(−k mn x + ωt)], where the proportionality depends on r <strong>and</strong> φ <strong>and</strong> describes<br />

how the pressure wave is distributed over the (r,φ) plane. The square <strong>of</strong> the wave vector k mn<br />

for the mode (m,n) is given by<br />

) 2 ( πqmn<br />

) 2<br />

− , (2.5)<br />

( ω<br />

kmn 2 = c a<br />

where m <strong>and</strong> n are mode numbers corresponding to the numbers <strong>of</strong> nodal diameters <strong>and</strong> nodal<br />

circles in the (r,φ) plane <strong>and</strong> q mn is defined by the boundary condition that there can be no flow<br />

in the radial direction at the cylinder walls.<br />

The (m,n) mode propagates (extends infinitely in the x direction, neglecting any losses) if<br />

k mn is real; otherwise the acoustic pressure decays exponentially with x <strong>and</strong> the mode is termed<br />

evanescent. The plane wave mode given by m = n = 0 will always propagate since q 00 = 0;<br />

higher modes will only propagate at frequencies above the cut<strong>of</strong>f frequency for the mode given<br />

by<br />

ω c = πq mnc<br />

. (2.6)<br />

a<br />

The higher mode with the lowest cut<strong>of</strong>f frequency is the (1, 0) mode, which has a single nodal<br />

plane. For this mode<br />

ω c = 1.83c<br />

a , (2.7)<br />

giving a cut<strong>of</strong>f frequency <strong>of</strong> approx. 10 kHz for a typical flute bore radius <strong>of</strong> 10 mm.<br />

Similar equations may be derived for conical horns although the cut<strong>of</strong>f frequency for each<br />

mode varies with radius along the length <strong>of</strong> the horn, <strong>and</strong> the wavefronts for the always-propagating<br />

mode are spherical rather than planar.<br />

Fortunately, in the bores <strong>of</strong> most wind instruments all higher modes are non-propagating at<br />

frequencies <strong>of</strong> musical interest. Higher modes are evoked in the vicinity <strong>of</strong> bore disturbances,<br />

such as discontinuities or tone holes, but decay quickly. The fields around tone holes may be

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

Saved successfully!

Ooh no, something went wrong!