Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction path analysis”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Using Membranes with Internal Microchannels to Prevent Drying-out during CO 2 Electrolysis

Scaling up CO 2 electrolysis is a vital aspect in the transition to manufacturing sustainable fuels and chemical compounds, satisfying the demand for chemicals and demand for storing renewable electricity. Depending on the employed catalyst, different products can be produced, such as carbon monoxide and ethylene, by applying a voltage on a CO 2 and H 2 O fed electrolyzer. CO is a desirable product according to techno-economic analysis, because it can be produced selectively using a silver catalyst and is a precursor for hydrocarbons in the Fischer–Tropsch process. In a state-of-the-art CO 2 electrolyzer, two electrodes are directly pressed against an ion-exchange membrane – this is called a zero-gap configuration. Therefore, the membrane is a crucial component for the system, since it has the role of providing a conductive medium between the electrodes. One of the challenges in CO 2 electrolysis is that water is consumed in the reaction. At high current density, this may cause the membrane's surface near the cathode to dry out, lowering efficiency or perhaps stopping the process entirely, since there is no longer a conductive medium. As a result, water management is critical for this process. In this work, we approach the drying-out challenge by studying the novel concept of a membrane with internal microchannels. These channels allow the circulation of water or an electrolyte inside the membrane, which reduces the water diffusion path and affects the membrane’s conductivity. The effects of channel geometry, location, and concentrations of electrolyte inside on water content, conductivity and overall performance are studied in a 2D COMSOL model. In addition, the effect of internal concentration of electrolyte on the membrane’s resistance, on the process performance and the K + cross-over to cathode side were investigated experimentally. Our modeling results prove that the presence of the channels can keep the membrane hydrated. The highest current densities are observed when the channel is closest to the cathode, and with smaller pores. Smaller pores are advantageous due to the trade-off between enhanced membrane conductivity and the lower conductivity of the liquid itself. If water is circulated in a large channel, it increases the membrane’s ionomer conductivity due to hydration but the overall conductivity is decreased since water is not highly conductive. Nonetheless, the results also show that a higher concentration of electrolyte inside the microchannels can significantly increase the total conductivity of the membrane, and therefore the energy efficiency of the process. These effects are most significant at higher current densities. In the experimental results, we’ve observed similar effects in terms of membrane conductivity and current density of the process – the higher the electrolyte concentration the higher the current density. Furthermore, a low concentration decreases the amount of potassium which crosses over to the cathode side, inhibiting salt deposition. We’ve concluded that a small channel, up to 90 µm wide, close to the pore with an electrolyte with a concentration of up to 10 mM could be very beneficial for the water management and energy efficiency of the process. This helps to keep the membrane hydrated at higher current densities, improves the conductivity of the membrane, and it doesn’t have significant impacts on the salt deposition.

Petrov, Kostadin Veselinov↗
DOE OSTI2020

Insight into the Catalytic Mechanism of GH11 Xylanase: Computational Analysis of Substrate Distortion Based on a Neutron Structure

The reaction mechanism of biomass decomposition by xylanases remains the subject of debate. Here, to clarify the mechanism, we investigated the glycosylation step of GH11 xylanase, an enzyme that catalyzes the hydrolysis of lignocellulosic hemicellulose (xylan). Making use of a recent neutron crystal structure, which revealed the protonation states of relevant residues, we used ab initio quantum mechanics/molecular mechanics (QM/MM) calculations to determine the detailed reaction mechanism of the glycosylation step. In particular, our focus is on the controversial question of whether or not an oxocarbenium ion intermediate is formed on the reaction pathway. The calculations support the validity of a basic retaining mechanism within a double-displacement scheme. The estimated free energy barrier of this reaction is ~18 kcal/mol with QM/MM-CCSD(T)/6-31(+)G**//MP2/6-31+G**/AMBER calculations, and the rate-determining step of the glycosylation is scission of the glycosidic bond after proton transfer from the acidic Glu177. The estimated lifetime of the oxocarbenium ion intermediate (on the order of tens of ps) and the secondary kinetic isotope effect suggest that there is no accumulation of this intermediate on the reaction path, although the intermediate can be transiently formed. In the enzyme–substrate (ES) complex, the carbohydrate structure of the xylose residue at the -1 subsite has a rather distorted (skewed) geometry, and this xylose unit at the active site has an apparent half-chair conformation when the oxocarbenium ion intermediate is formed. The major catalytic role of the protein environment is to orient residues that take part in the initial proton transfer. Because of a fine alignment of catalytic residues, the enzyme can accelerate the glycosylation reaction without paying a reorganization energy penalty.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Tailoring the Interfacial Interactions of Porous Polymer Membranes to Accelerate Atomic Layer Deposition: The Latent Path to Antifouling Membranes

