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 4 mm 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.     

Fabric dependence of quasi-waves in anisotropic porous media

Assessment of bone loss and osteoporosis by ultrasound systems is based on the speed of sound and broadband ultrasound attenuation of a single wave. However, the existence of a second wave in cancellous bone has been reported and its existence is an unequivocal signature of poroelastic media. To account for the fact that ultrasound is sensitive to microarchitecture as well as bone mineral density (BMD), a fabric-dependent anisotropic poroelastic wave propagation theory was recently developed for pure wave modes propagating along a plane of symmetry in an anisotropic medium. 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 mixed wave modes along an arbitrary direction in anisotropic porous media called quasi-waves. It was found that differences between phase and group velocities are due to the anisotropy of the bone microarchitecture, and that the experimental wave velocities are more accurately predicted by the poroelastic model when the fabric tensor variable is taken into account. This poroelastic wave propagation theory represents an alternative for bone quality assessment beyond BMD.     

Generalization of Biot’s equations with allowance for shear relaxation of a fluid

On the basis of the generalized variational principle for dissipative continuum mechanics, a system of generalized Biot’s equations is derived to describe the wave propagation in a two-phase porous permeable medium in the presence of shear relaxation in the pore-filling fluid. It was shown that the inclusion of shear viscoelasticity of the fluid leads to the appearance of two transverse modes in addition to two longitudinal modes described by the Biot theory. One of the transverse modes is an acoustic mode, whereas the other is a diffusion mode characterized by the linear frequency dependence of phase velocity and attenuation coefficient in the low-frequency region.     

Theoretical and experimental study of the ultrasonic attenuation in bovine cancellous bone

In this study a theoretical approach for the estimation of ultrasonic attenuation is proposed. The approach combines two models which take into account both absorption and scattering. Attenuation due to absorption is studied by using the Biot’s analytical model whereas that due to scattering is described by means of a generalized weak scattering model which is formulated for binary mixtures. The scattering model takes account of the density fluctuation of the porous medium in addition to the propagation velocity fluctuation. For the calculation of the attenuation coefficient due to absorption, experimental values have been used to link size of pores to porosity. The theoretical results have been compared with experimental data obtained on bovine cancellous bone samples filled with water. Using an immersion acoustic transmission method, the ultrasonic attenuation has been measured at a frequency range between 0.1 and 1.0 MHz for 12 bovine cancellous bone samples with a porosity range between 40% and 70%. The prediction of attenuation with this model appears to correspond more closely to its experimentally observed behavior. This study indicates that scattering is the predominant mechanism which is responsible for attenuation in trabecular bone. Furthermore, it shows that the density fluctuations contribute significantly to the phenomenon of attenuation and cannot thus be neglected.     

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.     

Time domain numerical modeling of wave propagation in 2D heterogeneous porous media

This paper deals with the numerical modeling of wave propagation in porous media described by Biot’s theory. The viscous efforts between the fluid and the elastic skeleton are assumed to be a linear function of the relative velocity, which is valid in the low-frequency range. The coexistence of propagating fast compressional wave and shear wave, and of a diffusive slow compressional wave, makes numerical modeling tricky. To avoid restrictions on the time step, the Biot’s system is splitted into two parts: the propagative part is discretized by a fourth-order ADER scheme, while the diffusive part is solved analytically. Near the material interfaces, a space–time mesh refinement is implemented to capture the small spatial scales related to the slow compressional wave. The jump conditions along the interfaces are discretized by an immersed interface method. Numerical experiments and comparisons with exact solutions confirm the accuracy of the numerical modeling. The efficiency of the approach is illustrated by simulations of multiple scattering.     

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.     

Propagation of two longitudinal waves in a cancellous bone with the closed pore boundary

Ultrasound propagation in cancellous bone (porous media) under the condition of closed pore boundaries was investigated. A cancellous bone and two plate-like cortical bones obtained from a racehorse were prepared. A water-immersion ultrasound technique in the MHz range and a three-dimensional elastic finite-difference time-domain (FDTD) method were used to investigate the waves. The experiments and simulations showed a clear separation of the incident longitudinal wave into fast and slow waves. The findings advance the evaluation of bones based on the two-wave phenomenon for in vivo assessment.     

