JOURNAL OF FLUTE ACOUSTICS

Volume 1 2006, Number 1, July 2006

 

Articles will be published without peer review if requested.  Submit copy for publication as an email attachment to <jellymuscles@adelphia.net>

 

Table of Contents

Review Article

The Physics of Flutes

by Nelson McAvoy

Introduction

            Modern woodwinds such as the clarinet, trumpet, saxophone, bagpipe, and flutes, were not invented, they evolved.  They became what they are by thousands of changes over generation with manufacturing and material advances to improve quality.    Today flute design continues.  There is still a lot of trial-and-error in design today, even though we completely understand the physics of instruments.  This is because mathematical solutions always assume simplified or idealized situations.  Luckily, the flute is simpler than other woodwinds. A flute is basically a tube open on both ends.  The sound made by other woodwinds is more complicated and less predictable because:

      So we will be able to make rather good predictions from calculations and knowledge of the physics of the flute that will be of value to designers, much more so than other woodwinds.   This book is geared to the flute designer and requires no knowledge of higher mathematics beyond trigonometry.  Other mathematics, such as complex variables will be explained as we go.   Knowledge of the physics of the flute is helpful for the designer mostly in a qualitative, instead of a quantitative way, in that it gives an appreciation of what quantities change due to modification without going through extensive calculations. 

For completeness and for those with engineering backgrounds, we have added appendices giving development of equations from basics.  We have chosen this format because many of the readers are very familiar with some of the design concepts and others are hearing them for the first time.  There is no other way to make the book complete and coherent for readers of both kinds. 

There is no way to understand how a flute works with out understanding how pressure and current are phase related and the consequence of this phase relation.  This, and an understanding of basic trigonometry are all the tools needed.  Readers that are not familiar with these concepts should go to Appendix A before delving into Chapter 2.


 

Chapter 1

Oscillator Concepts

   An oscillator is a device of repetitive motion.  Musical instruments are oscillators.   They are as old as human history.  But for the last two centuries they have become sophisticated and stylized.  Their workings have been explained completely by Newtonian physics.  The history of the last century has been completely shaped by a new family of oscillators-electronic oscillators.  Not until the extensive development of electronic engineering was musical instrument engineering formalized.  These concepts were not introduced into books on flute engineering until (Fletcher and Rossing 1997) as based on the formalism first presented by (Kinsler and Frey 1942).  This is why the nomenclature, perspective, and concepts used now in flute engineering are those used in electronic oscillators.   This formalized treatment makes for easier understanding of musical instruments, especially flutes.

  Before getting into flute design, it is helpful to have a perspective of oscillator concepts and history.

 When stored energy switches from potential (stationary) energy,  force x distance, to kinetic (motion) energy, ½ mv2,  it is called an oscillator, or resonator. 

Some examples of oscillators are:

1.      A weight bobbing on a spring.

2.      A bouncing ball.

3.      Molecules.  For example, the Nitrogen of an (ammonia molecule) vibrating back and forth from one side of the middle of the triangle made by the three H atoms.  Ammonia is blown over the embouchure of a flute tube to make the most accurate tone there is.  It drifts by one cycle in 1013 years.  Charles Towens got the Nobel Prize for this flute, it is called the ammonia maser.  The flute has been taken on space journeys and checked when it returned to measure the Einstein prediction that it would “age”.  It did not.  It just kept on tooting at the same note. 

4.      A small tube with mirrors on the end and filled with fluorescent helium and neon.  This like any flute oscillates at multiples of  c/2L, in the 632,8 nmeter range (orange light).

5.      A clock.  Any clock, like the tuning fork watches made in the mid-20th century by the Hamilton Watch Co (it hums).  A pendulum clock.  Today’s watches have a transistor in series with a coil and capacitor. 

6.      A musical instrument.

7.      Water swishing back and forth in a bathtub.

8.      A radio transmitter or the many electronic oscillators in computers. 

9.      A suspension bridge, sometimes when it is not designed right.

10.  A boat when it is rocking back and forth.

11.  A microwave oven.

12.  An ordinary playground swing.

 

A resonator that runs by itself is called a perpetual motion machine.  There is no such gadget. Oscillators all have universal aspects and in their design and study, engineers use well defined nomenclature.  The commonly used ones are:

Table 1 gives some examples of oscillators and their characteristics.

 

Oscillator L/D* Frequency

Range

Potential

Energy

Kinetic

Energy

Quality

Factor,Q

Feed

Mechanism

Remarks
Weight on

 a spring

L 1-20,000

Hertz

Spring con-

stant x distance

50 to 1000 Push v is speed
Basketball L 5 Hz Height x

weight

10 Magic

Johnson

 
NH3 maser D 1.2 GHz Spring con-

stant x distance

10,000 Molecular beam

over embouchure

 
Wind up

Watch

L 1 Hz Torque x angle 50 Spring-ratchet  is angular

speed

Electronic

oscillator

L 60 Hz to

1 THz

100-1000 Transistor C-capacitance

L-inductance

Microwave

oven

D 3 to 10 GHz 1000 Electron beam

embouchure

Z is impedance of space
Poorly designed bridge D 1 Hz Same as violin string Same as violin string 30 wind  
Rocking boat L 1 Hz Torque x angle 20 Waves and

wind

 
Flute D 3 octaves 80 embouchure 2 ends open
Woodwind D 3 octaves 80 Vibrating end One end closed
He-Ne laser D .6 wave-

length

1000 Plasma induced amplifying Red-orange

light

Harmonica L 4 octaves Rod bend Rod vibrate 100 embouchure  
String Instrument D 4 octaves Stretched string Vibrating

string

10,000 Bow or pluck Ring time indicates Q

                        Table 1 Examples of oscillators.  L and D stand for lumped and distributed elements.

 

 

 

Chapter 2

Elements of Flute Design

 

    Appendix A familiarizes you with the description of sound waves and how they are described mathematically.   The Appendix explains how waves are described in trigonometric functions sine, cosines, and tangents in a pipe.  It also explains how waves are described with imaginary numbers to keep track of phase relations.  If you are familiar with these concepts, we can start right in describing the impedance of a sound wave in a pipe.  Impedance is the ratio of pressure to current .  It was shown in Appendix A that a plane wave traveling in space has an impedance of .   We can think of u as the actual speed of a wind in meters/second.  Of course in a sound wave it is the speed of the wind going back and forth in phase with the pressure.  The impedance in a pipe either very long, or terminated with an absorbent material such as cotton wool,is

                                                                        .                                                                (2.1)  

where S is the area of the pipe.  When a sound wave is launched down a pipe  is the density of air in kilograms per cubic meterat room temperature or 1.293 at .  c is the speed of sound which is  meters/sec. u is the additional speed of the molecules in the air just as the pressure p is the additional pressure beyond atmospheric pressure of 14.7 and T is in Celsius degrees, usually 25.  At sea level and atmospheric pressure and 77F, it is about 344.86 M/s.  See Appendix B for details.  If S is the area of the pipe and u is speed in meters per second, then uS is the speed per unit volume.  So in a pipe, if then the pressure and current are always in phase.  With no reflection in the pipe, their description in time and space are exactly the same.   More specifically Equation 2.1 is,

                                                                                                      (2.1a)    

