Issue 
Acta Acust.
Volume 7, 2023



Article Number  16  
Number of page(s)  9  
Section  Musical Acoustics  
DOI  https://doi.org/10.1051/aacus/2023007  
Published online  12 May 2023 
Scientific Article
Amplitudedependent modal coefficients accounting for localized nonlinear losses in a timedomain integration of woodwind model
^{1}
Buffet Crampon, 5 rue Maurice Berteaux, 78711 ManteslaVille, France
^{2}
Aix Marseille Univ, CNRS, Centrale Marseille, LMA, 13013 Marseille, France
^{*} Corresponding author: szwarcberg@lma.cnrsmrs.fr
Received:
3
February
2023
Accepted:
23
March
2023
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. [Journal of Sound and Vibration 528 (2022) 116892] to the study of forced systems. It is now implemented for selfoscillating systems. The employed method extends the definition of the input impedance to the nonlinear domain by adding a dependance on the RMS acoustic velocity at a geometric discontinuity. The poles and residues resulting from the modal decomposition are fitted as a function of this velocity. Thus, the pressureflow relation defined 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. [PhD thesis, Université du Maine (2004)] and Dalmont and Frappé [Journal of the Acoustical Society of America 122(2) (2007) 1173–1179] were carried out. Simulations show that the model reproduces these experimental results qualitatively and quantitatively.
Key words: Musical acoustics / Nonlinear losses / Reed instruments / Modal decomposition
© The Author(s), Published by EDP Sciences, 2023
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
Sound synthesis by modal decomposition is based on an input impedance measurement, which captures the passive acoustical response of a real instrument. This method thus has the advantage of capturing the acoustical subtleties that differentiate two clarinets by means of a single measurement. Other methods, such as delay lines [1], waveguides [2, 3] or spatial discretization [4, 5] would require more effort on the geometrical description of the resonator to render such acoustical subtleties.
Furthermore, sound synthesis by modal decomposition requires little RAM compared to such other methods, which can be beneficial for embedded devices. Indeed, in modal decomposition synthesis, the temporal integration scheme requires to keep only a few previous iterations (precisely two for [1]) to compute a new one. For the delay line synthesis, it is necessary to keep in memory all the temporal iterations during a round trip, i.e. during 2LF_{s}/c_{0} iterations for an academic cylindrical resonator without lateral holes. Although this is not a problem for most embedded processors, fine modeling of viscothermal losses in the waveguide formalism also requires additional computations involving fractional derivatives [6]. Finally, the ODE formalism is particularly well suited to the bifurcation analysis, which is not the case of the PDE [6] or of the DDE formalisms [7].
The input impedance is a measure of the linear frequency response of the resonator, for a lowamplitude excitation. However, in the case of woodwind instrument playing, the measured acoustic pressure and velocity in the resonator can be very high. For instance, [8] measure a pressure in the mouthpiece of 4 kPa, i.e. 163 dB for a monochromatic wave. According to [9] of Chap. 8.4.5, from an acoustic speed of roughly 1 m/s, jet separation phenomena appear. At the level of a geometrical discontinuity, such as the open end of a pipe [10] or a lateral orifice [11], vortices are observed. A part of the kinetic energy of the jet is absorbed by these vortices and is dissipated as heat by friction. These losses are accounted for through a nonlinear relationship derived from Bernoulli’s law, as demonstrated by the implementation in waveguide modeling or delay lines of Refs. [7, 12–14].
Following [11], the works from Atig et al. [12, 15] and Dalmont et al. [13, 16] model localized nonlinear losses through a resistive impedance Z_{t}. This impedance is in series with the radiation impedance Z_{R}, and depends on the amplitude of the acoustic velocity v_{0} at the location of the discontinuity:
(1)where Z_{c} = ρ_{0}c_{0}/S is the characteristic impedance of the medium for plane waves, and c_{d} ∈ [0; 3] is a parameter depending on the geometry of the termination. The larger the radius of curvature at the output, the smaller this coefficient is. The coefficient c_{d} is difficult to predict theoretically, and easier to determine experimentally [12]. As shown by Atig et al. [12], nonlinear losses have a significant influence on the playing range of a clarinet, hence models of sound production in woodwinds should include such effects.
Models of woodwind instruments taking into account localized nonlinear losses at the end of the resonator have never been developped in the framework of modal synthesis. This is the main contribution of this work.
First, a method proposed by Diab et al. [17] to account for nonlinear losses by modal decomposition is presented. It is then applied to a cylindrical tube with nonlinear losses located at the open end. From this resonator, a selfoscillating clarinetlike system is defined and simulated by time integration. These simulations are finally compared to the experimental results published by Atig et al. [12] and Dalmont and Frappé [13].
2 Method of modal decomposition accounting for localized nonlinear losses
The definition of a “nonlinear impedance” (i.e. adjusting the impedance so that the link between acoustic flow and pressure stays valid in nonlinear conditions) opens the way to consider localized nonlinear losses by modal decomposition. In Diab et al. [17], the evolution of the surface impedance of perforated plates is computed by temporal simulation for a broadband excitation, with respect to the RMS amplitude of the acoustic velocity at the level of the hole, noted v_{RMS}. From the decomposition of the impedance in the linear domain into a sum of N modes, two methods are studied for increasing values of v_{RMS}: the interpolation of the impedance, on one hand, and the regression of the poles and complex residues (s_{n}, C_{n}) on the other hand. This second method will be considered hereafter.
Nonlinear effects are taken into account by allowing the poles and residues to vary with respect to v_{RMS}. In the Laplace domain, the modal decomposition of the input impedance is written:
(2)where s is the Laplace variable and •* denotes the complex conjugation operation. The relationship between poles and residues and v_{RMS} is assumed to be a rational fraction [17]. In the present work, the approximation is limited to a polynomial of degree N_{p}:
(4)where the index •^{(0)} refers to the linear part of the modal coefficient. In Diab et al. [17], the regression coefficients of the poles and residues are then adjusted on a nonlinear surface impedance model Z_{NL} (v_{RMS}), given by Laly et al. [18]. It is worth noting that in the present paper, the poles and residues are fitted on the dimensioned quantity of the RMS velocity. Therefore, the unit of coefficients and are expressed in .
The last step consists in computing v_{RMS}, given by:
(5)where v_{0} (t) = u_{0}/S is the acoustic velocity at the discontinuity, of cross section S. This definition is rather impractical in the framework of a timedomain simulation. In order to obtain a state formalism, it is replaced by:
(6)which is equivalent. By considering the righthand side of equation (6) as an input for the ODE solver, can be computed, hence v_{RMS}.
3 Application to a cylindrical tube with nonlinear losses at the open end
The impedance regression technique presented by Diab et al. [17] is applied to the case of a closedopen cylinder of radius R and length L, which is a similar case study to Ref. [12]. Viscothermal losses are taken into account in the propagation:
(7)where η = 3 × 10^{−5} s^{1/2}. The boundary condition at x = L is given by the radiation impedance Z_{R}, which includes in series the nonlinear impedance Z_{t} defined by equation (1):
(8)and jk = s/c_{0}, Δl ≈ 0.6R ([9], Chap. 12.6.1.3), K_{NL} = 4c_{d}/(3π). In the following, v_{RMS}(L, t) will be denoted as v_{RMS} for reading comfort. The input impedance is defined by:
(9)introducing the dimensionless notation z_{•} = Z_{•}/Z_{c}. A linear input impedance is also defined, such that
(10)which is equivalent to z_{in} for v_{RMS} = 0.
Figure 1 illustrates the evolution of z_{in} (computed through equation (9) and equation (10), in which the Laplace variable is substituted by j2πf) as a function of the RMS velocity at the open end. It can be observed that taking into account nonlinear losses has a consequence on the amplitude, but not on the frequency of the resonance peaks. Figure 1b shows that this impact is particularly accentuated on the first resonance peak: its amplitude is 43% lower for v_{RMS} = 24 m/s than in the case without nonlinear losses. The difference is reduced to 30% for the second peak, and to 16% for the sixth peak. Near the antiresonances, the negative relative discrepancies are high, since the impedance values are all close to zero. The absolute deviation remains very low. From the lack of frequency variation of the resonance peaks with v_{RMS}, we could anticipate that nonlinear losses at the end of the pipe would rather impact dynamics than intonation.
Figure 1 Input impedance of a cylinder z_{in} (computed through Eq. (9) and Eq. (10)) of length L = 64 cm and radius R = 8 mm for different values of the acoustic RMS velocity at the open end. The value of c_{d} has been set to 13/9, which corresponds to a radius of curvature of 0.3 mm, according to [29]. From top to bottom, from left to right: modulus of the input impedance; relative gap between nonlinear and linear definition of z_{in}; detail view on the first peak; detail view on the fourth peak. 
3.1 Modal decomposition of the nonlinear input impedance
In order to take into account nonlinear losses in temporal simulation, the modal coefficients of the input impedance must be determined. In the same way that the nonlinear input impedance depends on the acoustic velocity at the open end, the coefficients C_{n} and s_{n} derived from the modal decomposition of z_{in} also depend on v_{RMS}. To formulate the relationship between the modal coefficients and v_{RMS}, a minor adaptation of the method presented in Ref. [7] of Chap. 5.5.3 is used. The definition of z_{R} used in [9] is substituted by the expression given by equation (8), taking into account nonlinear losses at the end of the tube.
3.1.1 Expression of the poles s_{n}
Equation (9) can also be written as:
The poles are the solutions of cosh[ΓL + h] = 0, i.e.:
The solutions of equation (12) give s_{n}, for different input values of v_{RMS}. They are plotted on Figure 2a. The consideration of nonlinear losses at the end of a cylindrical pipe has almost no influence on , which is the resonance angular frequency of peak n. However, a larger increases . This shift remains almost the same regardless of the index of the pole. However, in terms of relative deviation, the ratio is decreasing as n increases. This is signaled by the amplitude of peaks of higher index appearing less affected by nonlinear losses, as shown in Figure 1b.
Figure 2 Graphical representation of the N = 8 first poles and residues. The values of v_{RMS} are linearly chosen between 0 m/s and 24 m/s, according to [29]. As in Figure 1, c_{d} = 13/9. (a) Real and imaginary parts of poles s_{n} (v_{RMS}), calculated by solving equation (12). Detailed view on s_{1}. (b) Real and imaginary parts of residues C_{n} (v_{RMS)}, calculated with equation (14). Detailed views on C_{1} and C_{4}. 
3.1.2 Expression of the residues C_{n}
It remains to calculate residue C_{n} (v_{RMS}). In the vicinity of a pole s_{n}, the denominator D of z_{in} can be written through a firstorder series expansion of cosh:
(13)introducing notation •′ = ∂•/∂s. According to Chaigne and Kergomard ([7], Chap. 5.5.3), after application of the residues theorem, the modal decomposition of the input impedance for the nth peak can be written as
The evolution of coefficients C_{n} with respect to v_{RMS} is plotted on Figure 2b. Real and imaginary parts of C_{n} slightly decrease when v_{RMS} increases. Following Diab et al. [17], it remains to fit the modal coefficients with respect to v_{RMS}, as in the example of equation (3).
3.1.3 Polynomial fitting of the modal coefficients
Figure 3a shows the mean relative fitting error for polynomials of different degrees, for the same data as in Figure 2. For linear regression, the mean relative error stays around 10^{−7} for each mode n, both for C_{n} and s_{n}. Further simulations reveal that choosing an excessive values of c_{d} tends to increase this error. For a tube with sharp edges at the open end, Atig et al. [12] estimated a maximum value of c_{d} = 2.8. For c_{d} = 5, for instance, the mean relative error of linear regression is around 10^{−5}. Although this error has increased, it remains very low.
Figure 3 Error due to the regression of modal coefficients C_{n} and s_{n} (a), and its consequence on the discrepancy between modal decomposition and definition (9) on z_{in} (b). (a) Mean relative error ε between fitted modal coefficient and its actual value X = {s_{n}, C_{n}}. Regression is based on N_{v} = 25 values of v_{RMS}. (b) L_{2}norm of the error over z_{in} for every frequencies, between definition given by equation (9) and modal decomposition approximation defined by equation (2). This error is plotted with respect to the degree of the fitting polynomial of C_{n} and s_{n}. 
Moreover, Figure 3b shows that choosing an polynomial regression degree larger than 1 has no consequence on the error over z_{in} between modal decomposition with fitted coefficients and the definition given by equation (9). This error is therefore inherent to the modal decomposition approximation (Eq. (2)) of the input impedance (Eq. (9)).
According to these previous results, in the rest of this article, the modal coefficients will both be fitted by a polynomial of order 1.
4 Application to the timeintegration simulation of a clarinetlike system
The cylindrical tube considered in Section 3 is now treated as the resonator of a clarinetlike instrument. The physical model governing the selfsustained oscillations of the instrument is presented in Section 4.1. It is then simulated by time integration, imposing a linear increase of the blowing pressure (crescendo) followed by a linear decrease (diminuendo). Details and results from the simulations are given in Section 4.2.
4.1 Physical model for the selfsustained oscillations
The selfsustained system consists of three main parts: reed dynamics, reed channel and resonator. This last subsystem is enriched compared to the classical modal decomposition formalism [19–21] by the addition of nonlinear losses at the end of the pipe.
4.1.1 Reed dynamics
The dynamic behavior of the reed is modeled for its dimensionless displacement x(t) by the following equation:
(15)where ω_{r} = 2π × 2200 rad/s is the reed resonance angular frequency [22] and q_{r} = 0.4 is the reed damping [23]. γ(t) = P_{m}(t)/p_{M} is the dimensionless blowing parameter and p(t) = p^{dim} (t)/p_{M} is the dimensionless pressure at the input of the resonator. The beating pressure p_{M} is the value of the blowing pressure P_{m} for which the reed closes the reed channel, in quasistatic regime. is the contact force function. In the following, a model of “ghost reed” [24] is considered, i.e. F_{c} = 0.
4.1.2 Reed channel
The characteristic of the flow u(t) through the reed channel is described by the following equation, according to [25]:
(16)where λ = 5.5 × 10^{−3}/c_{0} is the reed flow parameter [26], ζ is the embouchure parameter, and γ(t) is the dimensionless blowing pressure. The operator [•]^{+} = (• + •)/2 refers to the positive part function. For the simulations, the absolute value functions are regularized to smooth the irregularities of the quasistatic flow characteristic, without altering its behavior, according to [27]. Thus, , with η = 0.001.
4.1.3 Resonator
The resonator is described in the frequency domain by its input impedance z_{in}, under its modal decomposition form given by equation (2). The modal poles and residues s_{n} and C_{n} are chosen to be linearly dependent on v_{RMS}, i.e.:
The modal decomposition of z_{in} allows to write the relation between the dimensionless pressure p and flow u in the temporal domain:
The values of s_{n} and C_{n} must be updated at each time step by computing v_{RMS} (L, t). The method for calculating the mean velocity at nonlinearity is detailed in the following section.
4.1.4 Computation of v_{RMS} at each time step
To compute v_{RMS} (L, t), the pressure at the termination L must be calculated. First, a linear problem is considered. The value of s_{n} is therefore chosen for v_{RMS} = 0 m/s. The modal components of the pressure are related to the total pressure at L through the following equation:
(18)and ξ ∈ [0, L] is the distance along the longitudinal axis of the resonator. For reading comfort, the expression p(t) will be reserved to denote p(0, t).
The dimensionless acoustic velocity at the open end is then calculated from its definition which is related to the pressure field through the dimensionless Euler equation:
The dimensionless acoustic velocity v(L, t) is thus obtained by numerical integration of the following equation:
The RMS velocity is finally obtained by double time integration of equation (20), according to equation (6). Since the equations of the system characterizing the selfsustained oscillations are dimensionless, it is necessary to resize v_{RMS}, because the modal coefficients and presented in equation (3) have been regressed from the dimensioned velocity. To do so, is multiplied by p_{M}/(ρ_{0}c_{0}). In the following simulations, a value of p_{M} = 8.5 kPa is set, according to experimental results from [12].
Regarding the calculation of v(L, t), it is not numerically guaranteed that the mean value of remains zero. Consequently, the integral may diverge. To avoid these divergence problems when integrating and v^{2}(L, t), a shortmemory term τ is added within each integral. Thus, v(L, t) is now defined by the following convolution product:
(21)which is written, in the Laplace () domain, as:
By going back to the time domain, equation (21) becomes:
The same method is applied for the computation of v_{RMS}:
4.1.5 Complete system of equations
The complete selfsustained system is governed by the equations defined in Sections 4.1.1, 4.1.2, 4.1.3 and 4.1.4. This article shows the resolution of this problem using an ode solver. The Cauchy problem writes, according to equation (28), as:
(28)where u(t) is computed by using equation (16), C_{k} (v_{RMS}) are computed with equation (14), and s_{k}(v_{RMS}) are computed with equation (12). The initial conditions will be set to 0 in the following simulations.
4.2 Simulation including nonlinear losses at the open end
The model including nonlinear losses at the open end as represented by the system of equation (28) is now simulated by time integration, using the solver ode45 from Matlab. The input data of the problem are based on experimental results from [12].
4.2.1 Simulation parameters
The parameters related to the resonator, the reed and the reed channel are presented in Table 1. Six values of c_{d} have been taken from [12], which are between c_{d} = 0 and c_{d} = 2.8. The value of ζ = 0.28 was calculated using data corresponding to a “loose embouchure” configuration. For the calculation of , the shortmemory term was set to the period of the first impedance peak, i.e. . The evolution of γ is first linearly ascending, from γ = 0 to γ = 3 in 8 s, then linearly descending, from γ = 3 to γ = 0 in 8 s. The total simulation time is therefore 16 s. Simulations were performed for N = 4 modes. A convergence study showed that the dynamic behavior of the system did not change for a higher number of modes.
Main constants used for the simulations.
4.2.2 Results
The bifurcation diagram of the input pressure p^{(dim)} = p · p_{M} over the blowing pressure P_{m} = γ · p_{M} is represented on Figure 4, with the parameters detailed in Section 4.2.1.
Figure 4 Bifurcation diagram of the L_{2} norm of p^{(dim)} with respect to P_{m}, for different values of c_{d}, in an ascending blowing pressure configuration (a) and a descending one (b). 
During a crescendo (Fig. 4a), the oscillation threshold is located near 5.7 kPa (γ_{th} ≈ 0.67). This threshold remains the same, whether nonlinear losses are taken into account (c_{d} ≠ 0) or not (c_{d} = 0). However, nonlinear losses have a significant influence on the extinction threshold , i.e. the blowing pressure from which the reed stops oscillating and is pressed completely against the mouthpiece. This influence is shown in Table 2. When c_{d} is increased (i.e. nonlinear losses increase), γ_{extup} diminishes. Similarly to the experiments conducted by Atig et al. [12], nonlinear losses have an important influence on the dynamic playing range of the musician. Figure 4 shows that when nonlinear losses are low, the range of stable oscillation amplitude that can be obtained is larger.
In the diminuendo phase, the inverse threshold is the same (γ_{inv} ≈ 0.89) for each geometry at the open end of the pipe. Around this threshold, the amplitude of the input pressure decreases slightly as the losses increase. This behavior is also observed in the experimental curves of Ref. [12] of Figures 12a, 12c and 12e.
This comparaison with the experimental results of Atig et al. is limited to a qualitative study. Parameters such as the reed resonance angular frequency ω_{r} and reed damping q_{r} have a strong influence on the dynamics of the system, as demonstrated by Silva et al. [28] for the oscillation threshold. The complete recalibration of the model on experimental results is beyond the scope of this article.
4.3 Comparison with experimental results from Dalmont and Frappé (2007)
An attempt to quantitatively validate the model of nonlinear losses is carried out. To do so, the values of the saturation threshold γ_{sat} measured by Dalmont and Frappé [13] are employed (see Fig. 5). This threshold corresponds to the blowing pressure γ for which p is maximum during a crescendo. The authors measured this threshold for different reed openings, which are transcribed here in different values of ζ. The dimensioning pressure values p_{M} used here correspond to the beating pressure estimated by the authors for a diminuendo. Furthermore, the measurements are performed on a cylindrical tube of dimensions L = 50 cm and R = 8 mm. The geometry at the opening has been estimated by the authors at c_{d} = 2.8.
Figure 5 Experimental results from [13] for the saturation threshold (·), compared with three methods: modal decomposition without nonlinear losses (+); modal decomposition including nonlinear losses (×); analytical solution provided by Raman model including nonlinear losses (*). 
The experimental data are compared to simulations using the same parameters L and R as in the experiment. Simulations are performed for two cases: the first one does not take into account nonlinear losses (c_{d} = 0), the other one takes them into account (c_{d} = 2.8). The simulations are performed for ascending ramps of blowing pressure varying from to in 10 s. The other parameters are the same as those given in Table 1.
The model presented in this paper is also compared to the Raman model including nonlinear losses proposed by Dalmont and Frappé [13]. In this model, viscothermal losses are simplified by a frequencyindependent coefficient α which replaces . The simplicity of this model allows the authors to obtain an analytical expression of the saturation threshold as a function of the parameter c_{d}.
In the present work, the value of α was first adjusted at , where is the angular frequency of the impedance peak supporting the oscillation. However, this value produced excessively high estimates of compared to experimental results, especially for high values of . To better match the experimental data, was chosen. It should be noted that the analytical expression of exhibits a linear dependence in , for low values of α. Thus, is extremely sensitive to the value of α in Raman’s model.
Figure 5 illustrates the evolution of the saturation threshold as a function of in the experimental case, through modal decomposition simulation, and through the analytical expression based on the Raman model. It appears first that the saturation threshold is overestimated when nonlinear losses are not taken into account (+). In comparison, the two models taking into account nonlinear losses produce results much closer to the experiment. This shows the importance of taking this phenomenon into account in wind instrument simulations.
Furthermore, for the analytical solution given by the Raman model (*), the evolution of follows a straight line whose slope depends on . Although the results are close to the experimental ones, the high sensitivity of to α as well as the linear evolution of reflect the excessive simplicity of the Raman model to describe the dynamics of a wind instrument.
Finally, the model presented in this paper (×) has a good agreement with the experimental data, except for very low values of the mouthpiece parameter (). For these low values of , the model overestimates compared to the experimental results. This overestimation may be related to the transient behavior at extinction caused by the temporal evolution of the control parameter . A characterization of the system by continuation could give results independent of the evolution rate of .
In conclusion, the correspondence between the results from the model presented in this article and the experimental data from [13] highlights its ability to describe the dynamic behavior of a simplified wind instrument at saturation.
5 Conclusion
This article develops the design of a sound synthesis model of a reed instrument by modal decomposition of the input impedance, taking into account viscothermal losses as well as nonlinear losses at the end of the resonator. The input impedance now depends on the RMS acoustic velocity at a geometric discontinuity (here, the open termination). Poles and residues resulting from the modal decomposition are fitted as a function of this velocity. In a physical model of wind instrument, the pressureflow relation defined by the resonator is then completed by new equations which account for this dependence with the velocity at the end of the pipe.
To evaluate the ability of the model to reproduce a real phenomenon, comparisons with the experimental results of [12] and [13] have been made. In the first case, simulations show a qualitatively similar behavior regarding the evolution of the extinction threshold depending on the geometry at the open end. In the second case, the model gives a good correspondence with the experimental results, in particular compared to a model without nonlinear losses.
This formalism could be promising in the sound synthesis of wind instruments in real time, which employs the modal decomposition formalism in order to keep a minimum of memory. The inclusion of nonlinear losses in a complete clarinet model could more accurately translate dynamic phenomena essential to the experience of the musician, as suggested by Dalmont and Frappé [13]: “any realistic model of the clarinet should include nonlinear losses in the side holes”.
Conflict of interest
The authors declare no conflict of interest.
Funding
This study has been supported by the French ANR LabCom LIAMFI (ANR16LCV200701).
Acknowledgments
The authors thank S. Maugeais and J.P. Dalmont for their valuable comments.
References
 P. Guillemain, J. Kergomard, T. Voinier: Realtime synthesis of clarinetlike instruments using digital impedance models. Journal of the Acoustical Society of America 118, 1 (2005) 483–494. [Google Scholar]
 G.P. Scavone: An acoustic analysis of singlereed woodwind instruments with an emphasis on design and performance issues and digital waveguide modeling techniques, PhD thesis, Stanford University, Stanford, CA, 1997. [Google Scholar]
 R. Mignot: Réalisation en guides d’ondes numériques stables d’un modèle acoustique réaliste pour la simulation en tempsréel d’instruments à vent. PhD thesis, Télécom ParisTech, 2009. [Google Scholar]
 S. Bilbao: Numerical sound synthesis: finite difference schemes and simulation in musical acoustics, John Wiley & Sons, 2009. [Google Scholar]
 A. Ernoult, J. Chabassier, S. Rodriguez, A. Humeau: Full waveform inversion for bore reconstruction of woodwindlike instruments. Acta Acustica 5 (2021) 47. [CrossRef] [EDP Sciences] [Google Scholar]
 S. Bilbao, J. Chick: Finite difference time domain simulation for the brass instrument bore. Journal of the Acoustical Society of America 134, 5 (2013) 3860–3871. [Google Scholar]
 P. Guillemain, J. Terroir: Digital synthesis models of clarinetlike instruments including nonlinear losses in the resonator, in the 9th International Conference on Digital Audio Effects, (2006). 83. [Google Scholar]
 B. Bergeot, A. Almeida, B. Gazengel, C. Vergez, D. Ferrand: Response of an artificially blown clarinet to different blowing pressure profiles. Journal of the Acoustical Society of America 135, 1 (2014) 479–490. [Google Scholar]
 A. Chaigne, J. Kergomard: Acoustics of musical instruments, Springer. 2016. [CrossRef] [Google Scholar]
 J.H.M. Disselhorst, L. Van Wijngaarden: Flow in the exit of open pipes during acoustic resonance. Journal of Fluid Mechanics 99, 2 (1980) 293–319. [Google Scholar]
 U. Ingård, S. Labate: Acoustic circulation effects and the nonlinear impedance of orifices. Journal of the Acoustical Society of America 22, 2 (1950) 211–218. [Google Scholar]
 M. Atig, J.P. Dalmont, J. Gilbert: Saturation mechanism in clarinetlike instruments, the effect of the localised nonlinear losses. Applied Acoustics 65, 12 (2004) 1133–1154. [CrossRef] [Google Scholar]
 J.P. Dalmont, C. Frappé: Oscillation and extinction thresholds of the clarinet: comparison of analytical results and experiments. Journal of the Acoustical Society of America 122, 2 (2007) 1173–1179. [Google Scholar]
 P.A. Taillard: Theoretical and experimental study of the role of the reed in clarinet playing. PhD thesis, Université du Maine, 2018. [Google Scholar]
 M. Atig, J.P. Dalmont, J. Gilbert: Termination impedance of openended cylindrical tubes at high sound pressure level. Comptes Rendus Mécanique 332, 4 (2004) 299–304. [CrossRef] [Google Scholar]
 J.P. Dalmont, C. Nederveen, V. Dubos, S. Ollivier, V. Méserette, E. Sligte: Experimental determination of the equivalent circuit of an open side hole: linear and non linear behaviour. Acta Acustica united with Acustica 88 (2002) 567–575. [Google Scholar]
 D. Diab, D. Dragna, E. Salze, M.A. Galland: Nonlinear broadband timedomain admittance boundary condition for duct acoustics. application to perforated plate liners. Journal of Sound and Vibration 528 (2022) 116892. [Google Scholar]
 Z. Laly, N. Atalla, S.A. Meslioui: Acoustical modeling of microperforated panel at high sound pressure levels using equivalent fluid approach. Journal of Sound and Vibration 427 (2018) 134–158. [Google Scholar]
 V. Dubos, J. Kergomard, A. Khettabi, J.P. Dalmont, D.H. Keefe, C.J. Nederveen: Theory of sound propagation in a duct with a branched tube using modal decomposition. Acta Acustica united with Acustica 85, 2 (1999) 153–169. [Google Scholar]
 F. Silva, V. Debut, J. Kergomard, C. Vergez, A. Deblevid, P. Guillemain: Simulation of single reed instruments oscillations based on modal decomposition of bore and reed dynamics, in Proceedings of the International Congress of Acoustics, 2007. [Google Scholar]
 P.A. Taillard, F. Silva, Ph Guillemain, J. Kergomard: Modal analysis of the input impedance of wind instruments. application to the sound synthesis of a clarinet. Applied Acoustics 141 (2018) 271–280. [CrossRef] [Google Scholar]
 V. Chatziioannou, M. Van Walstijn: Estimation of clarinet reed parameters by inverse modelling. Acta Acustica united with Acustica 98, 4 (2012) 629–639. [CrossRef] [Google Scholar]
 F. Silva, Ph Guillemain, J. Kergomard, B. Mallaroni, A.N. Norris: Approximation formulae for the acoustic radiation impedance of a cylindrical pipe. Journal of Sound and Vibration 322, 1–2 (2009) 255–263. [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. [CrossRef] [Google Scholar]
 T.A. Wilson, G.S. Beavers: Operating modes of the clarinet. Journal of the Acoustical Society of America 56, 2 (1974) 653–658. [Google Scholar]
 J.P. Dalmont, J. Gilbert, S. Ollivier: Nonlinear characteristics of singlereed instruments: Quasistatic volume flow and reed opening measurements. Journal of the Acoustical Society of America 114, 4 (2003) 2253–2262. [Google Scholar]
 T. Colinot: Numerical simulation of woodwind dynamics: investigating nonlinear sound production behavior in saxophonelike instruments. PhD thesis, Laboratoire de Mécanique et d'Acoustique [Marseille], November, 2020. [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, 5 (2008) 3284–3295. [Google Scholar]
 M. Atig: Nonlinéarité acoustique localisée à l’extrémité ouverte d’un tube. Mesure, modélisation et application aux instruments à vent. PhD thesis, Université du Maine, 2004. [Google Scholar]
Cite this article as: Szwarcberg N. Colinot T. Vergez C. & Jousserand M. 2023. Amplitudedependent modal coefficients accounting for localized nonlinear losses in a timedomain integration of woodwind model. Acta Acustica, 7, 16.
All Tables
All Figures
Figure 1 Input impedance of a cylinder z_{in} (computed through Eq. (9) and Eq. (10)) of length L = 64 cm and radius R = 8 mm for different values of the acoustic RMS velocity at the open end. The value of c_{d} has been set to 13/9, which corresponds to a radius of curvature of 0.3 mm, according to [29]. From top to bottom, from left to right: modulus of the input impedance; relative gap between nonlinear and linear definition of z_{in}; detail view on the first peak; detail view on the fourth peak. 

In the text 
Figure 2 Graphical representation of the N = 8 first poles and residues. The values of v_{RMS} are linearly chosen between 0 m/s and 24 m/s, according to [29]. As in Figure 1, c_{d} = 13/9. (a) Real and imaginary parts of poles s_{n} (v_{RMS}), calculated by solving equation (12). Detailed view on s_{1}. (b) Real and imaginary parts of residues C_{n} (v_{RMS)}, calculated with equation (14). Detailed views on C_{1} and C_{4}. 

In the text 
Figure 3 Error due to the regression of modal coefficients C_{n} and s_{n} (a), and its consequence on the discrepancy between modal decomposition and definition (9) on z_{in} (b). (a) Mean relative error ε between fitted modal coefficient and its actual value X = {s_{n}, C_{n}}. Regression is based on N_{v} = 25 values of v_{RMS}. (b) L_{2}norm of the error over z_{in} for every frequencies, between definition given by equation (9) and modal decomposition approximation defined by equation (2). This error is plotted with respect to the degree of the fitting polynomial of C_{n} and s_{n}. 

In the text 
Figure 4 Bifurcation diagram of the L_{2} norm of p^{(dim)} with respect to P_{m}, for different values of c_{d}, in an ascending blowing pressure configuration (a) and a descending one (b). 

In the text 
Figure 5 Experimental results from [13] for the saturation threshold (·), compared with three methods: modal decomposition without nonlinear losses (+); modal decomposition including nonlinear losses (×); analytical solution provided by Raman model including nonlinear losses (*). 

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.