Finite difference time domain modeling of phase grating diffusion

In this paper, a method for modeling diffusion caused by non-smooth boundary surfaces in simulations of room acoustics using finite difference time domain (FDTD) technique is investigated. The proposed approach adopts the well-known theory of phase grating diffusers to efficiently model sound scattering from rough surfaces. The variation of diffuser well-depths is attained by nesting allpass filters within the reflection filters from which the digital impedance filters used in the boundary implementation are obtained. The presented technique is appropriate for modeling diffusion at high frequencies caused by small surface roughness and generally diffusers that have narrow wells and infinitely thin separators. The diffusion coefficient was measured with numerical experiments for a range of fractional Brownian diffusers.

A Spherical Array Approach for Simulation of Binaural Impulse Responses using the Finite Difference Time Domain Method

The finite difference time domain (FDTD) method has direct applications in musical instrument modeling, simulation of environmental acoustics, room acoustics and sound reproduction paradigms, all of which benefit from auralization. However, rendering binaural impulse responses from simulated
data is not straightforward to accomplish as the calculated pressure at FDTD grid nodes does not contain any directional information. This paper addresses this issue by introducing a spherical array to capture sound pressure on a finite difference grid, and decomposing it into a plane-wave density
function. Binaural impulse responses are then constructed in the spherical harmonics domain by combining the decomposed grid data with free field head-related transfer functions. The effects of designing a spherical array in a Cartesian grid are studied, and emphasis is given to the relationships
between array sampling and the spatial and spectral design parameters of several finite-difference
schemes.

Comparison of Wavelet Transform and Time-Domain Analysis of Second Trimester Uterine Artery Doppler Waveforms in Screening for Pre-Eclampsia

The aim of this study was to compare time-domain waveform analysis of second-trimester uterine artery Doppler using the resistance index (RI) with waveform analysis using a mathematical tool known as wavelet transform for the prediction of pre-eclampsia (PE). This was a retrospective, nested case-cohort study of 336 women, 37 of whom subsequently developed PE. Uterine artery Doppler waveforms were analysed using both RI and waveform analysis. The utility of these indices in screening for PE was then evaluated using receiver operating characteristic curves. There were significant differences in uterine artery RI between the PE women and those with normal pregnancy outcome. After wavelet analysis, significant difference in the mean amplitude in wavelet frequency band 4 was noted between the 2 groups. The sensitivity for both Doppler RI and frequency band 4 for the detection of PE at a 10% false-positive rate was 45%. This small study demonstrates the application of wavelet transform analysis of uterine artery Doppler waveforms in screening for PE. Further prospective studies are needed in order to clearly define if this analytical approach to waveform analysis may have the potential to improve the detection of PE by uterine artery Doppler screening.

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.

Crack detection of steel girders using Brillouin optical time domain analysis

This article presents the results from an experimental program designed to evaluate the performance of a system consisting of a readout unit and a ribbon type Fiber Optic Sensor (FOS) based on Brillouin Optical Time Domain Analysis (BOTDA). The system is intended for the detection of cracks as well as the monitoring of long-term performance for steel bridge girders. The program consisted of introducing a crack at the center of a 3-m-long steel beam and monitoring its progression using static loading tests performed at ambient and sub-zero temperatures. For sensor lengths similar to those used in the field, the resonant frequency shifts per unit increase in crack width were found to decrease from 114 MHz/mm at ambient temperature (~25C) to 65 MHz/mm at -10C. Results also revealed nonlinearity and variability, which can be attributed to an incompatibility between the settings of the laser pump in the readout unit and the sensor length. Significant losses were detected along the bonded segments of the sensor and were attributed to the presence of ripples along the sensor. These undulations worsen with a reduction in temperature and are induced by the bonding procedure as well as the slack provided in the plastic sleeves containing the splices.

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.

Numerical analysis of bodyworn UHF antenna systems

Bodyworn antennas are found in a wide range of medical, military and personal communication applications, yet reliable communication from the surface of the human body still presents a range of engineering challenges. At UHF and microwave frequencies, bodyworn antennas can suffer from reduced efficiency due to electromagnetic absorption in tissue, radiation pattern fragmentation and variations in feed-point impedance. The significance and nature of these effects are system specific and depend on the operating frequency, propagation environment and physical constraints on the antenna itself. This paper describes how numerical electromagnetic modelling techniques such as FDTD (finite-difference time-domain) can be used in the design of bodyworn antennas. Examples are presented for 418 MHz, 916 .5 MHz and 2 . 45 GHz, in the context of both biomedical signalling and wireless personal-area networking applications such as the Bluetooth(TM)* wireless technology.

Periodic FDTD analysis of a 2-D leaky-wave planar antenna based on dipole frequency selective surfaces

A periodic finite-difference time-domain (FDTD) analysis is presented and applied for the first time in the study of a two-dimensional (2-D) leaky-wave planar antenna based on dipole frequency selective surfaces (FSSs). First, the effect of certain aspects of the FDTD modeling in the modal analysis of complex waves is studied in detail. Then, the FDTD model is used for the dispersion analysis of the antenna of interest. The calculated values of the leaky-wave attenuation constants suggest that, for an antenna of this type and moderate length, a significant amount of power reaches the edges of the antenna, and thus diffraction can play an important role. To test the validity of our dispersion analysis, measured radiation patterns of a fabricated prototype are presented and compared with those predicted by a leaky-wave approach based on the periodic FDTD results.

FDTD analysis of close-coupled 418 MHz radiating devices for human biotelemetry

This paper describes the finite-difference time-domain (FDTD) analysis of antenna-body interaction effects occurring when chest-mounted 418 MHz radio transmitters are used for medical telemetry applications. Whole-body software models (homogeneous, layered and tissue-segmented) were developed for an adult male subject. Using an electrically small (300 mm(2)) planar loop antenna, calculated radiation efficiencies ranged between 33.5% and 39.2% for a whole-body model, and between 60.7% and 66.1% for a torso; radiation patterns were found to be largely independent of model composition. The computed radiation efficiency for a 21.5 kg phantom representing a six-year-old female was within 1.1 dB of measured results (actual body mass 28 kg) and well-correlated azimuthal radiation patterns were noted.

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.

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.

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.

Manipulating plasmon modes in tip-enhanced Raman microscopy

The influence of polarization on the plasmon modes excited in tip-enhanced near-field optical microscopy have been investigated using the Finite Difference Time Domain Method. Analysis of the calculated results have laid particular emphasis on the ability to align local field enhancements with the orientation of molecules in order to optimize Raman signals, with particular relevance to recent experimental work on carbon nanotubes.