SEARCH · Engineering Papers

Results for “generalized linear model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Two-Dimensional Array Beam Scanning Via Externally and Mutually Injection Locked Coupled Oscillators

Some years ago, Stephan proposed an approach to one dimensional (linear) phased array beam steering which requires only a single phase shifter. This involves the use of a linear array of voltage-controlled electronic oscillators coupled to nearest neighbors. The oscillators are mutually injection locked by controlling their coupling and tuning appropriately. Stephan's approach consists of deriving two signals from a master oscillator, one signal phase shifted with respect to the other by means of a single phase shifter. These two signals are injected into the end oscillators of the array. The result is a linear phase progression across the oscillator array. Thus, if radiating elements are connected to each oscillator and spaced uniformly along a line, they will radiate a beam at an angle to that line determined by the phase gradient which is, in turn, determined by the phase difference between the injection signals.The beam direction is therefore controlled by adjusting this phase difference. Recently, Pogorzelski and York presented a formulation which facilitates theoretical analysis of the above beam steering technique. This was subsequently applied by Pogorzelski in analysis of two dimensional beam steering using perimeter detuning of a coupled oscillator array. The formulation is based on a continuum model in which the oscillator phases are represented by a continuous function satisfying a partial differential equation of diffusion type. This equation can be solved via the Laplace transform and the resulting solution exhibits the dynamic behavior of the array as the beam is steered. Stephan's beam steering technique can be similarly generalized to two-dimensional arrays in which the beam control signals are applied to the oscillators on the perimeter of the array. In this paper the continuum model for this two-dimensional case is developed and the dynamic solution for the corresponding aperture phase function is obtained. The corresponding behavior of the resulting far-zone radiation pattern is displayed as well.

Pogorzelski, Ronald J.↗

NASA NTRS2010

Multi-Model Assessment of the Factors Driving the Ozone Evolution Over the 21st Century

The evolution of ozone from 1960 to 2100 is examined in simulations from fourteen chemistry-climate models. There is general agreement among the models at the broadest levels, with all showing column ozone decreasing at all latitudes from 1960 to around 2000, then increasing at all latitudes over the first half of the 21 st century (21 C), and latitudinal variations in the rate of increase and date of return to historical values. In the second half of the century, ozone is projected to carry on increasing, level off or even decrease depending on the latitude, resulting in variable dates of return to historical values at latitudes where column ozone has declined below those levels. Separation into partial column above and below 20 hPa reveals that these latitudinal differences are almost completely due to differences in the lower stratosphere. At all latitudes, upper stratospheric ozone increases throughout the 21 C and returns to 1960 levels before the end of the century, although there is a spread among the models in dates that ozone returns to historical values. Using multiple linear regression the upper stratospheric ozone increase comes from almost equal contributions due to decrease in halogens and cooling from increased greenhouse gas concentrations. The evolution of lower stratospheric ozone differs with latitude. In the tropical lower stratosphere an increase in tropical upwelling causes a steady decrease in ozone through the 21C, and total column ozone does not return to 1960 levels in all models. In contrast, lower stratospheric and total column ozone in middle and high latitudes increases during the 21 C and returns to 1960 levels. For all models there is an earlier return for ozone to historical levels in the northern hemisphere. This is thought to be due to interhemispheric differences in transport.

Oman, L.↗

NASA NTRS2021

Variabilities in Pm2.5 and Black Carbon Surface Concentrations Reproduced by Aerosol Optical Properties Estimated by In-Situ Data, Ground Based Remote Sensing and Modeling

