A stratified model to predict dispersion in trabecular bone

A generalized fractional-order elastodynamic theory for non-local attenuating media

This study presents a generalized elastodynamic theory, based on fractional-order operators, capable of modelling the propagation of elastic waves in non-local attenuating solids and across complex non-local interfaces. Classical elastodynamics cannot capture hybrid field transport processes that are characterized by simultaneous propagation and diffusion. The proposed continuum mechanics formulation, which combines fractional operators in both time and space, offers unparalleled capabilities to predict the most diverse combinations of multiscale, non-local, dissipative and attenuating elastic energy transport mechanisms. Despite the many features of this theory and the broad range of applications, this work focuses on the behaviour and modelling capabilities of the space-fractional term and on its effect on the elastodynamics of solids. We also derive a generalized fractional-order version of Snell's Law of refraction and of the corresponding Fresnel's coefficients. This formulation allows predicting the behaviour of fully coupled elastic waves interacting with non-local interfaces. The theoretical results are validated via direct numerical simulations.

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.

Artificial neural network to estimate micro-architectural properties of cortical bone using ultrasonic attenuation: A 2-D numerical study

The goal of this study is to estimate micro-architectural parameters of cortical porosity such as pore diameter (φ), pore density (ρ) and porosity (ν) of cortical bone from ultrasound frequency dependent attenuation using an artificial neural network (ANN). First, heterogeneous structures with controlled pore diameters and pore densities (mono-disperse) were generated, to mimic simplified structure of cortical bone. Then, more realistic structures were obtained from high resolution CT scans of human cortical bone. 2-D finite-difference time-domain simulations were conducted to calculate the frequency-dependent attenuation in the 1-8 MHz range. An ANN was then trained with the ultrasonic attenuation at different frequencies as the input feature vectors while the output was set as the micro-architectural parameters (pore diameter, pore density and porosity). The ANN is composed of three fully connected dense layers with 24, 12 and 6 neurons, connected to the output layer. The dataset was trained over 6000 epochs with a batch size of 16. The trained ANN exhibits the ability to predict the micro-architectural parameters with high accuracy and low losses. ANN approaches could potentially be used as a tool to help inform physics-based modelling of ultrasound propagation in complex media such as cortical bone. This will lead to the solution of inverse-problems to retrieve bone micro-architectural parameters from ultrasound measurements for the non-invasive diagnosis and monitoring osteoporosis.

Velocity dispersion and backscatter in marrow-filled and water-filled trabecular bone samples in vitro

The phase velocity and the backscatter coefficient were measured in 28 bovine femoral trabecular bone samples filled with marrow and water in vitro from 0.2 to 0.6 MHz. The phase velocities decreased approximately linearly with increasing frequency and the average dispersion rate of -34 ms^-1MHz^-1 in the marrow-filled samples was higher than that of -42 ms^-1MHz^-1 in the water-filled samples. The backscatter coefficients exhibited nonlinear, monotonically increasing dependences on the frequency and the average value of the exponent n = 2.92 (frequency dependence) in the marrow-filled samples was higher than the value of n = 2.79 in the water-filled samples.

Micro-scale finite element modeling of ultrasound propagation in aluminum trabecular bone-mimicking phantoms: A comparison between numerical simulation and experimental results

The present study investigated the accuracy of micro-scale finite element modeling for simulating broadband ultrasound propagation in water-saturated trabecular bone-mimicking phantoms. To this end, five commercially manufactured aluminum foam samples as trabecular bone-mimicking phantoms were utilized for ultrasonic immersion through-transmission experiments. Based on micro-computed tomography images of the same physical samples, three-dimensional high-resolution computational samples were generated to be implemented in the micro-scale finite element models. The finite element models employed the standard Galerkin finite element method (FEM) in time domain to simulate the ultrasonic experiments. The numerical simulations did not include energy dissipative mechanisms of ultrasonic attenuation; however, they expectedly simulated reflection, refraction, scattering, and wave mode conversion. The accuracy of the finite element simulations were evaluated by comparing the simulated ultrasonic attenuation and velocity with the experimental data. The maximum and the average relative errors between the experimental and simulated attenuation coefficients in the frequency range of 0.6-1.4 MHz were 17% and 6% respectively. Moreover, the simulations closely predicted the time-of-flight based velocities and the phase velocities of ultrasound with maximum relative errors of 20 m/s and 11 m/s respectively. The results of this study strongly suggest that micro-scale finite element modeling can effectively simulate broadband ultrasound propagation in water-saturated trabecular bone-mimicking structures.

