Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Redes de rega de campos de golfe

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Studies on the universal parameters and onset of chaos in dissipative systems

It has become clear over the last few years that many deterministic dynamical systems described by simple but nonlinear equations with only a few variables can behave in an irregular or random fashion. This phenomenon, commonly called deterministic chaos, is essentially due to the fact that we cannot deal with infinitely precise numbers. In these systems trajectories emerging from nearby initial conditions diverge exponentially as time evolves)and therefore)any small error in the initial measurement spreads with time considerably, leading to unpredictable and chaotic behaviour The thesis work is mainly centered on the asymptotic behaviour of nonlinear and nonintegrable dissipative dynamical systems. It is found that completely deterministic nonlinear differential equations describing such systems can exhibit random or chaotic behaviour. Theoretical studies on this chaotic behaviour can enhance our understanding of various phenomena such as turbulence, nonlinear electronic circuits, erratic behaviour of heart and brain, fundamental molecular reactions involving DNA, meteorological phenomena, fluctuations in the cost of materials and so on. Chaos is studied mainly under two different approaches - the nature of the onset of chaos and the statistical description of the chaotic state.

Quantified real constraint solving using modal intervals with applications to control

Les restriccions reals quantificades (QRC) formen un formalisme matemàtic utilitzat per modelar un gran nombre de problemes físics dins els quals intervenen sistemes d'equacions no-lineals sobre variables reals, algunes de les quals podent ésser quantificades. Els QRCs apareixen en nombrosos contextos, com l'Enginyeria de Control o la Biologia. La resolució de QRCs és un domini de recerca molt actiu dins el qual es proposen dos enfocaments diferents: l'eliminació simbòlica de quantificadors i els mètodes aproximatius. Tot i això, la resolució de problemes de grans dimensions i del cas general, resten encara problemes oberts. Aquesta tesi proposa una nova metodologia aproximativa basada en l'Anàlisi Intervalar Modal, una teoria matemàtica que permet resoldre problemes en els quals intervenen quantificadors lògics sobre variables reals. Finalment, dues aplicacions a l'Enginyeria de Control són presentades. La primera fa referència al problema de detecció de fallades i la segona consisteix en un controlador per a un vaixell a vela.

A comparison of two methods for developing the linearization of a shallow-water model

We develop the linearization of a semi-implicit semi-Lagrangian model of the one-dimensional shallow-water equations using two different methods. The usual tangent linear model, formed by linearizing the discrete nonlinear model, is compared with a model formed by first linearizing the continuous nonlinear equations and then discretizing. Both models are shown to perform equally well for finite perturbations. However, the asymptotic behaviour of the two models differs as the perturbation size is reduced. This leads to difficulties in showing that the models are correctly coded using the standard tests. To overcome this difficulty we propose a new method for testing linear models, which we demonstrate both theoretically and numerically. © Crown copyright, 2003. Royal Meteorological Society

Cyclic extrusion of a lava dome based on a stick-slip mechanism

Lava dome eruptions are sometimes characterised by large periodic fluctuations in extrusion rate over periods of hours that may be accompanied by Vulcanian explosions and pyroclastic flows. We consider a simple system of nonlinear equations describing a 1D flow of lava extrusion through a deep elastic dyke feeding a shallower cylindrical conduit in order to simulate this short-period cyclicity. Stick-slip conditions depending on a critical shear stress are assumed at the wall boundary of the cylindrical conduit. By analogy with the behaviour of industrial polymers in a plastic extruder, the elastic dyke acts like a barrel and the shallower cylindrical portion of the conduit as a die for the flow of magma acting as a polymer. When we applied the model to the Soufrière Hills Volcano, Montserrat, for which the key parameters have been evaluated from previous studies, cyclic extrusions with periods from 3 to 30 h were readily simulated, matching observations. The model also reproduces the reduced period of cycles observed when a major unloading event occurs due to lava dome collapse.

A new method for locating the global optimum : application of the cutting angle method to molecular structure prediction

Many problems in chemistry depend on the ability to identify the global minimum or maximum of a function. Examples include applications in chemometrics, optimization of reaction or operating conditions, and non-linear least-squares analysis. This paper presents the results of the application of a new method of deterministic global optimization, called the cutting angle method (CAM), as applied to the prediction of molecular geometries. CAM is shown to be competitive with other global optimization techniques for several benchmark molecular conformation problem. CAM is a general method that can also be applied to other computational problems involving global minima, global maxima or finding the roots of nonlinear equations.