Because of the increased temporal and spatial resolutions of the sensors onboard recently launched satellites, satellite-based surface aerosol concentration, which is usually estimated from the aerosol optical depth (AOD), is expected to become a strategic tool for air quality studies in the future. By further exploring the relationships of aerosol concentrations and their optical properties using ground observations, the accuracies of these products can be improved. Here, we analyzed collocated observations of surface mass concentrations of fine particulate matter (PM2.5) and black carbon (BC), as well as columnar aerosol optical properties from a sky radiometer and aerosol extinction profiles obtained by multi-axis differential optical absorption spectroscopy (MAX-DOAS), during the 2019–2020 period. We focused the analyses on a daily scale, emphasizing the role of the ultraviolet (UV) spectral region. Generally, the correlation between the AOD of the fine fraction (i.e., fAOD) and the PM2.5 surface concentration was moderately strong, regardless of considerations of boundary layer humidity and altitude. In contrast, the fAOD of the partial column below 1 km, which was obtained by combining sky radiometer and MAX-DOAS retrievals, better reproduced the variability of the PM2.5 and resulted in a linear relationship. In the same manner, we demonstrated that the absorption AOD of the fine fraction (fAAOD) of the partial column was related to the variability of the BC concentration. Analogous analyses based on aerosol products from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) confirmed these findings and highlighted the importance of the shape of the aerosol profile. Overall, our results indicated a remarkable consistency among the retrieved datasets, and between the datasets and MERRA-2 products. These results confirmed the well-known sensitivity to aerosol absorption in the UV spectral region; they also highlighted the efficacy of combined MAX-DOAS and sky radiometer observations.

Alessandro Damiani↗

NASA NTRS1995

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗

NASA NTRS2007

Zero-Gravity Locomotion Simulators: New Ground-Based Analogs for Microgravity Exercise Simulation

Maintaining health and fitness in crewmembers during space missions is essential for preserving performance for mission-critical tasks. NASA's Exercise Countermeasures Project (ECP) provides space exploration exercise hardware and monitoring requirements that lead to devices that are reliable, meet medical, vehicle, and habitat constraints, and use minimal vehicle and crew resources. ECP will also develop and validate efficient exercise prescriptions that minimize daily time needed for completion of exercise yet maximize performance for mission activities. In meeting these mission goals, NASA Glenn Research Center (Cleveland, OH, USA), in collaboration with the Cleveland Clinic (Cleveland, Ohio, USA), has developed a suite of zero-gravity locomotion simulators and associated technologies to address the need for ground-based test analog capability for simulating in-flight (microgravity) and surface (partial-gravity) exercise to advance the health and safety of astronaut crews and the next generation of space explorers. Various research areas can be explored. These include improving crew comfort during exercise, and understanding joint kinematics and muscle activation pattern differences relative to external loading mechanisms. In addition, exercise protocol and hardware optimization can be investigated, along with characterizing system dynamic response and the physiological demand associated with advanced exercise device concepts and performance of critical mission tasks for Exploration class missions. Three zero-gravity locomotion simulators are currently in use and the research focus for each will be presented. All of the devices are based on a supine subject suspension system, which simulates a reduced gravity environment by completely or partially offloading the weight of the exercising test subject s body. A platform for mounting treadmill is positioned perpendicularly to the test subject. The Cleveland Clinic Zero-g Locomotion Simulator (ZLS) utilizes a pneumatic subject load device to apply a near constant gravity-replacement load to the test subject during exercise, and is currently used in conjunction with the General Clinical Research Center for evaluating exercise protocols using a bedrest analog. The enhanced ZLS (eZLS) at NASA Glenn Research Center features an offloaded treadmill that floats on a thin film of air and interfaces to a force reaction frame via variably-compliant isolators, or vibration isolation system. The isolators can be configured to simulate compliant interfaces to the vehicle, which affects mechanical loading to crewmembers during exercise, and has been used to validate system dynamic models for new countermeasures equipment designs, such as the second International Space Station treadmill slated for use in 2010. In the eZLS, the test subject and exercise device can be pitched at the appropriate angle for partial gravity simulations, such as lunar gravity (1/6th earth gravity). On both the eZLS and the NASA-Johnson Space Center standalone ZLS installed at the University of Texas Medical Branch in Galveston, Texas, USA, the subject's body weight relative to the treadmill is controlled via a linear motor subject load device (LM-SLD). The LM-SLD employs a force-feedback closed-loop control system to provide a relatively constant force to the test subject during locomotion, and is set and verified for subject safety prior to each session. Locomotion data were collected during parabolic flight and on the eZLS. The purpose was to determine the similarities and differences between locomotion in actual and simulated microgravity. Subjects attained greater amounts of hip flexion during walking and running during parabolic flight. During running, subjects had greater hip range of motion. Trunk motion was significantly less on the eZLS than during parabolic flight. Peak impact forces, loading rate, and impulse were greater on the eZLS than during parabolic while walking with a low external load (EL) and rning with a high EL. Activation timing differences existed between locations in all muscles except for the rectus femoris. The tibialis anterior and gluteus maximus were active for longer durations on the eZLS than in parabolic flight during walking. Ground reaction forces were greater with the LM-SLD than with bungees during eZLS locomotion. While the eZLS serves as a ground-based analog, researchers should be aware that subtle, but measurable, differences in kinematics and leg musculature activities exist between the environments. Aside from space applications, zero-gravity locomotion simulators may help medical researchers in the future with development of rehabilitative or therapeutic protocols for injured or ill patients. Zero-gravity locomotion simulators may be used as a ground-based test bed to support future missions for space exploration, and eventually may be used to simulate planetary locomotion in partial gravity environments, including the Moon and Mars. Figure: Zero-gravity Locomotion Simulator at the Cleveland Clinic, Cleveland, Ohio, USA

