## Physical Wave Propagation Modeling for Real-Time Synthesis of Natural Sounds (2002)

### Cached

### Download Links

Citations: | 10 - 3 self |

### BibTeX

@TECHREPORT{Essl02physicalwave,

author = {Georg Essl},

title = {Physical Wave Propagation Modeling for Real-Time Synthesis of Natural Sounds},

institution = {},

year = {2002}

}

### OpenURL

### Abstract

This thesis proposes banded waveguide synthesis as an approach to real-time sound synthesis based on the underlying physics. So far three main approaches have been widely used: digital waveguide synthesis, modal synthesis and finite element methods. Digital waveguide synthesis is efficient and realistic and captures the complete dynamics of the underlying physics but is restricted to instruments that are well-described by the one-dimensional string equation. Modal synthesis is efficient and realistic yet abandons complete dynamical description and hence cannot used for certain types of performance interactions like bowing. Finite element methods are realistic and capture the behavior of the constituent physical equations but on current commodity hardware does not perform in real-time. Banded waveguides offer efficient simulations for cases for which modal synthesis is appropriate but traditional digital waveguide synthesis is not applicable. The key realization is that the dynamic behavior of traveling waves, which is being used in waveg-uide synthesis, can be applied to individual modes and that the efficient computational

### Citations

484 |
Auditory scene analysis
- Bregman
- 1990
(Show Context)
Citation Context ...question has not been thoroughly addressed in this thesis: the question of a quantitative measure of perceptual quality of a sound simulation. Perceptual models have a history of research behind them =-=[26, 130]-=-, though only recently in the context of digital speech coding for wireless telephony [4, 11, 164] led attempts at using them to arrive at quantitative measures. Another field where perceptual models ... |

341 |
An Introduction to Multigrid Methods
- Wesseling
- 1992
(Show Context)
Citation Context ...M being the width of the band [156, pp. 50-55]. Higher dimensional and more complex domains and meshes generally don’t form well-structured matrices, which increases the complexity. Multigrid methods =-=[220, 22, 190, 217, 204]-=- address this problem by solving on a coarse scale and adding refinement steps. Other possibilities are matrix preconditioning techniques [77, pp. 532ff] and manipulation (sparsing) [177]. In any case... |

304 |
An Introduction to the Psychology of Hearing
- Moore
- 1997
(Show Context)
Citation Context ...question has not been thoroughly addressed in this thesis: the question of a quantitative measure of perceptual quality of a sound simulation. Perceptual models have a history of research behind them =-=[26, 130]-=-, though only recently in the context of digital speech coding for wireless telephony [4, 11, 164] led attempts at using them to arrive at quantitative measures. Another field where perceptual models ... |

266 |
The Computer Music Tutorial
- Roads
- 1996
(Show Context)
Citation Context ... Synthesis Methods Sound synthesis methods can be classified into various groups. The previously used methods with direct relation to this work are generally classified as physical modeling synthesis =-=[165]-=- 7 . Though Roads [165] classifies many separate categories of methods within the field, many are conceptually very closely related. For instance McIntyre, Schumacher, and Woodhouse synthesis [122], K... |

204 | 3D Sound for Virtual Reality and Multimedia
- Begault
- 1994
(Show Context)
Citation Context ...tion of a delay line uses a circular buffer [55, p. 227][194, p. 258-260], something that had been realized by Schroeder and others while studying synthetic reverberation and room-acoustic simulation =-=[105, 75, 76, 12, 175, 90, 28]-=-. In the remainder of the thesis I will refer to all of these categories simply as waveguide synthesis 8 . Another line of physical modeling are methods which are concerned with direct discrete simula... |

190 | Advanced Engineering Mathematics - Kreyszig - 1988 |

161 | Graphical Modeling and Animation of Brittle Fracture - O’BRIEN, HODGINS - 1999 |

159 |
Numerical Methods for Partial Differential Equations
- Ames
- 1977
(Show Context)
Citation Context ...ty, the matrices arising in solid mechanics finite element methods have a direct mass-spring-damper interpretation [222] and finite elements can be seen as generalization of finite difference methods =-=[222, 5, 201]-=-. 8 Julius Smith often uses the term waveguide digital synthesis to indicate that this is a digital synthesis method, throughout this thesis all methods are digital in the sense that are all numerical... |

