Three-dimensional finite element modeling of guided ultrasound wave propagation in intact and healing long bones

The use of guided waves has recently drawn significant interest in the ultrasonic characterization of bone aiming at supplementing the information provided by traditional velocity measurements. This work presents a three-dimensional finite element study of guided wave propagation in intact and healing bones. A model of the fracture callus was constructed and the healing course was simulated as a three-stage process. The dispersion of guided modes generated by a broadband 1-MHz excitation was represented in the time-frequency domain. Wave propagation in the intact bone model was first investigated and comparisons were then made with a simplified geometry using analytical dispersion curves of the tube modes. Then, the effect of callus consolidation on the propagation characteristics was examined. It was shown that the dispersion of guided waves was significantly influenced by the irregularity and anisotropy of the bone. Also, guided waves were sensitive to material and geometrical changes that take place during healing. Conversely, when the first-arriving signal at the receiver corresponded to a nondispersive lateral wave, its propagation velocity was almost unaffected by the elastic symmetry and geometry of the bone and also could not characterize the callus tissue throughout its thickness. In conclusion, guided waves can enhance the capabilities of ultrasonic evaluation.     

The effect of boundary conditions on guided wave propagation in two-dimensional models of healing bone

Guided wave propagation has recently drawn significant interest in the ultrasonic characterization of bone. In this work, we present a two-dimensional computational study of ultrasound propagation in healing bones aiming at monitoring the fracture healing process. In particular, we address the effect of fluid loading boundary conditions on the characteristics of guided wave propagation, using both time and time-frequency (t-f) signal analysis techniques, for three study cases. In the first case, the bone was assumed immersed in blood which occupied the semi-infinite spaces of the upper and lower surfaces of the plate. In the second case, the bone model was assumed to have the upper surface loaded by a 2mm thick layer of blood and the lower surface loaded by a semi-infinite fluid with properties close to those of bone marrow. The third case, involves a three-layer model in which the upper surface of the plate was again loaded by a layer of blood, whereas the lower surface was loaded by a 2mm layer of a fluid which simulated bone marrow. The callus tissue was modeled as an inhomogeneous material and fracture healing was simulated as a three-stage process. The results clearly indicate that the application of realistic boundary conditions has a significant effect on the dispersion of guided waves when compared to simplified models in which the bone's surfaces are assumed free.     

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.     

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

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

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.     

On ultrasound waves guided by bones with coupled soft tissues: A mechanism study and in vitro calibration

The influence of soft tissues coupled with cortical bones on precision of quantitative ultrasound (QUS) has been an issue in the clinical bone assessment in conjunction with the use of ultrasound. In this study, the effect arising from soft tissues on propagation characteristics of guided ultrasound waves in bones was investigated using tubular Sawbones phantoms covered with a layer of mimicked soft tissue of different thicknesses and elastic moduli, and an in vitro porcine femur in terms of the axial transmission measurement. Results revealed that presence of soft tissues can exert significant influence on the propagation of ultrasound waves in bones, leading to reduced propagation velocities and attenuated wave magnitudes compared with the counterparts in a free bone in the absence of soft tissues. However such an effect is not phenomenally dependent on the variations in thickness and elastic modulus of the coupled soft tissues, making it possible to compensate for the coupling effect regardless of the difference in properties of the soft tissues. Based on an in vitro calibration, this study proposed quantitative compensation for the effect of soft tissues on ultrasound waves in bones, facilitating development of high-precision QUS.     

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.     

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.     

Ultrasonically determined thickness of long cortical bones: Three-dimensional simulations of in vitro experiments

It was reported in a previous study that simulated guided wave axial transmission velocities on two-dimensional (2D) numerically reproduced geometry of long bones predicted moderately real in vitro ultrasound data on the same bone samples. It was also shown that fitting of ultrasound velocity with simple analytical model yielded a precise estimate (UTh) for true cortical bone thickness. This current study expands the 2D bone model into three dimensions (3D). To this end, wave velocities and UTh were determined from experiments and from time-domain finite-difference simulations of wave propagation, both performed on a collection of 10 human radii (29 measurement sites). A 3D numerical bone model was developed with tuneable fixed material properties and individualized geometry based on X-ray computed tomography reconstructions of real bones. Simulated UTh data were in good accordance (root-mean-square error was 0.40 mm; r(2)=0.79, p<0.001) with true cortical thickness, and hence the measured phase velocity can be well estimated by using a simple analytical inversion model also in 3D. Prediction of in vitro data was improved significantly (by 10% units) and the upgraded bone model thus explained most of the variability (up to 95% when sites were carefully matched) observed in in vitro ultrasound data.     

