Samuel Callé

Exploration of trabecular bone nonlinear elasticity using time-of-flight modulation

Bone tissue contains microcracks that may affect its mechanical properties as well as the whole trabecular structure. The relationship between crack density and bone strength is nevertheless poorly understood. Linear ultrasound techniques being almost insensitive to the level of damage, we propose a method to measure acoustic non- linearity in trabecular bone using time-of-flight modulation (TOFM) measurements. Ultrasonic short bursts times-of- flight (TOF) are modulated as a result of nonlinear interaction with a low-frequency (LF) wave in the medium. TOF variations are directly related to elastic modulus variations. Classical and nonclassical nonlinear parameters beta, delta, and alpha can be derived from these measurements. The method was validated in materials with classical, quadratic, nonlinear elasticity. In dense trabecular bone region, TOFM related to classical, quadratic, nonlinear elasticity as a function of the LF pressure exhibits tension-compression asymmetry. The TOFM amplitude measured in dense areas of trabecular bone is almost one order of magnitude higher than in a low-density area, but the linear parameters show much smaller variations: 5% for ultrasound propagation velocity and 100% for broadband ultrasonic attenuation (BUA). In high-density trabecular bone regions, beta depends on the LF pressure amplitude and can reach 400 at 50 kPa.

A two-dimensional pseudospectral model for time reversal and nonlinear elastic wave spectroscopy

One way to characterize metallic materials in the presence of defects like dislocation networks is to measure their large dynamic nonlinear elastic response. In this numerical study, a new method combining the nonlinear elastic wave spectroscopy (NEWS) method with a time reversal (TR) process is proposed. This method, called NEWS-TR, uses nonlinear analysis as a pretreatment of time reversal and then consists of retrofocusing only nonlinear components on the defect position. A two-dimensional pseudospectral time domain algorithm is developed here to validate the NEWS-TR method as a potential technique for damage location. Hysteretic nonlinear behavior of the materials being studied is introduced using the Preisach-Mayergoyz model. Moreover, in order to extend this solver in two dimensions, the Kelvin notation is used to modify the elastic coefficient tensor. Simulations performed on a metallic sample show the feasibility and value of the NEWS-TR methodology for microdamage imaging. Retrofocusing quality depends on different parameters such as the filtering method used to keep only nonlinear components and the nonlinear effect measured. In harmonic generation, pulse inversion filtering seems to be a more appropriate filtering method than classical harmonic filtering for most defect positions, mainly because of its ability to filter all fundamental components.

Temporal analysis of tissue displacement induced by a transient ultrasound radiation force

One of the stress sources that can be used in dynamic elastography imaging methods is the acoustic radiation force. However, displacements of the medium induced by this stress field are generally not fully understood in terms of spatial distribution and temporal evolution. A model has been developed based on the elastodynamic Greens function describing the different acoustic waves generated by focused ultrasound. The function is composed of three terms: two far-field terms, which correspond to a purely longitudinal compression wave and a purely transverse shear wave, and a coupling near-field term which has a longitudinal component and a transverse component. For propagation distances in the shear wavelength range, the predominant term is the near field term. The displacement duration corresponds to the propagation duration of the shear wave between the farthest source point and the observation point. This time therefore depends on the source size and the local shear modulus of the tissue. Evolution of the displacement/time curve profile, which is directly linked to spatial and temporal source profiles, is computed at different radial distances, for different durations of force applications and different shear elastic coefficients. Experimental results performed with an optical interferometric method in a homogeneous tissue-mimicking phantom agreed with the theoretical profiles.

Nonlinear elastodynamics in micro-inhomogeneous solids observed by head-wave based dynamic acoustoelastic testing

Dynamic acoustoelastic testing provides a more complete insight into the acoustic nonlinearity exhibited by micro-inhomogeneous media like granular and cracked materials. This method consists of measuring time of flight and energy modulations of pulsed ultrasonic waves induced by a low-frequency standing wave. Here pulsed ultrasonic head waves were employed to assess elastic and dissipative nonlinearities in a region near the surface of a solid. Synchronization of the ultrasound pulse sequence with the low-frequency excitation provided instantaneous variations in the elastic modulus and the attenuation as functions of the instantaneous low-frequency strain. Weak quadratic elastic nonlinearity and no dissipative nonlinearity were detected in duralumin. In limestone, distinction between tensile and compressive behaviors revealed an asymmetry in the acoustic nonlinearity and hysteresis in both the elastic modulus and the attenuation variations. Measured nonlinear acoustical parameters are in good agreement with values obtained by different techniques. Reversible acoustically induced conditioning modified the acoustic nonlinearity both quantitatively and qualitatively. It reduced tension-compression asymmetry, suggesting a nonequilibrium modification of the sources of acoustic nonlinearity. Additionally to the metrology of the acoustic nonlinearity, head wave based dynamic acoustoelastic testing may be a useful tool to monitor changes in the microstructure or the accumulation of damage in solids.

Monitoring trabecular bone microdamage using a dynamic acousto-elastic testing method

Dynamic acousto-elastic testing (DAET) is based on the coupling of a low-frequency (LF) acoustic wave and high-frequency ultrasound (US) pulses (probing wave). It was developed to measure US viscoelastic and dissipative non-linearity in trabecular bone. It is well known that this complex biphasic medium contains microdamage, even when tissues are healthy. The purpose of the present study was to assess the sensitivity of DAET to monitor microdamage in human calcanei. Three protocols were therefore performed to investigate the regional heterogeneity of the calcaneus, the correlation between DAET measurements and microdamage revealed by histology, and DAET sensitivity to mechanically induced fatigue microdamage. The non-linear elastic parameter β was computed for all these protocols. The study demonstrated the presence of high viscoelastic and dissipative non-linearity only in the region of the calcaneus close to the anterior talocalcaneal articulation (region of high bone density). Protocols 1 and 2 also showed that most unsorted calcanei did not naturally exhibit high non-linearity, which is correlated with a low level of microcracks. Nevertheless, when microdamage was actually present, high levels of US non-linearity were always found, with characteristic non-linear signatures such as hysteresis and tension/compression asymmetry. Finally, protocol 3 demonstrated the high sensitivity of DAET measurement to fatigue-induced microdamage.

