Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution species”

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 271 records · Page 15

Electromagnetic studies of redox systems for energy storage

Both chromium and iron couples were studied on various electrode surfaces in acidic perchlorate solution by using rotating ring-disk techniques. It was found that chloride which forms inner sphere coordination complexes with the redox species enhances the electrode kinetics dramatically. The effects of lead underpotential deposition and surface alloy formation on the kinetics of the chromium couple on gold were studied using both linear sweep voltammetry and potential step techniques. The lad underpotential deposition was found to slow down the kinetics of the reduction of the Cr species on gold surfaces although increase the hydrogen overvoltage. The effect on the chromium kinetics can be explained in terms of principally a double layer effect. The underpotential deposition lead species with its positive charge results in a decrease in the concentration of the Cr species at the electrode surface. Similar phenomena were also observed with bismuth underpotential deposition on gold for the iron couple.

Wu, C. D.↗

Stability of a laminar premixed supersonic free shear layer with chemical reactions

The stability of a two-dimensional compressible supersonic flow in the wake of a flat plate is discussed. The fluid is a multi-species mixture which is undergoing finite rate chemical reactions. The spatial stability of an infinitesimal disturbance in the fluid is considered. Numerical solutions of the eigenvalue stability equations for both reactive and nonreactive supersonic flows are presented and discussed. The chemical reactions have significant influence on the stability behavior. For instance, a neutral eigenvalue is observed near the freestream Mach number of 2.375 for the nonreactive case, but disappears when the reaction is turned on. For reactive flows, the eigenvalues are not very dependent on the free stream Mach number.

Menon, S.↗

Nonlinear upper hybrid drift waves for a longitudinal electric field perpendicular to a uniform magnetic field in the Vlasov-Maxwell approximation

Upper hybrid drift waves are found as a special solution to a Vlasov-Maxwell plasma which has a longitudinal electric field and a perpendicular uniform magnetic field. A single-species plasma with a constant-density mobile neutralizing background supports spatially varying disturbances that oscillate at the upper hybrid frequency. The general functional dependences of the electric field, the plasma number density, and the one-particle distribution function for the special case are found from more general Vlasov-Maxwell equations invariant under a Lie group point transformation. The one-particle distribution function for the plasma is a function of the Liouville invariant, which is the energy in the generalized Bernstein-Greene-Kruskal (BGK) reference frame, and the momentum in the drift direction.

Abraham-Shrauner, B.↗

Enhanced thermochemical nonequilibrium computations of flow around the aeroassist flight experiment vehicle

A three-dimensional explicit, finite-rate, shock-capturing numerical algorithm is used to calculate thermochemical nonequilibrium flowfields about the Aeroassist Flight Experiment vehicle at one of its flight trajectory points. The full Navier-Stokes equations and an eleven species chemical model with the latest reaction rates are incorporated into the code. Results are compared against experimental data and other numerical solutions. The effects of changes to the physical model on the stagnation line and base region flow and surface quantities on the forebody are investigated.

Palmer, Grant↗

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.↗

HANA - A three-dimensional Navier-Stokes code for chemically reacting flows

A recently developed computer program for analyzing chemically reacting flows is described. The program uses a diagonal implicit algorithm to solve the equations of viscous flow for a gas in thermochemical nonequilibrium. The diagonal scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The program uses multiple zones of grids patched together and includes a k-epsilon turbulence model. Solutions are presented for hypersonic flows of air, hypersonic flows of Martian atmospheric gases, combusting flows of hydrogen/air mixtures, and combusting flows of methane/air mixtures.

Imlay, Scott T.↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗

Uniqueness of a solution of a steady state photochemical problem: Applications to Mars

Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or nonzero velocities should not be less than the number of chemical elements in the system, and the boundary conditions for species given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is shown that the problem of the CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have a unique solution, if only CO2 and H2O densities are given at the lower boundary, and the remaining boundary conditions are fluxes. Two examples of models of this type are discussed. Two models of the photochemistry of the Martian atmosphere, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape ratio of 1.2 x 10(exp 8)/cu cm/s is given at 240 km and is balanced with the total hydrogen escape rate within an uncertainty of 1% for both models. Both models fit the measured O2 and CO mixing ratios, the O3 abundance, and the O2 1.27-micrometer dayglow almost within the uncertainties of the measured values, though the model without nitrogen chemistry fits better. The importance of nitrogen chemistry in the lower and middle atmosphere of Mars depends on a fine balance between production of NO and N in the upper atmosphere which is not known within the required accuracy.

Krasnopolsky, Vladimir A.↗

Fluorescence Studies of Protein Crystal Nucleation