各向异性渗流条件下弹性波的传播特征

研究了弹性波在非均匀裂纹孔隙介质中的传播特性,建立了各向异性喷射流模型.当弹性波通过裂纹孔隙介质时,由于波的扰动及裂纹和孔隙几何结构的不一致,导致在裂纹内部及裂纹与周边孔隙之间同时存在着流体压力梯度.此时的弹性波波动响应中包含着裂纹内连通性特征和背景孔隙渗透率信息.流体的动态流动过程使得介质的等效弹性参数为复数(非完全弹性),并且具有频率依赖性.当弹性波为低频和高频极限时,介质为完全弹性;当处于中间频段时,波有衰减和频率依赖.裂纹孔隙介质的各向异性连通性(渗透率)对应着各向异性特征频率(当渗流长度等于非均匀尺度时的弹性波频率),波的传播受到裂纹内连通性的影响.在一定频段内,随着裂纹厚度的增加,将出现第二峰值,峰值大小同时受到裂纹厚度和半径的影响.     

Ultrasonic pulse waves in cancellous bone analyzed by finite-difference time-domain methods

The trabecular frame of cancellous bone has a high degree of porosity, anisotropy and inhomogeneity. The propagation of ultrasonic waves in cancellous bone is significantly affected by the trabecular structure. In this paper, two two-dimensional finite-difference time-domain (FDTD) methods, which were the popular viscoelastic FDTD method for a viscoelastic medium and Biot's FDTD method for a fluid-saturated porous medium, have been applied to numerically analyze the ultrasonic pulse waves propagating through bovine cancellous bone in the directions parallel and perpendicular to the trabecular alignment. The Biot's fast and slow longitudinal waves, which were identified in previous experiments for the propagation parallel to the trabecular orientation, could be analyzed using Biot's FDTD method rather than the viscoelastic FDTD method. For the single wave propagation in the perpendicular direction, on the other hand, the viscoelastic FDTD result was found to be in more good agreement with the experimental result.     

Numerical and experimental study on the wave attenuation in bone--FDTD simulation of ultrasound propagation in cancellous bone

In cancellous bone, longitudinal waves often separate into fast and slow waves depending on the alignment of bone trabeculae in the propagation path. This interesting phenomenon becomes an effective tool for the diagnosis of osteoporosis because wave propagation behavior depends on the bone structure. Since the fast wave mainly propagates in trabeculae, this wave is considered to reflect the structure of trabeculae. For a new diagnosis method using the information of this fast wave, therefore, it is necessary to understand the generation mechanism and propagation behavior precisely. In this study, the generation process of fast wave was examined by numerical simulations using elastic finite-difference time-domain (FDTD) method and experimental measurements. As simulation models, three-dimensional X-ray computer tomography (CT) data of actual bone samples were used. Simulation and experimental results showed that the attenuation of fast wave was always higher in the early state of propagation, and they gradually decreased as the wave propagated in bone. This phenomenon is supposed to come from the complicated propagating paths of fast waves in cancellous bone.     

曲线坐标系下孔隙介质圆柱纵向表面波的传播特性

为了研究孔隙介质圆柱纵向表面波的传播规律,分析其频散和衰减特性,在正交曲线坐标系下建立了表面波的频散方程,通过数值计算得到频散曲线,将纵向导波最低模态与表面波进行对比,并分析了曲率半径及孔隙参数对表面波频散和衰减的影响。结果表明,当频率足够大时,导波最低模态的频散曲线与表面波近似;在同一频率下,表面波的相速度随曲率半径的增大而增大,随孔隙度的增大而减小;表面波的衰减随孔隙度的增大而增大。研究结果为开展孔隙介质圆柱结构的超声无损评价提供了一定的理论参考。     

Propagation of longitudinal elastic waves in a fluid-saturated porous medium with spherical inclusions

