Fabric dependence of quasi-waves in anisotropic porous media

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.

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.

ULTRASONIC WAVE IS DETERMINED BY FABRIC TENSOR: AN APPLICATION TO CALCANEUS

The fabric tensor is a good measure for determining the mechanical properties of cancellous bone. Ultrasound is one method used to measure these mechanical properties. Ultrasound-generated speed of sound (SOS) measures the mechanical properties of cancellous bone. Thus, in this paper, we started with the fact that the fast wave in poroelastic theory is identical to the bulk wave velocity. We then formulate the equation for the fast wave in terms of fabric tensor for the calcaneus. The formulation in this paper is simpler than previously published results and will be easy to use in future experiments.

Enhanced correlation between quantitative ultrasound and structural and mechanical properties of bone using combined transmission-reflection measurement

Quantitative ultrasound (QUS) is capable of predicting the principal structural orientation of trabecular bone; this orientation is highly correlated with the mechanical strength of trabecular bone. Irregular shape of bone, however, would increase variation in such a prediction, especially under human in vivo measurement. This study was designed to combine transmission and reflection modes of QUS measurement to improve the prediction for the structural and mechanical properties of trabecular bone. QUS, mechanical testing, and micro computed tomography (μCT) scanning were performed on 24 trabecular bone cubes harvested from a bovine distal femur to obtain the mechanical and structural parameters. Transmission and reflection modes of QUS measurement in the transverse and frontal planes were performed in a confined 60° angle range with 5° increment. The QUS parameters, attenuation (ATT) and velocity (UV), obtained from transmission mode, were normalized to the specimen thickness acquired from reflection mode. Analysis of covariance showed that the combined transmission-reflection modes improved prediction for the structural and Young's modulus of bone in comparison to the traditional QUS measurement performed only in the medial-lateral orientation. In the transverse plane, significant improvement between QUS and μCT was found in ATT vs bone surface density (BS/BV) (p < 0.05), ATT vs trabecular thickness (Tb.Th) (p < 0.01), ATT vs degree of anisotropy (DA) (p < 0.05), UV vs trabecular bone number (Tb.N) (p < 0.05), and UV vs Tb.Th (p < 0.001). In the frontal plane, significant improvement was found in ATT vs structural model index (SMI) (p < 0.01), ATT vs bone volume fraction (BV/TV) (p < 0.01), ATT vs BS/BV (p < 0.001), ATT vs Tb.Th (p < 0.001), ATT vs DA (p < 0.001), and ATT vs modulus (p < 0.001), UV vs SMI (p < 0.01), UV vs BV/TV (p < 0.05), UV vs BS/BV (p < 0.05), UV vs Tb.Th (p < 0.01), UV vs trabecular spacing (p < 0.05), and UV vs modulus (p < 0.01). These data suggested that the combined transmission-reflection QUS method is capable of providing information more relevant to the structural and mechanical properties of trabecular bone.

Changes of elastic constants and anisotropy patterns in trabecular bone during disuse-induced bone loss assessed by poroelastic ultrasound

Currently, the approach most widely used to examine bone loss is the measurement of bone mineral density (BMD) using dual X-ray absorptiometry (DXA). However, bone loss due to immobilization creates changes in bone microarchitecture, which in turn are related to changes in bone mechanical function and competence to resist fracture.Unfortunately, the relationship between microarchitecture and mechanical function within the framework of immobilization and antiresorptive therapy has not being fully investigated. The goal of the present study was to investigate the structure–function relationship in trabecular bone in the real-world situations of a rapidly evolving osteoporosis(disuse), both with and without antiresorptive treatment. We evaluated the structure–function relationship in trabecular bone after bone loss (disuse-induced osteoporosis)and bisphosphonate treatment (antiresorptive therapy using risedronate) in canine trabecular bone using lCT and ultrasound wave propagation. Microstructure values determined from lCT images were used into the anisotropic poroelastic model of wave propagation in order to compute the apparent elastic constants (EC) and elastic anisotropy pattern of bone. Immobilization resulted in a significant reduction in trabecular thickness (Tb.Th) and bone volume fraction (BV/TV), while risedronate treatment combined with immobilization exhibited a lesser reduction in Tb.Th and BV/TV, suggesting that risedronate treatment decelerates bone loss, but it was unable to fully stop it. Risedronate treatment also increased the tissue mineral density (TMD), which when combined with the decrease in Tb.Th and BV/TV may explain the lack of significant differences invBMD in both immobilization and risedronate treated groups. Interestingly, changes inapparent EC were much stronger in the superior–inferior (SI) direction than in the medial–lateral (ML) and anterior–posterior (AP) anatomical directions, producing changes in elastic anisotropy patterns. When data were pooled together, vBMD was able to explain 58% of ultrasound measurements variability, a poroelastic wave propagation analytical model (i.e., BMD modulated by fabric directionality) was able to predict 81%of experimental wave velocity variability, and also explained 91% of apparent EC and changes in elastic anisotropy patterns. Overall, measurements of vBMD were unable to distinguish changes in apparent EC due to immobilization or risedronate treatment.However, anisotropic poroelastic ultrasound (PEUS) wave propagation was able to distinguish functional changes in apparent EC and elastic anisotropy patterns due to immobilization and antiresorptive therapy, providing an enhanced discrimination of anisotropic bone loss and the structure–function relationship in immobilized and risedronate-treated bone, beyond vBMD.