General solutions for some classes of interacting two field kinks

In this work we present some classes of models whose the corresponding two coupled first-order nonlinear equations can be put into a linear form, and consequently be solved completely. In these cases the so-called trial orbit method is completely unnecessary. We recall that some physically important models as, for instance, the problem of tiling a plane with a network of defects and polymer properties are in this class of models. (c) 2005 Elsevier B.V. All rights reserved.

Dynamics of Bose-Einstein condensates in cigar-shaped traps

The Gross-Pitaevskii equation for a Bose-Einstein condensate confined in an elongated cigar-shaped trap is reduced to an effective system of nonlinear equations depending on only one space coordinate along the trap axis. The radial distribution of the condensate density and its radial velocity are approximated by Gaussian functions with real and imaginary exponents, respectively, with parameters depending on the axial coordinate and time. The effective one-dimensional system is applied to a description of the ground state of the condensate, to dark and bright solitons, to the sound and radial compression waves propagating in a dense condensate, and to weakly nonlinear waves in repulsive condensate. In the low-density limit our results reproduce the known formulas. In the high-density case our description of solitons goes beyond the standard approach based on the nonlinear Schrodinger equation. The dispersion relations for the sound and radial compression waves are obtained in a wide region of values of the condensate density. The Korteweg-de Vries equation for weakly nonlinear waves is derived and the existence of bright solitons on a constant background is predicted for a dense enough condensate with a repulsive interaction between the atoms.

Numerical study of the spherically symmetric Gross-Pitaevskii equation in two space dimensions

The Gross-Pitaevskii equation for Bose-Einstein condensation (BEC) in two space dimensions under the action of a harmonic oscillator trap potential for bosonic atoms with attractive and repulsive interparticle interactions was numerically studied by using time-dependent and time-independent approaches. In both cases, numerical difficulty appeared for large nonlinearity. Nonetheless, the solution of the time-dependent approach exhibited intrinsic oscillation with time iteration which is independent of space and time steps used in discretization.

Variational calculation of positronium-helium-atom scattering length

A basis-set calculation scheme for S-waves Ps-He elastic scattering below the lowest inelastic threshold was formulated using a variational expression for the transition matrix. The scheme was illustrated numerically by calculating the scattering length in the electronic doublet state: a=1.0±0.1 a.u.

Soliton propagation in a medium with Kerr nonlinearity and resonant impurities: A variational approach

The soliton propagation in a medium with Kerr nonlinearity and resonant impurities was studied by a variational approach. The existence of a solitary wave was shown within the framework of a combined nonintegrable system composed of one nonlinear Schrödinger and a pair of Bloch equations. The analytical solution which was obtained, was tested through numerical simulations confirming its solitary wave nature.

Array of Bose-Einstein condensates under time-periodic Feshbach-resonance management

The investigation of the dynamics of a discrete soliton in an array of Bose-Einstein condensates under the action of a periodically time-modulated atomic scattering length [Feshbach-resonance management (FRM)] was discussed. The slow and rapid modulations, in comparison with the tunneling frequency were considered. An averaged equation, which was a generalized discrete nonlinear Schrödinger equation, including higher-order effective nonlinearities and intersite nonlinear interactions was derived in the case of the rapid modulation. It was demonstrated that the modulations of sufficient strength results in splitting of the soliton by direct simulations.

Optimal Boussinesq model for shallow-water waves interacting with a microstructure

In this paper, we consider the propagation of water waves in a long-wave asymptotic regime, when the bottom topography is periodic on a short length scale. We perform a multiscale asymptotic analysis of the full potential theory model and of a family of reduced Boussinesq systems parametrized by a free parameter that is the depth at which the velocity is evaluated. We obtain explicit expressions for the coefficients of the resulting effective Korteweg-de Vries (KdV) equations. We show that it is possible to choose the free parameter of the reduced model so as to match the KdV limits of the full and reduced models. Hence the reduced model is optimal regarding the embedded linear weakly dispersive and weakly nonlinear characteristics of the underlying physical problem, which has a microstructure. We also discuss the impact of the rough bottom on the effective wave propagation. In particular, nonlinearity is enhanced and we can distinguish two regimes depending on the period of the bottom where the dispersion is either enhanced or reduced compared to the flat bottom case. © 2007 The American Physical Society.