Search NASA⌕ Search

SEARCH · Search NASA

Results for “gas solution 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 181 records · Page 10

Extension of the Viscous Collision Limiting Direct Simulation Monte Carlo Technique to Multiple Species

There are many flows fields that span a wide range of length scales where regions of both rarefied and continuum flow exist and neither direct simulation Monte Carlo (DSMC) nor computational fluid dynamics (CFD) provide the appropriate solution everywhere. Recently, a new viscous collision limited (VCL) DSMC technique was proposed to incorporate effects of physical diffusion into collision limiter calculations to make the low Knudsen number regime normally limited to CFD more tractable for an all-particle technique. This original work had been derived for a single species gas. The current work extends the VCL-DSMC technique to gases with multiple species. Similar derivations were performed to equate numerical and physical transport coefficients. However, a more rigorous treatment of determining the mixture viscosity is applied. In the original work, consideration was given to internal energy non-equilibrium, and this is also extended in the current work to chemical non-equilibrium.

Liechty, Derek S.↗

Theory of Dust Voids in Plasmas

Dusty plasmas in a gas discharge often feature a stable void, i.e., a dust-free region inside the dust cloud. This occurs under conditions relevant to both plasma processing discharges and plasma crystal experiments. The void results from a balance of the electrostatic and ion drag forces on a dust particle. The ion drag force is driven by a flow of ions outward from an ionization source and toward the surrounding dust cloud, which has a negative space charge. In equilibrium the force balance for dust particles requires that the boundary with the dust cloud be sharp, provided that the particles are cold and monodispersive. Numerical solutions of the one-dimensional nonlinear fluid equations are carried out including dust charging and dust-neutral collisions, but not ion-neutral collisions. The regions of parameter space that allow stable void equilibria are identified. There is a minimum ionization rate that can sustain a void. Spatial profiles of plasma parameters in the void are reported. In the absence of ion-neutral collisions, the ion flow enters the dust cloud's edge at Mach number M = 1. Phase diagrams for expanding or contracting voids reveal a stationary point corresponding to a single stable equilibrium void size, provided the ionization rate is constant. Large voids contract and small voids expand until they attain this stationary void size. On the other hand, if the ionization rate is not constant, the void size can oscillate. Results are compared to recent laboratory and microgravity experiments.

Goree, J.↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

Hypersonic Aerodynamic Characteristics of a Proposed Single-Stage-To-Orbit Vehicle

The hypersonic aerodynamic characteristics of a winged body concept representing a candidate single- stage-to-orbit vehicle which features wing tip fin controllers and elevon/body flap control surfa'Fs are predicted at points along a nominal trajectory for Mach numbers from 5 to 27 and angles of attack from 19 to 32 degrees. Predictions are derived from surface properties based on flow solvers for inviscid and viscous, laminar flows acting as a perfect gas, as a gas in chemical equilibrium and as a gas in chemical non- equilibrium. At a Mach number of 22, the lateral aerodynamic characteristics of the vehicle are determined based on an inviscid analysis at side slip angles of 2 and 4 degrees and 32 degrees angle of attack; a viscous analysis was carried out to determine the effect of gas chemistry model on surface pressure and to determine the incremental aerodynamics for control surface deflections. The results show that the longitudinal pitch characteristics of the baseline configuration, i.e., zero control surface deflections, are significantly altered by real gas chemistry at angles of attack greater than 30 degrees and Mach numbers greater than 9; and, that aerodynamics derived from inviscid solutions are of sufficient accuracy for preliminary analysis. Also, it is shown that a Mach number of 22, the choice of gas chemistry model has a large impact on surface pressure levels at highly localized regions on the vehicle and that the vehicle can be trimmed at control surface deflections less than 11 degrees.

Weilmuenster, K. James↗

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗

A flux-limited treatment for the conductive evaporation of spherical interstellar gas clouds

