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.

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.

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.

An investigation of the changes in the natural frequency of a Pile affected by Scour

Scour around bridge foundations is one of the leading causes of bridge failure. Up until recently, the monitoring of this phenomenon was primarily based around using underwater instrumentation to monitor the progression of scour holes as they develop around foundation systems. Vibration-based damage detection techniques have been used to detect damage in bridge beams. The application of these vibration based methods to the detection of scour has come to the fore in research in recent years. This paper examines the effect that scour has on the frequency response of a driven pile foundation system, similar to those used to support road and rail bridges. The effect of scour on the vibration characteristics of the pile is examined using laboratory and field testing. It is clear that there is a very clear reduction in the natural frequency of the pile as the severity of scour increases. It is shown that by combining state-of-the-art geotechnical techniques with relatively simple finite element modelling approaches, it is possible to accurately predict the natural frequency of the pile for a given scour depth. Therefore, the paper proposes a method that would allow the estimation of scour depth for a given observed pile frequency.

Modelling underwater acoustic noise as a tool for coastal management

Dissertação de Mestrado, Gestão da Água e da Costa, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2009

A time frequency approach to blind deconvolution in multipath underwater channels

Blind deconvolution is studied in the underwater acoustic channel context, by time-frequency (TF) processing. The acoustic propagation environment is modelled by ray tracing and mathematically described by a multipath propagation channel. Representation of the received signal by means of a signal-dependent TF distribution (radially Gaussian kernel distribution) allowed to visualize the resolved replicas of the emitted signal, while signi cantly attenuating the inherent interferences of classic quadratic TF distributions. The source signal instantaneous frequency estimation was the starting point for both source and channel estimation. Source signature estimation was performed by either TF inversion, based on the Wigner-Ville distribution of the received signal, or a subspace- -based method. The channel estimate was obtained either via a TF formulation of the conventional matched- lter, or via matched- - ltering with the previously obtained source estimate. A shallow water realistic scenario is considered, comprising a 135-m depth water column and an acoustic source located at 90-m depth and 5.6-km range from the receiver. For the corresponding noiseless simulated data, the quality of the best estimates was 0.856 for the source signal, and 0.9664 and 0.9996 for the amplitudes and time-delays of the impulse response, respectively. Application of the proposed deconvolution method to real data of the INTIMATE '96 sea trial conduced to source and channel estimates with the quality of 0.530 and 0.843, respectively. TF processing has proved to remove the typical ill-conditioning of single sensor deterministic deconvolution techniques.

Analysis of underwater acoustic data via 3-D segmentation

The analysis of seabed structure is important in a wide variety of scientific and industrial applications. In this paper, underwater acoustic data produced by bottom-penetrating sonar (Topas) are analyzed using unsupervised volumetric segmentation, based on a three dimensional Gibbs-Markov model. The result is a concise and accurate description of the seabed, in which key structures are emphasized. This description is also very well suited to further operations, such as the enhancement and automatic recognition of important structures. Experimental results demonstrating the effectiveness of this approach are shown, using Topas data gathered in the North Sea off Horten, Norway.

Estimating depth and explosive charge weight for an extremely shallow underwater explosion of the ROKS Cheonan sinking in the Yellow Sea

We estimated the detonation depth and net explosive weight for a very shallow underwater explosion using cutoff frequencies and spectral analysis. With detonation depth and a bubble pulse the net explosive weight for a shallow underwater explosion could simply be determined. The ray trace modeling confirms the detonation depth as a source of the hydroacoustic wave propagation in a shallow channel. We found cutoff frequencies of the reflection off the ocean bottom to be 8.5 Hz, 25 Hz, and 43 Hz while the cutoff frequency of the reflection off the free surface to be 45 Hz including 1.01 Hz for the bubble pulse, and also found the cutoff frequency of surface reflection to well fit the ray-trace modeling. We also attempted to corroborate our findings using a 3D bubble shape modeling and boundary element method. Our findings led us to the net explosive weight of the underwater explosion offshore of Baengnyeong-do for the ROKS Cheonan sinking to be approximately 136 kg TNT at a depth of about 8 m within an ocean depth of around 44 m. © 2015 Elsevier B.V.

Measuring underwater noise with high endurance surface and underwater autonomous vehicles

This paper describes the results of AcousticRobot'13 - a noise measurement campaign that took place off the Portuguese Coast in May 2013, using two high endurance autonomous vehicles capable of silent operation (an underwater glider and an autonmomous sailing vessel) equipped with hydrophones, and a moored hydrophone that served as reference. We show that the autonomous vehicles used can provide useful measurements of underwater noise, and describe the main advantages and shortcomings that became evident during the campaign.

Geostatistical Prediction of Ocean Outfall Plume Characteristics Based on an Autonomous Underwater Vehicle

Geostatistics has been successfully used to analyze and characterize the spatial variability of environmental properties. Besides giving estimated values at unsampled locations, it provides a measure of the accuracy of the estimate, which is a significant advantage over traditional methods used to assess pollution. In this work universal block kriging is novelty used to model and map the spatial distribution of salinity measurements gathered by an Autonomous Underwater Vehicle in a sea outfall monitoring campaign, with the aim of distinguishing the effluent plume from the receiving waters, characterizing its spatial variability in the vicinity of the discharge and estimating dilution. The results demonstrate that geostatistical methodology can provide good estimates of the dispersion of effluents that are very valuable in assessing the environmental impact and managing sea outfalls. Moreover, since accurate measurements of the plume’s dilution are rare, these studies might be very helpful in the future to validate dispersion models.