An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.

Compact Routing Messages in Self-Healing Trees

Existing compact routing schemes, e.g., Thorup and Zwick [SPAA 2001] and Chechik [PODC 2013], often have no means to tolerate failures, once the system has been setup and started. This paper presents, to our knowledge, the first self-healing compact routing scheme. Besides, our schemes are developed for low memory nodes, i.e., nodes need only O(log2 n) memory, and are thus, compact schemes.
We introduce two algorithms of independent interest: The first is CompactFT, a novel compact version (using only O(log n) local memory) of the self-healing algorithm Forgiving Tree of Hayes et al. [PODC 2008]. The second algorithm (CompactFTZ) combines CompactFT with Thorup-Zwick’s treebased compact routing scheme [SPAA 2001] to produce a fully compact self-healing routing scheme. In the self-healing model, the adversary deletes nodes one at a time with the affected nodes self-healing locally by adding few edges. CompactFT recovers from each attack in only O(1) time and ∆ messages, with only +3 degree increase and O(log∆) graph diameter increase, over any sequence of deletions (∆ is the initial maximum degree).
Additionally, CompactFTZ guarantees delivery of a packet sent from sender s as long as the receiver has not been deleted, with only an additional O(y log ∆) latency, where y is the number of nodes that have been deleted on the path between s and t. If t has been deleted, s gets informed and the packet removed from the network.

Numerical simulations of fluid-structure interactions in single-reed mouthpieces

Most single-reed woodwind instrument models rely on a quasistationary approximation to describe the relationship between the volume flow and. the pressure difference across the reed channel. Semiempirical models based on the quasistationary approximation are very useful in explaining the fundamental characteristics of this family of instruments such as self-sustained oscillations and threshold of blowing pressure. However, they fail at explaining more complex phenomena associated with the fluid-structure interaction during dynamic flow regimes, such as the transient and steady-state behavior of the system as a function. of the mouthpiece geometry. Previous studies have discussed the accuracy of the quasistationary approximation but the amount of literature on the subject is sparse, mainly due to the difficulties involved in the measurement of dynamic flows in channels with an oscillating reed. In this paper, a numerical technique based on the lattice Boltzmann method and a finite difference scheme is proposed in order to investigate the characteristics of fully coupled fluid-structure interaction in single-reed mouthpieces with different channel configurations. Results obtained for a stationary simulation with a static reed agree very well with those predicted by the literature based on the quasistationary approximation. However, simulations carried out for a dynamic regime with dn oscillating reed show that the phenomenon associated with flow detachment and reattachment diverges considerably frorn the theoretical assumptions. Furthermore, in the case of long reed channels, the results obtained for the vena contracta factor are in significant disagreement with those predicted by theory. For short channels, the assumption of constant vena contracta was found to be valid for only 40% of the duty cycle. (c) 2007 Acoustical Society of America.

Room acoustics simulation using 3-D compact explicit FDTD schemes

This paper presents methods for simulating room acoustics using the finite-difference time-domain (FDTD) technique, focusing on boundary and medium modeling. A family of nonstaggered 3-D compact explicit FDTD schemes is analyzed in terms of stability, accuracy, and computational efficiency, and the most accurate and isotropic schemes based on a rectilinear grid are identified. A frequency-dependent boundary model that is consistent with locally reacting surface theory is also presented, in which the wall impedance is represented with a digital filter. For boundaries, accuracy in numerical reflection is analyzed and a stability proof is provided. The results indicate that the proposed 3-D interpolated wideband and isotropic schemes outperform directly related techniques based on Yee's staggered grid and standard digital waveguide mesh, and that the boundary formulations generally have properties that are similar to that of the basic scheme used.

On the numerical solution of the 2d wave equation with compact fdtd schemes

This paper discusses compact-stencil finite difference time domain (FDTD) schemes for approximating the 2D wave equation in the context of digital audio. Stability, accuracy, and efficiency are investigated and new ways of viewing and interpreting the results are discussed. It is shown that if a tight accuracy constraint is applied, implicit schemes outperform explicit schemes. The paper also discusses the relevance to digital waveguide mesh modelling, and highlights the optimally efficient explicit scheme.

On the derived category of a regular toric scheme

Let X be a quasi-compact scheme, equipped with an open covering by affine schemes U s = Spec A s . A quasi-coherent sheaf on X gives rise, by taking sections over the U s , to a diagram of modules over the coordinate rings A s , indexed by the intersection poset S of the covering. If X is a regular toric scheme over an arbitrary commutative ring, we prove that the unbounded derived category of quasi-coherent sheaves on X can be obtained from a category of Sop-diagrams of chain complexes of modules by inverting maps which induce homology isomorphisms on hyper-derived inverse limits. Moreover, we show that there is a finite set of weak generators, one for each cone in the fan S. The approach taken uses the machinery of Bousfield–Hirschhorn colocalisation of model categories. The first step is to characterise colocal objects; these turn out to be homotopy sheaves in the sense that chain complexes over different open sets U s agree on intersections up to quasi-isomorphism. In a second step it is shown that the homotopy category of homotopy sheaves is equivalent to the derived category of X.