Atomic layer deposition (ALD) is a powerful strategy to engineer hybrid organic-inorganic membranes with emergent functionalities. The combination of atomic-level thickness control, a wide materials palette, and unprecedented conformality allows the physiochemical properties (e.g., hydrophilicity) of mesoporous polymer membranes to be precisely tuned. The nucleation of ALD material growth on polymer surfaces relies on Lewis acid-base interactions and remains an overlooked motif with tremendous potential to accelerate ALD nucleation and growth. Strategies to enhance these interactions could enable desirable properties such as antifouling behavior to be imparted on inert polymer surfaces that lack the necessary functional groups for ALD nucleation. Here, in this study, we demonstrate that the reactivity of polyacrylonitrile (PAN) membranes toward ALD metal oxide (MO) precursors with Lewis acid characteristics is enhanced by introducing strong Lewis base functional groups (amidoxime: Am) on the PAN backbone (Am-PAN). The resulting Lewis acid-base interactions accelerate the MO nucleation in Am-PAN and reduce the number of deposition cycles required to achieve hydrophilicity compared with the untreated PAN membrane. Unveiling the reaction mechanism, in situ Fourier transform infrared (FTIR) spectroscopy measurements established enhanced interaction dynamics between the ALD MO precursors and the Am-PAN membrane, unlike the PAN membrane. Spectroscopic ellipsometry and thermogravimetric analysis measurements revealed higher MO loadings in Am-PAN membranes compared to PAN membranes for the same number of ALD cycles. We found that strong Lewis acid-base interactions accelerated the ALD for a range of materials including Al 2 O 3 , TiO 2 , SnO 2 , and ZnO. More broadly, our work demonstrates that tailoring metal-precursor-polymer interactions is a powerful strategy to accelerate and modulate the ALD. We used this design strategy to fabricate Al 2 O 3 -Am-PAN hybrid membranes that showed 2-fold higher antifouling capability compared to pristine PAN membranes prepared with an equivalent number of Al 2 O 3 ALD cycles. Our approach expands the scope of design options for fouling-resistant porous hybrid inorganic-organic membranes and may ultimately reduce the operational costs of water treatment.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Space Shuttle Columbia Reconstruction Effort: An Analysis of the Thermal Protection System.

Upon the breakup of the Space Shuffle Columbia, which occurred during peak heating, at an altitude of 200,000 feet, and traveling at Mach 19, the daunting responsibility was placed upon National Aeronautics and Space Administration (NASA) to determine the root cause of the accident. One of the foremost tasks was to begin locating Columbia's debris, which was scattered acr ss a 2,000 square miles area of Texas and into western Louisiana. Each potential shuttle debris piece found was categorized and the longitude and latitude noted. This provided crucial data to the vehicle's breakup trajectory path and aided in locating key shuffle components. Since the airframe and tiles would be the primary focus during the reconstruction effort at the Kennedy Space center in Florida, packing and shipping was crucial to the preservation of the shape and condition of each recovered object. Over 83,000 pieces of Columbia, in a five-month timeframe, were located. The recovered tiles still bonded to structure indicated failures modes due to heating and aerodynamic loads. Also, different types of tile bond failure modes were observed, with adhesion failures at the Koropon layer or material failure at the tile's densification layer. Numerous tiles had the identification markings burned off, slumping of the Reaction-Cured Glass (RCG) coating, and of the Ll-900 silica base material. The conditions in which the tiles were found resulted in developing methods to aid in the identification of the tile. Correctly identifying each tile was critical to understanding the heating patterns of the tiles in relationship to the vehicle, since it was known that the lower left wing received the initial external tank foam impact. With the aid of multiple disciplines, the clues found in the TPS debris, lead to the root cause of the breakup of Columbia.