Principal trabecular structural orientation predicted by quantitative ultrasound is strongly correlated with μFEA determined anisotropic apparent stiffness

The microarchitecture and alignment of trabecular bone adapts to the particular mechanical milieu applied to it. Due to this anisotropic mechanical property, measurement orientation has to be taken into consideration when assessing trabecular bone quality and fracture risk prediction. Quantitative ultrasound (QUS) has demonstrated the ability in predicting the principal structural orientation (PSO) of trabecular bone. Although the QUS prediction for PSO is very close to that of μCT, certain angle differences still exist. It remains unknown whether this angle difference can induce significant differences in mechanical properties or not. The objective of this study was to evaluate the mechanical properties in different PSOs predicted using different methods, QUS and μCT, thus to investigate the ability of QUS as a means to predict the PSO of trabecular bone noninvasively. By validating the ability of QUS to predict the PSO of trabecular bone, it is beneficial for future QUS applications because QUS measurements in the PSO can provide information more correlated with the mechanical properties than with other orientations. In this study, seven trabecular bone balls from distal bovine femurs were used to generate finite element models based on the 3-dimensional μCT images. Uniaxial compressive loading was performed on the bone ball models in the finite element analysis (FEA) in six different orientations (three anatomical orientations, two PSOs predicted by QUS and the longest vector of mean intercept length (MIL) tensor calculated by μCT). The stiffness was calculated based on the reaction force of the bone balls under loading, and the von Mises stress results showed that both the mechanical properties in the PSOs predicted by QUS are significantly higher than the anatomical orientations and comparatively close to the longest vector of MIL tensor. The stiffness in the PSOs predicted by QUS is also highly correlated with the stiffness in the MIL tensor orientation (ATTmax vs. MIL, R(2) = 0.98, p < 001; UVmax vs. MIL, R(2) = 0.92, p < 001). These results were validated by in vitro mechanical testing on the bone ball samples. This study demonstrates that the PSO of trabecular bone predicted by QUS has an equally strong apparent stiffness with the orientation predicted by μCT.

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.

Tortuosity and the Averaging of Microvelocity Fields in Poroelasticity

The relationship between the macro- and microvelocity fields in a poroelastic representative volume element (RVE) has not being fully investigated. This relationship is considered to be a function of the tortuosity: a quantitative measure of the effect of the deviation of the pore fluid streamlines from straight (not tortuous) paths in fluid-saturated porous media. There are different expressions for tortuosity based on the deviation from straight pores, harmonic wave excitation, or from a kinetic energy loss analysis. The objective of the work presented is to determine the best expression for tortuosity of a multiply interconnected open pore architecture in an anisotropic porous media. The procedures for averaging the pore microvelocity over the RVE of poroelastic media by Coussy and by Biot were reviewed as part of this study, and the significant connection between these two procedures was established. Success was achieved in identifying the Coussy kinetic energy loss in the pore fluid approach as the most attractive expression for the tortuosity of porous media based on pore fluid viscosity, porosity, and the pore architecture. The fabric tensor, a 3D measure of the architecture of pore structure, was introduced in the expression of the tortuosity tensor for anisotropic porous media. Practical considerations for the measurement of the key parameters in the models of Coussy and Biot are discussed. In this study, we used cancellous bone as an example of interconnected pores and as a motivator for this study, but the results achieved are much more general and have a far broader application than just to cancellous bone.

Microarchitecture and bone quality in the human calcaneus: local variations of fabric anisotropy

