Arxiv:1708.06571V1 [Nlin.AO] 22 Aug 2017 the Generating Unit, Which Is the Pressure Reservoir Be- Coupling: Neath the Pipe at a Basically Constant Rate
Total Page:16
File Type:pdf, Size:1020Kb
SAW16, dated September 28, 2018 Synchronization of organ pipes 1, 2 1 Jakub Sawicki, ∗ Markus Abel, and Eckehard Schöll 1Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstraße 36, 10623 Berlin, Germany 2Department of Physics and Astronomy, Potsdam University, Karl-Liebknecht-Straße 24, 14476 Potsdam Germany We investigate synchronization of coupled organ pipes. Synchronization and reflection in the organ lead to undesired weakening of the sound in special cases. Recent experiments have shown that sound interaction is highly complex and nonlinear, however, we show that two delay-coupled Van-der-Pol oscillators appear to be a good model for the occurring dynamical phenomena. Here the coupling is realized as distance-dependent, or time-delayed, equivalently. Analytically, we investigate the synchronization frequency and bifurcation scenarios which occur at the boundaries of the Arnold tongues. We successfully compare our results to experimental data. PACS numbers: 05.45.Xt, 05.45.-a Introduction - The physics of organ pipes is an in- we investigate the bifurcation scenarios in the context of terdisciplinary topic where many fields of science meet. two delay-coupled Van-der-Pol oscillators as a represen- It is highly interesting as it includes elements of non- tation of the system of two coupled organ pipes, such linear dynamical system theory [4, 6, 11], aeroacoustic as in the experimental setup of Bergweiler et al. [5]. modeling [13] and synchronization theory [17]. The fo- In extension of previous work, we study the dependence cus of these different research areas is the “queen of in- of Arnold tongues under variation of the time delay τ struments” which captivates through the grandeur of her and the coupling strength κ, to explore how undesired sight and majesty of her sound. Here, we investigate synchronization or chaotic behavior can be avoided. We an interesting nonlinear effect: organ pipes close to each compare our results to experimental measurements of the other synchronize. Recent studies have been of exper- synchronization under variation of the pipe distance and imental nature as well as theoretical [1, 2, 7, 8]. For find a qualitative coincidence of the nonmonotonic mod- musical purposes, synchronization of sound might be de- ulation of the shape of the Arnold tongue which is in sired or not: it might stabilize the pitch of special organ contrast to the linear boundaries of the Arnold tongues pipes as a favorable effect, whereas sound weakening, as for systems with undelayed coupling. observed in the prospect of an organ, is highly undesired In Section II we introduce two delay-coupled Van der [9, 18, 20]. This weakening occurs as an amplitude mini- Pol oscillators as a simple model of coupled organ pipes. mum due to destructive interaction between pipes during In section II we apply two analytical methods to get a the actuating of the swell box, where the pipes stand close better understanding of the synchronization phenomena. to each other. The central part of the paper is Sec. III, where we present A qualitative understanding of the nonlinear mecha- the analytical results. In Sec. IV we compare these results nisms is obtained following the arguments of [2, 7]: a with acoustic experiments. We conclude with Sec. V. single organ pipe can be described as a self-sustained os- cillator, where the oscillating unit consists of the jet, or “air sheet” which exits at the pipe mouth. The resonator I. A MODEL OF COUPLED ORGAN PIPES is of course the pipe body; there, sound waves emitted at the labium (i.e., the sharp edge in the upper part of To obtain a deeper insight into the synchronization the pipe’s opening) travel up and down and can trigger a phenomena of two coupled organ pipes we model the regular oscillation of the air sheet. Energy is supplied by pipes by Van der Pol oscillators with delayed cross- arXiv:1708.06571v1 [nlin.AO] 22 Aug 2017 the generating unit, which is the pressure reservoir be- coupling: neath the pipe at a basically constant rate. Experimen- 2 tal and numerical investigations by Abel et al. [2] yield x¨i + !i xi µ x_ i f_(xi) + κxj(t τ) = 0; (1) the conclusion that an organ pipe can be approximated − − − h i satisfactorily by a Van der Pol oscillator. Whereas the where i; j = 1; 2. These equations represent a harmonic work of Fischer [7] focuses on the nonlinearities in sound oscillator with an intrinsic angular frequency !i, supple- generation and their effect on the synchronization prop- mented with linear and nonlinear damping of strength erties, in this paper we investigate the effect of the finite µ > 0. The nonlinear damping can be described by the distance of two coupled pipes which in turn is reflected nonlinear function by a delay in the coupling function. More specifically, γ f(x ) = x3; (2) i 3 i where γ is the anisochronicity parameter and f_(xi) = ∗ corresponding author: [email protected] 2 γxi x_i. The coupling strength in Eq. (1) is κ, and the 2 II. ANALYTIC APPROACHES width of sync-region A. Method of averaging The method of averaging (quasiharmonic reduction) describes weakly nonlinear oscillations in terms of slowly varying amplitude and phase. For µ = 0 the uncoupled system reduces to the harmonic 2 oscillator x¨i + !i xi = 0 with solution ͙ xi = Ri sin(!it + φi); (4) with constant amplitude R and phase φ . For 0 < µ 1 i i ω-sync we look for a solution in the form Eq. (4) but assume that the amplitude R 0 and the phase φ are time- i ≥ i dependent functions: xi = Ri(t) sin(!it + φi(t)); (5) x_ i = Ri(t)!i cos(!it + φi(t)): where terms involving the slowly varying functions R_i, Φ_i are neglected. Without loss of generality, we choose !2 = FIG. 1: Angular frequency Ω of oscillator x1 (dark 1. For small µ we use the method of averaging, assuming green circles), and oscillator x2 (light yellow circles) that the product µτ is small, and Taylor expand R (t τ) i − versus the detuning ∆ of the oscillators. Full (empty) and φi(t τ) in the following way: circles correspond to symmetric (non-symmetric) initial − conditions. Parameters: ! = 1, µ = 0:1, γ = 1, κ = 0:4, τ 2 2 R (t τ) = R (t) τR_ (t) + R¨ (t) + :::; (6) τ = 0:1π. i − i − i 2 i We introduce the phase difference (t) = φ1(t) φ2(t). ~ 2t − Defining a new time scale t = µ we find the equations coupling delay is τ. Since for synchronization the fre- which describe the system (1) on a slow time scale: quency difference of the two oscillators is important, we ~ 2 introduce the detuning parameter ∆ R by R_ (t~) = R (t~) 1 γR1=2(t) 2 1=2 1=2 − 4 κR2=1(t~) sin( (t~) + τ); (7) 2 2 ∓ !1 = !2 + µ∆: (3) ~ _ ~ R1(t) ~ (t) = ∆ + κ ~ cos( (t) τ) − R2(t) − ~h R2(t) ~ In Fig. 1 we display frequency locking as obtained by nu- ~ cos( (t) + τ) : (8) − R1(t) merical simulation of Eq. (1) with symmetric initial con- i ditions . The observed angular frequencies Ω are plotted For simplicity, in the following we omit the tilda . The ∼ versus the detuning ∆ of two Van der Pol oscillators. A method of averaging together with truncation of the pronounced synchronization region and a sharp transi- Taylor expansion in τ reduces the infinite-dimensional tion to synchronization is observed. Within the synchro- problem to a finite-dimensional problem by assuming nization region, only for small ∆ the in-phase synchro- that the product µτ is small. This key step enables us j j nized solution is observed (lower branch), while for larger to handle the original delay differential equation as a ∆ the anti-phase synchronized solution (upper branch) system of ordinary differential equations [19, 22]. We j j is obtained even for symmetric initial conditions (full cir- now have two dynamical equations (7) for the amplitudes cles). Note that for small ∆ the anti-phase synchronized R1 and R2, and one equation (8) for the phase difference j j solution is also observed for any non-symmetric initial (t) which is also called slow phase. The latter equation conditions (empty circles). Our purpose is to analyze is a generalized Adler equation [3] and contains the main the synchronization frequency, the width of the synchro- features of synchronization. nization region, the phase difference in the synchronized state, its stability, and the bifurcation scenarios which Generalized Adler equation - The equilibria of the occur at the boundaries. Adler equation correspond to the locking of phase and 3 3 frequency, since the difference between the phases is con- B ×10-3 stant. To investigate the stability and bifurcation sce- E nario of such fixed points we take a closer look at the 6 generalized Adler equation(8), written in general form: width of 4 sync-region _(t) = ∆ + κq( (t)) (9) N(A)- G(iωs) − 2 where the averaged forcing term q( (t)) is the 2π- NONLINEAR LINEAR periodic function -2 A q( (t)) = R1(t) cos [ (t) τ] R2(t) cos [ (t) + τ] : R2(t) − − R1(t) -4 (10) -6 FIG. 3: Block diagram of a nonlinear feedback control The generalized Adler equation (8) is necessary for the loop consisting of a linear transfer function G(i!s) and calculation of the Arnold tongue (see below), which is one 0 1 2 3 4 5 6 anonlineardescribingfunctionN(A). A stands for the of the main characteristics of synchronization in nonlin- C feedback signal and E for the reference input, which is ear systems. For further analysis it is useful to eliminate ×10-2 added at the summation point (circle).