Issue 
Acta Acust.
Volume 4, Number 6, 2020



Article Number  27  
Number of page(s)  12  
Section  Musical Acoustics  
DOI  https://doi.org/10.1051/aacus/2020026  
Published online  11 December 2020 
Scientific Article
Minimal blowing pressure allowing periodic oscillations in a simplified reed musical instrument model: BouasseBenade prescription assessed through numerical continuation
^{1}
Laboratoire d’Acoustique de l’Université du Mans, UMR CNRS 6613, 72085 Le Mans, France
^{2}
Laboratoire Manceau de Mathématiques – Le Mans Université, 72085 Le Mans, France
^{3}
Aix Marseille Univ, CNRS, Centrale Marseille, LMA, UMR 7031, France
^{*} Corresponding author: joel.gilbert@univlemans.fr
Received:
19
March
2020
Accepted:
6
November
2020
A reed instrument model with N acoustical modes can be described as a 2N dimensional autonomous nonlinear dynamical system. Here, a simplified model of a reedlike instrument having two quasiharmonic resonances, represented by a four dimensional dynamical system, is studied using the continuation and bifurcation software AUTO. Bifurcation diagrams of equilibria and periodic solutions are explored with respect to the blowing mouth pressure, with focus on amplitude and frequency evolutions along the different solution branches. Equilibria and periodic regimes are connected through Hopf bifurcations, which are found to be direct or inverse depending on the physical parameters values. Emerging periodic regimes mainly supported by either the first acoustic resonance (first register) or the second acoustic resonance (second register) are successfully identified by the model. An additional periodic branch is also found to emerge from the branch of the second register through a perioddoubling bifurcation. The evolution of the oscillation frequency along each branch of the periodic regimes is also predicted by the continuation method. Stability along each branch is computed as well. Some of the results are interpreted in terms of the ease of playing of the reed instrument. The effect of the inharmonicity between the first two impedance peaks is observed both when the amplitude of the first is greater than the second, as well as the inverse case. In both cases, the blowing pressure that results in periodic oscillations has a lowest value when the two resonances are harmonic, a theoretical illustration of the BouasseBenade prescription.
© J. Gilbert et al., Published by EDP Sciences, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
An important goal of the acoustics of wind instruments is to understand key components of intonation and also the ease of playing. From the physics modelling point of view, it is interesting to study the main variables that control the playing frequency (for intonation) and the minimum mouth pressure to achieve autooscillations (for ease of playing). It is assumed that part of the musician’s judgement of ease of playing of a note is inversely related to the sounding resistance represented by the threshold blowing pressure. Support for this hypothesis is offered by measurements on the performing properties of saxophone reeds by [1]: a significant correlation was found between the softhard scale on which the sounding resistance of different reeds was judged by saxophonists and the threshold pressures measured in the mouths of the performers. A useful overview of the acoustics of reed and lip wind instruments can be found in books such as [2–6].
It is often commented that the flaring bore of brass instruments are designed such that the input impedance are, as close as possible, harmonically related. While this alignment is said to be important for intonation, it is also likely to determine the oscillation threshold and therefore improve the ease of playing. Here, the necessity of an alignment in a harmonic series is called the BouasseBenade prescription because of what Benade wrote in his famous book [2], or in [7]: The usefulness of the harmonically related air column resonances in fostering stable oscillations sustained by a reedvalve was first pointed out by the French physicist Henri Bouasse in his book “Instruments à vents” [8]. In order to illustrate this prescription, a horn was designed to provide an air column whose resonance frequencies (frequencies of maximum input impedance) were chosen to avoid all possible integer relations between them, called “tacet horn” in [7]. The purpose of this instrument is to deliberately make the conditions for oscillation unfavorable.
The effect of inharmonicity of the two first resonance frequencies on both tone colour and ease of playing have been examined experimentally on alto saxophone fingerings during a project for the design of microinterval systems [9]. An decreased harmonicity by extending the bore of the Tintignac carnyx improves its ease of playing [10, 11]. The harmonicity of resonances is also necessary for proper intonation when a reed instrument is expected to play in upper registers, and is often used as a target in optimisation problems [12–16]. This paper focuses on the assessment of the BouasseBenade prescription on a model of reed musical instruments. The goal is to study the influence of the inharmonicity on the playing frequency and on the minimum mouth pressure required to achieve autooscillation in the first register. However, it appears that this mouth pressure cannot always be determined by a study of small amplitude oscillations only. On the contrary, a complete bifurcation diagram, including all periodic branches with the blowing pressure as the continuation parameter, needs to be computed. To achieve this, a simplified model of a reed instrument is derived from a generic model that is valid for both reed and brass instruments, and constitutes a simplified version of the problem. The reed is modeled as a simple spring [17], only two acoustic resonances are taken into account [18, 19], and the nonlinear coupling between the reed and the acoustic resonances through the incoming flow is reduced to a polynomial expansion (Kergomard in [20]). This may be considered as the simplest model of reed instruments that includes inharmonicity. Furthermore, this simple model helps isolate the effects of the main parameters without the added complications that arise when considering real instruments.
Inharmonicity Inh between the acoustic resonances f_{res1} and f_{res2} is defined as the deviation from harmonicity: Inh = f_{res2}/(2f_{res1}) − 1. Therefore once the first resonance frequency f_{res1} and the Inharmonicity Inh are known, the second resonance frequency is fixed through the relation f_{res2} = 2f_{res1} (1 + Inh). Note that if the resonances are exactly harmonic (Inh = 0) the problem can be solved analytically and two bifurcation diagrams have been obtained (Figs. 8 and 10 of [18]).
In Section 2 of this paper the theoretical background, and particularly the equations of the elementary model of reed instruments, are briefly presented. The behaviour of the elementary model at the stability threshold of the equilibrium position, and the nature of the Hopf bifurcations, are discussed in the second part of this section. Section 3 documents the procedure used to calculate bifurcation diagrams using a continuation method, after having reformulated the two equations of the model into a set of four firstorder ODE equations. The influence of the inharmonicity on the bifurcation diagrams is shown and discussed in Section 4. The section is divided into two parts: the first study assumes that the amplitude of the first resonance is larger than the second resonance (Z_{1} > Z_{2}), for which preliminary results have been presented in [21], and the second study the opposite condition (Z_{2} > Z_{1}) is considered. In order to link the BouasseBenade prescription to the ease of playing experienced by musicians, bifurcation diagrams are analysed with respect to the minimal mouth pressure necessary to achieve oscillation. Additionally, the effect of inharmonicity is also considered.
2 Theoretical background
2.1 Elementary acoustical model
The model presented and used in the present publication is labeled as elementary because a number of major simplifications are made in deriving it (see for example Hirschberg in [20, 22]). The vibrating reeds or lips are modeled as a linear onedegreeoffreedom oscillator. The upstream resonances of the player’s windway are neglected, as is nonlinear propagation of sound in the air column of the instrument. Wall vibrations are also ignored. Despite these simplifications, the elementary model is capable of reproducing many of the important aspects of performances by human players on realistic reed and brass instruments (see [4, 5]). The model is based on a set of three equations, which have to be solved simultaneously to predict the nature of the sound radiated by the instrument. These three constituent equations of the model are presented hereafter. Besides the control parameters defining the embouchure of the players, including the reed or lips parameters and the mouth pressure P_{m}, and the input impedance of the wind instrument, there are three variables in the set of the forthcoming three equations as a function of the time t: , the reed or lipopening height, the pressure in the mouthpiece of the instrument, and the volume flow entering the instrument.
In order to describe the vibrating reeds or lips, the first of the three constituent equations of the elementary model is:(1)
In this equation, which describes the reeds or lips as a onedegreeoffreedom (1DOF) mechanical oscillator, the symbols ω_{r}, Q_{r}, h_{o} and μ represent the angular reed resonance frequency, the quality factor of the reed resonance, the value of the reed or lipopening height at rest, and the effective mass per unit area of the reed or lips respectively. These quantities are parameters of the model, which are either constant (in a stable note) or changing slowly in a prescribed way (in a music performance). Note that if μ is positive, an increase of the pressure difference will imply a closing of the reed or lips aperture. It is called the “inwardstriking” model, used mainly for reed instruments. If μ is negative, an increase of the pressure difference will imply an opening of the reed or lips aperture. It is called the “outward striking” model, used preferably for brass (lip reed) instruments.
The second constituent equation describes the relationship between pressure and flow velocity in the reed channel:(2)where the square root originates from the Bernoulli equation, and the positive part of the reed or lips aperture implies that the volume flow vanishes when the reed or lips are closed.
The third and last constituent equation describes the relationship between flow and pressure in the instrument mouthpiece. It is written in its frequency domain form by using the input impedance Z(ω) of the wind instrument:(3)
Other than the difference of sign of μ between inwardstriking reed instruments model and outward striking brass instruments model, there is another difference between these two subfamilies of wind instruments. The control parameter ω_{r} of vibrating lips varies a lot, over four octaves, to get the entire tessitura of a given brass instrument. On the other hand the ω_{r} associated to reeds is more fixed (slightly varying because of the lower lip of the clarinet or saxophone player) and most of the time very large compared to the playing frequencies. This justifies a lowfrequency approximation of the elementary model: ω_{r} is assumed infinite and the reed undamped. In other words, the reed is reduced to its stiffness only and the set of three equations becomes a set of two equations as follows:(4)
When the mouth pressure is too high, the reed can be blocked against the lay of the mouthpiece. Then the closure pressure defined by is the minimal mouth pressure for which the reed remains closed in the static regime ( becomes equal to 0). By using this closure pressure, a dimensionless mouth pressure γ can be defined: γ = P_{m}/P_{M}. It is convenient to define another dimensionless parameter, a dimensionless reed height at rest:(5)where is the characteristic impedance for plane wave inside the resonator of input cross section S, ρ is the air density and c is the sound velocity.
In the following, the nonlinear equation of the model is approximated by its thirdorder Taylor series around the equilibrium position defined by , , and .
The approximated nonlinear equation can be written in the following dimensionless form (see for example Kergomard in [20]):(6)with , , and . The value of the dimensionless reed height at rest ζ is chosen to be equal to 0.1.
It is this elementary lowfrequency model for reed instruments which is studied in the present paper. If a nonbeating reed is assumed which is typically obtained for a dimensionless mouth pressure γ lower than 0.5, the third order approximation of the flow rates is appropriate. The elementary model based on the set of two equations has to be solved to predict the nature of the sound radiated by the instrument. Low amplitude solutions for a few specific cases are reviewed in the following subsection.
2.2 Small amplitude behaviour
The equilibrium position is the trivial permanent (steady) regime corresponding to silence. Sound can happen if the equilibrium position becomes unstable. For a lossless cylindrical air column, it becomes unstable for a specific value of γ which is γ_{thr} = 1/3. If losses are taken into account, then γ_{thr} is a bit higher (see [5]). If the losses are very important, the threshold value γ_{thr} can reach 1 and the reed channel is closed at equilibrium. In this case the equilibrium remains stable for any value of γ. Hence no sound can be produced. For an extensive analysis of stability of the equilibrium position with an experimental comparison for cylindrical air columns, see [23] and [24].
The step beyond the above linear stability analysis is the study of the small oscillations around the threshold. It has been done first by [25] and then extended by analysing the nature of the bifurcation at the threshold by [26] which can be direct or inverse Hopf bifurcation. The results are displayed in Figure 1 as a 2D map where the xaxis is C the third coefficient of the Taylor expansion Equation (6), and the yaxis is 1/Z_{2} − 1/Z_{1}, the difference between the admittance amplitude between the two first resonances (assumed to be harmonic, the ratio between their frequencies, being equal to 2).
Figure 1 Diagram showing the regions where the bifurcation is direct as well as those regions where it is inverse. The xaxis shows the values of the third coefficient of the Taylor expansion Equation (6), and the yaxis is the difference between the admittance amplitude between the two first resonances (assumed to be harmonic) Y_{2} −Y_{1} = 1/Z_{2} − 1/Z_{1}. The hatched region is for a direct bifurcation, and the unhatched region for an inverse bifurcation. Adapted from [26]. 
In our specific case, the coefficient C is negative. Then, for a specific negative C value, following an imaginary vertical line coming from an infinite positive value of 1/Z_{2} − 1/Z_{1} (second resonance peak absent like for the cylindrical tube) the bifurcation is direct. It becomes inverse in a particular point for a particular positive value of 1/Z_{2} − 1/Z_{1} not far from zero: 1/Z_{2} − 1/Z_{1} = −2B^{2}/(3C). And when 1/Z_{2} − 1/Z_{1} becomes negative, and whatever how 1/Z_{2} < 1/Z_{1} (it means whatever Z_{2} > Z_{1}) is, the bifurcation becomes and stays direct. Properties of small amplitude oscillations of the singlereed woodwind instruments near the oscillation threshold have been investigated more recently by using analytical formulae with explicit dependence on the physical parameters of the instrument and the instrumentist allowing to determine the bifurcation point, the nature of the bifurcation, the amplitude of the first harmonics and the oscillation frequency [27]. Apart from a few very simplified cases, such as a clarinetlike model with a lossless cylindrical tube ([28], Kergomard in [20, 5, 29]), or by taking into account losses independent of frequency, sometimes called Raman model [30], the equations are not tractable analytically, and the bifurcation diagrams can not be easilly obtained.
The simplest nontrivial resonator that can be studied, is a resonator having two quasiharmonic resonance frequencies f_{res1} and f_{res2}. This kind of resonator can be obtained in practice in the midle and high ranges of the first register of saxophone (see for example Figs. 17 and 12 of [9]). Bifurcation diagrams have been analytically calculated in [18] in the restrictive case of perfect harmonicity between the two resonances. In the following sections, this kind of resonator but with a nonzero inharmonicty Inh is analysed.
3 Typical bifurcation diagram obtained by continuation method
To overcome the difficulties of the analytical analysis of small amplitude oscillations near thresholds, and to get results for any inharmonicity value arbitrarily far from the oscillation threshold, simulation techniques in time domain are often used. An alternative method is possible. A nice way to have an overview of the dynamics over small and large amplitudes is to use the bifurcation diagram representation. Very few of them can be obtained analytically (see the previous subsection). It is possible to obtain bifurcation diagrams numerically for a large range of situations by using continuation methods, such as in the AUTO software [31] or MANLAB software [32] for example. In order to use AUTO technique in the following section, the elementary model has to be mathematically reformulated in a set of firstorder ODE equations.
The principle of continuation is to seek solution branches of a nonlinear algebraic system rather than solution points. A solution branch is an 1Dcurve in a space whose axes are an unknown to the problem and a parameter of interest called a bifurcation parameter. In the following the dimensionless mouth pressure γ is chosen as the bifurcation parameter. It provides more information than a set of solution points obtained for successive values of the bifurcation parameter. Branches of static and periodic solutions are computed numerically hereafter using the software AUTO, freely available online [33].
The model analysed in this paper is a nonlinear dynamical system. In order to obtain a nonlinear algebraic system in which numerical continuation can be applied, some additional work may be required. For instance, for continuing periodic solutions of a dynamical system, a discretisation is necessary to come down to an algebraic system. Many approaches are possible among which a timedomain discretisation of the (unknown) solution over one (unknown) period. The unknowns of the resulting nonlinear algebraic system are the sampled values of the periodic solution and the period. The time discretisation implemented in AUTO is called orthogonal collocation and relies on the use of Lagrange polynomials. The stability of each solution is also assessed. Stability is a very important information for the interpretation of the bifurcation diagram since only stable solutions are observable. Stability of both equilibria and periodic solutions is found through a linearization of the system of equations around the solution considered. The solution is stable if and only if the real parts of all the eigenvalues of a matrix characteristic of the linearized system are negative. This matrix is the Jacobian matrix if the solution considered is an equilibrium, and the socalled monodromy matrix if the solution considered is periodic. Stability of a solution along a branch is an output of AUTO. For comprehensive details about continuation of static/periodic solutions using AUTO, please refer to [31].
In order to use the AUTO technique, the input impedance equation (Eq. (3)) is reformulated by a sum of individual acoustical resonance modes in the frequency domain, and then translated them in the time domain. There are two ways to manage that: sum of real modes (see for example [34]), sum of complex modes (see for example [35]). These two ways of approximating the input impedance in the frequency domain lead to two different sets of firstorder equations with two different X vectors. In the present paper the real mode representation of the input impedance Z is used.
The modalfitted input impedance with N resonance modes, is written as follows:(7)where the nth resonance is defined by three real constants, the amplitude Z_{n}, the dimensionless quality factor Q_{n} and the angular frequency ω_{n}.
Translation of Equation (7) in the time domain and reconstruction of p(t) from real modal components p_{n}, such that the acoustical pressure is , results in a second order ODE for each p_{n}:(8)
Taking into account the other equation of the elementary model, the time derivative of the volume flow nonlinear equation (Eq. (6)), the previous set of N second order ODE (Eq. (8)) can be rewritten by using the following expression of :(9)
Then the equations can be put into a statespace representation , where F is a nonlinear vector function, and X the state vector having 2N real components defined as follows:(10)
In practice, because our paper is dedicated to a two quasi harmonic resonance instrument, the state space representation is based on the state vector of four real components , and the nonlinear vector function F can be written as:(11)
Before discussing extensively bifurcation diagrams for different values of inharmonicity and for different configurations of relative amplitudes between Z_{1} and Z_{2} of the two resonances, let us begin by showing and discussing typical elements of a bifurcation diagram. Figure 2 has been obtained by choosing Z_{1} = 1.5Z_{2} and two harmonic resonances (i.e. Inh = 0). The values of the modal parameters of the two resonance’s air column given in Table 1 are inspired from [9] and [19]. The main plot displays the continuation results obtained with AUTO, whereas the six smaller plots above correspond to timedomain simulations of the same system between t = 0 and t = 0.5 s, for different values of γ pointed by numbers. Time integration is performed with an ordinary differential equation solver, namely ode15s from the Matlab ODE Suite.
Figure 2 Bifurcation diagram and time domain simulations of the two harmonic resonance air column (parameter’s values in Tab. 1 with Z_{2} = Z_{1}/1.5 or with Z_{2} = 0, and Inh = 0) with respect to the control parameter γ. Upper plots: six time domain simulations of the dimensionless pressure p = p_{1} + p_{2} calculated between t = 0 and t = 0.5 s for γ = 0.32, γ = 0.36, γ = 0.40 (two simulations with different initial conditions), γ = 0.45 and γ = 0.53. The dimensionless pressure of the plots numbered from 1 to 3 (respectively 4 to 6) is displayed between −0.3 and +0.3 (respectively −1.2 and +1.2). Lower plot: Maximum of the absolute value of the periodic solution p over one period with respect to γ. The branch in green (respectively in black) corresponds to the case Z_{2} = Z_{1}/1.5 (respectively Z_{2} = 0), and illustrates an inverse (respectively direct) Hopf bifurcation scenario. Stable (respectively unstable) solutions are plotted with thick (respectively thin) lines. For each scenario, the Hopf bifurcation point where the equilibrium becomes unstable, is noted γ_{thr}. In the case of an inverse bifurcation the subcrital threshold γ_{subthr} is highlighted with a vertical dashed line. 
Values of the modal parameters of the two resonance’s air column.
The main plot displays maxp, the maximum of the absolute value of pressure in the mouthpiece over one period with respect to the blowing pressure γ. While it is not highlighted here, the horizontal line maxp = 0 corresponds to the equilibrium solution. Below a certain critical value of γ (namely γ < γ_{thr1}), the equilibrium is stable as illustrated by the three time domain simulations calculated for γ = 0.32, γ = 0.36 and γ = 0.40. For initial conditions chosen around the equilibrium, these oscillating solutions decay in time back to the (stable) equilibrium. It is worth noting that the decaying transient lasts all the longer as the value of γ is approaching the critical value γ_{thr1}. When γ = γ_{thr1}, the equilibrium becomes unstable and a branch of periodic solution emerges from the equilibrium. This branch is represented in green on the main plot of Figure 2: it first goes backward in terms of γ and is unstable (thin line), then after a turning point (also called a fold) goes forward and is stable (thick line). This scenario is called an inverse Hopf bifurcation and the value γ = γ_{subthr} the subcritical threshold (see, e.g. [36]).
As explained above, the bifurcation point γ = γ_{thr1} is reached when the real part of one eigenvalue of the jacobian matrix crosses the imaginary axis. The imaginary part of the eigenvalue concerned gives the angular frequency of the emerging periodic solution. In the present case, it is close to ω_{1}. Hence the periodic solution is classified as “first register” or fundamental regime. If the angular frequency of the emerging periodic solution had been close to ω_{2}, the periodic solution would have been classified as “second register” or octave regime. Note that the frequency of the periodic solution along the green branch is not locked at any value but is modified according to the nonlinearity. This is exemplified and discussed in the next section. Two time domain simulations are shown with γ = 0.4 and γ = 0.45 and reveal that the solution is repelled from the equilibrium and converges toward a periodic solution. Note that in the case of γ = 0.4 the choice of the initial condition is crucial since two stable solutions exist: the equilibrium (plot number 3 in Fig. 2) and the periodic solution (plot number 4). A thorough look at the time domain simulation would reveal that maxp deduced from the steadystate (periodic) regime is equal to the ordinate of the green curve at the corresponding value of γ.
The black curve corresponds to emerging branch of periodic solutions in the case where only one acoustic resonance is considered (Z_{2} = 0). In that case, the amplitude maxp is simply a squareroot shaped function of the bifurcation parameter γ in the neighbourhood of the threshold. The thick line denotes a stable periodic solution. Such a scenario is called a direct Hopf bifurcation. Just above the Hopf bifurcation point (γ = γ_{thr1}), the direct scenario leads to stable periodic oscillations with infinitely small amplitudes. Sounds can be played with the nuance pianissimo. On the contrary, in the case of an inverse bifurcation, stable periodic oscillations found just above the Hopf bifurcation point have finite amplitude. Playing with the pianissimo nuance is no longer possible.
For pedagogical purposes, the bifurcation diagram is limited here to the neighborhood of one Hopf bifurcation point, coming from the value γ = γ_{thr1}. However, it will be shown in the next section that for other values of γ, another Hopf bifurcation point is found as well as other bifurcations of the periodic branches.
4 Effects of the inharmonicity. Results and discussion
4.1 Large first resonance amplitude
The discussion is initiated by analysing the case corresponding to Z_{1} slightly higher than Z_{2} (in practice Z_{1}/Z_{2} = 3/2). Three bifurcation diagrams corresponding to Inh = 0, Inh = 0.02 and Inh = 0.04 are shown in Figure 3 (remember that a semi tone corresponds to 0.059).
Figure 3 Three bifurcation diagrams of the two quasiharmonic resonance air column (with Z_{1} > Z_{2}, in practice Z_{1}/Z_{2} = 3/2; parameters values in Tab. 1) with respect to the control parameter γ. From top to bottom: Inh = 0 (a), Inh = 0.02 (b) and Inh = 0.04 (c). Each case is described with two plots. Upper plot: Maximum of the absolute value of the periodic solution with respect to γ. The branch in green (respectively in red, and in blue) corresponds to the fundamental regime, the standard Helmholtz motion (respectively the octave regime, and the inverted Helmholtz motion fundamental regime). Note that a black curve corresponding to a direct Hopf bifurcation is branched at γ = γ_{thr1}. This fundamental regime corresponds to an air column having only one resonance at the frequency f_{res1}. Lower plot: Frequency with respect to γ. The frequency branch in green (respectively in red) corresponds to the fundamental frequency of the fundamental regime (respectively the octave regime, frequency divided by 2). The reference dashed horizontal lines are the reference frequencies: f_{res1} and f_{res2}/2. 
The results shown in Figure 3 for the case Inh = 0 are qualitatively consistent with the one published in [18] (see in particular its Fig. 8). Note that the continuation method gives an additional information: the stability nature of the periodic oscillations.
In Figure 3a the bifurcation diagram shows two branches coming from the equilibrium position:
The first branch originates from the linear threshold γ = γ_{thr1}, associated to the first resonance f_{res1}, originating through an inverse Hopf bifurcation. This fundamental regime, or first register regime, is a standard Helmholtz motion according to [18]. The branch is unstable and then becomes stable at the limit point at γ = γ_{subthr} (subcritical threshold). Compared to the case of a single mode (black curve), important differences are observed, including the nature of the bifurcation.
The second branch originates from the linear threshold γ = γ_{thr2} and is associated to the second resonance f_{res2}, originating through a direct Hopf bifurcation. Note that γ_{thr2} is above γ_{thr1}, because Z_{1} is bigger than Z_{2}. This branch which would correspond to the octave regime, or second register regime, is not observable in practice, because the periodic solutions are unstable.
The nature of the bifurcation of the two branches originating from the linear thresholds γ = γ_{thr1} and γ = γ_{thr2} is in agreement with the publication of [26].
There is a third branch which originates from the unstable octave branch, thanks to a period doubling bifurcation. This branch which would correspond to another fundamental regime (the inverted Helmholtz motion according to [18]) is unstable.
The associated lower plot shows the frequency of the periodic oscillations corresponding to the branches of the bifurcation diagram. In particular the frequency of the fundamental regime (green curve) is almost locked to the value f_{res1} = f_{res2}/2 for any value of γ.
For an inharmonicity of 0.02 (Fig. 3b). The bifurcation diagram is quite close to the one with Inh = 0. However two things are pointed out. First, at the threshold γ = γ_{thr1} the Hopf bifurcation has become direct as it can be predicted theoretically [37]. Second, again there are periodic oscillations for values of the mouth pressure γ under γ = γ_{thr1} until a new value γ = γ_{subthr} which is a bit larger than the one of the case Inh = 0. This is due to the occurence of two folds (limit points on the solution branch) corresponding to saddlenode bifurcations. Note that the frequency of the fundamental regime (green curve) is not locked at the value f_{res1} anymore but is partially pulled towards the value f_{res2}/2, which is reasonable. If the inharmonicity was negative, the same kind of results would have been obtained, the frequency being pulled towards f_{res2}/2 lower than f_{res1}.
For an inharmonicity of 0.04 (Fig. 3c). Now the branch coming from the threshold γ = γ_{thr1}, corresponding to the fundamental regime, looks like a classical branch associated to the direct Hopf bifurcation, there is no γ_{subthr} anymore, since the folds noted in the previous case have disappeared, γ = γ_{thr1} is now the threshold of oscillation. In fact, when the inharmonicity increases, the dynamics of the system behaves more and more like the dynamics of a singleresonance system. The frequency of the fundamental regime comes from the threshold value f_{thr1} at the direct Hopf bifurcation point, and then is partially pulled toward the value f_{res2}/2. Note that in Figure 3a curve corresponding to a direct Hopf bifurcation is branched at γ = γ_{thr1}, this curve corresponds to an air column having only one resonance at the frequency f_{res1}.
Under certain circumstances, for instance when the inharmonicity is high enough, a branch of quasiperiodic solutions may emerge from a Neimark–Sacker bifurcation (often refered as a Hopf bifurcation for a periodic regime). Above this bifurcation point, the periodic branch still exists but it becomes unstable. Such bifurcation has not been encountered in this work, but it has been observed experimentally with a modified saxophone played in the medium range of its tessitura [9], simulated by [38], extensively studied in [19], and it has been studied with continuation on a toy model of saxophone in Section 3 of [39].
The above analysis illustrates significant things because of the inverse Hopf bifurcation (cases Inh = 0 and Inh = 0.02):

