SEARCH · Engineering Papers
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This document discloses, among other things, a system and method for implementing an algorithm to determine pose, velocity, acceleration or other navigation information using feature tracking data. The algorithm has computational complexity that is linear with the number of features tracked.
Long-duration space exploration missions and sustained lunar or Martian surface operations present greater demands for multifunctional and reusable materials. By scaling down the amount of material to be launched from Earth, both mission cost and risk can be reduced. In this regard, leveraging in-space manufacturing capabilities with reusable feedstock materials is an attractive option, as it will allow for articles to be generated on demand, utilized, and then recycled for additional use. NASA’s Enabling Sustained Presence Using Recyclables (ESPUR) project aims to develop reusable materials using polymer-coated microparticles that are bonded via reversible Diels-Alder reactions, where only modest heat is needed to trigger the reverse reaction and enable reuse. For proof-of-concept demonstration, research is currently focused on the fabrication of epoxy microparticles that contain a copoly(carbonate urethane) coating with maleimide and furan functionalities. Here, we discuss the integration of computational materials modeling approaches to help navigate the large design space in this development effort. We perform molecular dynamics (MD) simulations with atomistic and coarse-grained models of the copolymer, which allow us to evaluate the effects of design parameters like the molecular weight and composition on the molecular interactions and chain dynamics. We show how the properties change with the reversible bonds. We also leverage discrete element method (DEM) simulations to assess how microparticle design parameters like the size ratio and volume fraction can be tuned to increase the packing density and number of microparticle contacts to improve the mechanical properties. Our results demonstrate how computational tools can be used in close collaboration with experimental efforts to accelerate material design.
One or the key objectives of the Applied Research Branch in the Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Allies Research Center is the accelerated introduction of highly parallel machines into a full operational environment. In this report we discuss the performance results obtained from the implementation of some computational fluid dynamics (CFD) applications on the Connection Machine CM-2 and the Intel iPSC/860. We summarize some of the experiences made so far with the parallel testbed machines at the NAS Applied Research Branch. Then we discuss the long term computational requirements for accomplishing some of the grand challenge problems in computational aerosciences. We argue that only massively parallel machines will be able to meet these grand challenge requirements, and we outline the computer science and algorithm research challenges ahead.
A spiral approach concept is proposed as a standard procedure for independent, commercial VTOL operations conducted in close proximity to CTOL operations. A guidance logic is developed for this VTOL application, although the results may be applicable to curved approaches by STOL or CTOL aircraft as well. The guidance concept attempts to maintain constant airspeed along a fixed radius nominal spiral. The presence of wind requires a continuous variation in bank angle and heading rate to remain on the desired path. Linear perturbation analysis is used to select satisfactory feedback gains for commanded bank angle, longitudinal acceleration and vertical speed. A wind estimator detects differences between the predicted and observed wind and uses the result to modify the nominal control. For 4-D guidance the nominal time must be computed as a function of turn angle, which requires the solution of an elliptic integral. Time-varying longitudinal accelerations are necessary for 4-D guidance when the observed wind differs from the predicted wind, or when wind shear is present. The logic and the linear feedback gains have been tested in a nonlinear simulation. Results have generally verified the performance predicted by linear analysis.
A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.
The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear. time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that nonlinear waves are generated in coronal holes by torsional Alfv\'{e}n waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the, fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature. and magnetic field geometry,) that will become available from the recently launched SOHO spacecraft.
The successful development and deployment of space nuclear power and propulsion technology requires a high-fidelity, efficient means of modeling ex-core radiation fields. Monte Carlo particle transport codes provide enable the highest level of fidelity in radiation field modeling but are inefficient for modeling ex-core radiation without modification. This paper details the development of an approach to ex-core radiation field modeling that maintains the fidelity of a Monte Carlo transport approach without sacrificing computational efficiency. This approach is based on recent improvements to an existing Monte Carlo transport acceleration technique known as surface source banking. As the transport of many particles is required to ensure adequate uncertainty in far field transport, analytical reconstruction techniques are employed to represent surface source banks as sampleable series of distributions in particle phase from which an arbitrarily large number of particles can be generated and simulated. All steps of this approach are incorporated into an AMA model to test both the underlying mathematics of the analytical source reconstruction process, as well as the relevance of the technique as a whole towards SNP applications. These tests demonstrated that this approach is capable of accelerating far-field radiation modeling in SNP-relevant scenarios by factors of at least nine over purely eigenvalue-based scenarios without significantly sacrificing simulation accuracy.
Recent interpretations of laser ranging for the LAGEOS satellite have rather conclusively established that the observed acceleration in the node of its orbit is just that expected to exist as a residual effect of the last deglaciation event which ended about 6000 years ago. The nontidal acceleration of rotation would be rather different than that observed if there were any significant melting of high latitude continental ice masses currently ongoing. The sensitivity of the expected nontidal acceleration to variations of several elements of the radial viscoelastic structure of the planet is explored using a new normal mode method for the computation of viscoelastic relaxation spectra. These calculations establish that the most important sensitivity is to variations in the mantle viscosity profile. Although the predicted nontidal acceleration does depend upon lithospheric thickness and on the elastic component of the radial structure, the dependence on these components of the structure is much weaker than it is upon mantle viscosity. The observed J sub 2 is therefore a particularly useful determinant of radial variations in the latter parameter.
We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.
The deployment of high data rate instruments onboard satellites with direct broadcast capability demands the implementation of powerful ground data processing systems. Conversely, to enable the widespread usage of direct broadcast data and the generation of subsequent data products, low cost data processing systems must be available. Reconfigurable computers based on Field Programmable Gate Arrays (FPGA) enable the implementation of ground data processing systems that lead to higher performance and lower cost while keeping a high degree of programmability and maintainability. This paper presents the conceptualization of the MODIS Adaptive Level One Accelerator (ALOA), a prototype under construction at NASA's Goddard Space Flight Center. The ALOA is a computer workstation augmented by a reconfigurable computer. The system will calibrate and geographically locate data from the EOS-AM1 MODerate resolution Imaging Spectroradiometer (MODIS) instrument at near real-time data rates.
The long-term performance of rechargeable battery cells has traditionally been verified through life-testing, a procedure that generally requires significant commitments of funding and test resources. In the situation of nickel hydrogen battery cells, which have the capability of providing extremely long cycle life, the time and cost required to conduct even accelerated testing has become a serious impediment to transitioning technology improvements into spacecraft applications. The utilization of computer simulations to indicate the changes in performance to be expected in response to design or operating changes in nickel hydrogen cells is therefore a particularly attractive tool in advanced battery development, as well as for verifying performance in different applications. Computer-based simulations of the long-term performance of rechargeable battery cells have typically had very limited success in the past. There are a number of reasons for the lack in progress in this area. First, and probably most important, all battery cells are relatively complex electrochemical systems, in which performance is dictated by a large number of interacting physical and chemical processes. While the complexity alone is a significant part of the problem, in many instances the fundamental chemical and physical processes underlying long-term degradation and its effects on performance have not even been understood. Second, while specific chemical and physical changes within cell components have been associated with degradation, there has been no generalized simulation architecture that enables the chemical and physical structure (and changes therein) to be translated into cell performance. For the nickel hydrogen battery cell, our knowledge of the underlying reactions that control the performance of this cell has progressed to where it clearly is possible to model them. The recent development of a relative generalized cell modelling approach provides the framework for translating the chemical and physical structure of the components inside a cell into its performance characteristics over its entire cycle life. This report describes our approach to this task in terms of defining those processes deemed critical in controlling performance over life, and the model architecture required to translate the fundamental cell processes into performance profiles.
A study was performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal of the study was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semi-analytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. In several cases this fixed mode approach resulted in very poor approximations of the stress sensitivities. Almost all of the original modes were required for an accurate sensitivity and for small numbers of modes, the accuracy was extremely poor. To overcome this poor accuracy, two semi-analytical techniques were developed. The first technique accounts for the change in eigenvectors through approximate eigenvector derivatives. The second technique applies the mode acceleration method of transient analysis to the sensitivity calculations. Both result in accurate values of the stress sensitivities with a small number of modes and much lower computational costs than if the vibration modes were recalculated and then used in an overall finite difference method.
The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear, time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that solitary-like waves are generated in coronal holes nonlinearly by torsional Alfven waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature, and magnetic field geometry) that will become available from the recently launched SOHO spacecraft.
On Shaft Data Acquisition System (OSDAS) is a rugged, compact, multiple-channel data acquisition computer system that is designed to record data from instrumentation while operating under extreme rotational centrifugal or gravitational acceleration forces. This system, which was developed for the Heritage Fuel Air Turbine Test (HFATT) program, addresses the problem of recording multiple channels of high-sample-rate data on most any rotating test article by mounting the entire acquisition computer onboard with the turbine test article. With the limited availability of slip ring wires for power and communication, OSDAS utilizes its own resources to provide independent power and amplification for each instrument. Since OSDAS utilizes standard PC technology as well as shared code interfaces with the next-generation, real-time health monitoring system (SPARTAA Scalable Parallel Architecture for Real Time Analysis and Acquisition), this system could be expanded beyond its current capabilities, such as providing advanced health monitoring capabilities for the test article. High-conductor-count slip rings are expensive to purchase and maintain, yet only provide a limited number of conductors for routing instrumentation off the article and to a stationary data acquisition system. In addition to being limited to a small number of instruments, slip rings are prone to wear quickly, and introduce noise and other undesirable characteristics to the signal data. This led to the development of a system capable of recording high-density instrumentation, at high sample rates, on the test article itself, all while under extreme rotational stress. OSDAS is a fully functional PC-based system with 48 channels of 24-bit, high-sample-rate input channels, phase synchronized, with an onboard storage capacity of over 1/2-terabyte of solid-state storage. This recording system takes a novel approach to the problem of recording multiple channels of instrumentation, integrated with the test article itself, packaged in a compact/rugged form factor, consuming limited power, all while rotating at high turbine speeds.
The complete output of the Synchronous Meteorological Satellite was recorded on one inch magnetic tape. A quality control subsystem tests cloud track vectors against four sets of criteria: (1) rejection if best match occurs on correlation boundary; (2) rejection if major correlation peak is not distinct and significantly greater than secondary peak; (3) rejection if correlation is not persistent; and (4) rejection if acceleration is too great. A cloud height program determines cloud optical thickness from visible data and computer infrared emissivity. From infrared data and temperature profile, cloud height is determined. A functional description and electronic schematics of equipment are given.
One of NASA's key goals is to increase the safety and reduce the cost of space transportation. Thus, a key element of NASA's new Integrated Space Transportation Plan is to develop new propulsion, structures, and operations for future generations of reusable launch vehicles (RLVs). As part of this effort to develop the next RLV, the ClCT Program's Computing, Networking, and Information Systems (CNIS) Project is developing and demonstrating collaborative software technologies that use the collective power of the NASA Grid to accelerate spacecraft design. One of these technologies, called AeroDB, automates the execution and monitoring of computational fluid dynamics (CFD) parameter studies on the NASA Grid. About the NASA Grid The NASA Grid, or Information Power Grid,. is being developed to leverage the distributed resources of NASA's many computers. instruments, simulators, and data storage systems. The goal is to use these combined resources to sdve difficult NASA challenges, such as iimulating the entire flight of a space vehicle from ascent to descent.To realize the vision of the NASA Grid, the CNIS Project is developing the software framework and protocols for building domain-specific environments and interfaces, new Grid services based on emerging industry standards, and advanced networking and computing testbeds to support new Grid-based applications such as AeroDB.
The convergence of inviscid and viscous hypersonic flow calculations using a two-dimensional flux-splitting code is accelerated by applying a Richardson-type overrelaxation method. Successful results are presented for various cases; and a 50 percent savings in computer time is usually achieved. An analytical formula for the overrelaxation factor is derived, and the performance of this scheme is confirmed numerically. Moreover, application of this overrelaxation scheme produces a favorable preconditioning for Wynn's epsilon-algorithm. Both techniques have been extended to viscous three-dimensional flows and applied to accelerate the convergence of the compressible Navier-Stokes code. A savings of 40 percent in computer time is achieved in this case.
Simulations of astronaut motions during extravehicular activity (EVA) tasks were performed using computational multibody dynamics methods. The application of computational dynamic simulation to EVA was prompted by the realization that physical microgravity simulators have inherent limitations: viscosity in neutral buoyancy tanks; friction in air bearing floors; short duration for parabolic aircraft; and inertia and friction in suspension mechanisms. These limitations can mask critical dynamic effects that later cause problems during actual EVA's performed in space. Methods of formulating dynamic equations of motion for multibody systems are discussed with emphasis on Kane's method, which forms the basis of the simulations presented herein. Formulation of the equations of motion for a two degree of freedom arm is presented as an explicit example. The four basic steps in creating the computational simulations were: system description, in which the geometry, mass properties, and interconnection of system bodies are input to the computer; equation formulation based on the system description; inverse kinematics, in which the angles, velocities, and accelerations of joints are calculated for prescribed motion of the endpoint (hand) of the arm; and inverse dynamics, in which joint torques are calculated for a prescribed motion. A graphical animation and data plotting program, EVADS (EVA Dynamics Simulation), was developed and used to analyze the results of the simulations that were performed on a Silicon Graphics Indigo2 computer. EVA tasks involving manipulation of the Spartan 204 free flying astronomy payload, as performed during Space Shuttle mission STS-63 (February 1995), served as the subject for two dynamic simulations. An EVA crewmember was modeled as a seven segment system with an eighth segment representing the massive payload attached to the hand. For both simulations, the initial configuration of the lower body (trunk, upper leg, and lower leg) was a neutral microgravity posture. In the first simulation, the payload was manipulated around a circular trajectory of 0.15 m radius in 10 seconds. It was found that the wrist joint theoretically exceeded its ulnal deviation limit by as much as 49. 8 deg and was required to exert torques as high as 26 N-m to accomplish the task, well in excess of the wrist physiological limit of 12 N-m. The largest torque in the first simulation, 52 N-m, occurred in the ankle joint. To avoid these problems, the second simulation placed the arm in a more comfortable initial position and the radius and speed of the circular trajectory were reduced by half. As a result, the joint angles and torques were reduced to values well within their physiological limits. In particular, the maximum wrist torque for the second simulation was only 3 N-m and the maximum ankle torque was only 6 N-m.