Investigation of an anisotropic tortuosity in a biot model of ultrasonic propagation in cancellous bone

Mathematical modeling of vancomycin release from Poly-L-Lactic Acid-Coated implants

This study aimed to develop a mathematical model to predict the release profile and antibacterial efficacy of a vancomycin delivery system integrated with poly(L-lactic acid)-coated bone implants specifically designed for bone plates. Using Fickian diffusion principles within an ANSYS-CFX computational fluid dynamic model, we validated the model against our in vitro vancomycin release and agar diffusion studies, as well as previously published in vivo data, confirming the reliability of the model. The model predictions demonstrated the effectiveness of the system in inhibiting bacterial growth in surrounding tissue with no observed toxicity, with a peak vancomycin concentration of 0.95 mg/ml at 6 hours, followed by a decrease to levels that remained effective for antibacterial activity. Furthermore, a sensitivity analysis revealed that the model is particularly sensitive to the half-life of vancomycin, with a maximum sensitivity index of 0.8, indicating its greater impact on the prediction accuracy than the diffusion coefficient, which has a maximum sensitivity index of 0.5. Therefore, precise input of vancomycin's half-life is critical for accurate predictions. These findings offer substantial support for the efficacy of the local delivery system as a promising therapeutic approach against implant-associated infections.

Prediction of osteoporotic degradation of tibia human bone at trabecular scale

A theoretical numerical model is proposed to predict patient dependent osteoporotic bone degradation. The model parameters are identified through a particle swarm optimization algorithm and based on individual patient high resolution peripherical quantitative computer tomography (HRpQCT) scan data. The degradation model is based on cellular activity initiated by the elastic strain energy developed in the bone microstructure through patient's body weight. The macro (organ scale) and meso (trabecular scale) scale analyses are carried out and predicted bone volume fraction and microstructure evolution are compared with in-vivo experimental bone degradation for four elderly women over a period of 10 years. A significant correlation (r > 0.9) is observed between the model predictions and in-vivo experiments in all cases with an average deviation error of 1.46%. The model can easily be extended to other patients and provide good predictions for different population categories such as ethnicity, gender, age, etc.

Transient Propagation of Longitudinal and Transverse Waves in Cancellous Bone: Application of Biot Theory and Fractional Calculus

In this paper, the influence of the transverse wave on sound propagation in a porous medium with a flexible structure is considered. The study is carried out in the time domain using the modified Biot theory obtained by the symmetry of the Lagrangian (invariance by translation and rotation). The viscous exchanges between the fluid and the structure are described by fractional calculus. When a sound pulse arrives at normal incidence on a porous material with a flexible structure, the transverse waves interfere with the longitudinal waves during propagation because of the viscous interactions that appear between the fluid and the structure. By performing a calculation in the Laplace domain, the reflection and transmission operators are derived. Their time domain expressions depend on the Green functions of the longitudinal and transverse waves. In order to study the effects of the transverse wave on the transmitted longitudinal waves, numerical simulations of the transmitted waves in the time domain by varying the characteristic parameters of the medium are realized whether the transverse wave is considered or not.

Influence of Higher Order Viscous and Thermal Effects on an Ultrasonic Wave Reflected from the First Interface of a Porous Material

Ultrasound propagation in porous materials involves some higher order physical parameters whose importance depends on the acoustic characteristics of the materials. This article concerns the study of the influence of two parameters recently introduced, namely, the viscous and thermal surfaces, on the acoustic wave reflected by the first interface of a porous material with a rigid structure. These two parameters describe the fluid/structure interactions in a porous medium during the propagation of the acoustic wave in the high-frequency regime. Both viscous and thermal surfaces are involved in Laurent expansion, which is limited to the dynamic tortuosity and compressibility to a higher order and corrects the visco-thermal losses. A sensitivity study is performed on the modulus of the reflection coefficient at the first interface as a function of frequency and on the waveforms reflected by the porous material in the time domain. The results of this study show that highly absorbent porous materials are the most sensitive to viscous and thermal surfaces, which makes the consideration of these two parameters paramount for the characterization of highly absorbent porous materials using the waves reflected from the first interface.

Bayesian inference of human bone sample properties using ultrasonic reflected signals