is the phase and determines what the pressure is at t=0 and x=0 and .  We can also describe the wave with out imaginary numbers, i.e. without j, as

                                                                      (2.1b)

or

                                     (2.1c)

Our next relation, Equation A28 from Appendix A is used extensively in flute design,

           .                                          (2.2)

It is the general transmission line equation.  It is applicable not only for sound waves in a pipe but for electro-magnetic waves in a cable, power transmission lines, and a wave traveling down a taught string.    Consider a pipe of constant area, S.  Choose a location down the pipe,.  Choose another location to the left of  , the “in” location or where.  The impedance “looking in” at x=0 is the of Equation 2.2.

We first illustrate the use of Equation 2.2.  With absorption material starting at .  There is no reflection and the load impedance at  will be  .  When substituting this for in Equation 2.2 ,  which is expected from Equation 2.1 if there is no reflection.

Figure 2.1  Pipe filled with absorbing material so there is no reflection of a sound wave launched for the left end. 

 

Open End Pipe

We now want to start talking about reflections in the pipe, i.e., situations where the pipe is terminated by something other than an absorbent material or is infinitely long.  Not only do we want to terminate it but we want to terminate is so that the pipe is a resonator, a flute.  In Chapter I  it was stated that all oscillators and resonate cavities are constructed so that the kinetic energy, current in our case, and the potential energy, pressure in our case, are out of phase and the energy switches back and fourth at the frequency of the flute.  A pipe open on both ends, as shown in Figure 2.2 satisfies this condition.  At a distance of 0.61a the pressure at the end of the pipe has become zero (we explain this later).  Therefore the load impedance, in Equation 2.2 so,

Figure 2.2  At a distance out from the end of an open pipe 0.61a, the pressure from inside the pipe has diminished to zero.  Actually, 0.6a is a fictitious number that assumes that it diminishes all at once. 

                                              (2.3)

From the definition of impedance as the ratio of pressure to current, Equation 2.3, is also

                                                    (2.4)

Figure A6 and an explanation in Appendix A shows that an imaginary quantity is out of phase with a real part.  Or another way to look at it is that  from the definition.   Now the  and , therefore, .  Note that p=0 when is and so on.  Or, in other words, when  is  because  This is the resonance condition, that the pipe has to be a multiple of a half wavelength. 

 We also want to show that the impedance of two pipes cascaded is the same as a single pipe of the same diameter.  If this were not true there would be something wrong with our equation for open pipe impedance.  This will be called series impedance, although it is not the sum of two impedances as one has with lump elements.  Sometimes we label this impedance as in distinction tofor  the “shunt” impedance which is the impedance of two branched pipes, i.e. parallel impedance. 

                        Figure 2.1. Impedance of a tube as a series of cascaded impedances.

 

If Equation 2.2 is valid, than the impedance at the location    is  .  Equastion 2.3 becomes ,

                                     (2.5)

Using as the new load impedance for the end of the pipe at , we have from Equation 2.2 for ,

.      (2.6)

Multiplying top and bottom by  gives,

    .                                   (2.7)

We note from the trigonometric identity in Appendix D that

.  Comparing this with Equation 2.6  results in

                                   (2.8)

which proves that the impedance of open pipes are additive if they have the same area S. 

We will have many occasions when the length is not more than the diameter.  I.e., when  and is a small angle.  Then  (if we use radians, not degrees)  For example a key hole of chimney height t is a tube of impedance .  In this situation the concept of a transmission line whose length is more than its diameter, where Equation 2.2 is valid, breaks down.  We will find out from experiments done by Keefe (1982) who measured the impedance of keys and key holes, where an experimentally determined   are valid concepts; where   for closed holes is for an open short tubes.  In other words can be used in the same manner as d is in Equations 2.3 through 3.8.  is a series impedance for a tone hole even though the main impedance element of a tone hole is a shunt (parallel) impedance which will be discussed under Branched Pipes below and depicted in Figure 2.6.

Closed End Pipe

                                             Figure 2.2  A closed end pipe of length l and area S =p a2 .

                       

The closed end of a pipe has the impedance of  because the pressure is high and the current is zero.  In this case Equation 2.2 becomes,

      (2.7)

For short pipes where  and ,

,                                                     (2.8)

where V is the volume of the closed pipe S and .     is called acoustic compliance or acoustic capacitance because a closed volume in acoustics behaves analogously to an electrical capacitance.  The impedance of an electrical capacitance is.  So the acoustical equivalent of an electrical capacitance, in phase and form is .  The analogy is not complete because you can not cascade two closed end pipes and you can cascade electrical capacitances.  It is a convenient concept in flutes because a closed key hole is just that, a very short pipe closed on one end and branched to the flute body.  We do not stack one closed key on top of another so that where the analogy with series electric capacitors brakes down is no concern to us.

Branched Pipes                                                   

The embouchure of a flute, the key holes, opened or closed, and the stopped left end are examples of branched pipes.  The concept of impedance implies a specific place and the direction in which you are “looking”.  For example, the impedance of a pipe that is conical in shape is not the same looking from the small end as when looking from the large end. In the cone, the pressure to current ratio will be different at each end.   When a pressure wave in a pipe branches, the current is split so  , looking from the  end.  The pressure, on the other hand, remains the same at the junction, so ,

.                                             (2.9)

In flutes we terminate the end of the branch pipe in two ways.  We open and close the keys.   The left (embouchure) end is closed.  The embouchure is an open branch, it is partly covered by the lips.  And, of course, the other end is always open. 

Figure 2.3 A branched pipe with lengths, , and;and cross sectional

areas, , and;and impedances, , and  at the ends.

 

Key Hole Impedance Measurement 

Key hole impedance measurement was made with an ingenious experimental arrangement by Dougles Keefe (J. Accoust. Soc. Am. Vol. 72, No. 3, Sept. 1982).  The experimental arrangement is shown in Figure 2.4.  The lowest order resonance has the pressure equal to zero at the tone hole being tested.  In that situation there is a small shortening of the pipe.  After resonance the plunger is backed off a measured amount to maintain resonance when the hole is eliminated by plug.  This represents a series impedance  where  is the change in the plunger length.  The a is for series (read length) change and (0) is for open tone hole, in contradistinction to s for shunt (parallel) and (C) is for a closed tone hole.  The shunt impedance is measured. without the adjustment of the plunger, by observing the length  at the next highest (symmetric) resonance.     

Figure 2.4 Experimental arrangement by Keefe for measuring the tone hole impedance.  Unfortunately for flute design, the pipe bore diameter was 40 mm and the tone hole was similar to that of a clarinet.  The height, t, was 6.4 mm and the key hole radius was 13.2 mm.   

 

Open Hole Asymmetric Mode