Micklos, Ann M.↗
DOE OSTI2020

Experimental and modeling analysis of detonation in circular arcs of the conventional high explosive PBX 9501

We examine the diffraction dynamics of a two-dimensional (2D) detonation in a circular arc of the conventional HMX-based, high performance, solid explosive PBX 9501, for which the detonation reaction zone length scale is estimated to be of the order of 100–150 µm. In this configuration, a steady propagating detonation will develop, sweeping around the arc with constant angular speed. We report on results from three PBX 9501 arc experiments, exploring the variation in linear speed on the inner and outer arc surfaces for the steady wave along with the structure of the curved detonation front, as a function of varying inner surface radius and arc thickness. Comparisons of the properties of the motion of the steady wave for each arc configuration are then made with a spatially-distributed PBX 9501 reactive burn model, calibrated to detonation performance properties in a 2D planar slab geometry. We show that geometry-induced curvature of the detonation near the inner arc surface has a significant effect on the detonation motion even for conventional high explosives. We also examine the detonation driving zone structure for each arc case, and thus the subsonic regions of the flow that determine the influence of the arc geometry on the detonation propagation. In addition, streamline paths and reaction progress isolines are calculated. Overall, we conclude that a common approximation for modeling conventional high explosive detonation, wherein the shock-normal detonation speed is assumed equal to the Chapman–Jouguet speed, can lead to significant errors in describing the speed at which the detonation propagates.

42 ENGINEERING↗
NASA NTRS2004

Characterizing the effects of droplines on target acquisition performance on a 3-D perspective display

The present study investigated the effects of droplines on target acquisition performance on a 3-D perspective display in which participants were required to move a cursor into a target cube as quickly as possible. Participants' performance and coordination strategies were characterized using both Fitts' law and acquisition patterns of the 3 viewer-centered target display dimensions (azimuth, elevation, and range). Participants' movement trajectories were recorded and used to determine movement times for acquisitions of the entire target and of each of its display dimensions. The goodness of fit of the data to a modified Fitts function varied widely among participants, and the presence of droplines did not have observable impacts on the goodness of fit. However, droplines helped participants navigate via straighter paths and particularly benefited range dimension acquisition. A general preference for visually overlapping the target with the cursor prior to capturing the target was found. Potential applications of this research include the design of interactive 3-D perspective displays in which fast and accurate selection and manipulation of content residing at multiple ranges may be a challenge.

Time and Motion Studies↗
NASA NTRS1995

Scramjet Performance Assessment Using Water Absorption Diagnostics (U)

Simultaneous multiple path measurements of temperature and H2O concentration will be presented for the AIMHYE test entries in the NASA Ames 16-Inch Shock Tunnel. Monitoring the progress of high temperature chemical reactions that define scramjet combustor efficiencies is a task uniquely suited to nonintrusive optical diagnostics. One application strategy to overcome the many challenges and limitations of nonintrusive measurements is to use laser absorption spectroscopy coupled with optical fibers. Absorption spectroscopic techniques with rapidly tunable lasers are capable of making simultaneous measurements of mole fraction, temperature, pressure, and velocity. The scramjet water absorption diagnostic was used to measure combustor efficiency and was compared to thrust measurements using a nozzle force balance and integrated nozzle pressures to develop a direct technique for evaluating integrated scramjet performance. Tests were initially performed with a diode laser tuning over a water absorption feature at 1391.7 nm. A second diode laser later became available at a wavelength near 1343.3 nm covering an additional water absorption feature and was incorporated in the system for a two-wavelength technique. Both temperature and mole fraction can be inferred from the lineshape analysis using this approach. Additional high temperature spectroscopy research was conducted to reduce uncertainties in the scramjet application. The lasers are optical fiber coupled to ports at the combustor exit and in the nozzle region. The output from the two diode lasers were combined in a single fiber, and the resultant two-wavelength beam was subsequently split into four legs. Each leg was directed through 60 meters of optical fiber to four combustor exit locations for measurement of beam intensity after absorption by the water within the flow. Absorption results will be compared to 1D combustor analysis using RJPA and nozzle CFD computations as well as to data from a nozzle metric balance measuring thrust and integrated pressure measurements along the length of the nozzle. Assessment of its value as a combustor performance evaluation tool will be conducted.

