Generalized eigenproblem of hybrid matrix for Floquet wave propagation in one-dimensional phononic crystals with solids and fluids

A method based on the solution to a generalized eigenproblem of hybrid matrix is presented for stable analysis of Floquet wave propagation in one-dimensional phononic crystals with solids and fluids. The method overcomes the numerical instability in the standard eigenproblem of transfer matrix, thus enabling Floquet waves to be determined reliably. The recursion relations of hybrid matrix for periodic multilayered structure of various solid and/or fluid phases are formulated. Dispersion relation and omnidirectional reflection for one-dimensional phononic crystals with solids and fluids are discussed. The frequency-thickness range of phononic bandgap is determined conveniently based on the Floquet wavenumbers.     

The use of an orthogonality relation for reducing the size of finite element models for 3D guided waves scattering problems

The scattering of guided waves by complex shaped defects in three-dimensional (3D) waveguides is considered. For such problems, analytical solutions do not exist, and modal decomposition techniques based on the establishment of the displacement and stress fields in the vicinity of the scatterer are quite heavy and complicated to perform. On the other hand, finite elements (FE)-based methods constitute a powerful way to obtain solutions, but they are known to be very memory consuming. This paper proposes a post-processing technique, based on a 3D orthogonality relation, to decompose a complex acoustic field produced by a scatterer and predicted by a 3D FE model, into plane waves, the amplitudes of which are quantified in the far field. This technique allows important reductions in the size of the FE models to be made. Two applications are presented to demonstrate the potential of this method. The first one concerns the scattering of the S₀ Lamb wave incident on a flat bottom circular hole. In this example, the amplitude of each mode is calculated via the orthogonality relation-based method, and compared to that obtained by simply monitoring the displacements at appropriate through-thickness positions. In the second application, the incident S₀ Lamb mode is converted into five modes scattered by a defect of complex geometry.     

A concise and efficient scattering matrix formalism for stable analysis of elastic wave propagation in multilayered anisotropic solids

This paper presents a concise and efficient scattering matrix formalism for stable analysis of elastic wave propagation in multilayered anisotropic solids. The formalism is capable of resolving completely the numerical instability problems associated with transfer matrix method, thereby obviating the extensive reformulation in its modified versions based on delta operator technique. In contrast to the earlier reflection matrix formalisms, all scattering matrices are obtained in a direct manner without invoking wave-propagator or scatterer operator concepts. Both local and global reflection and transmission matrices corresponding to scatterings in two and more layers are derived. The derivation of global scattering matrices in terms of the local ones is carried out concisely based on physical arguments to provide better insights into scattering mechanism. Another formulation which is even more succinct is also devised for obtaining the global scattering matrices directly from eigensolutions. The resultant expressions and algorithm are terse, efficient and convenient for implementation.     

Asymptotic Analysis of Vertical Branch-Cut Integral of Shear Waves in a Fluid-Filled Borehole Utilizing the Steepest-Descent Method

We obtain an asymptotic solution to the vertical branch-cut integral of shear waves excited by an impulsive pressure point source in a fluid-filled borehole, by taking the effect of the infinite singularity of the Hankel functions related to shear waves in the integrand at the shear branch point into account and using the method of steepest-descent to expand the vertical branch-cut integral of shear waves. It is theoretically proven that the saddle point of the integrand is located at ks-i/z, where ks and z are the shear branch point and the offset. The continuous and smooth amplitude spectra and the resonant peaks of shear waves are numerically calculated from the asymptotic solution. These asymptotic results are generally in agreement with the numerical integral results. It is also found by the comparison and analysis of two results that the resonant factor and the effect of the normal and leaking mode poles around the shear branch point lead to the two-peak characteristics of the amplitude spectra of shear waves in the resonant peak zones from the numerical integral calculations.     

Resonant acoustic scattering by a finite linear grating of elastic shells

Theoretical and experimental studies of the acoustic scattering by a finite linear grating of elastic cylindrical shells are performed. It is observed that a resonant interaction takes place at low frequency when the shells are very close to each other. This phenomenon can be clearly associated to the Scholte-Stoneley wave that propagates around a single shell. It is shown that each resonance of the Scholte-Stoneley wave is split up into N resonances when N shells compose the grating.     

