Search NASASearch

Engineering topics

Whalen, R. T.

Publications and source records attributed to Whalen, R. T..

The X-ray attenuation characteristics and density of human calcaneal marrow do not change significantly during adulthood

Changes in the material characteristics of bone marrow with aging can be a significant source of error in measurements of bone density when using X-ray and ultrasound imaging modalities. In the context of computed tomography, dual-energy computed techniques have been used to correct for changes in marrow composition. However, dual-energy quantitative computed tomography (DE-QCT) protocols, while increasing the accuracy of the measurement, reduce the precision and increase the radiation dose to the patient in comparison to single-energy quantitative computed tomography (SE-QCT) protocols. If the attenuation properties of the marrow for a particular bone can be shown to be relatively constant with age, it should be possible to use single-energy techniques without experiencing errors caused by unknown marrow composition. Marrow was extracted by centrifugation from 10 mm thick frontal sections of 34 adult cadaver calcanei (28 males, 6 females, ages 17-65 years). The density and energy-dependent linear X-ray attenuation coefficient of each marrow sample were determined. For purposes of comparing our results, we then computed an effective CT number at two GE CT/i scan voltages (80 and 120 kVp) for each specimen. The coefficients of variation for the density, CT number at 80 kVp and CT number at 120 kVp were each less than 1%, and the parameters did not change significantly with age (p > 0.2, r2 < 0.02, power > 0.8 where the minimum acceptable r2 = 0.216). We could demonstrate no significant gender-associated differences in these relationships. These data suggest that calcaneal bone marrow X-ray attenuation properties and marrow density are essentially constant from the third through sixth decades of life.

NASA Discipline Musculoskeletal

Calcaneal loading during walking and running

PURPOSE: This study of the foot uses experimentally measured kinematic and kinetic data with a numerical model to evaluate in vivo calcaneal stresses during walking and running. METHODS: External ground reaction forces (GRF) and kinematic data were measured during walking and running using cineradiography and force plate measurements. A contact-coupled finite element model of the foot was developed to assess the forces acting on the calcaneus during gait. RESULTS: We found that the calculated force-time profiles of the joint contact, ligament, and Achilles tendon forces varied with the time-history curve of the moment about the ankle joint. The model predicted peak talocalcaneal and calcaneocuboid joint loads of 5.4 and 4.2 body weights (BW) during walking and 11.1 and 7.9 BW during running. The maximum predicted Achilles tendon forces were 3.9 and 7.7 BW for walking and running. CONCLUSIONS: Large magnitude forces and calcaneal stresses are generated late in the stance phase, with maximum loads occurring at approximately 70% of the stance phase during walking and at approximately 60% of the stance phase during running, for the gait velocities analyzed. The trajectories of the principal stresses, during both walking and running, corresponded to each other and qualitatively to the calcaneal trabecular architecture.

NASA Center ARC

Physical Activity and Bone Density in Women

A mathematical model of bone density regulation as a function of the daily tissue "effective" stress has been derived. Using the model, the influence of daily activity in the form of a daily loading history has been related to bone density of the calcaneus. The theory incorporates a stress exponent m to account for differences in the importance of magnitude and number of load cycles experienced during daily activity. We have derived a parameter from the model, the "Bone Density Index" (BDI). We have developed a method of collecting daily habitual loading histories using an insole force sensor interfaced to a portable digital data logger carried in a fanny pack. Our goal for this study was to determine a stress exponent, m, relating GRFz history to Calcaneal Bone Mineral Density (CBMD).

Bowley, Susan M.

Reconstruction algorithm for polychromatic CT imaging: application to beam hardening correction

This paper presents a new reconstruction algorithm for both single- and dual-energy computed tomography (CT) imaging. By incorporating the polychromatic characteristics of the X-ray beam into the reconstruction process, the algorithm is capable of eliminating beam hardening artifacts. The single energy version of the algorithm assumes that each voxel in the scan field can be expressed as a mixture of two known substances, for example, a mixture of trabecular bone and marrow, or a mixture of fat and flesh. These assumptions are easily satisfied in a quantitative computed tomography (QCT) setting. We have compared our algorithm to three commonly used single-energy correction techniques. Experimental results show that our algorithm is much more robust and accurate. We have also shown that QCT measurements obtained using our algorithm are five times more accurate than that from current QCT systems (using calibration). The dual-energy mode does not require any prior knowledge of the object in the scan field, and can be used to estimate the attenuation coefficient function of unknown materials. We have tested the dual-energy setup to obtain an accurate estimate for the attenuation coefficient function of K2 HPO4 solution.