155 |
The Physics of Musical Instruments
- Fletcher, Rossing
- 1998
(Show Context)
Citation Context ...braphone bar (Measured), and simulated results for a low-resolution (Simulated 1) and high-resolution mesh (Simulated 2). The vertical lines located at 1, 4, and 10 show the tuning ratios reported in =-=[70]-=-. . . . . . . . 119 6.3 These figures show a round weight being dropped onto two different surfaces. The surface shown in (a) is rigid while the one shown in (b) is more compliant. . . . . . . . . . .... |

152 |
Theoretical Acoustics
- Morse, Ingard
- 1968
(Show Context)
Citation Context ...ly presented starting with the ideal wave-equation which describes the perfectly elastic string under tension and oscillations in an air-tube or related wave-guiding structure for small displacements =-=[132, 194, 188]-=-: ∂2y ∂t2 = c2 ∂2y ∂x2 (2.1) 10 Note that modal means different things to different communities. The word modal was chosen over spectral because the later usually is used in the Computer Music communi... |

144 |
Can one hear the shape of a drum
- Kac
- 1966
(Show Context)
Citation Context ...ng dynamics is fixed, then indeed the resulting modes directly depend on the geometry of the object. This direct connection was highlighted by Kac’s famous question “Can we hear the shape of a drum?” =-=[97]-=-. While this question has been answered with “no” for specific nonsmooth constructions [57], the question remains open for most and very general classes of shapes (like all smooth boundaries). Of cour... |

120 |
A geometrical theory of Diffraction
- Keller
(Show Context)
Citation Context ...ehavior of traveling wave rays. This is similar to using principles of geometric optics to describe general optical wave phenomena asymptotically, as Keller did in his geometric theory of diffraction =-=[102]-=-. Keller and Rubinow showed how, using their approach, circular membranes and elliptical domains can be solved. They also showed how whispering gallery modes and bouncing ball modes can be derived on ... |

110 |
Quantum Mechanics
- Schiff
- 1955
(Show Context)
Citation Context ... to conclude, that all this should not work. Surprisingly though it has repeatedly noted, that even for the long wavelength solutions (that is the low eigenmodes), the asymptotic method performs well =-=[104, 178]-=-. In specific cases, like the d’Alembert solution for the wave equation, the “asymptotic” traveling wave Ansatz is in deed exact.sCHAPTER 3. THEORY OF PROPAGATION MODELING 45 3.8.3 Discrete Asymptotic... |