One of the most powerful and versatile methods for studying molecules in solution is fluorescence. Crystallization typically takes place in a concentrated solution environment, whereas fluorescence typically has an upper concentration limit of approximately 1 x 10(exp -5)M, thus intrinsic fluorescence cannot be employed, but a fluorescent probe must be added to a sub population of the molecules. However the fluorescent species cannot interfere with the self-assembly process. This can be achieved with macromolecules, where fluorescent probes can be covalently attached to a sub population of molecules that are subsequently used to track the system as a whole. We are using fluorescence resonance energy transfer (FRET) to study the initial solution phase self-assembly process of tetragonal lysozyme crystal nucleation, using covalent fluorescent derivatives which crystallize in the characteristic P432121 space group. FRET studies are being carried out between cascade blue (CB-lys, donor, Ex 376 nm, Em 420 nm) and lucifer yellow (LY-lys, acceptor, Ex 425 nm, Em 520 nm) asp101 derivatives. The estimated R0 for this probe pair, the distance where 50% of the donor energy is transferred to the acceptor, is approximately 1.2 nm, compared to 2.2 nm between the side chain carboxyls of adjacent asp101's in the crystalline 43 helix. The short CB-lys lifetime (approximately 5 ns), coupled with the large average distances between the molecules ((sup 3) 50 nm) in solution, ensure that any energy transfer observed is not due to random diffusive interactions. Addition of LY-lys to CB-lys results in the appearance of a second, shorter lifetime (approximately 0.2 ns). Results from these and other ongoing studies will be discussed in conjunction with a model for how tetragonal lysozyme crystals nucleate and grow, and the relevance of that model to microgravity protein crystal growth

Pusey, Marc L.↗

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

An improved algorithm for the modeling of vapor flow in heat pipes

A heat pipe vapor flow algorithm suitable for use in codes on microcomputers is presented. The incompressible heat pipe vapor flow studies of Busse are extended to incorporate compressibility effects. The Busse velocity profile factor is treated as a function of temperature and pressure. The assumption of a uniform saturated vapor temperature determined by the local pressure at each cross section of the pipe is not made. Instead, a mean vapor temperature, defined by an energy integral, is determined in the course of the solution in addition to the pressure, saturation temperature at the wall, and the Busse velocity profile factor. For alkali metal working fluids, local species equilibrium is assumed. Temperature and pressure profiles are presented for several cases involving sodium heat pipes. An example for a heat pipe with an adiabatic section and two evaporators in sequence illustrates the ability to handle axially varying heat input. A sonic limit plot for a short evaporator falls between curves for the Busse and Levy inviscid sonic limits.

Tower, Leonard K.↗

Nonequilibrium thermo-chemical calculations using a diagonal implicit scheme

A recently developed computer program for hypersonic vehicle flow analysis is described. The program uses a diagonal implicit algorithm to solve the equations of viscous flow for a gas in thermochemical nonequilibrium. The diagonal scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The program uses multiple zones of grids patched together and includes radiation wall and rarefied gas boundary conditions. Solutions are presented for hypersonic flows of air and hydrogen air mixtures.

Imlay, Scott T.↗

UV-Resistant Non-Spore-Forming Bacteria From Spacecraft-Assembly Facilities

Four species of non-spore-forming bacteria collected from clean-room surfaces in spacecraft-assembly facilities could survive doses of ultraviolet (UV) radiation that would suffice to kill most known cultivable bacterial species. In a previous study, high UV resistance was found in spores of the SAFR-032 strain of Bacillus pumilus, as reported in "Ultraviolet- Resistant Bacterial Spores," NASA Tech Briefs, Vol. 31, No. 9 (September 2007), page 94. These studies are parts of a continuing effort to understand the survival of hardy species of bacteria under harsh conditions, and develop means of sterilizing spacecraft to prevent biocontamination of Mars that could in turn interfere with future life detection missions. The four species investigated were Arthrobacter sp. KSC_Ak2i, Microbacterium schleiferi LMA_AkK1, Brevundimonas diminuta KSC_Ak3a, and Sphingomonas trueperi JSC_Ak7-3. In the study, cells of these species were mixed into Atacama Desert soil (to elucidate the shadowing effect of soil particles) and the resulting mixtures were tested both in solution and in a desiccated state under simulated Martian atmospheric and UV conditions. The UV-survival indices of Arthrobacter sp. and Microbacterium schleiferi were found to be comparable to those of Bacillus pumilus spores.

Venkateswaran, Kasthuri↗

Evidence for Ammonia-bearing Species on the Uranian Satellite Ariel Supports Recent Geologic Activity