Cavolowsky, John A.↗
DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Pu(IV) quantification via visible–near-infrared absorption spectroscopy: tackling interferences using D-optimal design and partial least squares

Here, this study presents a novel analytical approach for quantifying Pu(IV) in glove box environments using fiber-optic-based visible–near-infrared absorption spectroscopy in combination with partial least squares regression (PLSR) and design of experiments. The method addresses significant challenges posed by overlapping spectral features arising from Nd(III), which is a common fission product impurity, and the speciation variability of Pu(IV) nitrato complexes in HNO 3 concentrations ranging from 2.5 to 11 M. A curated training set consisting of data from 20 samples was developed via D-optimal design to enable robust PLSR model calibration for Pu(IV) using the near-infrared band near 1050 nm. The training set was acquired from samples in cuvettes with a 1-cm path length and was used to build the PLSR model. The robustness of the model was validated with data collected using a dip probe with a 1-cm path length and varying Pu(IV) concentrations. The strong performance of the model indicates good model transfer from cuvette to dip probe and highlights the potential for in situ measurements and online monitoring of reactions in a crystallization reactor vessel. The results demonstrate that this combined spectroscopic and chemometric approach can accurately and simultaneously quantify Pu(IV) and HNO 3 , thereby offering a promising tool for real-time monitoring in process environments.

Actinide↗
DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Chondrites and the Protoplanetary Disk, Part 3

Contents include the following: Ca-, Al-Rich Inclusions and Ameoboid Olivine Aggregates: What We Know and Don t Know About Their Origin. Aluminium-26 and Oxygen Isotopic Distributions of Ca-Al-rich Inclusions from Acfer 214 CH Chondrite. The Trapping Efficiency of Helium in Fullerene and Its Implicatiion to the Planetary Science. Constraints on the Origin of Chondritic Components from Oxygen Isotopic Compositions. Role of Planetary Impacts in Thermal Processing of Chondrite Materials. Formation of the Melilite Mantle of the Type B1 CAIs: Flash Heating or Transport? The Iodine-Xenon System in Outer and Inner Portions of Chondrules from the Unnamed Antarctic LL3 Chondrite. Nucleosynthesis of Short-lived Radioactivities in Massive Stars. The Two-Fluid Analysis of the Kelvin-Helmholtz Instability in the Dust Layer of a Protoplanetary Disk: A Possible Path to the Planetesimal Formation Through the Gravitational Instability. Shock-Wave Heating Model for Chonodrule Formation: Heating Rate and Cooling Rate Constraints. Glycine Amide Hydrolysis with Water and OH Radical: A Comparative DFT Study. Micron-sized Sample Preparation for AFM and SEM. AFM, FE-SEM and Optical Imaging of a Shocked L/LL Chondrite: Implications for Martensite Formation and Wave Propagation. Infrared Spectroscopy of Chondrites and Their Components: A Link Between Meteoritics and Astronomy? Mid-Infrared Spectroscopy of CAI and Their Mineral Components. The Origin of Iron Isotope Fractionation in Chondrules, CAIs and Matrix from Allende (CV3) and Chainpur (LL3) Chondrites. Protoplanetary Disk Evolution: Early Results from Spitzer. Kinetics of Evaporation-Condensation in a Melt-Solid System and Its Role on the Chemical Composition and Evolution of Chondrules. Oxygen Isotope Exchange Recorded Within Anorthite Single Crystal in Vigarano CAI: Evidence for Remelting by High Temperature Process in the Solar Nebula. Chondrule Forming Shock Waves in Solar Nebula by X-Ray Flares. Organic Globules with Anormalous Nitrogen Isotopic Compositions in the Tagish Lake Meteorite: Products of Primitive Organic Reactions. Yet Another Chondrule Formation Scenario. CAIs are Not Supernova Condensates. Microcrystals and Amorphous Material in Comets and Primitive Meteorites: Keys to Understanding Processes in the Early Solar System. A Nearby Supernova Injected Short-lived Radionuclides into Our Protoplanetary Disk. REE+Y Systematics in CC and UOC Chondrules. Meteoritic Constraints on Temperatures, Pressures, Cooling Rates, Chemical Compositions, and Modes of Condensation in the Solar Nebula. The I-Xe Record of Long Equilibration in Chondrules from the Unnamed Antarctic Meteorite L3/LL3. Early Stellar Evolution.