Modeling of piezoelectric transducers with combined pseudospectral and finite-difference methods.

A new hybrid finite-difference (FD) and pseudospectral (PS) method adapted to the modeling of piezoelectric transducers (PZTs) is presented. The time-dependent equations of propagation are solved using the PS method while the electric field induced in the piezoelectric material is determined through a FD representation. The purpose of this combination is to keep the advantages of both methods in one model: the adaptability of FD representation to model piezoelectric elements with various geometries and materials, and the low number of nodes per wavelength required by the PS method. This approach is implemented to obtain an accurate algorithm to simulate the propagation of acoustic waves over large distances, directly coupled to the calculation of the electric field created inside the piezoelectric material, which is difficult with classical algorithms. These operations are computed using variables located on spatially and temporally staggered grids, which attenuate Gibbs phenomenon and increase the algorithms accuracy. The two-dimensional modeling of a PZT plate excited by a 50 MHz sinusoidal electrical signal is performed. The results are successfully compared to those obtained using the finite-element (FE) algorithm of ATILA software with configurations spatially and temporally adapted to the FE requirements. The cost efficiency of the FD-PS time-domain method is quantified and verified.

Dynamic acoustoelastic testing of weakly pre-loaded unconsolidated water-saturated glass beads

Dynamic acoustoelastic testing is applied to weakly pre-loaded unconsolidated water-saturated glass beads. The gravitational acceleration produces, on the probed beads, a static stress of order 130 Pa, thus the granular medium is close to the jamming transition. A low-frequency (LF) acoustic wave gently disturbs the medium, inducing successively slight expansion and compaction of the granular packing expected to modulate the number of contacts between beads. Ultrasound (US) pulses are emitted simultaneously to dynamically detect the induced modification of the granular skeleton. US propagation velocity and attenuation both increase when the LF pressure increases. The quadratic nonlinear elastic parameter β, related to the pressure dependence of US propagation velocity, was measured in the range 60-530 if water-saturated glass beads are considered as an effective medium. A dynamic modification of US scattering induced by beads is proposed to modulate US attenuation. Complex hysteretic behaviors and tension-compression asymmetry are also observed and analyzed by time-domain and spectral analyses. Furthermore acoustic nonlinearities are measured in cases of quasi-static and dynamic variations of the LF wave amplitude, providing quantitatively similar acoustic nonlinearities but qualitatively different.

Modeling nonlinear viscoelasticity in dynamic acoustoelasticity

Dynamic Acoustoelastic Testing (DAET) has been developed to non-invasively assess the nonlinear viscoelastic properties of fluids and solids. We propose a phenomenologically motivated model for harmonic regime to identify nonlinear viscoelastic parameters from DAET experiments. The nonlinear elastic and viscous parameters are derived from the real and imaginary components of the Taylor series expansion of the complex longitudinal modulus. The model is applied to Newtonian fluids that exhibit classical elastic nonlinearities and glass beads saturated with water that exhibit nonclassical viscoelastic nonlinearities. Hysteresis, asymmetry, and DC offset are well accounted for in the model.

Single-element ultrasonic transducer modeling using a hybrid FD-PSTD method.

In a recent publication [E. Filoux, S. Callé, D. Certon, M. Lethiecq, F. Levassort, Modeling of piezoelectric transducers with combined pseudospectral and finite-difference methods, J. Acoust. Soc. Am. 123 (6) (2008) 4165-4173], a new finite-difference/pseudospectral time-domain (FD-PSTD) algorithm was presented and used to model the generation of acoustic waves by a piezoelectric resonator and their propagation in the structure and the surrounding water. In this paper, the model has been extended to simulate the two-dimensional behaviour of a complete single-element transducer, composed of the resonator, a backing and a front matching layer. This further version of the model takes into account the mechanical loss in materials, and enables the calculation of electrical impedance, which is a characteristic of high interest to optimize the performance of ultrasonic transducers. The impedance curves of a PZT [URL: http://www.ferroperm-piezo.com (last viewed 04/2008); B. Jaffe, R.S. Roth, S. Marzullo, Piezoelectric properties of lead zirconate-lead titanate solid-solution ceramics, J. Appl. Phys. 25 (1954) 809-810] plate-based high-frequency transducer, with a 50 MHz thickness resonant frequency, were compared to those of a KLM model [R. Krimholtz, D.A. Leedom, G.L. Matthei, New equivalent circuit for elementary piezoelectric transducers, Electron. Lett. 6 (1970) 398-399] in the one-dimensional case. The acoustical properties were also found to be in good agreement with those obtained using the finite element (FE) method of ATILA software in two-dimensional configuration.

Optimised properties of high frequency transducers based on curved piezoelectric thick films obtained by pad printing process

Abstract Abstract A high frequency transducer for medical imaging (25 MHz) was fabricated using a pad printing process to deposit a curved lead zirconate titanate (PZT) thick film on electroded backing (porous PZT). This piezoelectric thick film was characterised, and a thickness coupling factor (47%) comparable with that of a bulk ceramic with similar composition was measured. This transducer was successfully modelled with a numerical tool previously published and specifically adapted to curved shapes. The experimental axial and lateral resolutions are 40 and 230 μm respectively. Moreover, the sensitivity is sufficiently high to consider this transducer to be integrated in an echographic system for high frequency imaging such as skin.