# Whistle

For other uses, see Whistle (disambiguation).
A party whistle.
A metal pea whistle.

An aerodynamic whistle (or call) is a simple aerophone, an instrument which produces sound from a stream of gas, most commonly air. It may be mouth-operated, or powered by air pressure, steam, or other means. Whistles vary in size from a small slide whistle or nose flute type to a large multi-piped church organ.

Whistles have been around since early humans first carved out a gourd or branch and found they could make sound with it. In prehistoric Egypt, small shells were used as whistles.[1] Many present day wind instruments are inheritors of these early whistles. With the rise of more mechanical power, other forms of whistles have been developed.

One characteristic of a whistle is that it creates a pure, or nearly pure, tone. There are many ways to create pure tones, but we restrict the descriptions here to what are called aerodynamic whistles. Strictly speaking, they are fluid mechanical whistles since they occur in gases, such as air or steam, as well as in liquids, such as water. The only difference between them is the fluid density and the sound speed.

The word aerodynamic whistle is used here since it is in common use. The conversion of flow energy to sound comes from an interaction between a solid material and a fluid stream. The forces in some whistles are sufficient to set the solid material in motion. Classic examples are Aeolian tones that result in galloping power lines,[2] or the Tacoma Narrows Bridge (galloping Gertie).[3] Other examples are circular disks set into vibration.[4]

The whistles described in this article are in a subclass where only the fluid is in motion and there is no significant dependent motion of the interacting solid. Depending on the geometry, there are two basic types of whistles: those that generate sound through oscillations of fluid mass flow, and those that generate sound through oscillations of the force applied to the surrounding medium.

## Early police whistles

Carved whalebone whistle dated 1821. 8 cm long.

Joseph Hudson set up J Hudson & Co in Birmingham, UK in 1870. With his younger brother James, he designed the 'Acme City' brass whistle. This became the first referee whistle used at association football matches during the 1878–79 Football Association Cup match between Nottingham Forest and Sheffield. Prior to the introduction of the whistle, handkerchiefs were used by the umpires to signal to the players.[5]

A police whistle being blown

In 1883 he began experimenting with pea-whistle designs that could produce an intense sound that could grab attention from over a mile away. His invention was discovered by accident, when he accidentally dropped his violin and it shattered on the floor. Observing how the discordant sound of the breaking strings travelled (trill effect), Hudson had the idea to put a pea in the whistle.[6] Prior to this, whistles were much quieter, and were only thought of as musical instruments or toys for children. After observing the problems that local police were having with effectively communicating with rattles,[7][8] he realised that his whistle designs could be used as an effective aid to their work.[9]

Hudson demonstrated his whistle to Scotland Yard and was awarded his first contract in 1884. Both Ratchet rattles and whistles were used to call for back-up in areas where neighborhood beats overlapped, and following their success in London, the whistle was adopted by most police in the United Kingdom.

This police whistle monopoly gradually made Hudson the largest whistle manufacturer in the world, supplying police forces and other general services everywhere. His whistle is still used by many forces worldwide. His design, was improved as the 'Acme Thunderer', the first ever pea whistle, which remains the most used whistle in the world; for train guards, dog handlers and police officers.From the 1880s and 90s, J. Hudson & Co began facing greater competition, as other whistle manufacturing companies were established, including W. Dowler & Sons, J. Barrall, R. A. Walton, H. A. Ward and A. De Courcy & Co. In 1987, Ron Foxcroft released the Fox 40 pealess whistle, designed to replace the pea whistle and be more reliable.

## Typical whistle sources and uses

There are numerous sources and uses of aerodynamic whistles. They can be used for hunting, controlling dogs, as toys, for sports events, as musical instruments, for police or military use, as ship or lighthouse fog horns, as steam whistles on trains, for scheduling and emergency use in industrial settings such as mines, in fluidic circuits, in domestic settings such at the tea pot whistle, for reducing foam bubbles, and for human whistling. They also occur as accidental byproducts of fluid flow such as supersonic jets, cavity resonances, whistling telephone wires, and idling circular saws.

## Types

Wilson, et al.,[10] in their study of human whistling (see below), pointed out the importance of including the symmetry or asymmetry of the unstable flow in addition to the feedback classes listed below. Because of the close relationship of flow symmetry to the sound field generated, their concept was included here as part of the sound source description (monopole -symmetric and dipole - asymmetric).

### The Monopole Whistle

Point Monopole

Whistles that generate sound through fluctuations of mass flow across a boundary are called monopole-like sources. The figure on the right is an example of a small sphere whose volume is oscillating. Any arbitrary fixed boundary drawn around the sphere will show a net mass flow across it. If the sphere is small enough relative to the sound wavelength it is emitting, it can be called a point monopole. For this type of source, the sound is emitted radially so the sound field is the same in every direction and decays with the inverse square of the distance. Real whistles are only approximations to this idealized model. most have boundaries around them such as the Hole Tone described below. Nevertheless, much can be learned about whistles with the useful form of the sound power equation for the monopole. Using the definitions below, it can be expressed as

${\displaystyle {{W}_{m}}={\frac {\pi {{\rho }_{0}}}{{c}_{0}}}{\widehat {{Q}^{2}}}S{{t}^{2}}{{U}^{4}}{{L}^{2}}}$

The variables U and L are considered characteristic of the source and their correct choice is important.

### The Dipole Whistle

Point Dipole

Whistles that generate sound through fluctuations of momentum or a force exerted on the surrounding medium are called dipole-like sources. The figure on the right is example of a small rigid sphere that is moving back and forth in a given direction. If the sphere is small with respect to the wavelength of the sound emitted, it can be called a point dipole. A force must be applied to the sphere in a specific direction to move it. The surrounding medium in the direction of motion is compressed to radiate sound, but the medium at right angles slides past the sphere and is not compressed.

This results in a non-uniform sound field, unlike the monopole whistle. Real whistles are only approximations to this idealized model. Although tuning forks are not whistles, they create sound fields that are very close to the idealized dipole model. Nevertheless, much can be learned about whistles with the useful form of the sound power equation for the dipole. Using the definitions below, it can be expressed as

${\displaystyle {{W}_{d}}={\frac {3\pi \rho {}_{0}}{c_{0}^{3}}}{{\widehat {F}}^{2}}S{{t}^{2}}{{U}^{6}}{{L}^{2}}}$

Once again, U and L must be chosen correctly.

## Feedback categories

Aerodynamic whistles rely on the instability of flows and some mechanism downstream to send the disturbance back to the origin to continue the instability. There are several ways that feedback can occur.[11]

### Category I

The sound from a Category I whistle is primarily a byproduct of source motion.[11] In every case, there is a back reaction of the medium on the source (resistive and reactive impedance). One example of a weak back reaction is a vibrating iron body in air. The densities are so different the back reaction is often ignored. Back reactions of air on an air source or water on a water source can be different. In many cases, say turbulent jets, the sound created is random and it is convenient to consider the sound to be merely a byproduct of the flow. In this category the back reaction is insufficient to strongly control source motion, so whistles are not in this category.

### Category II

The back reaction of the medium is an determinant of source motion.[11] In many important cases, linear thinking (small cause = small effect) is fallacious. Unstable fluid motion or the sound generated by it can feedback to the source and control it, much like audio feedback squeal. The basic requirements for a feedback controlled system are: (1) a source of steady power; (2) an amplification mechanism that can convert the steady power to time varying power; (3) a disturbance which supplies the oscillations to be amplified; (4) a means of generating sound or other oscillatory fluid motion; and (5) a means for feedback of that oscillatory motion as a disturbance to the input of the amplifier. Whistles are in this category. There are several ways to describe the feedback process.

#### Class 1

The feedback is essentially incompressible; the speed of sound, although finite, is sufficiently large that it can be considered infinite.[11] This action may be called near field or hydrodynamic feedback. There are a number Class I devices. The feedback that causes a stick in a water stream to vibrate, or a flag to wave, is due to hydrodynamic feedback.

#### Class II

Whistle Feedback Loop

The feedback is compressible and depends on the speed of sound.[11] This action may be called intermediate field, quasi-compressible feedback. A well known example is the hole tone (described below) where the feedback distance of a compressible (sound) wave is very small as compared to the wavelength of the sound.

#### Class III

The feedback is compressible and depends on the speed of sound.[11] This may be called far field or acoustic feedback. the feedback distance of a compressible wave can be an appreciable fraction of the wavelength of the sound. An example is the flute.

The figure on the right shows a block diagram of these feedback mechanisms. All aerodynamic whistles operate under one of the classes.

## Stages

Nonlinear feedback stages

There are feedback loops associated with many whistle operations and they are non-linear.[12] Because of the nonlinearity, there can be several possible conditions for a given flow speed or geometry. Which of these is dominant is determined by the gain of the unstable flow at a particular frequency and whether the feedback is constructive or destructive.

Early studies have used the term stage to describe the possible states of feedback as shown schematically in the figure on the right. As the flow speed increases (Reynolds number - Re) the frequency slowly climbs (nearly constant Strouhal number -St) but then the frequency jumps up abruptly to a higher stage. As the flow speed is later decreased the frequency slowly decreases but then jumps down abruptly to a lower stage. This pattern is called a hysteresis loop.

At any particular flow speed it is possible for one of several loops to be dominant, depending on how that speed was achieved. In a number of the whistles described here, Stage I is associated with the development of a single vortex in the distance between initiation of the flow instability and initiation of the feedback signal.

Higher stages are associated with more vortices in that distance, hinting that this distance can be an important characteristic dimension. In several whistles, three stages have been identified (edge tone). By blowing hard in some musical wind instruments Stage I jumps to Stage II; it is called overblowing.

