Search

Localized nonlinear losses are taken into account in a simple Raman clarinet model.The complete system is expressed as an iterated map, enabling to study the stability of the different playing regimes. A parametric study is carried out with respect to three major parameters: blowing pressure, embouchure and nonlinear losses coefficient.The model exhibits the well-known effect of reducing the extinction threshold when nonlinear losses increase.However, the stability analysis also shows that increasing nonlinear losses increases the oscillation threshold.

The Laboratory of Mechanics and Acoustics in Marseilles (France) was created in 1941, under the name of Centre de Recherches Scientifiques, Industrielles et Maritimes (CRSIM). But it was actually issued from the French Naval Research Center created in Toulon by the French Navy to work on submarine detection during World War I. LMA is therefore the result of a long and quite amazing story with several moves and even more name changes. It benefited from all these events and is today established in a new campus with large facilities specially designed for its latest research activities. This article presents the story in some details, summarize the evolution of the research domains through all these years and finally gives a description of the LMA today.

A cantilever beam under axial flow, confined or not, is known to develop self-sustained oscillations at sufficiently large flow velocities. In recent decades, the analysis of this archetypal system has been mostly pursued under linearized conditions, to calculate the critical boundaries separating stable from unstable behavior. However, nonlinear analysis of the self-sustained oscillations ensuing flutter instabilities are considerably rarer. Here we present a simplified one-dimensional nonlinear model describing a cantilever beam subjected to confined axial flow, for generic axial profiles of the fluid channels. In particular, we explore how the shape of the confinement walls affects the dynamics of the system. To simplify the problem, we consider symmetric channels with plane walls in either converging or diverging configurations. The beam is modeled in a modal framework, while bulk-flow equations, including singular head-loss terms, are used to model the flow-structure coupling forces. The dynamics of the system are first analyzed through linear stability analysis to assess the stabilizing/destabilizing effects of the channel walls configuration. Subsequently, we develop a systematic nonlinear analysis based on the continuation of periodic solutions. The harmonic balance method is used in conjunction with the asymptotic numerical method to calculate branches of periodic solutions. The continuation-based methods are used to investigate bifurcations with respect to both the reduced flow velocity and the channel slope parameter. From the results presented, we illustrate how continuationbased approaches and bifurcation analysis provide an efficient tool to analyze the nonlinear behavior of flow-induced vibration problems, particularly when reduced/simplified models are available., Comment: 22nd International Conference of Numerical Analysis and Applied Mathematics, Sep 2024, Heraklion, Greece

This study investigates the role of localized nonlinear losses in the register hole on the production of second-register notes. First, an experiment is conducted to study the ability of a register hole to produce second register. A cylindrical tube is drilled with holes of increasing diameter. Five are at the same level as the register hole of a B-flat clarinet, and five are at the same level as the thumb hole. Participant clarinetists are then asked to play with constant control parameters. At the beginning of each measurement, all holes are closed. The operator then opens randomly one of the ten holes.The resulting register is noted. The experiment is replicated numerically by time integration of two different models. The first is the state-of-the-art model based on the modal decomposition of the input impedance of the resonator. The second accounts for localized nonlinear losses in the register hole, through the model from Dalmont and Nederveen (2002). These losses are handled through a variable modal coefficients method. For the first model, simulations never produce second register, for any of the open holes. For the second, the proportion of second-register production is close to the experiment for upstream holes, but remains at zero for downstream holes.

This article develops the design of a sound synthesis model of a woodwind instrument by modal decomposition of the input impedance, taking into account viscothermal losses as well as localized nonlinear losses at the end of the resonator. This formalism has already been applied by Diab et al. (2022) to the study of forced systems. It is now implemented for self-oscillating systems. The employed method extends the denition of the input impedance to the nonlinear domain by adding a dependance on the RMS acoustic velocity at a geometric discontinuity. The poles and residuals resulting from the modal decomposition are interpolated as a function of this velocity. Thus, the pressure-ow relation dened by the resonator is completed by new equations which account for the dependence with the velocity at the end of the tube. To assess the ability of the model to reproduce a real phenomenon, comparisons with the experimental results of Atig et al. (2004) and Dalmont et al. (2007) were carried out. Simulations show that the model reproduces these experimental results qualitatively and quantitatively.

