Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Rates”

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 775 records · Page 43

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

Characterizing and Engineering Trehalose Biosynthesis Pathways to Improve Acetate Tolerance in Pseudomonas putida KT2440

Alkaline pretreatment can depolymerize lignin into a diverse mixture of monomers and smaller compounds, enabling bioconversion of these compounds to value-added chemicals in engineered bacterial hosts such as Pseudomonas putida KT2440. Acetate and salts comprise approximately 8% and 32% of the lignin stream, respectively. High concentrations of these components can cause cell death in P. putida, so lignin-derived streams must be fed in lower amounts to decrease stress. Investigation of the mechanisms of bacterial tolerance to chemicals in lignin-rich substrates is therefore required to overcome this obstacle in lignin bioconversion. Previous work has shown that trehalose, an endogenously produced disaccharide, is involved in tolerance to heat, cold, and osmotic stress, as well as implicated in acid stress in other bacteria. As such, trehalose biosynthesis presents a useful target to increase acetate and osmotic tolerance in P. putida. This work investigates the native TreSA, TreSB, and TreY/TreZ trehalose biosynthesis pathways in this bacterial strain as well as the OtsAB pathway of Sphingobium sp. SYK-6. Overexpression of the OtsAB pathway in P. putida increased tolerance to acetate relative to wild-type. Additionally, functional knockouts of treSA, treSB, treY, and treZ in P. putida showed that interruption of any trehalose biosynthesis pathway, and, in the case of TreY/TreZ, interruption of the pathway at any point, significantly hindered growth in acetate as well as on glucose and LB broth. Each trehalose biosynthesis pathway was then engineered for overexpression in P. putida to determine whether overexpression of these genes conferred additional tolerance to acetate. The TreSA, TreSB, OtsA, and OtsB enzymes were also expressed and purified to determine the reaction rates and kinetics of these pathways. The hypothetical glycoside hydrolase family 15 protein, encoded in the gene located in between otsB and otsA in the Sphingobium operon, was purified as well to investigate its potential involvement in the OtsAB pathway. The findings of this work illuminated the mechanisms of trehalose biosynthesis in P. putida and identified genetic targets to improve strain tolerance to common components of alkaline-pretreated lignin streams. These findings could ultimately improve yields of desired products from lignin in engineered strains of P. putida.

alkaline pretreatment↗

Criticality Experiments to Reduce Compensating Errors in Plutonium Nuclear Data

Compensating errors between nuclear data observables in a library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors in nuclear data. A new criticality experiment, described in this work, was designed with the specific target nuclear data of 239 Pu fission, inelastic scattering, elastic scattering, capture, nu-bar, and prompt fission neutron spectrum (PFNS). This work will focus on the design and execution of the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). The criticality experiment includes two different configurations with very different geometries: one is cube-like to minimize neutron leakage while the other is slab-like to maximize leakage. Having these two widely varying configurations allows the scattering sensitivities of 239 Pu to the neutron multiplication factor to be greatly changed while minimally impacting the other cross section sensitivities. Both configurations utilize the Pu ZPPR (Zero Power Physics Reactor) plates as fuel. The experiments were designed using a D-Optimality criteria, which is an optimization method minimizing the log-determinant of the adjusted nuclear data covariance for the target reactions. These experiments include not only inference of k eff , as done in all critical benchmark experiments, but several other responses as well, such as neutron multiplication measurements and reaction rate ratios. After analysis of the measured data is complete, adjustment of nuclear data will be performed to assess whether the new experimental data successfully reduced compensating errors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗

Implementation and testing of the GPT sensitivities in TRIPOLI-4

Computing the sensitivity of given reactor parameters with respect to nuclear data is mandatory in key applications such as the assessment of reactor safety or core design. In recent years, the development of Monte Carlo methods capable of estimating sensitivity coefficients based on the Standard or Generalized Perturbation Theory (SPT and GPT, respectively) has attracted intensive research efforts, especially in view of having reference golden standard simulations to be compared to faster but approximate deterministic calculations. Following our previous implementation and validation of the SPT functionalities in the Monte Carlo TRIPOLI-4, devoted to the sensitivity of the k eigenvalue, in this work we detail the development of GPT methods for ratios of reaction rates, and validate them against a few relevant benchmark results. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimizing group structures using hierarchical division

Creating group structures with few groups that give low errors is a difficult problem in reactor analysis. In recent years, automated optimization techniques have been applied to this task. We continue this trend by applying the hierarchical division algorithm to generate optimized group structures that minimize a cost function. At each stage, the algorithm adds a single group boundary to an existing group structure, dividing one group into two to increase the resolution of the group structure. The location of the added boundary is the one that gives the lowest error over all possible new boundary locations. Our implementation requires a beginning group structure, a set of candidate new boundary locations, and a set of reference reaction rates. As a proof of concept, we used WIMS-69 as the initial group structure, XMAS-172 as the ending group structures, and a 344-group reference group structure. Testing on two simple, homogenized reactor problems, we found that hierarchical division was able to reduce the error by a factor of around 5 with an increase of only 15% in the number of groups. Because hierarchical division can get stuck in local minima, it often reaches a plateau in its error reduction capability as many groups are added. Nevertheless, we find hierarchical division has strong potential to make good group structures into great group structures at a modest increase in computational cost. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Scalable algorithms for domain decomposition of Monte Carlo neutral particle transport simulations on unstructured mesh

