SEARCH · Search NASA
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.
Zero-RK is a software package that simulates chemically reacting systems using sparse, preconditioned, adaptive matrix methods to achieve orders-of-magnitude reduction in simulation time while maintaining accurate results.
Zero-RK is a software package that simulates chemically reacting systems using sparse, preconditioned, adaptive matrix methods to achieve orders-of-magnitude reduction in simulation time while maintaining accurate results
Zero-RK is a software package that simulates chemically reacting systems using sparse, preconditioned, adaptive matrix methods to achieve orders-of-magnitude reduction in simulation time while maintaining accurate results.
Zero-RK is a software package that simulates chemically reacting systems using sparse, preconditioned, adaptive matrix methods to achieve orders-of-magnitude reduction in simulation time while maintaining accurate results.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS), as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Gamma Technologies, LLC (GT or Participante), to incorporate the ability to access LLNL chemical kinetics technologies while using GT-SUITE, GT’s market leading engine simulation software. At the end of the project, LLNL has released Zero-RK version 3.5 with zero- and one-dimensional (0-D and 1-D) solver functionality that interfaces with GT’s GT-SUITE v2023 and later releases. GT has tested its product to assure their customers that the interface can provide reduction in chemistry solution time for detailed chemistry simulations. The process has also positioned GT to easily benefit from future improvements of the Zero-RK suite of tools developed under the DOE Vehicle Technologies Office.
Explore the source record for details and available documents.
Slide deck for research highlight summaries for Thrust 2.
The successful commercialization of electrochemical energy-conversion systems hinges on a deeper understanding of electrocatalytic reaction kinetics. Despite extensive research, a key descriptor that characterizes electrolyte effects on reaction kinetics remains elusive. Here, surface charge in electrical double layers (EDLs) is introduced as a descriptor for electrolyte-dependent kinetics. The surface charge is calculated with a continuum EDL model parameterized by density-functional theory. The model is validated by reproducing the anomalously low slope of Pt(111) in Parsons-Zobel plots. Strong correlations are observed between calculated surface charge and experimental kinetic currents for hydrogen evolution, oxygen reduction, and CO 2 -reduction reactions across various pH levels and cationic species. These correlations can be either promotional or inhibitory, depending on solute-intermediate interactions. In acidic media, incorporating adsorbate charge captures specific adsorption effects in oxygen reduction reaction. In conclusion, these findings establish surface charge density as a key descriptor for electrolyte-dependent kinetics, which will guide the design of the electrode/electrolyte interface.
In this work, kinetic measurements for silica-supported Pt and Pt-Mo catalysts were collected in vapor-phase acetic acid hydrodeoxygenation by varying the hydrogen partial pressure between 18 and 72 kPa and acetic acid partial pressure between 7 and 18 kPa at 423-473 K. At all testing conditions, the Pt-Mo catalyst was more active and selective. In addition, the apparent activation energy for Pt-Mo of 76±3 kJ/mol was lower than that of 84±4 kJ/mol for Pt. The apparent reaction orders were also different. The order in hydrogen of 0.8±0.1 for Pt-Mo changed to zero at higher hydrogen partial pressures while that for Pt remained constant at 0.6±0.1. A near-zero order in acetic acid for Pt-Mo changed to -2.1 at higher acetic acid pressures while that for Pt remained constant at -2.9±0.3. These differences in reaction kinetics as well as in selectivity trends with changes in temperature and feed composition indicated a change in the reaction mechanism for Pt-Mo. The catalysts were characterized with hydrogen temperature programmed desorption, oxygen temperature programmed oxidation and transmission electron microscopy with energy-dispersed X-ray spectroscopy elemental mapping. Mo was present in the form of subnanometer-size clusters on the surface of Pt nanoparticles. Both Pt and Pt-Mo catalysts were stable under the reaction conditions for 10 h, and the size and structure of Pt and Pt-Mo particles remained mostly unchanged, without coke accumulation. Density functional theory calculations show that surface acetate is not a major reaction intermediate on both Pt and Pt-Mo and, instead, C-OH bond splitting with the formation of acetyl is the first reaction step in hydrodeoxygenation. The activation energy for this step is dramatically lower on Pt-Mo. Furthermore, the activity of acetyl on Pt-Mo is different. Unlike on pure Pt, the reaction of C-O bond splitting becomes exothermic with a lower activation energy on Pt-Mo, and the reaction of C-C bond splitting, in contrast, becomes endothermic with a higher activation energy, explaining the experimentally observed higher activities and selectivities of Pt-Mo. In addition, the calculations demonstrate that a pair of Pt-Mo surface atoms acts as a single active site where Mo serves as a preferential binding anchor for O atoms. The presence of Mo atoms changes the structure and energy of adsorbed and reacting surface species.
Distribution of relaxation time (DRT) is used to interpret electrochemical impedance spectroscopy (EIS) for proton exchange membrane (PEM) water electrolyzers, with an attempt to separate overlapped relaxation processes in Nyquist plots. By varying operating conditions and catalyst loadings, four main relaxation peaks arising from EIS can be identified and successfully separated from low to high frequencies as (P1) mass transport, (P2) oxygen evolution reaction kinetics, (P3) reaction kinetics (with faster time constant than P2), and (P4) ionic transport. Here, the shape, height, and frequency of the DRT peaks change with different membrane electrode assembly (MEA) configurations. Electron microscopy reveals distinct features from the cross-sectioned MEAs which verify critical DRT results in that increasing the iridium (Ir)-anode loading from 0.2 mgIr/cm 2 to 1.5 mgIr/cm 2 reduces kinetic losses due to higher site-access; a thick and compacted anode, however, also triggers higher ohmic resistances from membrane/catalyst layer hydration and increases transport losses due to longer ionomer pathways. DRT provides higher resolution to EIS for deconvoluting processes with different relaxation times and the quantification of DRT peaks improves the accounting of total losses from each process.
he reaction of ground state methylidyne (CH) with water vapor (H 2 O) is theoretically re-investigated using high-level coupled cluster computations in combination with semi-classical transition state theory (SCTST) and two-dimensional master equation simulations. Insertion of CH into a H–O bond of H 2 O over a submerged barrier via a well-skipping mechanism yielding solely H and CH 2 O is characterized. The reaction kinetics is effectively determined by the formation of a pre-reaction van der Waals complex (PRC, HC—OH 2 ) and its subsequent isomerization to activated CH 2 OH in competition with PRC re-dissociation. The tunneling effects are found to be minor, while variational effects in the PRC → CH 2 OH step are negligible. The calculated rate coefficient k ( T ) is nearly pressure-independent, but strongly depends on temperature with pronounced down-up behavior: a high value of 2 x 10 -10 cm 3 s -1 at 50 K, followed by a fairly steep decrease down to 8 × 10 -12 cm 3 s -1 at 900 K, but increasing again to 5 × 10 -11 cm 3 s -1 at 3500 K. Over the T-range of this work, k ( T ) can be expressed as: k ( T , P = 0) = 2.31 × 10 -11 ( T /300 K) -1.615 exp(-38.45/ T ) cm 3 s -1 for T = 50-400 K k ( T , P = 0) = 1.15 × 10 -12 ( T /300 K) 0.8637 exp(892.6/ T ) cm 3 s -1 for T = 400-1000 K k ( T , P = 0) = 4.57 × 10 -15 ( T /300 K) 3.375 exp(3477.4/ T ) cm 3 s -1 for T = 1000–3500 K.
The ability of two-photon lithography (TPL) to print deterministic nanoporous 3D structures is highly valuable for many applications. However, it is challenging to print such structures rapidly due to the proximity effects that cause closely spaced features to enlarge and merge. A key challenge is the limited understanding of the origins and spatiotemporal dynamics of the long-range proximity effects that extend beyond the optical focal spot. Here, we empirically investigate these proximity effects in serial TPL using custom-made acrylate-based photoresists. We demonstrate that the complex spatiotemporal dynamics of the proximity effects can be explained through the kinetics of the reaction–diffusion photopolymerization mechanisms that underlie the curing process. Specifically, we show that the proximity effects arise due to comparable timescales of reaction and diffusion of oxygen. Furthermore, we demonstrate that long-range proximity effects can be suppressed by introducing phenolic inhibitors, which reduce oxygen consumption through non-inhibiting side reactions and promote faster termination of curing reactions. These insights enabled improving the linewidths from >400 to 260 nm during printing of nanoporous 3D woodpiles at scanning speeds of 50 mm s −1 . Thus, the knowledge generated here can be applied to deterministically tune the proximity effects and enable rapid printing of high-resolution nanoporous 3D structures.
Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.