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.     

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.     

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 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.     

Attenuation coefficient estimation using experimental diffraction corrections with multiple interface reflections

The ultrasonic attenuation coefficient of a fluid or solid is an acoustic parameter routinely estimated for the purpose of materials characterization and defect/disease detection. This paper describes a broadband attenuation coefficient estimation technique that combines two established estimation approaches. The key elements of these two approaches are: (1) the use of magnitude spectrum ratios of front surface, first back surface, and second back surface reflections from interfaces of materials with plate-like geometries, and (2) the use of an experimental diffraction correction approach to avoid diffraction losses. The combined estimation approach simplifies the attenuation coefficient estimation process by eliminating the need to explicitly make diffraction corrections or calculate reflection/transmission coefficients. The approach yields estimates of the attenuation coefficient, reflection coefficient, and material density. Models, experimental procedures, and signal analysis procedures, which support implementation of the approach, are presented. Attenuation coefficient and reflection coefficient estimates are presented for water and solid samples with estimates based on measurements made with multiple transducers.     

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.     

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.     

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.     

Ultrasonic material characterization using large-aperture PVDF receivers

This work describes the use of a large-aperture PVDF receiver in the measurement of liquid density and composite material elastic constants. The density measurement of several liquids is obtained with accuracy of 0.2% using a conventional NDE emitter transducer and a 70-mm-diameter, 52-μm P(VDF-TrFE) membrane with gold electrodes. The determination of the elastic constants is based on the phase velocity measurement. Diffraction can lead to errors around 1% in velocity measurement when using alternatively the conventional pair of ultrasonic transducers (1-MHz frequency and 19-mm-diameter) operating in through-transmission mode, separated by a distance of 100 mm. This effect is negligible when using a pair of 10-MHz, 19-mm-diameter transducers. Nevertheless, the dispersion at 10 MHz can result in errors of about 0.5%, when measuring the velocity in composite materials. The use of an 80-mm diameter, 52-μm-thick PVDF membrane receiver practically eliminates the diffraction effects in phase velocity measurement. The elastic constants of a carbon fiber reinforced polymer were determined and compared with the values obtained by a tensile test.     

Parameter measurement of the cylindrically curved thin layer using low-frequency circumferential Lamb waves

The characteristic parameters of a cylindrically curved thin layer include its elastic constants, thickness and curved radius. A layer is considered thin if the echoes from the front and back surfaces of the layer cannot be separated in the time domain, and/or that the wave arrivals corresponding to longitudinal and shear wave part cannot be identified in the time or space domain. This paper describes a low-frequency circumferential Lamb wave method to characterize those parameters of a cylindrically curved thin layer. The technique is based on the measurement of circumferential Lamb wave phase velocity and the unknown parameter is estimated through minimizing the mean square error obtained by comparing theoretical and experimental phase velocities. The sensitivity and accuracy of the proposed technique to different parameters are analyzed. Using the present technique, a cylindrically curved thin layer with thickness down to ten percent of the longitudinal wavelength can be successfully measured with an average relative error less than two-percent in our experiment.     

Excitations of thermoelastic waves in plates by a pulsed laser

The method of the eigenfunction expansion, also known as the expansion in normal modes, is employed to study numerically the axisymmetric excitation of the thermoelastic waves in plates by a pulsed laser. This method gives a systematic treatment and allows one to investigate not only the quasistatic and dynamic thermoelastic responses of pulsed photothermal deformation on the time scale of 1 s, but also the thermoelastic generation of longitudinal, transverse, and surface acoustic waves in thick materials, as well as the excitations of the Rayleigh-Lamb wave modes in thin plates. The formalism is particularly suitable for waveform analyses of the excitations of transient Lamb waves in thin plates because one needs only to calculate the contributions of several lower eigenmodes. The numerical technique provides a quantitative tool for the experimental determination of material properties, especially the mechanical and elastic properties of free-standing films and thicker sheet materials by thermoelastic detection.     

In-line ultrasonic monitoring of semi-solid magnesium die casting process

Ultrasonic velocity and attenuation measurements in AZ91D magnesium (Mg) alloy with dendritic, rosette and globular microstructures were performed at elevated temperatures using a non-contact laser-ultrasonic technique. It was found that the ultrasonic velocity in the globular microstructure and the ultrasonic attenuation in the dendritic microstructure are the highest among the three microstructures. An ultrasonic clad steel buffer rod sensor embedded in the die has been used to monitor the semi-solid die casting process in-line for the AZ91D Mg alloy. This probe monitored the completion of the die filling, the release of the pressure, the opening of the die, part detachment, solidification of the part, the averaged temperature of the die and the part.     

Estimation of critical and viscous frequencies for Biot theory in cancellous bone