On the one hand, there may be a minimum value γ = γ_{subthr} lower than γ_{thr1} above which there are stable periodic oscillations. This particular value γ_{subthr} can be thought of as a quantitative characterisation of the ease of playing. In Figure 3 it is shown that the lowest value of γ_{subthr} is obtained when the two resonances are perfectly harmonic (Inh = 0). If it is assumed that a lower γ_{subthr} corresponds to an instrument easier to play, then it suggests the reed instrument considered is the easiest to play when Inh = 0. In a way that is a theoretical illustration of the BouasseBenade prescription. The threshold of oscillation, equal to γ_{subthr} for low inharmonicities, and equal to γ_{thr}1 for higher inharmonicities, is displayed in Figure 4. The minimum of the threshold of oscillation correspond to Inh = 0.
Figure 4 Minimum value of the mouth pressure γ (green line) corresponding to a stable periodic solution (fundamental regime) with respect to the inharmonicity I_{n}h (case Z_{1} > Z_{2}; in practice Z_{1}/Z_{2} = 3/2). Linear threshold γ_{thr1} (blue dashed line).
On the other hand, the stable periodic oscillations which appear for γ slightly above γ_{subthr} can have fundamental frequencies significantly different from f_{thr1} = f_{res1} because of the effect of the second resonance which controls partially the intonation of the fundamental regime. This study highlights the intrinsic limitation of the linear stability analysis: it should be only considered to assess the stability of the equilibrium. Conclusions concerning the existence of periodic solutions can only be provided through a nonlinear analysis, analitycally in specific cases or with tools like AUTO otherwise.
In addition an animation showing the evolution of the bifurcation diagram as a function of the inharmonicity increasing from Inh = −0.05 to Inh = +0.05 is available from the link of footnote^{1}. Most of the illustrations displayed in the figures are corresponding to a positive inharmonicity Inh, but the animation and the Figure 4 illustrate the fact that the behaviour is qualitavely the same for negative values of Inh.
In order to illustrate the bifurcation diagram (Fig. 3), it is interesting to do simulations by solving the equation in the time domain (sounds available from the links of footnotes^{2}, ^{3}). Figures 5 and 6 show a signal corresponding to an inharmonicity Inh = 0.040 and Z_{1} = 1.5Z_{2}:
In Figure 5 the control parameter increases linearly from γ = 0.43 to γ = 0.50 (crescendo). Because the branch is coming from a direct Hopf bifurcation in the bifurcation diagram, the amplitude of the oscillation is a smoothly increasing mathematical function with respect to γ. Therefore, in the time domain simulation, the amplitude of the signal (fundamental regime) is increasing smoothly, as it is with a resonator having only one resonance f_{res1}.
In Figure 6 the control parameter decreases slowly from γ = 0.55 to γ = 0.50 (decrescendo). Because of the chosen initial conditions, the periodic regime obtained is corresponding to the upper octave, but when γ reaches the value 0.53, the branch coming from γ = γ_{thr2} becomes unstable, and then the periodic solution jumps on the first branch one octave below, the stable branch coming from γ = γ_{thr1} (fundamental regime).
Figure 5 Signal (dimensionless acoustical pressure) obtained by simulation in time domain with an inharmonicity Inh = 0.04 and Z_{1} = 1.5Z_{2} (like in Fig. 3c). The dimensionless mouth pressure is printed in black, and is increasing linearly from γ = 0.43 (constant before t = 2 s) to γ = 0.50 (constant after t = 9 s). 
Figure 6 Signal (dimensionless acoustical pressure) obtained by simulation in time domain with an inharmonicity Inh = 0.04 and Z_{1} = 1.5Z_{2} (like in Fig. 3c). The dimensionless mouth pressure is printed in black, and is decreasing linearly from γ = 0.55 (constant before t = 2 s) to γ = 0.50 (constant after t = 9 s). 
4.2 Large second resonance amplitude
The discussion continues by analysing the case corresponding to Z_{1} slightly lower than Z_{2} (in practice Z_{2}/Z_{1} = 3/2). Three bifurcation diagrams corresponding to Inh = 0, Inh = 0.015 and Inh = 0.03 are shown Figure 7.
Figure 7 Three bifurcation diagrams of the two quasiharmonic resonance air column (with Z_{1} < Z_{2}, in practice Z_{2}/Z_{1} = 3/2) with respect to the control parameter γ. From top to bottom: Inh = 0 (a), Inh = 0.015 (b) and Inh = 0.030 (c). Each case is described with two plots. Upper plot (a): Maximum of the absolute value of the periodic solution with respect to γ. The branch in red (respectively in blue, and in green) corresponds to the octave regime (respectively the fundamental regime associated to the standard Helmholtz motion, and the inverted Helmholtz motion fundamental regime). Lower plot (c): Frequency with respect to γ. The frequency branch in red (respectively in blue) corresponds to the fundamental frequency of the octave regime (respectively the fundamental regime). The frequency of the red branch has been divided by 2 for sake of clarity. The reference dashed horizontal lines are the reference frequencies: f_{res1} and f_{res2}/2. 
The results shown in Figure 7 for the case Inh = 0 are qualitatively consistent with the one published in [18] (see in particular its Fig. 10).
In Figure 7a the bifurcation diagram shows two branches coming from the equilibrium position:
On the lefthand side, the first branch originates from the linear threshold γ = γ_{thr2}, associated to the second resonance f_{res2}, originating through a direct bifurcation. This octave regime is stable until a period doubling bifurcation point, and then becomes unstable. At the bifurcation point, there is an emerging branch corresponding to a fundamental regime. It is a standard Helmholtz motion according to [18]. This fundamental regime is unstable until a turning point (a fold) corresponding to a minimum value of γ = γ_{subthr} where the periodic oscillations become stable. Note that the threshold of oscillation of the fundamental regime γ = γ_{subthr} is significantly lower than the value γ_{thr1} predicted by the linear stability analysis.
The second branch originates from the linear threshold γ = γ_{thr1}, associated to the first resonance f_{res1}, originating through a direct bifurcation. Note that γ_{thr1} is bigger than γ_{thr2}, because Z_{1} is lower than Z_{2}. This branch which would correspond to a second fundamental regime is not observable in practice, because the periodic solutions are unstable. This branch would correspond to the inverted Helmholtz motion according to [18].
The associated lower curve shows the frequency of the periodic oscillations corresponding to the branches of the bifurcation diagram. In particular the frequency of the stable fundamental regime (blue curve) is close to the value f_{res1} = f_{res2}/2 for any value of γ.
For an inharmonicity of 0.015 (Fig. 7b). The bifurcation diagram is qualitatively quite close to the one with Inh = 0. Two things are now pointed out. Once again there are periodic oscillations for dimensionless mouth pressure γ values below γ = γ_{thr2} < γ_{thr1} until a new value γ = γ_{subthr} which is a bit bigger than the one in the case of Inh = 0. Note that the frequency of the fundamental regime (blue curve) is surprisingly close to the value f_{res1}, the fundamental frequency is not much pulled towards the value f_{res2}/2. Note that the range of γ where there are two stable periodic regimes, octave and the standard Helmholtz motion fundamental regime, is larger: from γ = γ_{thr2} to the value of γ where the period doubling bifurcation point occurs.
For an inharmonicity of 0.03 (Fig. 7c). Again the bifurcation diagram is qualitatively quite close to the ones with Inh = 0 and Inh = 0.015. The minimum pressure of fundamental periodic oscillations γ_{subthr} (on the blue branch keeps increasing with inharmonicity, and becomes higher than γ = γ_{thr2}.
The above discussion illustrates significant things:

There may be a minimum value γ = γ_{subthr} lower than γ_{thr2} < γ_{thr1} where there are stable periodic oscillations. This particular value γ_{subthr} can be chosen as a kind of quantitative characterisation of the ease of playing. In Figure 8 it is shown again (as in Fig. 4) that the lowest value of γ_{subthr} is obtained when the two resonances are perfectly harmonic (Inh = 0). If it is assumed that a lower γ_{subthr} corresponds to an instrument easier to play, then it suggests that the reed instrument considered to be the easiest to play when Inh = 0. In a way, even if Z_{1} < Z_{1}, again that is a theoretical illustration of the BouasseBenade prescription. The threshold of oscillation is displayed in Figure 8: the minimum of the threshold of oscillation is corresponding to Inh = 0.
Figure 8 Minimum value of the mouth pressure γ (blue line) corresponding to a stable periodic solution (fundamental regime) with respect to the inharmonicity Inh (case Z_{1} < Z_{2}; in practice Z_{2}/Z_{1} = 3/2). Linear threshold γ_{thr1} (green dashed line).
The stable periodic oscillations which appear for γ slightly bigger than γ_{subthr} have fundamental frequencies quite close to f_{res1}.
It is worth emphasising that, whatever the inharmonicity, the fundamental regime does never come from the first threshold γ = γ_{thr1}, but comes through a perioddoubling bifurcation point attached to the octave branch. A naive analysis of the time domain simulations (at least with Inh = 0) would probably suggest that the fundamental regime emerges from the equilibrium trough an inverse Hopf bifurcation, but this is not correct. It is also worth noting that a linear stability analysis (LSA) of the equilibrium is useless here to give some hints about the oscillation behaviour of the model.
Note as well that sometimes, there are several stable regimes (equilibrium position and periodic regime, or two periodic regimes) for a given value of γ. For such cases, the stable regime reached is the consequence of the initial conditions.
In addition, an animation showing the evolution of the bifurcation diagram as a function of the inharmonicity increasing from Inh = −0.05 to Inh = +0.05 is available from the link of footnote^{4}. Most of the illustrations displayed in the figures correspond to a positive inharmonicity Inh, but the animation and the Figure 8 illustrate the fact that the behaviour is qualitavely the same for negative values of Inh. Unlike Figure 4, it can be noted that the plot is slightly asymetric with respect to the vertical axis Inh = 0.
In order to illustrate the bifurcation diagrams (Fig. 7), it is interesting to do simulations by solving the equation in the time domain (sound available from the link of footnote^{5}). Figure 9 shows a signal corresponding to an inharmonicity Inh = 0.015 and Z_{1} = Z_{2}/1.5. The control parameter increases slowly from γ = 0.38 (just below the period doubling bifurcation) to γ = 0.45 (crescendo). Because of the chosen initial conditions, the periodic regime obtained corresponds to the octave, but when γ reaches the value 0.39, the branch coming from γ = γ_{th2} becomes unstable, and then the periodic solution jumps to the only stable branch one octave below (branch coming from the perioddoubling bifurcation).
Figure 9 Signal (dimensionless acoustical pressure) obtained by simulation in time domain with an inharmonicity Inh = 0.015 and Z_{1} = Z_{2}/1.5 (like in Fig. 7b). The dimensionless mouth pressure is printed in black, and is increasing linearly from γ = 0.38 (constant before t = 2 s) to γ = 0.45 (constant after t = 9 s). 
5 Conclusion
Bifurcation diagrams of a basic reed instrument modeled by two quasiharmonic resonances have been computed by using a continuation method (AUTO software), where the mouth pressure is the control parameter. Some of the mouth pressure thresholds results are interpreted in terms of the ease of playing of the reed instrument. When there is an inverse Hopf bifurcation (perfect harmonicity) or of a double fold after a direct Hopf bifurcation (moderate inharmonicity), there may be a minimum value γ = γ_{subthr} lower than γ_{thr1} for which periodic stable oscillations can be observed. This value γ_{subthr} may be considered as a quantitative characterisation of the ease of playing. It has been shown that the lowest value of γ_{subthr} is obtained when the two resonances are harmonic, harmonicity equal to 2. This is a theoretical illustration of the BouasseBenade prescription [2, 8]). Even if a few AUTO simulations using other parameter’s values than the one used in the present study have been done, a large set of other tests should be done with many other parameter’s values to verify that the conlusions of the present paper are robust.
An interesting direction for future work could include experimental validation, particularly using the modified saxophone used in [9]. There, the saxophone was modified by the addition of two closed side tubes on the neck. Movable pistons are used to change the volume of the side tubes, which results in a shift of the resonant frequencies. As explained in this publication, it is possible to choose particular positions and volumes of the closed tubes to ensure a control of the inharmonicity. This strategy is used with fingerings corresponding to the middle and high ranges of the first regime of the saxophone, where its input impedance consists essentially of two resonances.
The results provided in the current manuscript depend on a physical model of reed instruments based on three strong approximations: the reed dynamics is ignored, only two acoustic resonances are taken into account, and the nonlinear equation describing the incoming volume flow is approximated by its third order Taylor series expansion. Therefore, the conclusions above cannot be directly extended to real instruments until further research is carried out on more complex models. At that point, many other interesting topics could be explored, such as sound production of lowpitched notes by conical reed instruments such as saxophones, oboes or bassoons. Replacing the inwardstriking reed model by an outwardstriking lip model is also planned in order to study nonlinear dynamics of brass instruments (preliminary results in [40]). More precisely, it is expected that numerical continuation could clarify their pedal note regime, recently simulated using a timedomain finite differences method in [41].
Conflict of interest
The authors declare no conflict of interest.
Acknowledgments
Authors acknowledge their colleagues Murray Campbell, JeanPierre Dalmont and Erik Petersen for fruitfull discussions. And authors wish to thank the reviewers, whose conscientious work has improved the clarity of this article and raised some interesting questions.
Animation showing the evolution of the bifurcation diagram as a function of the inharmonicity increasing from Inh = −0.05 to Inh = +0.05 (case Z_{1} > Z_{2}), corresponding to Figure 3, in http://perso.univlemans.fr/~jgilbert/output_Z1_sup_Z2_stab.webm.
Time domain simulation corresponding to Figure 5 in http://perso.univlemans.fr/~jgilbert/Inh0p040_10s.wav.
Time domain simulation corresponding to Figure 6 in http://perso.univlemans.fr/~jgilbert/gamma0p55a0p50.wav.
Animation showing the evolution of the bifurcation diagram as a function of the inharmonicity increasing from Inh = −0.05 to Inh = +0.05 (case Z_{1} < Z_{2}), corresponding to Figure 7, in http://perso.univlemans.fr/~jgilbert/output_Z2_sup_Z1_stab.webm.
Time domain simulation corresponding to Figure 9 in http://perso.univlemans.fr/~jgilbert/Inh0p015_Z2supZ1_10s.wav.
References
 J.F. Petiot, P. Kersaudy, G. Scavone, S. Mac Adams, B. Gazengel: Investigations of the relationships between perceived qualities and sound parameters of saxophone reeds. Acustica United With Acta Acustica 103 (2017) 812–829. [CrossRef] [Google Scholar]
 A.H. Benade: Fundamentals of musical acoustics, 2nd ed. Dover, 1990. [Google Scholar]
 M. Campbell, C. Greated: The Musician’s Guide to Acoustics. Oxford University Press, 1989. [Google Scholar]
 N.H. Fletcher, T.D. Rossing: The Physics of Musical Instruments, 2nd ed. Springer, 1998. [CrossRef] [Google Scholar]
 A. Chaigne, J. Kergomard: Acoustics of Musical Instruments. Springer, 2016. [CrossRef] [Google Scholar]
 M. Campbell, J. Gilbert, A. Myers: The Science of Brass Instruments. Springer, 2020. [Google Scholar]
 A.H. Benade, D.J. Gans: Sound production in wind instruments. Annals of the New York Academy of Science 155 (1968) 247–263. [CrossRef] [Google Scholar]
 H. Bouasse: Instruments à vent tomes I et II. Delagrave, Paris, 1929; repr. Librairie Scientifique et Technique Albert Blanchard, Paris, 1986. [Google Scholar]
 J.P. Dalmont, B. Gazengel, J. Gilbert, J. Kergomard: Some aspects of tuning and clean intonation in reed instruments. Applied Acoustics 461 (1995) 19–60. [CrossRef] [Google Scholar]
 J. Gilbert, E. Brasseur, J.P. Dalmont, C. Maniquet: Acoustical evaluation of the Carnyx of Tintignac. Proceedings of Acoustics 2012, Nantes, 2012. [Google Scholar]
 D.M. Campbell, J. Gilbert, P. Holmes: Seeking the sound of ancient horns. ASA Meeting, Boston, 2017. [Google Scholar]
 W. Kausel: Optimization of brasswind instruments and its application in bore reconstruction. Journal of New Music Research 30 (2001) 69–82. [CrossRef] [Google Scholar]
 A. Braden, M. Newton, D.M. Campbell: Trombone bore optimization based on input impedance targets. Journal of the Acoustical Society of America 125 (2009) 2404–2412. [CrossRef] [Google Scholar]
 D. Noreland, J. Kergomard, F. Laloë, C. Vergez, P. Guillemain, A. Guilloteau: The logical clarinet: Numerical optimization of the geometry of woodwind instruments. Acta Acustica United With Acustica 99 (2013) 615–628. [CrossRef] [Google Scholar]
 W.L. Coyle, P. Guillemain, J. Kergomard, J.P. Dalmont: Predicting playing frequencies for clarinets: A comparison between numerical simulations and simplified analytical formulas. Journal of the Acoustical Society of America 138 (2015) 2770–2781. [CrossRef] [Google Scholar]
 R. Tournemenne, J.F. Petiot, B. Talgorn, M. Kokkolaras, J. Gilbert: Sound simulation based design optimization of brass wind instruments. Journal of the Acoustical Society of America 145 (2019) 3795–3804. [CrossRef] [Google Scholar]
 M.E. McIntyre, R.T. Schumacher, J. Woodhouse: On the oscillations of musical instruments. Journal of the Acoustical Society of America 74 (1983) 1325–1345. [CrossRef] [Google Scholar]
 J.P. Dalmont, J. Gilbert, J. Kergomard: Reed instruments, from small to large amplitude periodic oscillations and the Helmholtz motion analogy. Acustica 86 (2000) 671–684. [Google Scholar]
 J.B. Doc, C. Vergez, S. Missoum: A minimal model of a singlereed instrument producing quasiperiodic sounds. Acta Acustica United With Acustica 100 (2014) 543–554. [CrossRef] [Google Scholar]
 A. Hirschberg, J. Kergomard, G. Weinreich: Mechanics of musical instruments. SpringerVerlag, Wien, Austria, 1995. [Google Scholar]
 J. Gilbert, S. Maugeais, C. Vergez: From the bifurcation diagrams to the ease of playing of reed musical instruments. A theoretical illustration of the BouasseBenade prescription? International Symposium on Musical Acoustics, Detmold, Germany, 2019. [Google Scholar]
 B. Fabre, J. Gilbert, A. Hirschberg: Modeling of Wind Instruments. Chapter 7 of Springer Handbook of Systematic Musicology. SpringerVerlag, 2018. [Google Scholar]
 T.A. Wilson, G.S. Beavers: Operating modes of the clarinet. Journal of the Acoustical Society of America 56 (1974) 653–658. [Google Scholar]
 F. Silva, J. Kergomard, C. Vergez, J. Gilbert: Interaction of reed and acoustic resonator in clarinetlike systems. Journal of the Acoustical Society of America 124 (2008) 3284–3295. [CrossRef] [PubMed] [Google Scholar]
 W.E. Worman: Selfsustained nonlinear oscillations of medium amplitude in clarinetlike systems, PhD Thesis. Case Western Reserve University, Cleveland, 1971. [Google Scholar]
 N. Grand, J. Gilbert, F. Laloë: Oscillation threshold of woodwind instruments. Acustica 83 (1997) 137–151. [Google Scholar]
 B. Ricaud, P. Guillemain, J. Kergomard, F. Silva, C. Vergez: Behavior of reed woodwind instruments around the oscillation threshold. Acta Acustica United With Acustica 95 (2009) 733–743. [CrossRef] [Google Scholar]
 C. Maganza, R. Caussé, F. Laloë: Bifurcations, period doubling and chaos in clarinet like systems. Europhysics Letters 1 (1986) 295–302. [CrossRef] [Google Scholar]
 P.A. Taillard, J. Kergomard, F. Laloë: Iterated maps for clarinetlike systems. Nonlinear Dynamics 62 (2010) 253–271. [CrossRef] [Google Scholar]
 J.P. Dalmont, J. Gilbert, J. Kergomard, S. Ollivier: An analytical prediction of the oscillation and extinction thresholds of a clarinet. Journal of the Acoustical Society of America 118 (2005) 3294–3305. [CrossRef] [PubMed] [Google Scholar]
 E.J. Doedel, A.R. Champneys, T.F. Fairgrieve, Yu.A. Kuznetsov, B. Sandstede, X.J. Wang: auto97: Continuation and bifurcation software for ordinary differential equations (with HomCont) User’s Guide. Concordia Univ. (1997). [Google Scholar]
 S. Karkar, B. Cochelin, C. Vergez: A highorder, purely frequency based harmonic balance formulation for continuation of periodic solutions: The case of nonpolynomial nonlinearities. Journal of Sound and Vibration 332 (2013) 968–977. [CrossRef] [Google Scholar]
 http://indy.cs.concordia.ca/auto/ [Google Scholar]
 V. Debut, J. Kergomard: Analysis of the selfsustained oscillations of a clarinet as a Van der Pol oscillator. International Congress on Acoustics, Kyoto, 2004. [Google Scholar]
 F. Silva, C. Vergez, P. Guillemain, J. Kergomard, V. Debut: MoReeSC: a framework for the simulation and analysis of sound production in reed and brass instruments. Acta Acustica United With Acustica 100 (2014) 126–138. [CrossRef] [Google Scholar]
 S.H. Strogatz: Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering (Studies in Nonlinearity), 2nd ed., Kindle, 2019. [Google Scholar]
 B. Gazengel: Caractérisation objective de la qualité de justesse, de timbre et d’émission des instruments à vent à anche simple, PhD Thesis. Université du Maine, 1994. [Google Scholar]
 B. Gazengel, J. Gilbert: From the measured input impedance to the synthesized pressure signal: application to the saxophone. Proceedings of the International Symposium on Musical Acoustics, Dourdan, July 2–6, 1995. [Google Scholar]
 L. Guillot, B. Cochelin, C. Vergez: A Taylor seriesbased continuation method for solutions of dynamical systems. Nonlinear Dynamics 98 (2019) 2827–2845. [CrossRef] [Google Scholar]
 V. Freour, H. Masuda, S. Usa, E. Tominaga, Y. Tohgi, B. Cochelin, C. Vergez: Numerical analysis and comparison of brass instruments by continuation. International Symposium on Musical Acoustics, Detmold, Germany, 2019. [Google Scholar]
 L. Velut, C. Vergez, J. Gilbert, M. Djahanbani: How well can Linear Stability Analysis predict the behaviour of an outward valve brass instrument model? Acustica United With Acta Acustica 103 (2016) 132–148. [CrossRef] [Google Scholar]