In this work, we present and analyze a new analytic solution for the saturated (flux-limited) thermal evaporation of a spherical cloud. This work is distinguished from earlier analytic studies by allowing the thermal conductivity to change continuously from a diffusive to a saturated form, in a manner usually employed only in numerical calculations. This closed form solution will be of interest as a computational benchmark. Using our calculated temperature profiles and mass-loss rates, we model the thermal evaporation of such a cloud under typical interstellar medium (ISM) conditions, with some restrictions. We examine the ionization structure of the cloud-ISM interface and evaluate column densities of carbon, nitrogen, oxygen, neon, and silicon ions toward the cloud. In accord with other investigations, we find that ionization equilibrium is far from satisfied under the assumed conditions. Since the inclusion of saturation effects in the heat flux narrows the thermal interface relative to its classical structure, we also find that saturation effects tend to lower predicted column densities.

Dalton, William W.↗

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗

Tables and charts of equilibrium normal shock and shock-tube solutions for helium-hydrogen mixtures with velocities to 70 km/sec

Equilibrium thermodynamic and flow properties are presented in tabulated and graphical form for moving, standing, and reflected normal shock waves into helium-hydrogen mixtures representative of proposed outer planet atmospheres. The volumetric compositions of these mixtures are 0.35He-0.65H2, 0.20He-0.80H2, and 0.05He-0.95H2. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular-weight ratio, isentropic exponent, velocity, and species mole fractions. Incident (moving) shock velocities are varied from 4 to 70 km/sec for a range of initial pressure of 5 N/sq m to 100 kN/sq m. The present results are applicable to shock-tube flows and to free-flight conditions for a blunt body at high velocities. A working chart illustrating idealized shock-tube performance with a 0.20He-0.80H2 test gas and heated helium driver gas is also presented.

Miller, C. G., III↗

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.↗

Evaporation and combustion of LOX under supercritical and subcritical conditions

The objective is to study the evaporation and combustion of LOX under supercritical and subcritical conditions both experimentally and theoretically. In the evaporation studies, evaporation rate and surface temperature were measured when LOX vaporizing in helium environments at pressures ranging from 5 to 68 atm. A Varian 3700 gas chromatograph was employed to measure the oxygen concentration above the LOX surface. For the combustion tests, high-magnification video photography was used to record direct images of the flame shape of a LOX/H2/He laminar diffusion flame. The gas composition in the post-flame region is also being measured with the gas sampling and chromatography analysis. These data are being used to validate the theoretical model. A comprehensive theoretical model with the consideration of the solubility of ambient gases as well as variable thermophysical properties was formulated and solved numerically to study the gasification and burning of LOX at elevated pressures. The calculated flame shape agreed reasonably well with the edge of the observed luminous flame surface. The effect of gravity on the flame structure of laminar diffusion flames was found to be significant. In addition, the predicted results using the flame-sheet model were compared with those based upon full equilibrium calculations (which considered the formation of intermediate species) at supercritical pressures. Except at the flame front where temperature exceeded 2,800 K, the flame-sheet and equilibrium solutions in terms of temperature distributions were in very close agreement. The temperature deviation in the neighborhood of the flame front is caused by the effect of high-temperature dissociation.

Yang, A. S.↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 1: Theory and numerical solution procedures

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 1 of a series of three reference publications that describe LENS, provide a detailed guide to its usage, and present many example problems. Part 1 derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved. The accuracy and efficiency of LSENS are examined by means of various test problems, and comparisons with other methods and codes are presented. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

A theoretical and observational analysis of He I triplet lines in class 1 Seyfert galaxy spectra