Source record↗
DOE OSTI2025

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗
DOE OSTI2021

Get on the BAND Wagon: a Bayesian framework for quantifying model uncertainties in nuclear dynamics

In this paper, we describe the Bayesian Analysis of Nuclear Dynamics (BAND) framework, a cyberinfrastructure that we are developing which will unify the treatment of nuclear models, experimental data, and associated uncertainties. We overview the statistical principles and nuclear-physics contexts underlying the BAND toolset, with an emphasis on Bayesian methodology's ability to leverage insight from multiple models. In order to facilitate understanding of these tools we provide a simple and accessible example of the BAND framework's application. Four case studies are presented to highlight how elements of the framework will enable progress on complex, far-ranging problems in nuclear physics. By collecting notation and terminology, providing illustrative examples, and giving an overview of the associated techniques, this paper aims to open paths through which the nuclear physics and statistics communities can contribute to and build upon the BAND framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2008

Tolerance and UQ4SIM: Nimble Uncertainty Documentation and Analysis Software

Ultimately, scientific numerical models need quantified output uncertainties so that modeling can evolve to better match reality. Documenting model input uncertainties and variabilities is a necessary first step toward that goal. Without known input parameter uncertainties, model sensitivities are all one can determine, and without code verification, output uncertainties are simply not reliable. The basic premise of uncertainty markup is to craft a tolerance and tagging mini-language that offers a natural, unobtrusive presentation and does not depend on parsing each type of input file format. Each file is marked up with tolerances and optionally, associated tags that serve to label the parameters and their uncertainties. The evolution of such a language, often called a Domain Specific Language or DSL, is given in [1], but in final form it parallels tolerances specified on an engineering drawing, e.g., 1 +/- 0.5, 5 +/- 10%, 2 +/- 10 where % signifies percent and o signifies order of magnitude. Tags, necessary for error propagation, can be added by placing a quotation-mark-delimited tag after the tolerance, e.g., 0.7 +/- 20% 'T_effective'. In addition, tolerances might have different underlying distributions, e.g., Uniform, Normal, or Triangular, or the tolerances may merely be intervals due to lack of knowledge (uncertainty). Finally, to address pragmatic considerations such as older models that require specific number-field formats, C-style format specifiers can be appended to the tolerance like so, 1.35 +/- 10U_3.2f. As an example of use, consider figure 1, where a chemical reaction input file is has been marked up to include tolerances and tags per table 1. Not only does the technique provide a natural method of specifying tolerances, but it also servers as in situ documentation of model uncertainties. This tolerance language comes with a utility to strip the tolerances (and tags), to provide a path to the nominal model parameter file. And, as shown in [1], having the ability to quickly mark and identify model parameter uncertainties facilitates error propagation, which in turn yield output uncertainties.

Kleb, Bil↗
NASA NTRS2015

Simulation of Liquid Rocket Engine Failure Propagation Using Self-Evolving Scenarios

Traditional probabilistic risk assessment approaches often require failure scenarios to be explicitly defined through event sequences that are then quantified as part of the integrated analysis. This approach becomes difficult when failure propagation paths change as a function of the system operation. Additionally, if the propagation paths represent interactions among even a modest number of components, the scenario count becomes combinatorially intractable. This paper presents an alternate approach for quantifying the probability of failure propagation in such a case. Rather than explicitly defining scenario sequences, simple physical models are created for each of the components. In this way, only the physical states and rules of component interactions must be defined, rather than event sequences for each individual scenario. Initiating failures are introduced into the system, either randomly or as defined by relative likelihood, and the failures cascade through the system via the interaction rules. This process is repeated using Monte Carlo methods and, as a result, the most probable scenarios “self-evolve” in terms of both sequence path and frequency. This approach was applied to failures occurring in the engine compartment of a space launch vehicle with four liquid rocket engines and four high-pressure helium tanks. Each engine was modeled with key components, such as turbomachinery, combustion chamber, propellant lines, and additional support systems. Three test cases were conducted with different high-energy engine failures. End results of interest included an additional engine-out failure and tank burst, which represent the loss-of-mission (LOM) and loss-of-crew (LOC) failure environments, respectively. Observations show that almost every scenario outcome is unique and that many scenarios involve complex chain reactions that are difficult to predict. This validates the usefulness of the modeling approach in assessing the overall risks to the crew during a launch vehicle abort.

