Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal diffusion”

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.

At least 685 records · Page 38

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

A 2D Kaleidoscope of Electron Heat Fluxes Driven by Auroral Electron Precipitation

Electron heat flux is an important value for ionospheric space weather modeling networks. Utilizing the 2D array of Time History of Events and Macroscale Interactions during Substorms all-sky-imager (ASI) observations, Gabrielse et al. (2021, https://doi.org/10.3389/fphy.2021.744298) described a new method that estimates the auroral scale sizes of intense precipitating electron energy fluxes and their mean energies during two substorms on 16 February 2010. These parameters in combination with SuperThermal Electron Transport code were used to develop a new methodology to calculate electron thermal fluxes from data inputs in 2D during one of the substorms at 09:40:00 UT across Canada and Alaska. To test the effect of various precipitation lifetimes on electron heat flux values, boxcar averages ranging from 0 to 900 s were applied to the ASI data. These data are then combined with the newly developed kinetic simulation to determine the thermal fluxes associated with the observed diffuse and discrete precipitation.

Electron precipitation in aurora↗

Cleaning of a thermal vacuum chamber with shrouds in place

In February, 1991, a failure of a rotary booster pump caused the diffusion pumps to backstream into a 10 ft x 15 ft thermal vacuum chamber. Concerns existed about the difficulty of removing and reinstalling the shrouds without causing leaks. The time required for the shroud removal was also of concern. These concerns prompted us to attempt to clean the chamber without removing the shrouds.

Bond, William R.↗

Cellular interface morphologies in directional solidification. III - The effects of heat transfer and solid diffusivity

The shape and stability of two-dimensional finite-amplitude cellular interfaces arising during directional solidification are compared for several solidification models that account differently for latent heat released at the interface, unequal thermal conductivities of melt and solid, and solute diffusivity in the solid. Finite-element analysis and computer-implemented perturbation methods are used to analyze the families of steadily growing cellular forms that evolve from the planar state. In all models a secondary bifurcation between different families of finite-amplitude cells exists that halves the spatial wavelength of the stable interface. The quantitative location of this transition is very dependent on the details of the model. Large amounts of solute diffusion in the solid retard the growth of large-amplitude cells.

Ungar, Lyle H.↗

Thermal Skin fabrication technology

Advanced fabrication techniques applicable to Thermal Skin structures were investigated, including: (1) chemical machining; (2) braze bonding; (3) diffusion bonding; and (4) electron beam welding. Materials investigated were nickel and nickel alloys. Sample Thermal Skin panels were manufactured using the advanced fabrication techniques studied and were structurally tested. Results of the program included: (1) development of improved chemical machining processes for nickel and several nickel alloys; (2) identification of design geometry limits; (3) identification of diffusion bonding requirements; (4) development of a unique diffusion bonding tool; (5) identification of electron beam welding limits; and (6) identification of structural properties of Thermal Skin material.

Milam, T. B.↗

Partially coherent scattering in stellar chromospheres. I - Effects on resonance line thermalization

The conditions necessary for partially coherent scattering to influence the thermalization depth are investigated, and the effects of partial redistribution in homogeneous slab atmospheres are approximated using several common resonance lines as examples. For electron densities above 10 exp 10/cu cm, it is concluded that even when coherent scattering dominates the escape process, the thermalization depths of the strongest resonance lines of H, Ca II, and Mg II agree roughly with the Doppler diffusion in frequency. No elastic scattering effects on the thermalization of Mg II h and k and Ca II H and K are found. At lower densities, such as for giant star chromospheres and QSO models, the results deviate strongly from the Doppler complete redistribution case.

Gayley, Kenneth G.↗

Pulsed-field-gradient measurements of time-dependent gas diffusion

Pulsed-field-gradient NMR techniques are demonstrated for measurements of time-dependent gas diffusion. The standard PGSE technique and variants, applied to a free gas mixture of thermally polarized xenon and O2, are found to provide a reproducible measure of the xenon diffusion coefficient (5.71 x 10(-6) m2 s-1 for 1 atm of pure xenon), in excellent agreement with previous, non-NMR measurements. The utility of pulsed-field-gradient NMR techniques is demonstrated by the first measurement of time-dependent (i.e., restricted) gas diffusion inside a porous medium (a random pack of glass beads), with results that agree well with theory. Two modified NMR pulse sequences derived from the PGSE technique (named the Pulsed Gradient Echo, or PGE, and the Pulsed Gradient Multiple Spin Echo, or PGMSE) are also applied to measurements of time dependent diffusion of laser polarized xenon gas, with results in good agreement with previous measurements on thermally polarized gas. The PGMSE technique is found to be superior to the PGE method, and to standard PGSE techniques and variants, for efficiently measuring laser polarized noble gas diffusion over a wide range of diffusion times. Copyright 1998 Academic Press.

NASA Discipline Environmental Health↗

Transport coefficients of partially ionized argon.

Properties of partially ionized Ar computed via Chapman-Enskog-Burnett expressions noting ambipolar diffusion coefficient, electron-atom momentum, electrical and thermal conductivity, etc

TRANSPORT COEFFICIENT↗

Diffusive shock acceleration - Comparison of a unified shock model to bow shock observations

A comparison is made between recent AMPTE/IRM observations of diffuse ions detected upstream of the earth's bow shock when the interplanetary magnetic field was nearly parallel to the solar wind direction and a known collisionless quasi-parallel shock model. These observations, which provide the proton spectrum for all of velocity space, give a direct measure of the shock acceleration efficiency and show how thermal solar winds are injected into the diffusive shock acceleration mechanism. The model accurately describes the proton spectrum; moreover, it predicts the shock structure, the complete particle spectrum, the relative velocity distributions of different ion species in the down stream region, and the enhancement of heavy ions over protons relative to the unshocked solar wind.

Ellison, Donald C.↗

Flight Computer Processing Avionics for Space Station Microgravity Experiments: A Risk Assessment of Commercial Off-the-Shelf Utilization

NASA/Marshall Space Flight Center (MSFC) is continually looking for ways to reduce the costs and schedule and minimize the technical risks during the development of microgravity programs. One of the more prominent ways to minimize the cost and schedule is to use off-the-shelf hardware (OTS). However, the use of OTS often increases the risk. This paper addresses relevant factors considered during the selection and utilization of commercial off-the-shelf (COTS) flight computer processing equipment for the control of space station microgravity experiments. The paper will also discuss how to minimize the technical risks when using COTS processing hardware. Two microgravity experiments for which the COTS processing equipment is being evaluated for are the Equiaxed Dendritic Solidification Experiment (EDSE) and the Self-diffusion in Liquid Elements (SDLE) experiment. Since MSFC is the lead center for Microgravity research, EDSE and SDLE processor selection will be closely watched by other experiments that are being designed to meet payload carrier requirements. This includes the payload carriers planned for the International Space Station (ISS). The purpose of EDSE is to continue to investigate microstructural evolution of, and thermal interactions between multiple dendrites growing under diffusion controlled conditions. The purpose of SDLE is to determine accurate self-diffusivity data as a function of temperature for liquid elements selected as representative of class-like structures. In 1999 MSFC initiated a Center Director's Discretionary Fund (CDDF) effort to investigate and determine the optimal commercial data bus architecture that could lead to faster, better, and lower cost data acquisition systems for the control of microgravity experiments. As part of this effort various commercial data acquisition systems were acquired and evaluated. This included equipment with various form factors, (3U, 6U, others) and equipment that utilized various bus structures, (VME, PC104, STD bus). This evaluation of hardware was performed in conjunction with a trade study that considered over twenty (20) different factors relevant to the selection of an optimum design approach. These factors included; safety, sizing and timing, radiation hardness and single event upset, power consumption, heat dissipation, size and volume, expected service life, maintainability, heritage, operating systems, requirements for software reuse, availability of compatible interface boards, relative cost, schedule, reliability, EMI/EMC factors, "hot swap" capability, standards for conduction cooling, I/O capabilities, unique carrier requirements and operating system considerations. The approach to evaluate Safety as part of this study included a review of the Preliminary Hazard Analysis (PHA) for each of the experiment designs and a determination of how each hazard could be addressed and eliminated when different processors were selected. This included evaluating various design approaches and trade-offs between fault tolerant designs and fail-safe designs in accordance with NSTS 1700.7B. This will include the results of radiation testing where available. Various operating systems, such as VxWorks, Linux, QNX, and Embedded NT are evaluated and the advantages and disadvantages of their utilization are also addressed. Design implementation strategies for the various operating systems are considered and discussed. This paper presents the results and recommendations from this trade study. Preliminary conclusions from this study are that safety concerns from lack or radiation testing on COTS equipment can be addressed by additional testing and design considerations, the PC104 bus provided adequate I/O for the SDLE and EDSE microgravity experiments, and PC104 bus components offered significant advantages over VME and cPCI for weight and space reductions.

Estes, Howard↗

Nonlinear Transient Thermal Analysis by the Force-Derivative Method

High-speed vehicles such as the Space Shuttle Orbiter must withstand severe aerodynamic heating during reentry through the atmosphere. The Shuttle skin and substructure are constructed primarily of aluminum, which must be protected during reentry with a thermal protection system (TPS) from being overheated beyond the allowable temperature limit, so that the structural integrity is maintained for subsequent flights. High-temperature reusable surface insulation (HRSI), a popular choice of passive insulation system, typically absorbs the incoming radiative or convective heat at its surface and then re-radiates most of it to the atmosphere while conducting the smallest amount possible to the structure by virtue of its low diffusivity. In order to ensure a successful thermal performance of the Shuttle under a prescribed reentry flight profile, a preflight reentry heating thermal analysis of the Shuttle must be done. The surface temperature profile, the transient response of the HRSI interior, and the structural temperatures are all required to evaluate the functioning of the HRSI. Transient temperature distributions which identify the regions of high temperature gradients, are also required to compute the thermal loads for a structural thermal stress analysis. Furthermore, a nonlinear analysis is necessary to account for the temperature-dependent thermal properties of the HRSI as well as to model radiation losses.

Balakrishnan, Narayani V.↗

Investigation of Aeronomic Processes on the Basis of Ionospheric Sporadic E

A method by means of which the characteristics of atmospheric turbulence can be determined is discussed. Turbulent diffusivity deduced from ionospheric sporadic E in different altitudes above the ionospheric stations at Bekescsaba and Juliusrah are given. Seasonal variations in different altitudes above Bekescsaba are shown. An attempt was made to determine the variation of the turbulent diffusivity during and after geomagnetic disturbances by computing the horizontal wind velocity from the thermal wind equation. It is argued that the change of the turbulent diffusivity below the turbopause is opposed to that above the turbopause.

Bencze, P.↗

Cosmic-ray-modified stellar winds. II - A perturbation approach

A perturbation method is developed to describe the modification of a stellar wind flow with termination shock by the Galactic cosmic rays. The perturbation parameter epsilon used in the analysis is the ratio of the galactic cosmic-ray pressure to the thermal gas pressure at a large distance from the star. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, particularly just upstream of the termination shock of the wind. A one-fluid polytropic model is used to describe the thermal stellar wind gas, with the coupling between the cosmic rays and the thermal gas being determined by an average hydrodynamical diffusion coefficient kappa describing the scattering of cosmic rays by magnetohydrodynamic waves traveling in the background flow. The analysis takes into account cosmic-ray modifications of the critical point of the wind, the thermal gas entropy constants on both sides of the shock, and the fluid velocity profile.

Ko, C. M.↗

Model calculations for diffuse molecular clouds

A steady state isobaric cloud model is developed. The pressure, thermal, electrical, and chemical balance equations are solved simultaneously with a simple one dimensional approximation to the equation of radiative transfer appropriate to diffuse clouds. Cooling is mainly by CII fine structure transitions, and a variety of heating mechanisms are considered. Particular attention is given to the abundance variation of H2. Inhomogeneous density distributions are obtained because of the attenuation of the interstellar UV field and the conversion from atomic to molecular hyrodgen. The effects of changing the model parameters are described and the applicability of the model to OAO-3 observations is discussed. Good qualitative agreement with the fractional H2 abundance determinations has been obtained. The observed kinetic temperatures near 80 K can also be achieved by grain photoelectron heating. The problem of the electron density is solved taking special account of the various hydrogen ions as well as heavier ones.

Glassgold, A. E.↗

Modeling Specular Exchange Between Concentric Cylinders in a Radiative Shielded Furnace

The objective of this research is to develop and validate mathematical models to characterize the thermal performance of a radiative shielded furnace, the University of Alabama in Huntsville (UAH) Isothermal Diffusion Oven. The mathematical models are validated against experimental data obtained from testing the breadboard oven in a terrestrial laboratory environment. It is anticipated that the validation will produce math models capable of predicting the thermal performance of the furnace over a wide range of operating conditions, including those for which no experimental data is available. Of particular interest is the furnace core temperature versus heater power parametric and the transient thermal response of the furnace. Application to a microgravity environment is not considered, although it is conjectured that the removal of any gravity dependent terms from the math models developed for the terrestrial application should yield adequate results in a microgravity environment. The UAH Isothermal Diffusion Oven is designed to provide a thermal environment that is conducive to measuring the diffusion of high temperature liquid metals. In addition to achieving the temperatures required to melt a sample placed within the furnace, reducing or eliminating convective motions within the melt is an important design consideration [1]. Both of these influences are reflected in the design of the furnace. Reducing unwanted heat losses from the furnace is achieved through the use of low conductivity materials and reflective shielding. As evidenced by the highly conductive copper core used to house the sample within the furnace, convective motions can be greatly suppressed by providing an essentially uniform thermal environment. An oven of this design could ultimately be utilized in a microgravity environment, presumably as a experiment payload. Such an application precipitates other design requirements that limit the resources available to the furnace such as power, mass, volume, and possibly even time. Through the experimental and numerical results obtained, the power requirements and thermal response time of the breadboard furnace are quantified.

Schunk, Richard Gregory↗

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