Perusek, Gail P.↗

NASA NTRS2015

A Self-Stabilizing Hybrid Fault-Tolerant Synchronization Protocol

This paper presents a strategy for solving the Byzantine general problem for self-stabilizing a fully connected network from an arbitrary state and in the presence of any number of faults with various severities including any number of arbitrary (Byzantine) faulty nodes. The strategy consists of two parts: first, converting Byzantine faults into symmetric faults, and second, using a proven symmetric-fault tolerant algorithm to solve the general case of the problem. A protocol (algorithm) is also present that tolerates symmetric faults, provided that there are more good nodes than faulty ones. The solution applies to realizable systems, while allowing for differences in the network elements, provided that the number of arbitrary faults is not more than a third of the network size. The only constraint on the behavior of a node is that the interactions with other nodes are restricted to defined links and interfaces. The solution does not rely on assumptions about the initial state of the system and no central clock nor centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. A mechanical verification of a proposed protocol is also present. A bounded model of the protocol is verified using the Symbolic Model Verifier (SMV). The model checking effort is focused on verifying correctness of the bounded model of the protocol as well as confirming claims of determinism and linear convergence with respect to the self-stabilization period.

Malekpour, Mahyar R.↗

NASA NTRS2016

Dynamic Characterization of Galfenol (Fe81.6Ga18.4)