Rocket Engine↗
DOE OSTI2020

Analysis of an all-solid state nanobattery using molecular dynamics simulations under an external electric field

Present Li-ion battery (LIB) technology requires strong improvements in performance, energy capacity, charging-time, and cost to expand their application to e-mobility and grid storage. Li-metal is one of the most promising materials to replace commercial anodes such as graphite because of its 10 times higher specific capacity. However, Li-metal has high reactivity with commercial liquid electrolytes; thus, new solid materials are proposed to replace liquid electrolytes when Li-metal anodes are used. We present a theoretical analysis of the charging process in a full nanobattery, containing a LiCoO 2 cathode, a Li 7 P2S 8 I solid-state electrolyte (SSE), a Li-metal anode as well as Al and Cu collectors for the cathode and anode, respectively. In addition, we added a Li 3 P/Li 2 S film as a solid electrolyte interphase (SEI) layer between the Li-anode and SSE. Thus, we focus this study on the SEI and SSE. We simulated the charging of the nanobattery with an external voltage by applying an electric field. We estimated temperature profiles within the nanobattery and analyzed Li-ion transport through the SSE and SEI. Here, we observed a slight temperature rise at the SEI due to reactions forming $PS_{3}^{–}$ and $P_{2}S_{7}^{4}$$^{–}$ fragments at the interfaces; however, this temperature profile changes due to the charging current under the presence of the external electric field ε = 0.75 V Å –1 . Without the external field, the calculated open-circuit voltage (OCV) was 3.86 V for the battery, which is within the range of values of commercial cobalt-based LIBs. This voltage implies a spontaneous fall of available Li-ions from the anode to the cathode (during discharge). The charge of this nanobattery requires overcoming the OCV plus an additional voltage that determines the charging current. Thus, we applied an external potential able to neutralize the OCV, plus an additional 1.6 V to induce the transport of Li + from the cathode up to the anode. Several interesting details about Li + transport paths through the SSE and SEI are discussed.

25 ENERGY STORAGE↗
NASA NTRS1993

Compressibility effects on rotor forces in the leakage path between a shrouded pump impeller and its housing

A modified approach to Childs' previous work on fluid-structure interaction forces in the leakage path between an impeller shroud and its housing is presented in this paper. Three governing equations consisting of continuity, path-momentum, and circumferential-momentum equations were developed to describe the leakage path inside a pump impeller. Radial displacement perturbations were used to solve for radial and circumferential force coefficients. In addition, impeller-discharge pressure disturbances were used to obtain pressure oscillation responses due to precessing impeller pressure wave pattern. Childs' model was modified from an incompressible model to a compressible barotropic-fluid model (the density of the working fluid is a function of the pressure and a constant temperature only). Results obtained from this model yielded interaction forces for radial and circumferential force coefficients. Radial and circumferential forces define reaction forces within the impeller leakage path. An acoustic model for the same leakage path was also developed. The convective, Coriolis, and centrifugal acceleration terms are removed from the compressible model to obtain the acoustics model. A solution due to impeller discharge pressure disturbances model was also developed for the compressible and acoustics models. The results from these modifications are used to determine what effects additional perturbation terms in the compressible model have on the acoustic model. The results show that the additional fluid mechanics terms in the compressible model cause resonances (peaks) in the force coefficient response curves. However, these peaks only occurred at high values of inlet circumferential velocity ratios greater than 0.7. The peak pressure oscillation was shown to occur at the wearing ring seal. Introduction of impeller discharge disturbances with n = 11 diametral nodes showed that maximum peak pressure oscillations occurred at nondimensional precession frequencies of f = 6.4 and f = 7.8 for this particular pump. Bolleter's results suggest that for peak pressure oscillations to occur at the wearing ring seal, the nondimensional excitation frequency should be on the order of f = 2.182 for n = 11. The resonances found in this research do not match the excitation frequencies predicted by Bolleter. At the predicted peak excitation frequencies given by Bolleter, the compressible model shows an attenuation of the pressure oscillations at the seal exit. The compressibility of the fluid does not have a significant influence on the model at low values of nondimensional excitation frequency. At high values of nondimensional frequency, the effects of compressibility become more significant. For the acoustic analysis, the convective, Coriolis, and centrifugal acceleration terms do affect the results to a limited extent for precession excitation and to a large extent for a pressure excitation when the fluid operates at relatively high Mach numbers.