Ultrasonic measurements on poroelastic slabs: Determination of reflection and transmission coefficients and processing for Biot input parameters

Ultrasonic reflection and transmission measurements on saturated, plane-parallel, poroelastic slabs are presented. A data processing technique is proposed to obtain poroelastic parameters from transmission measurements. A special experimental data acquisition and processing technique is applied to minimize the finite beam effects of the transducers. This technique yields results which can be compared with poroelastic plane-wave theory. Results of normal- and oblique-incidence measurements are presented, and transmission data are processed to yield wave speeds, sample thickness, angle of incidence, tortuosity, and permeability. The results show good agreement with independent measurements, and they are subsequently used as input for a forward modeling of the complete transmitted and reflected waveforms utilizing Biot theory. The agreement between recorded and modeled signals is good, both in time and frequency domain.     

Perturbation Analysis on Guided Waves in a Fluid-Filled Borehole Surrounded by a Cubic Crystal Anisotropic Medium

The perturbation method is employed to analyse the guided waves in a borehole surrounded by a cubic crystal medium for the first time. The cubic crystal medium is regarded as a reference unperturbed isotropic state added to the perturbation. The dispersion characteristics of Stoneley wave, pseudo-Rayleigh wave, flexural wave, and screw wave are investigated in detail. It is found that dispersion of the guided waves excited by monopole and dipole sources does not depend on the azimuth of the source, whereas the dispersion of screw wave excited by quadrupole source is significantly related to the azimuth of the source. Screw waves propagated along different azimuth in the borehole can be split. This is different from screw waves in transversely isotropic media （hexagonal crystal）, which have been widely studied.     

Normal modes of a poroelastic plate and their relation to the reflection and transmission coefficients

A poroelastic plate that obeys the Biot theory is considered. Compact new forms of its reflection and transmission coefficients, similar to those of the resonance scattering theory for an elastic plate, are derived. A numerical comparison of the reflection coefficient modulus with the plate normal modes, at low frequency, shows that a study of the reflection or transmission coefficient does not provide the same kind of information on the poroelastic plate than an investigation of guided leaky waves propagation.     

Minimization of finite beam effects in the determination of reflection and transmission coefficients of an elastic layer

A special data acquisition technique was applied to determine the acoustic plane-wave reflection and transmission properties of a plane-parallel aluminum plate. In this technique, the reflected and transmitted wavefield along a plane or line normal to the reflected or transmitted wave vector is recorded at equidistant receiver positions. The obtained traces are subsequently added up in the temporal domain to satisfy plane-wave conditions, thus effectively removing the effect of the limited beam of commonly used transducers. The agreement between plane-wave theory and experiment was found to be excellent, both in the temporal and in the frequency domain.     

Phase Transition in Acoustic Localization in a Soft Medium Permeated with Air Bubbles

Propagation of an acoustic wave in a soft medium permeated with air bubbles is theoretically investigated by using a self-consistent approach. The soft medium is assumed to be viscoelastic to estimate the effect of acoustic absorption on the acoustic localization in such a medium. The oscillation phases of bubbles are examined by employing a phase diagram method. A collective oscillation of the bubbles is observed once the acoustic localization occurs, which is known as a phenomenon of ＇phase transition ＇, and such a phenomenon persists as we manually increase the viscosity factor of the soft medium. Therefore it is proven that the phenomenon of phase transition may serve as a unique criterion to effectively identify acoustic localization in a bubbly soft medium even in the presence of viscosity, and the directions of the phase vectors help to determine the extent of localization. This is of practical significance for experimental research studying the acoustic localization in such a medium, for which the presence of viscosity generally causes great ambiguity in distinguishing the effects of localization and acoustic absorption.     

Impact of localized inhomogeneity on the surface-wave velocity and bulk-wave reflection in solids

The effect of a weak surface, near-surface and interfacial inhomogeneity on the frequency dependence of the surface wave velocity and of the SH (shear horizontal) wave reflectivity in isotropic elastic media is studied analytically and numerically. The inhomogeneity is modeled as an infinite planar layer with continuously varying properties. Weak inhomogeneity may markedly affect the dispersion of the Rayleigh velocity and especially of the reflectivity. It is demonstrated how this effect, particularly pronounced at high frequency, depends on the extent of inhomogeneity. The material data for damaged and ideal concrete and several simple examples of inhomogeneity profiles are utilized for the numerical calculations based on the Peano expansion. The use of explicit low- and high-frequency approximations is also exemplified. Among these, simple WKB asymptotics are shown to be particularly helpful for the Rayleigh velocity in the case of a prominent inhomogeneity attached to the surface and for the reflection on weak interfaces.     