Cite this article as: Gilbert J, Maugeais S & Vergez C. 2020. Minimal blowing pressure allowing periodic oscillations in a simplified reed musical instrument model: BouasseBenade prescription assessed through numerical continuation. Acta Acustica, 4, 27.
All Tables
All Figures
Figure 1 Diagram showing the regions where the bifurcation is direct as well as those regions where it is inverse. The xaxis shows the values of the third coefficient of the Taylor expansion Equation (6), and the yaxis is the difference between the admittance amplitude between the two first resonances (assumed to be harmonic) Y_{2} −Y_{1} = 1/Z_{2} − 1/Z_{1}. The hatched region is for a direct bifurcation, and the unhatched region for an inverse bifurcation. Adapted from [26]. 

In the text 
Figure 2 Bifurcation diagram and time domain simulations of the two harmonic resonance air column (parameter’s values in Tab. 1 with Z_{2} = Z_{1}/1.5 or with Z_{2} = 0, and Inh = 0) with respect to the control parameter γ. Upper plots: six time domain simulations of the dimensionless pressure p = p_{1} + p_{2} calculated between t = 0 and t = 0.5 s for γ = 0.32, γ = 0.36, γ = 0.40 (two simulations with different initial conditions), γ = 0.45 and γ = 0.53. The dimensionless pressure of the plots numbered from 1 to 3 (respectively 4 to 6) is displayed between −0.3 and +0.3 (respectively −1.2 and +1.2). Lower plot: Maximum of the absolute value of the periodic solution p over one period with respect to γ. The branch in green (respectively in black) corresponds to the case Z_{2} = Z_{1}/1.5 (respectively Z_{2} = 0), and illustrates an inverse (respectively direct) Hopf bifurcation scenario. Stable (respectively unstable) solutions are plotted with thick (respectively thin) lines. For each scenario, the Hopf bifurcation point where the equilibrium becomes unstable, is noted γ_{thr}. In the case of an inverse bifurcation the subcrital threshold γ_{subthr} is highlighted with a vertical dashed line. 