Nonlinear attenuation and dispersion in human calcaneus in vitro: statistical validation and relationships to microarchitecture

Through-transmission measurements were performed on 30 human calcaneus samples in vitro. Nonlinear attenuation and dispersion measurements were investigated by estimating 95% confidence intervals of coefficients of polynomial expansions of log magnitude and phase of transmission coefficients. Bone mineral density (BMD) was measured with dual x-ray absorptiometry. Microarchitecture was measured with microcomputed tomography. Statistically significant nonlinear attenuation and nonzero dispersion were confirmed for a clinical bandwidth of 300-750 kHz in 40%-43% of bone samples. The mean linear coefficient for attenuation was 10.3 dB/cm MHz [95% confidence interval (CI): 9.0-11.6 dB/cm MHz]. The mean quadratic coefficient for attenuation was 1.6 dB/cm MHz(2) (95% CI: 0.4-2.8 dB/cm MHz(2)). Nonlinear attenuation provided little information regarding BMD or microarchitecture. The quadratic coefficient for phase (which is related to dispersion) showed moderate correlations with BMD (r = -0.65; 95% CI: -0.82 to -0.36), bone surface-to-volume ratio (r = 0.47; 95% CI: 0.12-0.72) and trabecular thickness (r = -0.40; 95% CI: -0.67 to -0.03). Dispersion was proportional to bone volume fraction raised to an exponent of 2.1 ± 0.2, which is similar to the value for parallel nylon-wire phantoms (2.4 ± 0.2) and supports a multiple-scattering model for dispersion.

Attenuated Fractional Wave Equations With Anisotropy

This paper develops new fractional calculus models for wave propagation. These models permit a different attenuation index in each coordinate to fully capture the anisotropic nature of wave propagation in complex media. Analytical expressions that describe power law attenuation and anomalous dispersion in each direction are derived for these fractional calculus models.

The finite element method for micro-scale modeling of ultrasound propagation in cancellous bone

Quantitative ultrasound for bone assessment is based on the correlations between ultrasonic parameters and the properties (mechanical and physical) of cancellous bone. To elucidate the correlations, understanding the physics of ultrasound in cancellous bone is demanded. Micro-scale modeling of ultrasound propagation in cancellous bone using the finite-difference time-domain (FDTD) method has been so far utilized as one of the approaches in this regard. However, the FDTD method accompanies two disadvantages: staircase sampling of cancellous bone by finite difference grids leads to generation of wave artifacts at the solid-fluid interface inside the bone; additionally, this method cannot explicitly satisfy the needed perfect-slip conditions at the interface. To overcome these disadvantages, the finite element method (FEM) is proposed in this study. Three-dimensional finite element models of six water-saturated cancellous bone samples with different bone volume were created. The values of speed of sound (SOS) and broadband ultrasound attenuation (BUA) were calculated through the finite element simulations of ultrasound propagation in each sample. Comparing the results with other experimental and simulation studies demonstrated the capabilities of the FEM for micro-scale modeling of ultrasound in water-saturated cancellous bone.

Multichannel instantaneous frequency analysis of ultrasound propagating in cancellous bone

An ultrasonic pulse propagating in cancellous bone can be separated into two waves depending on the condition of the specimen. These two waves, which are called the fast wave and the slow wave, provide important information for the diagnosis of osteoporosis. The present study proposes to utilize a signal processing method that extracts the instantaneous frequency (IF) of waveforms from multiple spectral channels. The instantaneous frequency was expected to be able to show detailed time-frequency properties of ultrasonic waves being transmitted through cancellous bone. The employed method, termed the multichannel instantaneous frequency (MCIF) method, showed robustness against background noise as compared to the IF that was directly derived from the original waveform. The extracted IF revealed that the frequency of the fast wave was affected by both the propagation distance within the specimen and the bone density, independently. On the other hand, the alternation of the center frequency of the originally transmitted wave did not produce proportional changes in the extracted IF values of the fast waves, suggesting that the fast wave IF mainly reflected the thickness of the specimens. These findings may provide the possibility of obtaining a more precise diagnosis of osteoporosis.

