Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium”

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 199 records · Page 11

Orientation and velocity dependence of the nonequilibrium partition coefficient

Monte Carlo simulations based on a Spin-1 Ising Model for binary alloys have been used to investigate the non-equilibrium partition coefficient (k(sub neq)) as a function of solid-liquid interface velocity and orientation. In simulations of Si with a second component k(sub neq) is greater in the [111] direction than the [100] direction in agreement with experimental results reported by Azlz et al. The simulated partition coefficient scales with the square of the step velocity divided by the diffusion coefficient of the secondary component in the liquid.

Beatty, K. M.↗

Nonlinear Dynamics and Nucleation Kinetics in Near-Critical Liquids

The objective of our study is to model the nonlinear behavior of a near-critical liquid following a rapid change of the temperature and/or other thermodynamic parameters (pressure, external electric or gravitational field). The thermodynamic critical point is manifested by large, strongly correlated fluctuations of the order parameter (particle density in liquid-gas systems, concentration in binary solutions) in the critical range of scales. The largest critical length scale is the correlation radius r(sub c). According to the scaling theory, r(sub c) increases as r(sub c) = r(sub 0)epsilon(exp -alpha) when the nondimensional distance epsilon = (T - T(sub c))/T(sub c) to the critical point decreases. The normal gravity alters the nature of correlated long-range fluctuations when one reaches epsilon approximately equal to 10(exp -5), and correspondingly the relaxation time, tau(r(sub c)), is approximately equal to 10(exp -3) seconds; this time is short when compared to the typical experimental time. Close to the critical point, a rapid, relatively small temperature change may perturb the thermodynamic equilibrium on many scales. The critical fluctuations have a hierarchical structure, and the relaxation involves many length and time scales. Above the critical point, in the one-phase region, we consider the relaxation of the liquid following a sudden temperature change that simultaneously violates the equilibrium on many scales. Below T(sub c), a non-equilibrium state may include a distribution of small scale phase droplets; we consider the relaxation of such a droplet following a temperature change that has made the phase of the matrix stable.

Patashinski, Alexander Z.↗

Emission Spectral Measurements in the Plenum of an Arc Jet Wind Tunnel

Arc jet wind tunnel facilities are used to evaluate thermal protection system materials for re-entry vehicles. The high speed, high temperature flowfield generated by the arc jet can simulate the extreme aerodynamic heating environment experienced during re-entry so that the survivability of heat shield materials and performance of various designs options can be tested. Although the re-entry heating environment can be approximated in the arc jet facility, the flowfield only partially simulates the actual re-entry flight conditions. Reynolds numbers are not matched so that surface shear stress distributions and mass transfer rates due to ablation or other mechanisms are not modeled correctly. Unlike flight conditions the arc freestream air is in non-equilibrium because of the rapid expansion that occurs in the supersonic nozzle. To properly study the actual re-entry flow environment, computational fluid dynamics, computational chemistry and radiation models must be used. Arc jet tunnel tests serve to validate these models. To perform accurate simulations inlet and boundary-conditions are needed, which come from measurements of the flowfield. The present study is concerned with measurements in the plenum region of an arc heater. In the past, conditions in the arc heater flowfield have been predicted using simulations since conventional measurement techniques could not be used in the harsh extremely high temperature environment. The present study is part of a recent push to utilize optical techniques to help better characterize the arc jet flowfields. Emission measurements have been made in the shock layer and the constrictor section of the arc heater to determine temperatures and species number densities. LIF measurements have been made in the free stream to determine temperature and velocity.

Donohue, Jim↗

High redshift quasars and high metallicities

A large-scale code called Cloudy was designed to simulate non-equilibrium plasmas and predict their spectra. The goal was to apply it to studies of galactic and extragalactic emission line objects in order to reliably deduce abundances and luminosities. Quasars are of particular interest because they are the most luminous objects in the universe and the highest redshift objects that can be observed spectroscopically, and their emission lines can reveal the composition of the interstellar medium (ISM) of the universe when it was well under a billion years old. The lines are produced by warm (approximately 10(sup 4)K) gas with moderate to low density (n less than or equal to 10(sup 12) cm(sup -3)). Cloudy has been extended to include approximately 10(sup 4) resonance lines from the 495 possible stages of ionization of the lightest 30 elements, an extension that required several steps. The charge transfer database was expanded to complete the needed reactions between hydrogen and the first four ions and fit all reactions with a common approximation. Radiative recombination rate coefficients were derived for recombination from all closed shells, where this process should dominate. Analytical fits to Opacity Project (OP) and other recent photoionization cross sections were produced. Finally, rescaled OP oscillator strengths were used to compile a complete set of data for 5971 resonance lines. The major discovery has been that high redshift quasars have very high metallicities and there is strong evidence that the quasar phenomenon is associated with the birth of massive elliptical galaxies.