Galfenol has the potential to transform the smart materials industry by allowing for the development of multifunctional, load-bearing devices. One of the primary technical challenges faced by this development is the very limited experimental data on Galfenol's frequency-dependent response to dynamic stress, which is critically important for the design of such devices. This report details a novel and precise characterization of the constitutive behavior of polycrystalline Galfenol (Fe81.6Ga18.4) under quasi-static (1 Hz) and dynamic (4 to 1000 Hz) stress loadings. Mechanical loads are applied using a high-frequency load frame. Quasi-static minor and major hysteresis loop measurements of magnetic flux density and strain are presented for constant electromagnet currents (0 to 1.1 A) and constant magnetic fields 0 to 14 kA/m (0 to 180 Oe). The dynamic stress amplitude for minor and major loops is 2.88 and 31.4 MPa (418 and 4550 psi), respectively. Quasi-static material properties closely match published values for similar Galfenol materials. Quasi-static actuation responses are also measured and compared to quasi-static sensing responses; the high degree of reversibility seen in the comparison is consistent with published measurements and modeling results. Dynamic major and minor loops are measured for dynamic stresses up to 31 MPa (4496 psi) and 1 kHz, and the bias condition resulting in maximum, quasi-static sensitivity. Eddy current effects are quantified by considering solid and laminated Galfenol rods. Three key sources of error in the dynamic measurements are accounted for: (1) electromagnetic noise in strain signals due to Galfenol's magnetic response, (2) error in load signals due to the inertial force of fixturing, and (3) phase misalignment between signals due to conditioning electronics. For dynamic characterization, strain error is kept below 1.2 percent of full scale by wiring two collocated gauges in series (noise cancellation) and through leadwire weaving. Inertial force error is kept below 0.41 percent by measuring the dynamic force in the specimen using a nearly collocated piezoelectric load washer. The phase response of all conditioning electronics is explicitly measured and corrected for. In general, as frequency is increased, the sensing response becomes more linear because of an increase in eddy currents. As frequency increases above approximately 100 Hz, the elbow in the strain-versus-stress response disappears as the active (soft) regime stiffens toward the passive (hard) regime. Under constant-field conditions, the loss factors of the solid rod peak between 200 and 600 Hz, rather than exhibiting a monotonic increase. Compared to the solid rod, the laminated rod exhibits much slower increases in hysteresis with frequency, and its quasi-static behavior extends to higher frequencies. The elastic modulus of the laminated rod decreases between 100 and 300 Hz; this trend is currently unexplained.

Dynamic sensing↗

NASA NTRS1992

A model for the characterization of the spatial properties in vestibular neurons

Quantitative study of the static and dynamic response properties of some otolith-sensitive neurons has been difficult in the past partly because their responses to different linear acceleration vectors exhibited no "null" plane and a dependence of phase on stimulus orientation. The theoretical formulation of the response ellipse provides a quantitative way to estimate the spatio-temporal properties of such neurons. Its semi-major axis gives the direction of the polarization vector (i.e., direction of maximal sensitivity) and it estimates the neuronal response for stimulation along that direction. In addition, the semi-minor axis of the ellipse provides an estimate of the neuron's maximal sensitivity in the "null" plane. In this paper, extracellular recordings from otolith-sensitive vestibular nuclei neurons in decerebrate rats were used to demonstrate the practical application of the method. The experimentally observed gain and phase dependence on the orientation angle of the acceleration vector in a head-horizontal plane was described and satisfactorily fit by the response ellipse model. In addition, the model satisfactorily fits neuronal responses in three-dimensions and unequivocally demonstrates that the response ellipse formulation is the general approach to describe quantitatively the spatial properties of vestibular neurons.

Non-NASA Center↗

NASA NTRS2014

Identifiability of Additive, Time-Varying Actuator and Sensor Faults by State Augmentation

Recent work has provided a set of necessary and sucient conditions for identifiability of additive step faults (e.g., lock-in-place actuator faults, constant bias in the sensors) using state augmentation. This paper extends these results to an important class of faults which may affect linear, time-invariant systems. In particular, the faults under consideration are those which vary with time and affect the system dynamics additively. Such faults may manifest themselves in aircraft as, for example, control surface oscillations, control surface runaway, and sensor drift. The set of necessary and sucient conditions presented in this paper are general, and apply when a class of time-varying faults affects arbitrary combinations of actuators and sensors. The results in the main theorems are illustrated by two case studies, which provide some insight into how the conditions may be used to check the theoretical identifiability of fault configurations of interest for a given system. It is shown that while state augmentation can be used to identify certain fault configurations, other fault configurations are theoretically impossible to identify using state augmentation, giving practitioners valuable insight into such situations. That is, the limitations of state augmentation for a given system and configuration of faults are made explicit. Another limitation of model-based methods is that there can be large numbers of fault configurations, thus making identification of all possible configurations impractical. However, the theoretical identifiability of known, credible fault configurations can be tested using the theorems presented in this paper, which can then assist the efforts of fault identification practitioners.

