Internal absorbing boundary conditions for closed-aperture wavefield decomposition in solid media with unknown interiors.

We present a method to create an internal numerical absorbing boundary within elastic solid media whose properties are largely unknown and use it to create the first wavefield separation method that retrieves all orders of outgoing elastic wavefield constituents for real data recorded on a closed free surface. The recorded data are injected into a numerical finite-difference (FD) simulation along a closed, transparent surface, and the new internal numerical absorbing boundary condition achieves high attenuation of the ingoing waves radiated from the injection surface. This internal wave absorption enables the data injection to radiate all outgoing waves for experimental domains that include arbitrary unknown scatterers in the interior. The injection-absorption-based separation scheme is validated using three-dimensional (3D) synthetic modeling and a real data experiment acquired using a 3D laser Doppler vibrometer on a granite rock. The wavefield separation method forms a key component of an elastic immersive wave experimentation laboratory, and the ability to numerically absorb ingoing scattered energy in an uncharacterized medium while still radiating the true outgoing energy is intriguing and may lead to other development and applications in the future.

Broadband acoustic invisibility and illusions.

Rendering objects invisible to impinging acoustic waves (cloaking) and creating acoustic illusions (holography) has been attempted using active and passive approaches. While most passive methods are inflexible and applicable only to narrow frequency bands, active approaches attempt to respond dynamically, interfering with broadband incident or scattered wavefields by emitting secondary waves. Without prior knowledge of the primary wavefield, the signals for the secondary sources need to be estimated and adapted in real time. This has thus far impeded active cloaking and holography for broadband wavefields. We present experimental results of active acoustic cloaking and holography without prior knowledge of the wavefield so that objects remain invisible and illusions intact even for broadband moving sources. This opens previously inaccessible research directions and facilitates practical applications including architectural acoustics, education, and stealth.

Closed-aperture unbounded acoustics experimentation using multidimensional deconvolution.

In physical acoustic laboratories, wave propagation experiments often suffer from unwanted reflections at the boundaries of the experimental setup. We propose using multidimensional deconvolution (MDD) to post-process recorded experimental data such that the scattering imprint related to the domain boundary is completely removed and only the Green's functions associated with a scattering object of interest are obtained. The application of the MDD method requires in/out wavefield separation of data recorded along a closed surface surrounding the object of interest, and we propose a decomposition method to separate such data for arbitrary curved surfaces. The MDD results consist of the Green's functions between any pair of points on the closed recording surface, fully sampling the scattered field. We apply the MDD algorithm to post-process laboratory data acquired in a two-dimensional acoustic waveguide to characterize the wavefield scattering related to a rigid steel block while removing the scattering imprint of the domain boundary. The experimental results are validated with synthetic simulations, corroborating that MDD is an effective and general method to obtain the experimentally desired Green's functions for arbitrary inhomogeneous scatterers.

Seismological Processing of Six Degree-of-Freedom Ground-Motion Data.

Recent progress in rotational sensor technology has made it possible to directly measure rotational ground-motion induced by seismic waves. When combined with conventional inertial seismometer recordings, the new sensors allow one to locally observe six degrees of freedom (6DOF) of ground-motion, composed of three orthogonal components of translational motion and three orthogonal components of rotational motion. The applications of such 6DOF measurements are manifold-ranging from wavefield characterization, separation, and reconstruction to the reduction of non-uniqueness in seismic inverse problems-and have the potential to revolutionize the way seismic data are acquired and processed. However, the seismological community has yet to embrace rotational ground-motion as a new observable. The aim of this paper is to give a high-level introduction into the field of 6DOF seismology using illustrative examples and to summarize recent progress made in this relatively young field. It is intended for readers with a general background in seismology. In order to illustrate the seismological value of rotational ground-motion data, we provide the first-ever 6DOF processing example of a teleseismic earthquake recorded on a multicomponent ring laser observatory and demonstrate how wave parameters (phase velocity, propagation direction, and ellipticity angle) and wave types of multiple phases can be automatically estimated using single-station 6DOF processing tools. Python codes to reproduce this processing example are provided in an accompanying Jupyter notebook.

Identifying multiply scattered wavepaths in strongly scattering and dispersive media.

