Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transition state”

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 847 records · Page 47

Irreversible Entropy Production in Two-Phase Mixing Layers

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.

Okongo, Nora↗

Turbulence in Supercritical O2/H2 and C7H16/N2 Mixing Layers

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.

Bellan, Josette↗

Update on Area Production in Mixing of Supercritical 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.

Okongo, Nora↗

Isotopic Ratios in Titan's Methane: Measurements and Modeling

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.

Nixon, C. A.↗

An Integrated Approach for Gear Health Prognostics

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.

He, David↗

Assimilating Remote Sensing Observations of Leaf Area Index and Soil Moisture for Wheat Yield Estimates: An Observing System Simulation Experiment

Observing system simulation experiments were used to investigate ensemble Bayesian state updating data assimilation of observations of leaf area index (LAI) and soil moisture (theta) for the purpose of improving single-season wheat yield estimates with the Decision Support System for Agrotechnology Transfer (DSSAT) CropSim-Ceres model. Assimilation was conducted in an energy-limited environment and a water-limited environment. Modeling uncertainty was prescribed to weather inputs, soil parameters and initial conditions, and cultivar parameters and through perturbations to model state transition equations. The ensemble Kalman filter and the sequential importance resampling filter were tested for the ability to attenuate effects of these types of uncertainty on yield estimates. LAI and theta observations were synthesized according to characteristics of existing remote sensing data, and effects of observation error were tested. Results indicate that the potential for assimilation to improve end-of-season yield estimates is low. Limitations are due to a lack of root zone soil moisture information, error in LAI observations, and a lack of correlation between leaf and grain growth.

Nearing, Grey S.↗

Cold Water Vapor in the Barnard 5 Molecular Cloud

After more than 30 yr of investigations, the nature of gas-grain interactions at low temperatures remains an unresolved issue in astrochemistry. Water ice is the dominant ice found in cold molecular clouds; however, there is only one region where cold ((is) approximately 10 K) water vapor has been detected-L1544. This study aims to shed light on ice desorption mechanisms under cold cloud conditions by expanding the sample. The clumpy distribution of methanol in dark clouds testifies to transient desorption processes at work-likely to also disrupt water ice mantles. Therefore, the Herschel HIFI instrument was used to search for cold water in a small sample of prominent methanol emission peaks. We report detections of the ground-state transition of o-H2O (J = 110-101) at 556.9360 GHz toward two positions in the cold molecular cloud, Barnard 5. The relative abundances of methanol and water gas support a desorption mechanism which disrupts the outer ice mantle layers, rather than causing complete mantle removal.

ISM: molecules↗

State Machine Modeling of the Space Launch System Solid Rocket Boosters

The Space Launch System is a Shuttle-derived heavy-lift vehicle currently in development to serve as NASA's premiere launch vehicle for space exploration. The Space Launch System is a multistage rocket with two Solid Rocket Boosters and multiple payloads, including the Multi-Purpose Crew Vehicle. Planned Space Launch System destinations include near-Earth asteroids, the Moon, Mars, and Lagrange points. The Space Launch System is a complex system with many subsystems, requiring considerable systems engineering and integration. To this end, state machine analysis offers a method to support engineering and operational e orts, identify and avert undesirable or potentially hazardous system states, and evaluate system requirements. Finite State Machines model a system as a finite number of states, with transitions between states controlled by state-based and event-based logic. State machines are a useful tool for understanding complex system behaviors and evaluating "what-if" scenarios. This work contributes to a state machine model of the Space Launch System developed at NASA Ames Research Center. The Space Launch System Solid Rocket Booster avionics and ignition subsystems are modeled using MATLAB/Stateflow software. This model is integrated into a larger model of Space Launch System avionics used for verification and validation of Space Launch System operating procedures and design requirements. This includes testing both nominal and o -nominal system states and command sequences.

model-based systems engineering↗

Instantons in Quantum Annealing: Thermally Assisted Tunneling Vs Quantum Monte Carlo Simulations

Recent numerical result (arXiv:1512.02206) from Google suggested that the D-Wave quantum annealer may have an asymptotic speed-up than simulated annealing, however, the asymptotic advantage disappears when it is compared to quantum Monte Carlo (a classical algorithm despite its name). We show analytically that the asymptotic scaling of quantum tunneling is exactly the same as the escape rate in quantum Monte Carlo for a class of problems. Thus, the Google result might be explained in our framework. We also found that the transition state in quantum Monte Carlo corresponds to the instanton solution in quantum tunneling problems, which is observed in numerical simulations.

