Prediction of Transonic LCOs using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic Lim
it Cycle Oscillations using high fidelity analysis. CFD based Harmonic Balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling HB methods to accurately capture Limit Cycle Oscillations (LCOs); this is achieved by defining a frequency updating procedure based on a coupled CFD/CSD Harmonic Balance formulation to find the LCO condition. A pitch/plunge aerofoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both LCO amplitude and frequency, while producing at least one order of magnitude reduction in computational time.

Prediction of limit cycle oscillations under uncertainty using a Harmonic Balance method

The Harmonic Balance method is an attractive solution for computing periodic responses and can be an alternative to time domain methods, at a reduced computational cost. The current paper investigates using a Harmonic Balance method for simulating limit cycle oscillations under uncertainty. The Harmonic Balance method is used in conjunction with a non-intrusive polynomial-chaos approach to propagate variability and is validated against Monte Carlo analysis. Results show the potential of the approach for a range of nonlinear dynamical systems, including a full wing configuration exhibiting supercritical and subcritical bifurcations, at a fraction of the cost of performing time domain simulations.

Prediction of Transonic Limit-Cycle Oscillations Using an Aeroelastic Harmonic Balance Method

This work proposes a novel approach to compute transonic limit-cycle oscillations using high-fidelity analysis. Computational-Fluid-Dynamics based harmonic balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling harmonic balance methods to accurately capture limit-cycle oscillations; this is achieved by defining a frequency-updating procedure based on a coupled computational-fluid-dynamics/computational-structural-dynamics harmonic balance formulation to find the limit-cycle oscillation condition. A pitch/plunge airfoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both limit-cycle oscillation amplitude and frequency while producing at least a one-order-of-magnitude reduction in computational time.

Modeling Frequency-Dependent Boundaries as Digital Impedance Filters in FDTD and K-DWM Room Acoustics Simulations

A new method for modeling-frequency-dependent boundaries in finite-difference time-domain (FDTD) and Kirchhoff variable digital waveguide mesh (K-DWM) room acoustics simulations is presented. The proposed approach allows the direct incorporation of a digital impedance filter (DIF) in the Multidimensional (2D or 3D) FDTD boundary model of a locally reacting surface. An explicit boundary update equation is obtained by carefully constructing a Suitable recursive formulation. The method is analyzed in terms of pressure wave reflectance for different wall impedance filters and angles of incidence. Results obtained from numerical experiments confirm the high accuracy of the proposed digital impedance filter boundary model, the reflectance of which matches locally reacting surface (LRS) theory closely. Furthermore a numerical boundary analysis (NBA) formula is provided as a technique for an analytic evaluation of the numerical reflectance of the proposed digital impedance filter boundary formulation.

A Phase Grating Approach to Modeling Surface Diffusion in FDTD Room Acoustics Simulations

In this paper, a method for modeling diffusive boundaries in finite difference time domain (FDTD) room acoustics simulations with the use of impedance filters is presented. The proposed technique is based on the concept of phase grating diffusers, and realized by designing boundary impedance filters from normal-incidence reflection filters with added delay. These added delays, that correspond to the diffuser well depths, are varied across the boundary surface, and implemented using Thiran allpass filters. The proposed method for simulating sound scattering is suitable for modeling high frequency diffusion caused by small variations in surface roughness and, more generally, diffusers characterized by narrow wells with infinitely thin separators. This concept is also applicable to other wave-based modeling techniques. The approach is validated by comparing numerical results for Schroeder diffusers to measured data. In addition, it is proposed that irregular surfaces are modeled by shaping them with Brownian noise, giving good control over the sound scattering properties of the simulated boundary through two parameters, namely the spectral density exponent and the maximum well depth.

Binaural Reproduction of Finite Difference Simulations Using Spherical Array Processing

