Physical Models of Wind Instruments : A Generalized Excitation Coupled with a Modular Tube Simulation Platform

To develop real-time simulations of wind instruments, digital waveguides filters can be used as an efficient representation of the air column. Many aerophones are shaped as horns which can be approximated using conical sections. Therefore the derivation of conical waveguide filters is of special interest. When these filters are used in combination with a generalized reed excitation, several classes of wind instruments can be simulated. In this paper we present the methods for transforming a continuous description of conical tube segments to a discrete filter representation. The coupling of the reed model with the conical waveguide and a simplified model of the termination at the open end are described in the same way. It turns out that the complete lossless conical waveguide requires only one type of filter.Furthermore, we developed a digital reed excitation model, which is purely based on numerical integration methods, i.e., without the use of a look-up table.

An overview of seismic hybrid testing of engineering structures

An overview of research on the development of the hybrid test method is presented. The maturity of the hybrid test method is mapped in order to provide context to individual research in the overall development of the test method. In the pseudo dynamic (PsD) test method, the equations of motion are solved using a time stepping numerical integration technique with the inertia and damping being numerically modelled whilst restoring force is physically measured over an extended timescale. Developments in continuous PsD testing led to the real-time hybrid test method and geographically distributed hybrid tests. A key aspect to the efficiency of hybrid testing is the substructuring technique where the critical structural subassemblies that are fundamental to the overall response of the structure are physically tested whilst the remainder of the structure whose response can be more easily predicted is numerically modelled. Much of the early research focused on developing the accuracy and efficiency of the test method, whereas more recently the method has matured to a level where the test method is applied purely as a dynamic testing technique. Developments in numerical integration methods, substructuring, experimental error reduction, delay compensation and speed of testing have led to a test method now in use as full-scale real-time dynamic testing method that is reliable, accurate, efficient and cost effective.

Optimal Scheduling Methods to Integrate Plug-In Electric Vehicles with the Power System: A Review

The introduction of the Tesla in 2008 has demonstrated to the public of the potential of electric vehicles in terms of reducing fuel consumption and green-house gas from the transport sector. It has brought electric vehicles back into the spotlight worldwide at a moment when fossil fuel prices were reaching unexpected high due to increased demand and strong economic growth. The energy storage capabilities from of fleets of electric vehicles as well as the potentially random discharging and charging offers challenges to the grid in terms of operation and control. Optimal scheduling strategies are key to integrating large numbers of electric vehicles and the smart grid. In this paper, state-of-the-art optimization methods are reviewed on scheduling strategies for the grid integration with electric vehicles. The paper starts with a concise introduction to analytical charging strategies, followed by a review of a number of classical numerical optimization methods, including linear programming, non-linear programming, dynamic programming as well as some other means such as queuing theory. Meta-heuristic techniques are then discussed to deal with the complex, high-dimensional and multi-objective scheduling problem associated with stochastic charging and discharging of electric vehicles. Finally, future research directions are suggested.

Molecular Distributions in the Inner Regions of Protostellar Disks

