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 paper describes the process for defining and evaluating the costs of alternative transportation node concepts for Space Station Freedom. A cost estimating methodology including costs of development, production, delivery, assembly, steady state operations, and mission-specific operations is described in detail. The primary drivers of transportation node costs are identified and discussed.
This paper presents the experimental and theoretical results for expanding methane and ethylene diffusion flames in microgravity. A small porous sphere made from a low-density and low-heat-capacity insulating material was used to uniformly supply fuel at a constant rate to the expanding diffusion flame. A theoretical model which includes soot and gas radiation is formulated but only the problem pertaining to the transient expansion of the flame is solved by assuming constant pressure infinitely fast one-step ideal gas reaction and unity Lewis number. This is a first step toward quantifying the effect of soot and gas radiation on these flames. The theoretically calculated expansion rate is in good agreement with the experimental results. Both experimental and theoretical results show that as the flame radius increases, the flame expansion process becomes diffusion controlled and the flame radius grows as gamma t. Theoretical calculations also show that for a constant fuel mass injection rate a quasi-steady state is developed in the region surrounded by the flame and the mass flow rate at any location inside this region equals the mass injection rate.
A modification of an implicit approximate-factorization finite-difference algorithm applied to partial differential equations is presented. This algorithm is applied to the two- and three-dimensional Euler equations in general curvilinear coordinates. The modification transforms the coupled system of equations into an uncoupled diagonal form that requires less computational work. For steady-state applications, the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. The diagonal algorithm reduces the storage requirement of the implicit solution process and therefore has an important effect on the application of implicit finite-difference schemes to vector processors. Results are presented for realistic two-dimensional transonic flow fields about airfoils. Computation costs are reduced 24-34%.
An advanced computer code is undergoing development for numerically simulating laser cooling and trapping of large numbers of atoms. The code is expected to be useful in practical engineering applications and to contribute to understanding of the roles that light, atomic collisions, background pressure, and numbers of particles play in experiments using laser-cooled and -trapped atoms. The code is based on semiclassical theories of the forces exerted on atoms by magnetic and optical fields. Whereas computer codes developed previously for the same purpose account for only a few physical mechanisms, this code incorporates many more physical mechanisms (including atomic collisions, sub-Doppler cooling mechanisms, Stark and Zeeman energy shifts, gravitation, and evanescent-wave phenomena) that affect laser-matter interactions and the cooling of atoms to submillikelvin temperatures. Moreover, whereas the prior codes can simulate the interactions of at most a few atoms with a resonant light field, the number of atoms that can be included in a simulation by the present code is limited only by computer memory. Hence, the present code represents more nearly completely the complex physics involved when using laser-cooled and -trapped atoms in engineering applications. Another advantage that the code incorporates is the possibility to analyze the interaction between cold atoms of different atomic number. Some properties that cold atoms of different atomic species have, like cross sections and the particular excited states they can occupy when interacting with each other and light fields, play important roles not yet completely understood in the new experiments that are under way in laboratories worldwide to form ultracold molecules. Other research efforts use cold atoms as holders of quantum information, and more recent developments in cavity quantum electrodynamics also use ultracold atoms to explore and expand new information-technology ideas. These experiments give a hint on the wide range of applications and technology developments that can be tackled using cold atoms and light fields. From more precise atomic clocks and gravity sensors to the development of quantum computers, there will be a need to completely understand the whole ensemble of physical mechanisms that play a role in the development of such technologies. The code also permits the study of the dynamic and steady-state operations of technologies that use cold atoms. The physical characteristics of lasers and fields can be time-controlled to give a realistic simulation of the processes involved such that the design process can determine the best control features to use. It is expected that with the features incorporated into the code it will become a tool for the useful application of ultracold atoms in engineering applications. Currently, the software is being used for the analysis and understanding of simple experiments using cold atoms, and for the design of a modular compact source of cold atoms to be used in future research and development projects. The results so far indicate that the code is a useful design instrument that shows good agreement with experimental measurements (see figure), and a Windows-based user-friendly interface is also under development.
The CMTB Thwacker hardware operates as a testbed analogue for the Flight Thwacker and Sieve components of CHIMRA, a device on the Curiosity Rover. The sieve separates particles with a diameter smaller than 150 microns for delivery to onboard science instruments. The sieving behavior of the testbed hardware should be similar to the Flight hardware for the results to be meaningful. The elastodynamic behavior of both sieves was studied analytically using the Rayleigh Ritz method in conjunction with classical plate theory. Finite element models were used to determine the mode shapes of both designs, and comparisons between the natural frequencies and mode shapes were made. The analysis predicts that the performance of the CMTB Thwacker will closely resemble the performance of the Flight Thwacker within the expected steady state operating regime. Excitations of the testbed hardware that will mimic the flight hardware were recommended, as were those that will improve the efficiency of the sieving process.
Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.
Numerical simulations of the MHD equations for a fully compressible, low-beta, resistive plasma are used to study the resonance absorption process for the heating of coronal active region loops. Comparisons with more approximate analytic models show that the major predictions of the analytic theories are, to a large extent, confirmed by the numerical computations. The simulations demonstrate that the dissipation occurs primarily in a thin resonance layer. Some of the analytically predicted features verified by the simulations are (a) the position of the resonance layer within the initial inhomogeneity; (b) the importance of the global mode for a large range of loop densities; (c) the dependence of the resonance layer thickness and the steady-state heating rate on the dissipation coefficient; and (d) the time required for the resonance layer to form. In contrast with some previous analytic and simulation results, the time for the loop to reach a steady state is found to be the phase-mixing time rather than a dissipation time. This disagreement is shown to result from neglect of the existence of the global mode in some of the earlier analyses. The resonant absorption process is also shown to behave similar to a classical driven harmonic oscillator.
A simplified analytical method for calculation of thermal response within a transpiration-cooled porous heat shield material in an intense radiative-convective heating environment is presented. The essential assumptions of the radiative and convective transfer processes in the heat shield matrix are the two-temperature approximation and the specified radiative-convective heatings of the front surface. Sample calculations for porous silica with CO2 injection are presented for some typical parameters of mass injection rate, porosity, and material thickness. The effect of these parameters on the cooling system is discussed.
The deviation from a Maxwellian velocity distribution caused by excitation of the J = 2 rotational level of parahydrogen (followed by radiative decay) in a molecular hydrogen gas is investigated. It is noted that inelastic collisions deplete the high-velocity tail of the distribution, while elastic collisions tend to refill it; the resulting steady-state distribution has a small depletion in the tail which slightly reduces the cooling rate of the gas. The elastic-collision transport is approximated by a continuous process, a simplified expression is obtained for the inelastic-collision operator, and the Boltzmann equation is then solved analytically. A correction to the inelastic-collision rate is evaluated, and it is found that the relative magnitude of the effect of tail depletion on the inelastic-collision rate is strongly temperature-dependent. The critical molecular density above which the J = 2 level deexcites by superelastic collisions is shown to be a weak function of temperature.
S-band time delay measurements were collected from the spacecrafts Helios A and B during three solar occultations in 1975/76 within heliocentric distances of about 3 and 215 solar radii in terms of range, Doppler frequency shift, and electron content. A description is given concerning some characteristic features of the methods of measurement and data processing. Typical data sets are discussed to probe the electron density distribution near the sun (west and east limb as well) including the outer and extended corona. Steady-state and dynamical aspects of the solar corona are presented and compared with earth-bound K-coronagraph measurements. Using a weighted least squares estimation 3 parameters of an average coronal electron density profile are derived in a preliminary analysis to yield electron densities of about 130 billion; 100 million; 7 million/cu m at r?3; 65; 215 solar radii. Transient phenomena are discussed and a velocity of propagation v approximately 900 km/s for plasma ejecta from a solar flare is determined from an extraordinary set of Helios B electron content measurements on April 30/May 1, 1976.
The magnetic reconnection process is studied here using incompressible MHD simulations with different inflow boundary conditions and different magnetic Reynolds numbers R(m). The angle between the magnetic separatrices is in steady state reconnection depends mainly on the normal magnetic field on the inflow boundary. In steady state nonuniform reconnection with large separatrix angles, field-aligned plasma jets appear slightly downstream of the magnetic separatrices. The field-aligned plasma jet are stronger when R(m) is larger. Each field-aligned plasma jet consists of two parts: a slow shock and a fast-mode compressional wave. The slow shock converts the magnetic energy into plasma kinetic energy by acceleration and heating. The fast-mode compressional wave decelerates the plasma to a smaller outflow speed and heats it further. Nearly all the magnetic energy flowing into the diffusion region is converted into other forms. The length and width of the diffusion region depend on the values of the reconnection rate, R(m), and the normal magnetic field on the inflow boundary.
The goals of our work on detached solidification have been to: Develop a complete understanding of all of the phenomena of detached solidification; Make it possible to achieve detached solidification reproducibly; Increase crystallographic perfection through detached solidification. Specifically, we aimed in this project to: Identify a system and develop methods that would allow viewing of the melt surface and convection in the melt during detached solidification in microgravity; Improve understanding of the origination and evolution of detachment through experiments and theoretical treatments; Achieve detachment on earth. The project resulted in 14 publications, 15 presentations, completion of 2 Ph.D. theses, and completion of 2 M.S. theses. Two additional papers are currently being reviewed for publication. Copies of most of the papers are attached as appendices. Among the accomplishments are: Achievement of detached solidification of InSb on earth and determination of the conditions favoring detachment on earth. Development of a new method for coating the interior of silica growth ampoules with transparent boron nitride, which yields the high contact angles for semiconductor melts that favor detachment. Development of a new coating for the interior of Pyrex ampoules yielding very high contact angles for water and molten organic compounds. Development of a material-balance model for steady detached solidification that provides greatly improved insight into the process, both in microgravity and on earth. Dimensionless parameters were found that clarify the role of operating conditions and physical properties. The reason for the occurrence of two steady states was clarified, along with their relative stability. Observation of periodic gas tubes during directional solidification of water and organic compounds. Gas bubbles did not propagate around the periphery to yield full detachment. Failure to obtain detachment is attributed to the ready plastic deformation of these materials, so that the frozen material continues to adhere to the ampoule wall during cooling.
A simplified analytical solution for thermal response of a transpiration-cooled porous heat-shield material in an intense radiative-convective heating environment is presented. Essential features of this approach are "two-flux method" for radiative transfer process and "two-temperature" assumption for solid and gas temperatures. Incident radiative-convective heatings are specified as boundary conditions. Sample results are shown using porous silica with CO2 transpiration and some parameters quantitatively show the effect on this transpiration cooling system. Summarized maps for mass injection rate, porosity and blowing correction factor for radiation are obtained in order to realize such a cooling system.
The presence of unthermalized photoelectrons in the sunlit polar cap leads to an enhanced ambipolar potential drop and enhanced upward ion acceleration. Observations in the topside ionosphere have led to the conclusion that large-scale electrostatic potential drops exist above the spacecraft along polar magnetic field lines connected to regions of photoelectron production. A kinetic approach is used for the O(+), H(+), and photoelectron (p) distributions, while a fluid approach is used to describe the thermal electrons (e) and self-consistent electric field (E(sub II)) electrons are allowed to carry a flux that compensates for photoelectron escape, a critical assumption. Collisional processes are excluded, leading to easier escape of polar wind particles and therefore to the formation of the largest potential drop consistent with this general approach. We compute the steady state electric field enhancement and net potential drop expected in the polar wind due to the presence of photoelectrons as a function of the fractional photoelectron content and the thermal plasma characteristics. For a set of low-altitude boundary conditions typical of the polar wind ionosphere, including 0.1% photoelectron content, we found a potential drop from 500 km to 5 R(sub E) of 6.5 V and a maximum thermal electron temperature of 8800 K. The reasonable agreement of our results with the observed polar wind suggests that the assumptions of this approach are valid.
The RIDGE effort continues the aegis of the earlier, NASA-sponsored, Isothermal Dendritic Growth Experiment (IDGE) series of experiments through the continued analysis of microgravity data acquired during these earlier space flights. The preliminary observations presented here demonstrate that there are significant differences between SCN and the more anisotropic PVA dendrites. The side branch structure becomes amplified only further behind the tip, and the interface shape is generally wider (i.e. more hyperbolic than parabolic) in PVA than in SCN. These characteristics are seen to affect the process of heat transport. Additionally, the dendrites grown during the fourth United States Microgravity Payload (USMP-4) exhibit time-dependent growth characteristics and may not always have reached steady-state growth during the experiment.
The atmospheric structure of Triton in the altitude range 25-150 kilometers shows an unexpectedly steep thermal gradient of 0.26 K per kilometer above 50 kilometer altitude, with a nearly isothermal profile below. The upper part of the profile can be explained by downward conduction of heat deposited by magnetospheric electrons and solar UV. However, the atmospheric temperature below 50 kilometers is too cold for identified radiative processes to dispose of the inferred heat flux (0.0012 erg per square centimeter per second) from the upper atmosphere. This implies that either the atmosphere is not in a steady state and/or an unidentified cooling mechanism is at work in the altitude range 25-50 kilometers. When extrapolated to the surface, the inversion results yield a pressure of 19.0 sup (+1.8) sub (-1.5), mubar, about 5mubar greater than that observed by Voyager.
A thermal evaluation of a composite tank wall design for a liquid hydrogen tank was performed in the present study. The primary focus of the current effort was to perform one-dimensional, temperature nonlinear, transient thermal analyses to determine the through-the-thickness temperature profiles. These profiles were used to identify critical points within the flight envelope that could have detrimental effects on the adhesive bondlines used in the construction of the tank wall. Additionally, this paper presents the finite element models, analysis strategies, and thermal analysis results that were determined for several vehicle flight conditions. The basic tank wall configuration used to perform the thermal analyses consisted of carbon-epoxy facesheets and a Korex honeycomb core sandwich that was insulated with an Airex cryogenic foam and an Alumina Enhanced Thermal Barrier (AETB-12). Nonlinear, transient thermal analyses were conducted using the ABAQUS finite element code. Tank wall models at a windward side location on the fuel tank were analyzed for three basic flight conditions: cold-soak (ground-hold), ascent, and re-entry. Additionally, three ambient temperature boundary conditions were applied to the tank wall for the cold-soak condition, which simulated the launch pad cooldown process. Time-dependent heating rates were used in the analyses of the ascent and reentry segments of the flight history along with temperature dependent material properties. The steady-state through-the-thickness temperature profile from the cold-soak condition was used as the initial condition for the ascent analyses. Results from the nonlinear thermal analyses demonstrated very good correlation with results from similar models evaluated by Northrop- Grumman using a different analysis tool. Wall through-the-thickness temperature gradients as a function of flight time were obtained for future incorporation into a full-scale thermostructural analysis to evaluate the adhesive bondlines. As a result of the thermal analyses conducted, a sufficient level of confidence was demonstrated in the thermal modeling and analysis capabilities of ABAQUS to warrant future use as a thermo-structural analysis tool to evaluate cryogenic tank wall designs.
The dynamics and energetics of intergranular crack growth along a flat grain boundary in aluminum is studied by a molecular-dynamics simulation model for crack propagation under steady-state conditions. Using the ability of the molecular-dynamics simulation to identify atoms involved in different atomistic mechanisms, it was possible to identify the energy contribution of different processes taking place during crack growth. The energy contributions were divided as: elastic energy, defined as the potential energy of the atoms in fcc crystallographic state; and plastically stored energy, the energy of stacking faults and twin boundaries; grain-boundary and surface energy. In addition, monitoring the amount of heat exchange with the molecular-dynamics thermostat gives the energy dissipated as heat in the system. The energetic analysis indicates that the majority of energy in a fast growing crack is dissipated as heat. This dissipation increases linearly at low speed, and faster than linear at speeds approaching 1/3 the Rayleigh wave speed when the crack tip becomes dynamically unstable producing periodic dislocation bursts until the crack is blunted.