The end termination was adjusted to insure symmetry around the tone hole by adjusting the termination assembly for minimum radiation out of the open tone hole The tone hole was then plugged and the cavity length was increased for resonance.  The fact that the length has to be increased from   to 2means that one effect of the open hole was to acoustically decrease  the tube impedance by

                                                 (2.10)

on the right of the key and

                                                 (2.11)

on the left of the key. This means that when designing a flute, the row of open holes will each decrease the length of the flute by   This is very non-consequential because the shunt effect of open holes dominate the sound created.   In contrast, when we get to the calculation of  for the closed hole asymmetric case below it is very consequential in the design because a row of closed holes above the open holes will effectively shorten the flute by   times the number of closed holes.  The theoretical calculation of this given by Keefe(1982) is,

                                      (2.12)

The resonance condition is, on each side of the center.  This means that .  This is the lowest order mode and is called an asymmetric mode because the pressure on one side of the center is out of phase with the pressure on the other.   Figure 2.5.

Figure 2.5 The asymmetric mode is the lowest frequency mode in a cavity closed on both ends and a minimum of radiation is transmitted out the key hole because the pressure is a minimum there.  The next highest mode is symmetric,  is adjusted to maximize the power out of the key hole and is about twice the frequency as the asymmetric mode. 

 

Open Hole Symmetric Mode

The next highest mode at double the frequency is a symmetric mode with a pressure maximum at the key hole as in Figure 2.5.  We will see that this arrangement is used to measure the branch (or shunt) effect of the hole.   Considering both arms of the cavity as branches looking into the open key hole, from Equation 2.9 we have

                                                       (2.13)

The left and right pipes have been adjusted to be the same so we have,

                                                    (2.14)

We can easily measure  and and ,,in meters and we have the effective key hole height, , to use in calculation of  the tuning effect resulting from opening a key hole.  Keefe also went through the rigorous theoretical treatment and came up with the calculation of as,

                                             (2.15)

He found good agreement but, unfortunately, he used key geometry for clarinet type keys.  We will use this when we start designing flutes.

 

 

 

Closed Hole Asymmetric Mode

Unlike the open hole asymmetric (series) parameter calculated above, the closed hole asymmetric (series) parameter is used to design the length of all flutes because all flutes have a row of closed holes that shorten the effective length.  The measurement method was identical to the open hole method.  The cavity was adjusted for resonance with the tone hole covered and sealed.  The tone hole was then plugged and the cavity length increased for resonance was measured.  The fact that the length has to be increased from   to 2means that one effect of the closed hole was to acoustically decrease  the tube length by and, consequently, the impedance by

                                                 (2.16)

on the right of the key and

                                                 (2.17)

on the left of the key.

The calculated value was,

                               (2.18)

 

Closed Hole Symmetric Mode

The closed hole symmetric mode measures the shunt impedance and it was found, as expected that the measured was the same as the physical depth of the closed hole, t. 

 from            (2.19)

 For closed key holes chimney high of 5 mm on wooden flutes and 2 mm on concert flutes, the shunt admittance 1/is  close to zero.  For example a concert flute with t = 2xmeters,  and , gives a     of j(1000)(2x) (7.7x)=1.6x.  One case where the  is a large enough to be significant is the shunt impedance of the cork end.  Notice that in this apparatus the lowest resonance frequency, the asymmetric mode, measures series impedance and the addition or subtraction of an extra pipe length.  In the apparatus’ next highest resonance  frequency, the symmetric mode, measures the shunt or branch effect. 

Figure 2.6.  A series impedance shortens the pipe, either closed or open hole situation.  The shunt impedance changes the effectiveness of a branch pipe.  The open hole branch is lengthened.  The closed hole branch is not changed. 

 

 

 

 

 

Flute examples

Example 1 Thin Walled Celtic Flute

 

Symbles used in Thin Wall Flute analysis (all distances in metes)

f        Frequency.  293.66 Hz for the low note and 329.63 Hz for the low note.

   c       Really should be  M/sec and ,T is in Celsius degrees.  We use 344.7   

            on this particular day.  Measured lengths change a millimeter or so for every degree change .

 

        Density of air at standard temperature and pressure.  1.21 Kg/.

   S       Area of pipe in

   t        Height or length of pipe.

      Length of open pipe used as a shunt. 

        Wavelength in meters.  .      

   wave number in reciprocal meters also = =5.35 or f=293.7 and 6.0 for f=329.63 and c=344.7.

   Q Energy stored / energy loss pre cycle.  .  are  at 3db off resonance.

       Decay time. .  . . has units of radians/sec. .

       Attenuation constant .  Units are meters/sec.

   a       Radius of the flute in meters.

   b       Radius of the branch in meters.

      Length of a substitute open pipe that extends from the center of the embouchure to the left.  It

            functions the same as the corked end of the flute.

     Distance from the center of the embouchure to the end of the flute.

        Distance beyond the right end of the flute, ,where the pressure is zero.

        Distance from the center of the lowest key to the end of the flute.

       Sometimes written as .  It is the substitute distance for the lowest key.  Abolish the last key and      

            replace it with a pipe from the center of the key long.

    the for the previous key location calculation.

      Distance from the substitute distance -.

 

As the first example we will calculate the geometry of a quality,  reputable Celtic flute that is not a complex geometry.  It is the 6 key hole flute of thinner walls than most Celtic flutes because it is made of plastic pipe instead of thicker walled wooden flutes turned on a lath.  Therefore the branch impedance for the key holes is for a very short length pipe of negotiable contribution.   You can look at this flute on  <http://home.earthlink.net/~life2all/dougswebspace/index.html> .

    First we do the substitute tube for the embouchure.  Figure 2.7 shows the branch.

 

            Figure 2.7  The embouchure for this flute is  a round 1 cm diameter hole that is 4.57 mm deep (t= 4.57x10M).  It is standard PVC pipe of 20.8 mm diameter.  The cork is 19 mm from the center line of the embouchure.   It has 6 finger holes in the key of D major.

 

    The cork end impedance has the form of a compliance (capacitance) because it is closed at one end and has the form as shown in Equation 2.19

                                      (2.20)                                      

The embouchure end of the branch is open and is of the form of an acoustic inertance (inductance) as in Equation 2.14,

                                                       (2.21)

where  

                                   (2.22)                          

The admittance  “looking in” to the left of the embouchure is

                                                (2.23)

where the reciprocal of  Equation 2.20 and 2.21 are the admittances of the embouchure and cork end.  The substitute tube replacing the left end as an open pipe of area =and length ?L is from,

.  This gives for the substitute end ?L,

                                               (2.24)

 

As shown in Figure 2.6, t=4.57mm, b=5mm, a=10.4mm, D=19mm, f=293.66Hz for .  We calculate the key chimney from Equation 2.22 obtain=10.9 mm.  For these short distances Equation 2.24 can simplify to

                                              (2.25)

Figure 2.8  The total open tube length for zero impedance on both ends is +.  The speed of sound, 344.7M/s, was calculated from Equations B14,B16 and B17 for a hot day F. 

 

