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.
The long-life deep-space missions associated with NASA's X2000 Advanced Flight Systems Program creates many unprecedented challenges for us.
Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit (CPU) on the order of ~30s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.
Advances in flood monitoring/forecasting have been constrained by the difficulty in estimating rainfall continuously over space (catchment-, national-, continental-, or even global-scale areas) and flood-relevant time scale. With the recent availability of satellite rainfall estimates at fine time and space resolution, this paper describes a prototype research framework for global flood monitoring by combining real-time satellite observations with a database of global terrestrial characteristics through a hydrologically relevant modeling scheme. Four major components included in the framework are (1) real-time precipitation input from NASA TRMM-based Multi-satellite Precipitation Analysis (TMPA); (2) a central geospatial database to preprocess the land surface characteristics: water divides, slopes, soils, land use, flow directions, flow accumulation, drainage network etc.; (3) a modified distributed hydrological model to convert rainfall to runoff and route the flow through the stream network in order to predict the timing and severity of the flood wave, and (4) an open-access web interface to quickly disseminate flood alerts for potential decision-making. Retrospective simulations for 1998-2006 demonstrate that the Global Flood Monitor (GFM) system performs consistently at both station and catchment levels. The GFM website (experimental version) has been running at near real-time in an effort to offer a cost-effective solution to the ultimate challenge of building natural disaster early warning systems for the data-sparse regions of the world. The interactive GFM website shows close-up maps of the flood risks overlaid on topography/population or integrated with the Google-Earth visualization tool. One additional capability, which extends forecast lead-time by assimilating QPF into the GFM, also will be implemented in the future.
This presentation presents experimental work conducted in the Icing Research Tunnel at NASA Glenn Research Center to characterize the velocity of large drops in the tunnel test section. Some icing spray clouds with large-sized drops were generated with Mod1 nozzles at low nozzle air pressure of 2 to 4 psig for various tunnel air speeds. Drop diameters and drop velocities were measured via high-resolution imaging with a Particle Imaging Particle Tracking Velocimetry probe developed by Artium Technologies. The probe was mounted at four different locations aligned with the centerline of the test section from near the end of the contraction to the constant height test section part of the tunnel. CFD analyses were performed. It showed that the probe head geometry affects the local air flow between the prongs and in front of the probe. Initial analysis of the air velocity data during the test also indicated that the probe mounting stand has blockage effect on the local tunnel air velocity measurement by a Pitot static probe affixed on the mounting plate. Those findings were later verified in the Icing Research Tunnel using a new Pitot probe design with a linear motion system to measure the local tunnel air velocity with and without the probe. As a result, additional drop trajectory simulations as airflow moving towards the probe head were run for a Langmuir-D 7-bin drop size distribution of a spray cloud with a nominally large value of medium volumetric diameter. The simulation results helped identify a critical drop-size threshold of 300 µm above which the velocities of larger drops are not affected by the adverse pressure gradient generated by the probe head due to their large drop inertia. From the dimensional analysis of the drop velocity measurement data obtained, it showed that at the tunnel test section reference location a generalized empirical correlation was developed for the non-dimensional drop velocity as a function of the non-dimensional corrected drop diameter independent of the actual spraybar pressure settings and the tunnel air speeds. The generalized curve-fit correlation showed that the drop velocity was universally asymptotic to about 86 percent of the corresponding tunnel air speed at the test section reference location for the largest drop diameter captured by the probe. Further evaluation of this correlation is recommended to assess its applicability for Supercooled Large Drop icing scaling applications in the Icing Research Tunnel.
This software provides a programming interface for automating data collection with a PhaseCam interferometer from 4D Technology, and distributing the image-processing algorithm across a cluster of general-purpose computers. Multiple instances of 4Sight (4D Technology s proprietary software) run on a networked cluster of computers. Each connects to a single server (the controller) and waits for instructions. The controller directs the interferometer to several images, then assigns each image to a different computer for processing. When the image processing is finished, the server directs one of the computers to collate and combine the processed images, saving the resulting measurement in a file on a disk. The available software captures approximately 100 images and analyzes them immediately. This software separates the capture and analysis processes, so that analysis can be done at a different time and faster by running the algorithm in parallel across several processors. The PhaseCam family of interferometers can measure an optical system in milliseconds, but it takes many seconds to process the data so that it is usable. In characterizing an adaptive optics system, like the next generation of astronomical observatories, thousands of measurements are required, and the processing time quickly becomes excessive. A programming interface distributes data processing for a PhaseCam interferometer across a Windows computing cluster. A scriptable controller program coordinates data acquisition from the interferometer, storage on networked hard disks, and parallel processing. Idle time of the interferometer is minimized. This architecture is implemented in Python and JavaScript, and may be altered to fit a customer s needs.
This paper reports on an initial assessment of using a Field-Programmable Gate Array (FPGA) computational device as a new tool for solving structural mechanics problems. A FPGA is an assemblage of binary gates arranged in logical blocks that are interconnected via software in a manner dependent on the algorithm being implemented and can be reprogrammed thousands of times per second. In effect, this creates a computer specialized for the problem that automatically exploits all the potential for parallel computing intrinsic in an algorithm. This inherent parallelism is the most important feature of the FPGA computational environment. It is therefore important that if a problem offers a choice of different solution algorithms, an algorithm of a higher degree of inherent parallelism should be selected. It is found that in structural analysis, an 'analog computer' style of programming, which solves problems by direct simulation of the terms in the governing differential equations, yields a more favorable solution algorithm than current solution methods. This style of programming is facilitated by a 'drag-and-drop' graphic programming language that is supplied with the particular type of FPGA computer reported in this paper. Simple examples in structural dynamics and statics illustrate the solution approach used. The FPGA system also allows linear scalability in computing capability. As the problem grows, the number of FPGA chips can be increased with no loss of computing efficiency due to data flow or algorithmic latency that occurs when a single problem is distributed among many conventional processors that operate in parallel. This initial assessment finds the FPGA hardware and software to be in their infancy in regard to the user conveniences; however, they have enormous potential for shrinking the elapsed time of structural analysis solutions if programmed with algorithms that exhibit inherent parallelism and linear scalability. This potential warrants further development of FPGA-tailored algorithms for structural analysis.
Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit(CPU)on the order of ~30 s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.
The National Aeronautics and Space Administration (NASA) intends to develop and maintain a human-lunar presence, requiring infrastructure on the Lunar South Pole. To keep pace with the continued growth in cislunar activity, there is national interest in developing a global lunar infrastructure network for power and Communication, Position, Navigation, and Timing (CPNT). Demonstrating autonomous construction capabilities is an infrastructure-enabling objective. This study explores several architecture trade studies to assess the feasibility of alternative lunar infrastructure concepts. 50 m towers, for instance, enable solar power generation up to 99% of a lunar year in addition to acting as wide-range communication and navigation surface relays. The infrastructure network considered in this study consisted of tower platforms, utilities hosted on the towers, and both crewed and robotic end users. High-level trades were analyzed, such as tower height, location, and distributed vs centralized power grids. Results indicated that even with minimal permanent surface assets (darkness survival power loads on the order of 5 kW or greater), at any of the five locations considered, taller towers resulted in lower infrastructure mass than shorter towers; up to 24% system mass savings were realized by reducing battery mass. Simultaneously, taller towers increased the range of communication coverage; a 50 m tower on the Connecting Ridge had 5 kg more structure and 76% more line-of-sight coverage to a two-meter-tall end user than a 10 m tower with the same location and user height. Leveraging a grid of tall towers for nominal power and communication would not only save up to 2,800 kg in battery mass but would be more maintainable and allow users to take advantage of generational upgrades.
The author has identified the following significant results. The use of clustering methods permits the development of relatively fast classification algorithms that could be implemented in an inexpensive computer system with limited amount of memory. Analysis of CCTs using these techniques can provide a great deal of detail permitting the use of the maximum resolution of LANDSAT imagery. Potential cases were detected in which the use of other techniques for classification using a Gaussian approximation for the distribution functions can be used with advantage. For jungle areas, channels 5 and 7 can provide enough information to delineate drainage patterns, swamp and wet areas, and make a reasonable broad classification of forest types.
A computational model for calculating the flow field, in the presence of solid particles, through a two dimensional compressor cascade is presented. The method treats the particle phase in the Lagrangian system and the fluid phase in the Eulerian system. The equations of momentum of the fluid phase are modified to account for the momentum exchange between the two phases. The resulting modified momentum and continuity equations are reduced to the conventional stream function-vorticity formulation. The analysis yields the change in the blade surface pressure distribution, the total pressure and velocity in the flow field. Based on this, it is possible to predict the decrease in compressor performance under the presence of particles.
The sounding rocket and satellite observations of space plasma waves within geospace in the frequency range from millihertz to megahertz are studied. Characteristic frequencies and source mechanisms of the plasma waves are described. The use of the Dynamic Explorer-1 Plasma Wave Instrument spectrograms to represent the plasma wave antenna and receiver system of geospace is examined. The ray tracing technique calculates the path of energy flow; the equations required for the analysis are presented. Cross-correlation of the wave electric and magnetic components provide data used to calculate the wave polarization, the direction of propagation, and the wave distribution function.
Chemiluminescent nightglow emission of the nitric oxide delta and gamma bands was measured from a sounding rocket launched on April 27, 1981, near local midnight. The integrated band emission rates for this near zenith observation above 205 km were less than 10 Rayleighs. The solar zenith angle was 127 deg. The branching ratio from the C2Pi state to the A2Sigma(+) state of NO was determined from comparison of the total emission rate of the delta band system to that of the gamma band system and found to be 0.30 + or - 0.06. The branching ratios within each of the band systems were found to be consistent with previous theoretical and experimental determinations. The vertical atomic nitrogen distribution, derived with the use of a model atmosphere, was found to have a peak density of 2.0 x 10 to the 7th atoms/cu cm at an altitude of 205 km. The analysis of these data indicate the presence of residual NO emission above 270 km at local midnight on the order of 1 Rayleigh of total band emission.
An error propagation model has been developed for multimodule computing systems in which the main parameters are the distribution functions of error propagation times. A digraph model is used to represent a multimodule computing system, and error propagation in the system is modeled by general distributions of error propagation times between all pairs of modules. Two algorithms are developed to compute systematically and efficiently the distributions of error propagation times. Experiments are also conducted to measure the distributions of error propagation times with the fault-tolerant microprocessor (FTMP). Statistical analysis of experimental data shows that the error propagation times in FTMP do not follow a well-known distribution, thus justifying the use of general distributions in the present model.
The ParaScope parallel programming environment, developed to support scientific programming of shared-memory multiprocessors, includes a collection of tools that use global program analysis to help users develop and debug parallel programs. This paper focuses on ParaScope's compilation system, its parallel program editor, and its parallel debugging system. The compilation system extends the traditional single-procedure compiler by providing a mechanism for managing the compilation of complete programs. Thus, ParaScope can support both traditional single-procedure optimization and optimization across procedure boundaries. The ParaScope editor brings both compiler analysis and user expertise to bear on program parallelization. It assists the knowledgeable user by displaying and managing analysis and by providing a variety of interactive program transformations that are effective in exposing parallelism. The debugging system detects and reports timing-dependent errors, called data races, in execution of parallel programs. The system combines static analysis, program instrumentation, and run-time reporting to provide a mechanical system for isolating errors in parallel program executions. Finally, we describe a new project to extend ParaScope to support programming in FORTRAN D, a machine-independent parallel programming language intended for use with both distributed-memory and shared-memory parallel computers.
Trace driven simulation is a well-established method of performance analysis for single processor computer systems. However, efficient and accurate memory address tracing for parallel computer systems is not well understood. In this paper we present a critical survey of recently implemented approaches to address tracing and highlight the issues specific to collection of traces for both shared and distributed memory parallel computers. These issues include potential distortion of the relative ordering of events by the address tracing activity, realistic interleaving of addresses generated by multiple processors, and I/O and storage problems associated with collecting traces for large parallel systems. The strengths and weaknesses of the parallel tracing approaches are described.
Flight flutter testing is an integral part of flight envelope clearance. This paper discusses advancements in several areas that are being investigated to improve efficiency and safety of flight test programs. Results are presented from recent flight testing of the F/A-18 Systems Research Aircraft. A wingtip excitation system was used to generate aeroelastic response data. This system worked well for many flight conditions but still displayed some anomalies. Wavelet processing is used to analyze the flight data. Filtered transfer functions are generated that greatly improve system identification. A flutter margin is formulated that accounts for errors between a model and flight data. Worst-case flutter margins are computed to demonstrate the flutter boundary may lie closer to the flight envelope than previously estimated. This paper concludes with developments for a distributed flight analysis environment and on-line health monitoring.
We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum exitation gap, gmin = O(n2(sup -n/2)), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.
Experimental, potentially real-time systems to detect floods and landslides related to heavy rain events are described. A key basis for these applications is high time resolution satellite rainfall analyses. Rainfall is the primary cause for devastating floods across the world. However, in many countries, satellite-based precipitation estimation may be the best source of rainfall data due to insufficient ground networks and absence of data sharing along many trans-boundary river basins. Remotely sensed precipitation from the NASA's TRMM Multi-satellite Precipitation Analysis (TMPA) operational system (near real-time precipitation at a spatial-temporal resolution of 3 hours and 0.25deg x 0.25deg) is used to monitor extreme precipitation events. Then these data are ingested into a macro-scale hydrological model which is parameterized using spatially distributed elevation, soil and land cover datasets available globally from satellite remote sensing. Preliminary flood results appear reasonable in terms of location and frequency of events, with implementation on a quasi-global basis underway. With the availability of satellite rainfall analyses at fine time resolution, it has also become possible to assess landslide risk on a near-global basis. Early results show that landslide occurrence is closely associated with the spatial patterns and temporal distribution of TRMM rainfall characteristics. Particularly, the number of landslides triggered by rainfall is related to rainfall climatology, antecedent rainfall accumulation, and intensity-duration of rainstorms. For the purpose of prediction, an empirical TMPA-based rainfall intensity-duration threshold is developed and shown to have skill in determining potential areas of landslides. These experimental findings, in combination with landslide surface susceptibility information based on satellite-based land surface information, form a starting point towards a potential operational landslide monitoring/warning system around the globe.