The equations of statistical equilibrium are solved for 11 levels of the helium atom that are most relevant to the production and radiative transfer of He I triplet lines. Solutions are obtained for several densities and temperatures over a run of optical depth in the 10,830-A line, and the results for emission-line strengths are plotted. Spectrophotometric scans and pertinent line strengths are presented for the class 1 Seyfert galaxies 3C 120, NGC 7469, and Mrk 618, as well as for Mrk 335, which exhibits a 'more normal' He I 5876/H-beta ratio. A comparison of the theoretically determined line intensities with the observational results verifies that the He I triplet line intensities are significantly enhanced by electron collisional excitation from 2(3)P and 2(3)S. It is suggested that the gas which gives rise to the broad He I lines in 3C 120, Mrk 618, and NGC 7469 is characterized by a number density of about 5 billion per cu cm, a temperature of approximately 15,000 K, and an optical depth of 100 to 120 in the 10,830-A line.

Feldman, F. R.↗

Boundary Layer Transition Correlations and Aeroheating Predictions for Mars Smart Lander

Laminar and turbulent perfect-gas air, Navier-Stokes computations have been performed for a proposed Mars Smart Lander entry vehicle at Mach 6 over a free stream Reynolds number range of 6.9 x 10(exp 6/m to 2.4 x 10(exp 7)m(2.1 x 10(exp 6)/ft to 7.3 x 10(exp 6)ft) for angles-of-attack of 0-deg, 11-deg, 16-deg, and 20-deg, and comparisons were made to wind tunnel heating data obtained at the same conditions. Boundary layer edge properties were extracted from the solutions and used to correlate experimental data on the effects of heat-shield penetrations (bolt-holes where the entry vehicle would be attached to the propulsion module during transit to Mars) on boundary-layer transition. A non-equilibrium Martian-atmosphere computation was performed for the peak heating point on the entry trajectory in order to determine if the penetrations would produce boundary-layer transition by using this correlation.

Hollis, Brian R.↗

Boundary Layer Transition Correlations and Aeroheating Predictions for Mars Smart Lander

Laminar and turbulent perfect-gas air, Navier-Stokes computations have been performed for a proposed Mars Smart Lander entry vehicle at Mach 6 over a free stream Reynolds number range of 6.9 x 10(exp 6)/m to 2.4 x 10(exp 7)/m (2.1 x 10(exp 6)/ft to 7.3 x 10(exp 6)/ft) for angles-of-attack of 0-deg, 11-deg, 16-deg, and 20-deg, and comparisons were made to wind tunnel heating data obtained a t the same conditions. Boundary layer edge properties were extracted from the solutions and used to correlate experimental data on the effects of heat-shield penetrations (bolt-holes where the entry vehicle would be attached to the propulsion module during transit to Mars) on boundary-layer transition. A non-equilibrium Martian-atmosphere computation was performed for the peak heating point on the entry trajectory in order to determine if the penetrations would produce boundary-layer transition by using this correlation.

Hollis, Brian R.↗

Equilibrium slab models of Lyman-alpha clouds

Solutions for the equilibrium configuration of a slab with ionizing radiation incident equally from both sides are explored. Radiation effects (photoionization, Ly-alpha photon trapping, and mock gravity) as well as external pressure and self gravity (with and without dark matter) are included. The general formalism is applied to structure growth on small scales at very high z due to mock gravity on dust. Emphasis is placed on the application of slab models at z of less than 5, particularly those that may correspond to Ly-alpha forest, Lyman limit, and damped Ly-alpha systems. The regime with a dominant outward force contributed by trapping of Ly-alpha photons is discussed. General expressions are given for the equilibrium, including dark matter, assuming various relationships between the density of the dark matter halo and the total gas column density.

Charlton, Jane C.↗

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗

Swept shock/boundary layer interaction experiments in support of CFD code validation

Research on the topic of shock wave/turbulent boundary-layer interaction was carried out during the past three years at the Penn State Gas Dynamics Laboratory. This report describes the experimental research program which provides basic knowledge and establishes new data on heat transfer in swept shock wave/boundary-layer interactions. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angle ranging from 10 deg to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a k-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.

Settles, G. S.↗