To check on acutal values of the of Figure 2.8, as compared with the calculated value of 47mm; we take the head off the flute.  becomes 125 mm.  Closing the end with the hand and blowing the lowest note; it is 4972 Hz.  This is a closed end, quarter wave resonator of length.  173.4-125=48.5mm= compared to 47mm for the low note in Equation 2.25.  We can open the end so the note is at a measured frequency of about 945 Hz so.  The substitute open embouchure length is 182.4-125-0.3x20.8=51.2 mm. If we take the next highest mode of the closed head joint, the for a measured frequency of 1475 Hz.  The is 175.3-125=50.3mm.   These notes are still of low enough frequency that the simply form, Equation 2.25 is valid.  It is a reasonable estimate that the effective would be a little more than the calculated value 47mm because of the lip coverage that changes as shown in Figure 2.8.  Finally for the total length of the flute adjusted for , as shown in Figure 2.8, that .  This is very close to the average of the other three measurements, i.e. (48.5+51.5+50.3)/3=50.1 mm. 

All of this, of course gives us a false accuracy because the air jet distance is different for different players and a temperature change of a few degrees will swing the real reading by a millimeter or two.  Blowing loudly will move the note by 2 Hz and change the slide setting to change the measured by a few millimeters.  Also, Figure 2.9 shows changes in lip  position, and consequently, change in in Equation 2.24.

Figure 2.9  Measurements were made by Fletcher (1975) of the jet length from the lips to the back side of the embouchure, by experienced flute players for notes of difference pitch.

 

 

 

Next we determine the position of the lowest key hole.  We need to determine D in Figure 2.10.  To do this we have to solve for  with the branch equation.  is the ficticious pipe end that replaces the open key and end branch.

                                                     (2.26)

From Equation 2.3

and.                               (2.27)

Equation 2.26 and 2.27 together give,

                                          (2.28)

using the trigonometry identity from Appendix  D that  we get

and finally from the general form for and , we get,

    (2.29)

Figure 2.10  In determining the location of the key hold for the E note, it can be seen that the substitute open end for the key-end branch plus the difference in wavelength of the two notes, is equal to D+ used in calculating whole length. 

 

Equation 2.29 is very much simplified (Fletcher 1998) by using  instead of  if the distances are short enough.  Equation 2.30 is the simplified version.  Note that Equations 2.29 and 2.30 differ only by the exclusion of which was added because of the need to separate .

                                             (2.30)

For our thin walled flute example we use t=3mm.  From Equation 2.22 we get =8.3mm. =64 mm. The diameters 2a and 2b are 20.8 and 8.13 mm for the key hole and pipe respectively.  These numbers into equation 2.30 give.  To check this answer note that from Figure 2.10.  My measured  is 94.8mm and d=.3x2a=6.04mm.  This gives =100.8mm.  calculated from Equation 2.30.  At higher frequencies up the scale, you may want to use Equation 2.29 instead of 2.30.  For example, even at this low frequency,

and .404/6=.067M or 67mm.  This used for in Equation 2.30 instead of 64 mm; we get .  Now which is closer to the measured value of=100.8mm. can now be considered as the effective length of the pipe with one hole open.  The procedure can be repeated to locate where the next key hole, the  key hole should be placed after choosing a diameter 2b for the hole.  This is shown in Figure 2.11.

Figure 2.11  Knowing the location of the lowest key hole, calculating the distance D to the next key hole for used the same procedure used to determine the lowest key hole.

The second hole up is a larger hole with a diameter of 2b=11.43 mm.  With a wall thickness of t= 3 mm, Equation 2.15 gives =7mm.  2a is still 20.8 mm and  is 36.1 mm.  Using Equation 2.30, the substitute distance becomes,

mm. 

So mm, i.e..  I measure 40.6 mm.  If we use Equation 2.29 instead of Equation 2.30 the substitution distance is a little different, but not worthwhile for a craftsman to bother with.  At higher frequencies, the difference begins to be pronounced.  For example, when we get to , , , and , a consequential difference. 

When fabricating a flute “from scratch” the key holes are made about 0.5 mm smaller than the diameter used in calculations.  The holes are then slightly reamed for final tuning.  When making the final tuning, like a haircut, you can no put it back, but you can get tricks from experienced flute makers, such as tapering the key hole toward the bottom or top of the flute.  The next day it will be out of tune because the temperature and barometric pressure, and the loudness you are playing has all changed. Another change in tuning might be a swollen lip that resulted for a hearing problem that caused you to stand up when someone told you to shut up.   Play out of tune or design in a slide tuner.  Of course that person may have been a fiddle player who was playing out of tune.  Even more of a reason to have a slide tuner in your flute. 

 

 

 

 

 


 

 

Example 2 Silver Concert Flute

 

        First we establish the substitute   left of the embouchure as shown in Figure 2.12.  The impedance “looking into the left end” is the reciprocal of the reciprocal sum of the impedance of the two branches, the open embouchure pipe and the closed cork-end pipe as described by Equation 2.9 and Figure 2.3.  The impedance “looking into the left end” is also the impedance of the factious substitute end pipe.  This is the typical branch situation shown in Figure 2.3.

Figure 2.12.  The embouchure of a modern concert flute is rectangular.  Even more so when the player’s lip coverage is included.  Both for the embouchure height and for the open end substitute pipe are calculated values.  For completeness, is decreased by=-0.1 mm and so is the right hand pipe.

 

The impedance of the embouchure chimney, being an open pipe, is from Equation 2.3,

 

.                                                   (2.31)

Ditto for the substitute pipe,

                                                    (2.32)

The cork-end pipe is closed so the impedance is, according to Equation 2.7 and Figure 2.2,

                             (2.33)

So from Equation 2.9 we have for the substitute pipe,

                                                           (2.34)

which gives Equation 2.24 for the substitute pipe,

                                            (2.35)

is always in the 30 to 50 mm range so we cannot let unless we don’t mind being a few mm too large.  Otherwise D, k, and t have to be in meters.  Note that D are in the range.  This is why the position of the cork, D, does not come into play until the second or third octave and becomes large.  For  we get from Equation 2.15,

 

                (2.36)

and for the  note at 261.6 Hz that gives k=4.76, becomes,

                  (2.31)  

this gives radians and .  For completeness, we should take off .  But this gives an accuracy better than the precession, as Thom McGonnigal, an old NASA communications expert used to day.  Figure 2.13 shows the predictions for higher frequencies, when the cork distance D comes into play.  It is similar to a more detailed graph in (Fletcher & Rossing 1998 page 542-3 Figures 16.26 and 16.27) .  

 

Figure 2.13  Plot of the equation for the substitution embouchure length,  , as a function of cork distance to the center of the embouchure, D, area of the flute at the embouchure, , effective area of the embouchure, the wave number, k, and the substitute shunt chimney height of the embouchure,. 

 

We next determine the influence of closed holes on the length of the flute for the bottom note at 261.6 Hz. 

