Conventional finite-difference or pseudo-spectral viscoacoustic modeling methods use uniform, rectangular grids to simulate wave propagation, leading to severe spurious diffractions at the surface topography. To overcome the problem, we mesh the velocity and quality factor (Q) models into curvilinear grids and map the models and the first-order velocity-stress viscoacoustic equation into an auxiliary

The radiated acoustic waves from impact pile driving produce high noise level into the water which may cause damage to marine mammals living close to the offshore construction location. In this paper, a linear, axisymmetric finite element (FE) model is applied to predict pile driving noise in the water. Measurement from bottom-mounted hydrophone deployed at a site 230 m from the source is used to validate

We derive a formula for the gradients of the total scattering cross-section (TSCS) with respect to positions of a set of cylindrical scatterers. The analytic form enhances modeling capability when combined with optimization algorithms and parallel computing. As application of the method, we consider a gradient-based minimization of TSCS for a set of cylindrical obstacles by incrementally repositioning

A method is presented for the closed-form evaluation of the singular and near-singular integrals arising in the Boundary Element Method solution of the Helmholtz equation. An error analysis is presented for the numerical evaluation of such integrals on a plane element, and used to develop a criterion for the selection of quadrature rules. The analytical approach is based on an optimized expansion of

Computational speeds are compared for modern implementations of three parabolic equation approximations. The split-step c0-insensitive model is 3.7 and 5.5 times faster than the finite-difference model OWWE and finite-element model RAM, respectively. Calculations are made between a source and receiver separated horizontally by 1000km at 600m depth near the minimum of sound speed in the deep ocean.

Propagation of Lamb waves in both laminated and functionally graded (FG) composites is analyzed. Anomalous dispersion at high frequencies is observed revealing substantial discrepancy in the high frequencies asymptotes for quasi-flexural and quasi-symmetric fundamental modes. The applied methodology is based on a variant of sextic complex formalism.

This paper collects the state of the art and the tremendous progress that has been made in hybrid modeling of aeroacoustic sound. Hybrid modeling is defined such that flow and acoustics are modeled separate and connected by an aeroacoustic model. The contributions will be classified with respect to the aeroacoustic models being developed, covering Lighthill’s analogy, Ffowcs Williams and Hawkings,

Crossing internal wave trains are commonly observed in continental shelf shallow water. In this paper, we study the effects of crossing internal wave structures on three-dimensional acoustic ducts with both theoretical and numerical approaches. We show that, depending on the crossing angle, acoustic energy, which is trapped laterally between internal waves of one train, can be either scattered, cross-ducted

Pile driving is used for constructing foundation supports for offshore structures. Underwater noise, induced by in-water pile driving, could adversely impact marine life near the piling location. Many studies have computed this noise in close ranges by using semi-analytical models and Finite Element Method (FEM) models. This work presents a Spectral Element Method (SEM) wave simulator as an alternative

The frequency-domain linearized Navier–Stokes equations (LNSEs) are used to describe the sound field of Helmholtz resonators and concentric perforated tube resonators in the presence of low Mach number flow. The numerical procedure of LNSEs method is performed in three steps, computational fluid dynamics (CFD) calculation, data transfer and acoustics calculation. The transmission loss predictions of

An acoustic metasurface (AM) is proposed to achieve a sound tunnel with switchable functionalities of asymmetric transmission and bidirectional sound isolation, which is composed of an array of steel containers filled with two inert gases, argon and xenon. Acoustic asymmetric transmission can be observed when the AM is placed vertically. The functionality of the tunnel can be switched to bidirectional

For various types of finite acoustic sources placed on an infinite cylindrical baffle, the pressure solution in cylindrical coordinates can be given by an infinite series of Inverse Fourier Integrals involving a singular quotient of Hankel functions. A hybrid method is introduced addressing these integrals’ singularity analytically and truncating their infinite integration range with predictable error

Kirchhoff approximations for the forward-scattering target strength of underwater objects are developed by combining Babinet’s principle and the Kirchhoff integral, where theoretical formulations and a numerical implementation are given in detail. The Kirchhoff approximation is found to be a high-frequency physical acoustic approximation. The forward-scattering target strength versus frequency and