Wideband and Isotropic Room Acoustics Simulation Using 2-D Interpolated FDTD Schemes

In this paper, a complete method for finite-difference time-domain modeling of rooms in 2-D using compact explicit schemes is presented. A family of interpolated schemes using a rectilinear, nonstaggered grid is reviewed, and the most accurate and isotropic schemes are identified. Frequency-dependent boundaries are modeled using a digital impedance filter formulation that is consistent with locally reacting surface theory. A structurally stable and efficient boundary formulation is constructed by carefully combining the boundary condition with the interpolated scheme. An analytic prediction formula for the effective numerical reflectance is given, and a stability proof provided. The results indicate that the identified accurate and isotropic schemes are also very accurate in terms of numerical boundary reflectance, and outperform directly related methods such as Yee's scheme and the standard digital waveguide mesh. In addition, one particular scheme-referred to here as the interpolated wideband scheme-is suggested as the best scheme for most applications.

Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

A comparison of nonstaggered compact FDTD schemes for the 3D wave equation

This paper aims at providing a better insight into the 3D approximations of the wave equation using compact finite-difference time-domain (FDTD) schemes in the context of room acoustic simulations. A general family of 3D compact explicit and implicit schemes based on a nonstaggered rectilinear grid is analyzed in terms of stability, numerical error, and accuracy. Various special cases are compared and the most accurate explicit and implicit schemes are identified. Further considerations presented in the paper include the direct relationship with other numerical approaches found in the literature on room acoustic modeling such as the 3D digital waveguide mesh and Yee's staggered grid technique.

Generation of a train of ultrashort pulses from a compact birefringent crystal array

A linear array of n calcite crystals is shown to allow the generation of a high contrast (> 10: 1) train of 2(n) high energy (> 100 mu J) pulses from a single ultrafast laser pulse. Advantage is taken of the pulse-splitting properties of a single birefringent crystal, where an incident laser pulse can be split into two pulses with orthogonal polarizations and equal intensity, separated temporally in proportion to the thickness of the crystal traversed and the difference in refractive indices of the two optic axes. In the work presented here an array of seven calcite crystals of sequentially doubled thickness is used to produce a train of 128 pulses, each of femtosecond duration. Readily versatile properties such as the number of pulses in the train and variable mark-space ratio are realized from such a setup. (c) 2007 Optical Society of America

Physical and numerical constraints in source modeling for finite difference simulation of room acoustics

In finite difference time domain simulation of room acoustics, source functions are subject to various constraints. These depend on the way sources are injected into the grid and on the chosen parameters of the numerical scheme being used. This paper addresses the issue of selecting and designing sources for finite difference simulation, by first reviewing associated aims and constraints, and evaluating existing source models against these criteria. The process of exciting a model is generalized by introducing a system of three cascaded filters, respectively, characterizing the driving pulse, the source mechanics, and the injection of the resulting source function into the grid. It is shown that hard, soft, and transparent sources can be seen as special cases within this unified approach. Starting from the mechanics of a small pulsating sphere, a parametric source model is formulated by specifying suitable filters. This physically constrained source model is numerically consistent, does not scatter incoming waves, and is free from zero- and low-frequency artifacts. Simulation results are employed for comparison with existing source formulations in terms of meeting the spectral and temporal requirements on the outward propagating wave.

Conical Horn Antenna with Spiral Phase Plate for Difference Pattern Generation

This paper describes the performance characteristics and experimental validation of a compact conical horn antenna with a dielectric cylinder spiral phase plate attached at its aperture. This performs the function of a spatial phase imprinting device creating a helical wave-front which results in a null in the far field radiation pattern of the antenna assembly.

Screened Hybrid Exact Exchange Correction Scheme for Adsorption Energies on Perovskite Oxides

The bond formation between an oxide surface and oxygen, which is of importance for numerous surface reactions including catalytic reactions, is investigated within the framework of hybrid density functional theory that includes nonlocal Fock exchange. We show that there exists a linear correlation between the adsorption energies of oxygen on LaMO3 (M = Sc–Cu) surfaces obtained using a hybrid functional (e.g., Heyd–Scuseria–Ernzerhof) and those obtained using a semilocal density functional (e.g., Perdew–Burke–Ernzerhof) through the magnetic properties of the bulk phase as determined with a hybrid functional. The energetics of the spin-polarized surfaces follows the same trend as corresponding bulk systems, which can be treated at a much lower computational cost. The difference in adsorption energy due to magnetism is linearly correlated to the magnetization energy of bulk, that is, the energy difference between the spin-polarized and the non-spin-polarized solutions. Hence, one can estimate the correction ...