With /2 =345.1/2x261.63=659.9 mm; deducting the effective embouchure end effect of 40.2 mm and the other end effect of .3 (19)= 5.7 mm; the calculated embouchure-end distance is 614mm, as shown in Figure 2.8.  To check these results, we inserted the headjoint into a plastic pipe of 19mm internal diameter.  We then cut the pipe off when it sounded , as shown in Figure 2.13.  The measured length from the embouchure center to the end was 614 mm in agreement with the calculated value. The concert flute measures 608 mm; it is 6 mm shorter than calculated.  This is due to the series impedance of distance  in Equation 2.18 of a concert flute with 13 closed holes of chimney height t= 2 mm; 2b=13.5 mm average and  2a = 19 mm.  From Equation 2.18 we have,

=/key.          (2.37)

 

Figure 2.13  Bohem flute head with 17 mm cork distance connected to a 19 mm i.d. schedule 80 plastic pipe.  The distance from the center of the embouchure to the end is 614 mm for a note of 261.3 Hz. The same distance for a typical silver concert flute is 608mm.

 

Example 3: Rudell and Rose 1835 Wooden Flute

Prior to about 1850 flutes were commonly designed with a taper in the lower part as shown in Figures 2.13.  These flutes are now used for traditional Celtic music  and authentic orchestration of the old masters. 

A tapered pipe has an input  impedance, analogous to the input impedance of a cylinder as in Equation 2.2, of

                                          (2.38)

where  Figure 2.12 shows how these terms are defined for the impedance looking into either end.

Figure 2.12  The impedance looking into a conic section of length L uses the same equation, however the terms are different.  S1 is always the end into which the impedance is defined.  L=x2 – x1 is negative when looking into the large end.   

 

We first wish to see how much shorter a  tapered flute is than a cylinder flute sounding the same note.  To do this we note that a conic section open on both ends will have an impedance on the end of zero, just as the cylindrical one does.  Therefore, setting in Equation 2.38 gives,

.                                         (2.39)

A more usable and enlightening form of this equation is,

                                              (2.40)

We get to Equation 2.40 from 2.39 by using the identity ; using the definition ; and noting that .  Because we used the boundary condition that =0 at the end, we have to take in the end effect that the pressure is only zero at a small distance  beyond the end, as we did with the cylinder in Equation 2.3.  So we have to rewrite Equation 2.40 such that L=,

                                                                                                              (2.41)

 Note the similarity between this equation and the one for a cylinder, which we repeat here.

                                                                  (2.3)

Therefore, a cylindrical substitute tube would be

                                                                  (2.38)

We measure, then make an actual substitute tube, Figure 2.14.  The difference in measured length of the two flutes was 24.2 mm as shown in Figure 2.13.  The calculated difference is,

Figure 2.13    Dementions of the end of a boxwood flute similar to one made in 1830 by Rudell and Rose of London, model # 742, owned and played by Chris Norman The measured difference in length of a straight pipe and the Rudall & Rose 1830 Model 742 was 24.2 mm.  Measurements were to ½ mm accuracy.            

.   

Figure 2.14  A substitute tube used to compare the calculated substitute length for the Rudell and Rose flute.  The headjoint used was not the one from the Rudell and Rose flute, but one of the similar flute made by Michael Chronnally,  of M&E Flutes.  Similar heads were tried with similar results. 

 

The dimensions for calculating are shown in Figure 2.13 and the frequency is 261.6 Hz.  This gives for c=343M/s, k=2p(261.6)=4.79.  So for this situation Equation 2.38 becomes;  k= .829.  You have to be careful because dividing k by k to get .829/4.79 = .172M(172mm) is the wrong answer.  The reason is that has an ambiguity as shown in Figure 2.15.  We have to chose the option other than  =172 mm because we know  has to be in the 470 mm range.  So we try the other option of k=(.829-p)=2.31; =.483 M.  The calculated difference between and L is 483-457=26mm.  This is about 2 mm off from the measured value of 24.2 as shown in Figure 2.13.  There is an end effect difference of 2 to 3 mm because of the difference in diameter of the end of the two flutes.  Remember the end effect as being =0.61a.  a is the radius.

Figure 2.15.  Even though the tangent of these two angles is the same, they have different arguments kL resulting in an ambiguity in the equation when solving for L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 


 

 

Appendix A

Mathematics of Oscillatory Motion

 

Below is the presentation my grandson made to run for Treasurer of his Fourth Grade Class.  What he said about him and money, I hope you will be able to say about you and sound waves after you finish this section:

          I think that I would be a good treasurer for a number of reasons. First I would like to say that if you do in fact find that this speech has either therefore convinced or changed your mind that will be a good treasurer I will like to say before hand thank you for this chance to tell you my reasons and for taking the time to listen.

            I think you will find these reasons rather good and will agree with me and maybe others that I can be a very good treasurer. Well first of all this may sound conceded  in fact I can grantee it will sound conceded but I promise you I am not conceded. So first up I say that I rock at keeping track of money. If I were to tell you about the times my mom ‘lost my money’ and it was right were I put it    would take up the school day. Maybe i exaggerated but man it happens a lot!

            Now to go to two. Adding money, subtracting money, doubling money, for Petes sake dividing money! I’m reminding you i told you i would say some things that will sound conceded but this is the truth and i don’t want to sound braggy. I am in the GT math and think that i can do the adding and subtracting needed. I don’t even know if we will even need to do any of this. But we might so thats why you need to know that i am in GT and probably can handle it.

            Now for the third. Last year the third grade teachers picked pupils to go for SCO meetings. I was one of those four kids. I mean yeah that was last year it may seem weird but why do i say this if i know that? Why? Because i am going to tell you what they said and what they do. Not the SCO together but the treasurer.

            “We have three hundred twelve dollars and thirty-four cents from our last lollipop sale.” is what they’d say. Maybe how much they owed someone or how much they lost in the lollipop sale for paying people to make them. See? They do a lot but i can handle it. i don’t see much more to say so i guess i have one more thing.

            Thank you so much for taking your time to listen to me. I hope that you listened to me and think i would be the right person for the job. I say thank you and remember my motto wish i convincingly came up with on the bus Monday afternoon: sometimes the best of all is the oddest of all. Thank you.

 

            Well, for Pete’s sake we have to learn to add, subtract, multiply and divide sound waves. 

 