In vitro estimation of fast and slow wave parameters of thin trabecular bone using space-alternating generalized expectation-maximization algorithm

In testing cancellous bone using ultrasound, two types of longitudinal Biot's waves are observed in the received signal. These are known as fast and slow waves and their appearance depend on the alignment of bone trabeculae in the propagation path and the thickness of the specimen under test (SUT). They can be used as an effective tool for the diagnosis of osteoporosis because wave propagation behavior depends on the bone structure. However, the identification of these waves in the received signal can be difficult to achieve. In this study, ultrasonic wave propagation in a 4mm thick bovine cancellous bone in the direction parallel to the trabecular alignment is considered. The observed Biot's fast and slow longitudinal waves are superimposed; which makes it difficult to extract any information from the received signal. These two waves can be separated using the space alternating generalized expectation maximization (SAGE) algorithm. The latter has been used mainly in speech processing. In this new approach, parameters such as, arrival time, center frequency, bandwidth, amplitude, phase and velocity of each wave are estimated. The B-Scan images and its associated A-scans obtained through simulations using Biot's finite-difference time-domain (FDTD) method are validated experimentally using a thin bone sample obtained from the femoral-head of a 30 months old bovine.

Correlations of group velocity, phase velocity, and dispersion with bone density in bovine trabecular bone

The present study investigated the correlations of the group velocity, the phase velocity, and the velocity dispersion with the apparent bone density in bovine trabecular bone in vitro. The phase velocity exhibited the negative dispersion, consistent with the behavior in human trabecular bone. The group and the phase velocities were found to increase with increasing apparent bone density, respectively, exhibiting similar high correlations of r=0.94 and 0.96. The negative dispersion rate exhibited a decreasing dependence on the apparent bone density, with a significant correlation of r=-0.86.

Determining attenuation properties of interfering fast and slow ultrasonic waves in cancellous bone

Previous studies have shown that interference between fast waves and slow waves can lead to observed negative dispersion in cancellous bone. In this study, the effects of overlapping fast and slow waves on measurements of the apparent attenuation as a function of propagation distance are investigated along with methods of analysis used to determine the attenuation properties. Two methods are applied to simulated data that were generated based on experimentally acquired signals taken from a bovine specimen. The first method uses a time-domain approach that was dictated by constraints imposed by the partial overlap of fast and slow waves. The second method uses a frequency-domain log-spectral subtraction technique on the separated fast and slow waves. Applying the time-domain analysis to the broadband data yields apparent attenuation behavior that is larger in the early stages of propagation and decreases as the wave travels deeper. In contrast, performing frequency-domain analysis on the separated fast waves and slow waves results in attenuation coefficients that are independent of propagation distance. Results suggest that features arising from the analysis of overlapping two-mode data may represent an alternate explanation for the previously reported apparent dependence on propagation distance of the attenuation coefficient of cancellous bone.

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.

Spatial backward planar projection in absorbing media possessing an arbitrary dispersion relation

Planar projection methods have been shown to rapidly relate fields between two planes. Such an approach is particularly useful for characterizing transducers, since only a single plane needs to be measured in order to characterize an entire field. The present work considers the same approach in the presence of an arbitrary dispersion relation. Unlike traditional methods that use Fourier solutions of the time-domain wave equation, the approach starts from a frequency-domain Helmholtz equation for waves in a dispersive medium. It is shown that a transfer function similar to that derived from time domain equations can be utilized. Both the forward- and backward-projection behaviors are examined and it is demonstrated that the approach is invariant to propagation direction.

Frequency dependence of average phase shift from human calcaneus in vitro