Ultrasonic wave propagation in heterogeneous solid media: theoretical analysis and experimental validation

This research deals with the ultrasonic characterization of thermal damage in concrete. This damage leads to the appearance of microcracks which then evolve in terms of volume rate and size in the material. The scattering of ultrasonic waves from the inclusions is present in this type of medium. The propagation of the longitudinal wave in the heterogeneous media is studied via a homogenization model that integrates the multiple scattering of waves. The model allows us to determine the phase velocity and the attenuation according to the elements which make the medium. Simulations adapted to the concrete are developed in order to test the responses of the model. These behaviors are validated by an experimental study: the measurements of phase velocity and attenuation are performed in immersion, with a comparison method, on a frequency domain which ranges from 160 kHz to 1.3 MHz. The analysis of different theoretical and experimental results obtained on cement-based media leads to the model validation, on the phase velocity behavior, in the case of a damage simulated by expanded polystyrene spheres in granular media. The application to the case of a thermally damaged concrete shows a good qualitative agreement for the changes in velocity and attenuation.     

Rapid computation of acoustic impulse scattering for rough penetrable surfaces

Chu’s theory for the impulse response of a point source to an isovelocity density contrast wedge [Chu D. Impulse response of density contrast wedge using normal coordinates. J Acoust Soc Am 1989;86:1883-96] enables wedge-assemblage rough surface scattering models to be extended to a broad range of penetrable seafloors, but is computationally intensive since it necessitates finding the multifold roots of a characteristic eigenvalue equation, and summing a power series, for each wedge apex. This present work considers the properties and relationships of the direct, reflected, and diffracted field components of a density contrast wedge. In particular, an analysis of the physical origin and behavior of diffractions associated with specular reflections of the source in the wedge faces leads to a simple extension of the Biot-Tolstoy theory [Biot MA, Tolstoy I. Formulation of wave propagation in infinite media by normal coordinates with an application to diffraction. J Acoust Soc Am 1957;29:381-91] to density contrast wedges with reflectivity ∣R∣ < 1, for wedge angles within the range 150 ? θ_w ? 210°, where the diffractions are predominantly associated with a single reflection in each wedge face. This facilitates rapid time domain calculations of acoustic bottom scattering and penetration for complex multilayered seafloors.     

Reflection of structural waves at a solid/liquid interface

This paper investigates the reflection characteristics of structural or guided waves in rods at a solid/liquid interface. Structural waves, whose wavelengths are much larger than the diameter of the rod, are described in a first approximation by classical one-dimensional wave theory. The reflection characteristics of such waves at a solid/liquid (melting) interface has been reported by two different ultrasonic measurement techniques: first, measuring the fast regression rate of a melting interface during the burning of metal rod samples in an oxygen-enriched environment, and second, monitoring the propagation of the solid/liquid interface during the slow melting and solidification of a rod sample in a furnace. The second work clearly shows that the major reflection occurs from the solid/liquid interface and not the liquid/gas interface as predicted by plane longitudinal wave reflectivity theory. The present work confirms this observation by reporting on the results of some specially designed experiments to identify the main interface of reflection for structural waves in rods. Hence, it helps in explaining the fundamental discrepancy between the reflection characteristics at a solid/liquid interface between low frequency structural waves and high frequency bulk waves, and confirms that the detected echo within a burning metallic rod clearly represents a reflection from the solid/liquid interface.     

Multi-focus design of underwater noise control linings based on finite element analysis

Varied, counter-demanding objectives in designing the underwater noise control linings are addressed using a finite element model based methodology. Four different kinds of designs are proposed to attend to diverse and conflicting requirements concerning echo reduction (ER) and transmission loss (TL) performance of these linings. In this regard a slightly modified hybrid type finite element based on the Pian and Tong (PT) formulation has been used to make the computational efforts less demanding as compared to the original one. The adequacy of this formulation has been shown by comparing its results with the analytical, finite element analysis based, and experimental results. Different unit cell representations for different types of distributions of air cavities on the linings are discussed with respect to their limitations and applicability. Effect of static pressure is studied by using a simplified technique which can be used to simulate deep sea testing environment. Performance variation of different designs is investigated under different water depths to study their applicability in such situations.     