Addition of Two Waves Traveling In the Same Direction

            A wave is the propagation or transmission of a pressure crest and rarefaction in a restorable or elastic medium.  An example is the transmission of a hump along a stretched wire.  Another is a water crest traveling on the surface of a pond.  Another is an electric field traveling as a radio wave.  These are all examples of waves that travel perpendicular to the direction of the pressure crest.  For example, in a water wave, the wave moves along the surface; but the motion of the water is up and down.   The same for a wave traveling along a taught wire.  The wave moves along the wire but the motion of the wire is up and down.  The force is perpendicular to the wire.   Unlike water waves or waves along a wire, sound waves are a pressure crest (and the associated rarefaction) along the direction of the traveling wave.  In order to understand how musical instruments work and to modify them, one has to be knowledgeable about waves and all aspects there of, pressure, current, phase, traveling waves, standing waves, boundary conditions, wave excitation, coupling, dampening, and resonators that generate waves.  Resonators are clocks, musical instruments, and molecular vibration.  In a larger sense, all measurements are made by resonators or repetitive motion.   Everything reduces to mass, length, and time, according to Einstein.  You really can not think of time without a clock, or distance without repetition of  yard sticks or mile stones or wavelengths.  You can not think of mass without thinking of bunches of things in grams, pounds, or pecks.  So before we go into describing musical instruments, we have to go into the details of describing waves and resonators.

             The mathematical construct used to describe waves is the trigonometric function, the sine and cosine wave. 

 Figure A1  Above:  Sine function in degrees.  Below: Cosine function in degrees.  Sin 30o = Cos 60o = 1/2.  Note also that the sine of an angle is the  sine of the remainder when  dividing  by 360o.  For example the remainder of 810 when divided by 360 is 90 and so sin(810o) = sin(90o) = 1 and sin40o = sin400o = 0.642787609686539326322643409907263.  In engineering, one uses only the first three decimals, 0.642.  This is also the cosine of 50o.  

Either the cosine or the sine wave can be used to define waves, they are the same thing with a different starting situation, called the phase.  Trigonometric functions have to be in some angular measure.  Before computer days, engineering books all had tables of trigonometric functions as in Figure 2, for convenience. These tables were worked out my Christian monks in the Middle Ages.

 

 

Figure A2 Table of trigonometric functions from the back of a surveyor's field note book.  I still use surveyor’s field books as rugged journals with 50% rag paper. 

 

Trigonometric functions used for repetitive motion and wave study, use radians as the unit of angular measure of choice, not degrees.  Because the circumference of a circle with a radius r is r, and there are  radians of angular measure in a circle, the angle in radians is the ratio of the cord to the radius of a circle.   Figure 3 shows a sine wave in both degrees and radians. 

Figure 3.  This is a plot of signal strength ,s, as a function of angle, a.  It is written as .  The lower graph shows the angel in radians and the upper one is in degrees.  There are 3.141x2 radians in 360 degrees.  If radians are used as the angular measurement the  sine of the angle approximately equals the angle for small angles.  

 

            The waves in Figures 1 and 3 do not describe any thing real, they are just mathematical functions. 

Let's first describe a real wave, say a water wave (Figure 4).  This wave is frozen in time and so is the skier.

Figure A4.  A water wave of wavelength and crest height of HMAX with a water skier located at 1200 .

 

If we want to describe where the skier is in terms of distance, instead of angle in degrees, we can say that,

.  If we want to describe where the skier is in terms of distance, instead of angle in radians, we can say that,

where is the angle in radians, in this case,

  radians.

The wave of Figure 4 along a direction X with time frozen is,

            or or                  

,

where  is the wavelength, x is the location on the wave where the height H is being described with maximum for a wave of maximum height of  .  The phase angel  or  gives the value of the wave when X=0.  Note that the .  If you wanted to express a wave, stationary in time, and in degrees instead of radians, it would be,

                                                ,

which is never used for many reasons; one is that, for small angles, the angle in radians is approximately equal to the sine of the angle, as shown in Figure 5.

                        Figure A5.  As the angle  becomes small and become equal.

 

This approximation is very handy with flutes design.  For example, a keyed hole is a closed or open pipe of short length X, therefore is a small angle and the sine of the angle is replaced by the angle.    is commonly denoted as k, the wave number and, finally, a pressure wave in the direction x in a pipe, frozen in time, is usually written as,

,

where p is the pressure in the pipe, a distance x from some reference point and is uniform across the pipe (except close to the walls).  In the three dimensions, this equation takes on the vector form,

,

With flutes we use mostly the one directional form.

            Keep in mind that this describes a “still” picture of a wave that starts at Po  amplitude when x is zero.  How do we describe a moving wave, a wave moving in the negative X direction (to the left) in time?  First we note that in Figure 1, the top wave is the same as the bottom one except the bottom one has moved (to the left) backwards by 90o, or in radian measure it would be.  So the top one would be written as 

The bottom would be written as

.

Therefore, you can see that if you want to describe the wave moving backwards in time, you can just change the phase element biger and biger in time.   In other words, would be twice as far to the left as  is.  So if we change the phase angle to

,

where ? is, so far, just a proportionality constant indicating how fast the wave moves backward in time.    must have units of radians (of angle) per second because t must be angle.  There is another unit commonly used in place of radians per second.  There are 2 p radians of angle for each circle of 360o.  So for our angle measurement we used, not degrees and not radians, but the number of whole cycles per second the wave has gone through, the above equation becomes,

,                                                       (A 0 )

where f is our old friend, frequency.  In other words, if a wave goes backwards by 20 cycles per second, then it has gone through an angle of (20)• (2p) radians of angle.  is commonly called the angular frequency of a wave and f is the linear frequency, or just frequency.  When the  wave has traveled backwards a distance of one wavelength, , in a time of 1/f seconds per cycle, then the speed of travel backwards is distance divided by  time, or , the speed of travel of the wave is.  Just one more thing to complete the description of the backward moving wave, we need a term to describe what the value of the wave pressure is when the time and x are both zero; in other words, the all important phase of the wave.  So we just add a term  for the initial condition and now our final description of a backward moving is,

                                                                        ,                                             (A.1)

and of course a forward moving wave would be,

                                                                                                                      (A.2)

where  is the initial condition for the backward moving wave and  is the phase of the forward moving wave, i.e.,  and  are the values when t and x = 0. Another way of thinking about the wave being forward or backward traveling is to remember that must be constant for p to stay constant. In Figue A4 if the skier is to stay at the same place on the wave, i.e. if   is to stay constant; then kx - t can stay constant .  This means that kx and  getting positive together.  Conversely, for  it is just the opposite.  When time increases, x must be getting smaller for to stay constant.   A wave cannot be moving forward and backward at the same time.  So either equation (A.1) describes the wave or equation (A.2) describes the wave.  Always look for the + or – sign between the position description, x, and the time description, t.

            These are perfectly good description of forward and backward moving waves but they are not used because they are awkward and time consuming when adding, multiplying and dividing waves and otherwise manipulating wave functions.  It is also very difficult to keep track of the phase of two waves in their sine and cosine form.  And you will soon find out that in flute analysis, or analysis of any resonator, we will be obsessed with phase.  A more convenient and rather ingenious way to represent the wave is to use the definition,

                                                .                                             (A.3)

j acts like and has the properties of .  Of course there is no such thing as , but ;; and; just as  does; because = -1;; and = - . 

Figure A6. The graphical depiction of a wave in the imaginary representation.  The vertical axis is the imaginary axis and the horizontal axis is the real axis.

 

In this representation the backward traveling wave becomes,

                                  (A.4)

and the forward traveling wave becomes,

                (A.5)