If dispersion in a medium is weak and approximately linear with frequency (over the experimental band of frequencies), then it can be shown that the constant term in a polynomial representation of phase shift as a function of frequency can produce errors in measurements of phase-velocity differences in through-transmission, substitution experiments. A method for suppressing the effects of the constant phase shift in the context of the single-wave-model was tested on measurements from 30 cancellous human calcaneus samples in vitro. Without adjustment for constant phase shifts, the estimated phase velocity at 500 kHz was 1516+/-6 m/s (mean+/-standard error), and the estimated dispersion was -24+/-4 m/s MHz (mean+/-standard error). With adjustment for constant phase shifts, the estimated mean velocity decreased by 4-9 m/s, and the estimated magnitude of mean dispersion decreased by 50%-100%. The average correlation coefficient between the measured attenuation coefficient and frequency was 0.997+/-0.0026 (mean+/-standard deviation), suggesting that the signal for each sample was dominated by one wave. A single-wave, linearly dispersive model conformed to measured complex transfer functions from the 30 cancellous-bone samples with an average root-mean-square error of 1.9%+/-1.0%.

Imaging of gaps in digital joints by measurement of ultrasound transmission using a linear array

In orthodontic dentistry for young subjects, it is important to assess the degree of growth of the jaw bones to determine the optimum time for treatment. The structure of the digital joint changes with age, with such changes correlating to the degree of bone growth (including jaw bones). There are two gaps in the digital joint of a young subject, one of which disappears with aging. In the present study, a method for noninvasive assessment of such change in the structure of a digital joint was examined, in which continuous-wave ultrasound is radiated to a digital joint by a single-element ultrasonic transducer. This continuous ultrasound, which passes through the digital joint, is received by a linear array ultrasonic probe situated opposite the transducer. The probe simultaneously realizes pulse-echo imaging and imaging of transmission ultrasound, which passes through the joint. Using this experimental apparatus, the existence and position of a gap can be detected clearly by imaging the transmission ultrasound on a pulse-echo image. In basic experiments, continuous-wave ultrasound generated by a planar or focused transducer was radiated to a gap between two acrylic bars, which simulated that in a digital joint; transmission ultrasound, which passed through the gap, was measured with a linear array probe. The basic experimental results showed that a gap with a width >0.4 mm is detectable and that the width at half maximum of the amplitude profile of the received transmission ultrasound that passed through the gap correlated with the width of the gap. Furthermore, in the preliminary in vivo experiments, transmission ultrasound that passed through two gaps in the case of a child was clearly imaged by the proposed method, and that which passed through only one gap in the case of an adult was also imaged. These results show the possibility for the use of the proposed method to noninvasively assess the change in the structure of a joint as a result of aging.

The dependencies of phase velocity and dispersion on volume fraction in cancellous-bone-mimicking phantoms

Frequency-dependent phase velocity was measured in eight cancellous-bone-mimicking phantoms consisting of suspensions of randomly oriented nylon filaments (simulating trabeculae) in a soft-tissue-mimicking medium (simulating marrow). Trabecular thicknesses ranged from 152 to 356 mum. Volume fractions of nylon filament material ranged from 0% to 10%. Phase velocity varied approximately linearly with frequency over the range from 300 to 700 kHz. The increase in phase velocity (compared with phase velocity in a phantom containing no filaments) at 500 kHz was approximately proportional to volume fraction occupied by nylon filaments. The derivative of phase velocity with respect to frequency was negative and exhibited nonlinear, monotonically decreasing dependence on volume fraction. The dependencies of phase velocity and its derivative on volume fraction in these phantoms were similar to those reported in previous studies on (1) human cancellous bone and (2) phantoms consisting of parallel nylon wires immersed in water.

Velocity dispersion in trabecular bone: influence of multiple scattering and of absorption