Ferland, Gary J.↗

Dynamical Evolution of a Coronal Streamer-Flux Rope System: 2. A Self-Consistent Non-Planar Magnetohydrodynamic Simulation

The dynamical response of a helmet streamer to a flux rope escape from the sub-photosphere is examined in a physically self-consistent manner within the approximation of axisymmetric three-dimensional magnetohydrodynamics (i.e., so-called '2 1/2 D'). In contrast to the previous planar analyses of Paper 1 (Wu, Guo, and Wang), the present study shows, with the inclusion of out-of-plane components of magnetic and velocity fields, that the magnetic configuration represents a helical flux rope instead of a planar bubble as shown in Paper 1. Because of this more physically-realistic configuration, we are able to examine the dynamical evolution of the helical flux rope's interaction with the helmet streamer. This process leads to the formation of two parts of the solar mass ejection: (i) the expulsion of the helmet dome due to eruption of this flux rope, and (ii) the flux rope's eruption itself. When this two-part feature propagates out to the interplanetary space, it exhibits all the physical characteristics of observed interplanetary magnetic clouds. These numerical simulations also show that the dynamical behavior of the streamer-flux rope system has three distinct states: (i) quasi-equilibrium, (ii) non-equilibrium, and (iii) eruptive state depending on the energy level of the flux rope.

Wu, S. T.↗

Theory of Metastable State Relaxation in a Gravitational Field for Non-Critical Binary Systems with Non-Conserved Order Parameter

A new mathematical ansatz is developed for solution of the time-dependent Ginzburg-Landau nonlinear partial differential equation describing metastable state relaxation in binary (solute+solvent) non-critical solutions with non-conserved scalar order parameter in presence of a gravitational field. It has been demonstrated analytically that in such systems metastability initiates heterogeneous solute redistribution which results in the formation of a non-equilibrium singly-periodic spatial solute structure in the new solute-rich phase. The critical radius of nucleation and the induction time in these systems are gravity-dependent. It has also been proved that metastable state relaxation in vertical columns of supersaturated non-critical binary solutions leads to formation of the solute concentration gradient. Analytical expression for this concentration gradient is found and analysed. It is concluded that gravity can initiate phase separation (nucleation or spinodal decomposition).

Izmailov, Alexander F.↗

Theory of Metastable State Relaxation for Non-Critical Binary Systems with Non-Conserved Order Parameter

A new mathematical ansatz for a solution of the time-dependent Ginzburg-Landau non-linear partial differential equation is developed for non-critical systems such as non-critical binary solutions (solute + solvent) described by the non-conserved scalar order parameter. It is demonstrated that in such systems metastability initiates heterogeneous solute redistribution which results in formation of the non-equilibrium singly-periodic spatial solute structure. It is found how the time-dependent period of this structure evolves in time. In addition, the critical radius r(sub c) for solute embryo of the new solute rich phase together with the metastable state lifetime t(sub c) are determined analytically and analyzed.

Izmailov, Alexander↗

Analysis of Flow From Arc-Jet Spectra

Materials testing is carried out at the JSC Atmospheric Reentry Materials and Structures Facility. A flow diagnostics program is under development to characterize the energy distribution in arc-jet flows used to simulate atmospheric entry. Heat transfer to materials depends on flow properties which includes the composition of and energy distribution among the atoms, ions, molecules, and molecular ions which make up the flow. This project involves analysis of shock layer characteristics from the radiation emitted and experimentally collected from near the front of the shock to near the surface of the material. The analysis has yielded information on relative populations of neutral molecules and molecular ions within the layer. In determining non-equilibrium temperatures within the layer, some insight into the spectral constants used to compute radiative emission has been gained.

Blackwell, H. E.↗

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗

The NASA Ames Hypersonic Combustor-Model Inlet CFD Simulations and Experimental Comparisons

Computations have been performed on a three-dimensional inlet associated with the NASA Ames combustor model for the hypersonic propulsion experiment in the 16-inch shock tunnel. The 3-dimensional inlet was designed to have the combustor inlet flow nearly two-dimensional and of sufficient mass flow necessary for combustion. The 16-inch shock tunnel experiment is a short duration test with test time of the order of milliseconds. The flow through the inlet is in chemical non-equilibrium. Two test entries have been completed and limited experimental results for the inlet region of the combustor-model are available. A number of CFD simulations, with various levels of simplifications such as 2-D simulations, 3-D simulations with and without chemical reactions, simulations with and without turbulent conditions, etc., have been performed. These simulations have helped determine the model inlet flow characteristics and the important factors that affect the combustor inlet flow and the sensitivity of the flow field to these simplifications. In the proposed paper, CFD modeling of the hypersonic inlet, results from the simulations and comparison with available experimental results will be presented.

