Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Michel Electron Selection with SPINE for DUNE Far Detector Simulation

Michel electrons are a valuable input for particle detector calibration due to their consistent kinetic energy distribution. This report details the evaluation of a Michel electron identification method's application to simulated data from the DUNE (Deep Underground Neutrino Experiment) far detector. This method, which relies on the neural network-based particle classification software SPINE (Scalable Particle Imaging with Neural Embeddings), was developed and calibrated using simulated data for the SBND (Short-Baseline Neutrino Detector) experiment before being applied to simulated DUNE data from a 1x2x6 subset of far detector modules.

Wilson, Dante [Colorado State U.]↗
DOE OSTI2024

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Transient Catalytic Reaction Analysis Through Signal Defragmentation

The Temporal Analysis of Products (TAP) pulse response technique provides valuable insights into catalytic function and reaction kinetics. However, complex fragmentation patterns in the TAP mass spectrometry signals can complicate precise quantification, particularly when analyzing transient gas flux data typical of TAP experiments. This work demonstrates a standard defragmentation method that deconvolves transient TAP signals while maintaining the temporal resolution of the experiment. First, the integrals of calibration gas fluxes are used to determine the fingerprint fragmentation pattern and construct a fragmentation matrix. This matrix is then used to defragment experimental flux data at each recorded time point via a non-negative least squares regression. The effectiveness of this method is demonstrated using virtual data and control experiments with a TAP reactor system. The defragmentation is then applied to the more complex propane dehydrogenation reaction on a chromia/alumina catalyst, which can contain up to ten significant gas species in the reactor outlet. Initial propane pulsing reveals an induction period during which propane is fully oxidized to CO2, followed by partial reduction to CO. Afterwards, there is a transition in chemistries towards coking and propylene production. Our example illustrates a practical method for the accurate determination of the time-dependent reactant/product concentrations and rates for a thorough analysis of the propane dehydrogenation kinetics. This approach can be broadly applied to any transient mass spectrometry experiment for a better understanding of catalyst-reaction dynamics.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Pauli potential formalism at finite temperature

At zero temperature, the Pauli potential—the functional derivative of the Pauli kinetic energy density functional—is the key to the accuracy of the orbital-free density functional theory (OFDFT) as it is supposed to capture all the effects associated with the Pauli exclusion principle. Here, we extend this concept to finite temperature by defining Pauli free energy and the modified Pauli free energy, both representing the natural generalizations of the Pauli term from zero- T to finite- T . We discuss their physical interpretation, the mathematical nuances, and the applicability, arguing that the modified Pauli potential should be used as an extension of the zero- T counterpart within the OF-DFT framework. Through analytical and numerical methods, we then analyze some of the exact properties concerning the modified Pauli free- and kinetic-energy terms and examine the temperature dependence of the modified Pauli potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electrochemically Enhanced Carbonate Precipitation into Building Materials: A Scalable Carbon Sequestration Strategy

Decarbonization goals across hard-to-abate industries have prompted an urgent need for advanced carbon capture and storage technologies. Sequestering CO2 into carbonate minerals is a scalable method of carbon management with the ability to produce value-added carbon negative materials from waste streams for the construction industry. Waste streams rich in Ca and Mg such as nickel mine tailings, iron/steel slag, and reverse osmosis brines can store 7.6 Mt CO2/year as minerals. Additionally, CO2 mineralization in acid-neutralization processes currently present in industrial waste treatment can eliminate associated CO2 emissions of lime processing by improving process circularity. The carbonate minerals formed from these waste sources are valuable as components of carbon-negative concrete, which have the potential to sequester 1.8 billion Mt of CO2/year. Electrochemical means of CO2 mineralization improves the kinetics of the thermodynamically favorable mineralization process, lessening or eliminating the high energy requirements of traditional methods. Here, we investigate benchtop scale electrochemical CO2 mineralization of alkaline mining waste, highlighting the effects of key constituents in mining waste on the mineralization process.

carbon capture↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2025

Online Bayesian State Estimation for Real-Time Monitoring of Growth Kinetics in Thin Film Synthesis

Rapid validation of newly predicted materials through autonomous synthesis requires real-time adaptive control methods that exploit physics knowledge, a capability that is lacking in most systems. Here, in this study, we demonstrate an approach to enable real-time control of thin film synthesis by combining in situ optical diagnostics with a Bayesian state estimation method. We developed a physical model for film growth and applied the direct filter (DF) method for real-time estimation of nucleation and growth rates during pulsed laser deposition (PLD). We validated the approach using simulated and experimental reflectivity data for WSe 2 growth and ultimately deployed the algorithm on an autonomous PLD system during the growth of 1T'-MoTe 2 . The DF robustly estimates growth parameters in real time at early stages of growth, down to 15% monolayer area coverage. This fusion of in situ diagnostics, data assimilation, and physical modeling opens new opportunities in adaptive control of synthesis trajectories toward desired material states.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Laser-plasma ion beam booster based on hollow-channel magnetic vortex acceleration

Laser-driven ion acceleration provides ultrashort, high-charge, low-emittance beams, which are desirable for a wide range of high-impact applications. Yet after decades of research, a significant increase in maximum ion energy is still needed. This paper introduces a quality-preserving staging concept for ultraintense ion bunches that is seamlessly applicable from the nonrelativistic plasma source to the relativistic regime. Full three-dimensional particle-in-cell simulations prove robustness and capture of a high-charge proton bunch, suitable for readily available and near-term laser facilities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Zero-Order Reaction Kinetics v. 3.6

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.