104 | Perceptual coding of digital audio
- Painter, Spanias
(Show Context)
Citation Context ...telephony [4, 11, 164] led attempts at using them to arrive at quantitative measures. Another field where perceptual models have found an application is music coding for compression (for a review see =-=[144]-=-). This work has let to a wealth of experimental data which can be used as reference for recorded music quality and in fact has been used to validate computational model measures by Thiede and Kabot [... |

97 |
Geometric Scattering Theory
- Melrose
- 1995
(Show Context)
Citation Context ...athematical texts include Fedoryuk [68], Mishchenko et al [129], Kravtsov and Orlov [110] and Kozlov and Treshchëv [109]. A related text in pure mathematics treating the scattering problem is Melrose =-=[128]-=-. 131sAPPENDIX A. GLOSSARY FOR ASYMPTOTICS TERMS 132 Asymptotics an Ansatz that approaches the exact solution in the limit of some variable. More precisely the generalized Poincaré definition of asymp... |

85 |
t a Lof Acoustics
- Kinsler, Coppens, et al.
- 1982
(Show Context)
Citation Context ...locity v. If we neglect viscous shear forces then the acoustic pressure, p, of the fluid adjacent to ds is given by p = zv · ˆnˆnˆn (6.1) where z = ρc is the fluid’s specific acoustic impedance. From =-=[106]-=-, the density of air at 20 ◦ C under one atmosphere of pressure is ρ = 1.21kg/m 3 , and the acoustic wave speed is c = 343m/s, giving z = 415Pa · s/m. Representing the pressure field over Ω requires s... |

82 | FoleyAutomatic: physically-based sound effects for interactive simulation and animation
- Doel, Kry, et al.
- 2001
(Show Context)
Citation Context ...nd memory operations and hence the complexity of this method is O(M) where M is the number of modes modeled. Models of this type have been proposed and used on a wide array of objects and instruments =-=[1, 219, 182, 210, 211, 208, 212]-=-. Modal shapes can be reconstructed by sampling the surface and hence reconstructing the modal shapes through an amplitude envelope as function of discrete space k. Then equation 2.5 becomes: y(n,k) =... |

72 | A Beam Tracing Approach to Acoustic Modeling for Interactive Virtual Environments
- FUNKHOUSER, CARLBOM, et al.
- 1998
(Show Context)
Citation Context ...various delays of time, having traveled down and up the valley and reflected at the side and bottom rocks. Propagation modeling of this kind is well known and widely used in simulating room acoustics =-=[74]-=-. As this propagation has a finite speed, air as a medium delays the arrival of sound . The listener does not really know or need to know the peculiarities or details of the sound traveling down the v... |

69 |
Real Sound Synthesis for Interactive Applications
- Cook
- 2002
(Show Context)
Citation Context ...o this thesis. They are either not physically informed or stochastic in nature. 2.1.1 Waveguide Synthesis There are many references which comprehensively review waveguide modeling at it current state =-=[188, 45, 187, 99, 8]-=-. Here I will only mention the key ideas and ideas in previous work that are relevant to this thesis. The derivation of the basic idea of waveguide synthesis is usually presented starting with the ide... |

66 | Bark and ERB bilinear transforms
- Smith, Abel
- 1999
(Show Context)
Citation Context ...], instrument body modeling [148],sCHAPTER 5. COMPARISON WITH ALTERNATIVE METHODS 107 reduction of dispersion errors in waveguide meshes [176], audio effects [66, 65] and warping to perceptual scales =-=[189]-=-. Otherwise it’s been studied in modeling room impulse responses [107], the approximation of delay systems [119, 118] and system identification [7, 67], filter design [120], and image compression [41]... |

61 | Synthesizing Sounds From Physically Based Motion
- O’BRIEN, COOK, et al.
- 2001
(Show Context)
Citation Context ...ference [60] and in its final version in an archival journal [61]. Finite element methods (section 6) as an alternative solution was investigated in collaboration with James F. O’Brien of UC Berkeley =-=[135]-=-. The generalization of the simulation method to higher dimensions and examples of Tabla and glass harmonica (sections 4.2 and 4.3) is in preparation for publication, though an unpublished manuscript ... |

58 |
On the oscillations of musical instruments
- McIntyre, Schumacher, et al.
- 1983
(Show Context)
Citation Context ...is [165] 7 . Though Roads [165] classifies many separate categories of methods within the field, many are conceptually very closely related. For instance McIntyre, Schumacher, and Woodhouse synthesis =-=[122]-=-, Karplus-Strong synthesis [100, 94, 99, 195] and Waveguide synthesis [185, 187, 188] are very closely related. In the case of Karplus-Strong and Waveguide synthesis, the difference can be seen as mer... |

58 |
Principles of digital waveguide models of musical instruments
- Smith
- 1998
(Show Context)
Citation Context ... (also referred to as scatterings) and discrete traveling wave approximations has been extensively studied by Bilbao [18], though the treatment of scattering in digital waveguides has long been known =-=[188]-=-. Hence one critical area for performance gain is the process of deciding when neglecting local scattering in a dynamical system is a good approximation or not (physically and perceptually). Finally, ... |

57 |
III, Techniques for digital filter design and system identification with application to the violin
- Smith
- 1983
(Show Context)
Citation Context ...ods within the field, many are conceptually very closely related. For instance McIntyre, Schumacher, and Woodhouse synthesis [122], Karplus-Strong synthesis [100, 94, 99, 195] and Waveguide synthesis =-=[185, 187, 188]-=- are very closely related. In the case of Karplus-Strong and Waveguide synthesis, the difference can be seen as mere physical interpretation of the structure. 6 The confusion that is part of the proce... |