Venkatapathy, E.↗

On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow

In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Sweby, P. K.↗

Computer Simulation of Motion of Heliospheric Termination Shock

A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.

Ratkiewicz, R. E.↗

Molecular Infrared Intensity Standards

Molecular rotational and vibration-rotation spectroscopy can be a powerful tool for determining abundances and physical conditions in gases over a wide range of molecular abundance, temperature and to some extent pressure. In this paper I discuss experimental requirements for obtaining absolute intensity standards for stable molecules as well as radical molecular species. I also review the application of absolute and relative molecular intensity information for determining "conditions" in equilibrium and non-equilibrium gases.

Chackerian, C.↗

Comparison of Experimental and Computational Aerothermodynamics of a 70-deg Sphere-Cone

Numerical solutions for hypersonic flows of carbon-dioxide and air around a 70-deg sphere-cone have been computed using an axisymmetric non-equilibrium Navier-Stokes solver. Freestream flow conditions for these computations were equivalent to those obtained in an experimental blunt-body heat-transfer study conducted in a high-enthalpy, hypervelocity expansion tube. Comparisons have been made between the computed and measured surface heat-transfer rates on the forebody and afterbody of the sphere-cone and on the sting which supported the test model. Computed forebody heating rates were within the estimated experimental uncertainties of 10% on the forebody and 15% in the wake except for within the recirculating flow region of the wake.

Hollis, Brian R.↗

Hypervelocity Heat-Transfer Measurements in an Expansion Tube

A series of experiments has been conducted in the NASA HYPULSE Expansion Tube, in both CO2 and air test gases, in order to obtain data for comparison with computational results and to assess the capability for performing hypervelocity heat-transfer studies in this facility. Heat-transfer measurements were made in both test gases on 70 deg sphere-cone models and on hemisphere models of various radii. HYPULSE freestream flow conditions in these test gases were found to be repeatable to within 3-10%, and aerothermodynamic test times of 150 microsec in CO2 and 125 microsec in air were identified. Heat-transfer measurement uncertainty was estimated to be 10-15%. Comparisons were made with computational results from the non-equilibrium Navier-Stokes solver NEQ2D. Measured and computed heat-transfer rates agreed to within 10% on the hemispheres and on the sphere-cone forebodies, and to within 10% in CO2 and 25% in air on the afterbodies and stings of the sphere-cone models.

Hollis, Brian R.↗

Low-temperature carbonate concretions in the Martian meteorite ALH84001: evidence from stable isotopes and mineralogy

The martian meteorite ALH84001 contains small, disk-shaped concretions of carbonate with concentric chemical and mineralogical zonation. Oxygen isotope compositions of these concretions, measured by ion microprobe, range from delta18O = +9.5 to +20.5 per thousand. Most of the core of one concretion is homogeneous (16.7 +/- 1.2 per thousand) and over 5 per thousand higher in delta18O than a second concretion. Orthopyroxene that hosts the secondary carbonates is isotopically homogeneous (delta18O = 4.6 +/- 1.2 per thousand). Secondary SiO2 has delta18O = 20.4 per thousand. Carbon isotope ratios measured from the core of one concretion average delta13C = 46 +/- 8 per thousand, consistent with formation on Mars. The isotopic variations and mineral compositions offer no evidence for high temperature (>650 degrees C) carbonate precipitation and suggest non-equilibrium processes at low temperatures (< approximately 300 degrees C).

NASA Center JSC↗

A Technique for Determining the Nozzle-Flow Properties of Air in an Equilibrium, Nonequilibrium, or Frozen State

One of the problems associated with the development of high-enthalpy hypersonic test facilities hos been the determination of nozzle flow properties. This note presents a technique for determining the flow properties of air in an equilibrium, non-equilibrium, or frozen state using two test section measurements. A knowledge of two stagnation properties is required before the technique can be applied. The stagnation pressure can be measured, and reference 1 offers a method for determining stagnation enthalpy. Results of the calculations using the method of reference 1 are presented in reference 2 in a chart (chart 21) which can be readily used to determine the stagnation enthalpy from measurements of stagnation pressure, mass flow of air through the tunnel, and throat diameter.

Hayman, Lovick O.↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