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.

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.

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.

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.

GALEX AND PAN-STARRS1 DISCOVERY OF SN IIP 2010aq: THE FIRST FEW DAYS AFTER SHOCK BREAKOUT IN A RED SUPERGIANT STAR

We present the early UV and optical light curve of Type IIP supernova (SN) 2010aq at z = 0.0862, and compare it to analytical models for thermal emission following SN shock breakout in a red supergiant star. SN 2010aq was discovered in joint monitoring between the Galaxy Evolution Explorer (GALEX) Time Domain Survey (TDS) in the NUV and the Pan-STARRS1 Medium Deep Survey (PS1 MDS) in the g, r, i, and z bands. The GALEX and Pan-STARRS1 observations detect the SN less than 1 day after the shock breakout, measure a diluted blackbody temperature of 31,000 +/- 6000 K 1 day later, and follow the rise in the UV/optical light curve over the next 2 days caused by the expansion and cooling of the SN ejecta. The high signal-to-noise ratio of the simultaneous UV and optical photometry allows us to fit for a progenitor star radius of 700 +/- 200R(circle dot), the size of a red supergiant star. An excess in UV emission two weeks after shock breakout compared with SNe well fitted by model atmosphere-code synthetic spectra with solar metallicity is best explained by suppressed line blanketing due to a lower metallicity progenitor star in SN 2010aq. Continued monitoring of PS1 MDS fields by the GALEX TDS will increase the sample of early UV detections of Type II SNe by an order of magnitude and probe the diversity of SN progenitor star properties.

Low-complexity iterative method of equalization for single carrier with cyclic prefix in doubly selective channels

Orthogonal frequency division multiplexing (OFDM) requires an expensive linear amplifier at the transmitter due to its high peak-to-average power ratio (PAPR). Single carrier with cyclic prefix (SC-CP) is a closely related transmission scheme that possesses most of the benefits of OFDM but does not have the PAPR problem. Although in a multipath environment, SC-CP is very robust to frequency-selective fading, it is sensitive to the time-selective fading characteristics of the wireless channel that disturbs the orthogonality of the channel matrix (CM) and increases the computational complexity of the receiver. In this paper, we propose a time-domain low-complexity iterative algorithm to compensate for the effects of time selectivity of the channel that exploits the sparsity present in the channel convolution matrix. Simulation results show the superior performance of the proposed algorithm over the standard linear minimum mean-square error (L-MMSE) equalizer for SC-CP.

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.

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.

Wavelet analysis of bender element signals

Accurate determination of shear wave arrival time using bender elements may be severely affected by a combination of near-field effects and reflected waves. These may mask the first arrival of the shear wave, making its detection difficult in the time domain. This paper describes an approach for measuring the shear wave arrival time through analysis of the output signal in the time-scale domain using a multi-scale wavelet transform. The local maxima lines of the wavelet transform modulus are observed at different scales, and all singularities are mathematically characterised, allowing subsequent detection of the singularity relating to the first arrival. Examples of the use of this approach on typical synthetic and experimental bender element signals are also supplied, and these results are compared with those from other time and frequency domain approaches. The wavelet approach is not affected by near-field effects, but instead is characterised by a relatively consistent singularity related to the shear wave arrival time, across a range of frequencies and noise levels.

A Physically-Constrained Source Model For FDTD Acoustic Simulation

The Finite Difference Time Domain (FDTD) method is becoming increasingly popular for room acoustics simulation. Yet, the literature on grid excitation methods is relatively sparse, and source functions are traditionally implemented in a hard or additive form
using arbitrarily-shaped functions which do not necessarily obey the physical laws of sound generation. In this paper we formulate
a source function based on a small pulsating sphere model. A physically plausible method to inject a source signal into the grid
is derived from first principles, resulting in a source with a near-flat spectrum that does not scatter incoming waves. In the final
discrete-time formulation, the source signal is the result of passing a Gaussian pulse through a digital filter simulating the dynamics of the pulsating sphere, hence facilitating a physically correct means to design source functions that generate a prescribed sound field.