The ability to extract information from scattered waves is usually limited to singly scattered energy even if multiple scattering might occur in the medium. As a result, the information in arrival times of higher-order scattered events is underexplored. This information is extracted using fingerprinting theory. This theory has never previously been applied successfully to real measurements, particularly when the medium is dispersive. The theory is used to estimate the arrival times and scattering paths of multiply scattered waves in a thin sheet using an automated scheme in a dispersive medium by applying an additional dispersion compensation method. Estimated times and paths are compared with predictions based on a sequence of straight ray paths for each scattering event given the known scatterer locations. Additionally, numerical modelling is performed to verify the interpretations of the compensated data. Since the source also acts as a scatterer in these experiments, initially, the predictions and the numerical results did not conform to the experimental observations. By reformulating the theory and the processing scheme and adding a source scatterer in the modelling, it is shown that predictions of all observed scattering events are possible with both prediction methods, verifying that the methods are both effective and practically achievable.

Compensating for source directivity in immersive wave experimentation.

A physical boundary mounted with active sources can cancel acoustic waves arriving at the boundary, and emit synthesized waves into the neighboring medium to fully control the acoustic wavefield in an experimental setup such as a water tank or air-filled cavity. Using the same principles, a physical experiment can be artificially immersed within an extended virtual (numerical) domain so that waves propagate seamlessly between the experimental setup and virtual domain. Such an immersive wave control experiment requires physical monopolar sources on the active boundary. However, real physical sources (e.g., piezoelectric transducers) project waves at middle-to-high sonic frequencies (e.g., 1-20 kHz) that do not fully conform to the theoretically required monopolar radiation pattern; if left uncorrected, this causes controlled wavefields to deviate from those desired in immersive experiments. A method is proposed to compensate for the non-monopole-like radiation patterns of the sources, and can be interpreted physically in terms of Huygens principle. The method is implemented as a pre-computation procedure that modifies the extrapolation Green's functions in the Kirchhoff-Helmholtz integral before the actual experiments take place. Two-dimensional finite-difference simulations show that the processing method can effectively suppress the undesired effect caused by non-monopolar active sources in immersive wave control experiments.

Characterization of interface reflection coefficients using a finite-difference injection technique.

In this paper, a numerical wave field injection technique for characterizing the reflection coefficient of a planar medium interface is proposed. By injecting recorded wave field quantities into a three-dimensional (3D) finite-difference calculation, two key objectives are addressed: first, the recorded wave field is separated into its incident and reflected constituents without the need of spatial Fourier transforms or a temporal separation of incident and reflected parts in the recorded data. Second, the separated constituents are independently extrapolated to the location of the reflecting interface to determine its reflecting properties. The methodology is experimentally validated on 3D laboratory data consisting of reflections from the water-air interface in a water tank and is shown to give accurate results for incidence angles of up to 60°.

Broadband cloaking and holography with exact boundary conditions.

Broadband cloaking and holography are achieved by creating an exact boundary condition on a surface enclosing an object or free space. A time-recursive, discrete version of the Kirchhoff-Helmholtz integral predicts the wavefield impinging on the surface, as well as its transmission through an arbitrary embedding or replacement medium. Surface source distributions proportional to the predicted wavefield cancel the incident waves and radiate the desired response. The fields inside and outside the surface can be controlled independently. A two-dimensional numerical example shows that cloaking and holography can be achieved to within numerical precision across the frequency range of the incident radiation.

Immersive experimentation in a wave propagation laboratory.

A wave propagation laboratory is proposed which enables the study of the interaction of broadband signals with complex materials. A physical experiment is dynamically linked to a numerical simulation in real time through transmitting and recording transducer surfaces surrounding the target. The numerical simulation represents an arbitrarily larger domain, allowing experiments to be performed in a total environment much greater than the laboratory experiment itself. Specific applications include the study of non-linear effects or wave propagation in media where the physics of wave propagation is not well understood such as the effect of fine scale heterogeneity on broadband propagating waves.

Exact wave field simulation for finite-volume scattering problems.

An exact boundary condition is presented for scattering problems involving spatially limited perturbations of arbitrary magnitude to a background model in generally inhomogeneous acoustic media. The boundary condition decouples the wave propagation on a perturbed domain while maintaining all interactions with the background model, thus eliminating the need to regenerate the wave field response on the full model. The method, which is explicit, relies on a Kirchhoff-type integral extrapolation to update the boundary condition at every time step of the simulation. The Green's functions required for extrapolation through the background model are computed efficiently using wave field interferometry.

Modeling of wave propagation in inhomogeneous media.

We present a methodology providing a new perspective on modeling and inversion of wave propagation satisfying time-reversal invariance and reciprocity in generally inhomogeneous media. The approach relies on a representation theorem of the wave equation to express the Green function between points in the interior as an integral over the response in those points due to sources on a surface surrounding the medium. Following a predictable initial computational effort, Green's functions between arbitrary points in the medium can be computed as needed using a simple cross-correlation algorithm.