Quantum Monte Carlo↗

Hybrid Differential Dynamic Programming with Stochastic Search

Differential dynamic programming (DDP) has been demonstrated as a viable approach to low-thrust trajectory optimization, namely with the recent success of NASA's Dawn mission. The Dawn trajectory was designed with the DDP-based Static/Dynamic Optimal Control algorithm used in the Mystic software.1 Another recently developed method, Hybrid Differential Dynamic Programming (HDDP),2, 3 is a variant of the standard DDP formulation that leverages both first-order and second-order state transition matrices in addition to nonlinear programming (NLP) techniques. Areas of improvement over standard DDP include constraint handling, convergence properties, continuous dynamics, and multi-phase capability. DDP is a gradient based method and will converge to a solution nearby an initial guess. In this study, monotonic basin hopping (MBH) is employed as a stochastic search method to overcome this limitation, by augmenting the HDDP algorithm for a wider search of the solution space.

Aziz, Jonathan↗

Hybrid Differential Dynamic Programming with Stochastic Search

Differential dynamic programming (DDP) has been demonstrated as a viable approach to low-thrust trajectory optimization, namely with the recent success of NASAs Dawn mission. The Dawn trajectory was designed with the DDP-based Static Dynamic Optimal Control algorithm used in the Mystic software. Another recently developed method, Hybrid Differential Dynamic Programming (HDDP) is a variant of the standard DDP formulation that leverages both first-order and second-order state transition matrices in addition to nonlinear programming (NLP) techniques. Areas of improvement over standard DDP include constraint handling, convergence properties, continuous dynamics, and multi-phase capability. DDP is a gradient based method and will converge to a solution nearby an initial guess. In this study, monotonic basin hopping (MBH) is employed as a stochastic search method to overcome this limitation, by augmenting the HDDP algorithm for a wider search of the solution space.

Aziz, Jonathan↗

RACORO Continental Boundary Layer Cloud Investigations: 1. Case Study Development and Ensemble Large-Scale Forcings

Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.

Boundary layers↗

Efficient Trajectory Propagation for Orbit Determination Problems

Regularized formulations of orbital motion apply a series of techniques to improve the numerical integration of the orbit. Despite their advantages and potential applications little attention has been paid to the propagation of the partial derivatives of the corresponding set of elements or coordinates, required in many orbit-determination scenarios and optimization problems. This paper fills this gap by presenting the general procedure for integrating the state-transition matrix of the system together with the nominal trajectory using regularized formulations and different sets of elements. The main difficulty comes from introducing an independent variable different from time, because the solution needs to be synchronized. The correction of the time delay is treated from a generic perspective not focused on any particular formulation. The synchronization using time-elements is also discussed. Numerical examples include strongly-perturbed orbits in the Pluto system, motivated by the recent flyby of the New Horizons spacecraft, together with a geocentric flyby of the NEAR spacecraft.

numerical methods↗

Advanced Symbolic Analysis Tools for Fault-Tolerant Integrated Distributed Systems

The project aims to develop advanced model-checking algorithms and tools to automate the verification of fault-tolerant distributed systems for avionics. We present a new method called Property-Directed K-Induction (PD-KIND) for synthesizing K-inductive invariants of state-transition systems. PD-KIND builds upon Satifiability Modulo Theories (SMT) to generalize Bradley's IC3 method and its variants. This method is implemented in a new tool called SALLY. Case studies show that PD-KIND can automatically verify fault-tolerant algorithms under a variety of fault models and that SALLY is competitive with other SMT-based model checkers.

Dutertre, Bruno↗

Contributions of Spherical Harmonics to Gravitational Moment