NASA Center ARC

Modeling of polychromatic attenuation using computed tomography reconstructed images

This paper presents a procedure for estimating an accurate model of the CT imaging process including spectral effects. As raw projection data are typically unavailable to the end-user, we adopt a post-processing approach that utilizes the reconstructed images themselves. This approach includes errors from x-ray scatter and the nonidealities of the built-in soft tissue correction into the beam characteristics, which is crucial to beam hardening correction algorithms that are designed to be applied directly to CT reconstructed images. We formulate this approach as a quadratic programming problem and propose two different methods, dimension reduction and regularization, to overcome ill conditioning in the model. For the regularization method we use a statistical procedure, Cross Validation, to select the regularization parameter. We have constructed step-wedge phantoms to estimate the effective beam spectrum of a GE CT-I scanner. Using the derived spectrum, we computed the attenuation ratios for the wedge phantoms and found that the worst case modeling error is less than 3% of the corresponding attenuation ratio. We have also built two test (hybrid) phantoms to evaluate the effective spectrum. Based on these test phantoms, we have shown that the effective beam spectrum provides an accurate model for the CT imaging process. Last, we used a simple beam hardening correction experiment to demonstrate the effectiveness of the estimated beam profile for removing beam hardening artifacts. We hope that this estimation procedure will encourage more independent research on beam hardening corrections and will lead to the development of application-specific beam hardening correction algorithms.

NASA Discipline Musculoskeletal

A new frame-based registration algorithm

This paper presents a new algorithm for frame registration. Our algorithm requires only that the frame be comprised of straight rods, as opposed to the N structures or an accurate frame model required by existing algorithms. The algorithm utilizes the full 3D information in the frame as well as a least squares weighting scheme to achieve highly accurate registration. We use simulated CT data to assess the accuracy of our algorithm. We compare the performance of the proposed algorithm to two commonly used algorithms. Simulation results show that the proposed algorithm is comparable to the best existing techniques with knowledge of the exact mathematical frame model. For CT data corrupted with an unknown in-plane rotation or translation, the proposed technique is also comparable to the best existing techniques. However, in situations where there is a discrepancy of more than 2 mm (0.7% of the frame dimension) between the frame and the mathematical model, the proposed technique is significantly better (p < or = 0.05) than the existing techniques. The proposed algorithm can be applied to any existing frame without modification. It provides better registration accuracy and is robust against model mis-match. It allows greater flexibility on the frame structure. Lastly, it reduces the frame construction cost as adherence to a concise model is not required.

NASA Center ARC

Prediction of human gait parameters from temporal measures of foot-ground contact

Investigation of the influence of human physical activity on bone functional adaptation requires long-term histories of gait-related ground reaction force (GRF). Towards a simpler portable GRF measurement, we hypothesized that: 1) the reciprocal of foot-ground contact time (1/tc); or 2) the reciprocal of stride-period-normalized contact time (T/tc) predict peak vertical and horizontal GRF, loading rates, and horizontal speed during gait. GRF data were collected from 24 subjects while they walked and ran at a variety of speeds. Linear regression and ANCOVA determined the dependence of gait parameters on 1/tc and T/tc, and prediction SE. All parameters were significantly correlated to 1/tc and T/tc. The closest pooled relationship existed between peak running vertical GRF and T/tc (r2 = 0.896; SE = 3.6%) and improved with subject-specific regression (r2 = 0.970; SE = 2.2%). We conclude that temporal measures can predict force parameters of gait and may represent an alternative to direct GRF measurements for determining daily histories of habitual lower limb loading quantities necessary to quantify a bone remodeling stimulus.

NASA Center ARC

Application of high-performance computing to numerical simulation of human movement

