Design Methodology for Real-Time FPGA-Based Sound Synthesis

Explicit finite difference (FD) schemes can realise highly realistic physical models of musical instruments but are computationally complex. A design methodology is presented for the creation of FPGA-based micro-architectures for FD schemes which can be applied to a range of applications with varying computational requirements, excitation and output patterns and boundary conditions. It has been applied to membrane and plate-based sound producing models, resulting in faster than real-time performance on a Xilinx XC2VP50 device which is 10 to 35 times faster than general purpose and DSP processors. The models have developed in such a way to allow a wide range of interaction (by a musician) thereby leading to the possibility of creating a highly realistic digital musical instrument.

Imaging domains in BaTiO3 single crystal nanostructures: comparing information from transmission electron microscopy and piezo-force microscopy

This article compares and contrasts information
obtained, using transmission electron microscopy (TEM)
and piezo-force microscopy (PFM), on domain configurations
adopted in single crystal lamellae of BaTiO3, that had
been cut directly from bulk using a focused ion beam
microscope with top and bottom surfaces parallel to
{100}pseudocubic. Both forms of imaging reveal domain
walls parallel to {110}pseudocubic, consistent with sets of 90
domains with dipoles oriented parallel to the two
\001[pseudocubic directions in the plane of the lamellae.
However, the domain width was observed to be dramatically
larger using PFM than it was using TEM. This suggests
significant differences in the surface energy densities
that drive the domain formation in the first place, that could
relate to differences in the boundary conditions in the two
modes of imaging (TEM samples are imaged under high
vacuum, whereas PFM imaging was performed in air).
Attempts were made to map local dipole orientations
directly, using a form of ‘vector’ PFM. However, information
inferred was largely inconsistent with the known
crystallography of the samples, raising concern about the
levels of care needed for accurate interpretation of PFM
images.

Simple and accurate analytical model of planar grids and high-impedance surfaces, comprising metal strips or patches

Simple analytical formulas are introduced for the grid impedance of electrically dense arrays of square patches and for the surface impedance of high-impedance surfaces based on the dense arrays of metal strips or square patches over ground planes. Emphasis is on the oblique-incidence excitation. The approach is based on the known analytical models for strip grids combined with the approximate Babinet principle for planar grids located at a dielectric interface. Analytical expressions for the surface impedance and reflection coefficient resulting from our analysis are thoroughly verified by full-wave simulations and compared with available data in open literature for particular cases. The results can be used in the design of various antennas and microwave or millimeter wave devices which use artificial impedance surfaces and artificial magnetic conductors (reflect-array antennas, tunable phase shifters, etc.), as well as for the derivation of accurate higher-order impedance boundary conditions for artificial (high-) impedance surfaces. As an example, the propagation properties of surface waves along the high-impedance surfaces are studied.

Analysis of compressive membrane action in concrete slabs

This paper summarises the results obtained from non-linear finite-element analysis (NLFEA) of a series of reinforced-concrete one-way slabs with various boundary conditions representative of a bridge deck slab strip in which compressive membrane action governs the structural behaviour. The application of NLFEA for the optimum analysis and design of in-plane restrained concrete slabs is explored. An accurate material model and various equation solution methods were assessed to find a suitable finite-element method for the analysis of concrete slabs in which arching action occurs. Finally, the results from the NLFEA are compared and validated with those from various experimental test data. Significantly, the numerical analysis was able to model the arching action that occurred as a result of external in-plane restraint at the supports and which enhanced the ultimate strength of the slab. The NLFEA gave excellent predictions for the ultimate load-carrying capacity and far more accurate predictions than those obtained using standard flexural or elastic theory.

Half-precessional dynamics of monsoon rainfall near the East African Equator

External climate forcings-such as long-term changes in solar insolation-generate different climate responses in tropical and high latitude regions(1). Documenting the spatial and temporal variability of past climates is therefore critical for understanding how such forcings are translated into regional climate variability. In contrast to the data-richmiddle and high latitudes, high-quality climate-proxy records from equatorial regions are relatively few(2-4), especially from regions experiencing the bimodal seasonal rainfall distribution associated with twice-annual passage of the Intertropical Convergence Zone. Here we present a continuous and well-resolved climate-proxy record of hydrological variability during the past 25,000 years from equatorial East Africa. Our results, based on complementary evidence from seismic-reflection stratigraphy and organic biomarker molecules in the sediment record of Lake Challa near Mount Kilimanjaro, reveal that monsoon rainfall in this region varied at half-precessional (similar to 11,500-year) intervals in phase with orbitally controlled insolation forcing. The southeasterly and northeasterly monsoons that advect moisture from the western Indian Ocean were strengthened in alternation when the inter-hemispheric insolation gradient was at a maximum; dry conditions prevailed when neither monsoon was intensified and modest local March or September insolation weakened the rain season that followed. On sub-millennial timescales, the temporal pattern of hydrological change on the East African Equator bears clear high-northern-latitude signatures, but on the orbital timescale it mainly responded to low-latitude insolation forcing. Predominance of low-latitude climate processes in this monsoon region can be attributed to the low-latitude position of its continental regions of surface air flow convergence, and its relative isolation from the Atlantic Ocean, where prominent meridional overturning circulation more tightly couples low-latitude climate regimes to high-latitude boundary conditions.

Hybridization of a class of "second order" models of traffic flow

Despite the simultaneous progress of traffic modelling both on the macroscopic and microscopic front, recent works [E. Bourrel, J.B. Lessort, Mixing micro and macro representation of traffic flow: a hybrid model based on the LWR theory, Transport. Res. Rec. 1852 (2003) 193–200; D. Helbing, M. Treiber, Critical discussion of “synchronized flow”, Coop. Transport. Dyn. 1 (2002) 2.1–2.24; A. Hennecke, M. Treiber, D. Helbing, Macroscopic simulations of open systems and micro–macro link, in: D. Helbing, H.J. Herrmann, M. Schreckenberg, D.E. Wolf (Eds.), Traffic and Granular Flow ’99, Springer, Berlin, 2000, pp. 383–388] highlighted that one of the most promising way to simulate efficiently traffic flow on large road networks is a clever combination of both traffic representations: the hybrid modelling. Our focus in this paper is to propose two hybrid models for which the macroscopic (resp. mesoscopic) part is based on a class of second order model [A. Aw, M. Rascle, Resurection of second order models of traffic flow?, SIAM J. Appl. Math. 60 (2000) 916–938] whereas the microscopic part is a Follow-the Leader type model [D.C. Gazis, R. Herman, R.W. Rothery, Nonlinear follow-the-leader models of traffic flow, Oper. Res. 9 (1961) 545–567; R. Herman, I. Prigogine, Kinetic Theory of Vehicular Traffic, American Elsevier, New York, 1971]. For the first hybrid model, we define precisely the translation of boundary conditions at interfaces and for the second one we explain the synchronization processes. Furthermore, through some numerical simulations we show that the waves propagation is not disturbed and the mass is accurately conserved when passing from one traffic representation to another.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Comment on “Electron transport through correlated molecules computed using the time-independent Wigner function: Two critical tests”

The many-electron-correlated scattering (MECS) approach to quantum electronic transport was investigated in the linear-response regime [I. Bâldea and H. Köppel, Phys. Rev. B 78, 115315 (2008). The authors suggest, based on numerical calculations, that the manner in which the method imposes boundary conditions is unable to reproduce the well-known phenomena of conductance quantization. We introduce an analytical model and demonstrate that conductance quantization is correctly obtained using open system boundary conditions within the MECS approach.

High pressures generated by laser driven shocks: applications to planetary physics

High power lasers are a tool that can be used to determine important parameters in the context of Warm Dense Matter, i.e. at the convergence of low-temperature plasma physics and finite-temperature condensed matter physics. Recent results concerning planet inner core materials such as water and iron are presented. We determined the equation of state, temperature and index of refraction of water for pressures up to 7 Mbar. The release state of iron in a LiF window allowed us to investigate the melting temperature near the inner core boundary conditions. Finally, the first application of proton radiography to the study of shocked material is also discussed.

Scaling of Superdomain Bands in Ferroelectric Dots

Bundles of 90° stripe domains have been observed to form into distinct groups, or bands, in mesoscale BaTiO3 single crystal dots. Vector piezoresponse force microscopy (PFM) shows that each band region, when considered as a single entity, possesses a resolved polarization that lies approximately along the pseudocubic direction; antiparallel alignment of this resultant polarization in adjacent bands means that these regions can be considered as 180° “superdomains.” For dots with sidewall dimensions below ~2 microns, Landau–Kittel like scaling in the width of these superdomains was observed, strongly suggesting that they form in response to lateral depolarizing fields. In larger dot structures, scaling laws break down. We have rationalized these observations by considering changes in the driving force for the adoption of equilibrium superdomain periodicities implied by Landau–Kittel-free energy models; we conclude that the formation of ordered bands of superdomains is a uniquely meso/nanoscale phenomenon. We also note that the superdomain bands found by PFM imaging in air contrast with the quadrant arrangements seen previously by Schilling et al. (Nano Lett., 9, 3359 (2009)) through transmission electron microscopy imaging in vacuum. The importance of the exact nature of the boundary conditions in determining the domain patterns that spontaneously form in nanostructures is therefore clearly implied.

Calculation of loss in conductors of a superconducting shielded slotted line and a coplanar waveguide

Finite conductivity in superconductors is taken into account by approximate boundary conditions imposed directly when deriving pair summatory equations, which are solved using the Galerkin method and the basis describing the edge singularity.

Comparison of Experimental PIV Data and CFD Simulations for Flow in a Diesel Particulate Filter Inlet Diffuser

Flow maldistribution of the exhaust gas entering a Diesel Particulate Filter (DPF) can cause uneven soot distribution during loading and excessive temperature gradients during the regeneration phase. Minimising the magnitude of this maldistribution is therefore an important consideration in the design of the inlet pipe and diffuser, particularly in situations where packaging constraints dictate bends in the inlet pipe close to the filter, or a sharp diffuser angle. This paper describes the use of Particle Image Velocimetry (PIV) to validate a Computational Fluid Dynamic (CFD) model of the flow within the inlet diffuser of a DPF so that CFD can be used with confidence as a tool to minimise this flow maldistribution. PIV is used to study the flow of gas into a DPF over a range of steady state flow conditions. The distribution of flow approaching the front face of the substrate was of particular interest to this study. Optically clear diffusing cones were designed and placed between pipe and substrate to allow PIV analysis to take place. Stereoscopic PIV was used to eliminate any error produced by the optical aberrations caused by looking through the curved wall of the inlet cone. In parallel to the experiments, numerical analysis was carried out using a CFD program with an incorporated DPF model. Boundary conditions for the CFD simulations were taken from the experimental data, allowing an experimental validation of the numerical results. The CFD model incorporated a DPF model, the cement layers seen in segmented filters and the intumescent matting that is commonly used to pack the filter into a metal casing. The mesh contained approximately 580,000 cells and used the realizable ?-e turbulence model. The CFD simulation predicted both pressure drop across the DPF and the velocity field within the cone and at the DPF face with reasonable accuracy, providing confidence in the use the CFD in future work to design new, more efficient cones.

Tree-code simulations of proton acceleration from laser-irradiated wire targets

Recent experiments using Terawatt lasers to accelerate protons deposited on thin wire targets are modeled with a new type of gridless plasma simulation code. In contrast to conventional mesh-based methods, this technique offers a unique capability in emulating the complex geometry and open-ended boundary conditions characteristic of contemporary experimental conditions. Comparisons of ion acceleration are made between the tree code and standard particle-in-cell simulations for a typical collisionless

Evolutionary dynamics of the spatial Prisoner's Dilemma with self-inhibition

In this paper we study the influence of interventions on self-interactions in a spatial Prisoner's Dilemma on a two-dimensional grid with periodic boundary conditions and synchronous updating of the dynamics. We investigate two different types of self-interaction modifications. The first type (FSIP) is deterministic, effecting each self-interaction of a player by a constant factor, whereas the second type (PSIP) performs a probabilistic interventions. Both types of interventions lead to a reduction of the payoff of the players and, hence, represent inhibiting effects. We find that a constant but moderate reduction of self-interactions has a very beneficial effect on the evolution of cooperators in the population, whereas probabilistic interventions on self-interactions are in general counter productive for the coexistence of the two different strategies. (C) 2011 Elsevier Inc. All rights reserved.

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of a VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions for stress resultant analysis.