Upchurch, Jason M.↗

NASA NTRS1999

Development of a 5-Component Balance for Water Tunnel Applications

The principal objective of this research/development effort was to develop a multi-component strain gage balance to measure both static and dynamic forces and moments on models tested in flow visualization water tunnels. A balance was designed that allows measuring normal and side forces, and pitching, yawing and rolling moments (no axial force). The balance mounts internally in the model and is used in a manner typical of wind tunnel balances. The key differences between a water tunnel balance and a wind tunnel balance are the requirement for very high sensitivity since the loads are very low (typical normal force is 90 grams or 0.2 lbs), the need for water proofing the gage elements, and the small size required to fit into typical water tunnel models. The five-component balance was calibrated and demonstrated linearity in the responses of the primary components to applied loads, very low interactions between the sections and no hysteresis. Static experiments were conducted in the Eidetics water tunnel with delta wings and F/A-18 models. The data were compared to forces and moments from wind tunnel tests of the same or similar configurations. The comparison showed very good agreement, providing confidence that loads can be measured accurately in the water tunnel with a relatively simple multi-component internal balance. The success of the static experiments encouraged the use of the balance for dynamic experiments. Among the advantages of conducting dynamic tests in a water tunnel are less demanding motion and data acquisition rates than in a wind tunnel test (because of the low-speed flow) and the capability of performing flow visualization and force/moment (F/M) measurements simultaneously with relative simplicity. This capability of simultaneous flow visualization and for F/M measurements proved extremely useful to explain the results obtained during these dynamic tests. In general, the development of this balance should encourage the use of water tunnels for a wider range of quantitative and qualitative experiments, especially during the preliminary phase of aircraft design.

Suarez, Carlos J.↗

NASA NTRS2025

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path↗

NASA NTRS1974

Effect of canard location and size on canard-wing interference and aerodynamic center shift related to maneuvering aircraft at transonic speeds

A generalized wind-tunnel model, typical of highly maneuverable aircraft, was tested in the Langley 8-foot transonic pressure tunnel at Mach numbers from 0.70 to 1.20 to determine the effects of canard location and size on canard-wing interference effects and aerodynamic center shift at transonic speeds. The canards had exposed areas of 16.0 and 28.0 percent of the wing reference area and were located in the chord plane of the wing or in a position 18.5 percent of the wing mean geometric chord above or below the wing chord plane. Two different wing planforms were tested, one with leading-edge sweep of 60 deg and the other 44 deg; both wings had the same reference area and span. The results indicated that the largest benefits in lift and drag were obtained with the canard above the wing chord plane for both wings tested. The low canard configuration for the 60 deg swept wing proved to be more stable and produced a more linear pitching-moment curve than the high and coplanar canard configurations for the subsonic test Mach numbers.

Gloss, B. B.↗

NASA NTRS2015

Generalized Fluid System Simulation Program (GFSSP) - Version 6

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

Majumdar, Alok↗

NASA NTRS1979

Effect of grazing flow on the acoustic impedance of Helmholtz resonators consisting of single and clustered orifices

A semiempirical fluid mechanical model is derived for the acoustic behavior of thin-walled single orifice Helmholtz resonators in a grazing flow environment. The incident and cavity sound fields are connected in terms of an orifice discharge coefficient whose values are determined experimentally using the two-microphone method. Measurements show that at high grazing flow speeds, acoustical resistance is almost linearly proportional to the grazing flow speed and almost independent of incident sound pressure. The corresponding values of reactance are much smaller and tend towards zero. For thicker-walled orifice plates, resistance and reactance were observed to be less sensitive to grazing flow as the ratio of plate thickness to orifice diameter increased. Loud tones were observed to radiate from a single orifice Helmholtz resonator due to interaction between the grazing flow shear layer and the resonator cavity. Measurements showed that the tones radiated at a Strouhal number equal to 0.26. The effects of grazing flow on the impedance of Helmholtz resonators consisting of clusters of orifices was also studied. In general, both resistance and reaction were found to be virtually independent of orifice relative spacing and number. These findings are valid with and without grazing flow.