Tip-enhanced Raman microscopy: practicalities and limitations

The feasibility of apertureless scanning near-field Raman microscopy, exploiting the local enhancement in Raman scattering in the vicinity of a silver or gold tip, was investigated. Using the finite difference time domain method we calculated the enhancement of electric field strength, and hence Raman scattering, achieved through the resonant excitation of local modes in the tip. By modelling the frequency-dependent dielectric response of the metal tip we were able to highlight the resonant nature of the tip-enhancement and determine the excitation wavelength required for the strongest electric field enhancement, and hence Raman scattering intensity, which occurs for the excitation of modes localized at the tip apex. It is demonstrated that a peak Raman enhancement of 10(7)-fold should be achievable with <5 nm spatial resolution. We show that surface-enhanced Raman scattering from carbon contamination on a silver or gold tip can be significant. However, we find for a tip of radius of curvature 20 nm that the Raman enhancement should decay totally within 20 nm from the tip. Hence withdrawal of the tip by this distance should lead to the disappearance of the tip-enhanced signal, leaving only that from carbon contamination on the tip itself and the intrinsic signal from the sample. Copyright (C) 2003 John Wiley Sons, Ltd.

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.

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.

Determination of Non-linear Damping Coefficients of bottom-hinged Oscillating Wave Surge Converters Using Numerical Free Decay Tests

Linear wave theory models are commonly applied to predict the performance of bottom-hinged oscillating wave surge converters (OWSC) in operational sea states. To account for non-linear effects, the additional input of coefficients not included in the model itself becomes necessary. In ocean engineering it is
common practice to obtain damping coefficients of floating structures from free decay tests. This paper presents results obtained from experimental tank tests and numerical computational fluid dynamics simulations of OWSC’s. Agreement between numerical and experimental methods is found to be very good, with CFD providing more data points at small amplitude rotations.
Analysis of the obtained data reveals that linear quadratic-damping, as commonly used in time domain models, is not able to accurately model the occurring damping over the whole regime of rotation amplitudes. The authors
conclude that a hyperbolic function is most suitable to express the instantaneous damping ratio over the rotation amplitude and would be the best choice to be used in coefficient based time domain models.

Optical coherence tomography for the diagnosis of neovascular age-related macular degeneration: a systematic review

The purpose is to study the diagnostic performance of optical coherence tomography (OCT) and alternative diagnostic tests for neovascular age-related macular degeneration (nAMD). Methods employed are as follows:systematic review and meta-analysis; Index test: OCT including time-domain (TD-OCT) and the most recently developed spectral domain (SD-OCT); comparator tests: visual acuity, clinical evaluation (slit lamp), Amsler chart, colour fundus photographs, infra-red reflectance, red-free images/blue reflectance, fundus autofluorescence imaging (FAF), indocyanine green angiography (ICGA), preferential hyperacuity perimetry (PHP), and microperimetry; reference standard: fundus fluorescein angiography. Databases searched included MEDLINE, MEDLINE In Process, EMBASE, Biosis, SCI, the Cochrane Library, DARE, MEDION, and HTA database. Last literature searches: March 2013. Risk of bias assessed using QUADAS-2. Meta-analysis models were fitted using hierarchical summary receiver operating characteristic (HSROC) curves. Twenty-two studies (2 abstracts and 20 articles) enrolling 2124 participants were identified, reporting TD-OCT (12 studies), SD-OCT (1 study), ICGA (8 studies), PHP (3 studies), Amsler grid, colour fundus photography and FAF (1 study each). Most studies were considered to have a high risk of bias in the patient selection (55%, 11/20), and flow and timing (40%, 8/20) domains. In a meta-analysis of TD-OCT studies, sensitivity and specificity (95% CI) were 88% (46–98%) and 78% (64–88%), respectively. There was insufficient information to undertake meta-analysis for other tests. TD-OCT is a sensitive test for detecting nAMD, although specificity was only moderate. Data on SD-OCT are sparse. Diagnosis of nAMD should not rely solely on OCT.