Search NASASearch

Engineering topics

Pennline, James A.

Publications and source records attributed to Pennline, James A..

Overview and Evaluation of a Computational Bone Physiology Modeling Toolchain and Its Application to Testing of Exercise Countermeasures

Prolonged microgravity exposure disrupts natural bone remodeling processes and can lead to a significant loss of bone strength, increasing injury risk during missions and placing astronauts at a greater risk of bone fracture later in life. Resistance-based exercise during missions is used to combat bone loss, but current exercise countermeasures do not completely mitigate the effects of microgravity. To address this concern, we present work to develop a personalizable, site-specific computational modeling toolchain of bone remodeling dynamics to understand and estimate changes in volumetric bone mineral density (BMD) in response to microgravity-induced bone unloading and in-flight exercise. The toolchain is evaluated against data collected from subjects in a 70-day bedrest study and is found to provide insight into the amount of exercise stimulus needed to minimize bone loss, quantitatively predicting post-study volumetric BMD of control subjects who did not perform exercise, and qualitatively predicting the effects of exercise. Results suggest that, with additional data, the toolchain could be improved to aid in developing customized in-flight exercise regimens and predict exercise effectiveness.

Schepelmann, Alexander

Overview and Evaluation of a Computational Bone Physiology Modeling Toolchain and Its Application to Testing of Exercise Countermeasures

Prolonged microgravity exposure disrupts natural bone remodeling processes and can lead to a significant loss of bone strength, increasing injury risk during missions and placing astronauts at a greater risk of bone fracture later in life. Resistance-based exercise during missions is used to combat bone loss, but current exercise countermeasures do not completely mitigate the effects of microgravity. To address this concern, we present work to develop a personalizable, site-specific computational modeling toolchain of bone remodeling dynamics to understand and estimate changes in volumetric bone mineral density (BMD) in response to microgravity-induced bone unloading and in-flight exercise. The toolchain is evaluated against data collected from subjects in a 70-day bed rest study and is found to provide insight into the amount of exercise stimulus needed to minimize bone loss, quantitatively predicting post-study volumetric BMD of control subjects who did not perform exercise, and qualitatively predicting the effects of exercise. Results suggest that, with additional data, the toolchain could be improved to aid in developing customized in-flight exercise regimens and predict exercise effectiveness.

Schepelmann, Alexander

Mapping Bone Mineral Density Obtained by Quantitative Computed Tomography to Bone Volume Fraction

Methods for relating or mapping estimates of volumetric Bone Mineral Density (vBMD) obtained by Quantitative Computed Tomography to Bone Volume Fraction (BVF) are outlined mathematically. The methods are based on definitions of bone properties, cited experimental studies and regression relations derived from them for trabecular bone in the proximal femur. Using an experimental range of values in the intertrochanteric region obtained from male and female human subjects, age 18 to 49, the BVF values calculated from four different methods were compared to the experimental average and numerical range. The BVF values computed from the conversion method used data from two sources. One source provided pre bed rest vBMD values in the intertrochanteric region from 24 bed rest subject who participated in a 70 day study. Another source contained preflight vBMD values from 18 astronauts who spent 4 to 6 months on the ISS. To aid the use of a mapping from BMD to BVF, the discussion includes how to formulate them for purpose of computational modeling. An application of the conversions would be used to aid in modeling of time varying changes in vBMD as it relates to changes in BVF via bone remodeling and/or modeling.

Pennline, James A.

Development of Bone Remodeling Model for Spaceflight Bone Physiology Analysis

Current spaceflight exercise countermeasures do not eliminate bone loss. Astronauts lose bone mass at a rate of 1-2% a month (Lang et al. 2004, Buckey 2006, LeBlanc et al. 2007). This may lead to early onset osteoporosis and place the astronauts at greater risk of fracture later in their lives. NASA seeks to improve understanding of the mechanisms of bone remodeling and demineralization in 1g in order to appropriately quantify long term risks to astronauts and improve countermeasures. NASA's Digital Astronaut Project (DAP) is working with NASA's bone discipline to develop a validated computational model to augment research efforts aimed at achieving this goal.

Microgravity

Development of the NASA Digital Astronaut Project Muscle Model

This abstract describes development work performed on the NASA Digital Astronaut Project Muscle Model. Muscle atrophy is a known physiological response to exposure to a low gravity environment. The DAP muscle model computationally predicts the change in muscle structure and function vs. time in a reduced gravity environment. The spaceflight muscle model can then be used in biomechanical models of exercise countermeasures and spaceflight tasks to: 1) develop site specific bone loading input to the DAP bone adaptation model over the course of a mission; 2) predict astronaut performance of spaceflight tasks; 3) inform effectiveness of new exercise countermeasures concepts.

Mathematical Models

A Computational Model for Simulating Spaceflight Induced Bone Remodeling