## Flow instability

Instability of laminar jets
Unstable water jet

Flow instability is the engine for whistles. It converts steady energy to time dependent energy. Conversion of laminar flow to turbulent flow is a well known example. Small disturbances to the laminar flow cause the transition.

An example is shown in the figure on the right with a water jet.[13] The laminar two-dimensional jet amplifies small disturbances at the orifice to generate a vortex street. For this case, the flow speed, in terms of Reynolds Number, was graphed against the disturbance frequency, in terms of the Strouhal Number for a variety of disturbance amplitudes to reveal the region of instability as shown in the figure on the left. The value of D in the figure represents the ratio of the lateral disturbance displacement to the nozzle width; the disturbances were minute.

The disturbance was temporal in the example, but can also be spatial. Transition to turbulence can occur over a rough surface or over an irregular shape, such as an aircraft spoiler. All whistle mechanisms described here are created by temporal disturbances that are in one of the three classes.

One important source of instability in a fluid is the presence of a velocity gradient or shear layer with an inflection point. Thus, fluid instability can be described as a three dimensional region with flow speed on one axis, disturbance amplitude on the second, and velocity profile on the third.

In a whistle, the instability starts at some point in the three dimensional region and then moves along some path in that region as the local variables change. This makes the comprehensive understanding of whistle instability mechanisms very difficult.

## Scaling

Whistles come in all shapes and sizes, but their operation can be unified through the concepts of dynamic and geometric similarity. Nature knows nothing of the specific measurement systems we use; it cares only about ratios between the various forces, time scales, and the several dimensions. To compare them we need to take into account the established ratios that are relevant to whistle operation.

Similarity is best exposed by determining a speed, U, that is characteristic of the dynamics, and a dimension, L, that is characteristic of the geometry. If these values are used in dimensionless numbers, such as those listed below, much understanding of the phenomenon can be achieved.

### Strouhal number

${\displaystyle St={\frac {f{{L}_{1}}}{U}},S{{t}_{a}}={\frac {f{{L}_{2}}}{{c}_{0}}}}$ The first number is the ratio of unsteady inertial forces to steady inertial forces. The number was named in honor of Vincenz Strouhal[14] who first deduced the relationship between the vortex shedding frequency around a cylinder and the flow speed. The characteristic variables were the cylinder diameter L1, and the speed of the flow over it U. He found the number to be reasonably constant over given a Reynolds number range.

This number permits relationships to be developed between the different sizes and speeds. Now the Strouhal Number can be derived directly from the dimensionless form of the mass continuity equation. This equation may be referred to as a fluid mechanical Strouhal number in comparison with the second version which may be referred to as the acoustical Strouhal number.