Comment on "Three-dimensional finite element modeling of guided ultrasound wave propagation in intact and healing long bones," [J. Acoust. Soc. Am. 121(6), 3907-3921 (2007)

Protopappas et al. performed finite element (FE) studies on the propagation of guided ultrasound waves in intact and healing long bones, and found that the dispersion of guided modes was significantly influenced by the irregularity and anisotropy of the bone. A time-frequency (t-f) method was applied to the obtained signals and several wave modes were identified. However, this technique was unable to quantify their observations and provide monitoring capabilities. One possible reason of this shortcoming may come from the inherent disadvantage of the t-f method. The objective of this comment is to demonstrate that it is necessary to combine other techniques with FE simulations for the extraction of significant quantitative ultrasonic features. Individual guided modes in an isotropic pipe have been theoretically examined using the normal mode expansion (NME) method, and many modes that are missed by the t-f analysis have been identified. It is concluded that in order to extract quantitative ultrasonic features, FE simulations should be supplemented by other techniques such as the NME.     

Numerical simulation of wave propagation in cancellous bone

Numerical simulation of wave propagation is performed through 31 3D volumes of trabecular bone. These volumes were reconstructed from high synchrotron microtomography experiments and are used as the input geometry in a simulation software developed in our laboratory. The simulation algorithm accounts for propagation into both the saturating fluid and bone but absorption is not taken into account. We show that 3D simulation predicts phenomena observed experimentally in trabecular bones : linear frequency dependence of attenuation, increase of attenuation and speed of sound with the bone volume fraction, negative phase velocity dispersion in most of the specimens, propagation of fast and slow wave depending on the orientation of the trabecular network compared to the direction of propagation of the ultrasound. Moreover, the predicted attenuation is in very close agreement with the experimental one measured on the same specimens. Coupling numerical simulation with real bone architecture therefore provides a powerful tool to investigate the physics of ultrasound propagation in trabecular structures.     

Use of Pulse-Echo Ultrasound Imaging in Transcranial Diagnostics of Brain Structures

The results are presented from computer simulations of acoustic pulse propagation in heterogenous media mimicking the human head in two-dimensional and three-dimensional geometries. In the three-dimensional experiment, the cranial bone is presented as a liquid layer with a speed of sound corresponding to that of longitudinal waves in the bone. In the two-dimensional experiment, both longitudinal and transverse waves are considered. Based on data obtained in the numerical experiments, the possibility of obtaining ultrasound images of point scatterers by compensating for aberrations introduced by cranial bones is studied. It is shown that even a simple time delay correction along straight rays greatly improves the quality of an ultrasound image obtained through a nonuniform-thickness solid layer.     

A numerical study on the propagation of Rayleigh and guided waves in cortical bone according to Mindlin's Form II gradient elastic theory

Cortical bone is a multiscale heterogeneous natural material characterized by microstructural effects. Thus guided waves propagating in cortical bone undergo dispersion due to both material microstructure and bone geometry. However, above 0.8 MHz, ultrasound propagates rather as a dispersive surface Rayleigh wave than a dispersive guided wave because at those frequencies, the corresponding wavelengths are smaller than the thickness of cortical bone. Classical elasticity, although it has been largely used for wave propagation modeling in bones, is not able to support dispersion in bulk and Rayleigh waves. This is possible with the use of Mindlin's Form-II gradient elastic theory, which introduces in its equation of motion intrinsic parameters that correlate microstructure with the macrostructure. In this work, the boundary element method in conjunction with the reassigned smoothed pseudo Wigner-Ville transform are employed for the numerical determination of time-frequency diagrams corresponding to the dispersion curves of Rayleigh and guided waves propagating in a cortical bone. A composite material model for the determination of the internal length scale parameters imposed by Mindlin's elastic theory is exploited. The obtained results demonstrate the dispersive nature of Rayleigh wave propagating along the complex structure of bone as well as how microstructure affects guided waves.     

Velocity field of wave-induced local fluid flow in double-porosity media

Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in the double-porosity medium,two porous skeletons are usually assumed,namely,host and inclusions.Of them,the volume ratio of inclusion skeletons is low.All previous studies have ignored the consideration of local fluid flow velocity field in inclusions,and therefore they can not completely describe the physical process of local flow oscillation and should not be applied to the situation where the fluid kinetic energy in inclusions cannot be neglected.In this paper,we analyze the local fluid flow velocity fields inside and outside the inclusion,rewrite the kinetic energy function and dissipation function based on the double-porosity medium model containing spherical inclusions,and derive the reformulated Biot-Rayleigh(BR)equations of elastic wave propagation based on Hamilton’s principle.We present simulation examples with different rock and fluid types.Comparisons between BR equations and reformulated BR equations show that there are significant differences in wave response characteristics.Finally,we compare the reformulated BR equations with the previous theories and experimental data,and the results show that the theoretical results of this paper are correct and effective.     

Transient acoustic wave propagation in rigid porous media: a time-domain approach

Wave propagation of acoustic waves in porous media is considered. The medium is assumed to have a rigid frame, so that the propagation takes place in the air which fills the material. The Euler equation and the constitutive relation are generalized to take into account the dispersive nature of these media. It is shown that the connection between the fractional calculus and the behavior of materials with memory allows time-domain wave equations, the coefficients of which are no longer frequency dependent, to be worked out. These equations are suited for direct and inverse scattering problems, and lead to the complete determination of the porous medium parameters.     

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.     

Analysis of the axial transmission technique for the assessment of skeletal status

Ultrasonic wave propagation in human cortical bone has been investigated in vitro using the so-called axial transmission technique. This technique, which relies on velocity measurement of the first arriving signal, has been used in earlier investigations to study bone status during fracture healing or osteoporosis. Two quasi-point-source elements, one transmitter and one receiver (central frequency 0.5 MHz), were used to generate a wide ultrasonic beam, part of which strikes the sample surface at the longitudinal critical angle, and to receive the signals reflected from the sample surface. The analysis of the field reflected from a fluid-solid interface for an incident spherical wave predicts the existence of a lateral wave propagating along the sample surface at a velocity close to the longitudinal velocity, in addition to the ordinary reflected wave and vibration modes. The transducer-sample and the transmitter-receiver distances were chosen such that the lateral wave is the first arriving signal. Validation of the measuring technique was performed on test materials and was followed by experiments on human cortical bones. Experimental results (arrival time and velocity) strongly suggest that the first detected signal corresponds to the lateral wave predicted by theory.     

About the transition frequency in Biot's theory

Biot's theory of wave propagation in porous media includes a characteristic frequency which is used to distinguish the low-frequency from the high-frequency range. Its determination is based on an investigation of fluid flow through different pore geometries on a smaller scale and a subsequent upscaling process. This idea is limited due to the assumptions made on the smaller scale. It can be enhanced for a general two-phase system by three properties: Inertia of the solid, elasticity of the solid, and frequency dependent corrections of the momentum exchange. They become important for highly porous media with liquids.     

Scholte波与含泥沙两相流介质属性关系的分析及仿真验证

基于四种超声悬浮液模型Urick, Urick-Ament, HT, Mcclements分析了Scholte波在两相流体与多孔介质固体界面处的传播特性. 结合各模型的复波数表达式建立含泥沙流体-多孔介质固体界面波特征方程, 分析了Scholte波速与两相流体积含量、粒径等介质属性的关系. 通过仿真实验获得界面波信号, 运用时延估计获得Scholte波速与泥沙含量、粒径的关系, 发现所得的波速与Urick-Ament和HT理论有相对好的一致性. 关键词： Scholte波 两相流体 多孔介质 泥沙含量     

Measurements of ultrasound velocity and attenuation in numerical anisotropic porous media compared to Biot’s and multiple scattering models

This article quantitatively investigates ultrasound propagation in numerical anisotropic porous media with finite-difference simulations in 3D. The propagation media consist of clusters of ellipsoidal scatterers randomly distributed in water, mimicking the anisotropic structure of cancellous bone. Velocities and attenuation coefficients of the ensemble-averaged transmitted wave (also known as the coherent wave) are measured in various configurations. As in real cancellous bone, one or two longitudinal modes emerge, depending on the micro-structure. The results are confronted with two standard theoretical approaches: Biot’s theory, usually invoked in porous media, and the Independent Scattering Approximation (ISA), a classical first-order approach of multiple scattering theory. On the one hand, when only one longitudinal wave is observed, it is found that at porosities higher than 90% the ISA successfully predicts the attenuation coefficient (unlike Biot’s theory), as well as the existence of negative dispersion. On the other hand, the ISA is not well suited to study two-wave propagation, unlike Biot’s model, at least as far as wave speeds are concerned. No free fitting parameters were used for the application of Biot’s theory. Finally we investigate the phase-shift between waves in the fluid and the solid structure, and compare them to Biot’s predictions of in-phase and out-of-phase motions.