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.
-Forced invariance- is a processing method that can subdue the expression of vegetation and enhance the expression of the underlying lithology in remotely sensed imagery.
This paper describes fundamental soils characterization work performed at NASA Langley Research Center in support of the Subsonic Rotary Wing (SRW) Aeronautics Program and the Orion Landing System (LS) Advanced Development Program (ADP). LS-DYNA(Registered TradeMark)1 soil impact model development and test-analysis correlation results are presented for: (1) a 38-ft/s vertical drop test of a composite fuselage section, outfitted with four blocks of deployable energy absorbers (DEA), onto sand, and (2) a series of impact tests of a 1/2-scale geometric boilerplate Orion capsule onto soil. In addition, the paper will discuss LS-DYNA contact analysis at the soil/structure interface, methods used to estimate frictional forces, and the sensitivity of the model to density, moisture, and compaction.
NASA's Energy and Water Cycle Study (NEWS) program fosters collaborative research towards improved quantification and prediction of water and energy cycle consequences of climate change. In order to measure change, it is first necessary to describe current conditions. The goal of the NEWS Water and Energy Cycle Climatology project is to develop "state of the global water cycle" and "state of the global energy cycle" assessments based on data from modern ground and space based observing systems and data integrating models. The project is a multiinstitutional collaboration with more than 20 active contributors. This presentation will describe results of the first stage of the water budget analysis, whose goal was to characterize the current state of the water cycle on mean monthly, continental scales. We examine our success in closing the water budget within the expected uncertainty range and the effects of forcing budget closure as a method for refining individual flux estimates.
A lightweight bulldozer blade prototype has been designed and built to be used as an excavation implement in conjunction with the NASA Chariot lunar mobility platform prototype. The combined system was then used in a variety of field tests in order to characterize structural loads, excavation performance and learn about the operational behavior of lunar excavation in geotechnical lunar simulants. The purpose of this effort was to evaluate the feasibility of lunar excavation for site preparation at a planned NASA lunar outpost. Once the feasibility has been determined then the technology will become available as a candidate element in the NASA Lunar Surface Systems Architecture. In addition to NASA experimental testing of the LANCE blade, NASA engineers completed analytical work on the expected draft forces using classical soil mechanics methods. The Colorado School of Mines (CSM) team utilized finite element analysis (FEA) to study the interaction between the cutting edge of the LANCE blade and the surface of soil. FEA was also used to examine various load cases and their effect on the lightweight structure of the LANCE blade. Overall it has been determined that a lunar bulldozer blade is a viable technology for lunar outpost site preparation, but further work is required to characterize the behavior in 1/6th G and actual lunar regolith in a vacuum lunar environment.
Direct to Earth return links are limited by the size and power of lander devices. A standard alternative is provided by a two-hops return link: a proximity link (from lander to orbiter relay) and a deep-space link (from orbiter relay to Earth). Although direct to Earth return links are limited by the size and power of lander devices, using an additional link and a proposed coding for relay channels, one can obtain a more reliable signal. Although significant progress has been made in the relay coding problem, existing codes must be painstakingly optimized to match to a single set of channel conditions, many of them do not offer easy encoding, and most of them do not have structured design. A high-performing LDPC (low-density parity-check) code for the relay channel addresses simultaneously two important issues: a code structure that allows low encoding complexity, and a flexible rate-compatible code that allows matching to various channel conditions. Most of the previous high-performance LDPC codes for the relay channel are tightly optimized for a given channel quality, and are not easily adapted without extensive re-optimization for various channel conditions. This code for the relay channel combines structured design and easy encoding with rate compatibility to allow adaptation to the three links involved in the relay channel, and furthermore offers very good performance. The proposed code is constructed by synthesizing a bilayer structure with a pro to graph. In addition to the contribution to relay encoding, an improved family of protograph codes was produced for the point-to-point AWGN (additive white Gaussian noise) channel whose high-rate members enjoy thresholds that are within 0.07 dB of capacity. These LDPC relay codes address three important issues in an integrative manner: low encoding complexity, modular structure allowing for easy design, and rate compatibility so that the code can be easily matched to a variety of channel conditions without extensive re-optimization. The main problem of half-duplex relay coding can be reduced to the simultaneous design of two codes at two rates and two SNRs (signal-to-noise ratios), such that one is a subset of the other. This problem can be addressed by forceful optimization, but a clever method of addressing this problem is via the bilayer lengthened (BL) LDPC structure. This method uses a bilayer Tanner graph to make the two codes while using a concept of "parity forwarding" with subsequent successive decoding that removes the need to directly address the issue of uneven SNRs among the symbols of a given codeword. This method is attractive in that it addresses some of the main issues in the design of relay codes, but it does not by itself give rise to highly structured codes with simple encoding, nor does it give rate-compatible codes. The main contribution of this work is to construct a class of codes that simultaneously possess a bilayer parity- forwarding mechanism, while also benefiting from the properties of protograph codes having an easy encoding, a modular design, and being a rate-compatible code.
Modeling of multicomponent trace contaminant adsorption on activated carbon for spacecraft life support applications presents challenges due to the very low contaminant concentrations in cabin air coupled with very strong adsorption of heavier components such as cyclic polydimethylsiloxanes. Both gas-phase and solid-phase micropore mass transfer resistances can be important, and competitive interactions between adsorbed components can result in behavior such as roll-up where lighter components displaced by heavier components have a higher concentration at the outlet than at the inlet. This paper describes the further development of a dynamic trace contaminant adsorption model based on Ideal Adsorbed Solution Theory (IAST). Alternative 1-dimensional and 2-dimensional versions of the model are described with different mass transfer resistance and driving force assumptions. Optimization of solution methods for speed and stability is also described. Mass balance equations are scaled to provide sensitivity over the wide range of component concentrations in cabin air. Model parameters are estimated from isotherm and breakthrough test data on single adsorbate and binary adsorbate systems. Multicomponent results are compared against predictions of a heritage model used extensively by NASA for trace contaminant control system (TCCS) carbon bed sizing. Finally, efforts to implement the model in a second software platform for wider access are described.
A method is suggested which determines the monthly mean longwave radiation at the surface over the oceans using currently available satellite data sets. An independent estimate of the frequency of occurrence of clouds or cloud top heights is not required. A climatological relationship between the LWCRF (longwave cloud radiation forcing) at the top of the atmosphere and surface may be substituted. This ratio is obtained from a GMC simulation but in principle could be obtained observationally, although global coverage is doubtful. A flow diagram shows the procedure used to obtain maps for the monthly mean net upward longwave flux at the surface.
When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.
This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.
This document outlines the tests performed to make aerodynamic force and torque measurements on the SOFIA wind tunnel model telescope. These tests were performed during the SOFIA 2 wind tunnel test in the 14 ft wind tunnel during the months of June through August 1994. The test was designed to measure the dynamic cross elevation moment acting on the SOFIA model telescope due to aerodynamic loading. The measurements were taken with the telescope mounted in an open cavity in the tail section of the SOFIA model 747. The purpose of the test was to obtain an estimate of the full scale aerodynamic disturbance spectrum, by scaling up the wind tunnel results (taking into account differences in sail area, air density, cavity dimension, etc.). An estimate of the full scale cross elevation moment spectrum was needed to help determine the impact this disturbance would have on the telescope positioning system requirements. A model of the telescope structure, made of a light weight composite material, was mounted in the open cavity of the SOFIA wind tunnel model. This model was mounted via a force balance to the cavity bulkhead. Despite efforts to use a 'stiff' balance, and a lightweight model, the balance/telescope system had a very low resonant frequency (37 Hz) compared to the desired measurement bandwidth (1000 Hz). Due to this mechanical resonance of the balance/telescope system, the balance alone could not provide an accurate measure of applied aerodynamic force at the high frequencies desired. A method of measurement was developed that incorporated accelerometers in addition to the balance signal, to calculate the aerodynamic force.
A method has been developed for real-time stability margin measurement calculations. The method relies on a tailored-forced excitation targeted to a specific frequency range. Computation of the frequency response is matched to the specific frequencies contained in the excitation. A recursive Fourier transformation is used to make the method compatible with real-time calculation. The method was incorporated into the X-38 nonlinear simulation and applied to an X-38 robustness test. X-38 stability margins were calculated for different variations in aerodynamic and mass properties over the vehicle flight trajectory. The new method showed results comparable to more traditional stability analysis techniques, and at the same time, this new method provided coverage that is more complete and increased efficiency.
In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.
A method of successive approximations is employed to solve the problem of heat transfer to a transparent gas flowing through a radiating-conducting tube with turbulent forced convection between the tube wall and the gas, and with energy generation in the wall. Emphasis is given to the effect of emissivity of the wall to the tube and gas temperature profiles.
The force required to pull 7t and 40t polyhedral boulders from the surface of an asteroid is simulated using the discrete element method considering the effects of microgravity, regolith cohesion and boulder acceleration. The connection between particle surface energy and regolith cohesion is estimated by simulating a cohesion sample tearing test. An optimal constant acceleration is found where the peak net force from inertia and cohesion is a minimum. Peak pulling forces can be further reduced by using linear and quadratic acceleration functions with up to a 40% reduction in force for quadratic acceleration.
Low-speed experimental force and data on a series of thin swept wings with sharp leading edges and leading and trailing-edge flaps are compared with predictions made using a linearized-theory method which includes estimates of vortex forces. These comparisons were made to assess the effectiveness of linearized-theory methods for use in the design and analysis of flap systems in subsonic flow. Results demonstrate that linearized-theory, attached-flow methods (with approximate representation of vortex forces) can form the basis of a rational system for flap design and analysis. Even attached-flow methods that do not take vortex forces into account can be used for the selection of optimized flap-system geometry, but design-point performance levels tend to be underestimated unless vortex forces are included. Illustrative examples of the use of these methods in the design of efficient low-speed flap systems are included.
A review of analytical methods used to analyze flexible vehicles with active controls is presented. Methods used to approximate and correct unsteady aerodynamic forces used in the analysis are discussed. Recent advances in the application of optimal methods to digital control law synthesis are presented. The use of active controls in an integrated design process is also discussed. Finally, the results of recent wind-tunnel studies aimed at demonstrating active control concepts and validating the analytical methods are presented.
Closed form analytical solutions are derived for the vibration and wave propagation of multisegmented elastic beams. Each segment is modeled as a Timoshenko beam with possible inclusion of material viscosity, elastic foundation and axial forces. Solutions are obtained by using transfer matrix methods. According to these methods formal solutions are first constructed which relate the deflection, slope, moment and shear force of one end of the individual segment to those of the other. By satisfying appropriate continuity conditions at segment junctions, a global 4x4 matrix results which relates the deflection, slope, moment and shear force of one end of the beam to those of the other. If any boundary conditions are subsequently invoked on the ends of the beam one gets the appropriate characteristic equation for the natural frequencies. Furthermore, by invoking appropriate periodicity conditions the dispersion relation for a periodic system is obtained. A variety of numerical examples are included.
Because the potential flow suction along the leading and side edges of a planform can be used to determine both leading- and side-edge vortex lift, the present investigation was undertaken to apply the vortex-lattice method to computing side-edge suction force for isolated or interacting planforms. Although there is a small effect of bound vortex sweep on the computation of the side-edge suction force, the results obtained for a number of different isolated planforms produced acceptable agreement with results obtained from a method employing continuous induced-velocity distributions. By using the method outlined, better agreement between theory and experiment was noted for a wing in the presence of a canard than was previously obtained.