SEARCH · Search NASA
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.
A computer program, PROSE, that produces nonlinear state equations from a simple topological description of an electrical or mechanical network is described. Unnecessary states are also automatically eliminated, so that a simplified terminal circuit model is obtained. The program also prints out the eigenvalues of a linearized system and the sensitivities of the eigenvalue of largest magnitude.
The results of a software reliability experiment conducted in a controlled laboratory setting are reported. The experiment was undertaken to gather data on software failures and is one in a series of experiments being pursued by the Fault Tolerant Systems Branch of NASA Langley Research Center to find a means of credibly performing reliability evaluations of flight control software. The experiment tests a small sample of implementations of radar tracking software having ultra-reliability requirements and uses n-version programming for error detection, and repetitive run modeling for failure and fault rate estimation. The experiment results agree with those of Nagel and Skrivan in that the program error rates suggest an approximate log-linear pattern and the individual faults occurred with significantly different error rates. Additional analysis of the experimental data raises new questions concerning the phenomenon of interacting faults. This phenomenon may provide one explanation for software reliability decay.
The objective of this program was to optimize and scale up the linear perfluoroalkylether stabilization process and to provide test data regarding the fluids' thermal oxidative stability in the presence of metal alloys. The stabilization of Fomblin Z-25 was scaled up to 300 g of fluid. The modified fluid was stable at 316 C in oxygen in the presence of M-50 alloy for more than 24 hrs but less than 40 hrs; the amount of volatiles produced after 24 hrs was 5.5 mg/g. In the presence of Ti(4Al,4Mn) alloy, under the above conditions, following an exposure of 24 hrs, the amount of volatiles formed was 6.2 mg/g; 56 hrs exposure yielded 13.9 mg/g. The commercial fluid at 288 C (in oxygen) in the presence of M-50 after 15 hrs of exposure decomposed extensively, 342 mg/g; in the presence of Ti(4Al,4Mn) alloy after only 8 hrs at 288 C, the amount of volatiles was 191 mg/g. Formulation of the commercial fluid with C2PN3 additive was not as effective as the stabilization processing. All the perfluoroalkylether fluids studied were stable in nitrogen at 343 C. The thermal oxidative stability in the absence of metal alloys varied, with Aflunox exhibiting the best behavior. All the fluids were degraded in oxygen at 316 C during 24 hrs exposure to Ti(4Al,4Mn) alloy with the exception of a perfluoroalkylether substituted triazine and the modified Z-25.
Minimum weight designs of helicopter rotor blades with constraints on multiple coupled flap-lag natural frequencies are studied. Constraints are imposed on the minimum value of the blade autorotational inertia to ensure sufficient rotary inertia to autorotate in case of engine failure and on stresses to guard against structural failure due to blade centrifugal forces. Design variables include blade taper ratio, dimensions of the box beam located inside the airfoil and magnitudes of nonstructural weights. The program CAMRAD is used for the blade modal analysis; the program CONMIN is used for the optimization. A linear approximation involving Taylor series expansion is used to reduce the analysis effort. The procedure contains a sensitivity analysis consisting of analytical derivatives for objective function and constraints on autorotational inertia and stresses. Central finite difference derivatives are used for frequency constraints. Optimal designs are obtained for both rectangular and tapered blades. Using this method, it is possible to design a rotor blade with reduced weight, when compared to a baseline blade, while satisfying all the imposed design requirements.
A spatial operator algebra for modeling the control and trajectory design of manipulation is discussed, with emphasis on its analytical formulation and implementation in the Ada programming language. The elements of this algebra are linear operators whose domain and range spaces consist of forces, moments, velocities, and accelerations. The effect of these operators is equivalent to a spatial recursion along the span of the manipulator. Inversion is obtained using techniques of recursive filtering and smoothing. The operator alegbra provides a high-level framework for describing the dynamic and kinematic behavior of a manipulator and control and trajectory design algorithms. Implementable recursive algorithms can be immediately derived from the abstract operator expressions by inspection, thus greatly simplifying the transition from an abstract problem formulation and solution to the detailed mechanization of a specific algorithm.
Boundary Force Method for Orthotropic Materials (BOFFO) computer program performs stress analyses of two-dimensional, linearly elastic, orthotropic or composite bodies with or without cracks, by use of boundary-force method. Easy to use because only boundaries of region of interest modeled mathematically by use of built-in mesh generator. Stresses computed at any specified point in body. Solutions involving stress-intensity factors and strain-energy-release rates computed for both mode-l and mixed-mode fracture problems. Written in FORTRAN V.
Using a comprehensive flight test database and a parameter identification software program produced at NASA Ames Research Center, a math model of the longitudinal aerodynamics of the Harrier aircraft was formulated. The identification program employed the equation error method using multiple linear regression to estimate the nonlinear parameters. The formulated math model structure adhered closely to aerodynamic and stability/control theory, particularly with regard to compressibility and dynamic manoeuvring. Validation was accomplished by using a three degree-of-freedom nonlinear flight simulator with pilot inputs from flight test data. The simulation models agreed quite well with the measured states. It is important to note that the flight test data used for the validation of the model was not used in the model identification.
This research program involves the study of single droplets and linear arrays of droplets in weakly-buoyant and non-bouyant environments. The primary purpose of the single droplet work was to (1) provide a data base from which to compare droplet array results and (2) to correlate the effects of buoyancy on flame shape. Traditionally convective effects in droplet combustion are represented in terms of the Reynolds number, Re, for forced convection and the Grashof number, Gr, for natural convection. Typically, corrections to the burning rate constant for convective effects are written in terms of Re or Gr(exp 1). The Stefan velocity is not included in these correlations, even though from purely physical reasons, one would expect it to be important, especially at higher burning rates. The flame distortion due to convective effects is less documented quantitatively. Kumagai and Isoda do predict flame shape in natural and forced convective flow fields. Their focus, however, was to predict the actual flame dimensions. Law and co-workers used reduced pressure, high oxidizer ambients to obtain spherical flames. This implies that buoyant flows were reduced at the low pressures, as indicated by a very small Grashof number. Ross et al, however, using scaling arguments showed that reducing the pressure does not have a large effect on the magnitude of the buoyant velocity. Struk et al showed elongated flame shapes during simulated (porous sphere) droplet combustion. The elongation of the flames was due to residual gravity levels aboard the reduced gravity aircraft on which the experiments were conducted. These flame shapes, as well as some data from the literature were interpreted based on a dimensionless grouping called the sphericity parameter, Sp. Sp is the ratio of a characteristic computed buoyant velocity to the Stefan velocity at the flame front. One purpose of the droplet arrays work is to extend the database and theories that exist for single droplets into the regime where droplet interactions are important. The eventual goal being to use the results of this work as inputs to models on spray combustion where droplets seldom burn individually; instead the combustion history of a droplet is strongly influenced by the presence of the neighboring droplets. Recently, Annamali and Ryan have summarized he current status of droplet array, cloud and spray combustion. A number of simplified theories led numerical studies of droplet vaporization/combustion where multiple droplet effects are present are now available. These theories all neglect the effect of buoyancy. Experimentally, most studies to date suffer the effects of buoyancy. It is the dominant transport mechanism in the problem. Only the works of Law and co-worker and more recently by Mikami et al were performed in an environment where buoyancy effects were small. Law and co-workers were limited to high oxygen index, low pressure ambient environments since there studies were conducted in normal gravity.
In this paper, the problem of random vibration of geometrically nonlinear MDOF structures is considered. The solutions obtained by application of two different versions of a stochastic linearization method are compared with exact (F-P-K) solutions. The formulation of a relatively new version of the stochastic linearization method (energy-based version) is generalized to the MDOF system case. Also, a new method for determination of nonlinear sti ness coefficients for MDOF structures is demonstrated. This method in combination with the equivalent linearization technique is implemented in a new computer program. Results in terms of root-mean-square (RMS) displacements obtained by using the new program and an existing in-house code are compared for two examples of beam-like structures.
The results of a Monte Carlo model for calculating fragment fluences and LET spectra are compared to data taken with 600 MeV/nucleon iron ions incident on an accelerator beamline configured for irradiation of biological samples, with no target and with 2, 5 and 8 cm of polyethylene. The model uses a multi-generation nuclear fragmentation code, coupled with a formulation of ionization energy loss based on the Bethe-Bloch equation. In the region where the data are reliable and the experimental acceptance is well understood, many of the features of the experimental spectra are well replicated by the model. To obtain good agreement with the experimental data, the model must allow for at least two generations of fragment production in the target.
Recent empirical studies have shown that multi-angle spectral data can be useful for predicting canopy height, but the physical reason for this correlation was not understood. We follow the concept of canopy spectral invariants, specifically escape probability, to gain insight into the observed correlation. Airborne Multi-Angle Imaging Spectrometer (AirMISR) and airborne Laser Vegetation Imaging Sensor (LVIS) data acquired during a NASA Terrestrial Ecology Program aircraft campaign underlie our analysis. Two multivariate linear regression models were developed to estimate LVIS height measures from 28 AirMISR multi-angle spectral reflectances and from the spectrally invariant escape probability at 7 AirMISR view angles. Both models achieved nearly the same accuracy, suggesting that canopy spectral invariant theory can explain the observed correlation. We hypothesize that the escape probability is sensitive to the aspect ratio (crown diameter to crown height). The multi-angle spectral data alone therefore may not provide enough information to retrieve canopy height globally
We report on airborne atmospheric pressure measurements using new fiber-based laser technology and the oxygen A-band at 765 nm. Remote measurements of atmospheric temperature and pressure are required for a number of NASA Earth science missions and specifically for the Active Sensing of CO2 Emissions Over Nights, Days, and Seasons (ASCENDS) mission. Accurate measurements of tropospheric CO2 on a global scale are very important in order to better understand its sources and sinks and to improve predictions on any future climate change. The ultimate goal of a CO2 remote sensing mission, such as ASCENDS, is to derive the CO2 concentration in the atmosphere in terms of mole fraction in unit of parts-per-million (ppmv) with regard to dry air. Therefore, both CO2 and the dry air number of molecules in the atmosphere are needed in deriving this quantity. O2 is a stable molecule and uniformly mixed in the atmosphere. Measuring the O2 absorption in the atmosphere can thus be used to infer the dry air number of molecules and then used to calculate CO2 concentration. With the knowledge of atmospheric water vapor, we can then estimate the total surface pressure needed for CO2 retrievals. Our work, funded by the ESTO IIP program, uses fiber optic technology and non-linear optics to generate 765 nm laser radiation coincident with the Oxygen A-band. Our pulsed, time gated technique uses several on- and off-line wavelengths tuned to the O2 absorption line. The choice of wavelengths allows us to measure the pressure by using two adjacent O2 absorptions in the Oxygen A-band. Our retrieval algorithm fits the O2 lineshapes and derives the pressure. Our measurements compare favorably with a local weather monitor mounted outside our laboratory and a local weather station.
A summary of computational and experimental aeroelastic and aeroservoelastic (ASE) results for the Semi-Span Super-Sonic Transport (S4T) wind-tunnel model is presented. A broad range of analyses and multiple ASE wind-tunnel tests of the S4T have been performed in support of the ASE element in the Supersonics Program, part of NASA's Fundamental Aeronautics Program. The computational results to be presented include linear aeroelastic and ASE analyses, nonlinear aeroelastic analyses using an aeroelastic CFD code, and rapid aeroelastic analyses using CFD-based reduced-order models (ROMs). Experimental results from two closed-loop wind-tunnel tests performed at NASA Langley's Transonic Dynamics Tunnel (TDT) will be presented as well.
We report on airborne atmospheric pressure measurements using new fiber-based laser technology and the oxygen A-band at 765 nm. Remote measurements of atmospheric temperature and pressure are required for a number of NASA Earth science missions and specifically for the Active Sensing of CO2 Emissions Over Nights, Days, and Seasons (ASCENDS) mission. Accurate measurements of tropospheric CO2 on a global scale are very important in order to better understand its sources and sinks and to improve predictions on any future climate change. The ultimate goal of a CO2 remote sensing mission, such as ASCENDS, is to derive the CO2 concentration in the atmosphere in terms of mole fraction in unit of parts-per-million (ppmv) with regard to dry air. Therefore, both CO2 and the dry air number of molecules in the atmosphere are needed in deriving this quantity. O2 is a stable molecule and uniformly mixed in the atmosphere. Measuring the O2 absorption in the atmosphere can thus be used to infer the dry air number of molecules and then used to calculate CO2 concentration. With the knowledge of atmospheric water vapor, we can then estimate the total surface pressure needed for CO2 retrievals. Our work, funded by the ESTO IIP program, uses fiber optic technology and non-linear optics to generate 765 nm laser radiation coincident with the Oxygen A-band. Our pulsed, time gated technique uses several on- and off-line wavelengths tuned to the O2 absorption line. The choice of wavelengths allows us to measure the pressure by using two adjacent O2 absorptions in the Oxygen A-band. Our retrieval algorithm fits the O2 lineshapes and derives the pressure. Our measurements compare favorably with a local weather monitor mounted outside our laboratory and a local weather station.
The Droplet Deformation and Breakup Model was used to predict deformation of droplets approaching the leading edge stagnation line of an airfoil. The quasi-steady model was solved for each position along the droplet path. A program was developed to solve the non-linear, second order, ordinary differential equation that governs the model. A fourth order Runge-Kutta method was used to solve the equation. Experimental slip velocities from droplet breakup studies were used as input to the model which required slip velocity along the particle path. The center of mass displacement predictions were compared to the experimental measurements from the droplet breakup studies for droplets with radii in the range of 200 to 700 mm approaching the airfoil at 50 and 90 m/sec. The model predictions were good for the displacement of the center of mass for small and medium sized droplets. For larger droplets the model predictions did not agree with the experimental results.
The Droplet Deformation and Breakup Model was used to predict deformation of droplets approaching the leading edge stagnation line of an airfoil. The quasi-steady model was solved for each position along the droplet path. A program was developed to solve the non-linear, second order, ordinary differential equation that governs the model. A fourth order Runge-Kutta method was used to solve the equation. Experimental slip velocities from droplet breakup studies were used as input to the model which required slip velocity along the particle path. The center of mass displacement predictions were compared to the experimental measurements from the droplet breakup studies for droplets with radii in the range of 200 to 700 mm approaching the airfoil at 50 and 90 m/sec. The model predictions were good for the displacement of the center of mass for small and medium sized droplets. For larger droplets the model predictions did not agree with the experimental results.
This paper describes a general purpose Coronagraph Performance Error Budget (CPEB) tool that we have developed under the NASA Exoplanet Exploration Program. The CPEB automates many of the key steps required to evaluate the scattered starlight contrast in the dark hole of a space-based coronagraph. It operates in 3 steps: first, a CodeV or Zemax prescription is converted into a MACOS optical prescription. Second, a Matlab program calls ray-trace code that generates linear beam-walk and aberration sensitivity matrices for motions of the optical elements and line-of-sight pointing, with and without controlled coarse and fine-steering mirrors. Third, the sensitivity matrices are imported by macros into Excel 2007 where the error budget is created. Once created, the user specifies the quality of each optic from a predefined set of PSDs. The spreadsheet creates a nominal set of thermal and jitter motions and combines them with the sensitivity matrices to generate an error budget for the system. The user can easily modify the motion allocations to perform trade studies.
Flow over idealized nozzle computed. Analysis and computer program calculate flow over idealized mixer nozzle. Nozzle idealized by unwrapping it so planform lies in z=O plane. Linearized compresible flow used to calculate flow and mixer shape given loading on mixer. End goal to achieve maximum amount of mixing downstream of mixer while retaining reasonable shape for mixer. Because analysis assumes linearized flow, calculation of effects of deep lobe penetration cannot be made using this program.