Synthetic aperture focusing for defect reconstruction in anisotropic media

Ultrasonic inspection plays an important role in numerous industrial fields. One of the prominent tasks with respect to quantitative nondestructive evaluation is the determination of location, shape, size and orientation of defects. In this respect, the synthetic aperture focusing technique (SAFT) has been successfully applied to isotropic materials over the years. In anisotropic media, however, its application suffers from several phenomena, which are the direction dependence of the ultrasonic velocities, the beam skewing effect and the modified transducer radiation characteristics. In this article, a SAFT imaging algorithm is presented which fully accounts for the nature of wave radiation and propagation within anisotropic materials. For three-dimensional defect reconstruction, the spatial dependence of the ultrasonic group velocities as well as the radiation characteristics of the transducer are exploited--respective algorithms have been implemented for orthotropic material symmetry. Tests have been performed on unidirectional composite material.     

A useful analytical description of the coefficients in an inhomogeneous wave decomposition of a symmetrical bounded beam

If a bounded beam is described using a superposition of infinite inhomogeneous waves, the values of the coefficients attributed to each inhomogeneous wave are found using a classical optimization procedure, whence it is impossible to describe the obtained values analytically. In this paper, we develop a new and easy to apply straightforward analytical method to find the appropriate values of the sought coefficients. Supplementary to its analytical and straightforward nature, the method proves to reduce the inherent instabilities found in the inhomogeneous wave decomposition.     

A New Method for Calculation of Single Seismic Phase of Cylindrically Multilayered Media Including Liquid Interlayer

A new method based on generalized reflection and transmission （R/T） coefficients method is proposed to calculate the single seismic phase （SSP） of cylindrically multilayered media including liquid interlayer. The use of normalization factors and normalized Lamé coefficients makes the algorithm stable numerically. Using the modified R/T matrices, we derive the iterative expressions of generalized R/T matrices, and by using the iterative relation we determine the SSP of each interface and the full waveforms. To show the superiority of this new approach for investigating of reflection and transmission properties of cylindrically multilayered media, we simulate the full waveforms and SSPs of cased hole model with annulus I （casing-cement interface） channelling （or, cross-flow）. The generalized reflection coefficient spectra and SSPs of interfaces obtained show the propagation mechanism of each component of full waveform clearly.     

Ultrasonic sizing of short surface cracks

This paper presents a method for ultrasonic sizing of surface cracks based on time domain and frequency domain Rayleigh wave near-field analysis. The procedure allows for the entire range of ratio of crack depth to Rayleigh wavelength a/λ to be covered with one single measurement. In the time domain the time-of-flight method was extended to cracks smaller than the wavelength by correlation of the time delay of the transmitted Rayleigh wave with the crack depth. In the frequency domain, the inverse scattering problem was solved by comparison of the measured scattering coefficients and central frequencies of the reflected and transmitted Rayleigh waves with theoretical curves. The sizing procedure was demonstrated experimentally with narrow slots and real fatigue cracks. The out-of-plane displacement component was measured pointwise in the scattered near field by means of laser interferometry. The determination of the scattering parameters in the near field was enabled by a procedure that allows for the Rayleigh wave to be separated from the other modes scattered at the defect. The experimental results showed good accuracy and repeatability down to the smallest available ratio of crack depth to Rayleigh wavelength a/λ = 0.15.     

Electromagnetic wave scattering by small bodies

A reduction of the Maxwell's system to a Fredholm second-kind integral equation with weakly singular kernel is given for electromagnetic (EM) wave scattering by one and many small bodies. This equation is solved asymptotically as the characteristic size of the bodies tends to zero. The technique developed is used for solving the many-body EM wave scattering problem by rigorously reducing it to solving linear algebraic systems, completely bypassing the usage of integral equations. An equation is derived for the effective field in the medium, in which many small particles are embedded. A method for creating a desired refraction coefficient is outlined.