In the text 
Figure 3 Three bifurcation diagrams of the two quasiharmonic resonance air column (with Z_{1} > Z_{2}, in practice Z_{1}/Z_{2} = 3/2; parameters values in Tab. 1) with respect to the control parameter γ. From top to bottom: Inh = 0 (a), Inh = 0.02 (b) and Inh = 0.04 (c). Each case is described with two plots. Upper plot: Maximum of the absolute value of the periodic solution with respect to γ. The branch in green (respectively in red, and in blue) corresponds to the fundamental regime, the standard Helmholtz motion (respectively the octave regime, and the inverted Helmholtz motion fundamental regime). Note that a black curve corresponding to a direct Hopf bifurcation is branched at γ = γ_{thr1}. This fundamental regime corresponds to an air column having only one resonance at the frequency f_{res1}. Lower plot: Frequency with respect to γ. The frequency branch in green (respectively in red) corresponds to the fundamental frequency of the fundamental regime (respectively the octave regime, frequency divided by 2). The reference dashed horizontal lines are the reference frequencies: f_{res1} and f_{res2}/2. 

In the text 
Figure 4 Minimum value of the mouth pressure γ (green line) corresponding to a stable periodic solution (fundamental regime) with respect to the inharmonicity I_{n}h (case Z_{1} > Z_{2}; in practice Z_{1}/Z_{2} = 3/2). Linear threshold γ_{thr1} (blue dashed line). 