Hersch, A. S.↗

NASA NTRS2001

Finite Element Methods for real-time Haptic Feedback of Soft-Tissue Models in Virtual Reality Simulators

We have applied the linear elastic finite element method to compute haptic force feedback and domain deformations of soft tissue models for use in virtual reality simulators. Our results show that, for virtual object models of high-resolution 3D data (>10,000 nodes), haptic real time computations (>500 Hz) are not currently possible using traditional methods. Current research efforts are focused in the following areas: 1) efficient implementation of fully adaptive multi-resolution methods and 2) multi-resolution methods with specialized basis functions to capture the singularity at the haptic interface (point loading). To achieve real time computations, we propose parallel processing of a Jacobi preconditioned conjugate gradient method applied to a reduced system of equations resulting from surface domain decomposition. This can effectively be achieved using reconfigurable computing systems such as field programmable gate arrays (FPGA), thereby providing a flexible solution that allows for new FPGA implementations as improved algorithms become available. The resulting soft tissue simulation system would meet NASA Virtual Glovebox requirements and, at the same time, provide a generalized simulation engine for any immersive environment application, such as biomedical/surgical procedures or interactive scientific applications.

Frank, Andreas O.↗

NASA NTRS1997

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗

NASA NTRS1989

Contact solution algorithms

Two algorithms for obtaining static contact solutions are described in this presentation. Although they were derived for contact problems involving specific structures (a tire and a solid rubber cylinder), they are sufficiently general to be applied to other shell-of-revolution and solid-body contact problems. The shell-of-revolution contact algorithm is a method of obtaining a point load influence coefficient matrix for the portion of shell surface that is expected to carry a contact load. If the shell is sufficiently linear with respect to contact loading, a single influence coefficient matrix can be used to obtain a good approximation of the contact pressure distribution. Otherwise, the matrix will be updated to reflect nonlinear load-deflection behavior. The solid-body contact algorithm utilizes a Lagrange multiplier to include the contact constraint in a potential energy functional. The solution is found by applying the principle of minimum potential energy. The Lagrange multiplier is identified as the contact load resultant for a specific deflection. At present, only frictionless contact solutions have been obtained with these algorithms. A sliding tread element has been developed to calculate friction shear force in the contact region of the rolling shell-of-revolution tire model.

Tielking, John T.↗

NASA NTRS2005

Chemical colloids versus biological colloids: a comparative study for the elucidation of the mechanism of protein fiber formation

Fiber formation from murine serum amyloid A1 (SAA) was compared to the linear aggregation and fiber formation of colloidal gold particles. Here we report the similarities of these processes. Upon incubation with acetic acid, SAA misfolds and adopts a new conformation, which we termed saa. saa apparently is less soluble than SAA in aqueous solution; it aggregates and forms nucleation units and then fibers. The fibers appear as a string of the nucleation units. Additionally, an external electric field promotes saa fiber formation. These properties of saa are reminiscent of colloidal gold formation from gold ions and one-dimensional aggregation of the gold colloids. Colloidal gold particles were also found to be capable of aggregating one-dimensionally under an electric field or in the presence of polylysine. These gold fibers resembled in structure that of saa fibers. In summary, protein aggregation and formation of fibers appear to follow the generalized principles derived in colloidal science for the aggregation of atoms and molecules, including polymers such as polypeptides. The analysis of colloidal gold formation and of one-dimensional aggregation provides a simple model system for the elucidation of some aspects of protein fiber formation.

Colloids↗