We can see from Figure A6 that what pressure we measure is .  It is important to get used to this representation because it is used extensively in disciplines where waves are described.   A second advantage of this representation is that you can solve problems graphically, especially when you are adding, dividing, and multiplying waves and want approximate answers.  Also, we will see when we describe waves in flutes, the real part of the ratio of pressure to current (impedance) in a flute must vanish to define the resonance condition.  As we will soon see, when describing oscillators, the current and pressure have to be 90° apart, i.e., out of phase. 

 

To help understand this representation, we will add the two waves of Figure A1, i.e. the sine and cosine waves of the same frequency.  This is done graphically in Figure 7.  Then we do it mathematically with the complex representation.

Figue A7  Addition of a sine and cosine wave of equal frequency

 

In Figure A7 the distances from the base line for sine and cosine waves (shown in black) are just added to give the wave for the sum of the two.  Note that at 45o, 225o, and 405o the arrows show the sum of two waves at their equal magnitude to give the peek of the summed wave. 

So now let’s calculate what we did graphically in Figure 7.  The two waves in Figure 1 are

                                                                                                      (A.5)                                                                ,                                                       (A.6)

Where the sine has been changed into a cosine by changing the phase by - 90 o. We can now put them into a complex form by adding a j sin ? term.  Equation A5 becomes,

                                            (A.7)

and Equation A6 becomes,

                                                                                    (A.8)

From algebra we have, we can now separate out the phases in Equations A7 and A8 as

                                                                                                (A.9)

Now the two waves can be  added in the complex forms as Equations A8 and A9 as,

                                                     or,                     (A.10)

The first term can be substituted to be and so equation (A10) becomes,

                                               

We now draw a picture of the two phase components in Figure 8

       

Figure A8  A graphic depiction of Equation A10 at x=0.  As the angle kx increases, the vector rotates in the counter-clockwise direction.  If you stop the rotation at some value of X and let t increase, the victor  rotates in the clockwise direction.  Just the opposite is true for a backward traveling wave.

 

From the diagram of Figue A8 and Equation A10,  we can see that the description of the sum of a sine wave and a cosine wave of equal amplitude Po can be re written as,

,

or,

.                                 (A10) 

 

   In review, two waves with the same frequency (and of course wavelength), of amplitude  P1  and  P2 and phase  traveling backwards,

 

                                                                       (A.11)

or,

            ,               (A.12)

where is the vector sum or 1 and   2 as shown in Figure A8.         

Another way to look at it is to add projections on the two axes in the usual way of vectors as,

                                                                                                                (A.13)

and

                                                                                                  (A.14)

and,

As shown in Figure A8.

 Figure A8   Phaser P, ?, for the addition of two waves with phasers P1, ?1 and P2, ?2  

 

 

Note that if P1=P 2  and ? 1 = ?  2 ± 180o

                                                                                         (A.15)  

In other words, two waves of equal strength traveling in the same direction out of phase will cancel each other.  This never happens in the real world.  Do you care to venture why?  How about conservation of energy?  You would have two sources in the same place.  Two things that occupy the same space are the same thing.  Cancellation of beams is common, but you have to have them reinforcing some place else.  It is worth doing a little experiment just to show this important concept.  Take two common lasers.  Set them so their beams cross.  At the intersection put a piece of partially reflecting glass at 45o to both as is shown in Figure 9. This results in two beams of the same wave length traveling in the same direction.  If you make a careful adjustment of the beams, one set will vanish and the other will be twice as strong. 

Figure A9  Two lasers shining on a partially reflecting plate.  Low cost helium-neon lasers at 0632.8 nM work fine.  If one beam is weakened because the two halves are out of phase, the other one is reinforced, otherwise, conservation of energy would be violated.  You can also do the experiment with one laser and two mirrors.

 

            It is interesting to note that a helium-neon laser is a gas filled tube with a mirror on each end.  One mirror is partially reflecting where the beam comes out.  It is a miniature flute only the gas builds up the wave inside.  But like a flute there are modes generated harmonically at frequencies spaced by c/2L i.e., (/2 = L).  For example a 3 cm long laser at a wavelength of 632.8 nM (6.328 x 10-4 cm) and a speed of light at 3 x 1010 cm/sec; gives a spacing in frequency of the laser light of  c/2L = (3 x 1010 cm/sec )/2 x 3 cm = 5 G Hz separation between the laser modes.  There are about three or four spikes in frequency of laser light centered around 632.8 nM (orange-red) light.  In other words, the flute is ½ wavelength long or 6 cm wavelength.   Therefore one gets a beat of 5 G Hz between radiation in the same laser beam with a detector.  So in summation the adding of two waves of the same frequency going in the same direction will be a wave with an amplitude and phase determined by the vector sum of the two phasers.

            Speaking of conservation of energy, a wave traveling in a pipe will not go on indefinitely without attenuation.  The pipe warms up and so does the air.  The pipe warms up due to the friction of the air rubbing the pipe.  The air warms up because of its viscosity (thickness).  We will have lots to say about this later in Appendix C, but for completeness a plane wave moving forward in a pipe is described by:

,                                                            (A15b)

describes the diminishing of the wave power as it goes down the pipe a distance x .  Often we express the diminishing in time and not position by using where ct=x and is used just as often, for example is the 3db bandwidth of an oscillator, we will learn when we get into flute losses and sound levels.   We mention this because authors use for both and as described here. 

 

Standing Waves in a Flutes

           

            With the phase equal zero at t=0 and x=0, we can describe two waves of equal ampliltued, A,moving “through each other” as,

                                                            .                                    (A16)

or,

                                                                                                              (A.17)

in the complex representation.  To put this in a usable form in the triganometric representation, use the identities from Appendix D,

and Equation A.16 becomes,

                                                                        .                                               (A.18)

To put Equation A.17 in a usable form in the complex representation, note that,

Equation A.17 can be written as,

.

Use the identity  from Appendix D to get .

The real part of this is exactly Equation (A.18), the triganometric represetation. 

Figure A17.  Depiction of a standing wave resulting from two traveling waves of equal amplitude passing through each other.  The dotted lines show the standing wave at different times, ,, and  .  Note that this pressure standing wave between two null points is what we have in a flute.

            The amplitude of the forward traveling wave and the backward traveling wave are seldom of equal magnatude due to attinuation and leakage from the ends.  So for the more realistic case, we rewrite Equation A.17 replacing the amplitude of the backward traveling wave from A to B, to wit,

                                                            .                                                        (A.19)

We put the real and imagionary parts together to form,

                                                .                                          (A.20) 

Compare the part of Equation (A.20) with the equation depicted in complex space in Figure A7.  We show it in Figure 18. 

                        Figure A18.  Shown is a depiction of Equation (A.20) in complex space.  Note

                        that a coeffecient of  ej?t rotates counterclockwise at an angular rate ?t .

 

When the reflected wave amplitude is smaller than the incident wave, then the standing wave pattern shown in Figure A17 changes to the one in Figure A19 as a depiction of Equation (A.20), or

                                                .                                   (A21) 

                        Figure A19. Depiction of a standing wave resulting from two traveling waves

                        of unequal amplitude passing through each other.  The dotted lines show the sand-

                        ing wave at different times, t1,t2, and t3 .  This is similar the the situation where the

                        two waves are of equal ampliltued, as shown in Figure A17. 

 

 