56 |
Splitting the unit delay: Tools for fractional delay filter design
- Laakso, Välimäki, et al.
- 1996
(Show Context)
Citation Context ...solution suffers. A filter modeling the fractional delay between the integer and actual length is necessary to accommodate for accurate tuning. This is achieved by using various filter design methods =-=[114]-=-. A number of modifications to the ideal waveguide have been studied to achieve more realism. These extensions are usually instrument-specific, for instance the bending stiffness of piano-strings [6] ... |

52 |
Physically informed sonic modeling (PhISM): synthesis of percussive sounds
- Cook
- 1997
(Show Context)
Citation Context ...hods. The final class of physical synthesis methods models the spectral response of the physical system. This can be achieved through additive sinusoidal modeling [182, 210], resonant filter modeling =-=[219, 42]-=-, model decomposition modeling [1]. These methods will be referred to as modal synthesis methods throughout this thesis 10 . Other methods of sound synthesis, as found in [165, 55, 198], are not revie... |

51 |
CORDIS-ANIMA: a modeling and simulation system for sound and image synthesis – The General formalism
- Cadoz, Luciani, et al.
- 1993
(Show Context)
Citation Context ...cerned with direct discrete simulation of the local dynamics responsible for sound generation. These include methods using finite differencing [174, 33, 20, 113, 35, 146], mass-spring-damper networks =-=[30, 71, 31]-=-, finite element methods [27, 21] and transmission-line methods for solving differential equations [214, 18] (which includes highly scattering digital waveguide simulations and digital waveguide meshe... |

45 |
PluckedString Models: From the Karplus-Strong Algorithm to Digital Waveguides and Beyond
- Karjalainen, Valimaki, et al.
- 1998
(Show Context)
Citation Context ...classifies many separate categories of methods within the field, many are conceptually very closely related. For instance McIntyre, Schumacher, and Woodhouse synthesis [122], Karplus-Strong synthesis =-=[100, 94, 99, 195]-=- and Waveguide synthesis [185, 187, 188] are very closely related. In the case of Karplus-Strong and Waveguide synthesis, the difference can be seen as mere physical interpretation of the structure. 6... |

41 | Physical modeling synthesis update
- Smith
- 1996
(Show Context)
Citation Context ...ods within the field, many are conceptually very closely related. For instance McIntyre, Schumacher, and Woodhouse synthesis [122], Karplus-Strong synthesis [100, 94, 99, 195] and Waveguide synthesis =-=[185, 187, 188]-=- are very closely related. In the case of Karplus-Strong and Waveguide synthesis, the difference can be seen as mere physical interpretation of the structure. 6 The confusion that is part of the proce... |

39 | Discrete Representation of Signals - Oppenheim, Johnson - 1972 |

38 |
o Rsponsive Input Devices and sound synthesis by simulation of instrumeatal mechanisms : The CORDIS system
- Cadoz, Luciani
(Show Context)
Citation Context ...cerned with direct discrete simulation of the local dynamics responsible for sound generation. These include methods using finite differencing [174, 33, 20, 113, 35, 146], mass-spring-damper networks =-=[30, 71, 31]-=-, finite element methods [27, 21] and transmission-line methods for solving differential equations [214, 18] (which includes highly scattering digital waveguide simulations and digital waveguide meshe... |

38 |
Computation of spectra with unequal resolution using the fast Fourier transform
- Oppenheim, Johnson
- 1971
(Show Context)
Citation Context ...g all-passes [94, 145, 213, 169].sCHAPTER 5. COMPARISON WITH ALTERNATIVE METHODS 106 The effect of all-passes replacing unit delays has long been studied and is classically known as frequency warping =-=[192, 141, 140]-=- and has been theoretically extended until today [216, 138]. To the same extent, the design of all-passes has been extensively studied [116]. An example design can be seen in Figure 5.2. Figure 5.2: A... |

33 |
III, “Efficient simulation of the reed-bore and bow-string mechanisms
- Smith
- 1986
(Show Context)
Citation Context ... Galilei [69, p. 331] and Chladni [40, 218]. Keller and Friedlander [101, 72] independly developed the theoretical foundation of bowed string action that is still widely used today in computer models =-=[186, 133]-=-. The friction characteristic used by them is depicted in Figure 2.8. This basic model was subject to a number of refinements taking into account the width of the bow [150] and more detailed knowledge... |

31 |
Numerical simulation of piano strings. i. a physical model for a struck string usin finite-difference methods
- Chaigne, Askenfelt
- 1994
(Show Context)
Citation Context ...other line of physical modeling are methods which are concerned with direct discrete simulation of the local dynamics responsible for sound generation. These include methods using finite differencing =-=[174, 33, 20, 113, 35, 146]-=-, mass-spring-damper networks [30, 71, 31], finite element methods [27, 21] and transmission-line methods for solving differential equations [214, 18] (which includes highly scattering digital wavegui... |

31 |
On the Sensations of
- Helmholtz
- 1877
(Show Context)
Citation Context ...omogeneity or scaling property and linearity between multiple different inputs, called additivity or superposition property [142, p. 18].sCHAPTER 2. BACKGROUND, PREVIOUS ART AND RELATED LITERATURE 21 =-=[86]-=-, Rayleigh [159] and Raman [158], though early mentionings can be traced back to Galileo Galilei [69, p. 331] and Chladni [40, 218]. Keller and Friedlander [101, 72] independly developed the theoretic... |

31 | Digital synthesis of plucked string and drum timbres - Karplus, Strong - 1983 |

31 |
Reducing the dispersion error in the digital waveguide mesh using interpolation and frequency warping techniques
- Savioja, Välimäki
- 2000
(Show Context)
Citation Context ...is approach has been used for warped linear prediction [117, 80], instrument body modeling [148],sCHAPTER 5. COMPARISON WITH ALTERNATIVE METHODS 107 reduction of dispersion errors in waveguide meshes =-=[176]-=-, audio effects [66, 65] and warping to perceptual scales [189]. Otherwise it’s been studied in modeling room impulse responses [107], the approximation of delay systems [119, 118] and system identifi... |

31 |
Synthesis of shape dependent sounds with physical modeling
- DOEL, PAI
- 1996
(Show Context)
Citation Context ...are well approximated as linear interactions are plucking, striking with hard beaters, and certain rough-surface frictions (those where the rough surface action is well-captured by periodic impulses) =-=[210, 209, 211, 163, 54, 208, 212, 42, 44]-=-. 2.2.2 Nonlinear or Friction Excitation: Stick-Slip For non-linear excitations, the independence of excitation from the state of the system does not hold. Hence, the force function of a finite elemen... |

30 |
Friction-induced vibration, chatter, squeal, and chaos. part ii: Dynamics and modeling
- Ibrahim
- 1994
(Show Context)
Citation Context ...[3, 87], propeller shaft noises of submarines [69], squeals from train wheels [84, 82], friction noises of insects [3] and more. This lead to a large corpus of research on friction-based oscillations =-=[3, 69, 92, 93, 203, 123, 87]-=- is concerned with the damping (i.e. control) of such sounds [83, 197, 85].sChapter 3 Theory of Propagation Modeling Using Banded Digital Waveguides Anybody else could have told me this in advance, bu... |

29 |
Theoretical Global Seismology
- Dahlen, Tromp
- 1998
(Show Context)
Citation Context ...pansions [78, 24, 89] [38, for a detailed review including the connection to rays in two-dimensional elasticity dynamics]. Asymptotic methods find applications in electromagnetics [23] and seismology =-=[50]-=-. 3.8.2 Previous Work: Numerical Traveling Waves Numerical simulation based on discretized traveling waves for the wave-equation was accidentally discovered by Karplus and Strong and was quickly recog... |

29 |
Extensions of the KarplusStrong plucked string algorithm
- Jaffe, Smith
- 1983
(Show Context)
Citation Context ...classifies many separate categories of methods within the field, many are conceptually very closely related. For instance McIntyre, Schumacher, and Woodhouse synthesis [122], Karplus-Strong synthesis =-=[100, 94, 99, 195]-=- and Waveguide synthesis [185, 187, 188] are very closely related. In the case of Karplus-Strong and Waveguide synthesis, the difference can be seen as mere physical interpretation of the structure. 6... |

28 |
A simplified approach to modeling dispersion caused by stiffness in strings and plates
- Duyne, Smith
- 1994
(Show Context)
Citation Context ...priately modeled phase delay response. This has strong connections to and is strongly motivated by work on modeling the stiffness of strings starting from a mixed bar/string-equation using all-passes =-=[94, 145, 213, 169]-=-.sCHAPTER 5. COMPARISON WITH ALTERNATIVE METHODS 106 The effect of all-passes replacing unit delays has long been studied and is classically known as frequency warping [192, 141, 140] and has been the... |

27 |
The Physics of the Violin
- Cremer
- 1984
(Show Context)
Citation Context ...l shapes is not necessary. These results have also been presented in [60, 61, 62] and will be discussed in section 4.1.1. The action of the bow has only been studied extensively when exciting strings =-=[46, 91]-=-. Numerical simulations typically use an efficient time-domain approach which follows from the constant phase delay characteristic of the ideal string equation [122] and this approach has been refined... |

27 |
Accurate dispersion simulation for piano strings
- Rocchesso, Scalcon
- 1996
(Show Context)
Citation Context ...priately modeled phase delay response. This has strong connections to and is strongly motivated by work on modeling the stiffness of strings starting from a mixed bar/string-equation using all-passes =-=[94, 145, 213, 169]-=-.sCHAPTER 5. COMPARISON WITH ALTERNATIVE METHODS 106 The effect of all-passes replacing unit delays has long been studied and is classically known as frequency warping [192, 141, 140] and has been the... |

27 | The sounds of physical shapes
- Doel, Pai
- 1998
(Show Context)
Citation Context ...ds as instances of finite element methods. The final class of physical synthesis methods models the spectral response of the physical system. This can be achieved through additive sinusoidal modeling =-=[182, 210]-=-, resonant filter modeling [219, 42], model decomposition modeling [1]. These methods will be referred to as modal synthesis methods throughout this thesis 10 . Other methods of sound synthesis, as fo... |

25 |
D.V.: Billiards: A genetic introduction to the dynamics of systems with impacts
- Kozlov, Treshchëv
- 1991
(Show Context)
Citation Context ...o its starting position. In this setting, billiard balls are said to be on periodic orbits, and the task is to find them. On the convex domain this has been extensively studied by Birkhoff and others =-=[109]-=-. Traveling wave methods are studied and used in a vast number of different fields. Connected theories have been developed in pure and applied mathematics with various motivations. Among them is appro... |

24 |
Wave and Scattering Methods for the Numerical Integration of Partial Differential Equations
- Bilbao
- 2004
(Show Context)
Citation Context ...thods using finite differencing [174, 33, 20, 113, 35, 146], mass-spring-damper networks [30, 71, 31], finite element methods [27, 21] and transmission-line methods for solving differential equations =-=[214, 18]-=- (which includes highly scattering digital waveguide simulations and digital waveguide meshes) 9 . Again, all these methods are in close relationship to each other. For instance there is a direct corr... |

23 | Circulant and Elliptic Feedback Delay Networks for Artificial Reverberation
- Rocchesso, Smith
- 1997
(Show Context)
Citation Context ... 25] 15 (see section 3.9.2). The connection between resonant response of a geometric enclosure and its shape has been realized in this context by Rocchesso and studied in the case of a rectangular box=-=[166, 170]-=- and sphere [167]. An alternate approach to multi-path feedback delay models can be found in the literature of physical models for string instruments. Until recently, the focus has been on studying ph... |

22 |
Asymptotic solution of eigenvalue problems
- Keller, Rubinow
- 1960
(Show Context)
Citation Context ...c approximations to dynamical systems, in particular optics and quantum mechanics. In joint work with Rubinow, he proposed that the method could also be used to find eigenvalues of dynamic structures =-=[104]-=-. The idea is that the dynamics of structures can be well approximated by studying the behavior of traveling wave rays. This is similar to using principles of geometric optics to describe general opti... |

21 |
The Missing Link: Modal Synthesis
- Adrien
- 1991
(Show Context)
Citation Context ...sis methods models the spectral response of the physical system. This can be achieved through additive sinusoidal modeling [182, 210], resonant filter modeling [219, 42], model decomposition modeling =-=[1]-=-. These methods will be referred to as modal synthesis methods throughout this thesis 10 . Other methods of sound synthesis, as found in [165, 55, 198], are not reviewed as they are not relevant to th... |