Biomechanical characterization and in vitro mechanical injury of elderly human femoral head cartilage: comparison to adult bovine humeral head cartilage.

OBJECTIVES: In vitro mechanical injury of articular cartilage is useful to identify events associated with development of post-traumatic osteoarthritis (OA). To date, many in vitro injury models have used animal cartilage despite the greater clinical relevance of human cartilage. We aimed to characterize a new in vitro injury model using elderly human femoral head cartilage and compare its behavior to that of an existing model with adult bovine humeral head cartilage. DESIGN: Mechanical properties of human and bovine cartilage disks were characterized by elastic modulus and hydraulic permeability in radially confined axial compression, and by Young's modulus, Poisson's ratio, and direction-dependent radial strain in unconfined compression. Biochemical composition was assessed in terms of tissue water, solid, and glycosaminoglycan (GAG) contents. Responses to mechanical injury were assessed by observation of macroscopic superficial tissue cracks and histological measurements of cell viability following single injurious ramp loads at 7 or 70%/s strain rate to 3 or 14 MPa peak stress. RESULTS: Confined compression moduli and Young's moduli were greater in elderly human femoral cartilage vs adult bovine humeral cartilage whereas hydraulic permeability was less. Radial deformations of axially compressed explant disks were more anisotropic (direction-dependent) for the human cartilage. In both cartilage sources, tissue cracking and associated cell death during injurious loading was common for 14 MPa peak stress at both strain rates. CONCLUSION: Despite differences in mechanical properties, acute damage induced by injurious loading was similar in both elderly human femoral cartilage and adult bovine humeral cartilage, supporting the clinical relevance of animal-based cartilage injury models. However, inherent structural differences such as cell density may influence subsequent cell-mediated responses to injurious loading and affect the development of OA.

Direct three-dimensional myocardial strain tensor quantification and tracking using zHARP.

Images of myocardial strain can be used to diagnose heart disease, plan and monitor treatment, and to learn about cardiac structure and function. Three-dimensional (3D) strain is typically quantified using many magnetic resonance (MR) images obtained in two or three orthogonal planes. Problems with this approach include long scan times, image misregistration, and through-plane motion. This article presents a novel method for calculating cardiac 3D strain using a stack of two or more images acquired in only one orientation. The zHARP pulse sequence encodes in-plane motion using MR tagging and out-of-plane motion using phase encoding, and has been previously shown to be capable of computing 3D displacement within a single image plane. Here, data from two adjacent image planes are combined to yield a 3D strain tensor at each pixel; stacks of zHARP images can be used to derive stacked arrays of 3D strain tensors without imaging multiple orientations and without numerical interpolation. The performance and accuracy of the method is demonstrated in vitro on a phantom and in vivo in four healthy adult human subjects.

Stiffness tomography by atomic force microscopy.

The atomic force microscope is a convenient tool to probe living samples at the nanometric scale. Among its numerous capabilities, the instrument can be operated as a nano-indenter to gather information about the mechanical properties of the sample. In this operating mode, the deformation of the cantilever is displayed as a function of the indentation depth of the tip into the sample. Fitting this curve with different theoretical models permits us to estimate the Young's modulus of the sample at the indentation spot. We describe what to our knowledge is a new technique to process these curves to distinguish structures of different stiffness buried into the bulk of the sample. The working principle of this new imaging technique has been verified by finite element models and successfully applied to living cells.

The clinical use of quantitative ultrasound (QUS) in the detection and management of osteoporosis

For the detection and management of osteoporosis and osteoporosis-related fractures, quantitative ultrasound (QUS) is emerging as a relatively low-cost and readily accessible alternative to dual-energy X-ray absorptiometry (DXA) measurement of bone mineral density (BMD) in certain circumstances. The following is a brief, but thorough review of the existing literature with respect to the use of QUS in 6 settings: 1) assessing fragility fracture risk; 2) diagnosing osteoporosis; 3) initiating osteoporosis treatment; 4) monitoring osteoporosis treatment; 5) osteoporosis case finding; and 6) quality assurance and control. Many QUS devices exist that are quite different with respect to the parameters they measure and the strength of empirical evidence supporting their use. In general, heel QUS appears to be most tested and most effective. Overall, some, but not all, heel QUS devices are effective assessing fracture risk in some, but not all, populations, the evidence being strongest for Caucasian females over 55 years old. Otherwise, the evidence is fair with respect to certain devices allowing for the accurate diagnosis of likelihood of osteoporosis, and generally fair to poor in terms of QUS use when initiating or monitoring osteoporosis treatment. A reasonable protocol is proposed herein for case-finding purposes, which relies on a combined assessment of clinical risk factors (CR.F) and heel QUS. Finally, several recommendations are made for quality assurance and control.

Comparison of aortic elasticity determined by cardiovascular magnetic resonance imaging in obese versus lean adults.