Monte Carlo simulations are the gold standard in particle transport simulations, for shielding and nuclear reactor simulations Full core simulations can require terabytes of memory for tracking pin-resolved (up to 50k pins), nuclide (200 tracked), axial (20-200 zones) reaction rates Domain decomposition, splitting the calculation and the memory storage, among several nodes is challenging as it involves communication of particles across domain boundaries We created a new MOOSE application, MaCaw, to study this problem Algorithms developed for ray tracing in MOOSE were adapted for simulating neutral particle transport Physics and tally routines in OpenMC, dynamically linked, are called from the application

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Embedding Neural Thermal Scattering (NeTS) Modules in SERPENT for Higher Fidelity Advanced Reactor Analysis

When a neutron born in fission thermalizes to the order of $k$ $B$ $T$, it’s de-Broglie wavelength and energy approach the order of inter-atomic spacing and elementary lattice oscillations, respectively. $S$($a,β,t$) or the scattering law, uuantify these temperature-dependent crystallographic contributions to total cross section (or reaction rate). In a Monte Carlo analysis, cumulative distribution functions (CDFs) of $S$($a,β,t$) are loaded to memory from “A Compact ENDF” (ACE) files for stochastically selecting thermal scattered neutron trajectories. In this work, novel neural thermal scattering (NeTS) modules for $S$($a,β,t$) CDFs are designed, trained, serialized and embedded within SERPENT using Python’s limited C-API for on-the-fly deployment of crystalline graphite $S$($a,β,t$) sampling. Torchscript tracing and Numba just-in-time (JIT) compilation streamline neural inference on NVIDIA GPUs with CUDA libraries. Demonstrations of bare sphere thermalization of fast and thermal sources show excellent agreement between embedded NeTS in SERPENT and MCNP. With an explicit model of the reactor, NeTS can predict on-the-fly changes in TREAT neutron spectra as a function of local temperature, which can serve to improve transient and accident predictions in a multiphysics analysis framework. This framework can be further extended to account on-the-fly for changes in local graphitic microstructure to scattering cross sections, and outlines a novel coupling of modern machine learning with state-of-the-art reactor physics methods.

97 MATHEMATICS AND COMPUTING↗

Investigation of Ring-heterogeneous Geometry Approximation for Efficient and Practical SFR Analysis

A novel geometry approximation scheme named ring-heterogeneous approximation for intermediate resolution multiphysics analysis of SFRs is presented. This approach particularly targets neutronics coupled with thermo-mechanics feedback, which dominate the passive safety of SFR designs. The ring-heterogeneous approximation casts an array of pins into volume-preserving bands of materials that constitute the pins. It allows to efficiently simulate SFRs with an unstructured mesh and also to explicitly consider the deformation of each material by retaining their heterogeneity. Verification results with single assembly and core problems of the ABTR benchmark demonstrated that the ring-heterogeneous approximation is a reasonable approximation for SFRs. For a 2D full core problem, the ring-heterogeneous approximation significantly reduced the number of elements by at least a factor of 70 from a fully-heterogeneous mesh while almost exactly preserving the reference fully-heterogeneous solutions with the RMS error of only 0.02% in the ring-wise fission reaction rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of SPND Material Candidates for Real-Time Neutron Spectrum Monitoring Poster

The future of the nuclear industry will be shaped by advanced reactor designs requiring innovative instrumentation for safe operation. A key area of innovation is the ability to accurately and reliably monitor the neutron energy spectrum during operations, which would also enhance the experimental capabilities of existing research reactors. To address this need, our research focuses on the potential use of Self-Powered Neutron Detectors (SPNDs) for real-time neutron spectrum monitoring. We investigated various candidate materials for SPNDs by evaluating their energy-dependent reaction rates. Promising materials were further modeled using MCNP to determine the free electron production per neutron across the neutron energy spectrum, approximating their energy-dependent signal sensitivity. Our findings identified Samarium (Sm) and Gadolinium (Gd) as top candidates for thermal neutrons, Hafnium (Hf) and Thulium (Tm) for epithermal neutrons, and Tantalum (Ta) for fast neutrons. These materials will undergo future experimental validation to confirm the modeled responses and assess their feasibility for real-time neutron spectrum measurements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

The CO2 mineralization potential of New Mexico basalts and associated critical element behavior and mobility.