The non-intrusiveness and low cost of ultrasonic interrogation is motivating the development of new means of detection of osteoporosis and other bone deficiencies. Bone is a porous media saturated with a viscous fluid and could thus be well characterized by the Biot model. The main purpose of this work is to present an in vitro methodology for the identification of the properties and structural parameters of the bone, adopting a statistical Bayesian inference technique using ultrasonic reflected signals at normal incidence. It is, in this respect, a companion paper to a previous work [J. Acoust. Soc. Am. 146, 3 (2019), pp. 1629-1640], where ultrasonic transmitted signals were considered. This approach allows the retrieval of some important parameters that characterize the bone structure and associated uncertainties. The method was applied to seven samples of bone extracted from femoral heads, immersed in water, and exposed to ultrasonic signals with a center frequency of ≈500 kHz. For all seven samples, signals at different sites were acquired to check the method robustness. The porosity, pore mean size and standard deviation, and the porous frame bulk modulus were all successfully identified using only ultrasonic reflected signals.

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.

Bayesian inference of a human bone and biomaterials using ultrasonic transmitted signals

Ultrasonic techniques could be good candidates to aid the assessment of osteoporosis detection, due to their non-intrusiveness and low cost. While earlier studies made use of the measured ultrasonic phase velocity and attenuation inside the bone, very few have considered an inverse identification of both the intrinsic pore microstructure and the mechanical properties of the bone, based on Biot's model. The main purpose of this work is to present an in vitro methodology for bone identification, adopting a statistical Bayesian inference technique using ultrasonic transmitted signals, which allows the retrieval of the identified parameters and their uncertainty. In addition to the bone density, Young's modulus and Poisson's ratio, the bone pore microstructure parameters (porosity, tortuosity, and viscous length) are identified. These additional microstructural terms could improve the knowledge on the correlations between bone microstructure and bone diseases, since they provide more information on the trabecular structure. In general, the exact properties of the saturating fluid are unknown (bone marrow and blood in the case of bone study) so in this work, the fluid properties (water) are identified during the inference as a proof of concept.

A method to locate spatial distribution of scattering centers from ultrasonic backscatter signal

The purpose of this work is to find a method to locate the scattering centers in spatial domain; by using this information, the mean scatter spacing (MSS) can be estimated, and the spatial information is the one-dimensional imaging of scattering centers. This paper presents a method that can locate the scattering centers in spatial domain robustly and automatically. By incorporating it with fast Fourier transformation, the MSS can be estimated. The three foremost processes, matched filtering, envelope extraction, and peak reconstruction, are incorporated in the authors' algorithm. Monte Carlo simulations demonstrate that the proposed method is a robust one to locate scattering centers in spatial domain, and has a better performance than spectrum-based MSS estimation techniques. Especially exploited in estimating MSS which varies from 0.6 to 1.2 mm in the range of human mean trabecular bone spacing, the proposed method shows great potential in medical use. Simple but widely used phantom experiments demonstrate that the proposed algorithm has the capacity to locate scattering centers in spatial domain.

Ultrasound temporal-spatial phase-interference in complex composite media; a comparison of experimental measurement and simulation prediction

The propagation of ultrasound through solid:liquid complex composite media such as cancellous bone suffers from a lack of a comprehensive understanding of the dependence upon density and structure. Assuming that a propagating ultrasound wave may be considered as an array of parallel sonic rays, we may determine the transit time of each by the relative proportion of the two constituents. A transit time spectrum (TTS) describes the proportion of sonic rays having a particular transit time between the minimum (tmin) and maximum (tmax) values; representing, for example, entire bone tissue and marrow respectively in the case of cancellous bone. Langton has proposed that the primary ultrasound attenuation mechanism in such media is phase-interference. The phase-interference of two or more ultrasound pulses detected at a phase-sensitive transducer has both temporal and spatial components. The temporal component is primarily dependent upon the transit time difference (dt) between the pulses and the propagating pulse-length (PL). The spatial component is primarily dependent upon the lateral separation (ds) of the detectedpulses of differing transit time and the lateral dimension of the ultrasound receive transducer aperture (dL). The aim of the paper was to explore these temporal and spatial dependencies through a comparison of experimental measurement and computer simulation in solid:liquid models of varying temporal and spatial complexity. Transmission measurements at nominal ultrasound frequencies of 1MHz and 5MHz were performed, thereby investigating the dependency upon period. The results demonstrated an overall agreement between experimental measurement and computer simulation of 87±16% and 85±12% for temporal and spatial components respectively. It is envisaged that a comprehensive understanding of ultrasound propagation through complex structures such as cancellous bone could provide an improved non-invasive tool for osteoporosis assessment.