Speed of sound measurements are widely used clinically to assess bone strength. Trabecular bone is an attenuating composite material in which negative values of velocity dispersion have been measured, this behavior remaining poorly explained physically. The aim of this work is to describe the ultrasonic propagation in trabecular bone modeled by infinite cylinders immersed in a saturating matrix, and to derive the physical determinants of velocity dispersion. A homogenization model accounting for the coupling of multiple scattering and absorption phenomena allows the computation of phase velocity and of dispersion while varying bone properties. The present model is adapted from the generalized self-consistent method developed in the work of Yang and Mal [(1994). "Multiple-scattering of elastic waves in a fiber-reinforced composite," J. Mech. Phys. Solids 42, 1945-1968]. It predicts negative values of velocity dispersion, in agreement with experimental results obtained in phantoms mimicking trabecular bone. In trabecular bone, mostly negative and also positive values of velocity dispersion are predicted by the model, which span within the range of values measured experimentally. Scattering effects are responsible for the negative values of dispersion, whereas the frequency dependence of the absorption coefficient in bone marrow and/or in the trabeculae results in an increase in dispersion, which may then become positive.

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.

Ultrasonic velocity dispersion in bovine cortical bone: an experimental study

Cortical bone quality is determinant in bone fragility and its ultrasonic evaluation has become possible in clinical practice. However, the interaction between a broadband ultrasonic pulse and this complex multiscale medium remains poorly understood. The frequency dependence of phase velocity, which may impact clinical measurements, has been sparsely investigated. Our objective is to evaluate the determinants of the frequency dependence of phase velocity in bovine femoral cortical bone samples using an in vitro ultrasonic transmission device. The apparent phase velocity varies quasilinearly on a 1 MHz restricted bandwidth around 4 MHz. After compensating for diffraction effects, significant differences in velocity dispersion are obtained according to the anatomical location. The microstructure of each sample is determined using an optical microscope, which allows assessing the dependence of dispersion on the type of bone microstructure. Mostly positive but also negative values of dispersion are measured. Negative dispersion is mainly obtained in samples constituted of mixed microstructure, which may be explained by phase cancellation effects due to the presence of different microstructures within the same sample. Dispersion is shown to be related to broadband ultrasonic attenuation values, especially in the radial direction. Results are compared with the local Kramers-Kronig relationships.

Negative dispersion in bone: the role of interference in measurements of the apparent phase velocity of two temporally overlapping signals

In this study the attenuation coefficient and dispersion (frequency dependence of phase velocity) are measured using a phase sensitive (piezoelectric) receiver in a phantom in which two temporally overlapping signals are detected, analogous to the fast and slow waves typically found in measurements of cancellous bone. The phantom consisted of a flat and parallel Plexiglas plate into which a step discontinuity was milled. The phase velocity and attenuation coefficient of the plate were measured using both broadband and narrowband data and were calculated using standard magnitude and phase spectroscopy techniques. The observed frequency dependence of the phase velocity and attenuation coefficient exhibit significant changes in their frequency dependences as the interrogating ultrasonic field is translated across the step discontinuity of the plate. Negative dispersion is observed at specific spatial locations of the plate at which the attenuation coefficient rises linearly with frequency, a behavior analogous to that of bone measurements reported in the literature. For all sites investigated, broadband and narrowband data (3-7 MHz) demonstrate excellent consistency. Evidence suggests that the interference between the two signals simultaneously reaching the phase sensitive piezoelectric receiver is responsible for this negative dispersion.

Prediction of negative dispersion by a nonlocal poroelastic theory

The objective of this work is to show that the negative dispersion of ultrasonic waves propagating in cancellous bone can be explained by a nonlocal version of Biot's theory of poroelasticity. The nonlocal poroelastic formulation is presented in this work and the exact solutions for one- and two-dimensional systems are obtained by the method of Fourier transform. The nonlocal phase speeds for solid- and fluid-borne waves show the desired negative dispersion where the magnitude of dispersion is strongly dependent on the nonlocal parameters and porosity. Dependence of the phase speed and attenuation is studied for both porosity and frequency variation. It is shown that the nonlocal parameter can be easily estimated by comparing the theoretical dispersion rate with experimental observations. It is also shown that the modes of Lamb waves show similar negative dispersion when predicted by the nonlocal poroelastic theory.

Predictions of the modified Biot-Attenborough model for the dependence of phase velocity on porosity in cancellous bone

