Ultrasonic wave propagation in stereo-lithographical bone replicas

Three-dimensional-printed replica models of bone for experimentally decoupling trabecular bone properties contribution to ultrasound propagation parameters

A detailed investigation of the relationship between ultrasonic (US) properties and trabecular bone microstructure is difficult because of the great variability in the bone loss process. The aim of this work was twofold. First, to verify by compressive tests that the three-dimensional (3D)-printer is able to produce precisely and repeatedly "bone replica models" of different size and density. Following, replicas of the original specimens with two different polymers and thinned trabeculae models were used to investigate US properties (speed of sound, SOS, and backscatter coefficient), aiming to deconvolute the influence of material properties on ultrasound characteristics. The results revealed that matrix material properties influence only the magnitude of the backscatter coefficient, whereas the characteristic undulated patterns are related to the trabecular structure. Simulation of perforation and thinning of cancellous bone, associated with bone loss, showed that SOS and mechanical properties were reduced perfectly linearly with apparent density when structure deteriorated. The 3D-printed bone replicas have the potential to enable systematic investigations of the influence of structure on both acoustical and mechanical properties and evaluate changes caused by bone loss. The development of replicas from materials with properties close to those of bone will permit quantitative conclusions for trabecular bone.

Evaluation of Structural Anisotropy in a Porous Titanium Medium Mimicking Trabecular Bone Structure Using Mode-Converted Ultrasonic Scattering

The mode-converted (longitudinal to transverse, L-T) ultrasonic scattering method was utilized to characterize the structural anisotropy of a phantom mimicking the structural properties of trabecular bone. The sample was fabricated using metal additive manufacturing from high-resolution computed tomography (CT) images of a sample of trabecular horse bone with strong anisotropy. Two focused transducers were used to perform the L-T ultrasonic measurements. A normal incidence transducer was used to transmit longitudinal ultrasonic waves into the sample, while the scattered transverse signals were received by an oblique incidence transducer. At multiple locations on the sample, four L-T measurements were performed by collecting ultrasonic scattering from four directions. The amplitude of the root mean square (rms) of the collected ultrasonic scattering signals was calculated for each L-T measurement. The ratios of rms amplitudes for L-T measurements in different directions were calculated to characterize the anisotropy of sample. The results show that the amplitude of L-T converted scattering is highly dependent on the direction of microstructural anisotropy. A strong anisotropy of the microstructure was observed, which coincides with simulation results previously published on the same structure as well as with the anisotropy estimated from the CT images. These results suggest the potential of mode-converted ultrasonic scattering methods to assess the anisotropy of materials with porous, complex structures, including trabecular bone.

Mechanisms of Interaction of Ultrasound With Cancellous Bone: A Review

Ultrasound is now a clinically accepted modality in the management of osteoporosis. The most common commercial clinical devices assess fracture risk from measurements of attenuation and sound speed in cancellous bone. This review discusses fundamental mechanisms underlying the interaction between ultrasound and cancellous bone. Because of its two-phase structure (mineralized trabecular network embedded in soft tissue-marrow), its anisotropy, and its inhomogeneity, cancellous bone is more difficult to characterize than most soft tissues. Experimental data for the dependencies of attenuation, sound speed, dispersion, and scattering on ultrasound frequency, bone mineral density, composition, microstructure, and mechanical properties are presented. The relative roles of absorption, scattering, and phase cancellation in determining attenuation measurements in vitro and in vivo are delineated. Common speed of sound metrics, which entail measurements of transit times of pulse leading edges (to avoid multipath interference), are greatly influenced by attenuation, dispersion, and system properties, including center frequency and bandwidth. However, a theoretical model has been shown to be effective for correction for these confounding factors in vitro and in vivo. Theoretical and phantom models are presented to elucidate why cancellous bone exhibits negative dispersion, unlike soft tissue, which exhibits positive dispersion. Signal processing methods are presented for separating "fast" and "slow" waves (predicted by poroelasticity theory and supported in cancellous bone) even when the two waves overlap in time and frequency domains. Models to explain dependencies of scattering on frequency and mean trabecular thickness are presented and compared with measurements. Anisotropy, the effect of the fluid filler medium (marrow in vivo or water in vitro), phantoms, computational modeling of ultrasound propagation, acoustic microscopy, and nonlinear properties in cancellous bone are also discussed.

Time-domain separation of interfering waves in cancellous bone using bandlimited deconvolution: simulation and phantom study

In through-transmission interrogation of cancellous bone, two longitudinal pulses ("fast" and "slow" waves) may be generated. Fast and slow wave properties convey information about material and micro-architectural characteristics of bone. However, these properties can be difficult to assess when fast and slow wave pulses overlap in time and frequency domains. In this paper, two methods are applied to decompose signals into fast and slow waves: bandlimited deconvolution and modified least-squares Prony's method with curve-fitting (MLSP + CF). The methods were tested in plastic and Zerdine(®) samples that provided fast and slow wave velocities commensurate with velocities for cancellous bone. Phase velocity estimates were accurate to within 6 m/s (0.4%) (slow wave with both methods and fast wave with MLSP + CF) and 26 m/s (1.2%) (fast wave with bandlimited deconvolution). Midband signal loss estimates were accurate to within 0.2 dB (1.7%) (fast wave with both methods), and 1.0 dB (3.7%) (slow wave with both methods). Similar accuracies were found for simulations based on fast and slow wave parameter values published for cancellous bone. These methods provide sufficient accuracy and precision for many applications in cancellous bone such that experimental error is likely to be a greater limiting factor than estimation error.

Estimation of fast and slow wave properties in cancellous bone using Prony's method and curve fitting

The presence of two longitudinal waves in poroelastic media is predicted by Biot's theory and has been confirmed experimentally in through-transmission measurements in cancellous bone. Estimation of attenuation coefficients and velocities of the two waves is challenging when the two waves overlap in time. The modified least squares Prony's (MLSP) method in conjuction with curve-fitting (MLSP + CF) is tested using simulations based on published values for fast and slow wave attenuation coefficients and velocities in cancellous bone from several studies in bovine femur, human femur, and human calcaneus. The search algorithm is accelerated by exploiting correlations among search parameters. The performance of the algorithm is evaluated as a function of signal-to-noise ratio (SNR). For a typical experimental SNR (40 dB), the root-mean-square errors (RMSEs) for one example (human femur) with fast and slow waves separated by approximately half of a pulse duration were 1 m/s (slow wave velocity), 4 m/s (fast wave velocity), 0.4 dB/cm MHz (slow wave attenuation slope), and 1.7 dB/cm MHz (fast wave attenuation slope). The MLSP + CF method is fast (requiring less than 2 s at SNR = 40 dB on a consumer-grade notebook computer) and is flexible with respect to the functional form of the parametric model for the transmission coefficient. The MLSP + CF method provides sufficient accuracy and precision for many applications such that experimental error is a greater limiting factor than estimation error.