McNenly, MatthewJ [Lawrence Livermore National Lab↗
DOE OSTI2024

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗
DOE OSTI2025

Catalytic water splitting kinetics and optimization

This project was motivated primarily by the opportunity of turning photocatalysis into a cost efficient hydrogen generation method for domestic energy and chemical production, and by the need to provide detailed mechanisms of photocatalytic reactions of broad scientific and technological interest. Specifically, we focused on the development of a new technique for analysis of (photo)catalytic reaction mechanisms to guide photocatalyst design. Our approach was based on applying velocity map imaging (VMI) to surface catalyzed reactions in order to identify photocatalytic reaction products, intermediates, and relate these to photocatalyst composition, dimensions, size and other characteristics. By coupling VMI with other spectroscopic and microstructural characterization, we aimed to reveal the features that limit performance of photocatalysis. As a first step to achieve this goal, we constructed a custom VMI microscope specifically designed for investigation of surface photocatalytic reactions and validated it’s design using a photodissociatinon of MnBr(CO) 5 adsorbed on a metallic surface. Using both 230 nm and 266 nm pulsed laser excitations, we observed that there are likely three main pathways for photodissociation of Mn(CO)5Br: from the gas phase, the surface, and ejected Mn complex from the surface increasing the gas phase concentration.

08 HYDROGEN↗
DOE OSTI2026

High temperature stabilization of ultrafine grain tungsten alloys through synergistic compositional complexities

Thermally-stabilized fine-grained microstructures in tungsten provide a pathway to harnessing enhanced properties such as a reduced ductile-to-brittle transition temperature and improved strength while mitigating the adverse effects of grain growth and recrystallization. Here, in this study, we employ a material design strategy that relies on grain boundary segregation in the nanocrystalline state driven by alloy thermodynamics balanced with impurity scavenging through in situ formation of kinetically-stabilizing metal carbides. A W-Ti-Cr alloy is designed through lattice Monte Carlo methods and subsequently synthesized in a single-phase nanocrystalline state, which upon annealing, evolves into an ultrafine grained microstructure containing chromium grain boundary segregation collectively with a dispersed TiX (X=C,O) phase through the reaction of titanium with carbon and oxygen impurities. In situ synchrotron X-ray diffraction experiments demonstrate that increased Ti promotes stabilization across a larger temperature range but with diminishing returns above 10 at.% Ti. Long-term stability was confirmed through the retention of the ultrafine grained microstructure for a total of 8 days (192 h) at 1300 °C without grain/carbide growth and/or recrystallization. Our results demonstrate that, through strategic tailoring of composition and microstructure, one can harness the benefits of both thermodynamic and kinetic stabilization mechanisms, opening pathways for future alloy formulations that expand the window of stability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗
DOE OSTI2026

Delayed Neutron Precursor Group Parameter and Spectra Generation from Fast Fission of 235U in SCALE

Delayed neutron precursor (DNP) group data are important for modeling reactor dynamics. Although the data for individual DNPs have been developed over time, the DNP group data present in the Evaluated Nuclear Data Files (ENDF) have not been updated in the past 20 years. In this work, we use SCALE to recreate the Godiva experiment that was used to generate the original DNP group structure for fast fission of 235U. However, each DNP is modeled using up-to-date data, and the results are then converted into a newly updated group structure. This conversion uses an iterative linear least squares solver to minimize chi-squared. The approaches used in this work also enable energy spectrum generation and uncertainty tracking. The method used in this paper for fast 235U fission DNP group structure updating can be applied to different energies and fissile nuclides. Demonstration of the uncertainty tracking in reactor kinetics and dynamics simulations is shown using point reactor kinetics simulations. Results show that there are data discrepancies between the International Atomic Energy Agency database and data used in ORIGEN, which are currently being fixed. Results also show that the proposed method for group spectra generation performs well.

Seifert, Luke [University of Illinois]↗
DOE OSTI2025

Kinetic Investigation by Batch Processing of Film Growth during the Electrodeposition of Titanium from Ethaline

Titanium electrodeposition as a technological advance has been out of reach due to the poor film quality obtained when using previously developed methods. This study shows that the electrodeposition of Ti is possible and feasible from the deep eutectic solvent ethaline. Moreover, this study shows through batch analysis that this deposition has temperature-dependent aspects that alter the film deposition kinetics upon modest heating. This ultraviolet–visible spectroscopic study is highlighting the more facile reduction of Ti 4+ to Ti 3+ and then to Ti 2+ at elevated temperatures. Understanding the underlying kinetics of electrodeposition opens up new venues for the application of this important metal.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evolution and Degradation Patterns of Electrochemical Cells Based on the Analysis of Interfacial Phenomena at Li Metal Anode/Electrolyte Interfaces

In this work, we report the results of a theoretical–computational analysis of the solid electrolyte interphase (SEI) growth and degradation dynamics occurring in lithium metal batteries during cycling. We use ab initio-kinetic Monte Carlo simulations to generate a synthetic data set, which is analyzed by machine learning methods. We aim to determine: (i) how modifications in interfacial interaction energies between solid electrolyte interphase (SEI) blocks and between Li ions and SEI facets impact the Coulombic efficiency (CE) of the battery and (ii) what factors, including reactions, microscopic transport, and other interfacial events, may lead to cell performance “failure” during prolonged charge and discharge cycles, signaled as a sharp decay in the CE over cycling. The demonstration of our approach is done on a cell including a Li metal surface interfacing with a previously introduced state-of-the-art electrolyte, and the idea can be applied to any electrochemical system. Outcomes include the identification of the leading chemical, physical, and structural variables causing cell failure and relating them to the electrolyte formulation, thus paving the way to future more refined analysis and electrolyte design.

batteries↗