The modified Biot-Attenborough (MBA) model for acoustic wave propagation in porous media has been found useful to predict wave properties in cancellous bone. The present study is aimed at applying the MBA model to predict the dependence of phase velocity on porosity in cancellous bone. The MBA model predicts a phase velocity that decreases nonlinearly with porosity. The optimum values for input parameters of the MBA model, such as compressional speed c(m) of solid bone and phase velocity parameter s(2), were determined by comparing the predictions with previously published measurements in human calcaneus and bovine cancellous bone. The value of the phase velocity parameter s(2)=1.23 was obtained by curve fitting to the experimental data for 53 human calcaneus samples only, assuming a compressional speed c(m)=2500 m/s of solid bone. The root-mean-square error (RMSE) of the curve fit was 15.3m/s. The optimized value of s(2) for all 75 cancellous bone samples including 22 bovine samples was 1.42 with a value of 55 m/s for the RMSE of the curve fit. The latter fit was obtained by using of a value of c(m)=3200 m/s. Although the MBA model relies on the empirical parameters determined from experimental data, it is expected that the model can be usefully employed as a practical tool in the field of clinical ultrasonic bone assessment.

Phase velocity and normalized broadband ultrasonic attenuation in Polyacetal cuboid bone-mimicking phantoms

The variations of phase velocity and normalized broadband ultrasonic attenuation (nBUA) with porosity were investigated in Polyacetal cuboid bone-mimicking phantoms with circular cylindrical pores running normal to the surface along the three orthogonal axes. The frequency-dependent phase velocity and attenuation coefficient in the phantoms with porosities from 0% to 65.9% were measured from 0.65 to 1.10 MHz. The results showed that the phase velocity at 880 kHz decreased linearly with porosity, whereas the nBUA increased linearly with porosity. This study provides a useful insight into the relationships between ultrasonic properties and porosity in bone at porosities lower than 70%.

Group velocity, phase velocity, and dispersion in human calcaneus in vivo

Commercial bone sonometers measure broadband ultrasonic attenuation and/or speed of sound (SOS) in order to assess bone status. Phase velocity, which is usually measured in frequency domain, is a fundamental material property of bone that is related to SOS, which is usually measured in time domain. Four previous in vitro studies indicate that phase velocity in human cancellous bone decreases with frequency (i.e., negative dispersion). In order to investigate frequency-dependent phase velocity in vivo, through-transmission measurements were performed in 73 women using a GE Lunar Achilles Insight commercial bone sonometer. Average phase velocity at 500 kHz was 1489 +/- 55 m/s (mean +/- standard deviation). Average dispersion rate was -59 +/- 52 m/sMHz. Group velocity was usually lower than phase velocity, as is expected for negatively dispersive media. Using a stratified model to represent cancellous bone, the reductions in phase velocity and dispersion rate in vivo as opposed to in vitro can be explained by (1) the presence of marrow instead of water as a fluid filler, and (2) the decreased porosity of bones of living (compared with deceased) subjects.

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

The modeling of ultrasonic propagation in cancellous bone is relevant to the study of clinical bone assessment. Historical experiments revealed the importance of both the viscous effects of bone marrow and the anisotropy of the porous microstructure. Of those propagation models previously applied to cancellous bone, Biot's theory incorporates viscosity, but has only been applied in isotropic form, while Schoenberg's anisotropic model does not include viscosity. In this paper we present an approach that incorporates the merits of both models, by utilizing the tortuosity, a key parameter describing pore architecture. An angle-dependent tortuosity for a layered structure is used in Biot's theory to generate the "Stratified Biot Model" for cancellous bone, which is compared with published bone data. While the Stratified Biot model was inferior to Schoenberg's model for slow wave velocity prediction, the proposed model improved agreement fast wave velocity at high propagation angles, particularly when sorted for porosity. An attempt was made to improve the fast wave agreement at low angles by introducing an angle-dependent Young's Modulus, which, while improving the agreement of predicted fast wave velocity at low angles, degraded agreement at high angles. In this paper the utility of the tortuosity in characterizing the architecture of cancellous bone is highlighted.

Numerical simulation of the dependence of quantitative ultrasonic parameters on trabecular bone microarchitecture and elastic constants