The distributions of molecules in the inner regions of a protostellar disk are presented. These were calculated using an uncoupled chemical/dynamical model, with a numerical integration of the vertical disk structure. A comparison between models with and without the effects of X-ray ionisation is made, and molecules are identified which are good tracers of the ionisation level in this part of the disk, notably CN and C_2H. In the region considered in this paper (r

A modified commercial DRIFTS cell for kinetically relevant operando studies of heterogeneous catalytic reactions

This paper discusses a number of checks that should be carried out to ensure that the kinetic and spectroscopic measurements made using a DRIFTS cell are meaningful. The observations reported here demonstrate how an appropriately modified commercial DRIFTS cell can provide pertinent kinetic information about both gaseous products and the related surface intermediates. The oxidation of CO with 02 was used as a test to assess the catalyst bed bypass by the reaction mixture. Full CO conversion was obtained after the light-off temperature in the case of the modified cell, contrary to the case of the original cell, for which 80% of the reaction mixture bypassed the catalyst bed. The water-gas shift reaction over a Pt/CeO2 catalyst was used as a model reaction to further characterize the behavior of the cell under reaction conditions. The catalyst bed was shown not to be a dead-zone and was purged in essentially the same time as that needed to purge the cell. The reaction chamber globally operated in a quasi plug-flow mode and the gas composition in the thin catalyst bed appears to be homogeneous when operated under differential conditions. The production of the gas-phase reaction product CO2 could be simultaneously followed both by mass spectrometry and DRIFTS, both techniques leading to identical results. Various IR bands integration methods were discussed to allow a precise and accurate determination of the surface concentration of adsorbates during isotopic exchange. (c) 2008 Elsevier B.V. All rights reserved.

Convergence acceleration of the doubly periodic Green's function for the analysis of thin wire arrays

A method is proposed to accelerate the evaluation of the Green's function of an infinite double periodic array of thin wire antennas. The method is based on the expansion of the Green's function into series corresponding to the propagating and evanescent waves and the use of Poisson and Kummer transformations enhanced with the analytic summation of the slowly convergent asymptotic terms. Unlike existing techniques the procedure reported here provides uniform convergence regardless of the geometrical parameters of the problem or plane wave excitation wavelength. In addition, it is numerically stable and does not require numerical integration or internal tuning parameters, since all necessary series are directly calculated in terms of analytical functions. This means that for nonlinear problem scenarios that the algorithm can be deployed without run time intervention or recursive adjustment within a harmonic balance engine. Numerical examples are provided to illustrate the efficiency and accuracy of the developed approach as compared with the Ewald method for which these classes of problems requires run time splitting parameter adaptation.

A heuristic model for the simulation of the deformation of elastic and spongy material for virtual reality applications

This paper presents a practical algorithm for the simulation of interactive deformation in a 3D polygonal mesh model. The algorithm combines the conventional simulation of deformation using a spring-mass-damping model, solved by explicit numerical integration, with a set of heuristics to describe certain features of the transient behaviour, to increase the speed and stability of solution. In particular, this algorithm was designed to be used in the simulation of synthetic environments where it is necessary to model realistically, in real time, the effect on non-rigid surfaces being touched, pushed, pulled or squashed. Such objects can be solid or hollow, and have plastic, elastic or fabric-like properties. The algorithm is presented in an integrated form including collision detection and adaptive refinement so that it may be used in a self-contained way as part of a simulation loop to include human interface devices that capture data and render a realistic stereoscopic image in real time. The algorithm is designed to be used with polygonal mesh models representing complex topology, such as the human anatomy in a virtual-surgery training simulator. The paper evaluates the model behaviour qualitatively and then concludes with some examples of the use of the algorithm.

Double-Electron Above-Threshold Ionization Resonances as Interference Phenomena

We report calculations of double ionization energy spectra and momentum distributions of laser-driven helium due to few-cycle pulses of wavelength 195 nm. The results are obtained from full-dimensional numerical integration of the two electron time-dependent Schr¨odinger equation. A momentum-space analysis of doubly ionizing wavepackets shows that the concentric-ring structure of above-threshold double ionization, together with the associated structure of peaks in the total kinetic energy spectrum, may be attributed to wavepacket interference effects, where at least two doubly-ionizing wavepackets from different recollision events populate the same spatial hemisphere.

Superluminal electromagnetic solitary waves in electron-positron plasmas

The propagation of an electromagnetic wave packet in an electron-positron plasma, in the form of coupled localized electromagnetic excitations, is investigated, from first principles. By means of the Poincare section method, a special class of superluminal localized nonlinear stationary solutions, existing along a separatrix curve, are proposed as intrinsic electromagnetic modes in a relativistic electron-positron plasma. The ratio of the envelope time scale to the carrier wave time scale of these envelope solitary waves critically depends on the carrier's phase velocity. In the strongly superluminal regime, v(ph)/c >> 1, the large difference between the envelope and carrier time scales enables us to carry out a multiscale perturbative analysis resulting in an analytical form of the solution envelope. The analytical prediction thus obtained is shown to be in agreement with the solution obtained via a direct numerical integration. Copyright (c) EPLA, 2012

Vibration based seismic damage identification procedure for a steel concentrically braced frame structure

An approach for seismic damage identification of a single-storey steel concentrically braced frame (CBF) structure is presented through filtering and double integration of a recorded acceleration signal. A band-pass filter removes noise from the acceleration signal followed by baseline correction being used to reduce the drift in velocity and displacement during numerical integration. The pre-processing achieves reliable numerical integration that predicts the displacement response accurately when compared to the measured lateral in-plane displacement of the CBF structure. The lateral displacement of the CBF structure is used to infer buckling and yielding of bracing members through seismic tests. The level of interstorey drift of the CBF during a seismic excitation allows the yield and buckling of the bracing members to be identified and indirectly detects damage based on exceedance of calculated displacement limits. The calculated buckling and yielding displacement threshold limits used to identify damage are demonstrated to accurately identify initial buckling and yielding in the bracing members.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Thermodynamic Properties of Dichloromethane, Bromochloromethane, and Dibromomethane under Elevated Pressure: Experimental Results and SAFT-VR Mie Predictions

The speeds of sound in dibromomethane, bromochloromethane, and dichloromethane have been measured in the temperature range from 293.15 to 313.15 K and at pressures up to 100 MPa. Densities and isobaric heat capacities at atmospheric pressure have been also determined. Experimental results were used to calculate the densities and isobaric heat capacities as the function of temperature and pressure by means of a numerical integration technique. Moreover, experimental data at atmospheric pressure were then used to determine the SAFT-VR Mie molecular parameters for these liquids. The accuracy of the model has been then evaluated using a comparison of derived experimental high-pressure data with those predicted using SAFT. It was found that the model provide the possibility to predict also the isobaric heat capacity of all selected haloalkanes within an error up to 6%.