The first version is used for dynamic similarity of the fluid motion in whistles while the second version is used for dynamic similarity of the acoustical motion in whistles. Many whistles, especially those with Class III feedback, require use of both numbers (see (Monopole-Dipole Whistles). The acoustical Strouhal number is essentially the Helmholtz number with the ${\displaystyle 2\pi }$ deleted.

### Mach number

${\displaystyle M={\frac {U}{{c}_{0}}}}$ It is the ratio of the steady speed to the speed of sound. The number was named in honor of Ernst Mach[15] who first studied (among other things) supersonic flow and shock waves.

This number describes the range between flows that can be considered incompressible and flows where compressibility is important. Now the Mach number can be derived directly from the dimensionless form of the momentum equation.

### Reynolds number

${\displaystyle \operatorname {Re} ={\frac {UL}{\upsilon }}}$ It is The ratio of the steady inertial forces to the steady viscous forces.

The number was named in honor of Osborne Reynolds,[16] an engineer who did pioneering studies on the transition of laminar to turbulent flow in pipes.

Now the Reynolds number can be derived directly from the dimensionless form of the momentum equation.

### Rossby number

${\displaystyle \operatorname {R} o={\frac {U}{{{f}_{0}}L}}}$ It is the ratio of linear velocity to tangential velocity for swirl flows. The frequency is characteristic of the rotation rate of the flow.

The number was named in honor of Carl-Gustav Rossby,[17] a meteorologist who first described the large scale motions of the atmosphere in terms of fluid mechanics. He described the jet stream, and his number was first used to describe the motion associated with the Coriolis force in the atmosphere.

Now the Rossby number can be derived directly from the dimensionless form of the momentum equation in curvilinear coordinates.

### Dimensionless force

${\displaystyle {\widehat {F}}={\frac {F}{{{\rho }_{0}}{{U}^{2}}{{L}^{2}}}}}$ The ratio of the actual dynamic force to the steady momentum..

### Dimensionless volumetric flow rate

${\displaystyle {\widehat {Q}}={\frac {Q}{U{{L}^{2}}}}}$ The ratio of the dynamic volumetric flow rate to the steady volumetric flow rate.

## Monopole-like whistles

In these whistles, the flow instability is symmetric, often resulting in periodic ring vortices and the sound generation is associated with fluctuations of volumetric (mass) flow rates. The sound field is as close to the directivity of an actual monopole source as local geometry allows.

### Hole tone (tea pot whistle, bird call)

Hole Tone geometry
Hole Tone frequencies

The steady flow from a circular orifice can be converted to an oscillatory flow by adding a downstream plate with a circular hole aligned with the orifice. Small perturbations in the flow feedback to the orifice to cause a variable volumetric flow rate through the downstream hole because of the symmetry of the feedback.

The disturbance in the jet is a symmetric vortex ring that moves at some speed slower than the mean jet speed until it encounters the hole and some fluid is forced through it, resulting in a monopole-like sound field in the half space outside. The oscillatory volumetric flow in the hole sends a wave back to the orifice to complete the feedback loop and causing a nearly pure tone.

The figure on the right shows a schematic of the geometry.

To invoke dynamic similarity,[18] the characteristic speed in a study was chosen to be the average speed of the jet at the orifice U (deduced from measured volumetric flow rate) and the characteristic length was chosen to be the orifice diameter δ. Tests were made at five spacing distances h/δ from the orifice. Two scaling laws were used; the Strouhal number was graphed as a function of Reynolds number. The results are shown in the figure on the right.

The frequency of the tone determined by how often a vortex encountered the hole while moving at some speed u less than the initial jet speed. Since the jet slowed as it proceeded toward the hole, the speed of the vortex slowed with it so the frequency and Strouhal number was greater at closer spacing. The Strouhal number data showed clearly the almost linear relation between frequency and initial jet speed. The number would be more constant if the actual jet speed at the hole could have been used as the characteristic speed. In four of the distances tested, there were jumps between Stage I and Stage II. The hysteresis loops are clear indications of the complex nature of the jet instability gain structure.

The uniformity of the measured sound field for this whistle confirmed its monopole-like nature. Measurements of speed dependence of the sound level showed it to be very close to ${\displaystyle {{U}^{4}}}$, further confirming the monopole nature of the source. At these speeds and spacings, the feedback was normally Class II, but reflecting surfaces as far away as three meters and with proper phasing, controlled the tone, converting the feedback to Class III.

At higher Reynolds Numbers, the flow became chaotic resulting in broad-band sound. The hole tone has been rediscovered in the form of the tea pot whistle.[19] They found that above a Reynolds number of 2000, hole tone operation occurred with symmetric vortex evolution and a constant Strouhal Number with Reynolds number. Comparison of their data with data in the figure suggests that the cylindrical enclosure between the two orifices raises the Strouhal number. There was no mention of frequency jumps. They noted that at lower speeds, the cylindrical volume responded as a Helmholtz resonator. Baron Rayleigh[20] was aware of this whistle; it was called the bird call then.

There seems to be evidence that events similar to the hole tone occur on aircraft landing gear covers with circular holes. In Australia, there is the Tenterfield Fox Whistle[21] and the Traditional Fox Whistle that appear to operate as hole tones.

### Corrugated pipe whistle

Corrugated pipe whistle

This whistle has dozens of popular names. Pipes with sinusoidal variations of radius are often created to permit bending. Steady flow through the pipe at low Reynolds numbers results in a fluctuating volumetric flow rate that generates a monopole-like sound field at the pipe exit. Examples of such pipes are shown in figure on the right.

The yellow plastic pipe is actually a child’s toy that sounds when the pipe is whirled around. The metal pipe shown was used in the Concorde cockpit to provide cooling air to the pilots, but its loud tone got it canceled. This whistle is similar in many respects to the hole tone, in particular the teapot whistle. It is subject to frequency jumps and hysteresis loops. There are numerous articles on the internet about this whistle and it has been studied.[22][23][24]

The characteristic speed is the mean flow through the pipe U and the characteristic length must be a multiple of the spacing between corrugations, nL, where n is an integer number and L is the distance between corrugations. At low speeds, the unstable interior flow needs to travel several corrugations to establish the feedback loop. As the speed increases, the loop can be established with fewer corrugations. Simple tests were performed on the yellow plastic tube.

The Strouhal number, ${\displaystyle {{S}_{t}}={\frac {{{f}_{n}}nL}{U}}}$ was used as the scaling factor. The highest frequency (7554 Hz) was found in the “overblown” condition and n was presumed as one corrugation. At the least flow rate, the frequency of 2452 Hz compared favorably to n=3. At intermediate flow rates, several non-harmonically related frequencies occurred simultaneously suggesting that several corrugations were involved in the sound generation. In the smaller metal tube, a predominant tone appeared at 6174 Hz and corresponded to n=2. A unique aspect of this whistle is that the internal flow carries both the unstable vortex downstream and the returning feedback signal upstream.

### Pipe tone (Pfeifenton)

An acoustical diode

The unique feature of this whistle is that the tone sounds only with flow through the orifice from outside; it is an acoustical diode. A cylindrical cavity with a small circular, square-edged hole at one end and totally open at the other is known to generate a tone when air is passed through it. It is subject to frequency jumps and hysteresis loops similar to the hole tone. There appear to be two stages and the feedback is likely Class II if the tube is lshort. The fundamental tone occurs near ${\displaystyle \lambda =4L}$, so one characteristic dimension is L, the length of the tube. The characteristic speed U is that of the flow through the hole.

A monopole-like sound field is generated by volumetric flow rate oscillations. Karthik[25] and Anderson[26][27][28] have studied this phenomenon and concluded that symmetric vortex shedding on the cavity side is the driving agency.

An example of this device is shown in the figure on the right; it had a 0.125 inch diameter hole, was 1.9 inches long, and 0.8 inches in diameter. The quarter wave resonance was calculated to be 1780 Hz, while the measured fundamental was 1625 Hz with detectable second and third harmonics. End corrections for radiation from the openings is needed to bring the two frequencies in consonance. To determine the end corrections, two additional dimensions are needed: the diameter of the orifice d1 and the diameter of the tube d2.

### Hartmann, Galton whistles (stem jet)

The Hartmann whistle

The previous whistles occur at low flow speeds, this whistle occurs at very high speeds. When a subsonic jet impinges on a cavity, jet instability becomes part of the feedback loop as with the hole tone. When a supersonic jet impinges on a cavity, shock wave instability becomes part of the feedback loop. The figure on the right is one example of this whistle. A cylindrical cavity with one end open and facing the supersonic circular jet will result in extremely intense sound. The shapes in the figure represent the shock/expansion cells within the jet. A related configuration, called the stem jet, has a central rod in the jet that extends to support and align the cavity. There are a number of other geometric variations, all of which operate in similar fashion.[29]

These devices have been studied,[30] and reviewed by Raman.[31] Here we look primarily at the Hartmann whistle. The shock cells of the jet interact with the shock in front of the cavity (the flow in the cavity being subsonic). Small symmetric disturbances in the jet stream are amplified as they proceed toward the cavity (similar in some respects to the hole tone) causing the shock in front of the cavity to oscillate. The shock front acts much like a piston source of high energy resulting in a monopole-like sound field. Again the volumetric flow is directional unlike the theoretical monopole.

The sound field may be similar to that created by oscillatory flow from a pipe, except for presence of the supersonic jet structure which can strongly modify the directivity. The original equation of Hartmann is shown below.

{\displaystyle {\begin{aligned}&{\frac {\lambda }{d}}=5.8+2.5\left\{{\frac {h}{d}}-\left(1+0.0041{{\left(P-0.9\right)}^{2}}\right)\right\}\\&{{S}_{t}}={\frac {fd}{{c}_{0}}}\approx 0.17\approx {\frac {fh}{U}}\\\end{aligned}}}

The diameter of the orifice and cavity is d, the distance between orifice and cavity is h, and the orifice pressure P was given in kilograms per meter squared. At the lower limit of h the second term disappears. In this case, the equation could have been reformatted in terms of the acoustical Strouhal Number as shown in the second equation above. The characteristic speed U at the nozzle is the sound speed ${\displaystyle {{c}_{0}}}$. It is interesting that the number is very close to that found by Strouhal for flow over a cylinder. There are two characteristic length scales. Nozzle diameter d characterizes the sound power while the separation distance h characterizes the frequency.

Comprehensive studies of this phenomenon[32][33] have shown that the position of the cavity is critical in creating sound. The process has hysteresis loops and the frequencies are related to multiples of the quarter wavelength resonance of the cavity. After reformatting Hartmann’s formula, and using the new formulation above, an equation for sound power can be written.

{\displaystyle {\begin{aligned}&{{W}_{h}}={\frac {\rho }{2}}{\frac {\pi {{d}^{2}}}{4}}{{\left(2\pi fa\right)}^{2}}{{c}_{0}}=A\rho {{f}^{2}}{{d}^{2}}{{h}^{2}}{{c}_{0}}=A{\frac {\rho }{{c}_{0}}}{{\left({\frac {fd}{{c}_{0}}}\right)}^{2}}c_{0}^{4}{{a}^{2}}\\&{{W}_{h}}=A{\frac {{\rho }_{0}}{{c}_{0}}}S{{t}^{2}}{{U}^{4}}{{L}^{2}}\\\end{aligned}}}

The Galton whistle

Since the characteristic speed U and sound speed are essentially the same, it can be rewritten as the second equation. This equation has the same structure as the one for the point monopole shown above. Although the amplitude factor A replaces the dimensionless volumetric flow rate in these equations, the speed dependence strongly confirms the monopole-like characteristics of the Hartmann whistle. A cousin to the Hartmann whistle is shown in the figure on the right, the Galton whistle. Here the cavity is excited by an annular jet which oscillates symmetrically around the sharp edges of the cavity. It appears to be a circular version of the edge tone (discussed below) in which the symmetry of the otherwise dipole source of the edge tone is converted to a monopole source.

Since it is highly likely that the oscillations are coherent around the periphery, there should be a fluctuating volumetric flow rate from the cavity with only a small net lateral force. Thus the source is yet another version of a monopole-like geometry; the volumetric flow rate is a cylindrical area between the jet and cavity.

### Rijke tube

The Rijke tube

There are a number of whistle phenomena in which heat plays a role. The temperature in a sound wave varies, but since it is so small normally it is common to neglect its effects. When amplification can occur a small variation can grow and have important influence on the sound field created. The most well known thermal whistle is the Rijke Tube.[34] Peter Rijke placed a heated a gauze material inside a vertical tube.

Originally, the gauze was heated with a Bunsen burner; later, a wire grid was heated electrically. The heat transferred to the air in the tube sets it into near half-wave resonance if the gauze is placed below the midpoint of the tube as shown in the figure on the right. There is no theoretically optimum position, as the wave speed upward is the sound speed ${\displaystyle {{c}_{0}}}$ plus u, the convection speed, while the wave speed downward is ${\displaystyle {{c}_{0}}}$ minus u. Without a convection flow, the midpoint and the lower tube end are the best locations for heat transfer. With convection, a compromise position halfway between the two points is normally chosen which depends on the amount of heat added. One characteristic length associated with frequency is tube length L.

Another characteristic length associated with sound power is αL, the position of the gauze. The characteristic speed must be the convection speed u at the heat source. For detailed study of the whistle, see Matveev.[35] Since the first mode resonance is about half-wave, the sound field emitted from the tube is from two in-phase monopole-like sources, one at either end. A gas flame inside a tube can drive resonance; it was called a singing flame. There is a reverse Rijke Tube when hot air passes through a cold grid.

### Sondhauss and Taconis tubes

The Sondhauss Tube is one of the early thermal tone generators; it was discovered in the glass blowing industry. A bulb with hot air is connected to one end of a tube that is at room temperature. When the cold tube is blown, tube acoustic oscillations occur. It was discussed by Baron Rayleigh in his Theory of Sound. This device is not considered a true whistle since the oscillations decay as the temperatures equalize.

In analyzing this tube, Rayleigh noted that if heat had been added at the point of highest density in the sound wave, and subtracted at the point of lowest density, the vibration would be encouraged. Another thermal effect is called the Taconis oscillation.[36][37] If a stainless steel tube has one side at room temperature and the other side in contact with liquid helium, spontaneous acoustic oscillations are observed. Again, the Sondhauss tube is not a true whistle.

### Human whistle

Human whistle sound

The number and variety of whistles created by humans is quite large, yet very little has been done to examine the physics of the process. There are three possible mechanisms: Helmholtz resonance, symmetric hole tone operation (monopole), or asymmetric edge tone operation (dipole).

Wilson and his colleagues[38] have simulated the human whistle by creating a cylinder 2.04 inches in diameter with a rounded orifice at one end that supplied a jet and another rounded orifice at the other end of the same diameter and on the same axis. The geometry was very similar to that of the tea pot whistle. After a number of tests at various speeds, orifice diameters, and orifice thicknesses, they concluded that the whistle was created by a Helmholtz resonance in the cylinder volume. There was enough data for one case in their study to calculate the Strouhal and Reynolds numbers. The results are shown in the figure on the right.

The Strouhal number was essentially constant over the limited speed range, suggesting hole tone operation with Class I or Class II feedback. Their work indicated symmetric unstable vortex flow as would be expected but there was no mention of stages. In the study by Henrywood, it was noted that Helmholtz resonance could occur at low speeds. The flexibility of the mouth suggests that although a hole tone feedback mechanism is highly likely, the possibility of Helmholtz resonances in the mouth cavity and asymmetric edge tone actions with the teeth are considered possible.

## Dipole-like whistles

In these whistles, the flow instability is asymmetric, often resulting in rows of alternate vortices and the sound generation is associated with fluctuations of applied force. The sound field is as close to a dipole source as local geometry allows.

### Aeolian tone

Cylinder vortex street
Aeolian Tone

The steady flow over a cylinder (or similar object) will generate vortex shedding and consequent sound. The early Greeks used this phenomenon to develop a harp and the sound was called an Aeolian tone after Aeolus, God of the Wind.

Whistling telephone wires, automobile radio antennae, certain automobile front grilles, and smoke stacks are other examples of this tone. At very low Reynolds numbers, the flow around a cylinder is stable, forming two fixed vortices behind it. As the speed increases, the flow, although laminar, becomes unstable and vortices are shed alternately.

Hydrodynamic feedback (Class I) influences the formation of new vortices and exerts a fluctuating force on the cylinder. The flow field is shown in the upper figure on the right (created by Gary Koopman). Theodor von Karman[39] identified and analyzed the flow behind objects like a cylinder and since then this special flow has been called the Karman vortex street. Vincenz Strouhal was the first to scientifically investigate the sound emitted by flow around a rigid cylinder. At low Reynolds numbers the tone was pure and the frequency was proportional to the steady flow speed U and inversely proportional to the cylinder diameter d.

For many applications, the first equation below is often used. A review of the literature[40] produced the figure on the right for the Strouhal number. At low Reynolds numbers the Strouhal number rises as inertial effects begin to dominate and then decays slightly at higher numbers. The second equation below is a best fit for the data between 1000<Re<100000.

{\displaystyle {\begin{aligned}&St={\frac {fd}{U}}\approx 0.2\\&St=.218{{\left({\frac {\operatorname {Re} }{1000}}\right)}^{-0.042}}\\\end{aligned}}}

Cylinder vortex development

It is surprising how often oscillatory flow phenomenon have Strouhal Numbers in this range. For shape comparison, the Strouhal number for an ellipse has been measured at 0.218, a cylinder at 0.188, a square at 0.160, and a triangle at 0.214. The characteristic dimension is that of the object lateral to the flow and the characteristic speed is that of the impinging flow.

The second equation suggests that the Strouhal Number is a weak negative function of Reynolds number. This suggests that the dynamic similarity approximation is reasonable. The fluctuating force exerted on the cylinder is a result of the flow circulation around it caused by the alternate vortex separation as suggested in the third figure. The fact that the vortices are not directly behind the cylinder suggests that the force vector has both a lift and drag component resulting in lift and drag dipoles.

An approximate way to relate the sound generated to the flow characteristics is to perturb the standard drag equation with velocity perturbations as shown in the upper equation below.(lift measurements for cylinders are generally not available). The upper equation is the modified drag equation with both drag component u and lift component v and the cross sectional area dL where d is the cylinder diameter and w is the length.

{\displaystyle {\begin{aligned}&F={{F}_{d}}+F_{d}^{'}+F_{l}^{'}={\frac {{{C}_{d}}{{\rho }_{0}}}{2}}{{\left(U+u'+v'\right)}^{2}}dw\\&{{W}_{d}}={\frac {3\pi {{\rho }_{0}}}{c_{0}^{3}}}{{S}^{2}}{{U}^{6}}dwC_{d}^{2}{\overline {{\left({\frac {{u}^{'}}{U}}\right)}^{2}}}\\&{{W}_{l}}={\frac {3\pi {{\rho }_{0}}}{c_{0}^{3}}}{{S}^{2}}{{U}^{6}}dwC_{d}^{2}{\overline {{\left({\frac {{v}^{'}}{U}}\right)}^{2}}}\\\end{aligned}}}

Aeolian Tone

Manipulation of the equation yields the lower two equations for the dipole sound power of both lift and drag. Each time a vortex is shed, the drag velocity fluctuation u has the same sign, but the lateral velocity fluctuation v, has opposite signs, since the vortex is shed on alternate sides. As a result, The drag dipole would be expected to be twice the frequency of the lift dipole.

Phillips[41] found the lateral velocity fluctuations were two orders of magnitude greater than the longitudinal, so the lift dipole is 20 dB above the drag dipole. He found the drag dipole did occur at twice the frequency of the lift dipole. At higher speeds, the vortex separation may not be correlated over the entire length of the cylinder resulting in multiple essentially independent dipole sources and lower sound power. The lower figure on the right shows the correlation coefficient as a function of distance along the cylinder and is from the Etkin, et al. study.

### Vortex meters

Vortex flow meter

Is there a use of Aeolian tone knowledge other than making musical instruments? There are several accurate flow meters that are designed to take advantage of the constancy of the Strouhal Number with Reynolds Number to provide a linear relationship between flow rate and measured frequency. They are called vortex meters.

A particular shaped object is placed within a pipe and a pressure sensor is embedded either in the pipe wall or in the inserted object. Although a number of shapes have been used, there are several that work well. The figure on the right shows data from a meter called the Deltameter.[42] The inserted shape was that of a trapezoid with the wider end facing upstream. The graph shows a nearly linear relationship with flow speed over nearly a 1000 to 1 range of Reynolds Numbers (12,530 to 1,181,000). Orifice plates as flow meters typically have a 5 to 1 speed range, while turbine meters may have up to 100 to 1 range.

It should be noted that three-dimensional (viscous) effects occur at low Reynolds numbers so dynamic similarity is not achieved there. At higher speeds the Strouhal Number is close to that for the Aeolian tone. The dependence of Strouhal Number on Reynolds Number for this confined geometry is slightly negative as was found also for the Aeolian tone.

### Mountain vortex tone

Mountain vortex street

Photographs taken from space have shown alternate arrangements of clouds around mountains; The figure on the right shows an example. Does this type of event create periodic sound? The NOAA laboratories in Boulder, Colorado, were tasked with detecting the extremely low frequency sound of nuclear tests. They detected one and by triangulation determined it was occurring in the Aleutian chain. As the figure suggests, it was vortex shedding from a volcanic cone. There are numerous satellite photos available on the web showing this phenomena in many places of the world.

These whistles generate enough sound to be detected but at frequencies below 1 Hz they are inaudible. Like the Aeolian tone, the feedback is Class I. Using Strouhal's number, it might be possible to estimate the wind velocities; they appear to be quite high.

### Trailing edge tone

Glider sound
Trailing edge tone
Trailing edge spectra

The boundary layer on the airfoil of a glider is laminar and vortex shedding similar to that of a cylinder occurs at the trailing edge. The sound can be a nearly pure tone.

The figure on the left shows a one-third octave band spectrum taken under a glider flyover; the tone is 15 dB above the broad band sound. The aircraft speed U was 51 m/s (170 ft/s) and the frequency was near 1400 Hz.

Based on a Strouhal Number of 0.20, the characteristic dimension δ was calculated to be near ¼ inch; the boundary layer thickness. A dipole sound field was created at the trailing edge due to the fluctuating force exerted on it.

At higher speeds on powered aircraft, the boundary layer on the airfoil is turbulent and more complex vortex shedding patterns have been observed. Since it is difficult to measure in flight, Hayden[43] made static tests.

The figure on the right shows an example. A boundary layer flow was created on both sides of a thin rigid flat plate which terminated with a square trailing edge. Note the nearly pure tone at 2000 Hz with a Strouhal number of 0.21 protruding above the turbulent sound spectrum. Once again the magic number of Strouhal appears. The characteristic speed was the mean speed of the jet, U and the characteristic dimension was chosen as the trailing edge thickness.t. The better characteristic dimension would have been the boundary layer thickness, but fortunately the two dimensions were almost the same. The measured sound field was clearly dipole-like (modified slightly by the plate presence).

The lower figure on the right shows a number of turbulent sound spectra measured at various speeds.[44] The frequencies were Strouhal number scaled with U and the sound levels were scaled with the dipole sound power rule of ${\displaystyle {{U}^{6}}}$ over a speed range of 3 to 1. The data fit was quite good, confirming dynamic similarity and the dipole model. The slight discrepancy in level and frequency overlap suggests that both the dimensionless force and the Strouhal number had weak dependence on the Reynolds number.

Another characteristic dimension is the airfoil chord. In these tests the jet width was sufficient to keep the vortex shedding coherent across it. On an airfoil there would be a correlation length less than the wingspan resulting in a number independent dipoles arrayed laterally. The sound power would be diminished somewhat. Since the dipole model is based on the time rate of change of the force, reduction of sound power might be accomplished by reducing that rate. One possible means would be for the opposite sides of the surface to gradually sense each other spatially prior to the trailing edge and thus reduce the rate at the edge. This might be done by a section of graduated porous or flexible materials.

### Circular saw whistle

An edge tone occurs when a jet impinges on a fixed surface. A trailing edge tone occurs when an exterior flow passes over a trailing edge. There is a whistle that is a combination of an edge tone and a trailing edge tone and might be called a wake-edge tone. It occurs in rotating circular saws under idling conditions and may be called the circular saw whistle. Under load conditions, blade vibration plays a role which is not addressed here.

There have been several studies of the fundamental sound generating mechanisms of this whistle.[45][46][47][48][49]

A drawing of typical blade construction is shown in the figure on the right. Research has shown that the sound field is dipole with the primary axis perpendicular to the blade plane. The sources are fluctuating forces acting on each cutting blade. Bies determined that the characteristic speed was the blade velocity, ${\displaystyle R\Omega }$, and the characteristic dimension was the tooth area. Other researchers used blade thickness as the characteristic dimension. Cho and Mote found that the Strouhal number ${\displaystyle St={\frac {fh}{U}}}$ was between 0.1 and 0.2 where h was the blade thickness. Poblete, et al., found Strouhal numbers between 0.12 and 0.18. If the edge tone is relevant, perhaps the characteristic dimension should be the gap between the blades.

The researchers deduced that the fluctuating force was proportional to ${\displaystyle {{U}^{2}}}$, but the sound power was found to vary from ${\displaystyle {{U}^{4.5}}}$ to ${\displaystyle {{U}^{6.0}}}$. If the measurement bandwidth is broad and the measurement distance is out of the near field, there are two dynamic factors (Strouhal number and dimensionless force), that can cause the exponent to be less than six. Both the Deltameter and hole tone data show the Strouhal number is a weak negative function of Reynolds number, which is squared in the sound power equation. This would result in a reduced speed exponent. This factor is not likely to explain the large reduction in exponent however.

The blade geometry was highly variable in the tests, so it is likely that the negative dependence of the dimensionless force on Reynolds number is the major factor. This whistle has two features that separate it from the other whistles described here. The first is that there are a multiplicity of these dipole sources arrayed around the periphery, that most likely are radiating at the same frequency, but incoherently. The second is that blade motion creates a steady, but rotating, pressure field at each blade. The rotating steady force creates a rotating dipole field which has an influence in the geometric near field. The feedback is Class I (hydrodynamic) and there is no indication that stages other than Stage 1 occur.

### Ring tone

Ring flow field
Ring sound

The word "ring" here refers to something akin to that worn on the finger and not the sound from a cellphone. The flow from a circular orifice impinging on a toroidal ring of the same diameter as the orifice will result in a tone; it is called a ring tone. It is similar to the hole tone described above except that because the plate was replaced by a ring a fundamental change in the resultant sound field occurs. Small disturbances at the ring feed back to the orifice to be amplified by the flow instability (Class I). The unstable flow creates a set of symmetric (ring) vortices that later impinge on the physical ring.

The passage of a vortex by the ring is shown schematically in the figure on the right in three steps. The flow vectors in the figure are merely suggestive of direction.

When two vortices are equidistant from the ring, one being beyond and the other approaching, the net circulation around the ring is zero; the null point for the flow oscillation. Each vortex creates a circular (ring) flow field whose axis varies slightly from the vertical as it passes. The figure suggests that the main component of the force on the physical ring is in the direction of the jet flow. If the vortex is a true ring (all parts are in phase), a dipole sound field directed along the jet axis is created.

The figure also suggests that there is a lateral component of force which can only be interpreted as a weak radial dipole. Experiments have been performed on the ring tone.[50] The lower figure on the right shows the relationship of frequency to Reynolds number. If the Strouhal number were graphed in lieu of the frequency, it would have shown that contours were reasonably constant similar to those for the hole tone. Close examination of the data in the figure showed a slight negative dependence of Strouhal number on Reynolds number.

It appears that this whistle has only two stages. The sound field was measured and clearly indicated a dipole whose axis was aligned with the jet axis. Since there were no reflecting surfaces near the source, the data also indicated that a weaker radial dipole component also existed. Such a field can only exist if there is a time delay at a distant point between each of the force components.

### Inaudible whistles

Most of the whistles described generate nearly pure tones that can be heard. The mountain tones discussed above are examples of tones that are inaudible because they are below the frequency range of humans. There are others whose sound levels are below the audible range of humans. For example, the vortex street behind a stick underwater might radiate at audible frequencies but not sufficiently to be heard by a Scuba diver. There are others that are both below audible frequencies and below audible levels.

An unstable water jet, similar to the one shown in the flow instability section above, was not disturbed deliberately, but was allowed to rise to a free water surface. On contact with the surface, a slight jet asymmetry caused an unsymmetrical raised surface that fed back to the jet origin and began a process that looked like a weak version of the flow instability figure. If the jet was not powered, but warmer than the surrounding fluid, it would rise and when encountering the surface would generate a similar feedback system.

Such a phenomenon was observed, but not photographed, in the Owens Valley of California. Early in the morning with no wind, thin clouds were observed to form above the valley. The distinction was that they were created alternately and moved in opposite directions away from a central location on the valley floor, suggesting the existence of an inaudible free convection whistle. The reason for including this type of whistle is that we tend to think that it is necessary for a forced jet flow to encounter a solid material to create a whistle. Perhaps it would be more correct to generalize the concept to a particular impedance mismatch rather than a solid object. The Hartmann whistle and the jet screech fits into this generalization. The concept also applies to any fluid motion as opposed to a strictly forced flow.

### Vortex whistle

The vortex whistle
Vortex Whistle frequencies

When the swirling flow within a pipe encounters the exit, it can become unstable. An example of the original system is shown in the figure on the left. The instability arises when there is a reversed flow on the axis.

The axis of rotation itself precesses around the pipe axis resulting in a rotating force at the pipe exit and results in a rotating dipole sound field. Studies of this whistle[51][52] have shown that dynamic similarity based on the pipe diameter d as the characteristic length scale, and inlet mean flow speed U as the characteristic speed was not achieved, as shown in the lower figure on the right. A more correct speed would be that characteristic of the swirl fd, where f is the precession (and sound) frequency, based on the Rossby number. To test the relevance of this new characteristic speed, the flow rate was increased and the frequency and level of the sound was measured. Using the dipole model, the calculated force was found to be nearly proportional to ${\displaystyle {{(fd)}^{2}}}$, confirming the correctness of the new characteristic speed.

Measurements showed that the vortex whistle was created by a rotating asymmetric vortex which created a rotating force vector in the plane of the exit and a rotating dipole sound field. The phenomenon of swirl instability has been shown to occur in other situations.[53] One was the flow separation on the upper side of delta-shaped airfoils of high speed aircraft (Concorde), The angle of attack of the leading edge resulted in a swirl flow that became unstable. Another is the flow within cyclone separators; the swirling flow there occurs in an annular region between two tubes. The flow reverses at the closed end of the outer tube and exits through the inner tube. Under certain conditions, the flow in the reversal region becomes unstable, resulting in a period rotating force on the outer tube.

Periodic vibration of a cyclone separator would indicate vortex instability. Large centrifugal fans sometimes use radial inlet blades that can be rotated to control the flow into the fan; they create a swirling flow. At near shutoff, where the swirl is very high, rotating blade stall of the fan blades occurs. Although not researched, it is highly likely that swirl instability is the cause. The feedback is clearly hydrodynamic (Class I) and there is no indication that more than one stage occurs.

### Swirl meter

Flow meter

The method of creating swirl in the vortex whistle was considered the cause for lack of dynamic similarity, so the swirl was created in a pipe with a contraction having swirl blades followed by an expansion to create the required axial backflow. This was the vortex whistle in a pipe. Measurements made with this geometry, are shown in the figure on the right. As can be seen, dynamic similarity was achieved with both air and water. This whistle became a flow meter called the swirlmeter. Its accuracy rivals that of the vortex shedding meters described above but has a higher pressure drop. The feedback is hydrodynamic (Class I) and only one stage was found.

### Edge tone

edge tone vortices
edge tone dipole-like flow
Wake edge tone

When a rectangular jet impinges on a sharp edged object such as a wedge, a feedback loop can be established resulting in a nearly pure tone. The figure on the right shows schematically the circulation of two vortices as they pass the wedge. This simple diagram suggests that there is a force applied to the wedge whose angle varies as the vortices pass.

As found in the Aeolian tone, the vertical component (lift) is large and results in a dipole-like sound field at the wedge (shown in the lower figure) and a much weaker horizontal component (drag) at twice the frequency (not shown). The drag component may contribute as part of the driving force for musical instruments (discussed below). A seminal study by Powell[54] of this phenomenon has exposed many details of the edge tone phenomenon. He showed that this whistle has three stages and the feedback loop was hydrodynamic (Class I). A semi-empirical equation for the frequency, developed by Curle,[55] when converted to Strouhal Number, is

${\displaystyle S{{t}_{n}}={\frac {fh}{U}}=\left[{\frac {4n+1}{8}}-{\frac {h}{60d}}\right]}$

This equation, applicable for ${\displaystyle {}^{h}\!\!\diagup \!\!{}_{d}\;>10}$, shows the mean speed U of the jet at the orifice as the characteristic speed and the distance h from orifice to the edge as the characteristic dimension. The integer n represents the various vortex modes. It also suggests that dynamic similarity is achieved to a first approximation; one deviation is that the speed at the wedge, which is less than that at the orifice, should be the characteristic speed. A weak negative Reynolds number effect is likely. The orifice width d also has some influence; it is related to vortex size and lateral correlation of the shedding process.

The presence of a dipole sound field and a periodic force proportional to ${\displaystyle {{U}^{2}}}$ was confirmed by Powell. Numerical simulations of the edge tone and extensive references can be found in a NASA report.[56] The lower figure on the right may be called a wake edge tone. If the preferred frequencies of the trailing edge instability match the preferred frequencies of the free edge tone, a stronger dipole sound should arise. There does not appear to be any research on this configuration.

### Shallow cavity tone

The study of sound generated by flow over cavities at high speed has been well funded by the federal government so a considerable amount of effort has been made. The problem relates to flow over aircraft cavities in flight such as bomb bays or wheel wells. Flow over a cavity in a surface can result in excitation of a feedback loop and nearly pure tones. Unlike the edge tone noted above, the cavity edge is typically square, but also can be an edge as part of a thin structural shell. Cavities can be separated into shallow or deep ones, the difference being that for deep cavities a Class III (acoustical) feedback path may be controlling. Shallow cavities are addressed here and are those in which the cavity length L is greater than the cavity depth D.

At high speeds U, the flow is turbulent and in some studies the speed can be supersonic and the sound generated can be quite high level. One study[57] has shown that several modes of oscillation (stages) can occur in a shallow cavity; the modes being related to the number of vortices in the distance L. For shorter cavities and lower Mach Numbers, there is a shear-layer mode, while for longer cavities and higher Mach Numbers there is a wake mode. The shear-layer mode is characterized well by the feedback process described by Rossiter. The wake mode is characterized instead by a large-scale vortex shedding with a Strouhal number independent of Mach number. There is an empirical equation for these data; it is called Rossiter’s formula.

Lee and others[58][59] have shown it in Strouhal number form as

${\displaystyle S{{t}_{n}}={\frac {{{f}_{n}}L}{U}}={\frac {n-\beta }{U\left({\frac {1}{{c}_{0}}}+{\frac {1}{{u}_{v}}}\right)}}}$

The bracketed term includes two feedback loop speeds; the downstream speed is the speed of the vortices u and the upstream speed is that of sound${\displaystyle {{c}_{0}}}$. The various modes are described by an integer n with an empirical delay constant β (near 0.25). The integer n is closely related to the number of vortices en route to the edge. It is clear from shadowgraphs that the fluctuating force near the downstream edge is the sound source. Since the Mach number of the flow can be appreciable, refraction makes it difficult to determine the major axis of the dipole-like sound field. The preferred frequencies in shallow cavities are different from those for the edge tone.

### Police whistle

police whistle operation

The title above is used since it is commonly used to describe whistles similar to those used by police in America and elsewhere. There are a number of whistles that operate in the same manner as the police whistle and there are number of whistles that are used by police elsewhere that do not operate in the same manner as the police whistle.[60] The London Metropolitan police use a linear whistle, more like a small recorder. Police whistles are commonly used by referees and umpires in sporting events.

The cross section of a common whistle is shown in the figure on the right. The cavity is a closed end cylinder (3/4 inch diameter), but with the cylinder axis lateral to the jet axis. The orifice is 1/16 inch wide and the sharp edge is 1/4 inch from the jet orifice. When blown weakly, the sound is mostly broad band with a weak tone. When blown more forcefully, a strong tone is established near 2800 Hz and adjacent bands are at least 20 dB down. If the whistle is blown yet more forcefully, the level of the tone increases and the frequency increases only slightly suggesting Class I hydrodynamic feedback and operation only in Stage I.

There does not appear to be any detailed research on police whistle operation. Considering the edge tone, noted above, one might expect several jumps in frequency, but none occur. This suggests that if multiple vortices exist in the unstable jet, they do not control.

The diagram on the right suggests a plausible explanation of whistle operation. Within the cavity is an off-center vortex. In the upper drawing, the vortex center is near the jet; the nearby cavity flow is slower and the pressure is less than atmospheric so the jet is directed into the cavity. When the jet moves toward the cavity an additional thrust is given to the interior vertical flow which then rotates around and back to the edge. At that point, the cavity flow and the local pressure are sufficient to force the jet to move away from the cavity.

An interior vortex of this type would explain why no frequency jumps occur. Since the excess fluid in the cavity must be discharged, the jet lateral movement must be considerably larger than that found in the edge tone; this is likely the reason for the high level sound. The flow over the edge results in an applied force and a dipole-like sound field. The characteristic speed must be U the jet exit speed. The characteristic dimension must be D, the cavity diameter.

The frequency of the sound is closely related to the rotation rate of the cavity vortex. With a frequency near 2800 Hz the interior rotation rate must be very high. It is likely that the Rossby number ${\displaystyle {\frac {U}{fD}}}$ would be a valuable dynamic similarity number. The Boatswain's pipe is similar to the police whistle except the cavity is spherical creating a more complex vortex.

### Pea Whistle/Referee's Whistle

A pea whistle is constructionally identical to a 'police whistle' but the chamber contains a small ball, known as the pea but usually a material such as plastic or hard rubber. When blown the pea moves chaotically in the chamber, interrupting and modulating the airflow to create a typical warbling/shrieking effect. Such whistles are traditionally used by association football referees and those of other games.

### Samba Whistles

Similar to pea whistles, samba whistles have a small ball or dowel to create the same sort of sound, but often also have two extensions either side of the chamber. None, one or both of these can be blocked to create a 'tri-tone' effect. The Apito de Samba is a traditional Portuguese example of a samba whistle.

### Levavasseur whistle

Levavasseur whistle

This whistle is essentially the police whistle turned into a torus, magnifying its sound making potential. A cross-section through the middle of the whistle is shown in the figure on the right.

An annular duct carries the fluid that creates the annular jet. The jet impinges on a sharp ended ring with two toroidal cavities on either side. In Levavasseur's patent,[61] a structure is added downstream of the annular opening to act as a coupling horn to direct the sound. The sound generated is very intense. It appears that no scientific study has been done to elucidate the detailed feedback mechanisms of its operation, although it is clear that this whistle has Class I feedback mechanism, similar to the police whistle.

The characteristic speed U is that of the annular jet. The characteristic dimension D is the cavity diameter and it appears that both cavities have similar dimensions. Again, the Rossby number V${\displaystyle {\frac {U}{fD}}}$ is likely to be a relevant dynamic number, since the operation of the inner cavity must be similar to that in the police whistle. It is likely that the vortex in the outer cavity is in anti-phase with the inner cavity to amplify jet displacement and thus the sound output.

### Screech tone

Strong tones can occur in both rectangular and circular jets when the pressure ratio is greater than the critical and the flow becomes supersonic on exit, resulting in a sequence of repetitive shock cells. These cells can be seen in the exhaust of rockets or jets operating with an afterburner. As with subsonic jets, these flows can be unstable.

In a rectangular jet, the instability can show as asymmetric cell distortions. The asymmetry sends waves back to the nozzle which sets up a Class III feedback loop and a strong periodic dipole sound field; it is called screech tone. Powell[62][63] first described the phenomenon and because of application to military aircraft and potential structural fatigue, much subsequent work has been done. The sound field is sufficiently intense for it to appear on a shadowgraph as shown in the figure on the right (from M.G. Davies) for a rectangular supersonic jet. The dipole nature of the source is clear by the phase reversal on either side of the jet. There is lateral motion of the shock cells that gives the dipole its axis.

Supersonic flows can be quite complex and some tentative explanations are available.[58][64] As with hole and ring tones, these jets can be sensitive to local sound reflecting surfaces.

The characteristic speed, U, is that in the exit plane, and the characteristic dimension L is the nozzle width, to which the cell dimensions are proportional. Circular supersonic jets also generate screech tones. In this case, however, there can be three modes of motion: symmetric (toroidal), asymmetric (sinuous), and helical.[65][66][67] These whistles are unlike the others listed above; the sound is generated without interaction with a solid; it is truly an aerodynamic whistle.

### Fluidic oscillators

Fluidic device
Strouhal number
Strouhal number

These devices are whistles that do not radiate sound, but are still aerodynamic whistles. The upper figure on the right shows the basic arrangement of one version of the device. The circle on the left is the fluid source (air or liquid). A jet is formed that either goes into the upper or lower channel.

The black lines are the feedback paths. If the fluid is in the lower channel, some fluid is fed back to the jet origin via the black tube and pushes the jet to the upper channel.

There has been considerable development of these devices from circuit switches that are immune to electromagnetic pulses to more modern uses.

One uniqueness of this whistle compared to the others described is that the length of the feedback path can be chosen arbitrarily. Although the channels are divided by a wedge shape, edge tone operation is avoided by the Coanda effect.[68] The second figure on the right shows results from one study[69] indicating a constant Strouhal number with Reynolds number. The data had been normalized to a reference value.

In another study[70] one set of their frequency data was recalculated in terms of Strouhal number and it was found to rise slowly and then be constant over a range of flow rates. Kim[71] found a similar result; the Strouhal number increased with Reynolds number and then stayed constant as shown in the lower figure on the right. Another uniqueness of this whistle is that the feedback is sufficiently strong that the jet is bodily diverted instead of depending on flow instability vortex development to control it. The geometry of the device suggests that it is essentially a dipole source that operates in Stage I with Class I (hydrodynamic) feedback.

## Monopole-dipole whistles

There are a number of whistles that possess the characteristics of both monopole and dipole sound sources. In several of the whistles described below, the driving source is dipole (generally an edge tone) and the responding source is a monopole (generally a tube or cavity in proximity to the dipole).

The fundamental difference of these whistles from those described above is that there are now two sets of characteristic variables. For the driving source, the characteristic speed is U, and the characteristic dimension is L1. For the responding source, the characteristic speed is ${\displaystyle {{c}_{0}}}$ and the characteristic dimension is L2, typically the corrected cavity depth or tube length. The non dimensional descriptors for each of these are the fluid mechanical Strouhal number and the acoustical Strouhal number. The tie between these two numbers is the commonality of the frequency.

### Jug whistle

Blowing over the edge of a jug or bottle can create a nearly pure tone of low frequency. The driving force is the flow over the jug edge so one might expect an edge tone dipole sound field. In this case, The curvature and roundness of the edge makes a strong edge tone unlikely. Any periodicity at the edge is likely submerged in the Class III feedback from the jug volume. The unsteady edge flow sets up a classical Helmholtz resonator response in which the interior geometry and the jug neck determines the resultant frequency. A resonance equation is shown below.[72]

{\displaystyle {\begin{aligned}&\cot \left(kL\right)={\frac {{{A}_{c}}\left({{L}_{o}}+{{\delta }_{e}}+{{\delta }_{i}}\right)}{{A}_{0}}}kL\\&kL=2\pi {\frac {fL}{{c}_{0}}}=2\pi {{S}_{t}}\\\end{aligned}}}

It is a transcendental equation where Ac is the cross sectional area of a cylindrical cavity of depth L. Ao is the area of the circular orifice of depth Lo, δe is the exterior end correction, δi is the interior end correction, and kL is the Helmholtz number (acoustical Strouhal Number with ${\displaystyle 2\pi }$ added). A cylindrical cavity nine inches deep and 4.25 inches in diameter was connected to a circular orifice 1.375 inches in diameter and 1.375 inches deep.[73] The measured frequency was close to 140 Hz. If the cavity acted as a one-quarter wavelength resonator, the frequency would have been 377 Hz; clearly not a longitudinal resonance.

The equation above indicated 146 Hz and the Nielsen equation[74] indicated 138 Hz. Clearly, the whistle was being driven by a cavity resonance. This is an example of a whistle being driven in edge tone fashion but the result is a monopole sound field.

### Deep cavity tone

Flow over a cavity that is considered deep can create a whistle similar to that over shallow cavities. Deep is generally distinguished from shallow by the cavity depth being greater than the width. There are two geometries that have been studied. The first geometry is flow exterior to the cavity such as on an aircraft.[75][76][77][78]

There are two characteristic dimensions (cavity width L associated with vortex development and cavity depth D associated with acoustical response). There are two characteristic speeds (flow speed U associated with vortex development and sound speed ${\displaystyle {{c}_{0}}}$ associated with cavity response). It was found that the feedback was Class III and the Strouhal Numbers ranging from 0.3 to 0.4 were associated with a single vortex pattern (Stage I) across the gap while Strouhal numbers ranging from 0.6 to 0.9 were associated with two vortices across the gap (Stage II).

The second geometry is flow in a duct with a side branch. Selamet and his colleagues[79][80][81] have made extensive studies of whistle phenomena in ducts with side branches that are closed at one end. For these studies, The cavity depth was L and D was the side branch diameter. The fluid mechanical Strouhal and acoustical Strouhal numbers were

{\displaystyle {\begin{aligned}&St={\frac {fD}{U}}\\&S{{t}_{a}}={\frac {f{{L}_{1}}}{{c}_{0}}}={\frac {\left(2n+1\right)}{4}}\\&{{L}_{1}}=L\left(1+\beta {\frac {D}{L}}\right)\\\end{aligned}}}

An arbitrary constant β was used to represent the impedance at the junction of the side branch with the duct. n was an integer representing the stage number. They noted that the Strouhal number remained constant with increase of speed.

### Pipe organ

The pipe organ is another example of a potentially dipole sound source being driven as a monopole source. An air jet is directed at a sharp edge setting up flow oscillations as in the edge tone. The edge is part of a generally cylindrical tube of length L. An example is shown in the figure on the right. The unstable jet drives fluid alternately into the tube and out. The streamlines clearly are distorted from those of the free edge tone. There is a stagnation point opposite the source. The dashed lines, colored in red, are those most strongly modified. The red streamlines in the tube are now augmented by the oscillatory flow in the tube, a superposition of resistive and reactive dipole flow and resistive acoustic flow.

The tube length determines whether the tube acoustic pressure or velocity is the dominant influence on the frequency of the tube. Simple models of organ pipe resonance is based on open-open pipe resonance (${\displaystyle \lambda ={}^{L}\!\!\diagup \!\!{}_{2}\;}$ but corrections must be made to take into account that one end of the pipe radiates into the surrounding medium and the other radiates through a slit with a jet flow. Boelkes and Hoffmann[82] have made measurements of end correction for open-open tubes and derived the relation δ=0.33D. This cannot be exact since the driving end is not open.

The radiation ±impedance at the driving end should move the tube toward a${\displaystyle {}^{\lambda }\!\!\diagup \!\!{}_{4}\;}$ condition, further lowering the frequency. Since there are two coupled systems, so there are two characteristics scales. For the pipe component, the characteristic dimension is L and the characteristic speed is ${\displaystyle {{c}_{0}}}$. For the edge tone component, the characteristic dimension is the orifice to edge distance h and the characteristic speed is that of the jet U. It would seem that the maximum oscillatory gain of the system would occur when the preferred pipe frequency matches the preferred edge tone frequency with suitable phase. This relationship expressed in terms of Strouhal Numbers is:

{\displaystyle {\begin{aligned}&S{}_{tc}={\frac {f{{L}_{1}}}{{c}_{0}}},{{S}_{te}}={\frac {fh}{U}}\\&{{S}_{tc}}={{S}_{te}}{\frac {{L}_{1}}{h}}M\\\end{aligned}}}

If dynamic similarity holds for both resonances, the latter equation suggests how organ pipes are scaled. The apparent simplicity of the equation hides important variable factors such as the effective pipe length ${\displaystyle {{L}_{1}}=L+{{\delta }_{1}}+{{\delta }_{2}}}$ where δ1 is correction for the open end and δ2 is the correction for the end near the jet. The jet disturbance (vortex) speed from orifice to edge will vary with mean speed U, edge distance h, and slit width d as suggested in the Edge Tone section.

The Strouhal relationship suggests that the jet Mach number and the ratio of effective pipe length to the edge distance are important in a first approximation. Normal pipe operation would be a monopole sound source in Stage I with Class III feedback.

### Flutes, recorders and piccolos

Native American Style flute

A number of musical instruments, other than the pipe organ, are based on the edge tone phenomenon, the most common of which are the flute, the piccolo (a small version of the flute), and the recorder. The flute can be blown lateral to the instrument or at the end as the other ones are. A native end blown flute is shown in the figure.

They are all subject to frequency jumps when overblown, suggesting the dipole-monopole relationship. The monopole aspects are relatively fixed. The characteristic dimension of the tube, L2, is fixed; the characteristic speed,${\displaystyle {{c}_{0}}}$, is fixed. The effective length of the tube is fixed since the radiation impedances at each end are fixed. Unlike the pipe organ, however, these instruments have side ports to change the resonance frequency and thus the acoustical Strouhal number.

The dipole aspects are also relatively fixed. The jet orifice dimension and the distance to the edge, h is fixed. Although the jet speed U can vary, the fluid mechanical Strouhal number is relatively constant and normally operates in Stage I. When there is phase coherent gain of the two aspects, they operate as Class III monopole sources. The efficiency of the monopole radiation is considerably greater than that of the dipole so the dipole pattern is noticed, The details of system gain and interaction between these two dynamic systems is yet to be fully uncovered. It is a testimony to the skills of early instrument makers that they were able to achieve the right port sizes and positions for a given note without scientific measurement instruments.

## Singing Sand

For many years, people have heard tones from the motion of sand, mostly on sand dunes. In many cases, the tune was akin to a whistle in that it could be at a single frequency. In most cases, the wind over the peak of a dune sets the sand on the dune into motion, and it is that motion that generates the sound.

In the recent past, a number of researchers have attempted to explain the origin of the sound.[83][84] One might consider this source an aerodynamic whistle since it is wind caused and one proposal is that the sand moves as an unstable fluidized sand/air mixture.

There are several other theories, but the aerodynamic whistle possibility is not correct since it can occur by simple mechanical motion of the sand. It is likely a stick/slip mechanism which would be a dipole like source rather than a volumetric (monopole) oscillation of the sand particles.

## References

1. ^ Arroyos, Rafael Pérez (2003). Egypt: Music in the Age of the Pyramids (1st ed.). Madrid: Centro de Estudios Egipcios. p. 28. ISBN 8493279617.
2. ^
3. ^ Tacoma Narrows Bridge (1940)
4. ^ Chanaud, Robert C. (1970). "Observations of Oscillatory Radial Flow between a Fixed Disk and a Free Disk". The Journal of the Acoustical Society of America. 47 (5B): 1471–2. doi:10.1121/1.1912065.
5. ^
6. ^ History of the Whistle
7. ^ Cross, David (2011-02-17). "On the Beat in Birmingham - Rules and regulations". BBC. Retrieved 11 March 2014. Police whistles came much later; the early Victorian constable would have carried a small wooden rattle.
8. ^ Taylor, J. "The Victorian Police Rattle Mystery" The Constabulary (2003) Archived February 18, 2010, at the Wayback Machine.
9. ^ ACME Whistles
10. ^ Wilson, T.A Beavers G. S. DeCoster M. A., Holger D. K., Regenfusd, M. D."Experiments on the fluid mechanics of whistling", J. Acoust. Soc. Am. 1, pp366-372, 1971
11. http://www.scientificamerican.com/article/aerodynamic-whistles/
12. ^ Powell, A. “On the Edge Tone”, Jo. Acoust. Soc. Am., 33, No. 4, 395-409 1961
13. ^ Chanaud,R.C.,MS Thesis University of California, Los Angeles, 1960
15. ^ Ernst Mach
16. ^ Osborne Reynolds
17. ^ Carl-Gustaf Rossby
18. ^ Chanaud, R.C., Powell, A. “Some Experiments Concerning the Hole and Ring Tone”, Jo. Acoust. Soc. Am. 37, No.5, pp. 902-911, 1965
19. ^ Henrywood,R.H., Agarwal1,A. "The aeroacoustics of a steam kettle",Phys. Fluids 25, 107101,2013
20. ^ Strutt, J.W. Baron Rayleigh, The Theory of Sound, MacMillan and Co. 1877
21. ^ http://www.tenterfieldfoxwhistle.net/factsheet.asp
22. ^ Nakiboğlu G, Rudenko O, Hirschberg A. "Aeroacoustics of the swinging corrugated tube: voice of the Dragon" J Acoust Soc Am. 31, 749-765,2012
23. ^ Rajavel, B., Prasad, M.G. "Parametric studies on acoustics of corrugated tubes using large eddy simulation (LES)", Noise Control Engineering Journal, 62(4),2014
24. ^ Lisa R., Taylor, M.E., "Experimental Study of the Acoustical Characteristics of Corrugated Tubing", Noise and Vibration Control Laboratory, Stevens Institute of Technology, Thesis 1994
25. ^ Karthik, B., Chakravarthy, S.R., Sujith, R.I, “Mechanism of pipe-tone excitation by flow through an orifice in a duct”, Int. Jo. Aeroacoustics, 7, Nos 3-4, 321-347, 2008
26. ^ Anderson, A.B.C. “Dependence of Pfeifenton Frequency on Pipe Length, Orifice Diameter and Gas Discharge Pressure”, Jo. Acoust. Soc. Am. 24, 675-681 1952
27. ^ Anderson, A.B.C., "Circular Orifice Number describing dependency of primary Pfeifenton frequency on differential pressure, gas density and orifice geometry" J. Acoust. Soc. Am., 25, 626-631, 1953
28. ^ Anderson, A.B.C., "Structure and velocity of the periodic Vortex Ring Flow Pattern of a Primary Pfeifenton Jet" J. Acoust. Soc. Am., 27, 1048-1053, 1955
29. ^ Steam whistle
30. ^ Hartmann, J.,Troll, B., "On a new method for the generation of sound waves", Phys. Rev., 20, 719-727, 1922
31. ^ Raman, G., Srinivasan, K., “The Power Resonance Tube: From Hartmann’s Discovery to Current Active Flow Control Applications”, Prog. Aerospace Sci., 45, 97-123, 2009
32. ^ Brun, E., Boucher, R.M.G., "Research in the acoustic air-jet generator: A new development", J. Acoust. Soc. Am., 29, 573-583, 1957
33. ^ Savoy, L.E., "Experiments with the Hartmann acoustic Generator", Engineering, 170, 99-100, 136-138, 1950
34. ^ Rijke tube
35. ^ Matveev, K., "Thermoacoustic instabilities in the Rijke Tube: Experiments and modeling." Thesis, Calif. Inst. Of Tech., 2003
36. ^ [Backhaus,S., Swift, G, "New Varieties of Thermoacoustic Engines", 9th International Congress on Sound and Vibration, 2002
37. ^ Thermoacoustics
38. ^ Wilson, T.A Beavers G.S., DeCoster M.A., Holger D.K., Regenfuss, M.D. "Experiments on the fluid mechanics of whistling", J. Acoust. Soc. Am. 1, 366-372, 1971
39. ^ von Karman, T. "Aerodynamics", McGraw-Hill, 1963
40. ^ Etkin, B., Ribner, H. "Canadian research on aerodynamic noise", Review 13, Institute of Physics, Univ. Toronto, 1958
41. ^ Phillips, O.M. “The Intensity of Aeolian Tones”, Jo. Fluid Mech. 1, 607-624 1956
42. ^ Rodely,A.E., White, D.F., Chanaud,R.C. “A Digital Flowmeter Without Moving Parts”, Amer. Soc. Mech. Engr., 65-WA/FM-6, 1965
43. ^ Hayden, R.E., Fox, H.L., Chanaud, R.C. “Some Factors Influencing Radiation from Flow Interaction with Edges of Finite Surfaces”, NASA CR-145073, 1976
44. ^ Rschevkin, S.N., “The Theory of Sound”, The MacMillan Company 1963
45. ^ Bies, D.A., "Circular saw aerodynamic noise", Jo. Sound & Vibration, 154(3), 495-513, 1992
46. ^ Martin, B.T., Bies, D.A., On aerodynamic noise from vortex shedding in rotating blades", Jo. Sound & Vibration, 155(2), 317-324, 1992<r
47. ^ Mote, C.D.,Jr., Wen Hua Zhu, "Aerodynamic far field noise in idling circular saws", Jo. Vibration, Stress, and Reliability Design', 106, 441-446, 1984
48. ^ Reiter, W.F., Jr., Keltie, R.F., "On the nature of idling noise of circular saw blades", Jo. Sound & Vibration, 44(4),531-543, 1976
49. ^ Poblete, V., Arenas, J.P., Rios, R., Millar, E. "Vibration and idling noise in commercial circular saws", Fifth Inter. Congress on Sound and vibration, 1997
50. ^ Chanaud, R.C., Powell, A. “Some Experiments Concerning the Hole and Ring Tone”, Jo. Acoust. Soc. Am. 37, (5), 902-911, 1965
51. ^ Vonnegut, B. “The Vortex Whistle”, Jo. Acoust. Soc. Am. 26, pp16-20, 1954
52. ^ Chanaud, R.C., “Experiments Concerning the Vortex Whistle”, Jo. Acoust. Soc. Am., 35, No. 7, 953-960, 1963
53. ^ Chanaud, R.C., “Observations of Oscillatory Motion in Certain Swirling Flows”, Jo. Fluid Mech. 21, 111-127, 1965
54. ^ Powell, A. “On the Edge Tone”, Jo. Acoust. Soc. Am., 33, No. 4, 395-409, 1961
55. ^ Curle, N. “The Mechanics of Edge Tones”, Proc. Roy. Soc. A231, 505, 1955
56. ^ Dougherty, B.L., O'Farrell, J.M., Numerical Simulation of the Edge Tone Phenomenon",NASA Contractor Report 4581,1994
57. ^ Rossiter, J.E. “Wind Tunnel Experiments on the Flow over Rectangular Cavities at Subsonic and Transonic speeds.”, Report 3438, Aeronautical Research Council (UK), 1964
58. ^ a b Lee, D.J, Lee,, I.C., Heo, D.N., Kim, Y.N., “Numerical Analysis of Aerodynamic Noise from Feedback Phenomena using Computational Aeroacoustics (CAA)”, Proc. 12th Asian Congress of Fluid Mechanics, August 2008
59. ^ Rowley, C.W., Colonius, T., Amit, J. "On self-sustained oscillations in two-dimensional compressible flow over rectangular cavities", Journal of Fluid Mechanics, 455 . 315-346. 2001
60. ^ Whistle
62. ^ Powell, A. "On the Mechanism of Choked Jet Noise”, Proc. Phys. Soc., B66, 1039-1056, 1953
63. ^ Powell, A. “The Reduction of Choked Jet Noise”, Proc. Phys. Soc., B67, 313-327, 1954
64. ^ Lin, D., Powell, A. “Symmetrical Oscillation Modes in Choked-jet Edge Tones and Screech from Rectangular Nozzles”, J. Acoust. Soc. Amer., 102(2), 1235-1238, 1997
65. ^ Davies, M.G., Oldfield, D.E.S., "Tones from a choked axisymmetric jet. I. Cell structure, eddy velocity and source locations”, Acta Acustica, 12, 257-276 1962
66. ^ Davies, M.G., Oldfield, D.E.S. "Tones from a choked axisymmetric jet. II. The self excited loop and mode of oscillations", Acta Acustica, 12, 267-277 1962
67. ^ Powell, A. “Observations of the Oscillation Modes of Choked Circular Jets”, J. Acoust. Soc. Amer., 92(5), 2823-2836, 1992
68. ^ Coandă effect
69. ^ Tesar, V., Peszynski, K., "Strangely Behaving Fluidic Oscillator" EPJ Web of Conferences, 45,01074,2013
70. ^ Gregory, J.W., Sullivan, J.P., Raman, G, Raghu,S. "Characterization of the Microfluidic Oscillator" AIAA J. 45(3), 568-576, 2007
71. ^ Kim, G, "A Study of Fluidic Oscillators as an alternative pulsed vortex generating jet actuator for flow separation" Masters Thesis, University of Manchester, 2011
72. ^ Chanaud, R.C. “Effect of Geometry on the Resonance Frequency of Helmholtz Resonators”, Jo. Sound & Vibration, 178, No. 3, 337-348, 1994
73. ^ Chanaud, R.C., Unpublished work
74. ^ Nielsen, A.K., "Acoustical resonators of circular cross section and with axial symmetry", Transactions of the Danish Academy of Technical Science, 1949
75. ^ Durgin, W.W., Graf, H.R. "Excited Acoustic resonance in a Deep Cavity: An Analytical Model", Synposium on Flow-induced Vibrtion and Noise, 7, 1992
76. ^ East, F., "Aerodynamically induced resonance in rectangular cavities", Jo. Sound and Vibration, 3, 277-287, 1966
77. ^ Forestier,N.,Jacquin, L., Geffroy, P., "The mixing layer over a deep cavity at high-subsonic speed", Jo. of Fluid Mechanics, 475, 101-145, 2003
78. ^ Rowley, C.W., Colonius, T., Basu, A.J. "On self-sustained oscillations in two-dimensional compressible flow over rectangular cavities", Jo. of Fluid Mechanics, 455, 315-346. 2001
79. ^ Selamet, A., Kurniawan, D., Knotts, B.D., “Whistles with a Generic Side Branch: Production and Suppression”, Jo. Sound & Vibration, 250(2), 277-298, 2002.
80. ^ Knotts, B.D., Selamet, A., “Suppression of Flow-acoustic Coupling in Side Branch Ducts by Interface Modification”, Jo. Sound & Vibration, 265, 025-1045, 2003
81. ^ Radvich, P.M., Selamet, A., Novak, J.M., “A Computational Approach for Flow-acoustic Coupling in Closed Side Branches”, Jo Acoust. Soc. Amer., 109(4), 1343-1351, 2001
82. ^ http://www.isb/ac./HS/JoP/index.html
83. ^ Andreotti, B., Bonneau, L., Clement, E., "Comment on solving the mystery of booming sand dunes by Nathalie M. Vriend, et. al. ",Geophysical Research Letters, 35, 2008
84. ^ Bohy. S., DuPont, S., Douahy, S., "Singing sand avalanches without dunes", Geophysical Research Letters, 39, 2012