Conventional, Bayesian, and Modified Prony's methods for characterizing fast and slow waves in equine cancellous bone

Conventional, Bayesian, and the modified least-squares Prony's plus curve-fitting (MLSP + CF) methods were applied to data acquired using 1 MHz center frequency, broadband transducers on a single equine cancellous bone specimen that was systematically shortened from 11.8 mm down to 0.5 mm for a total of 24 sample thicknesses. Due to overlapping fast and slow waves, conventional analysis methods were restricted to data from sample thicknesses ranging from 11.8 mm to 6.0 mm. In contrast, Bayesian and MLSP + CF methods successfully separated fast and slow waves and provided reliable estimates of the ultrasonic properties of fast and slow waves for sample thicknesses ranging from 11.8 mm down to 3.5 mm. Comparisons of the three methods were carried out for phase velocity at the center frequency and the slope of the attenuation coefficient for the fast and slow waves. Good agreement among the three methods was also observed for average signal loss at the center frequency. The Bayesian and MLSP + CF approaches were able to separate the fast and slow waves and provide good estimates of the fast and slow wave properties even when the two wave modes overlapped in both time and frequency domains making conventional analysis methods unreliable.

Fast characterization of two ultrasound longitudinal waves in cancellous bone using an adaptive beamforming technique

The received signal in through-transmission ultrasound measurements of cancellous bone consists of two longitudinal waves, called the fast and slow waves. Analysis of these fast and slow waves may reveal characteristics of the cancellous bone that would be good indicators of osteoporosis. Because the two waves often overlap, decomposition of the received signal is an important problem in the characterization of bone quality. This study proposes a fast and accurate decomposition method based on the frequency domain interferometry imaging method with a modified wave transfer function that uses a phase rotation parameter. The proposed method accurately characterized the fast and slow waves in the experimental study, and the residual intensity, which was normalized with respect to the received signal intensity, was less than -20 dB over the bone specimen thickness range from 6 to 15 mm. In the simulation study, the residual intensity was less than -20 dB over the specimen thickness range from 3 to 8 mm. Decomposition of a single received signal takes only 5 s using a laptop personal computer with a single central processing unit. The proposed method has great potential to provide accurate and rapid measurements of indicators of osteoporosis in cancellous bone.

A deconvolution method for deriving the transit time spectrum for ultrasound propagation through cancellous bone replica models

The acceptance of broadband ultrasound attenuation for the assessment of osteoporosis suffers from a limited understanding of ultrasound wave propagation through cancellous bone. It has recently been proposed that the ultrasound wave propagation can be described by a concept of parallel sonic rays. This concept approximates the detected transmission signal to be the superposition of all sonic rays that travel directly from transmitting to receiving transducer. The transit time of each ray is defined by the proportion of bone and marrow propagated. An ultrasound transit time spectrum describes the proportion of sonic rays having a particular transit time, effectively describing lateral inhomogeneity of transit times over the surface of the receiving ultrasound transducer. The aim of this study was to provide a proof of concept that a transit time spectrum may be derived from digital deconvolution of input and output ultrasound signals. We have applied the active-set method deconvolution algorithm to determine the ultrasound transit time spectra in the three orthogonal directions of four cancellous bone replica samples and have compared experimental data with the prediction from the computer simulation. The agreement between experimental and predicted ultrasound transit time spectrum analyses derived from Bland–Altman analysis ranged from 92% to 99%, thereby supporting the concept of parallel sonic rays for ultrasound propagation in cancellous bone. In addition to further validation of the parallel sonic ray concept, this technique offers the opportunity to consider quantitative characterisation of the material and structural properties of cancellous bone, not previously available utilising ultrasound.

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.

Transient ultrasound propagation in porous media using Biot theory and fractional calculus: application to human cancellous bone