The numerical determination of dispersion relations in periodic materials via the finite element method is a difficult task in most standard codes. Here, we propose a novel technique which allows the computation of these band structures from local elemental subroutines in contrast with existing methods which impose Bloch boundary conditions on the global arrays. The proposed local approach is thus

This paper presents a new method known as laser-assisted, modified Helmholtz Equation Least Squares method-based reconstruction to determine vibro-acoustic quantities of an arbitrarily shaped vibrating structure subject to a noncontact acoustic excitation. Unlike the traditional near-field acoustical holography (NAH) that relies on the acoustic pressures measured in the near-field, this new methodology

The acoustic propagation problem in the ocean is modeled via the wide angle parabolic equation with a Neumann to Dirichlet map bottom boundary condition. An environment consisting of the water column, a sediment layer and the semi-infinite sub-bottom region is considered. The derivatives of a new cost function with respect to the unknown environmental parameters are calculated analytically via the

In shallow-water waveguide, the normal mode information of acoustic signals plays an important role in localization research. In this paper, an extraction algorithm of modal intensity distribution based on horizontal circular array is proposed. By using the horizontal wavenumber of each mode, the modal domain beamforming is used to estimate the modal intensity distribution of the acoustic signal. The

Algorithms that enable acoustic sensing vehicles to autonomously map acoustic fields are being developed, but some require an analytical representation of the data collected from the acoustic field with specific continuity properties. One optimal control technique that has these requirements and has been used in autonomous acoustic sensing algorithms is Pontryagin’s Maximum Principle (PMP). Given that

A new spectrally optimized physical dispersion relation preserving scheme has been introduced to solve computational acoustics and aeroacoustics problems, accurately. The derived fourth-order accurate scheme has significant spectral resolution and physical dispersion relation preserving (DRP) nature. The scheme displays neutral stability at high CFL numbers. The developed scheme has an ability to add

Adiabatic approximation (AP) combined with perturbation theory gives a fast normal-mode solution of temporal coherence for sound field in a two-dimensional deep water with time-varying random internal waves. Internal waves induced mode changes are deduced using the first-order perturbation theory [C. T. Tindle, L. M. O’Driscoll and C. J. Higham, Coupled mode perturbation theory of range dependence

The bottom parameters of the deep ocean are difficult to obtain through in situ measurement. These parameters demonstrate a significant physical meaning for predicting sound field accurately. Thus, geoacoustic inversion is required. An acoustic experiment was performed on a reliable acoustic path (RAP) in the Philippine Sea in 2013. A single bottom-moored hydrophone was deployed as the receiver, and

The Kirchhoff approximation (K-A) to calculate the acoustic backscattering of a complex structure can be evaluated using a discretized version of its surface (i.e, a mesh). From the computational viewpoint, the most handy approach is the one based on flat facets. However, in the high-frequency range, where the K-A provides good agreement and is therefore applicable, it requires a mesh with such a large

We propose a fast 2D full-waveform forward and inverse scattering solver to reconstruct ultrasonic speed and attenuation properties. The forward method is based on the volume integral equation, accelerated by the extended Born approximation and biconjugate-gradient fast Fourier transform method. The inverse method is based on the extended contrast source inversion. The research results show that our

The source locations of axisymmetric modes of screech tones are numerically investigated. Fourth-order optimized compact scheme and fourth-order Runge–Kutta method are used to solve the 2-D axisymmetric Euler equations. The screech tone is successfully reproduced, and the change in wavelength with respect to jet Mach number shows good agreement with the experimental data. At various low supersonic

Acoustic scattering-field reconstruction of structures with arbitrary shape is the research basis of the scattering characteristics for underwater targets. Firstly, using boundary element method (BEM) and acoustic radiation modes (ARMs) solution in fluid domain, it is proved that the scattering pressures can be expressed by ARMs. Secondly, the acoustic field distribution modes (AFDMs) are constructed

Reverberation level (RL) is modeled in a shallow water environment with an underwater sound channel and internal solitary wave (ISW) using the coupled mode method. Numerical RL result based on one-way coupled mode shows an abnormal increase when a source is located near the channel axis and the ISW is located far from the source. The abnormal increase is analyzed by using a two-mode approach (assuming