Curent Impedence and Admittance

            As with the design of any oscillator, it is the ratio of pressure to current that we have to deal with because their ratio must be imaginary, i.e. out of phase, and the energy switches from one to the other at the resonant frequency, as explained above.  We get the relation between pressure and current from Newtan’s law, the force equals the mass times the acceleration.  Figure 20   shows our arrangement for considering this.

               Figure A20.    S is the area of the flute and ? is 1.21 Kgm/M3 and force is in Newtons.  It is pressure times area.

 

 

From  we get ,gives the all important,

                                                                        .                                                               (A22)

c is called the impedance of the medium in which the sound wave is traveling.  p is in Newtons per square meter and U is the “speed up” caused by the pressure p in meters per second.  .  For air at room conditions, ?= 1.2 Kgm/M3 (about 2 lbs. per cubic yard, heavy stuff) and c= 343 M/sec,giving an impedance of 415 Kgm/M sec.  For flutes of cross-sectional area S, we want to use the characteristic impedance of the flute which would be, for a one-way traveling wave; i.e. the pressure and current are in phase.

                                                                                                                 (A.23)

The units of U are ?x/?t, the speed of the sliver of Figure A20.  The units of u = ?x S/?t is the speed of the volume unit ?V=S ?x.  Some chose to refer to the characteristic impedance of a pipe as Zo =?c/S but I chose to use Zo = ?c= 415 Kgm/M sec and keep ?c/S as a symbol to emphasize the importance of pipe area, especially for those designing conical bore flutes.  If we go through the same process for a backward traveling wave, the current is out of phase with the pressure, to give,

                                                                                                              (A.24)

            Finally, we can get to the impedance of a standing wave in a pipe.  Equation (A.19) is what we used for the standing wave pressure in a pipe with a forward wave of amplitude A and a weakened reflected wave amplitude B, we repeat,

                                                            .                                                        (A.19)

The impedance of the standing wave will be,

                                                                               (A.25)

or from Equation A.19 and A.24, the impedance becomes                                                                                                                                            .                                                      (A.26)

If we write this equation for x=0, it is commonly labeled because it is the impedance that you “see” when “looking into a pipe from x=0 and it is,

                                                                        .                                                     (A.27)

For clarity and convention, we will relabel Equation (A.26) replacing x = L, indicating that is the impedance of a place down the tube a distance L.  So Equation (A.26) is rewritten as,

                                                            .                                                    (A.26)

If we change Equations A 26 and A.27 into the form,

                                                                        .

This allows for the elimination of B/A and gives an impedance for a pipe that is independent of a signal strength.  In other words, it gives an expression for the impedance “looking” into the pipe that is dependent on the pipe geometry alone.  This is the, all important equation,

                                                .                                                  (A.28)

This is the equation used for transmission of power on high voltage transmission lines at a 60 Hz frequency (5000 Km wavelength).  It is the equation used for sending and receiving signals from an antenna along a coaxial cable; and it is the equation for sending sound waves along a pipe.   The early work on sound in a pipe done by Helmhertz and followed up by others was not expressed in the form of Equation A.28.  It was not until formalization of electronics and the work of Lawrence Kinsler and Austin Frey(1942); Neville Fletcher; Thomas Rossing(1998); Arther Benade; and John Coltman and others, using that formalism; that the engineering was expressed in an intuitively obvious way.  The resonance frequency  occurs at  the value of kL, for which the real part of Equation A.28 vanishes.  For the flute we do not need to bother going through separating out the real and imaginary parts because one can see two cases where this occurs. The first is the resonance condition for clarinets, saxophones, bagpipes, and most woodwinds with reeds or lips closing on end.  When one end of the pipe is closed thenbecause the current has to be zero.  Rewriting Equation A.28 to see what vanishes, and temperedly neglecting the loss part,

                                               

                                                                                     

                                                                                                                 (A.29)

you can see that when ,

                                                                                                          (A.30)

       Our second condition, the one for flutes, that provides for to be imaginary, i.e., for the real part of Equation (A.28) to disappear is for the pipe to be open on the end, giving .  In this case  of Equation (A.28) becomes

                                                                                                           (A.33) which becomes the classical impedance for a flute open on both ends of effective length .  when or any multiples of half wavelengths. All other woodwinds are constrained by the boundary condition or Equation A30.  In this case when and or . 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 


Appendix B

Losses and Sound Output

The Speed of Sound

The two single names that loom first in contribution to physics are Newton and Clark Maxwell; then a close third and forth are Albert Einstein and Schrödinger.  Maxwell noses out the other two because all electricity and magnetism are described by his four famous equations.  The tedious measurements of Joule and Avogadro in the first half of the 19th century brought on the kinetic theory of gasses.  Avogadro, an exocentric Italian, could not sleep at night so, instead of counting sheep, he counted the number of molecules in a box.  His box was 22.4 liters and he counted  molecules.  The weight of the molecules in the box was always the same as the molecular weight of the gas.  In other words, when he had molecules of helium, the box weighed 4 grams.  When he had water vapor in the box, it weighed 18 grams.  Avagadro earned a lot of money by going around to strangers and saying, “I bet 200 farthings that I know how many molecules are in this box.”  The trick was to put one mole of gas in the box.  He became rich that way.  The reason he knew this was because of Maxwill’s 4 famous postulates (no less famous that his 4 postulates describing electricity and magnetism).  In any enclosed gas at a fixed temperature the molecules are moving around in a purely random manor.

  • The molecules are very small compared with the space between them.
  • All molecules are equally spaced, on the average.
  • The molecules exert no force on each other, except during collision.
  • The impact of the molecules when they do collide is perfectly elastic. 

Amazing!  This simple model, along with Newton’s laws took all the mystery out of the discipline of thermal and sound properties of gasses.  We start our description of the speed of sound by defining the capacity for gas to hold heat, the specific heat, where n is the number of moles, not molecules and is a small amount of heat in joules.  C is not a well defined quantity.  We can see why by considering the thermally insulated container shown in Figure B1.  The container is designed to maintain a constant pressure by having a weighted piston on top.  The specific heat of a container of gas that cannot expand, , is less than that for a container that is allowed to expand and keep the pressure constant, .  The reason is that the expanding gas does mechanical work on the piston.  Conservation of energy dictates that,

.                                                    (B1)

This equation says that the total heat  goes partly into raising the temperature of the gas at constant volume, plus the energy that goes into the change in pressure and volume, mechanical work.  We rewrite this as

.                                         (B2)

If is a lot smaller than we can say that the change in due to pressure, plus the change in due to volume, is the total change.  Including this into Equation 2B gives,

                                (B3)

When not adding heat and , if we push the piston to change the volume by , then the temperature change in the container is

                                (B4)

Figure 2B summarizes these equations that result in the equation needed for the speed of sound, to wit,

                                                     (B5)

for air is measured to be,