Cao, Nhai The↗
NASA NTRS2010

Physics-Based Methods of Failure Analysis and Diagnostics in Human Space Flight

The Integrated Health Management (IHM) for the future aerospace systems requires to interface models of multiple subsystems in an efficient and accurate information environment at the earlier stages of system design. The complexity of modern aeronautic and aircraft systems (including e.g. the power distribution, flight control, solid and liquid motors) dictates employment of hybrid models and high-level reasoners for analysing mixed continuous and discrete information flow involving multiple modes of operation in uncertain environments, unknown state variables, heterogeneous software and hardware components. To provide the information link between key design/performance parameters and high-level reasoners we rely on development of multi-physics performance models, distributed sensors networks, and fault diagnostic and prognostic (FD&P) technologies in close collaboration with system designers. The main challenges of our research are related to the in-flight assessment of the structural stability, engine performance, and trajectory control. The main goal is to develop an intelligent IHM that not only enhances components and system reliability, but also provides a post-flight feedback helping to optimize design of the next generation of aerospace systems. Our efforts are concentrated on several directions of the research. One of the key components of our strategy is an innovative approach to the diagnostics/prognostics based on the real time dynamical inference (DI) technologies extended to encompass hybrid systems with hidden state trajectories. The major investments are into the multiphysics performance modelling that provides an access of the FD&P technologies to the main performance parameters of e.g. solid and liquid rocket motors and composite materials of the nozzle and case. Some of the recent results of our research are discussed in this chapter. We begin by introducing the problem of dynamical inference of stochastic nonlinear models and reviewing earlier results. Next, we present our analytical approach to the solution of this problem based on the path integral formulation. The resulting algorithm does not require an extensive global search for the model parameters, provides optimal compensation for the effects of dynamical noise, and is robust for a broad range of dynamical models. In the following Section the strengths of the algorithm are illustrated illustrated by inferring the parameters of the stochastic Lorenz system and comparing the results with those of earlier research. Next, we discuss a number of recent results in application to the development of the IHM for aerospace system. Firstly, we apply dynamical inference approach to a solution of classical three tank problems with mixed unknown continuous and binary parameters. The problem is considered in the context of ground support system for filling fuel tanks of liquid rocket motors. It is shown that the DI algorithm is well suited for successful solution of a hybrid version of this benchmark problem even in the presence of additional periodic and stochastic perturbation of unknown strength. Secondly, we illustrate our approach by its application to an analysis of the nozzle fault in a solid rocket motor (SRM). The internal ballistics of the SRM is modelled as a set of one-dimensional partial differential equations coupled to the dynamics of the propellant regression. In this example we are specifically focussed on the inference of discrete and continuous parameters of the nozzle blocking fault and on the possibility of an application of the DI algorithm to reducing the probability of "misses" of an on-board FD&P for SRM. In the next section re-contact problem caused by first stage/upper stage separation failure is discussed. The reaction forces imposed on the nozzle of the upper stage during the re-contact and their connection to the nozzle damage and to the thrust vector control (TVC) signal are obtained. It is shown that transient impact induced torquean be modelled as a response of an effective damped oscillator. A possible application of the DI algorithm to the inference of damage parameters and predicting fault dynamics ahead of time using the actuator signal is discussed. Finally, we formulate Bayesian inferential framework for development of the IHM system for in-flight structural health monitoring (SHM) of composite materials. We consider the signal generated by piezoelectric actuator mounted on composite structure generating elastic waves in it. The signal received by the sensor is than compared with the baseline signal. The possibility of damage inference is discussed in the context of development of the SHM.

Smelyanskiy, Vadim N.↗