The Frenkel-Biot theory is used to study the propagation of a longitudinal harmonic wave of the first kind in an isotropic porous matrix with inclusions contrasting in elastic properties and hydrodynamic permeability. The generation of elastic waves of the second kind at the boundaries of inclusions is taken into account. The effective wave number of the longitudinal wave is calculated using the equations of multiple scattering theory. The characteristic size of inhomogeneities is assumed to be much greater than the size of pores. The parameters of the model used for calculations correspond to sandstone with centimeter-scale inhomogeneities. The presence of such inhomogeneities is typical of sedimentary rocks. Calculations show that, in the frequency range of acoustic logging, the effective attenuation factor of the longitudinal wave may noticeably exceed the attenuation factors of longitudinal waves of the first kind in both matrix and inclusions. From the results obtained, it follows that, when studying the propagation of elastic waves in fluid-saturated porous media, it is necessary to take into account the hydrodynamic effects associated with the filtration overflows that arise at the boundaries of inhomogeneities.     

Ultrasonic wave propagation in cancellous and cortical bone: prediction of some experimental results by Biot's theory.

Pulse transmission ultrasound was used to determine the longitudinal wave speed along the direction of trabecular alignment in 32 water-saturated anisotropic tibial bovine cancellous bone samples and in one cortical bone sample also from the bovine tibia. These results are compared to published ultrasound wave-speed data obtained from bovine femoral specimens. Nonlinear regression was used to fit Biot's theory to the data. The correlation coefficient for regression analysis between the experimental ultrasound velocities and the velocities predicted by Biot's theory was r = 0.78.     

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₂, 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₂ = 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.3 m/s. The optimized value of s₂ 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.     

Estimation of critical and viscous frequencies for Biot theory in cancellous bone

The use of Biot theory for modelling ultrasonic wave propagation in porous media involves the definition of a "critical frequency" above which both fast and slow compressional waves will, in principle, propagate. Critical frequencies have been evaluated for healthy and osteoporotic cancellous bone filled with water or marrow, using data from the literature. The range of pore sizes in bone gives rise to a critical frequency band rather than a single critical frequency, the mean of which is lower for osteoporotic bone than normal bone. However, the critical frequency is a theoretical concept and previous researchers considered a more realistic "viscous frequency" above which both fast and slow waves may be experimentally observed. Viscous frequencies in bone are found to be several orders of magnitude greater than calculated critical frequencies. Whereas two waves may well be observed at all ultrasonic frequencies for water-filled cancellous bone at 20 degrees C, it is probable megahertz frequencies would be needed for observation of two waves in vivo.     

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.     

A multiscale poromicromechanical approach to wave propagation and attenuation in bone

Ultrasonics is an important diagnostic tool for bone diseases, as it allows for non-invasive assessment of bone tissue quality through mass density–elasticity relationships. The latter are, however, quite complex for fluid-filled porous media, which motivates us to develop a rigorous multiscale poromicrodynamics approach valid across the great variety of different bone tissues. Multiscale momentum and mass balance, as well as kinematics of a hierarchical double porous medium, together with Darcy’s law for fluid flow and micro–poro-elasticity for the solid phase of bone, give access to the so-called dispersion relation, linking the complex wave numbers to corresponding wave frequencies. Experimentally validated results show that 2.25 MHz acoustical signals transmit healthy cortical bone (exhibiting a low vascular porosity) only in the form of fast waves, agreeing very well with experimental data, while both fast and slow waves transmit highly osteoporotic as well as trabecular bone (exhibiting a large vascular porosity). While velocities and wavelengths of both fast and slow waves, as well as attenuation lengths of slow waves, are always monotonously increasing with the permeability of the bone sample, the attenuation length of fast waves shows a minimum when considered as function of the permeability.     

Theory of sound propagation in superfluid-filled porous media

The theory of sound propagation in macroscopically isotropic and homogeneous porous media saturated with superfluid ⁴He has been developed neglecting all damping processes. The case when the normal fluid component is locked inside a porous medium by viscous forces is investigated in detail. It is shown that in this case one shear wave and two longitudinal, fast and slow, waves exist. Fast wave as well as slow wave is accompanied with temperature oscillations. The velocities of these waves are obtained.