The local variability of microarchitecture of human trabecular calcaneus bone is investigated using high-resolution micro-computed tomography (µCT) scanning. The fabric tensor is employed as the measure of the microarchitecture of the pore structure of a porous medium. It is hypothesized that a fabric tensor-dependent poroelastic ultrasound approach will more effectively predict the data variance than will porosity alone. The specific aims of the present study are as follows: (1) to quantify the morphology and local anisotropy of the calcaneus microarchitecture with respect to anatomical directions; (2) to determine the interdependence, or lack thereof, of microarchitecture parameters, fabric, and volumetric bone mineral density (vBMD); and (3) to determine the relative ability of vBMD and fabric measurements in evaluating the variance in ultrasound wave velocity measurements along orthogonal directions in the human calcaneus. Our results show that the microarchitecture in the analyzed regions of human calcanei is anisotropic, with a preferred alignment along the posterior-anterior direction. Strong correlation was found between most scalar architectural parameters and vBMD. However, no statistical correlation was found between vBMD and the fabric components, the measures of the pore microstructure orientation. Therefore, among the parameters usually considered for cancellous bone (ie, classic histomorphometric parameters such as porosity, trabecular thickness, number and separation), only fabric components explain the data variance that cannot be explained by vBMD, a global mass measurement, which lacks the sensitivity and selectivity to distinguish osteoporotic from healthy subjects because it is insensitive to directional changes in bone architecture. This study demonstrates that a multidirectional, fabric-dependent poroelastic ultrasound approach has the capability of characterizing anisotropic bone properties (bone quality) beyond bone mass, and could help to better understand anisotropic changes in bone architecture using ultrasound.

Prediction of trabecular bone principal structural orientation using quantitative ultrasound scanning

Bone has the ability to adapt its structure in response to the mechanical environment as defined as Wolff's Law. The alignment of trabecular structure is intended to adapt to the particular mechanical milieu applied to it. Due to the absence of normal mechanical loading, it will be extremely important to assess the anisotropic deterioration of bone during the extreme conditions, i.e., long term space mission and disease orientated disuse, to predict risk of fractures. The propagation of ultrasound wave in trabecular bone is substantially influenced by the anisotropy of the trabecular structure. Previous studies have shown that both ultrasound velocity and amplitude is dependent on the incident angle of the ultrasound signal into the bone sample. In this work, seven bovine trabecular bone balls were used for rotational ultrasound measurement around three anatomical axes to elucidate the ability of ultrasound to identify trabecular orientation. Both ultrasound attenuation (ATT) and fast wave velocity (UV) were used to calculate the principal orientation of the trabecular bone. By comparing to the mean intercept length (MIL) tensor obtained from μCT, the angle difference of the prediction by UV was 4.45°, while it resulted in 11.67° angle difference between direction predicted by μCT and the prediction by ATT. This result demonstrates the ability of ultrasound as a non-invasive measurement tool for the principal structural orientation of the trabecular bone.

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.

Mixture theory-based poroelasticity as a model of interstitial tissue growth

This contribution presents an alternative approach to mixture theory-based poroelasticity by transferring some poroelastic concepts developed by Maurice Biot to mixture theory. These concepts are a larger RVE and the subRVE-RVE velocity average tensor, which Biot called the micro-macro velocity average tensor. This velocity average tensor is assumed here to depend upon the pore structure fabric. The formulation of mixture theory presented is directed toward the modeling of interstitial growth, that is to say changing mass and changing density of an organism. Traditional mixture theory considers constituents to be open systems, but the entire mixture is a closed system. In this development the mixture is also considered to be an open system as an alternative method of modeling growth. Growth is slow and accelerations are neglected in the applications. The velocity of a solid constituent is employed as the main reference velocity in preference to the mean velocity concept from the original formulation of mixture theory. The standard development of statements of the conservation principles and entropy inequality employed in mixture theory are modified to account for these kinematic changes and to allow for supplies of mass, momentum and energy to each constituent and to the mixture as a whole. The objective is to establish a basis for the development of constitutive equations for growth of tissues.

Role of structural anisotropy of biological tissues in poroelastic wave propagation

Ultrasound waves have a broad range of clinical applications as a non-destructive testing approach in imaging and in the diagnoses of medical conditions. Generally, biological tissues are modeled as an homogenized equivalent medium with an apparent density through which a single wave propagates. Only the first wave arriving at the ultrasound probe is used for the measurement of the speed of sound. However, the existence of a second wave in tissues such as cancellous bone has been reported and its existence is an unequivocal signature of Biot type poroelastic media. To account for the fact that ultrasound is sensitive to microarchitecture as well as density, a fabric-dependent anisotropic poroelastic ultrasound (PEU) propagation theory was recently developed. Key to this development was the inclusion of the fabric tensor - a quantitative stereological measure of the degree of structural anisotropy of bone - into the linear poroelasticity theory. In the present study, this framework is extended to the propagation of waves in several soft and hard tissues. It was found that collagen fibers in soft tissues and the mineralized matrix in hard tissues are responsible for the anisotropy of the solid tissue constituent through the fabric tensor in the model.