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.
The calculated transition state structure of 1-iodobutane reacting with DABCO is shown. The kinetic rates of this reaction in several organic solvents and temperatures have been measured along with those of iodomethane reacting with DABCO.
Electrocatalysts encapsulated by nanoscopic overlayers can control the rate of redox reactions at the outer surface of the overlayer or at the buried interface between the overlayer and the active catalyst, leading to complex behavior in the presence of two competing electrochemical reactions. This study investigated oxide encapsulated electrocatalysts (OECs) comprised of iridium (Ir) thin films coated with an ultrathin (2–10 nm thick) silicon oxide (SiO x ) or titanium oxide (TiO x ) overlayer. The performance of SiO x |Ir and TiO x |Ir thin film electrodes towards the oxygen evolution reaction (OER) and Fe(II)/Fe(III) redox reactions were evaluated. An improvement in selectivity towards the OER was observed for all OECs. Overlayer properties, namely ionic and electronic conductivity, were assessed using a combination of electroanalytical methods and molecular dynamics simulations. SiO x and TiO x overlayers were found to be permeable to H 2 O and O 2 such that the OER can occur at the MO x |Ir (M = Ti, Si) buried interface, which was further supported with molecular dynamics simulations of model SiO 2 coatings. In contrast, Fe(II)/Fe(III) redox reactions occur to the same degree with TiO x overlayers having thicknesses less than 4 nm as bare electrocatalyst, while SiO x overlayers inhibit redox reactions at all thicknesses. This observation is attributed to differences in electronic transport between the buried interface and outer overlayer surface, as measured with through-plane conductivity measurements of wetted overlayer materials. These findings reveal the influence of oxide overlayer properties on the activity and selectivity of OECs and suggest opportunities to tune these properties for a wide range of electrochemical reactions.
Due to an accidentally large s-wave scattering length, in a relatively wide range of energy, neutrons are approximately described by the nonrelativistic conformal field theory of unitarity fermions, perturbed by one relevant and an infinite number of irrelevant operators. We develop a formalism which provides a nonperturbative definition of local operators in that nonrelativistic conformal field theory. We compute the scattering-length and effective-range corrections to the two-point functions of primary charge-three operators using the technique of conformal perturbation theory. These calculations allow us to find the first corrections to the scale-invariant behavior of the rate of nuclear reactions with three neutrons in the final state in the regime when the neutrons have small relative momenta.
Metabolic flux analysis (MFA) is a valuable tool for quantifying cellular phenotypes and to guide plant metabolic engineering. By introducing stable isotopic tracers and employing mathematical models, MFA can quantify the rates of metabolic reactions through biochemical pathways. Recent applications of isotopically nonstationary MFA (INST–MFA) to plants have elucidated nonintuitive metabolism in leaves under optimal and stress conditions, described coupled fluxes for fast–growing algae, and produced a synergistic multi–organ flux map that is a first in MFA for any biological system. Finally, these insights could not be elucidated through other approaches and show the potential of INST–MFA to correct an oversimplified understanding of plant metabolism.
Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.
TAPsolver is a python package, which automates TAP simulation and analysis routines. TAPsolver is built around the python packages FEniCS and Dolfin-Adjoint, which help take advantage of model adjoints to provide automatic derivatives. TAPsolver is flexible, with reaction mechanisms and rate constants that can be set through input files that allow users to take advantage of the different functionalities, which include sensitivity analyses, parameter optimization and uncertainty quantification.
SAND2022-12609 O PlasMod0 is a 0-dimensional plasma modeling which approximates a plasma as a well-stirred reactor. The model uses two temperatures, one for electrons and one for ions and neutrals. This model was originally developed for laser produced plasma in single-particle aerosol mass spectrometers. As a result, there are some specialized capabilities that include spherical expansion into vacuum at a user-specified rate and photon reactions. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Savannah River National Laboratory (SRNL) has completed a work scope in support of the Savannah River Remediation (SRR) Nitric-Glycolic flowsheet. Glycolic acid will replace formic acid as a reducing agent in the pre-processing of High-Level Waste sludge performed at the Defense Waste Processing Facility (DWPF). Glycolic acid will perform the same function as formic acid, namely reduction of mercury and adjustment of feed rheology and melter oxidation/reduction potential, while significantly reducing the potential for hydrogen generation in DWPF processing. Development testing has demonstrated glycolic acid virtually eliminates hydrogen production in the pre-processing steps. The Nitric-Glycolic flowsheet significantly improves DWPF’s ability to address one of the Savannah River Site’s key challenges - the incorporation of effluent received from the Salt Waste Processing Facility (SWPF). SWPF will deliver significant effluent volume to DWPF, resulting in a concurrent increase in DWPF effluent returned to the Concentration, Storage and Transfer Facilities (CSTF). This work scope demonstrates that glycolate can be destroyed under the Nitric-Glycolic flowsheet prior to returning the DWPF recycle stream to CSTF. To avoid potential flammability issues due to thermolysis of glycolate in CSTF, SRR tasked SRNL to quantify and mitigate glycolate returns via DWPF’s recycle stream. The development of a strategy for glycolic mitigation was initiated with a system’s engineering workshop. Various chemical and/or physical solutions for how and where to destroy glycolate were considered – consistent with DWPF’s operational capabilities and process requirements. The workshop identified chemical oxidation of glycolate within the DWPF Recycle Collection Tank (RCT) as the most promising option. (The RCT collects offgas condensate from pre-processing, vitrification, and other unit operations performed in DWPF and is the singular return vessel delivering recycle effluent back to CSTF.) A series of down-select experimental tests were performed that showed sodium permanganate was the most effective oxidant. This work was performed to identify an appropriate operational protocol for DWPF implementation. Building upon the downselect experiments, the work targeted an improved description of suitable RCT conditions prior to the strike (pH, etc.), permanganate addition rate, and destruction reaction kinetics. Data collected were also used to determine the best expression of reaction stoichiometry and the effective minimum permanganate addition required.
The objective of this work is to examine the consequences of a residual water on change in condition of cladding and SNF during dry storage. Small amounts of residual water may remain in SNF dry storage canisters after the fuel assemblies are transferred to the canister in the spent fuel pool (SFP), the canister is removed from the SFP and drained or pumped, vacuum dried in several pressure reducing steps, and backfilled with helium. This residual water from incomplete drying could degrade fuel rod cladding, oxidize exposed fuel pellets, and cause cladding failure resulting in SNF fines’ release in canister cavity. If residual water content is high enough, it could corrode internal components inside the canister, e.g., fuel basket, and neutron absorber plates, but CNWRA work indicated that corrosion of the internal component is expected to be negligible even with 55 moles of residual water. Considering this, focus of this work to determine partitioning of residual water towards fuel and cladding oxidation. Residual water will decompose by radiolysis and create reactive oxygen species, which, at sufficient high temperatures will oxidize cladding, and UO 2 exposed by clad breaches. If a significant fraction of UO 2 in a fuel rod is oxidized to U 3 O 8 , swelling of the fuel pellets will rupture the cladding and release fuel particles, contaminating the inside of the canister and complicating retrieval of fuel assemblies, transport and ultimate disposal. The report describes an evaluation of the potential cladding and exposed SNF pellets due to residual water. An integrated quantitative approach was used to estimate the effects. The approach consisted of (i) the evolution of environmental conditions, i.e., temperature and radiation field, in the cask, (ii) the physicochemical processes that affect degradation of the materials within the cask, and (iii) the potential damage states that may result from these conditions and processes. The assessment considers the following set of cask environmental conditions: (i) the quantity of residual water after drying, (ii) internal temperature fields and their intensity decrease over time, and (iii) the strength of the internal radiation field. The quantity of residual water will affect the mass of cladding, SNF, and other components that are degraded by oxidation. Temperature controls the rate of the oxidation reactions, the phase composition of water, and, to a lesser extent, the rate at which water vapor is radiolyzed. The strength of the radiation field controls the radiolysis rate of residual water. These processes, i.e., time evolution of temperature fields and their effects on radiolysis of residual water, fuel and cladding oxidations, were modeled using the integrated approach.
Los Alamos National Laboratory developed an ultra-high-throughput screening platform to engineer custom biocatalysts that enhance the rate of chemical reactions critical in pharmaceuticals, renewable energy, and environmental cleanup. Current methods to find biocatalysts are slow. This platform scans genetic variations to optimize a single enzyme or microbial cell to generate a product efficiently. It selects rare mutations needed for biocatalyst optimization orders of magnitude faster than current screening methods. A custom sensor reporter gene circuit causes cells to fluoresce when they are making the target product. When coupled to flow cytometry, a million biocatalyst variants can be screened in hours.
The Zel’dovich-von Neumann-Doering (ZND) profile of a detonation wave is derived. Two basic material properties are required: i. An equation of state (EOS) for a partly burned explosive; P(V, e, λ). ii. A burn rate for the reaction progress variable; d/dtλ = R(V, e, λ). For a steady planar detonation wave the reactive flow PDEs can be reduced to ODEs. The detonation wave profile can be determined from an ODE plus algebraic equations for points on the partly burned detonation loci with a specified wave speed. Furthermore, for the CJ detonation speed, the end of the reaction zone is sonic. A solution to the reactive flow equations can be constructed with a rarefaction wave following the detonation wave profile. This corresponds to an underdriven detonation wave, and the rarefaction is know as a Taylor wave.
Determining kinetic parameters such as activation energy (Ea), pre-exponential factor (A), rate constants, and reaction order (n) via thermal analysis techniques is important to material synthesis and fabrication, industrial production, biomedicals, energy storage, catalysis, etc. Various kinetic methods based on the non-isothermal techniques have been developed to obtain the kinetic parameters. In this review, we focus on introducing various kinetic methods in non-isothermal analysis and their applications in different reactions, including thermal decomposition, solid-state phase transformation, crystallization, thermal ignition, curing process, etc.
While detailed chemical kinetic models have been successful in representing rates of chemical reactions in continuum scale computational fluid dynamics (CFD) simulations, applying the models in simulations for engineering device conditions is computationally prohibitive. To reduce the cost, data-driven methods, e.g., autoencoders, have been used to construct reduced chemical kinetic models for CFD simulations. Despite their success, data-driven methods rely heavily on training data sets and can be unreliable when used in out-of-distribution (OOD) regions (i.e., when extrapolating outside of the training set). In this paper, we present an enhanced autoencoder model for combustion chemical kinetics with uncertainty quantification to enable the detection of model usage in OOD regions, and thereby creating an OOD-aware autoencoder model that contributes to more robust CFD simulations of reacting flows.Here, we first demonstrate the effectiveness of the method in OOD detection in two well-known datasets, MNIST and Fashion-MNIST, in comparison with the deep ensemble method, and then present the OOD-aware autoencoder for reduced chemistry model in syngas combustion.
The invention pertains to a PPA from a microorganism belonging to the family Halobacteriaceae (HPPA), for example, a PPA from Haloferax volcanii. The HPPA provided by the invention is soluble, thermostable and active at high concentrations of salt and/or organic solvent. An embodiment of the invention provides a method of increasing the rate of a reaction by adding an HPPA to the reaction mixture, wherein the reaction produces PPi, for example, an enzymatic reaction, and wherein the reaction is carried out at moderately high temperature and/or low water activity. Further embodiments of the invention provide an assay to detect the PPi released during a reaction which produces PPi by adding an HPPA to convert the PPi in to Pi and measuring the resultant Pi. The invention further pertains to an assay to monitor a reaction which produces PPi in the presence or the absence of an HPPA.
TAPsolver is a python package, which automates TAP simulation and analysis routines. TAPsolver is built around the python packages FEniCS and Dolfin-Adjoint, which help take advantage of model adjoints to provide automatic derivatives. TAPsolver is flexible, with reaction mechanisms and rate constants that can be set through input files that allow users to take advantage of the different functionalities, which include sensitivity analyses, parameter optimization and uncertainty quantification.
Solvent molecules interact with reactive species and alter the rates and selectivities of catalytic reactions by orders of magnitude. Specifically, solvent molecules can modify the free energies of liquid phase and surface species via solvation, participating directly as a reactant or co-catalyst, or competitively binding to active sites. These effects carry consequences for reactions relevant for the conversion of renewable or recyclable feedstocks, the development of distributed chemical manufacturing, and the utilization of renewable energy to drive chemical reactions. First, we describe the quantitative impact of these effects on steady-state catalytic turnover rates through a rate expression derived for a generic catalytic reaction (A → B), which illustrates the functional dependence of rates on each category of solvent interaction. Second, we connect these concepts to recent investigations of the effects of solvents on catalysis to show how interactions between solvent and reactant molecules at solid–liquid interfaces influence catalytic reactions. Furthermore, this discussion demonstrates that the design of effective liquid phase catalytic processes benefits from a clear understanding of these intermolecular interactions and their implications for rates and selectivities.
In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.
We describe sensitivities of elementary steps as salient parameters in determining the rate determining character of elementary steps within a reaction network and develop a formalism wherein the overall composite reaction is described by an apparent rate-determining step that is a sensitivity-weighted average of the elementary steps that comprise the reaction. Reaction parameters—apparent reaction orders, apparent enthalpy and entropy of activation—of the composite reaction network are determined within the framework we develop by application of transition state theory to the apparent rate-determining step. From this formalism we develop methods for determination of surface coverages by measuring only reaction orders, novel methods for discrimination between proposed mechanisms, a proof that the kinetic degrees of rate control sum to unity, and a proof for the relationship between fractional coverages and thermodynamic degrees of rate control. Two complementary formalisms for identifying rate-limiting transition states are broadly employed in chemical kinetics—De Donder relations based on measuring the thermodynamic driving forces of elementary steps during reaction and degrees of rate control based on knowledge of the kinetics of elementary reaction steps. Here, the formalism developed herein unifies these two strategies to elucidate in the most general case a relationship between the thermodynamic driving forces of elementary steps captured by the reversibilities and the kinetic and thermodynamic degrees of rate control.