We have examined the feasibility of using massively-parallel and vector-processing supercomputers to solve large-scale optimization problems for human movement. Specifically, we compared the computational expense of determining the optimal controls for the single support phase of gait using a conventional serial machine (SGI Iris 4D25), a MIMD parallel machine (Intel iPSC/860), and a parallel-vector-processing machine (Cray Y-MP 8/864). With the human body modeled as a 14 degree-of-freedom linkage actuated by 46 musculotendinous units, computation of the optimal controls for gait could take up to 3 months of CPU time on the Iris. Both the Cray and the Intel are able to reduce this time to practical levels. The optimal solution for gait can be found with about 77 hours of CPU on the Cray and with about 88 hours of CPU on the Intel. Although the overall speeds of the Cray and the Intel were found to be similar, the unique capabilities of each machine are better suited to different portions of the computational algorithm used. The Intel was best suited to computing the derivatives of the performance criterion and the constraints whereas the Cray was best suited to parameter optimization of the controls. These results suggest that the ideal computer architecture for solving very large-scale optimal control problems is a hybrid system in which a vector-processing machine is integrated into the communication network of a MIMD parallel machine.

NASA Discipline Musculoskeletal

Dynamic Leg Exercise Improves Tolerance to Lower Body Negative Pressure

These results clearly demonstrate that dynamic leg exercise against the footward force produced by LBNP substantially improves tolerance to LBNP, and that even cyclic ankle flexion without load bearing also increases tolerance. This exercise-induced increase of tolerance was actually an underestimate, because subjects who completed the tolerance test while exercising could have continued for longer periods. Exercise probably increases LBNP tolerance by multiple mechanisms. Tolerance was increased in part by skeletal muscle pumping venous blood from the legs. Rosenhamer and Linnarsson and Rosenhamer also deduced this for subjects cycling during centrifugation, although no measurements of leg volume were made in those studies: they found that male subjects cycling at 98 W could endure 3 Gz centrifugation longer than when they remained relaxed during centrifugation. Skeletal muscle pumping helps maintain cardiac filling pressure by opposing gravity-, centrifugation-, or LBNP-induced accumulation of blood and extravascular fluid in the legs.

Watenpaugh, D. E.

Simulation of Hypo-and Hyper-Gravity Locomotion

We have proposed a novel method of using air pressure to apply an external axial force to the body, coincident with the body's center of mass, that has the potential of enabling Earth-equivalent musculoskeletal forces in space. By reversing the direction of the pressure differential, walking and running at reduced musculoskeletal levels are possible on Earth. We hypothesize that hypo- and hyper-g walking and running can be adequately simulated by this method of loading since gait is primarily governed by gravity acting at the center of mass and only secondarily by its distributed action on limb segments. We report here results from simulated hyper-g walking and running.

Whalen, R. T.

Quantification of Cyclic Ground Reaction Force Histories During Daily Activity in Humans

Theoretical models and experimental studies of bone remodeling suggest that bone density and structure are influenced by local cyclic skeletal tissue stress and strain histories. Estimation of long-term loading histories in humans is usually achieved by assessment of physical activity level by questionnaires, logbooks, and pedometers, since the majority of lower limb cyclic loading occurs during walking and running. These methods provide some indication of the mechanical loading history, but fail to consider the true magnitude of the lower limb skeletal forces generated by various daily activities. These techniques cannot account for individual gait characteristics, gait speed, and unpredictable high loading events that may influence bone mass significantly. We have developed portable instrumentation to measure and record the vertical component of the ground reaction force (GRFz) during normal daily activity. This equipment allows long-term quantitative monitoring of musculoskeletal loads, which in conjunction with bone mineral density assessments, promises to elucidate the relationship between skeletal stresses and bone remodeling.

Breit, G. A.

Determinants of Long Bone Structural Properties

The objective of our research is to determine whether a non-invasive determination of long bone cross-sectional areal properties using only the mineral component of bone accurately predicts the true structural properties. In this study section properties of a whole long bone were compared using two methods: (1) special analysis of bone densitometry data, and (2) experimental determination of flexural rigidities from bone surface strain measurements during controlled loading.

Cleek, T. M.

Mechanical properties of particulate composites based on a body-centered-cubic Mg-Li alloy containing boron

The effect of substituting the Mg metal in Mg-B composites by a Mg-14 wt pct Li solid solution on the ductility of the resulting composite was investigated using elastic modulus measurements on the P/M composite material prepared with a dispersion of B particles (in a vol pct range of 0-30) in a matrix of Mg-14 wt pct Li-1.5 wt pct Al. It was found that the elastic modulus of the composites increased rapidly with increasing boron, with specific stiffness values reaching about two times that of most structural materials. The values of the compression and tensile strengths increased significantly with boron additions. Good tensile ductility was achieved at the level of 10 vol pct B. However, at 20 vol pct B, the Mg-Li composite exhibited only limited tensile ductility (about 2 percent total elongation).

Whalen, R. T.