We investigated whether ammonia-rich constituents are present on the surface of the Uranian moon Ariel by analyzing 32 near-infrared reflectance spectra collected over a wide range of sub-observer longitudes and latitudes. We measured the band areas and depths of a 2.2 μm feature in these spectra, which has been attributed to ammonia-bearing species on other icy bodies. Ten spectra display prominent 2.2 μm features with band areas and depths >2σ. We determined the longitudinal distribution of the 2.2 μm band, finding no statistically meaningful differences between Ariel’s leading and trailing hemispheres, indicating that this band is distributed across Ariel’s surface. We compared the band centers and shapes of the five Ariel spectra displaying the strongest 2.2 μm bands to laboratory spectra of various ammonia-bearing and ammonium-bearing species, finding that the spectral signatures of the Ariel spectra are best matched by ammonia-hydrates and flash frozen ammonia-water solutions. Our analysis also revealed that four Ariel spectra display 2.24 μm bands (>2σ band areas and depths), with band centers and shapes that are best matched by ammonia ice. Because ammonia should be efficiently removed over short timescales by ultraviolet photons, cosmic rays, and charged particles trapped in Uranus’ magnetosphere, the possible presence of this constituent supports geologic activity in the recent past, such as emplacement of ammonia-rich cryolavas and exposure of ammonia-rich deposits by tectonism, impact events, and mass wasting.

Richard J. Cartwright↗

Monte Carlo simulation of the classical two-dimensional one component plasma

Monte Carlo simulation, lattice dynamics in the harmonic approximation, and solution of the hypernetted chain equation were used to study the classical two-dimensional one component plasma. The system consists of a single species of charged particles immersed in a uniform neutralizing background. The particles interact via a l/r potential, where r is the two dimensional separation. Equations of state were calculated for both the liquid and solid phases. Results of calculation of the thermodynamic functions and one and two particle correlation functions are presented.

Gann, R. C.↗

Probing the presently tenuous link between comets and the origin of life

The possibilities of using millimeter-wave technology to probe the subsurface processes of comets to investigate links between cometary materials and the origins of life are explored. It is noted that current theories hold that the necessities for life to begin comprise a fairly uniform temperature, the presence of a solvent to give materials mobility, and the presence of atoms which can form long chains of molecules. Consideration is given to two cometary nuclei models: a core with an equal amount of liquid water and lunar material, and a nucleus with equal amounts of frozen water ice and lunar material. Solutions to the radiative transfer equation for the two models are presented to characterize identifiable emissions using radiometric spectrometer instrumentation on a spacecraft. Particular species such as OH, CN, HCN, and glycine are expected to be detectable if present.

Hobbs, R. W.↗

A Priori Analysis of Flamelet-Based Modeling for a Dual-Mode Scramjet Combustor

An a priori investigation of the applicability of flamelet-based combustion models to dual-mode scramjet combustion was performed utilizing Reynolds-averaged simulations (RAS). For this purpose, the HIFiRE Direct Connect Rig (HDCR) flowpath, fueled with a JP-7 fuel surrogate and operating in dual- and scram-mode was considered. The chemistry of the JP-7 fuel surrogate was modeled using a 22 species, 18-step chemical reaction mechanism. Simulation results were compared to experimentally-obtained, time-averaged, wall pressure measurements to validate the RAS solutions. The analysis of the dual-mode operation of this flowpath showed regions of predominately non-premixed, high-Damkohler number, combustion. Regions of premixed combustion were also present but associated with only a small fraction of the total heat-release in the flow. This is in contrast to the scram-mode operation, where a comparable amount of heat is released from non-premixed and premixed combustion modes. Representative flamelet boundary conditions were estimated by analyzing probability density functions for temperature and pressure for pure fuel and oxidizer conditions. The results of the present study reveal the potential for a flamelet model to accurately model the combustion processes in the HDCR and likely other high-speed flowpaths of engineering interest.

Quinlan, Jesse R.↗

BLAYER User Guide

Estimation of boundary layer edge or thickness is an essential aspect of utilizing high fidelity aerothermodynamic flow solutions. For viscous flow over a flat plate, the boundary layer thickness is usually defined as the distance from the wall at which the flow velocity (zero at the wall) has increased through some profile to 99% or 99.5% of the free stream velocity. In real gas flows over aerospace vehicles, a total enthalpy profile is employed for edge calculations. The importance of knowing boundary layer thickness and its relation to possible transition from laminar to turbulent flow is illustrated by Space Shuttle mission STS-114, when the appearance of tile gap filler protrusions on the nose was deemed serious enough that a crew member was called upon to exit the vehicle and remove the potential threat prior to its descent from orbit. The BLAYER utility used as part of post-processing Orion/MPCV flow solutions is a generalization of the earlier BLAYER_RESULTS utility developed following the STS-107 Columbia accident. BLAYER handles a variable number of gas species (not just 5), allows for more than one temperature and for optional extra flow quantities, and omits optional handling of Orbiter tile datasets.

momentum thicknes↗