A scalar gravitational potential function expressed as a series of spherical harmonics frequently serves as the basis for a model of an astronomical body's gravitational field. The contribution of a generic spherical harmonic to gravitational gradient is expressed as a dyadic, which is then used to obtain an analytical expression in vector-dyadic form for the contribution to the moment of gravitational forces about the mass center of a small body such as a spacecraft. The expression developed for a harmonic's contribution to gravitational gradient can be applied in areas beyond the scope of the paper; for example, gravitational gradient plays an important role in the state propagation matrix and the state transition matrix that are used in spacecraft trajectory targeting and Kalman filtering. Additionally, it can be employed in numerical simulations of orbit determination based on measurements obtained with a gradiometer in low-Earth orbit. Contributions of spherical harmonics to gravitational moment may be of interest in connection with attitude control of a spacecraft in the vicinity of a body with an irregular shape, such as an asteroid. Normalized spherical harmonic coefficients up to degree and order 10 are obtained for the asteroid 216 Kleopatra and used in numerical evaluations of contributions to gravitational moment.

Roithmayr, Carlos M.↗

Mars Optimal Aerobrake Maneuver Estimation

Mars science satellites often perform orbit changes to obtain different measurements, ground tracks or relay operations. Large reductions in semi-major axis and eccentricity can be done efficiently using the atmospheric drag, a.k.a aerobrake. Aerobraking is one of the most challenging planetary orbit maneuvers in terms of planning and operations. The most important consideration for aerobraking is maintaining the spacecraft's periapsis within an allocated atmospheric density corridor, which is accomplished by raising or lowering periapsis through one or a series of very small and short maneuvers. These maneuvers must be performed as efficiently as possible due to propellant constraints. Work herein details a fast and accurate method to calculate the required impulsive velocity changes in the orbit to guarantee that the spacecraft remains in a prescribed density corridor. The method makes use of the orbit's state transition matrix to map the solution space around the reference orbit. It evaluates the most efficient maneuver epochs to target a given periapsis change with a linear optimal control for single or multiple maneuvers. A fast calculation of the maneuver allows for a more comprehensive evaluation of the trade space, and the selected maneuver may be re-targeted later with a higher- fidelity model. Comparisons against fully propagated models and direct method optimizations demonstrate the new method's performance.

Sarli, Bruno V.↗

X-Ray Spectral and Flux Variability of the Microquasar GRS 1758-258 on Timescales from Weeks to Years

We present the spectral and timing evolution of the persistent black hole X-ray binary GRS 1758-258 based on almost 12 years of observations using the Rossi X-ray Timing Explorer Proportional Counter Array. While the source was predominantly found in the hard state during this time, it entered the thermally dominated soft state seven times. In the soft state GRS 1758-258 shows a strong decline in flux above 3 keV rather than the pivoting flux around 10 keV more commonly shown by black hole transients. In its 3-20 keV hardness intensity diagram, GRS 1758-258 shows a hysteresis of hard and soft state fluxes typical for transient sources in outburst. The RXTE-PCA and RXTE-ASM long-term light curves do not show any orbital modulations in the range of 2 to 30 d. However, in the dynamic power spectra significant peaks drift between 18.47d and 18.04d for the PCA data, while less significant signatures between 19d and 20d are seen for the ASM data as well as for the Swift/BAT data. We discuss different models for the hysteresis behavior during state transitions as well as possibilities for the origin of the long term variation in the context of a warped accretion disk.

Maria Hirsch↗

Adaptive Stress Testing of Trajectory Predictions in Flight Management Systems

To find failure events and their likelihoods in flight-critical systems, we investigate the use of an advanced black-box stress testing approach called adaptive stress testing. We analyze a trajectory predictor from a developmental commercial flight management system which takes as input a collection of lateral waypoints and en-route environmental conditions. Our aim is to search for failure events relating to inconsistencies in the predicted lateral trajectories. The intention of this work is to find likely failures and report them back to the developers so they can address and potentially resolve shortcomings of the system before deployment. To improve search performance, this work extends the adaptive stress testing formulation to be applied more generally to sequential decision-making problems with episodic reward by collecting the state transitions during the search and evaluating at the end of the simulated rollout. We use a modified Monte Carlo tree search algorithm with progressive widening as our adversarial reinforcement learner. The performance is compared to direct Monte Carlo simulations and to the cross-entropy method as an alternative importance sampling baseline. The goal is to find potential problems otherwise not found by traditional requirements-based testing. Results indicate that our adaptive stress testing approach finds more failures and finds failures with higher likelihood relative to the baseline approaches.

adaptive stress testing↗