Finite-difference numerical simulation of ultrasound propagation in complex media such as cancellous bone represents a fertile alternative to analytical approaches because it can manage the complex 3D bone structure by coupling the numerical computation with 3D numerical models of bone microarchitecture obtained from high-resolution imaging modalities. The objective of this work was to assess in silico the sensitivity of ultrasound parameters to controlled changes of microarchitecture and variation of elastic constants. The simulation software uses a finite-difference approach based on the Virieux numerical scheme. An incident plane wave was propagated through a volume of bone of approximately 5 x 5 x 8 mm(3). The volumes were reconstructed from high-resolution micro-computed tomography data. An iterative numerical scenario of "virtual osteoporosis" was implemented using a dedicated image processing algorithm in order to modify the initial 3D microstructures. Numerical computations of wave propagation were performed at each step of the process. The sensitivity to bone material properties was also tested by changing the elastic constants of bone tissue. Our results suggest that ultrasonic variables (slope of the frequency-dependent attenuation coefficient and speed of sound) are mostly influenced by bone volume fraction. However, material properties and structure also appear to play a role. The impact of modifications of the stiffness coefficients remained lower than the variability caused by structural variations. This study emphasizes the potential of numerical computations tools coupled to realistic 3D structures to elucidate the physical mechanisms of interaction between ultrasound and bone structure and to assess the sensitivity of ultrasound variables to different bone properties.

Anomalous negative dispersion in bone can result from the interference of fast and slow waves

The goal of this work was to show that the apparent negative dispersion of ultrasonic waves propagating in bone can arise from interference between fast and slow longitudinal modes, each exhibiting positive dispersion. Simulations were carried out using two approaches: one based on the Biot-Johnson model and one independent of that model. Results of the simulations are mutually consistent and appear to account for measurements from many laboratories that report that the phase velocity of ultrasonic waves propagating in cancellous bone decreases with increasing frequency (negative dispersion) in about 90% of specimens but increases with frequency in about 10%.

Effects of frequency-dependent attenuation and velocity dispersion on in vitro ultrasound velocity measurements in intact human femur specimens

Numerous studies have shown that ultrasonic velocity measured in bone provides a good assessment of osteoporotic fracture risk. However, a lack of standardization of signal processing techniques used to compute the speed of sound (SOS) complicates the comparison between data obtained with different commercial devices. In this study, 38 intact femurs were tested using a through-transmission technique and SOS determined using different techniques. The resulting difference in measured SOS was determined as functions of the attenuation and the velocity dispersion. A numerical simulation was used to explain how attenuation and dispersion impact two different SOS measurements (group velocity, velocity based on the first zero crossing of the signal). A new method aimed at compensating for attenuation was devised and led to a significant reduction in the difference between SOS obtained with both signal processing techniques. A comparison between SOS and X-ray density measurements indicated that the best correlation was reached for SOS based on the first zero crossing apparently because it used a marker located in the early part of the signal and was less sensitive to multipath interference. The conclusion is that first zero crossing velocity may be preferred to group velocity for ultrasonic assessment at this potential fracture site.

Three-dimensional simulation of ultrasound propagation through trabecular bone structures measured by synchrotron microtomography

Three-dimensional numerical simulations of ultrasound transmission were performed through 31 trabecular bone samples measured by synchrotron microtomography. The synchrotron microtomography provided high resolution 3D mappings of bone structures, which were used as the input geometry in the simulation software developed in our laboratory. While absorption (i.e. the absorption of ultrasound through dissipative mechanisms) was not taken into account in the algorithm, the simulations reproduced major phenomena observed in real through-transmission experiments in trabecular bone. The simulated attenuation (i.e. the decrease of the transmitted ultrasonic energy) varies linearly with frequency in the MHz frequency range. Both the speed of sound (SOS) and the slope of the normalized frequency-dependent attenuation (nBUA) increase with the bone volume fraction. Twenty-five out of the thirty-one samples exhibited negative velocity dispersion. One sample was rotated to align the main orientation of the trabecular structure with the direction of ultrasonic propagation, leading to the observation of a fast and a slow wave. Coupling numerical simulation with real bone architecture therefore provides a powerful tool to investigate the physics of ultrasound propagation in trabecular structures. As an illustration, comparison between results obtained on bone modelled either as a fluid or a solid structure suggested the major role of mode conversion of the incident acoustic wave to shear waves in bone to explain the large contribution of scattering to the overall attenuation.

