Issue 
Acta Acust.
Volume 5, 2021



Article Number  57  
Number of page(s)  12  
Section  Musical Acoustics  
DOI  https://doi.org/10.1051/aacus/2021049  
Published online  22 December 2021 
Scientific Article
Minimal blowing pressure allowing periodic oscillations in a model of bass brass instruments
^{1}
Laboratoire d’Acoustique de l’Université du Mans (LAUM), UMR CNRS 6613, Institut d’Acoustique  Graduate School (IAGS), CNRS, Le Mans Université, 72085 Le Mans, France
^{2}
Aix Marseille Univ, CNRS, Centrale Marseille, LMA, UMR 7031, 13625 AixenProvence Cedex, France
^{3}
Laboratoire Manceau de Mathématiques – Le Mans Université, 72085 Le Mans, France
^{*} Corresponding author: remi.matteoli@univlemans.fr
Received:
6
September
2021
Accepted:
23
November
2021
In this study, an acoustic resonator – a bass brass instrument – with multiple resonances coupled to an exciter – the player’s lips – with one resonance is modelled by a multidimensional dynamical system, and studied using a continuation and bifurcation software. Bifurcation diagrams are explored with respect to the blowing pressure, in particular with focus on the minimal blowing pressure allowing stable periodic oscillations and the associated frequency. The behaviour of the instrument is first studied close to a (non oscillating) equilibrium using linear stability analysis. This allows to determine the conditions at which an equilibrium destabilises and as such where oscillating regimes can emerge (corresponding to a sound production). This approach is useful to characterise the ease of playing of a brass instrument, which is assumed here to be related – as a first approximation – to the linear threshold pressure. In particular, the lower the threshold pressure, the lower the physical effort the player has to make to play a note [The Science of Brass Instruments. SpringerVerlag, 2021]. Cases are highlighted where periodic solutions in the bifurcation diagrams are reached for blowing pressures below the value given by the linear stability analysis. Thus, bifurcation diagrams allow a more indepth analysis. Particular attention is devoted to the first playing regime of bass brass instruments (the pedal note and the ghost note of a tuba in particular), whose behaviour qualitatively differs from a trombone to a euphonium for instance.
Key words: Musical acoustics / Brass instruments / Ease of playing / Dynamical systems / Bifurcation diagrams
© R. Mattéoli et al., Published by EDP Sciences, 2021
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
One main goal of the acoustics of wind instruments is to describe and quantify the intonation and ease of playing of an instrument. From the physics viewpoint, it is interesting to model the coupled system formed by the musician and the instrument. Of particular interest is the influence of the musician’s control parameters on the oscillation frequency (linked to the intonation), and the minimum mouth pressure required to achieve autooscillations (related to the ease of playing). Indeed, it is assumed here that the musician’s feeling of ease of playing partly relies on the threshold blowing pressure: the higher the latter, the higher the physical effort the player has to make to play a note [1]. In practice, the musician can play several notes without depressing any valves in the case of a tuba, or moving the slide in the case of a trombone. These playing regimes are called the natural notes (B♭1, B♭2, F3, B♭3, D4, F4,... in the case of a trombone or a euphonium for instance), and their frequencies are close to the resonance frequencies of the instrument as a whole, except for the lowest note playable (B♭1).
Several tools are available to infer the oscillation frequencies and the threshold blowing pressures from a model of the system. Linear stability analysis allows to predict the behaviour of a dynamical system in the vicinity of its equilibrium solutions. The method consists in linearising the system and analysing its eigenvalues to determine whether or not the equilibrium solutions are stable. This method has already been applied to physical models of musical instruments – brass instruments – in [2]. Linear stability analysis is, in essence, only valid close to the equilibrium solutions. As such, this method alone might not explain the existence of periodic regimes, as highlighted for instance in [2] for the first regime of a saxhorn.
Alternatively, bifurcation analysis gives access to an extensive knowledge of permanent regimes (namely equilibrium and periodic regimes) of the dynamical system, as well as their stability when relevant. Bifurcation diagrams ideally represent all the families of equilibrium and periodic solutions with respect to one parameter of interest. In this paper, we investigate bifurcation diagrams of a physical model of a brass instrument to understand some aspects of its behaviour.
We consider here a simple model: the player’s lips are modeled through the “outwardstriking valve” [1, 3, 4], a one degreeoffreedom system. Also, the nonlinear sound propagation in the instrument’s bore which accounts for the “brassy sounds” [1] at high sound levels is neglected [5, 6]. The airflow blown by the musician into the instrument is described by a nonlinear algebraic equation. Some of the functions in this equation being nonsmooth, a regularised form of this equation is used for numerical reasons.
In Section 2, we describe the considered model and numerical methods, namely linear stability analysis and numerical continuation. In Sections 3 and 4, we use linear stability analysis and bifurcation diagrams, respectively, to investigate the minimal blowing pressure and playing frequencies of the trombone. Eventually, we focus in Section 5 on the first regime of two bass brass instruments (the trombone and the euphonium), mainly thanks to bifurcation diagrams. This highlights a major difference in the behaviour of these two instruments.
2 Theoretical background
2.1 Generic brass model
This subsection details the brassinstrument model considered throughout the article. Brass instruments as a whole can be described through both linear and nonlinear mechanisms. More precisely, a localised nonlinear element (the lips’ valve effect, i.e. the airflow modulation caused by the lips’ vibration) excites a passive linear acoustic multimode element (the musical instrument, usually characterised by its input impedance in the frequency domain). The latter acoustic resonator exerts, in turn, a retroaction on the former mechanical resonator. Such musical instruments are selfsustained oscillators: they generate an oscillating acoustic pressure (the note played) from a static overpressure in the player’s mouth (the blowing pressure).
The brass instrument coupled to the player is then described by a system of three equations. Because it relies on major simplifications [1, 4], this model is often referred to as “elementary”. More precisely, the three equations link the lipopening height h(t), the pressure in the mouthpiece p(t) and the volume airflow entering the instrument u(t), which will be the three independent variables of interest in this paper. They are all represented schematically in Figure 1. Several control parameters (i.e. controlled by the musician) are also involved: this includes in particular the blowing mouth pressure p_{m} and the lips’ resonance frequency .
Figure 1 Schematic representation of the player’s mouth coupled to the mouthpiece of a brass instrument. 
First, the vibrating lips of the musician are described by a onedegreeoffreedom damped oscillator:
where and are the angular resonance frequency and the quality factor of the lips, respectively, μ is the lips’ mass per unit area, and H is the lipopening height at rest.
The Bernoulli’s theorem applied between the mouth and the mouthpiece and taking into account the pressure drop caused by the presence of turbulence in the mouthpiece leads to the following equation:
with w the lipopening width (considered as constant) and ρ the air density. Here, h^{+} = max(h, 0) accounts for the fact that the lips cannot physically interpenetrate: as soon as the lips touch (h = 0), the volume flow is forced to zero. The sign function sgn accounts for the possibility of air flowing from the instrument into the player’s mouth.
Finally, the acoustic input impedance Z(ω) of the resonator is described in the Fourier domain as the ratio between the pressure and the volume flow at the input of the instrument. This provides a last link between the mouth pressure and the volume airflow:
with ω the angular frequency. Figure 2 shows the modulus and phase of Z with respect to the angular frequency ω, for a bass trombone (A. Courtois Legend AC502).
Figure 2 Modulus (top) and phase (bottom) of the input impedance of a bass trombone (A. Courtois Legend AC502) with respect to frequency. Blue curve: measured impedance; dashed orange curve: modal fit function with N = 15 modes; green curve: error between measured and fitted impedances, defined by Z_{fit} − Z_{exp} and arg(Z_{fit}) − arg(Z_{exp}) respectively as regards the module and phase. 
Several limitations of this generic model have already been pointed out in previous works [1, 2]. Firstly, it is worth noting that the instrument is described with a linear model of its input impedance: the nonlinear propagation of sound in the instrument itself is not taken into account, preventing any description of the characteristic “brassy sound” of a brass instrument at high sound levels. However, since we focus here on the emergence of autooscillations, this effect can be neglected. Furthermore, the harmonics created by the nonlinear distortion of the wavefront propagating in the instrument are known not to be reflected back to the player when they reach the bell [7].
Secondly, lips are described as a onedegreeoffreedom damped oscillator, which is also a strong assumption. As a matter of fact, both experiments using artificial lips [3] and invivo measurements [8] demonstrated that the lips’ oscillation mechanism relies mainly on a pair of lips’ resonance frequencies. Furthermore, only the first resonance of the lips reflects the “outwardstriking swinging door” mechanism described in [1], hence the single degree of freedom. In any case, the onedegreeoffreedom damped oscillator prevents any description of the “buzzing phenomenon”, which corresponds to the lips’ vibration without the instrument.
Thirdly, the lips’ opening area S in the airflow (Eq. (2)) is assumed to be proportional to h, and is thus written S = w × h with w a constant. This assumption was only partially validated in [9], depending on which regime is studied as well as on the sound level. In particular, a law S ∝ h^{2} would be more adequate when studying the first regimes like the pedal note.
2.2 Numerical considerations
In this paper, we aim at determining the minimal blowing pressures for which each periodic regime (natural notes B♭1, B♭2, F3, B♭3, D4, F4, …) is observable, as well as the corresponding playing frequencies.
In this subsection, the numerical methods allowing for the investigation of the periodic solutions of model {(1) ∪ (2) ∪ (3)} are presented, as well as the modal representation of the input impedance required for the practical implementation of these methods.
2.2.1 Numerical methods
Several numerical methods are available to investigate the influence of a control parameter – such as p_{m} or , as they are the most obvious parameters the player can change – on the dynamics of the system. Linear stability analysis [2] consists in studying the stability of the equilibrium solutions to small perturbations, obtained by zeroing all time derivatives. This allows one to determine the threshold value of the control parameter at which the equilibrium solution destabilises, which means that at least one of the eigenvalues of the system exhibits a positive real part, and where a smallamplitude periodic regime (stable or unstable) can emerge. When a smallamplitude periodic regime indeed emerges, such a change in the system’s behaviour is called a Hopf bifurcation [10].
In the case of brass instruments, both p_{m} and can be considered as control parameters. Figure 4 (blue curve) represents the threshold pressure predicted by the linear stability analysis with respect to . This analysis also provides the threshold frequency, that is to say the frequency of the periodic regime that emerges when the equilibrium solution destabilises. This is represented in the bottom plot in Figure 4.
One of the main advantages of this method is its straightforward implementation. However, it gives very little information on the oscillation regime of the model far from the socalled Hopf bifurcation point at which the equilibrium loses its stability. This method has already been used in [2], and is applied here to the trombone (see blue curves in Fig. 4).
To know more about the oscillating solution arbitrarily far from the bifurcation point (i.e. the point at which the equilibrium solution becomes unstable), another approach consists in numerically solving the whole system, thanks to an ODE solver [2]. In doing so, both the transient and the stationary parts of the solution are obtained for any value of the control parameter. Nevertheless, this approach becomes tedious and unsuitable in the context of the systematic investigation of the influence of the control parameter on the oscillation regime.
Continuation methods, on the other hand, are more suitable to gain access to a more extensive view of all the oscillating regimes of the system. It consists in computing the waveform of the oscillating solution for successive values of the control parameter. The waveform corresponding to a new value of the control parameter is then deduced from previously computed waveforms, through a predictor/corrector algorithm. The behaviour of an oscillatory solution of the system with respect to a control parameter such as p_{m} can then be assessed by plotting bifurcation diagrams, which are shown in Section 4. This approach is implemented in several softwares, such as AUTO [11–15] which is used in this publication. However, continuation methods require the system to be written in the form , with certain smoothness properties on F. Therefore, some work has yet to be done on the equations presented in Section 2.1, which is done in the following.
2.2.2 Input impedance
Since the continuation method requires the equations to be written in the time domain, equation (3) cannot be used as it is. An analytical form of the input impedance Z is therefore required to perform an inverse Fourier transform of (3). This quantity can be quite easily measured, and is represented in Figure 2 (blue curve). This measured impedance can then be numerically fitted, in the frequency domain, by a sum of N individual acoustical resonance modes of the following form [16]:
Here, are the modal parameters of the modal decomposition. For numerical reasons, it has been chosen to convert this set of parameters (A_{n}, ω_{n}, ξ_{n}) into another one (C_{n}, s_{n}) so that the fitted input impedance writes (z* is the complex conjugate of z):
where Z_{c} is the characteristic input impedance of the resonator defined as with S_{e} the input crosssectional area, and are defined as follows:
The new coefficients (C_{n}, s_{n}) also verify the relation , since one switched from a 3real parameters description with (A_{n}, ω_{n}, ξ_{n}) to a 2complex parameters description with (C_{n}, s_{n}).
The fitted impedance Z_{fit}(ω) is plotted (orange dashed curve) in Figure 2.
Reinjecting equation (5) in equation (3) and applying an inverse Fourier transform (for details, see Appendix A) leads to the following expression in the time domain [13, 17]:
The mouthpiece pressure p can then be written as . Eventually, the system {(1) ∪ (2) ∪ (7)} is now written in the form , with:
the state vector, so that h = X_{1}, , and . Taking the real and imaginary part of the N equation (7) yields 2N real equations, so that the nonlinear function F is defined as:
with .
2.2.3 Regularisation of the volume airflow
For the continuation method implemented in AUTO to work correctly, F has to be a vector function. In this respect, noting that sgn(x) = x/x for x ≠ 0 and , equation (2) is regularised based on the regularisation used in Section 2 of [18]:
or equivalently,
in terms of the components of the state vector X, with η the regularisation parameter, which is fixed to 10^{−6} in the following. p_{0} is defined in a similar way as the closure pressure for woodwind instruments [14]: , except that ω_{0} is chosen close to the first resonance frequency of the resonator. Indeed, the choice of generally considered for woodwind instruments is not suitable in the case of brass instruments where is no longer constant. In practice, the choice is not ω_{0} = ω_{res,1} (where are the resonance angular frequencies of the resonator), but rather ω_{0} = ω_{res,4}/4, because the fourth resonance of instruments having the same tube length (a trombone and a euphonium for instance) appears to be quite constant, contrary to the first resonance frequency which can vary up to 8 semitones between a trombone and euphonium.
In the following, the system equation (9) with u(X) given by equation (11) are processed numerically in a dimensionless form, which is detailed in Appendix B. Typical bifurcation diagrams will be shown in Section 4.
3 Lips’ parameters and linear stability analysis
The results of the linear stability analysis detailed in Section 3.1 suggest that the lips’ parameter can be tuned so as to obtain a given desired oscillatory regime. The values of the other lips’ parameters are difficult to measure and thus hardly found in the literature. In a first subsection, we consider the same parameters as in [2], which are detailed in Table 1. In order to be more consistent with experimental results, we introduce, in a second subsection; a frequencydependent lipopening at rest, based on experimental data from [4].
3.1 Constant lipopening height at rest
In the case of constant lips’ parameters, the results of the linear stability analysis of the system are represented in Figure 4 by the blue Ushaped patterns. As described in Section 2.2, these results give the values of the mouth pressure p_{m} at which the equilibrium solution destabilises (the threshold mouth pressures), as well as the related frequency at which the oscillating solution emerges (the threshold frequencies). As a matter of fact, the multiple Ushaped patterns on the top plot of Figure 4 reflect the fact that for a given configuration of the resonator (i.e. without moving the slide nor activating any trigger), the trombone player can play several notes called “natural notes” just by changing the lips’ resonance frequency (see for instance [1]) represented in Figure 3, from the B♭1 (58 Hz) to B♭4 (466 Hz) approximately. For each natural note or regime n to emerge, there is an optimal lips’ resonance frequency which corresponds to the minimum of a Ushaped pattern. This is associated with an optimal threshold mouth pressure and an optimal threshold frequency . The notation has been chosen to be as consistent as possible with [2]: in this article, the equilibrium optimal threshold value (i.e. the threshold inferred from the linear stability analysis) of a quantity q in the nth regime was written . Here, the threshold quantities are obtained either using linear stability analysis or through a bifurcation analysis, which will be defined and addressed in Section 4. Therefore, it has been chosen to add an extra superscript “eq” for “equilibrium”, referring to the linear stability analysis, or “per” for “periodic solutions”, referring to the analysis of bifurcation diagrams. Furthermore, it can be noticed that the patterns on the bottom plot of Figure 4 are always above the line whose equation is , which is a characteristic of the outwardstriking valve model: the instrument always plays a note slightly above the lips’ resonance frequency.
Figure 3 Natural notes playable by a trombone, i.e. without moving the slide nor activating any trigger. From left to right: B♭1 (pedal note), B♭2, F3, B♭3, D4, F4, A♭4 (naturally a bit flat compared to an equal tempered scale), B♭4. 
Figure 4 Results of the linear stability analysis with the constant, linear and inverse Hprofile described on Figure 5. Top and bottom plots represent respectively the threshold mouth pressure and the threshold frequency versus the lips’ resonance frequency. Blue: H = h_{0}, orange: , green: . Circles point out the minima of each Ushaped pattern on the top plot. Red points represent the experimental values of minimal threshold pressures given by [19]. The black dashed line on the bottom plot represents . 
3.2 Frequencydependent lipopening height at rest
In order to assess the validity of the results presented above, a comparison is drawn with the minimal threshold mouth pressures estimated experimentally in [19], for the exact same trombone as the one considered in this paper. In Figure 4, the comparison between the minima of the blue Ushaped patterns and the red dots (representing the experimental measures of the minimal threshold mouth pressures) highlights significant differences. In particular, the results of the linear stability analysis, on the one hand, and of the experiment, on the other hand, are quite different: compared to the experimental data, threshold pressures computed numerically are too low for the first regimes, and too high for the last regimes. To overcome this issue, additional information coming from experiments in [4] is considered. Indeed, the lips’ opening height at rest H is measured to decrease monotonously with . The experimental points are plotted in black on the top plot of Figure 5, together with two regression functions which aim at fitting these points. The inverse function (k is given in Table 1) appears here to better fit the experimental points than the linear regression . Therefore, the former Hprofile will be retained hereafter, even if the case H = h_{0} is still given alongside for sake of comparison until the end of Section 4.
Figure 5 Considered Hprofiles versus lips’ resonance frequency. Blue line corresponds to the constant value of H = h_{0} from [2], black points correspond to experimental data from [4], and the orange and green curves correspond to the regression functions chosen to fit the experimental points, respectively and . 
4 Bifurcation diagrams and minimal blowing pressure
In this section several typical bifurcation diagrams of a trombone are shown and discussed in detail. The additional information obtained with the continuation method on the threshold mouth pressure is also discussed, leading to new conclusions compared to the case of the linear stability analysis.
4.1 Case of the third regime
In this section the example of the third regime of a trombone is considered (F3, third Ushaped pattern on top plot of Fig. 4), as it shows a variety of different behaviours of the system.
The left plot of Figure 6 shows the bifurcation diagrams obtained by continuation for a constant Hprofile H = h_{0} and for three different but close values of . For each branch of periodic solution, there is a global minimum of the mouth pressure (identified by the vertical dotted lines) below which the regime does not exist anymore. This minimum can be lower than or equal to the linear threshold mouth pressure given by the linear stability analysis in Section 3. As for this latter method, an optimal threshold mouth pressure ^{1} which can be lower than or equal to is introduced, as well as the associated optimal threshold frequency and the associated optimal lips’ resonance frequency . For instance, in the case = 175.09 Hz (green curve) the observed bifurcation is a direct Hopf [10] and the minimal blowing pressure is the same as the one given by the linear stability analysis. On the contrary, in the case (blue curve) and (orange curve), a small stable portion at the very beginning of the branch is noticed, thus rigorously classifying the bifurcation as a direct Hopf. However, since this stable portion only exists for a narrow range of p_{m} it would be very difficult to observe it experimentally, hence we choose to neglect it. In the following, this configuration will be referred to as “almostinverse Hopf” instead. In the case of this almostinverse Hopf bifurcation, the solution emerging from the equilibrium is unstable (if we neglect the small stable portion at the very beginning) and stabilises further on the branch at a value of threshold mouth pressure lower than at the equilibrium. For , the actual minimal blowing pressure is not only below the linear threshold pressure but also below the minimal blowing pressure given by the case .
Figure 6 Top and bottom plots represent respectively the maximum amplitude of the periodic oscillation branches and their oscillation frequency versus the blowing pressure. Left: case H = h_{0}; right: case . Blue curves: case = 165.84 Hz (left) or 167.75 Hz (right); orange curves: case = 164.00 Hz (left) or 165.17 Hz (right); green curves: case = 175.09 Hz, which corresponds to an arbitrary value at which the Hopf bifurcations are direct. The line thickness indicates whether the branch portion is stable (thick line) or unstable (thin line). The vertical dashed lines locate the value of the minimal threshold mouth pressure of each curve. 
However, it is worth noting that the linear stability analysis provides a rather accurate first guess of the note which is the easiest to play. Indeed, the frequency interval between the easiest note according to the linear stability analysis (corresponding to the minimal blowing pressure of the blue curve on Fig. 6) and according to the bifurcation diagrams (corresponding to the minimal blowing pressure of the orange curve on Fig. 6) is less than a semitone.
Eventually, it is reassuring to note that even if the variable Hprofile has a significant impact on the threshold pressure of the system, it has nonetheless qualitatively no influence on its threshold frequencies as shown on bottom plot of Figure 4.
For the interested reader, an animation showing the evolution of the bifurcation diagrams as a function of the lips’ resonance frequency is provided online.^{2} It especially shows how the bifurcation can go from almostinverse Hopf to direct Hopf through the merging of the two stable portions of the branch into a single one.
4.2 Minimal blowing pressures
The investigation performed in Section 4.1 for two values of is repeated for a wide range of lips’ resonance frequency. For each value of , the global minimal blowing pressure is extracted from the bifurcation diagrams, and these threshold values are compared with the results of linear stability analysis. The results of both methods are represented in Figure 7, which clearly highlights the added value of the continuation approach. In particular, the fact the first regimes (B♭1, B♭2 and F3) emerge through inverse Hopf bifurcations results in an actual threshold pressure which is lower than the one predicted by the linear stability analysis. Nevertheless, it is worth noting that the system only exhibit direct Hopf bifurcations for regimes ≥4. This means that the linear stability analysis is sufficient to characterise the threshold mouth pressures for these regimes. Nevertheless, it is worth noting that similar investigations performed for a trumpet (which is very similar to a trombone considering its bore, except that its pitch is an octave higher) in [13] highlighted an inverse Hopf bifurcation for the fourth regime. Also, the behaviour observed here for the trombone has also been observed in the case of the euphonium or the saxhorn, except for the first regime which will be addressed in Section 5.
Figure 7 Top and bottom plots represent respectively the threshold pressures and threshold frequencies given by the linear stability analysis (blue) and the continuation method (orange) versus the lips’ frequency. Left: case H = h_{0}; right: case . The circles identify the value of corresponding to a local minimum of blowing pressure. Horizontal dotted lines on the bottom plots locate the values of the acoustical resonances of the resonator; the black dashed line on the bottom plots represents . 
As in Section 4.1, it is reassuring here to note that the variable Hprofile has qualitatively no influence on the system’s behaviour discussed above, in particular its threshold frequencies as shown on bottom plots of Figure 7.
Among all the regimes that can be played by the musician according to Figures 3 and 7, the first one is especially worth focusing on. Indeed, it exhibits a very different behaviour depending on whether the instrument has a predominantlycylindrical bore profile (case of the trombone or the trumpet for instance) or a predominantlyexpanding bore profile (case of the tubas or the flugelhorn); see for instance Figure 7.34 of [1]. In this respect, the last section of this paper focuses on the study of the first regime of the trombone (predominantlycylindrical bore profile) and the euphonium (predominantlyexpanding bore profile).
5 Bifurcation diagrams, case of the first regime
5.1 Case of the trombone
Similarly to the case of the third regime (see Sect. 4.1), the bifurcation diagrams of the trombone show that the first regime generally emerges through an inverse Hopf bifurcations (see Fig. 8). As already discussed in Section 4.2, this leads to an optimal threshold pressure smaller than . However, it is worth noting that, similarly to other regimes, only one note is accessible to the musician since only one branch of the bifurcation diagrams is stable. The corresponding note is called the “pedal note” (B♭1, lowest note playable in first position and roughly one octave below the second regime B♭2, see also Fig. 3).
Figure 8 Top and bottom plots represent respectively the maximum amplitude of the periodic solutions branches and their oscillation frequency versus the blowing pressure in the case of a trombone, for = 54.86 Hz. The line thickness indicates whether the branch portion is stable (thick line) or unstable (thin line). 
The optimal playing frequency of the pedal note predicted by the bifurcation diagrams is about a semitone lower than the optimal linear threshold frequency given by the linear stability analysis used in [2]. Therefore, the analysis of the bifurcation diagrams adds a significant information as regards the playing frequency, even though – as previously discussed in Section 4.2 – the linear stability analysis already provides a rather good estimation of the playing frequency.
5.2 Case of a bass brass instruments with predominantlyexpanding bore profile
For the regimes 2–8, the bifurcation diagrams of the euphonium (not shown here) are very similar to the trombone’s ones, in terms of the direct or inverse nature of the Hopf bifurcations. However, the euphonium exhibits a specific behaviour for the first regime, which is qualitatively different from the one observed for the trombone.
In the case of the euphonium, the left part of Figure 10 shows typical bifurcation diagrams, obtained for a value of lips’ resonance frequency in the range of the first regime. Contrary to the trombone, these diagrams display two stable regimes corresponding to two distinct periodic solution branches, thus resulting in two different playing regimes accessible to the musician. The note with the lowest oscillation frequency is called – as for the trombone – the “pedal note” (B♭1, lowest note playable without any valve pushed, see also Fig. 9), whereas the other one is referred to as the “ghost note”^{3} (see Fig. 9).
Figure 9 Natural notes playable by a euphonium, i.e. without depressing any valves. From left to right: B♭1 (pedal note, red), D2 (ghost note, green), B♭2, F3, B♭3, D4, F4, A♭4 (naturally a bit flat compared to an equal tempered scale), B♭4. 
On the one hand, the pedal note (red note on Fig. 9) is a natural note, which is as familiar as any other regime for both trombone and euphonium players. On the other hand, the ghost note (green note on Fig. 9) is a regime exclusively accessible to bass brass instruments with predominantlyexpanding bore profile as a whole, such as the euphonium. Even though the ghost note is not much used in a musical context by euphonium players, it can be played by any tuba player and was brought to light for the first time in [2] using linear stability analysis. In particular, it was the only note playable for values of close to the first acoustic resonance of the resonator according to this method. As a matter of fact, since the branch corresponding to the pedal note (see the left plots of Fig. 10) does not emerge from the equilibrium solution (whereas the branch of the ghost note does), this regime could not be found by linear stability analysis which gives trustworthy information only in the vicinity of an equilibrium. Furthermore, the ghost note can be linked to one of the “factitious notes” mentioned in [20], where it is referred to as a “phenomenon whereby a player can sound a note intermediate between notes low in the series of natural notes” and sounds a perfect fourth above the pedal note in the case of an E♭tuba.
Figure 10 Typical bifurcation diagrams (left) and timedomain simulations (right) in the case of the euphonium’s first regime, obtained for p_{m} = 1.8 kPa and (PedalNote) = 49 Hz. Left plots represent the maximum amplitude of the periodic oscillation branches versus the blowing pressure (top) and the frequency of the corresponding periodic solutions versus the blowing pressure (bottom). The line thickness indicates whether the branch portion is stable (thick line) or unstable (thin line). Right plots represent the time evolution of the acoustic pressure at two points of the bifurcation diagrams above, given by the intersection of the vertical dotted line and the branches. Top: timedomain simulation of the ghost note (blue point on the bifurcation diagrams); bottom: timedomain simulation of the pedal note (orange point on the bifurcation diagrams). The horizontal dotted line represents the value of p_{m}. 
To illustrate the uncommon behaviour of the bifurcation diagrams shown in Figure 10 (left) and (right) shows results of timedomain simulations. The first simulation (left) corresponds to the blue point on the ghost note branch of the bifurcation diagrams. Its instantaneous frequency is 69.5 Hz, which is coherent with the oscillation frequency displayed on the top right plot in Figure 10. The second timedomain simulation (right) corresponds to the orange point on the pedal note branch of the bifurcation diagrams. Its instantaneous frequency is 60.6 Hz, which is coherent with the oscillation frequency displayed on the top left bottom in Figure 10.
In the case of the pedal note, it is worth noting that over one period, p(t) becomes higher than p_{m} (horizontal dotted line), which is physically questionable as it was never observed experimentally. This phenomenon might be considered as a weakness of the model used, caused by the raw modelling of the contact between the lips in equation (2). To overcome this limitation, a contact force between the lips – currently unavailable in literature – could be added in this equation. However, it is worth mentioning that such a contact force has already been introduced in the case of reed instruments [21].
5.3 Assessment of the frequency interval between pedal note and ghost note
It is possible to qualitatively check the relevance, in terms of intonation, of these two regimes predicted by the bifurcation analysis. In particular, one can compare the frequency intervals between the “easiesttoplay” notes predicted by the bifurcation diagrams on the one hand, and by analysing recordings of a euphonium player on the other hand. In the case of a euphonium and considering the Ushaped patterns in Figure A.1, the “easiesttoplay” notes are those with the lowest threshold mouth pressure. The corresponding oscillation frequencies are = 60.6 Hz for the pedal note, = 75.6 Hz for the ghost note, and = 120.4 Hz for the B♭2 (in the bifurcation diagrams, we have actually ). Therefore, it has been found that the frequency interval between the pedal note and the ghost note inferred from the bifurcation analysis (75.6/60.6 ≈ 1.25 = 5/4) is about 2 tones (major third), whereas the recording of a euphonium player^{4} leads to a high major third between the pedal note and the ghost note, which is reasonably close. Furthermore, the frequency interval between the pedal note and the second regime obtained in the bifurcation diagrams is found to be close to 12 semitones (octave), as one would expect between a B♭1 and a B♭2.
Figure A.1 Top and bottom plots represent respectively the threshold pressures and threshold frequencies given by the linear stability analysis (blue) and the continuation method (orange) versus the lips’ resonance frequency in the case of a euphonium. The circles identify the value of corresponding to a local minimum of blowing pressure. In the case of the second regime, the easiest note playable is given by the minimum of the continuation method (orange curve), since its minimal threshold mouth pressure is lower than the one given by the linear stability analysis (blue curve). The first two regimes on the left of each plot represent respectively the pedal note and the ghost note. 
Eventually, it is worth highlighting that the ghost note does exist for other tuba types than the euphonium, such as the E♭tuba for instance (see Footnote 4).
6 Conclusion
Most results in this study highlight the usefulness of linear stability analysis to understand various nearthreshold behaviours of a complete nonlinear model of brass instrument applied to trombone and euphonium, as well as the ability of bifurcation diagrams to quickly give valuable information about their periodic solutions.
The ease of playing of an instrument is chosen to be assessed based on the minimal threshold pressures. A good estimate can be obtained using linear stability analysis. At the same time, the associated threshold oscillation frequency is inferred, which is directly linked to the pitch of the note. For the considered elementary model of brass instrument, this method provides information on all the periodic regimes accessible to a trombone player for a given position of the slide, and on almost all periodic regimes accessible for a given fingering of a euphonium.
However, linear stability analysis turns out to be unable to describe the pedal note of the euphonium. In this respect, the analysis of bifurcation diagrams are conveniently introduced. It allows to illustrate the existence of the pedal regime in the form of a branch of periodic solutions far from equilibrium (hence the failure of linear stability analysis to describe it) separated from the one describing the ghost note, which itself emerges from the destabilisation of the equilibrium solution. Moreover, compared to the linear stability analysis, bifurcation diagrams give a more precise prediction of the threshold mouth pressures in the event of an inverse Hopf bifurcation. However, it is worth highlighting that the minimal threshold mouth pressure of each periodic regimes for a given slide position or fingering in the case of a trombone or a euphonium – except the pedal note of the latter – given by the bifurcation diagrams is qualitatively the same as given by the linear stability analysis.
If the ghost note seems to exist on every model of tubas, the interval between the ghost note and the first regime appears to vary from a minor third to a perfect fourth depending on the tuba’s pitch. Future work include a more quantitative study on a wide range of tubas going from the B♭contrabass tuba to the B♭baritone horn (an octave higher), in order to determine whether or not the frequency interval between the pedal note and the ghost note remains constant regardless of the considered bass brass instrument with predominantlyexpanding bore profile.
Moreover, it is worth noting that the pedal note also exhibits different behaviours for the high brass instruments, depending on the bore profile of the instrument. For example, the pedal note appears to be more “slotted” (meaning that the player has the sensation of a welldefined pitch) on a flugelhorn (predominantlyexpanding bore profile) than on a trumpet (predominantlycylindrical bore profile).
One important limitation of the present study lies in the difficulty to estimate the lips’ parameters used in the model. Indeed, it would be reasonable to think that lips’ parameters would vary from playing a low note to playing a high note, not to mention from playing a contrabass tuba to playing an alto trombone for instance. Yet, these lips’ parameters being difficult to measure both on an artificial mouth or directly on the player, they were kept the same regardless the instrument or the note played. Even though the obtained results look reasonable, i.e. consistent with the musicians’™ experience, in vivo measurements of lips parameters during musical performance would be very valuable. Furthermore, even if the system dynamics does not seem to depend sensitively on the model’s parameters such as the lipopening height at rest H, it may be worth checking the impact of a finer modelling of various elements on the ease of playing and pitch: for instance, considering a mouth section proportional to H^{2} (instead of only H in the present study), taking a second resonance frequency of the lips into account, or adding a contact force between the lips when they are about to touch each other.
Appendix A
Inverse Fourier transform of the input impedance
This appendix is based on the demonstration found in Section 6.3 of [22].
Quantities in the Fourier domain are written using capital letters, whereas quantities in the time domain are written using small letters.
Starting from equation (5), we have P(ω) = Z_{fit}(ω)U(ω) by definition of the input impedance. Since Z_{fit} represents the impulse response of the resonator, the inverse Fourier transform of Z_{fit} has to be real. Thus, Z_{fit} is hermitiansymmetric, that is to say . The pressure in the mouthpiece in the time domain p can then be inferred by taking the inverse Fourier transform of P:
with ω′ = −ω. U being also hermitiansymmetric, we have:
where p_{n} is the complex component of the pressure in the mouthpiece p, whose dynamical behaviour is linked to C_{n} and s_{n}:
Appendix B
Dimensionless variables and equations
From a numerical point of view, it is reasonable to make equations dimensionless in order to prevent the variables from varying on a toolarge scale, especially when using the continuation method in AUTO. In all that follows, most dimensionless quantity are written with a tilde “~” unless otherwise stated.
The dimensionless variables of the system , the dimensionless time , the dimensionless mouth pressure γ, the dimensionless lips’ resonance frequency and the dimensionless modal parameters (κ_{n}, σ_{n}) are defined as follows:
where h_{0} is given in Table 1, is the quarter of the fourth resonance angular frequency of the resonator (chosen as so because the fourth resonance frequency appears to be quite similar between instruments of the same tube length, such as a trombone and a euphonium for instance), so that p_{0} can be defined in a similar way as the closure pressure for woodwind instruments [18]: . Given the new set of variables, the state vector is now:
so equation (9) is now written as:
with now written as:
and eventually the dimensionless lipopening height at rest, defined similarly to [14].
Conflict of interest
The authors declare that they have no conflicts of interest in relation to this article.
Acknowledgments
We wish to thank Emmanuel Brasseur and Julia Mourier for useful discussions and for their help in this work.
See Section 3.1 for details about this notation.
Animation showing the evolution of the bifurcation diagrams as a function of the lips’ resonance frequency increasing from 134.2 Hz to 182.1 Hz in the case of the constant Hprofile H = h_{0} at http://perso.univlemans.fr/~rmatte/bd_dyn_AC_basse_reg_reg3.avi or in the case of the inverse Hprofile at http://perso.univlemans.fr/~rmatte/bd_dyn_AC_basse_Hinv_reg_reg3.avi.
We chose here to use the term employed by [22] to qualify this note being between a minor third and a perfect fourth higher than the pedal note.
Recording of a euphonium player at http://perso.univlemans.fr/~rmatte/bbeuphonium_reg16.mp4 and of an E♭tuba player at http://perso.univlemans.fr/~rmatte/ebtuba_reg16.mp4 playing successively an ascending series and a descending series of the six first natural notes, going through the ghost note when descending.
References
 M. Campbell, J. Gilbert, A. Myers: The Science of Brass Instruments. SpringerVerlag, 2021. [Google Scholar]
 L. Velut, C. Vergez, J. Gilbert, M. Djahanbani: How well can linear stability analysis predict the behavior of an outward valve brass instrument model?. Acta Acustica United with Acustica 103, 1 (2017) 132–148. [Google Scholar]
 J. Cullen, J. Gilbert, M. Campbell: Brass instruments: Linear stability analysis and experiments with an artificial mouth. Acta Acustica United with Acustica 86, 4 (2000) 704–724. [Google Scholar]
 S.J. Elliott, J.M. Bowsher: Regeneration in brass wind instruments. Journal of Sound and Vibration 83, 2 (1982) 181–217. [Google Scholar]
 H. Berjamin, B. Lombard, C. Vergez, E. Cottanceau: Timedomain numerical modeling of brass instruments including nonlinear wave propagation, viscothermal losses, and lips vibration. Acta Acustica United with Acustica 103 (2017) 117–131. [Google Scholar]
 R. Msallam, S. Dequidt, R. Caussé, S. Tassart: Physical model of the trombone including nonlinear effects. Application to the sound synthesis of loud tones. Acta Acustica United with Acustica 86 (2000) 725–736. [Google Scholar]
 A. Hirschberg, J. Gilbert, R. Msallam, A.P.J. Wijnands: Shock waves in trombones. Journal of the Acoustical Society of America 99, 3 (1995) 1754–1758. [Google Scholar]
 M.J. Newton, M. Campbell, J. Gilbert: Mechanical response measurements of real and artificial brass players lips. The Journal of the Acoustical Society of America 123, 1 (2008) EL14–EL20. [Google Scholar]
 S. Bromage, M. Campbell, J. Gilbert: Open areas of vibrating lips in trombone playing. Acta Acustica United with Acustica 96, 4 (2010) 603–613. [Google Scholar]
 Y.A. Kuznetsov: Elements of Applied Bifurcation Theory. SpringerVerlag, 2004. [Google Scholar]
 B. Cochelin, C. Vergez: A high order purely frequencybased harmonic balance formulation for continuation of periodic solutions. Journal of Sound and Vibration 324, 1–2 (2009) 243–262. [Google Scholar]
 E. Doedel, A. Champneys, T. Fairgrieve, Y. Kuznetsov, B. Sandstede, X. Wang: AUTO 97: Continuation and bifurcation software for ordinary differential equations (with HomCont). Tech. report. [Google Scholar]
 V. Fréour, L. Guillot, H. Masuda, S. Usa, E. Tominaga, Y. Tohgi, C. Vergez, B. Cochelin: Numerical continuation of a physical model of brass instruments: Application to trumpet comparisons. The Journal of the Acoustical Society of America 148, 2 (2020) 748–758. [Google Scholar]
 J. Gilbert, S. Maugeais, C. Vergez: Minimal blowing pressure allowing periodic oscillations in a simplified reed musical instrument model: BouasseBenade prescription assessed through numerical continuation. Acta Acustica 4, 6 (2020) 27. [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, 4 (2013) 968–977. [Google Scholar]
 F. Ablitzer: Peakpicking identification technique for modal expansion of input impedance of brass instruments. Acta Acustica 5 (2021) 53. [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, 1 (2014) 126–138. [Google Scholar]
 T. Colinot, L. Guillot, C. Vergez, P. Guillemain, J.B. Doc, B. Cochelin: Influence of the “Ghost Reed” Simplification on the Bifurcation Diagram of a Saxophone Model. Acta Acustica United with Acustica 105, 6 (2019) 1291–1294. [Google Scholar]
 J. Gilbert, L. Leblanc, C. Vergez: L’analyse de stabilité linéaire pour évaluer la facilité d’émission des cuivres – étude comparative de trombones ténor et basse, in Proc. 14 e Congrès Français d’Acoustique, Le Havre, France. 2018. [Google Scholar]
 T. Herbert, A. Myers, J. Wallace (Eds.): The Cambridge Encyclopedia of Brass Instruments. Cambridge University Press, Cambridge, UK New York, NY Port Melbourne New Delhi Singapore, 2019. [Google Scholar]
 A. Munoz Arancon, B. Gazengel, J.P. Dalmont: In vivo and in vitro characterization of single cane reeds, in Proceedings of the Stockholm Music Acoustics Conference, Roberto Bresin & Anders Askenfelt, Editors, KTH Royal Institute of Technology, Stockholm. 2013. [Google Scholar]
 F. Silva: Émergence Des AutoOscillations Dans Un Instrument de Musique à Anche Simple. Theses, Université de Provence – AixMarseille I, 2009. [Google Scholar]
Cite this article as: Mattéoli R. Gilbert J. Vergez C. Dalmont J.P. Maugeais S, et al. 2021. Minimal blowing pressure allowing periodic oscillations in a model of bass brass instruments. Acta Acustica, 5, 57.
All Tables
All Figures
Figure 1 Schematic representation of the player’s mouth coupled to the mouthpiece of a brass instrument. 

In the text 
Figure 2 Modulus (top) and phase (bottom) of the input impedance of a bass trombone (A. Courtois Legend AC502) with respect to frequency. Blue curve: measured impedance; dashed orange curve: modal fit function with N = 15 modes; green curve: error between measured and fitted impedances, defined by Z_{fit} − Z_{exp} and arg(Z_{fit}) − arg(Z_{exp}) respectively as regards the module and phase. 

In the text 
Figure 3 Natural notes playable by a trombone, i.e. without moving the slide nor activating any trigger. From left to right: B♭1 (pedal note), B♭2, F3, B♭3, D4, F4, A♭4 (naturally a bit flat compared to an equal tempered scale), B♭4. 

In the text 
Figure 4 Results of the linear stability analysis with the constant, linear and inverse Hprofile described on Figure 5. Top and bottom plots represent respectively the threshold mouth pressure and the threshold frequency versus the lips’ resonance frequency. Blue: H = h_{0}, orange: , green: . Circles point out the minima of each Ushaped pattern on the top plot. Red points represent the experimental values of minimal threshold pressures given by [19]. The black dashed line on the bottom plot represents . 

In the text 
Figure 5 Considered Hprofiles versus lips’ resonance frequency. Blue line corresponds to the constant value of H = h_{0} from [2], black points correspond to experimental data from [4], and the orange and green curves correspond to the regression functions chosen to fit the experimental points, respectively and . 

In the text 
Figure 6 Top and bottom plots represent respectively the maximum amplitude of the periodic oscillation branches and their oscillation frequency versus the blowing pressure. Left: case H = h_{0}; right: case . Blue curves: case = 165.84 Hz (left) or 167.75 Hz (right); orange curves: case = 164.00 Hz (left) or 165.17 Hz (right); green curves: case = 175.09 Hz, which corresponds to an arbitrary value at which the Hopf bifurcations are direct. The line thickness indicates whether the branch portion is stable (thick line) or unstable (thin line). The vertical dashed lines locate the value of the minimal threshold mouth pressure of each curve. 

In the text 
Figure 7 Top and bottom plots represent respectively the threshold pressures and threshold frequencies given by the linear stability analysis (blue) and the continuation method (orange) versus the lips’ frequency. Left: case H = h_{0}; right: case . The circles identify the value of corresponding to a local minimum of blowing pressure. Horizontal dotted lines on the bottom plots locate the values of the acoustical resonances of the resonator; the black dashed line on the bottom plots represents . 

In the text 
Figure 8 Top and bottom plots represent respectively the maximum amplitude of the periodic solutions branches and their oscillation frequency versus the blowing pressure in the case of a trombone, for = 54.86 Hz. The line thickness indicates whether the branch portion is stable (thick line) or unstable (thin line). 

In the text 
Figure 9 Natural notes playable by a euphonium, i.e. without depressing any valves. From left to right: B♭1 (pedal note, red), D2 (ghost note, green), B♭2, F3, B♭3, D4, F4, A♭4 (naturally a bit flat compared to an equal tempered scale), B♭4. 

In the text 
Figure 10 Typical bifurcation diagrams (left) and timedomain simulations (right) in the case of the euphonium’s first regime, obtained for p_{m} = 1.8 kPa and (PedalNote) = 49 Hz. Left plots represent the maximum amplitude of the periodic oscillation branches versus the blowing pressure (top) and the frequency of the corresponding periodic solutions versus the blowing pressure (bottom). The line thickness indicates whether the branch portion is stable (thick line) or unstable (thin line). Right plots represent the time evolution of the acoustic pressure at two points of the bifurcation diagrams above, given by the intersection of the vertical dotted line and the branches. Top: timedomain simulation of the ghost note (blue point on the bifurcation diagrams); bottom: timedomain simulation of the pedal note (orange point on the bifurcation diagrams). The horizontal dotted line represents the value of p_{m}. 

In the text 
Figure A.1 Top and bottom plots represent respectively the threshold pressures and threshold frequencies given by the linear stability analysis (blue) and the continuation method (orange) versus the lips’ resonance frequency in the case of a euphonium. The circles identify the value of corresponding to a local minimum of blowing pressure. In the case of the second regime, the easiest note playable is given by the minimum of the continuation method (orange curve), since its minimal threshold mouth pressure is lower than the one given by the linear stability analysis (blue curve). The first two regimes on the left of each plot represent respectively the pedal note and the ghost note. 

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.