The vascular properties of large vessels in the obese have not been adequately studied. We used cardiovascular magnetic resonance imaging to quantify the cross-sectional area and elastic properties of the ascending thoracic and abdominal aorta in 21 clinically healthy obese young adult men and 25 men who were age-matched lean controls. Obese subjects had greater maximal cross-sectional area of the ascending thoracic aorta (984 +/- 252 vs 786 +/- 109 mm(2), p <0.01) and of the abdominal aorta (415 +/- 71 vs 374 +/- 51 mm(2), p <0.05). When indexed for height the differences persisted, but when indexed for body surface area, a significant difference between groups was found only for the maximal abdominal aortic cross-sectional area. The obese subjects also had decreased abdominal aortic elasticity, characterized by 24% lower compliance (0.0017 +/- 0.0004 vs 0.0021 +/- 0.0005 mm(2)/kPa/mm, p <0.01), 22% higher stiffness index beta (6.0 +/- 1.5 vs 4.9 +/- 0.7, p <0.005), and 41% greater pressure-strain elastic modulus (72 +/- 25 vs 51 +/- 9, p <0.005). At the ascending thoracic aorta, only the pressure-strain elastic modulus was different between obese and lean subjects (85 +/- 42 vs 65 +/- 26 kPa, respectively; p <0.05), corresponding to a 31% difference-but arterial compliance and stiffness index were not significantly different between groups. In clinically healthy young adult obese men, obesity is associated with increased cross-sectional aortic area and decreased aortic elasticity.

Assessing dystrophies and other muscle diseases at the nanometer scale by atomic force microscopy.

AIM: Atomic force microscopy nanoindentation of myofibers was used to assess and quantitatively diagnose muscular dystrophies from human patients. MATERIALS & METHODS: Myofibers were probed from fresh or frozen muscle biopsies from human dystrophic patients and healthy volunteers, as well as mice models, and Young's modulus stiffness values were determined. RESULTS: Fibers displaying abnormally low mechanical stability were detected in biopsies from patients affected by 11 distinct muscle diseases, and Young's modulus values were commensurate to the severity of the disease. Abnormal myofiber resistance was also observed from consulting patients whose muscle condition could not be detected or unambiguously diagnosed otherwise. DISCUSSION & CONCLUSION: This study provides a proof-of-concept that atomic force microscopy yields a quantitative read-out of human muscle function from clinical biopsies, and that it may thereby complement current muscular dystrophy diagnosis.

An objective anisotropic elastic plastic model and Algorithm applicable to bone mechanics

The implicit projection algorithm of isotropic plasticity is extended to an objective anisotropic elastic perfectly plastic model. The recursion formula developed to project the trial stress on the yield surface, is applicable to any non linear elastic law and any plastic yield function.A curvilinear transverse isotropic model based on a quadratic elastic potential and on Hill's quadratic yield criterion is then developed and implemented in a computer program for bone mechanics perspectives. The paper concludes with a numerical study of a schematic bone-prosthesis system to illustrate the potential of the model.

Permeability estimates from poro-elastic dispersion analyses of sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Traitement des fractures de l'enfant par enclouage centro-médullaire élastique stable (ECMES) [Treatment of pediatric fractures by intramedullary stable elastic pinning]

During the past decade several new techniques for the treatment of children's fractures respecting the specificity of the growing bone have been described. The goal of all these techniques was to mechanically stabilise the fracture however to preserve a certain instability of the fracture gap itself inducing early callus formation and subsequent consolidation. The dynamic external fixation as well as the elastic stable intramedullary pinning have become accepted means in the treatment of long bone fractures in the paediatric age group. We report our experience of the last seven years with the intramedullary pinning of 105 fractures. Eighty-four were fractures of the femur, 9 of the humerus, 8 of the forearm, and a further 4 of the tibial shaft. The intramedullary elastic pinning represents a simple technique which supports or even enhances the natural process of fracture healing of the growing bone. The method is not very invasive, is cost effective, and allows short hospitalisation. Early physical activity is guaranteed due to early consolidation of the fracture. Complications are rare and the final orthopedic and cosmetic outcome is excellent.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Constraints on the permeability structure of alluvial aquifers from the poro-elastic inversion of multi-frequency P-wave velocity logs

We have explored the possibility of obtaining first-order permeability estimates for saturated alluvial sediments based on the poro-elastic interpretation of the P-wave velocity dispersion inferred from sonic logs. Modern sonic logging tools designed for environmental and engineering applications allow one for P-wave velocity measurements at multiple emitter frequencies over a bandwidth covering 5 to 10 octaves. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and typical emitter frequencies ranging from approximately 1 to 30 kHz, the observable velocity dispersion should be sufficiently pronounced to allow one for reliable first-order estimations of the permeability structure. The corresponding predictions have been tested on and verified for a borehole penetrating a typical surficial alluvial aquifer. In addition to multifrequency sonic logs, a comprehensive suite of nuclear and electrical logs, an S-wave log, a litholog, and a limited number laboratory measurements of the permeability from retrieved core material were also available. This complementary information was found to be essential for parameterizing the poro-elastic inversion procedure and for assessing the uncertainty and internal consistency of corresponding permeability estimates. Our results indicate that the thus obtained permeability estimates are largely consistent with those expected based on the corresponding granulometric characteristics, as well as with the available evidence form laboratory measurements. These findings are also consistent with evidence from ocean acoustics, which indicate that, over a frequency range of several orders-of-magnitude, the classical theory of poro-elasticity is generally capable of explaining the observed P-wave velocity dispersion in medium- to fine-grained seabed sediments