A temporal model based on the Biot theory is developed to describe the transient ultrasonic propagation in porous media with elastic structure, in which the viscous exchange between fluid and structure are described by fractional derivatives. The fast and slow waves obey a fractional wave equation in the time domain. The solution of Biot's equations in time depends on the Green functions of each of the waves (fast and slow), and their fractional derivatives. The reflection and transmission operators for a slab of porous materials are derived in the time domain, using calculations in the Laplace domain. Their analytical expressions, depend on Green's function of fast and slow waves. Experimental results for slow and fast waves transmitted through human cancellous bone samples are given and compared with theoretical predictions.

Mechanical competence of bone: a new parameter to grade trabecular bone fragility from tortuosity and elasticity

With the elderly population increase, osteoporosis and its consequences have become not just a health issue but also a serious economic burden. The trabecular bone structure plays a very important role for the bone quality and mechanical competence of the scaffold. Currently, it is claimed that the trabecular microarchitecture understanding can improve the fracture risk prediction above 65%. Several parameters seem to be correlated providing structural details of the trabecular bone network. However, the tortuosity of the trabeculae has not yet been systematically taken into account and its contribution has not been fully investigated and understood. In this paper, we discuss the relationship between the trabecular tortuosity, connectivity, volume fraction, and elasticity, and provide a unified parameter to estimate the mechanical competence of the structure. It is shown that the trabecular network tortuosity presents high linear correlation with the other parameters and that the trabeculae tend to get aligned in the direction where the structure is mostly submitted to stress, corresponding to higher stiffness orientation. This new parameter will help to integrate the relevant information of bone microarchitecture quality and assess more directly the real trabecular fragility in osteoporotic patients.

Cancellous bone fast and slow waves obtained with Bayesian probability theory correlate with porosity from computed tomography

A Bayesian probability theory approach for separating overlapping ultrasonic fast and slow waves in cancellous bone has been previously introduced. The goals of this study were to investigate whether the fast and slow waves obtained from Bayesian separation of an apparently single mode signal individually correlate with porosity and to isolate the fast and slow waves from medial-lateral insonification of the calcaneus. The Bayesian technique was applied to trabecular bone data from eight human calcanei insonified in the medial-lateral direction. The phase velocity, slope of attenuation (nBUA), and amplitude were determined for both the fast and slow waves. The porosity was assessed by micro-computed tomography (microCT) and ranged from 78.7% to 94.1%. The method successfully separated the fast and slow waves from medial-lateral insonification of the calcaneus. The phase velocity for both the fast and slow wave modes showed an inverse correlation with porosity (R(2) = 0.73 and R(2) = 0.86, respectively). The slope of attenuation for both wave modes also had a negative correlation with porosity (fast wave: R(2) = 0.73, slow wave: R(2) = 0.53). The fast wave amplitude decreased with increasing porosity (R(2) = 0.66). Conversely, the slow wave amplitude modestly increased with increasing porosity (R(2) = 0.39).

Relationships of quantitative ultrasound parameters with cancellous bone microstructure in human calcaneus in vitro

Ultrasound parameters (attenuation, phase velocity, and backscatter), bone mineral density (BMD), and microarchitectural features were measured on 29 human cancellous calcaneus samples in vitro. Regression analysis was performed to predict ultrasound parameters from BMD and microarchitectural features. The best univariate predictors of the ultrasound parameters were the indexes of bone quantity: BMD and bone volume fraction (BV/TV). The most predictive univariate models for attenuation, phase velocity, and backscatter coefficient yielded adjusted squared correlation coefficients of 0.69-0.73. Multiple regression models yielded adjusted correlation coefficients of 0.74-0.83. Therefore attenuation, phase velocity, and backscatter are primarily determined by bone quantity, but multiple regression models based on bone quantity plus microarchitectural features achieve slightly better predictive performance than models based on bone quantity alone.

The 25th Anniversary of BUA for the Assessment of Osteoporosis: Time for a New Paradigm?

The measurement of broadband ultrasonic attenuation (BUA) in cancellous bone at the calcaneus for the assessment of osteoporosis was first described within this journal 25 years ago. It was recognized in 2006 by Universities UK as being one of the ‘100 discoveries and developments in UK Universities that have changed the world’ over the past 50 years. In 2008, the UK's Department of Health also recognized BUA assessment of osteoporosis in a publication highlighting 11 projects that have contributed to ‘60 years of NHS research benefiting patients’. The BUA technique has been extensively clinically validated and is utilized worldwide, with at least seven commercial systems currently providing calcaneal BUA measurement. However, there is still no fundamental understanding of the dependence of BUA upon the material and structural properties of cancellous bone. This review aims to provide an ‘engineering in medicine’ perspective and proposes a new paradigm based upon phase cancellation due to variation in propagation transit time across the receive transducer face to explain the non-linear relationship between BUA and bone volume fraction in cancellous bone.