In the text 
Figure 5 Signal (dimensionless acoustical pressure) obtained by simulation in time domain with an inharmonicity Inh = 0.04 and Z_{1} = 1.5Z_{2} (like in Fig. 3c). The dimensionless mouth pressure is printed in black, and is increasing linearly from γ = 0.43 (constant before t = 2 s) to γ = 0.50 (constant after t = 9 s). 

In the text 
Figure 6 Signal (dimensionless acoustical pressure) obtained by simulation in time domain with an inharmonicity Inh = 0.04 and Z_{1} = 1.5Z_{2} (like in Fig. 3c). The dimensionless mouth pressure is printed in black, and is decreasing linearly from γ = 0.55 (constant before t = 2 s) to γ = 0.50 (constant after t = 9 s). 

In the text 
Figure 7 Three bifurcation diagrams of the two quasiharmonic resonance air column (with Z_{1} < Z_{2}, in practice Z_{2}/Z_{1} = 3/2) with respect to the control parameter γ. From top to bottom: Inh = 0 (a), Inh = 0.015 (b) and Inh = 0.030 (c). Each case is described with two plots. Upper plot (a): Maximum of the absolute value of the periodic solution with respect to γ. The branch in red (respectively in blue, and in green) corresponds to the octave regime (respectively the fundamental regime associated to the standard Helmholtz motion, and the inverted Helmholtz motion fundamental regime). Lower plot (c): Frequency with respect to γ. The frequency branch in red (respectively in blue) corresponds to the fundamental frequency of the octave regime (respectively the fundamental regime). The frequency of the red branch has been divided by 2 for sake of clarity. The reference dashed horizontal lines are the reference frequencies: f_{res1} and f_{res2}/2. 

In the text 
Figure 8 Minimum value of the mouth pressure γ (blue line) corresponding to a stable periodic solution (fundamental regime) with respect to the inharmonicity Inh (case Z_{1} < Z_{2}; in practice Z_{2}/Z_{1} = 3/2). Linear threshold γ_{thr1} (green dashed line). 

In the text 
Figure 9 Signal (dimensionless acoustical pressure) obtained by simulation in time domain with an inharmonicity Inh = 0.015 and Z_{1} = Z_{2}/1.5 (like in Fig. 7b). The dimensionless mouth pressure is printed in black, and is increasing linearly from γ = 0.38 (constant before t = 2 s) to γ = 0.45 (constant after t = 9 s). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.