This presentation reports on laboratory experiments assessing the carbon dioxide (CO₂) mineralization potential of New Mexico basalts, including trachybasalt and basalt lithologies. Using simulated groundwater and controlled conditions, the study evaluates the kinetics of CO₂ sequestration, mineral saturation indices, and secondary mineral formation. It also investigates the mobility of critical and trace elements such as Ni, Zn, Cu, and Li during basalt dissolution. Results indicate that New Mexico basalts can effectively sequester CO₂ through carbonate formation, with grain size and mineral composition influencing reaction rates and elemental behavior. These findings support the viability of geologic carbon storage in mafic volcanic terrains and inform strategies for critical element recovery.

36 MATERIALS SCIENCE↗

High-fidelity Pebble Bed Reactor Depletion Based on Pebble Tracking Transport in Griffin

The pebble tracking transport (PTT) method is a high-fidelity, heterogeneous deterministic transport technique for pebble bed reactor analysis. It discretizes the broad-group neutron transport equation in space and angle with the discontinuous finite element and the discrete ordinates method, and utilizes various solving techniques, including mesh sweeping and diffusion acceleration, to provide pebble- wise reaction rates. This work presents the extension of the PTT method to enable fuel depletion capability in the Griffin code. We discuss the implementation details of the PTT-based high-fidelity depletion where isotope inventory of all individual pebbles is tracked through pre-determined pebble flow paths in the core. The implementation is verified with a generic pebble bed reactor model. Some preliminary equilibrium core results are included. Future works are also discussed.

97 - MATHEMATICS AND COMPUTING↗

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock tube/laser absorption measurement of the rate constant of the reaction: H 2 O 2 + CO 2 2OH + CO 2

We address the role of the linear mixing rule in the kinetics of the H 2 O 2 decomposition system by reporting the rate constant for H 2 O 2 + M = 2OH + M (M = Ar and CO 2 ) in the temperature range of 1087-1234 K at low pressures in a mixture of 20% CO 2 in Argon. The reaction rate constant was inferred from H 2 O concentrations monitored by using a laser-absorption spectroscopy-based water diagnostic. To the best of our knowledge, this is the first measurement of the rate constant of this reaction in a mixture to be reported in literature. A significant discrepancy was found between the rate constants derived using the traditional linear mixing rule and the reduced pressure linear mixing rule. Furthermore, this discrepancy can have serious implications on the predictive accuracy of these kinetic models, especially under conditions relevant to the operation of supercritical CO 2 (sCO 2 ) power cycles that rely on oxy-fuel combustion in a working fluid comprised almost entirely of CO 2 .

20 FOSSIL-FUELED POWER PLANTS↗

Methods and compositions to increase the rate of ligation reactions catalyzed by a sortase

In various embodiments methods and compositions to increase the rate and/or activity of ligation reactions catalyzed by a sortase. In certain embodiments reagents are provided that comprises, inter alia, a polypeptide comprising an amino terminal polyglycine sequence comprising at least three contiguous Gly residues (e.g., a triglycine sequence, a tetraglycine sequence, a pentaglycine sequence, etc.) followed by a peptide that is to be ligated to a moiety followed by a sequence comprising the catalytic domain of a Sortase A enzyme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Particle-Scale Reduction Analysis of CuFeMnO 4 with Hydrogen for Chemical Looping Combustion

In this work CuFeMnO 4 (copper iron manganese oxide) oxygen carrier was characterized using DSC/TGA, in-situ XRD, and SEM-EDS to gain a clear elucidation of the chemical looping combustion (CLC) reactions with H2 which is a component of synthesis gas. A reaction model which best described the experimental TGA reaction data was selected from the analysis. Intrinsic rate of reaction parameters obtained from the study could be used for designing a reactor for scale-up. A model with combined contributions from nucleation and growth model and dimensional phase boundary model was found to exhibit the best correlation with the experimental data. The analysis provided validation of intrinsic reaction rates and other rate parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anti-Arrhenius behavior of electron transfer reactions in molecular dimers

Rates of chemical reactions typically accelerate as the temperature rises, following the Arrhenius law. However, electron transfer reactions may exhibit weak temperature dependence or counterintuitive behavior, known as anti-Arrhenius behavior, wherein reaction rates decrease as temperature increases. Solvent reorganization energy and torsion-induced changes in electronic couplings could contribute to this unusual behavior, but how each contributes to the overall temperature dependence is unclear. One can decelerate the charge recombination process in photogenerated radical pairs or charge-separated states by harnessing this often-overlooked phenomenon. This means that we could achieve long-lived radical pairs without relying on conventional cooling. Using a series of homo molecular dimers, we showed that the degree of torsional hindrance dictates temperature-dependent torsion-induced changes in electronic coupling and, therefore, charge recombination rates. The overall temperature dependence is controlled by how changes in electronic coupling and the temperature-dependent solvent reorganization energy contribute to the rates of charge recombination. Our findings pave the way for rationally designing molecules that exhibit anti-Arrhenius behavior to slow down charge recombination, opening possibilities for applications in energy-related and quantum information technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