Influence of cancellous bone microstructure on two ultrasonic wave propagations in bovine femur: an in vitro study

The influence of cancellous bone microstructure on the ultrasonic wave propagation of fast and slow waves was experimentally investigated. Four spherical cancellous bone specimens extracted from two bovine femora were prepared for the estimation of acoustical and structural anisotropies of cancellous bone. In vitro measurements were performed using a PVDF transducer (excited by a single sinusoidal wave at 1 MHz) by rotating the spherical specimens. In addition, the mean intercept length (MIL) and bone volume fraction (BV/TV) were estimated by X-ray micro-computed tomography. Separation of the fast and slow waves was clearly observed in two specimens. The fast wave speed was strongly dependent on the wave propagation direction, with the maximum speed along the main trabecular direction. The fast wave speed increased with the MIL. The slow wave speed, however, was almost constant. The fast wave speeds were statistically higher, and their amplitudes were statistically lower in the case of wave separation than in that of wave overlap.

Cancellous bone analysis with modified least squares Prony's method and chirp filter: phantom experiments and simulation

The presence of two longitudinal waves in porous media is predicted by Biot's theory and has been confirmed experimentally in cancellous bone. When cancellous bone samples are interrogated in through-transmission, these two waves can overlap in time. Previously, the Modified Least-Squares Prony's (MLSP) method was validated for estimation of amplitudes, attenuation coefficients, and phase velocities of fast and slow waves, but tended to overestimate phase velocities by up to about 5%. In the present paper, a pre-processing chirp filter to mitigate the phase velocity bias is derived. The MLSP/chirp filter (MLSPCF) method was tested for decomposition of a 500 kHz-center-frequency signal containing two overlapping components: one passing through a low-density-polyethylene plate (fast wave) and another passing through a cancellous-bone-mimicking phantom material (slow wave). The chirp filter reduced phase velocity bias from 100 m/s (5.1%) to 69 m/s (3.5%) (fast wave) and from 29 m/s (1.9%) to 10 m/s (0.7%) (slow wave). Similar improvements were found for 1) measurements in polycarbonate (fast wave) and a cancellous-bone-mimicking phantom (slow wave), and 2) a simulation based on parameters mimicking bovine cancellous bone. The MLSPCF method did not offer consistent improvement in estimates of attenuation coefficient or amplitude.

Ultrasonic wave propagation in stereo-lithographical bone replicas

Predictions of a modified anisotropic Biot-Allard theory are compared with measurements of pulses centered on 100 kHz and 1 MHz transmitted through water-saturated stereo-lithographical bone replicas. The replicas are 13 times larger than the original bone samples. Despite the expected effects of scattering, which is neglected in the theory, at 100 kHz the predicted and measured transmitted waveforms are similar. However, the magnitude of the leading negative edge of the waveform is overpredicted, and the trailing parts of the waveforms are not predicted well. At 1 MHz, although there are differences in amplitudes, the theory predicts that the transmitted waveform is almost a scaled version of that incident in conformity with the data.

Decomposition of two-component ultrasound pulses in cancellous bone using modified least squares prony method--phantom experiment and simulation

Porous media such as cancellous bone often support the simultaneous propagation of two compressional waves. When small bone samples are interrogated in through-transmission with broadband sources, these two waves often overlap in time. The modified least-squares Prony's (MLSP) method was tested for decomposing a 500 kHz-center-frequency signal containing two overlapping components: one passing through a polycarbonate plate (to produce the "fast" wave) and another passing through a cancellous-bone-mimicking phantom (to produce the "slow" wave). The MLSP method yielded estimates of attenuation slopes accurate to within 7% (polycarbonate plate) and 2% (cancellous bone phantom). The MLSP method yielded estimates of phase velocities accurate to within 1.5% (both media). The MLSP method was also tested on simulated data generated using attenuation slopes and phase velocities corresponding to bovine cancellous bone. Throughout broad ranges of signal-to-noise ratio (SNR), the MLSP method yielded estimates of attenuation slope that were accurate to within 1.0% and estimates of phase velocity that were accurate to within 4.3% (fast wave) and 1.3% (slow wave).