The ghost note is a natural note which can be played exclusively on bass brass instruments with a predominantly-expanding bore profile such as tubas, euphoniums or saxhorns. It stands between the pedal note-the lowest natural note playable, or first regime-and the instrument's second regime. However, if the interval between the pedal note and the second regime remains close to an octave regardless of the instrument, the interval between the pedal note and the ghost note vary from a minor third to a perfect fourth. References about this note are very scarce, and it is not commonly known among tuba players.This study shows that an elementary brass model describing the player coupled to the instrument is capable of bringing both the ghost and the pedal note to light. Here, we adopt a dynamical systems point of view and perform a bifurcation analysis using a software of numerical continuation. The numerical results provided in terms of frequency intervals between pedal note and ghost note are compared with frequency intervals experimentally inferred from recordings of seven different types of tuba, each of them being played by two professional tuba players., Comment: arXiv admin note: text overlap with arXiv:2112.08751

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 [Campbell et al., 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 in-depth 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.

A time-domain numerical modeling of brass instruments is proposed. On one hand, outgoing and incoming waves in the resonator are described by the Menguy-Gilbert model, which incorporates three key issues: nonlinear wave propagation, viscothermal losses, and a variable section. The non-linear propagation is simulated by a TVD scheme well-suited to non-smooth waves. The fractional derivatives induced by the viscothermal losses are replaced by a set of local-in-time memory variables. A splitting strategy is followed to couple optimally these dedicated methods. On the other hand, the exciter is described by a one-mass model for the lips. The Newmark method is used to integrate the nonlinear ordinary differential equation so-obtained. At each time step, a coupling is performed between the pressure in the tube and the displacement of the lips. Finally, an extensive set of validation tests is successfully completed. In particular, self-sustained oscillations of the lips are simulated by taking into account the nonlinear wave propagation in the tube. Simulations clearly indicate that the nonlinear wave propagation has a major influence on the timbre of the sound, as expected. Moreover, simulations also highlight an influence on playing frequencies, time envelopes and on the playability of the low frequencies in the case of a variable lips tension.

Simple models of clarinet instruments based on iterated maps have been used in the past to successfully estimate the threshold of oscillation of this instrument as a function of a constant blowing pressure. However, when the blowing pressure gradually increases through time, the oscillations appear at a much higher value than what is predicted in the static case. This is known as bifurcation delay, a phenomenon studied in [1] for a clarinet model. In numerical simulations the bifurcation delay showed a strong sensitivity to numerical precision., Comment: Stockholm Music Acoustics Conference, Stockholm : Su\`ede (2013). arXiv admin note: substantial text overlap with arXiv:1407.3547

Simple models of clarinet instruments based on iterated maps have been used in the past to successfully estimate the threshold of oscillation of this instrument as a function of a constant blowing pressure. However, when the blowing pressure gradually increases through time, the oscillations appear at a much higher value, called dynamic oscillation threshold, than what is predicted in the static case. This is known as bifurcation delay, a phenomenon studied in [1,2] for a clarinet model. In particular the dynamic oscillation threshold is predicted analytically when the blowing pressure is linearly increased. However, the mouth pressure cannot grow indefinitely. During a note attack, after an increasing phase, the musician stabilizes the mouth pressure. In the present work, the analytical prediction of the dynamic oscillation threshold is extended to a situations in which the mouth pressure approaches a steady state pressure according to an exponential time profile. The predictions still show a good agreement with simulation of the simple clarinet-model. This situation is compared in terms of dynamic oscillation bifurcation., Comment: International Symposium on Musical Acoustics, Le Mans : France (2014)

Since they correspond to a jump from a given note to another one, the mouth pressure thresholds leading to regime changes are particularly important quantities in flute-like instruments. In this paper, a comparison of such thresholds between an artificial mouth, an experienced flutist and a non player is provided. It highlights the ability of the experienced player to considerabily shift regime change thresholds, and thus to enlarge its control in terms of nuances and spectrum. Based on recent works on other wind instruments and on the theory of dynamic bifurcations, the hypothe- sis is tested experimentally and numerically that the dynamics of the blowing pressure influences regime change thresholds. The results highlight the strong influence of this parameter on thresholds, suggesting its wide use by experienced musicians. Starting from these observations and from an analysis of a physical model of flute-like instruments, involving numerical continuation methods and Floquet stability analysis, a phenomenological modelling of regime change is proposed and validated. It allows to predict the regime change thresholds in the dynamic case, in which time variations of the blowing pressure are taken into account.

Flautas de chinos - prehispanic chilean flutes played during ritual celebrations in central Chile - are known to produce very particular beating sounds, the so-called sonido rajado. Some previous works have focused on the spectral analysis of these sounds, and on the input impedance of the complex resonator. However, the beating sounds origin remains to be investigated. Throughout this paper, a comparison is provided between the characteristics of both the sound produced by flautas de chinos and a synthesis sound obtained through time-domain simulation of the jet-drive model for flute-like instruments. Jet-drive model appears to be able to produce quasiperiodic sounds similar to sonido rajado. Finally, the analysis of the system dynamics through numerical continuation methods allows to explore the production mechanism of these quasiperiodic regimes., Comment: Stockholm Music Acoustics Conference, Stockholm : Sweden (2013)

This paper proposes a methodology to map the various acoustic regimes of wind instruments. The maps can be generated in a multi-dimensional space consisting of design, control parameters, and initial conditions. The bound- aries of the maps are obtained explicitly in terms of the parameters using a support vector machine (SVM) classifier as well as a dedicated adaptive sam- pling scheme. The approach is demonstrated on a simplified clarinet model for which several maps are generated based on different criteria. Examples of computation of the probability of occurrence of a specific acoustic regime are also provided. In addition, the approach is demonstrated on a design optimization example for optimal intonation.

Using an artificial mouth with an accurate pressure control, the onset of the pressure oscillations inside the mouthpiece of a simplified clarinet is studied experimentally. Two time profiles are used for the blowing pressure: in a first set of experiments the pressure is increased at constant rates, then decreased at the same rate. In a second set of experiments the pressure rises at a constant rate and is then kept constant for an arbitrary period of time. In both cases the experiments are repeated for different increase rates. Numerical simulations using a simplified clarinet model blown with a constantly increasing mouth pressure are compared to the oscillating pressure obtained inside the mouthpiece. Both show that the beginning of the oscillations appears at a higher pressure values than the theoretical static threshold pressure, a manifestation of bifurcation delay. Experiments performed using an interrupted increase in mouth pressure show that the beginning of the oscillation occurs close to the stop in the increase of the pressure. Experimental results also highlight that the speed of the onset transient of the sound is roughly the same, independently of the duration of the increase phase of the blowing pressure., Comment: 14 pages

This paper presents an analysis of the effects of noise and precision on a simplified model of the clarinet driven by a variable control parameter. When the control parameter is varied the clarinet model undergoes a dynamic bifurcation. A consequence of this is the phenomenon of bifurcation delay: the bifurcation point is shifted from the static oscillation threshold to an higher value called dynamic oscillation threshold. In a previous work [8], the dynamic oscillation threshold is obtained analytically. In the present article, the sensitivity of the dynamic threshold on precision is analyzed as a stochastic variable introduced in the model. A new theoretical expression is given for the dynamic thresholds in presence of the stochastic variable, providing a fair prediction of the thresholds found in finite-precision simulations. These dynamic thresholds are found to depend on the increase rate and are independent on the initial value of the parameter, both in simulations and in theory., Comment: 14 pages

In this paper, we extend the method proposed by Cochelin and Vergez [A high order purely frequency-based harmonic balance formulation for continuation of periodic solutions, Journal of Sound and Vibration, 324 (2009) 243-262] to the case of non-polynomial nonlinearities. This extension allows for the computation of branches of periodic solutions of a broader class of nonlinear dynamical systems. The principle remains to transform the original ODE system into an extended polynomial quadratic system for an easy application of the harmonic balance method (HBM). The transformation of non-polynomial terms is based on the differentiation of state variables with respect to the time variable, shifting the nonlinear non-polynomial nonlinearity to a time-independent initial condition equation, not concerned with the HBM. The continuation of the resulting algebraic system is here performed by the asymptotic numerical method (high order Taylor series representation of the solution branch) using a further differentiation of the non-polynomial algebraic equation with respect to the path parameter. A one dof vibro-impact system is used to illustrate how an exponential nonlinearity is handled, showing that the method works at very high order, 1000 in that case. Various kinds of nonlinear functions are also treated, and finally the nonlinear free pendulum is addressed, showing that very accurate periodic solutions can be computed with the proposed method.

The tone hole geometry of a clarinet is optimized numerically. The instrument is modeled as a network of one dimensional transmission line elements. For each (non-fork) fingering, we first calculate the resonance frequencies of the input impedance peaks, and compare them with the frequencies of a mathematically even chromatic scale (equal temperament). A least square algorithm is then used to minimize the differences and to derive the geometry of the instrument. Various situations are studied, with and without dedicated register hole and/or enlargement of the bore. With a dedicated register hole, the differences can remain less than 10 musical cents throughout the whole usual range of a clarinet. The positions, diameters and lengths of the chimneys vary regularly over the whole length of the instrument, in contrast with usual clarinets. Nevertheless, we recover one usual feature of instruments, namely that gradually larger tone holes occur when the distance to the reed increases. A fully chromatic prototype instrument has been built to check these calculations, and tested experimentally with an artificial blowing machine, providing good agreement with the numerical predictions.

Self-sustained musical instruments (bowed string, woodwind and brass instruments) can be modeled by nonlinear lumped dynamical systems. Among these instruments, flutes and flue organ pipes present the particularity to be modeled as a delay dynamical system. In this paper, such a system, a toy model of flute-like instruments, is studied using numerical continuation. Equilibrium and periodic solutions are explored with respect to the blowing pressure, with focus on amplitude and frequency evolutions along the different solution branches, as well as "jumps" between periodic solution branches. The influence of a second model parameter (namely the inharmonicity) on the behaviour of the system is addressed. It is shown that harmonicity plays a key role in the presence of hysteresis or quasi-periodic regime. Throughout the paper, experimental results on a real instrument are presented to illustrate various phenomena, and allow some qualitative comparisons with numerical results.

Reed instruments are modeled as self-sustained oscillators driven by the pressure inside the mouth of the musician. A set of nonlinear equations connects the control parameters (mouth pressure, lip force) to the system output, hereby considered as the mouthpiece pressure. Clarinets can then be studied as dynamical systems, their steady behavior being dictated uniquely by the values of the control parameters. Considering the resonator as a lossless straight cylinder is a dramatic yet common simplification that allows for simulations using nonlinear iterative maps. In this paper, we investigate analytically the effect of a time-varying blowing pressure on the behavior of this simplified clarinet model. When the control parameter varies, results from the so-called dynamic bifurcation theory are required to properly analyze the system. This study highlights the phenomenon of bifurcation delay and defines a new quantity, the dynamic oscillation threshold. A theoretical estimation of the dynamic oscillation threshold is proposed and compared with numerical simulations., Comment: 14 pages

This paper focuses on the oscillation threshold of single reed instruments. Several characteristics such as blowing pressure at threshold, regime selection, and playing frequency are known to change radically when taking into account the reed dynamics and the flow induced by the reed motion. Previous works have shown interesting tendencies, using analytical expressions with simplified models. In the present study, a more elaborated physical model is considered. The influence of several parameters, depending on the reed properties, the design of the instrument or the control operated by the player, are studied. Previous results on the influence of the reed resonance frequency are confirmed. New results concerning the simultaneous influence of two model parameters on oscillation threshold, regime selection and playing frequency are presented and discussed. The authors use a numerical continuation approach. Numerical continuation consists in following a given solution of a set of equations when a parameter varies. Considering the instrument as a dynamical system, the oscillation threshold problem is formulated as a path following of Hopf bifurcations, generalizing the usual approach of the characteristic equation, as used in previous works. The proposed numerical approach proves to be useful for the study of musical instruments. It is complementary to analytical analysis and direct time-domain or frequency-domain simulations since it allows to derive information that is hardly reachable through simulation, without the approximations needed for analytical approach.

The concept of nonlinear modes is applied in order to analyze the behavior of a model of woodwind reed instruments. Using a modal expansion of the impedance of the instrument, and by projecting the equation for the acoustic pressure on the normal modes of the air column, a system of second order ordinary differential equations is obtained. The equations are coupled through the nonlinear relation describing the volume flow of air through the reed channel in response to the pressure difference across the reed. The system is treated using an amplitude-phase formulation for nonlinear modes, where the frequency and damping functions, as well as the invariant manifolds in the phase space, are unknowns to be determined. The formulation gives, without explicit integration of the underlying ordinary differential equation, access to the transient, the limit cycle, its period and stability. The process is illustrated for a model reduced to three normal modes of the air column.

Schlieren flow visualizations of transverse oscillations of jets submitted to an acoustic perturbation are analyzed in this paper. The aim is to estimate the shape and the position of the median line of the jet. Two methods for image processing are proposed, based on complementary approaches : inter image and intra image analysis. Synthesized images of an oscillating jet are used to validate each method and compare their performances in the case of noisy pictures. Illustrations are then shown on real laminar and turbulent jets. The results obtained using both methods are very close, showing their reliability. Applications investigated in this paper are focused on the estimation of the convection velocity of perturbations along the jet, and the influence of the Reynolds number and of the channel geometry upstream the jet formation.

Sound emergence in clarinetlike instruments is investigated in terms of instability of the static regime. Various models of reed-bore coupling are considered, from the pioneering work of Wilson and Beavers ["Operating modes of the clarinet", J. Acoust. Soc. Am. 56, 653--658 (1974)] to more recent modeling including viscothermal bore losses and vena contracta at the reed inlet. The pressure threshold above which these models may oscillate as well as the frequency of oscillation at threshold are calculated. In addition to Wilson and Beavers' previous conclusions concerning the role of the reed damping in the selection of the register the instrument will play on, the influence of the reed motion induced flow is also emphasized, particularly its effect on playing frequencies, contributing to reduce discrepancies between Wilson and Beavers' experimental results and theory, despite discrepancies still remain concerning the pressure threshold. Finally, analytical approximations of the oscillating solution based on Fourier series expansion are obtained in the vicinity of the threshold of oscillation. This allows to emphasize the conditions which determine the nature of the bifurcation (direct or inverse) through which the note may emerge, with therefore important consequences on the musical playing performances.

Combinig the harmonic balance method (HBM) and a continuation method is a well-known technique to follow the periodic solutions of dynamical systems when a control parameter is varied. However, since deriving the algebraic system containing the Fourier coefficients can be a highly cumbersome procedure, the classical HBM is often limited to polynomial (quadratic and cubic) nonlinearities and/or a few harmonics. Several variations on the classical HBM, such as the incremental HBM or the alternating frequency/time domain HBM, have been presented in the literature to overcome this shortcoming. Here, we present an alternative approach that can be applied to a very large class of dynamical systems (autonomous or forced) with smooth equations. The main idea is to systematically recast the dynamical system in quadratic polynomial form before applying the HBM. Once the equations have been rendered quadratic, it becomes obvious to derive the algebraic system and solve it by the so-called ANM continuation technique. Several classical examples are presented to illustrate the use of this numerical approach., Comment: PACS numbers: 02.30.Mv, 02.30.Nw, 02.30.Px, 02.60.-x, 02.70.-c

This paper considers a "striking outwards" reed coupled to a resonator. This expression, due to Helmholtz, is not discussed here : it corresponds to the most common model of a lip-type valve, when the valve is assumed to be a one degree of freedom oscillator. The presented work is an extension of the works done by Wilson and Beavers (1974), Tarnopolsky (2000). The range of the playing frequencies is investigated. The first results are analytical : when no losses are present in the resonator, it is proven that the ratio between the threshold frequency and the reed resonance frequency is found to be necessarily within the interval between unity and the square root of 3. This is a musical sixth. Actually the interval is largely smaller, and this is in accordance with e.g. the results by Cullen et al.. The smallest blowing pressure is found to be directly related to the quality factor of the reed. Numerical results confirm these statements, and are discussed in comparison with previous ones by Cullen et al (2000).

This paper investigates the sound production in a system made of a bore coupled with a reed valve. Extending previous work (Debut, 2004), the input impedance of the bore is projected on the modes of the air column. The acoustic pressure is therefore calculated as the sum of modal components. The airrrflow blown into the bore is modulated by reed motion, assuming the reed to be a single degree of freedom oscillator. Calculation of self-sustained oscillations controlled by time-varying mouth pressure and player's embouchure parameter is performed using ODE solvers. Results emphasize the par ticipation of the whole set of components in the mode locking process. Another impor tant feature is the mutual innnfluence of reed and bore resonance during growing blowing pressure transients, oscillation threshold being altered by the reed natural frequency and the reed damping. Steady-state oscillations are also investigated and compared with results given by harmonic balance method and by digital sound synthesis.

This article proposes a characterisation of the double-reed in quasi-static regimes. The non-linear relation between the pressure drop $\Delta p$ in the double-reed and the volume flow crossing it $q$ is measured for slow variations of these variables. The volume flow is determined from the pressure drop in a diaphragm replacing the instrument's bore. Measurements are compared to other experimental results on reed instrument exciters and to physical models, revealing that clarinet, oboe and bassoon quasi-static behavior relies on similar working principles. Differences in the experimental results are interpreted in terms of pressure recovery due to the conical diffuser role of the downstream part of double-reed mouthpieces (the staple).

The harmonic balance method (HBM) was originally developed for finding periodic solutions of electronical and mechanical systems under a periodic force, but has later been adapted to self-sustained musical instruments. Unlike time-domain methods, this frequency-domain method does not capture transients and so is not adapted for sound synthesis. However, its independence of time makes it very useful for studying every periodic solution of the model, whether stable or unstable without care of initial conditions. A computer program for solving general problems involving nonlinearly coupled exciter and resonator, Harmbal, has been developed based on the HBM. The method as well as convergence improvements and continuations facilities are thorougly presented and discussed in the present paper. Application of the method is demonstrated on various problems related to a common model of the clarinet: a reed modelled as a simple spring with and without mass and damping, a nonlinear coupling and a cubic simplification of it, and a cylindrical bore with or without dissipation and dispersion as well as a bore formed as a stepped cone.