The dependencies of phase velocity and dispersion on trabecular thickness and spacing in trabecular bone-mimicking phantoms

Frequency-dependent phase velocity was measured in trabecular-bone-mimicking phantoms consisting of two-dimensional arrays of parallel nylon wires (simulating trabeculae) with thicknesses ranging from 152 to 305 microm and spacings ranging from 700 to 1000 microm. Phase velocity varied approximately linearly with frequency over the range from 400 to 750 kHz. Dispersion was characterized by the slope of a linear least-squares regression fit to phase velocity versus frequency data. The increase in phase velocity (compared with that in water) at 500 kHz was approximately proportional to the (1) square of trabecular thickness, (2) inverse square of trabecular spacing, and (3) volume fraction occupied by nylon wires. The first derivative of phase velocity with respect to frequency was negative and exhibited nonlinear, monotonically decreasing dependencies on trabecular thickness and volume fraction. The dependencies of phase velocity and its first derivative on volume fraction in the phantoms were consistent with those reported in trabecular bone.

Comparison of measurements of phase velocity in human calcaneus to Biot theory

Biot's theory for elastic propagation in porous media has previously been shown to be useful for modeling the dependence of phase velocity on porosity in bovine cancellous bone in vitro. In the present study, Biot's theory is applied to measurements of porosity-dependent phase velocity in 53 human calcanea in vitro. Porosity was measured using microcomputed tomography for some samples (n = 23) and estimated based on bone mineral densitometry for the remaining samples (n = 30). The phase velocity at 500 kHz was measured in a water tank using a through-transmission technique. Biot's theory performed well for the prediction of the dependence of sound speed on porosity. The trend was quasilinear, but both the theory and experiment show similar slight curvature. The root mean square error (RMSE) of predicted versus measured sound speed was 15.8 m/s.

Comparison of acoustic characteristics predicted by Biot's theory and the modified Biot-Attenborough model in cancellous bone

Biot's theory and the modified Biot-Attenborough (MBA) model are applied to predict the dependences of acoustic characteristics on frequency and on porosity in cancellous bone. The phase velocities and the attenuation coefficients predicted by both theories are compared with the experimental data of bovine cancellous bone specimens published in the literature. Biot's theory successfully predicts the dependences of the phase velocity on frequency and on porosity in cancellous bone, whereas a significant discrepancy is observed between the predicted and the measured attenuation coefficients. The MBA model agrees well with the frequency and the porosity dependences of the phase velocity and the attenuation coefficient experimentally measured in bovine bones. Although the MBA model relies on phenomenological parameters derived from the experimental data, its approach to cancellous bone can be usefully employed in the field of clinical ultrasonic bone assessment.

Acoustic wave propagation in bovine cancellous bone: application of the Modified Biot-Attenborough model

Acoustic wave propagation in bovine cancellous bone is experimentally and theoretically investigated in the frequency range of 0.5-1 MHz. The phase velocity, attenuation coefficient, and broadband ultrasonic attenuation (BUA) of bovine cancellous bone are measured as functions of frequency and porosity. For theoretical estimation, the Modified Biot-Attenborough (MBA) model is employed with three new phenomenological parameters: the boundary condition, phase velocity, and impedance parameters. The MBA model is based on the idealization of cancellous bone as a nonrigid porous medium with circular cylindrical pores oriented normal to the surface. It is experimentally observed that the phase velocity is approximately nondispersive and the attenuation coefficient linearly increases with frequency. The MBA model predicts a slightly negative dispersion of phase velocity linearly with frequency and the nonlinear relationships of attenuation and BUA with porosity. The experimental results are in good agreement with the theoretical results estimated with the MBA model. It is expected that the MBA model can be usefully employed in the field of clinical bone assessment for the diagnosis of osteoporosis.