Predictions of angle dependent tortuosity and elasticity effects on sound propagation in cancellous bone

The anisotropic pore structure and elasticity of cancellous bone cause wave speeds and attenuation in cancellous bone to vary with angle. Previously published predictions of the variation in wave speed with angle are reviewed. Predictions that allow tortuosity to be angle dependent but assume isotropic elasticity compare well with available data on wave speeds at large angles but less well for small angles near the normal to the trabeculae. Claims for predictions that only include angle-dependence in elasticity are found to be misleading. Audio-frequency data obtained at audio-frequencies in air-filled bone replicas are used to derive an empirical expression for the angle-and porosity-dependence of tortuosity. Predictions that allow for either angle dependent tortuosity or angle dependent elasticity or both are compared with existing data for all angles and porosities.

Propagation of two longitudinal waves in human cancellous bone: an in vitro study

The ultrasonic wave propagation of fast and slow waves was investigated in vitro in 35 cubic cancellous bone specimens extracted from human femoral heads. Measurements were performed in three orthogonal directions using home-made PVDF transducers excited by a single sinusoidal wave at 1 MHz. The apparent density of the specimens was measured. Two separated fast and slow waves were clearly observed in 16 specimens, mainly in the main load direction. The waveforms and the sound speeds of fast and slow waves were similar to the reported data in bovine bone. The group of specimens in which the two waves were observed did not exhibit statistically higher apparent density than the rest of the specimens, but did exhibit statistically higher acoustic anisotropy ratio. The speeds in the main load direction were higher than those in the other direction. The fast and slow wave speeds were in good agreement with Biot's model, showing an increase with bone volume fraction (BV/TV). The ratio of peak amplitudes of the fast and slow waves nonlinearly increased as a function of BV/TV. These results open interesting perspective for acoustic assessment of cancellous bone micro-architecture and especially anisotropy that might lead to an improved assessment of bone strength.

Analysis of frequency dependence of ultrasonic backscatter coefficient in cancellous bone

The ultrasonic scattering mechanism in cancellous bone is investigated theoretically and a model describing the frequency dependence of ultrasonic scattering from cancellous bone is presented. The ultrasonic backscatter coefficient (BSC) of bovine tibiae, human calcanei in vitro and in vivo, were measured and discussed. The data of BSC were also fitted by polynomial. The results demonstrate that BSC is a nonlinear function of frequency and increases with frequency. A good agreement was obtained between BSC values from theory and experiment. Also, the high correlation coefficient between BSC and bone mineral density was obtained, r=0.85+/-0.07 (mean+/-SD) (n=15, p<0.001). Based on the values of BSC, the status of cancellous bone and the degree of osteoporotic fracture risk may be assessed.

Interference between wave modes may contribute to the apparent negative dispersion observed in cancellous bone

Previous work has shown that ultrasonic waves propagating through cancellous bone often exhibit a linear-with-frequency attenuation coefficient, but a decrease in phase velocity with frequency (negative dispersion) that is inconsistent with the causality-imposed Kramers-Kronig relations. In the current study, interfering wave modes similar to those observed in bone are shown to potentially contribute to the observed negative dispersion. Biot theory, the modified Biot-Attenborogh model, and experimental results are used to aid in simulating multiple-mode wave propagation through cancellous bone. Simulations entail constructing individual wave modes exhibiting a positive dispersion using plausible velocities and amplitudes, and then summing the individual modes to create mixed-mode output wave forms. Results of the simulations indicate that mixed-mode wave forms can exhibit negative dispersion when analyzed conventionally under the assumption that only one wave is present, even when the individual interfering waves exhibit positive dispersions in accordance with the Kramers-Kronig relations. Furthermore, negative dispersion is observed when little or no visual evidence of interference exists in the time-domain data. Understanding the mechanisms responsible for the observed negative dispersion could aid in determining the true material properties of cancellous bone, as opposed to the apparent properties measured using conventional data analysis techniques.

Effects of structural anisotropy of cancellous bone on speed of ultrasonic fast waves in the bovine femur