An overview of an initial development of a model of bone loss due to skeletal unloading in weight bearing sites is presented. The skeletal site chosen for the initial application of the model is the femoral neck region because hip fractures can be debilitating to the overall performance health of astronauts. The paper begins with the motivation for developing such a model of the time course of change in bone in order to understand the mechanism of bone demineralization experienced by astronauts in microgravity, to quantify the health risk, and to establish countermeasures. Following this, a general description of a mathematical formulation of the process of bone remodeling is discussed. Equations governing the rate of change of mineralized bone volume fraction and active osteoclast and osteoblast are illustrated. Some of the physiology of bone remodeling, the theory of how imbalance in remodeling can cause bone loss, and how the model attempts to capture this is discussed. The results of a preliminary validation analysis that was carried out are presented. The analysis compares a set of simulation results against bone loss data from control subjects who participated in two different bed rest studies. Finally, the paper concludes with outlining the current limitations and caveats of the model, and planned future work to enhance the state of the model.

Pennline, James A.

Evaluating Daily Load Stimulus Formulas in Relating Bone Response to Exercise

Six formulas representing what is commonly referred to as "daily load stimulus" are identified, compared and tested in their ability to relate skeletal mechanical loading to bone maintenance and osteogenic response. Particular emphasis is placed on exercise- induced skeletal loading and whether or not the formulas can adequately capture the known experimental observations of saturation of continuous cyclic loading, rest insertion between repetitions (cycles), recovery of osteogenic potential following saturation, and multiple shorter bouts versus a single long bout of exercise. To evaluate the ability of the formulas to capture these characteristics, a set of exercise scenarios with type of exercise bout, specific duration, number of repetitions, and rest insertion between repetitions is defined. The daily load values obtained from the formulas for the loading conditions of the set of scenarios is illustrated. Not all of the formulas form estimates of daily load in units of stress or in terms of strain at a skeletal site due to the loading force from a specific exercise prescription. The comparative results show that none of the formulas are able to capture all of the experimentally observed characteristics of cyclic loading. However, the enhanced formula presented by Genc et al. does capture several characteristics of cyclic loading that the others do not, namely recovery of osteogenic potential and saturation. This could be a basis for further development of mathematical formulas that more adequately approximates the amount of daily stress at a skeletal site that contributes to bone adaptation.

Pennline, James A.

The Digital Astronaut Project Bone Remodeling Model

Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, particularly in the lower extremities such as the proximal femur: (1) The most commonly used countermeasure against bone loss has been prescribed exercise, (2) However, current exercise countermeasures do not completely eliminate bone loss in long duration, 4 to 6 months, spaceflight, (3,4) leaving the astronaut susceptible to early onset osteoporosis and a greater risk of fracture later in their lives. The introduction of the Advanced Resistive Exercise Device, coupled with improved nutrition, has further minimized the 4 to 6 month bone loss. But further work is needed to implement optimal exercise prescriptions, and (5) In this light, NASA's Digital Astronaut Project (DAP) is working with NASA physiologists to implement well-validated computational models that can help understand the mechanisms of bone demineralization in microgravity, and enhance exercise countermeasure development.

Femoral Neck

Predicting Bone Mechanical State During Recovery After Long-Duration Skeletal Unloading Using QCT and Finite Element Modeling

During long-duration missions at the International Space Station, astronauts experience weightlessness leading to skeletal unloading. Unloading causes a lack of a mechanical stimulus that triggers bone cellular units to remove mass from the skeleton. A mathematical system of the cellular dynamics predicts theoretical changes to volume fractions and ash fraction in response to temporal variations in skeletal loading. No current model uses image technology to gather information about a skeletal site s initial properties to calculate bone remodeling changes and then to compare predicted bone strengths with the initial strength. The goal of this study is to use quantitative computed tomography (QCT) in conjunction with a computational model of the bone remodeling process to establish initial bone properties to predict changes in bone mechanics during bone loss and recovery with finite element (FE) modeling. Input parameters for the remodeling model include bone volume fraction and ash fraction, which are both computed from the QCT images. A non-destructive approach to measure ash fraction is also derived. Voxel-based finite element models (FEM) created from QCTs provide initial evaluation of bone strength. Bone volume fraction and ash fraction outputs from the computational model predict changes to the elastic modulus of bone via a two-parameter equation. The modulus captures the effect of bone remodeling and functions as the key to evaluate of changes in strength. Application of this time-dependent modulus to FEMs and composite beam theory enables an assessment of bone mechanics during recovery. Prediction of bone strength is not only important for astronauts, but is also pertinent to millions of patients with osteoporosis and low bone density.

Chang, Katarina L.

Simulating Bone Loss in Microgravity Using Mathematical Formulations of Bone Remodeling

Most mathematical models of bone remodeling are used to simulate a specific bone disease, by disrupting the steady state or balance in the normal remodeling process, and to simulate a therapeutic strategy. In this work, the ability of a mathematical model of bone remodeling to simulate bone loss as a function of time under the conditions of microgravity is investigated. The model is formed by combining a previously developed set of biochemical, cellular dynamics, and mechanical stimulus equations in the literature with two newly proposed equations; one governing the rate of change of the area of cortical bone tissue in a cross section of a cylindrical section of bone and one governing the rate of change of calcium in the bone fluid. The mechanical stimulus comes from a simple model of stress due to a compressive force on a cylindrical section of bone which can be reduced to zero to mimic the effects of skeletal unloading in microgravity. The complete set of equations formed is a system of first order ordinary differential equations. The results of selected simulations are displayed and discussed. Limitations and deficiencies of the model are also discussed as well as suggestions for further research.