Due to its efficiency and simplicity, the finite-difference time-domain method is becoming a popular choice for solving wideband, transient problems in various fields of acoustics. So far, the issue of extracting a binaural response from finite difference simulations has only been discussed in the context of embedding a listener geometry in the grid. In this paper, we propose and study a method for binaural response rendering based on a spatial decomposition of the sound field. The finite difference grid is locally sampled using a volumetric array of receivers, from which a plane wave density function is computed and integrated with free-field head related transfer functions, in the spherical harmonics domain. The volumetric array is studied in terms of numerical robustness and spatial aliasing. Analytic formulas that predict the performance of the array are developed, facilitating spatial resolution analysis and numerical binaural response analysis for a number of finite difference schemes. Particular emphasis is placed on the effects of numerical dispersion on array processing and on the resulting binaural responses. Our method is compared to a binaural simulation based on the image method. Results indicate good spatial and temporal agreement between the two methods.

Discrete-Time Conserved Quantities for Damped Oscillators

Numerical sound synthesis is often carried out using the finite difference time domain method. In order to analyse the stability of the derived models, energy methods can be used for both linear and nonlinear settings. For Hamiltonian systems the existence of a conserved numerical energy-like quantity can be used to guarantee the stability of the simulations. In this paper it is shown how to derive similar discrete conservation laws in cases where energy is dissipated due to friction or in the presence of an energy source due to an external force. A damped harmonic oscillator (for which an analytic solution is available) is used to present the proposed methodology. After showing how to arrive at a conserved quantity, the simulation of a nonlinear single reed shows an example of an application in the context of musical acoustics.

A developed stage of Alfven wave phase mixing

Alfven wave phase mixing is an extensively studied mechanism for dissipating wave energy in an inhomogeneous medium. It is common in the vast majority of phase mixing papers to assume that even though short scale lengths and steep gradients develop as a result of phase mixing, nonlinear wave coupling does not occur. However, weakly nonlinear studies have shown that phase mixing generates magnetoacoustic modes. Numerical results are presented which show the nonlinear generation of magnetosonic waves by Alfven wave phase mixing. The efficiency of the effect is determined by the wave amplitude, the frequency of the Alfven waves and the gradient in the background Alfven speed. Weakly nonlinear theory has shown that the amplitude of the fast magnetosonic wave grows linearly in time. The simulations presented in this paper extend this result to later times and show saturation of the fast magnetosonic component at amplitudes much lower than that of the Alfven wave. For the case when Alfven waves are driven at the boundary, simulating photospheric footpoint motion, a clear modulation of the saturated amplitude is observed. All the results in this paper are for a low amplitude (less than or equal to 0.1), single frequency Alfven wave and a uniform background magnetic field in a two dimensional domain. For this simplified geometry, and with a monochromatic driver, we concluded that the nonlinear generation of fast modes has little effect on classical phase mixing.

Reduced-complexity iterative equalization for severe time-dispersive MIMO channels

In this paper, a reduced-complexity soft-interference-cancellation minimum mean-square-error.(SIC-MMSE) iterative equalization method for severe time-dispersive multiple-input-multiple-output (MIMO) channels is proposed. To mitigate the severe time dispersiveness of the channel, a single carrier with cyclic prefix is employed, and the equalization is per-formed in the frequency domain. This simplifies the challenging problem of equalization in MIMO channels due to both the intersymbol interference (ISI) and the coantenna interference (CAI). The proposed iterative algorithm works in two stages. The first stage estimates the transmitted frequency-domain symbols using a low-complexity SIC-MMSE equalizer. The second stage converts the estimated frequency-domain symbols in the time domain and finds their means and variances to incorporate in the SIC-MMSE equalizer in the next iteration. Simulation results show the bit-/symbol-error-rate performance of the SIC-MMSE equalizer, with and without coding, for various modulation schemes.

Radiowave propagation from a tissue-implanted source at 418 MHz and 916.5 MHz