Ultrasonic waves in cancellous bone change dramatically depending on its structural complexity. One good example is the separation of an ultrasonic longitudinal wave into fast and slow waves during propagation. In this study, we examined fast wave propagation in cancellous bone obtained from the head of the bovine femur, taking the bone structure into consideration. We investigated the wave propagation perpendicular to the bone axis and found the two-wave phenomenon. By rotating the cylindrical cancellous bone specimen, changes in the fast wave speed due to the rotation angle then were observed. In addition to the ultrasonic evaluation, the structural anisotropy of each specimen was measured by X-ray micro-computed tomography (CT). From the CT images, we obtained the mean intercept length (MIL), degree of anisotropy (DA), and angle of insonification relative to the trabecular orientation. The ultrasonic and CT results showed that the fast wave speed was dependent on the structural anisotropy, especially on the trabecular orientation and length. The fast wave speeds always were higher for propagation parallel to the trabecular orientation. In addition, there was a strong correlation between the DA and the ratio between maximum and minimum speeds (V(max)/V(min)) (R(2) = 0.63).

Simplified inverse filter tracking algorithm for estimating the mean trabecular bone spacing

Ultrasonic backscatter signals provide useful information relevant to bone tissue characterization. Trabecular bone microstructures have been considered as quasi-periodic tissues with a collection of regular and diffuse scatterers. This paper investigates the potential of a novel technique using a simplified inverse filter tracking (SIFT) algorithm to estimate mean trabecular bone spacing (MTBS) from ultrasonic backscatter signals. In contrast to other frequency-based methods, the SIFT algorithm is a time-based method and utilizes the amplitude and phase information of backscatter echoes, thus retaining the advantages of both the autocorrelation and the cepstral analysis techniques. The SIFT algorithm was applied to backscatter signals from simulations, phantoms, and bovine trabeculae in vitro. The estimated MTBS results were compared with those of the autoregressive (AR) cepstrum and quadratic transformation (QT) . The SIFT estimates are better than the AR cepstrum estimates and are comparable with the QT values. The study demonstrates that the SIFT algorithm has the potential to be a reliable and robust method for the estimation of MTBS in the presence of a small signal-to-noise ratio, a large spacing variation between regular scatterers, and a large scattering strength ratio of diffuse scatterers to regular ones.

The measurement of broadband ultrasonic attenuation in cancellous bone--a review of the science and technology

The measurement of broadband ultrasonic attenuation (BUA) in cancellous bone at the calcaneus was first described in 1984. The assessment of osteoporosis by BUA has recently been recognized by Universities UK, within its EurekaUK book, as being one of the "100 discoveries and developments in UK Universities that have changed the world" over the past 50 years, covering the whole academic spectrum from the arts and humanities to science and technology. Indeed, BUA technique has been clinically validated and is utilized worldwide, with at least seven commercial systems providing calcaneal BUA measurement. However, a fundamental understanding of the dependence of BUA upon the material and structural properties of cancellous bone is still lacking. This review aims to provide a science- and technology-orientated perspective on the application of BUA to the medical disease of osteoporosis.

Application of Biot's theory to ultrasonic characterization of human cancellous bones: determination of structural, material, and mechanical properties

This paper is devoted to the experimental determination of distinctive macroscopic structural (porosity, tortuosity, and permeability) and mechanical (Biot-Willis elastic constants) properties of human trabecular bones. Then, the obtained data may serve as input parameters for modeling wave propagation in cancellous bones using Biot's theory. The goal of the study was to obtain experimentally those characteristics for statistically representative group of human bones (35 specimens) obtained from a single skeletal site (proximal femur). The structural parameters were determined using techniques devoted to the characterization of porous materials: electrical spectroscopy, water permeametry, and microcomputer tomography. The macroscopic mechanical properties, Biot-Willis elastic constants, were derived based on the theoretical consideration of Biot's theory, micromechanical statistical models, and experimental results of ultrasonic studies for unsaturated cancellous bones. Our results concerning structural parameters are consistent with the data presented by the other authors, while macroscopic mechanical properties measured within our studies are situated between the other published data. The discrepancies are mainly attributed to different mechanical properties of the skeleton frame, due to strong structural anisotropy varying from site to site. The results enlighten the difficulty to use Biot's theory for modeling wave propagation in cancellous bone, implying necessity of individual evaluation of input parameters.