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.
Integrated program management process provides management visual assistance through three interrelated charts - system model that identifies each function to be performed, matrix that identifies personnel responsibilities for these functions, process chart that breaks down the functions into discrete tasks.
The importance of the rigorous treatment of photon angular momentum in molecular-collision processes occurring in the presence of intense radiation is investigated. An alternate approximate treatment, which essentially neglects the angular momentum coupling between the photon and the molecular degrees of freedom by averaging over the angular dependence of the interaction matrix elements, is presented and applied to a model calculation. The degeneracy-averaged results of this calculation compare remarkably well with the results of a rigorous calculation, from which we conclude (with reservation) that the explicit consideration of photoangular momentum coupling in molecular-collision problems is unnecessary.
The research effort on the control of large space structures at the University of Houston has concentrated on the mathematical theory of finite-element models; identification of the mass, damping, and stiffness matrix; assignment of damping to structures; and decoupling of structure dynamics. The objective of the work has been and will continue to be the development of efficient numerical algorithms for analysis, control, and identification of large space structures. The major consideration in the development of the algorithms has been the large number of equations that must be handled by the algorithm as well as sensitivity of the algorithms to numerical errors.
The mathematical background for combining and connecting linear, multi-input, multi-output subsystem models into an overall system model is provided. Several examples of subsystem configurations are examined in detail. A description of a MATRIX (sub x) command file to aid in the process of combining and connecting these subsystem models is contained.
This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.
A systematic approach for the identification of terrain media such as vegetation canopy, forest, and snow-covered fields is developed using the optimum polarimetric classifier. The covariance matrices for various terrain cover are computed from theoretical models of random medium by evaluating the scattering matrix elements. The optimal classification scheme makes use of a quadratic distance measure and is applied to classify a vegetation canopy consisting of both trees and grass. Experimentally measured data are used to validate the classification scheme. Analytical and Monte Carlo simulated classification errors using the fully polarimetric feature vector are compared with classification based on single features which include the phase difference between the VV and HH polarization returns. It is shown that the full polarimetric results are optimal and provide better classification performance than single feature measurements.
In this talk we shall discuss algorithms and CAD tools for the design and analysis of structures for high performance applications using advanced composite materials. An extensive mathematical theory for optimal structural (e.g., shape) design was developed over the past thirty years. Aspects of this theory have been used in the design of components for hypersonic vehicles and thermal diffusion systems based on homogeneous materials. Enhancement of the design methods to include optimization of the microstructure of the component is a significant innovation which can lead to major enhancements in component performance. Our work is focused on the adaptation of existing theories of optimal structural design (e.g., optimal shape design) to treat the design of structures using advanced composite materials (e.g., fiber reinforced, resin matrix materials). In this talk we shall discuss models and algorithms for the design of simple structures from composite materials, focussing on a problem in thermal management. We shall also discuss methods for the integration of active structural controls into the design process.
In this paper we study the robustness with respect to stability of the closed-loop system with collocated rate sensor using LQG (mean square rate) optimized compensators. Our main result is that the transmission zeros of the compensator are precisely the structure modes when the actuator/sensor locations are 'pinned' and/or 'clamped': i.e., motion in the direction sensed is not allowed. We have stability even under parameter mismatch, except in the unlikely situation where such a mode frequency of the assumed system coincides with an undamped mode frequency of the real system and the corresponding mode shape is an eigenvector of the compensator transfer function matrix at that frequency. For a truncated modal model - such as that of the NASA LaRC Phase Zero Evolutionary model - the transmission zeros of the corresponding compensator transfer function can be interpreted as the structure modes when motion in the directions sensed is prohibited.
Atmospheric ionizing radiation (AIR) produces chemically active radicals in biological tissues that alter the cell function or result in cell death. The AIR ER-2 flight measurements will enable scientists to study the radiation risk associated with the high-altitude operation of a commercial supersonic transport. The ER-2 radiation measurement flights will follow predetermined, carefully chosen courses to provide an appropriate database matrix which will enable the evaluation of predictive modeling techniques. Explicit scientific results such as dose rate, dose equivalent rate, magnetic cutoff, neutron flux, and air ionization rate associated with those flights are predicted by using the AIR model. Through these flight experiments, we will further increase our knowledge and understanding of the AIR environment and our ability to assess the risk from the associated hazard.
This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.
The durability of polymer matrix composites exposed to harsh environments is a major concern. Surface degradation and damage are observed in polyimide composites used in air at 125 to 300 C. It is believed that diffusion of oxygen into the material and oxidative chemical reactions in the matrix are responsible. Previous work has characterized and modeled diffusion behavior, and thermogravimetric analyses (TGAs) have been carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. However, the model developed using these data was not able to capture behavior seen in isothermal tests, especially those of long duration. A test program that focuses on lower temperatures and makes use of isothermal tests was undertaken to achieve a better understanding of the degradation reactions under use conditions. A new low-cost technique was developed to collect chemical degradation data for isothermal tests lasting over 200 hr in the temperature range 125 to 300 C. Results indicate complex behavior not captured by the previous TGA tests, including the presence of weight-adding reactions. Weight gain reactions dominated in the 125 to 225 C temperature range, while weight loss reactions dominated beyond 225 C. The data obtained from isothermal tests was used to develop a new model of the material behavior. This model was able to fully capture the behavior seen in the tests up to 275 C. Correlation of the current model with both isothermal data at 300 C and high rate TGA test data is mediocre. At 300 C and above, the reaction mechanisms appear to change. Attempts (which failed) to measure non-oxidative degradation indicate that oxidative reactions dominate the degradation at low temperatures. Based on this work, long term isothermal testing in an oxidative atmosphere is recommended for studying the degradation behavior of this class of materials.
NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP#5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 gigahertz attenuation from the disdrometer and the 20 gigahertz time-series as well as to directly measure the 40 gigahertz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data. In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer-derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.
NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 GHz attenuation from the disdrometer and the 20 GHz timeseries as well as to directly measure the 40 GHz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data.In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.
The primary contribution of this paper is an evaluation of the potential value of a tailored trajectory management (TTM) capability for uncrewed aircraft (UA) operators that is proactive in detecting conflicts and developing trajectory-based solutions for UA prior to air traffic control (ATC) performing conflict resolution. The experiment matrix is composed of one baseline simulation that models current air traffic operations without such a capability and four test simulations with different configurations of such a capability. In each simulation, five UA operations into Fort Worth Alliance airport were modeled in the presence of recorded tracks for about 4700 flights on January 18, 2022. The analysis focused on the extent to which such a capability was able to preclude an event that could spike UA operator workload. More specifically, in this study, the emulated UA operator TTM capability for multi-vehicle regional air cargo operations was able to reduce the number of instances of concurrent UA conflicts in the modeled ATC conflict resolution timeframe of 8 minutes or less from three to as low as one.
Missions to Mars will differ from all previous human spaceflight missions in that the crew of astronauts will be required to operate in an Earth-independent manner due to the long communication delays, making the decision on crew size critical to mission success. To support National Aeronautics and Space Administration (NASA) decision makers, NASA’s Engineering and Safety Center (NESC) developed a quantitative methodology for assessing the number of crew necessary to meet primary mission objectives and respond to unforeseen failures based on modeling of NASA’s crew training qualifications. We present results of a custom-built model of a Mars crew qualification and responsibility matrix (CQRM) and discuss the implications on trade space analysis of Mars mission crew size.
Missions to Mars will differ from all previous human spaceflight missions in that the crew of astronauts will be required to operate in an Earth-independent manner due to the long communication delays, making the decision on crew size critical to mission success. To support National Aeronautics and Space Administration (NASA) decision makers, NASA’s Engineering and Safety Center (NESC) developed a quantitative methodology for assessing the number of crew necessary to meet primary mission objectives and respond to unforeseen failures based on modeling of NASA’s crew training qualifications. We present results of a custom-built model of a Mars crew qualification and responsibility matrix (CQRM) and discuss the implications on trade space analysis of Mars mission crew size.
Extending the work of Goree and Gross (1979), solutions are given for a two-dimensional region of unidirectional fibers embedded in an elastic matrix whose initial flaw may take the form of a transverse notch, a rectangular cutout, or a circular hole. Subsequent flaw-induced damage is generated by remote stresses acting parallel to the fibers. For the case of such ductile matrix composites as boron/aluminum, present results indicate that both longitudinal matrix yielding and transverse notch extension must be included in order for the model to agree with experimental results. Little difference is found for the three types of initial damage considered. In all cases, the presence of additional damage changes the nature of stress distribution through the unbroken fibers.
Pade Approximation with Scaling, (PAWS) and Scaling Taylor Exponential Matrix (STEM) computer programs are software tools for design and validation. Provide flexible, user-friendly, language-based interface for input of Markov mathematical methods describing behaviors of fault-tolerant computer systems. Markov models include both recovery from faults via reconfiguration and behaviors of such systems when faults occur. PAWS and STEM produce exact solutions of probability of system failure and provide conservative estimate of number of significant digits in solution. Written in PASCAL and FORTRAN.