Tissue-implanted ultra-high frequency (UHF) radio devices are being employed in both humans and animals for telemetry and telecommand applications, This paper describes the experimental measurement and electromagnetic modeling of propagation from 418-MHz and 916.5-MHz sources placed in the human vagina. Whole-body homogeneous and semi-segmented software models were constructed using data from the Visible Human Project. Bodyworn radiation efficiencies for a vaginally placed 418-MHz source were calculated using finite-difference time-domain and ranged between 1.6% and 3.4% (corresponding to net body losses of between 14.7 and 18.0 dB), Greater losses were encountered at 916.5 MHz, with efficiencies between 0.36% and 0.46% (net body loss ranging between 23.4 and 24.4 dB), Practical measurements were in good agreement with simulations, to within 2 dB at 418 MHz and 3 dB at 916.5 MHz. The degree of tissue-segmentation for whole-body models was found to have a minimal effect on calculated azimuthal radiation patterns and bodyworn radiation efficiency, provided the region surrounding the implanted source was sufficiently detailed.

Quality factor assessment of subwavelength resonant cavities at FIR frequencies

Subwavelength resonators at FIR are presented and studied. The structures consist of 1D cavities formed between a metallized (silver) surface and a metamaterial surface comprising a periodic array of silver patches on a silver-backed silicon substrate. The concept derives from recent discoveries of artificial magnetic conductors (AMC). By studying the currents excited on the metamaterial surface by a normally incident plane wave, the nature of the emerging resonant phenomena and the physical mechanism underlying the AMC operation are investigated. Full wave simulations, based on finite element method and time-domain transmission line modelling technique, have been carried out to demonstrate the effective AMC boundary condition and prove the possibilities for subwavelength cavities. The quality factor of the resonant cavities is assessed as a function of the cavity profile. It is demonstrated that the quality factor drops to about 1/8 of the half-wavelength value for lambda/8 resonant cavity.

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.

Formulation of Locally Reacting Surfaces in FDTD/K-DWM Modelling of Acoustic Spaces

In this paper, we present new methods for constructing and analysing formulations of locally reacting surfaces that can be used in finite difference time domain (FDTD) simulations of acoustic spaces. Novel FDTD formulations of frequency-independent and simple frequency-dependent impedance boundaries are proposed for 2D and 3D acoustic systems, including a full treatment of corners and boundary edges. The proposed boundary formulations are designed for virtual acoustics applications using the standard leapfrog scheme based on a rectilinear grid, and apply to FDTD as well as Kirchhoff variable digital waveguide mesh (K-DWM) methods. In addition, new analytic evaluation methods that accurately predict the reflectance of numerical boundary formulations are proposed. numerical experiments and numerical boundary analysis (NBA) are analysed in time and frequency domains in terms of the pressure wave reflectance for different angles of incidence and various impedances. The results show that the proposed boundary formulations structurally adhere well to the theoretical reflectance. In particular, both reflectance magnitude and phase are closely approximated even at high angles of incidence and low impedances. Furthermore, excellent agreement was found between the numerical boundary analysis and the experimental results, validating both as tools for researching FDTD boundary formulations.

Time-resolved studies of ultrafast wavepacket dynamics in hydrogen molecules

Recent advances in the study of quantum vibrations and rotations in the fundamental hydrogen molecules are reported. Using the deuterium molecules (D-2(+) and D-2) as exemplars, the application of ultrafast femtosecond pump-probe experiments to study the creation and time-resolved imaging of coherent nuclear wavepackets is discussed. The ability to study the motion of these fundamental molecules in the time-domain is a notable milestone, made possible through the advent of ultrashort intense laser pulses with durations on sub-vibrational (and sub-rotational) timescales. Quantum wavepacket revivals are characterised for both vibrational and rotational degrees of freedom and quantum models are used to provide a detailed discussion of the underlying ultrafast physical dynamics for the specialist and non-specialist alike. (C) 2009 Elsevier B.V. All rights reserved.

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.