SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Databases of transitional states obtained from Direct Numerical simulations (DNS) of temporal, supercritical mixing layers for two species systems, O2/H2 and C7H16/N2, are analyzed to elucidate species-specific turbulence aspects and features of fluid disintegration.
Databases of transitional states obtained from Direct Numerical Simulations (DNS) of temporal, supercritical mixing layers for two species systems, 02/H2 and C7Hle/N2, are analyzed to elucidate species-specific turbulence aspects.
The hydrogen- and self-broadened line widths of 116 (sup 14)NH(sub 3) ground state transitions with J,K = 1,0 to 10,10 have been measured at 0.0006 cm(sup -1) resolution using a Bruker spectrometer between 40 to 210 cm(sup -1). These experimental widths have been reproduced to 2.4% and 11% respectively using an heuristically derived expression of the form....
Single planar Langmuir probes and fiber optic probes are used to concurrently measure the plasma properties and neutral density variation in a 30cm diameter ion engine discharge chamber, from the immediate vicinity of the keeper to the near grid plasma region. The fiber optic probe consists of a collimated optical fiber recessed into a double bore ceramic tube fitted with a stainless steel light-limiting window. The optical fiber probe is used to measure the emission intensity of excited neutral xenon for a small volume of plasma, at various radial and axial locations. The single Langmuir probes, are used to generate current-voltage characteristics at a total of 140 spatial locations inside the discharge chamber. Assuming a maxwellian distribution for the electron population, the Langmuir probe traces provide spatially resolved measurements of plasma potential, electron temperature, and plasma density. Data reduction for the NSTAR TH8 and TH15 throttle points indicates an electron temperature range of 1 to 7.9 eV and an electron density range of 4e10 to le13 cm(sup -3), throughout the discharge chamber, consistent with the results in the literature. Plasma potential estimates, computed from the first derivative of the probe characteristic, indicate potential from 0.5V to 11V above the discharge voltage along the thruster centerline. These values are believed to be excessively high due to the sampling of the primary electron population along the thruster centerline. Relative neutral density profiles are also obtained with a fiber optic probe sampling photon flux from the 823.1 nm excited to ground state transition. Plasma parameter measurements and neutral density profiles will be presented as a function of probe location and engine discharge conditions. A discussion of the measured electron energy distribution function will also be presented, with regards to variation from pure maxwellian. It has been found that there is a distinct primary population found along the thruster centerline, which causes estimates of electron temperature, electron density, and plasma potential, to err on the high side, due this energetic population. Computation of the energy distribution fimction of the plasma clearly indicates the presence of primaries, whose presence become less obvious with radial distance from the main discharge plume.
Models for large eddy simulation (LES) are assessed on a database obtained from direct numerical simulations (DNS) of supercritical binary-species temporal mixing layers. The analysis is performed at the DNS transitional states for heptane/nitrogen, oxygen/hydrogen and oxygen/helium mixing layers. The incorporation of simplifying assumptions that are validated on the DNS database leads to a set of LES equations that requires only models for the subgrid scale (SGS) fluxes, which arise from filtering the convective terms in the DNS equations. Constant-coefficient versions of three different models for the SGS fluxes are assessed and calibrated. The Smagorinsky SGS-flux model shows poor correlations with the SGS fluxes, while the Gradient and Similarity models have high correlations, as well as good quantitative agreement with the SGS fluxes when the calibrated coefficients are used.
Density functional theory calculations at the B3LYP/6-31+G** level were employed to characterize the critical points for adducts, isomers, products, and intervening transition states for the reactions between benzene and the ethynyl (C2H) or cyano (CN) radicals. Both addition reactions were found to have no barriers in their entrance channels, making them efficient at the low temperature and pressure conditions that prevail in the haze-forming region of Titan's atmosphere as well as in the dense interstellar medium (ISM). The dominant products are ethynylbenzene (C6H5C2H) and cyanobenzene (C6H5CN). Hydrogen abstraction reactions were also characterized but found to be non-competitive. Trajectory calculations based on potentials fit to about 600 points calculated at the ROMP2/6-31+G** level for each interaction surface were used to determine reaction rates. The rates incorporated any necessary corrections for back reactions as ascertained from a multiwell treatment used to determine outcome distributions over the range of temperatures and pressures pertinent to Titan and the ISM and are in good agreement with the limited available experimental data.
A computational-simulation study has been presented of effects of perturbation wavelengths and initial Reynolds numbers on the transition to turbulence of a heptane/nitrogen mixing layer at supercritical pressure. The governing equations for the simulations were the same as those of related prior studies reported in NASA Tech Briefs. Two-dimensional (2D) simulations were performed with initially im posed span wise perturbations whereas three-dimensional (3D) simulations had both streamwise and spanwise initial perturbations. The 2D simulations were undertaken to ascertain whether perturbations having the shortest unstable wavelength obtained from a linear stability analysis for inviscid flow are unstable in viscous nonlinear flows. The goal of the 3D simulations was to ascertain whether perturbing the mixing layer at different wavelengths affects the transition to turbulence. It was found that transitions to turbulence can be obtained at different perturbation wavelengths, provided that they are longer than the shortest unstable wavelength as determined by 2D linear stability analysis for the inviscid case and that the initial Reynolds number is proportionally increased as the wavelength is decreased. The transitional states thus obtained display different dynamic and mixture characteristics, departing strongly from the behaviors of perfect gases and ideal mixtures.
Product channels for the self-reaction of the resonance-stabilized allyl radical, C3H5 + C3H5, have been studied with isomeric specificity at temperatures from 300-600 K and pressures from 1-6 Torr using time-resolved multiplexed photoionization mass spectrometry. Under these conditions 1,5-hexadiene was the only C6H10 product isomer detected. The lack of isomerization of the C6H10 product is in marked contrast to the C6H6 product in the related C3H3 + C3H3 reaction, and is due to the more saturated electronic structure of the C6H10 system. The disproportionation product channel, yielding allene + propene, was also detected, with an upper limit on the branching fraction relative to recombination of 0.03. Analysis of the allyl radical decay at 298 K yielded a total rate coefficient of (2.7 +/- 0.8) x 10(exp -11) cu cm/molecule/s, in good agreement with pre.vious experimental measurements using ultraviolet kinetic absorption spectroscopy and a recent theoretical determination using variable reaction coordinate transition state theory. This result provides independent indirect support for the literature value of the allyl radical ultraviolet absorption cross-section near 223 nm.
A particular periodic orbit in the Earth-Sun circular restricted three body problem is shown to have the characteristics needed for a ballistic lunar capture transfer. An injection from a circular parking orbit into the periodic orbit serves as an initial guess for a targeting algorithm. By targeting appropriate parameters incrementally in increasingly complicated force models and using precise derivatives calculated from the state transition matrix, a reliable algorithm is produced. Ballistic lunar capture trajectories in restricted four body systems are shown to be able to be produced in a systematic way.
This paper discusses a method used for the systematic improvement of NASA s Lunar Surface Systems avionics architectures in the area of reliability and fault-tolerance. This approach utilizes an integrated system model to determine the effects of component failure on the system s ability to provide critical functions. A Markov model of the potential degraded system modes is created to characterize the probability of these degraded modes, and the system model is run for each Markov state to determine its status (operational or system loss). The probabilistic results from the Markov model are first produced from state transition rates based on NASA data for heritage failure rate data of similar components. An additional set of probabilistic results are created from a representative set of failure rates developed for this study, for a variety of component quality grades (space-rated, mil-spec, ruggedized, and commercial). The results show that careful application of redundancy and selected component improvement should result in Lunar Surface Systems architectures that exhibit an appropriate degree of fault-tolerance, reliability, performance, and affordability.
Motivated by large eddy simulation (LES) modeling of supercritical turbulent flows, transitional states of databases obtained from direct numerical simulations (DNS) of binary-species supercritical temporal mixing layers were examined to understand the subgrid-scale dissipation, and its variation with filter size. Examination of the DSN-scale domain- averaged dissipation confirms previous findings that, out of the three modes of viscous, temperature and species-mass dissipation, the species-mass dissipation is the main contributor to the total dissipation. The results revealed that the percentage of species-mass by total dissipation is nearly invariant across species systems and initial conditions. This dominance of the species-mass dissipation is due to high-density-gradient magnitude (HDGM) regions populating the flow under the supercritical conditions of the simulations; such regions have also been observed in fully turbulent supercritical flows. The domain average being the result of both the local values and the extent of the HDGM regions, the expectations were that the response to filtering would vary with these flow characteristics. All filtering here is performed in the dissipation range of the Kolmogorov spectrum, at filter sizes from 4 to 16 times the DNS grid spacing. The small-scale (subgrid scale, SGS) dissipation was found by subtracting the filtered-field dissipation from the DNS-field dissipation. In contrast to the DNS dissipation, the SGS dissipation is not necessarily positive; negative values indicate backscatter. Backscatter was shown to be spatially widespread in all modes of dissipation and in the total dissipation (25 to 60 percent of the domain). The maximum magnitude of the negative subgrid- scale dissipation was as much as 17 percent of the maximum positive subgrid- scale dissipation, indicating that, not only is backscatter spatially widespread in these flows, but it is considerable in magnitude and cannot be ignored for the purposes of LES modeling. The Smagorinsky model, for example, is unsuited for modeling SGS fluxes in the LES because it cannot render backscatter. With increased filter size, there is only a modest decrease in the spatial extent of backscatter. The implication is that even at large LES grid spacing, the issue of backscatter and related SGS-flux modeling decisions are unavoidable. As a fraction of the total dissipation, the small-scale dissipation is between 10 and 30 percent of the total dissipation for a filter size that is four times the DNS grid spacing, with all OH cases bunched at 10 percent, and the HN cases spanning 24 30 percent. A scale similarity was found in that the domain-average proportion of each small-scale dissipation mode, with respect to the total small-scale dissipation, is very similar to equivalent results at the DNS scale. With increasing filter size, the proportion of the small-scale dissipation in the dissipation increases substantially, although not quite proportionally. When the filter size increases by four-fold, 52 percent for all OH runs, and 70 percent for HN runs, of the dissipation is contained in the subgrid-scale portion with virtually no dependence on the initial conditions of the DNS. The indications from the dissipation analysis are that modeling efforts in LES of thermodynamically supercritical flows should be focused primarily on mass-flux effects, with temperature and viscous effects being secondary. The analysis also reveals a physical justification for scale-similarity type models, although the suitability of these will need to be confirmed in a posteriori studies.
The environmental signatures for habitability are not necessarily biosignatures, even though on Earth, they are definitive proof of habitability. It is the constant overprint of the chemical signatures of life that makes it difficult to recognize the chemical and physical properties of a potentially habitable environment as distinct from an inhabited one. Mars Science Laboratory (MSL) will soon embark on a mission to Mars to assess its past or present habitability, so it is useful to examine how we measure habitability on Earth and prepare for how that approach may differ for Mars. This exercise includes: (a) articulation of fundamental assumptions about habitability, (b) an inventory of factors that affect habitability, (c) development of metrics, measurement approach and implementation, and (d) a new classification scheme for planetary habitability that goes beyond the binary "yes" or "no." There may be dozens of factors that affect habitability and they can be weighted as a function of specific environment. However a robotic, in situ investigation even on Earth has constraints that prevent the measurement of every environmental factor, so metrics must be reduced to the most relevant subset, given available time, cost, technical feasibility and scientific importance. Many of the factors could be measured with a combination of orbital data and the MSL payload. We propose that, at a minimum, a designation of high habitability potential requires the following conditions be met: (a) thermally stable with respect to extremes and frequency of fluctuation, (b) has more than one energy source, (c) sufficient chemical diversity to make compounds with covalent and hydrogen bonding, (d) can moderate ionizing radiation enough to allow a stable or evolving pool of organic molecules, (e) must have water or other high quality polar solvent, (f) must be able to renew chemical resources (e.g., plate tectonics, volcanism or something else we haven't envisioned). A measurement approach we have taken to measure habitability on Earth is : 1. Study remote sensing data, maps, etc. 2. Decide how big an area to measure. 3. Determine the spatial sampling rate. 4. Determine the temporal sampling rate. 5. Determine the order of measurements 6. Decide where to begin measurements 7. Select locations at field site and proceed While science drives each of the steps, there are additional constraints, e.g., technical, time, cost, safety (risk). This approach is also executable on Mars. Measurement of past habitability is more challenging both for Earth and Mars where access to the past means subsurface access and confrontation with unknowns about preservation of the martian past. Some environments preserve evidence of past habitability better than others, and this is where selection of the landing site to maximize the preservation potential of habitability indicators will be key. Mars presents an opportunity to discover transitional states between habitable or not, and we offer a ranking scale for planetary habitability with Mars as the second test subject: CLASS ONE Uninhabitable and likely has never been so CLASS TWO Has a high potential but no confirmed observation of life (as defined above) CLASS THREE Inhabited (we find life) 3-A Globally inhabited 3-B Primitive life; early in its evolution, but not yet globally established 3-C Exists only in refugia -- planet heading toward class four CLASS FOUR Post-habitable (there once was life, but now it's gone) MSL provides an opportunity to carefully investigate the habitability of at least one site on Mars and it will reveal much about the possible states of planetary habitability
Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.
This report presents a study of dissipation (irreversible production of entropy) in three-dimensional, temporal mixing layers laden with evaporating liquid drops. The purpose of the study is to examine the effects of evaporating drops on the development of turbulent features in flows. Direct numerical simulations were performed to analyze transitional states of three mixing layers: one without drops, and two that included drops at different initial mass loadings. Without drops, the dissipation is essentially due to viscous effects. It was found that in the presence of drops, the largest contribution to dissipation was made by heating and evaporation of the drops, and that at large length scales, this contribution is positive (signifying that the drops reduce turbulence), while at small scales, this contribution is negative (the drops increase turbulence). The second largest contribution to dissipation was found to be associated with the chemical potential, which leads to an increase in turbulence at large scales and a decrease in turbulence at small scales. The next smaller contribution was found to be that of viscosity. The fact that viscosity effects are only third in order of magnitude in the dissipation is in sharp contrast to the situation for the mixing layer without the drops. The next smaller contribution - that of the drag and momentum of the vapor from the drops - was found to be negative at lower mass loading but to become positive at higher mass loading.
This report presents a study of numerical simulations of mixing layers developing between opposing flows of paired fluids under supercritical conditions, the purpose of the study being to elucidate chemical-species-specific aspects of turbulence. The simulations were performed for two different fluid pairs O2/H2 and C7H16/N2 at similar reduced initial pressures (reduced pressure is defined as pressure divided by critical pressure). Thermodynamically, O2/H2 behaves more nearly like an ideal mixture and has greater solubility, relative to C7H16/N2, which departs strongly from ideality. Because of a specified smaller initial density stratification, the C7H16/N2 layers exhibited greater levels of growth, global molecular mixing, and turbulence. However, smaller density gradients at the transitional state for the O2/H2 system were interpreted as indicating that locally, this system exhibits enhanced mixing as a consequence of its greater solubility and closer approach to ideality. These thermodynamic features were shown to affect entropy dissipation, which was found to be larger for O2/H2 and concentrated in high-density-gradient-magnitude regions that are distortions of the initial density-stratification boundary. In C7H16/N2, the regions of largest dissipation were found to lie in high-density-gradient-magnitude regions that result from mixing of the two fluids.
The focus of this research is on supercritical C7H16/N2 and O2/H2 mixing layers undergoing transitions to turbulence. The C7H16/N2 system serves as a simplified model of hydrocarbon/air systems in gas-turbine and diesel engines; the O2/H2 system is representative of liquid rocket engines. One goal of this research is to identify ways of controlling area production to increase disintegration of fluids and enhance combustion in such engines. As used in this research, "area production" signifies the fractional rate of change of surface area oriented perpendicular to the mass-fraction gradient of a mixing layer. In the study, a database of transitional states obtained from direct numerical simulations of the aforementioned mixing layers was analyzed to investigate global layer characteristics, phenomena in regions of high density-gradient magnitude (HDGM), irreversible entropy production and its relationship to the HDGM regions, and mechanisms leading to area production.
The existence of methane in Titan's atmosphere (approx. 6% level at the surface) presents a unique enigma, as photochemical models predict that the current inventory will be entirely depleted by photochemistry in a timescale of approx 20 Myr. In this paper, we examine the clues available from isotopic ratios (C-12/C-13 and D/H) in Titan's methane as to the past atmosphere history of this species. We first analyze recent infrared spectra of CH4 collected by the Cassini Composite Infrared Spectrometer, measuring simultaneously for the first time the abundances of all three detected minor isotopologues: (13)CH4, (12)CH3D, and (13)CH3D. From these we compute estimates of C-12/C-13 = 86.5 +/- 8.2 and D/H = (1.59 +/- 0.33) x 10(exp -4) , in agreement with recent results from the Huygens GCMS and Cassini INMS instruments. We also use the transition state theory to estimate the fractionation that occurs in carbon and hydrogen during a critical reaction that plays a key role in the chemical depletion of Titan's methane: CH4 + C2H yields CH3 + C2H2. Using these new measurements and predictions we proceed to model the time evolution of C-12/C-13 and D/H in Titan's methane under several prototypical replenishment scenarios. In our Model 1 (no resupply of CH4), we find that the present-day C-12/C-13 implies that the CH4 entered the atmosphere 60-1600 Myr ago if methane is depleted by chemistry and photolysis alone, but much more recently-most likely less than 10 Myr ago-if hydrodynamic escape is also occurring. On the other hand, if methane has been continuously supplied at the replenishment rate then the isotopic ratios provide no constraints, and likewise for the case where atmospheric methane is increasing, We conclude by discussing how these findings may be combined with other evidence to constrain the overall history of the atmospheric methane.
In this paper, an integrated approach for gear health prognostics using particle filters is presented. The presented method effectively addresses the issues in applying particle filters to gear health prognostics by integrating several new components into a particle filter: (1) data mining based techniques to effectively define the degradation state transition and measurement functions using a one-dimensional health index obtained by whitening transform; (2) an unbiased l-step ahead RUL estimator updated with measurement errors. The feasibility of the presented prognostics method is validated using data from a spiral bevel gear case study.