Pennline, James A.

Modeling and Grid Generation of Iced Airfoils

SmaggIce Version 2.0 is a software toolkit for geometric modeling and grid generation for two-dimensional, singleand multi-element, clean and iced airfoils. A previous version of SmaggIce was described in Preparing and Analyzing Iced Airfoils, NASA Tech Briefs, Vol. 28, No. 8 (August 2004), page 32. To recapitulate: Ice shapes make it difficult to generate quality grids around airfoils, yet these grids are essential for predicting ice-induced complex flow. This software efficiently creates high-quality structured grids with tools that are uniquely tailored for various ice shapes. SmaggIce Version 2.0 significantly enhances the previous version primarily by adding the capability to generate grids for multi-element airfoils. This version of the software is an important step in streamlining the aeronautical analysis of ice airfoils using computational fluid dynamics (CFD) tools. The user may prepare the ice shape, define the flow domain, decompose it into blocks, generate grids, modify/divide/merge blocks, and control grid density and smoothness. All these steps may be performed efficiently even for the difficult glaze and rime ice shapes. Providing the means to generate highly controlled grids near rough ice, the software includes the creation of a wrap-around block (called the "viscous sublayer block"), which is a thin, C-type block around the wake line and iced airfoil. For multi-element airfoils, the software makes use of grids that wrap around and fill in the areas between the viscous sub-layer blocks for all elements that make up the airfoil. A scripting feature records the history of interactive steps, which can be edited and replayed later to produce other grids. Using this version of SmaggIce, ice shape handling and grid generation can become a practical engineering process, rather than a laborious research effort.

Vickerman, Mary B.

SmaggIce Version 1.8

SmaggIce version 1.8 is a set of software tools for geometrical modeling of, and generation of grids that conform to, both clean and iced airfoils. Ice shapes, especially those that include rough surfaces, pose difficulty in generating high-quality grids that are essential for predicting airflows by use of computational fluid dynamics. SmaggIce version 1.8 contains software tools needed to overcome this difficulty. For a given airfoil, it allows the user to define the flow domain, decompose the domain into blocks, generate grids, merge gridded blocks, and control the density and smoothness of each grid. Among the unique features of version 1.8 is a thin C-shaped block, called a "viscous sublayer block," which is wrapped around an iced airfoil and its wake line and serves as a means to generate highly controlled grids near the rough ice surface. Users can modify block boundary shapes using control points of non-uniform rational B-spline (NURBS) curves. Concave ice regions can be smoothed during geometrical modeling or creation of the viscous sublayer block.

Vickerman, Mary B.

Preparing and Analyzing Iced Airfoils

SmaggIce version 1.2 is a computer program for preparing and analyzing iced airfoils. It includes interactive tools for (1) measuring ice-shape characteristics, (2) controlled smoothing of ice shapes, (3) curve discretization, (4) generation of artificial ice shapes, and (5) detection and correction of input errors. Measurements of ice shapes are essential for establishing relationships between characteristics of ice and effects of ice on airfoil performance. The shape-smoothing tool helps prepare ice shapes for use with already available grid-generation and computational-fluid-dynamics software for studying the aerodynamic effects of smoothed ice on airfoils. The artificial ice-shape generation tool supports parametric studies since ice-shape parameters can easily be controlled with the artificial ice. In such studies, artificial shapes generated by this program can supplement simulated ice obtained from icing research tunnels and real ice obtained from flight test under icing weather condition. SmaggIce also automatically detects geometry errors such as tangles or duplicate points in the boundary which may be introduced by digitization and provides tools to correct these. By use of interactive tools included in SmaggIce version 1.2, one can easily characterize ice shapes and prepare iced airfoils for grid generation and flow simulations.

Vickerman, Mary B.

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations That Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + definite integral of g(x, t)F(t,y(t))dt with limits between 0 and 1,0 less than or equal to x les than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integral equations arise, e.g., when one applied Green's function techniques to nonlinear two-point boundary value problems of the form y "(x) =f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and y(l) = y(sub l), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trepezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal rule, thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations that Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram

Cray performance data from five benchmarks

The five benchmark programs discussed in TM-88956, February 1987, were run on the CRAY X-MP/24 under different operating systems and compilers. Performance data is reported for runs under early versions of UNICOS and CFT77. The most recent data includes a system of configuration for a X-MP hardware upgrade. Performance figures for the Y-MP are shown for comparison. Differences in the figures are analyzed and discussed.

Pennline, James A.

A comparison of five benchmarks

Five benchmark programs were obtained and run on the NASA Lewis CRAY X-MP/24. A comparison was made between the programs codes and between the methods for calculating performance figures. Several multitasking jobs were run to gain experience in how parallel performance is measured.

Huss, Janice E.