The use of Biot theory for modelling ultrasonic wave propagation in porous media involves the definition of a "critical frequency" above which both fast and slow compressional waves will, in principle, propagate. Critical frequencies have been evaluated for healthy and osteoporotic cancellous bone filled with water or marrow, using data from the literature. The range of pore sizes in bone gives rise to a critical frequency band rather than a single critical frequency, the mean of which is lower for osteoporotic bone than normal bone. However, the critical frequency is a theoretical concept and previous researchers considered a more realistic "viscous frequency" above which both fast and slow waves may be experimentally observed. Viscous frequencies in bone are found to be several orders of magnitude greater than calculated critical frequencies. Whereas two waves may well be observed at all ultrasonic frequencies for water-filled cancellous bone at 20 degrees C, it is probable megahertz frequencies would be needed for observation of two waves in vivo.     

Simultaneous determination of the ultrasound velocity and the thickness of solid plates from the analysis of thickness resonances using air-coupled ultrasound

A method that combines transmission of air-coupled ultrasound pulses through solid plates and amplitude and phase spectral analysis is presented. In particular, the method analyzes the first thickness resonance of the plates. The purpose is to determine, simultaneously, velocity and attenuation coefficient of the ultrasounds in the material and the thickness of the plate. This is especially useful when thickness can not be measured independently. The method is successfully applied to soft membranes, biological samples and FRP composites.     

Analysis of ultrasonic attenuation in particle-reinforced plastics by a differential scheme

Attenuation of ultrasonic longitudinal waves in some particle-reinforced polymer composites is studied theoretically by a micromechanical model based on a differential (incremental) scheme. A set of differential equations is established by which the attenuation spectrum of the composite can be computed from the known properties of viscoelastic matrix and elastic particles. For a composite reinforced with glass particles with radius 0.15 mm, the proposed scheme is shown to predict the attenuation in better agreement with the foregoing experimental results than the previous simplistic independent scattering model. Based on this scheme, the dependence of the longitudinal attenuation spectrum of a particulate polymer composite on the wavelength-to-particle radius ratio and the particle volume fraction is examined in detail. It is then shown theoretically that the attenuation of the composite decreases monotonically with the particle volume fraction when the particle radius is sufficiently small compared to the incident wavelength, while it shows non-monotonic particle-fraction dependence when the ratio of the particle radius to the wavelength is larger. To examine this theoretical finding from an experimental point of view, the longitudinal attenuation in a glass-particle-reinforced polyester composite with particle radius 0.0225 mm is measured for different particle volume fractions. The measured attenuation characteristics are shown to support the qualitative features of the theoretical prediction.     

Reflection and mode conversion of surface acoustic waves studied by scanning acoustic force microscopy

We present measurements of the reflection and mode conversion of surface acoustic waves (SAWs) by scanning acoustic force microscopy (SAFM). The SAFM offers a unique combination of high lateral resolution and high sensitivity towards acoustic modes of all polarizations. Since a SAW mixing experiment of two waves can be performed even if the amplitude difference between both waves is 40 dB, wavefields of extremely small amplitudes can be investigated. Using SAFM, the reflection of SAWs from a metallic wedge is investigated with submicron lateral resolution. We are able to identify two reflected wave modes, a Love and a non-coupling Rayleigh mode, by measuring their phase velocities. Received: 4 December 2000 / Accepted: 6 December 2000 / Published online: 9 February 2001     

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.     

Shock waves generated by confined XeCl excimer laser ablation of polyimide

We investigate shock waves generated by excimer laser ablation of sheet polyimide confined in water. The velocities of the ablation-induced pressure waves in the water are determined by an optical probe system. We measure supersonic velocities up to a few hundred microns away from the irradiated surface, indicating the formation of shock waves. We use these velocities to calculate the corresponding pressures. They are already in the kbar range at fluences comparable to the threshold of ablation. The shock pressure varies as the square root of the incident laser fluence, a behavior that is explained by the rapid heating of the confined gaseous products of ablation.The initially planar shock waves propagate, become spherical, and decay within a few hundred microns in the surrounding water to acoustic waves. During spherical expansion the shock pressure drops as the inverse of the square of the propagation distance.The shock waves generated may be relevant in explaining photoacoustic damage observed in biological tissue after excimer-ablation at corresponding irradiances. They may also be important in material processing applications of excimer laser ablation of polymers as they can lead to plastic deformation.     

Resonant ultrasound spectroscopy and homogeneity in polycrystals

Resonant ultrasound spectroscopy (RUS) is capable of determining the bulk elastic properties of a solid from its characteristic vibration frequencies, given the dimensions, density and shape of the sample. The model used for extracting values of the elastic constants assumes perfect homogeneity, which can be approximated by average-isotropic polycrystals. This approximation is excellent in the small grain regime assumed for most averaging procedures, but for real samples with indeterminate grain size distributions, it is not clear where the approximation breaks down. RUS measurements were made on pure copper samples where the grain size distribution was changed by progressive heat treatments in order to find a quantitative limit for the loss of homogeneity. It is found that when a measure of the largest grains is 15% of the sample’s smallest dimension, the deviation in RUS fits indicates elastic inhomogeneity.