Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic model data”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

This work presents a semi-detailed kinetic model to address the pyrolytic gas-phase reactivity of volatiles formed during thermal degradation of polyethylene (PE). The model builds on a validated multi-step condensed-phase model and employs validated lumping approaches. Short-chain compounds are modelled with high detail, while long-chain ones are described by surrogate species representative of diesel-cuts (NC16H32) and waxes (NC30H60). The reactivity of short chains is described through the comprehensive CRECK kinetic model, updated to align C5-C7 olefins based on recent literature experimental data. Due to the lack of experimental data for longer olefins, their reactivity is modeled by analogy to the shorter ones, ensuring an asymptotic behavior with increasing carbon numbers. The semi-detailed model is validated through experimental data on PE pyrolysis, assuming an instantaneous mixing of the inert inlet flow with released volatiles, followed by a segregated plug-flow behavior. Validation across different reactor setups confirms the model’s capability to predict detailed product distributions. Despite minor discrepancies, the proposed model effectively captures experimental trends. Further work will address modelling the reactivity in oxygen-containing environments.

kinetics↗
DOE OSTI2024

Material analysis of metal hydrides for bulk hydrogen storage

The United States Department of Energy (DOE) Hydrogen Storage Engineering Center of Excellence (HSECoE), 2009–2014, was established to evaluate the viability of material-based hydrogen storage for light duty vehicles. The HSECoE concluded that there were no metal hydrides that simultaneously met the DOE targets for gravimetric and volumetric capacity, and rate of onboard refueling. However, it was noted by the DOE that potential applications existed for stationary hydrogen storage systems. An application of particular interest is using a metal-hydride to store a hydrogen supply for fuel cell backup power to a 5 MWe data center. A generalized numerical model that coupled mass and energy conservation, chemical kinetics and thermodynamics was developed to simulate the discharge process for a given metal-hydride. The model was applied to an interstitial and a complex metal-hydride. In conclusion, it was found that the interstitial metal-hydride, which is discharged by a pressure swing is preferrable in a practical sense to the complex metal-hydride, which requires a temperature swing discharge process.

08 HYDROGEN↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Solid-State Transformation of an Additive Manufactured Inconel 625 Alloy at 700 °C

Inconel 625, a nickel-based superalloy, has drawn much attention in the emerging field of additive manufacturing (AM) because of its excellent weldability and resistance to hot cracking. The extreme processing condition of AM often introduces enormous residual stress (hundreds of MPa to GPa) in the as-fabricated parts, which requires stress-relief heat treatment to remove or reduce the internal stresses. Typical residual stress heat treatment for AM Inconel 625, conducted at 800 °C or 870 °C, introduces a substantial precipitation of the δ phase, a deleterious intermetallic phase. In this work, we used synchrotron-based in situ scattering and diffraction methods and ex situ electron microscopy to investigate the solid-state transformation of an AM Inconel 625 at 700 °C. Our results show that while the δ phase still precipitates from the matrix at this temperature, its precipitation rate and size at a given time are both smaller when compared with their counterparts during typical heat treatment temperatures of 800 °C and 870 °C. A comparison with thermodynamic modeling predictions elucidates these experimental findings. Our work provides the rigorous microstructural kinetics data required to explore the feasibility of a promising lower-temperature stress-relief heat treatment for AM Inconel 625. The combined methodology is readily extendable to investigate the solid-state transformation of other AM alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The kinetics of SARS-CoV-2 infection based on a human challenge study

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Integrated modeling framework (FVCOM-ICM/Toxi) to simulate the fate and transport of polychlorinated biphenyls in urban estuaries–Case study for Puget Sound, WA

Puget Sound is an urban estuary that exhibits persistent polychlorinated biphenyls (PCBs) contamination despite years of remediation efforts. For robust management actions, determining the timing, location, and magnitude of PCB loading sources and transport pathways using field data alone is challenging due to complex water circulation and biogeochemical processes. Here, this study aims to develop an integrated modeling framework that couples complex estuarian circulation with biogeochemical processes and associated interactions with PCB kinetics. The model simulates PCB accumulation in the lower tropic food web, demonstrating PCB intrusion into primary producers and its biomagnification in pelagic consumers. The PCB data from a new field survey was used to calibrate/validate the new PCB modeling framework for Puget Sound. The software program of the modeling framework is available to the user community for the applications of toxic contaminants transport in marine waters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A mechanism-based quantitative multi-scale framework for investigating irradiation hardening of tungsten at low temperature

We present the development and validation of a mechanism-based multi-scale modeling framework to quantitatively link the irradiation defect evolution kinetics at the microscopic scale to the macroscopic yield strength and flow stress evolution of tungsten irradiated at low temperature (0.08 < T/T m < 0.18, where T m is the melting point and is equal to 3673 K). The mechanism-based strength model, proposed as the superposition of thermal softening and modified dispersed barrier hardening, is developed to understand the underlying strengthening mechanisms. The thermal softening exponent as defined in the Johnson-Cook model is obtained by fitting unirradiated yield strength at different temperatures. A set of irradiation-induced defect kinetics equations with the thermodynamics parameters derived from the atomic calculations reported in the literature is used to determine the densities and sizes of defect clusters at the meso-scale, and the predicted irradiation damage characteristics are in reasonable quantitative agreement with experimental data from literature. The effects of irradiation condition (temperature and irradiation dose) and test temperature on yield strength are quantitatively predicted and compared with experimental measurements from the literature. Here, the predicted irradiation hardening and strain hardening are compared with experimental data as model validation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Decoupling electrode kinetics to elucidate reaction mechanisms in alkaline water electrolysis

Alkaline water electrolysis (AWE) presents key advantages, including reduced material costs, enhanced operational stability, and compatibility with non-precious metal catalysts, positioning it as a scalable route for hydrogen production. In this study, we introduce a minimally invasive single-cell configuration incorporating a reference electrode via diaphragm extension to form an internal ion channel. This setup, combined with an interfaced potentiostat and auxiliary electrometer, enables real-time, independent monitoring of anode and cathode behavior, offering high-resolution electrochemical diagnostics. While it is well established that the hydrogen evolution reaction (HER) exhibits sluggish kinetics in alkaline media, our study reveals that this limitation persists even in practical AWE systems where nickel-based substrates are used as electrodes. This observation is supported by both experimental data and voltage breakdown modeling. Arrhenius-type analysis reveals that localized electric fields induced by catalysts shift the reaction kinetics from classical Butler–Volmer behavior toward a Marcus-like regime, where interfacial molecular dynamics and bimolecular charge transfer dominate. We propose a semi-empirical model and a surficial reaction mechanism to describe these dynamics. This work underscores the critical need for cathode innovation and provides a rational framework for designing advanced catalysts and electrode architectures to optimize AWE performance.

08 HYDROGEN↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2021

Incorporating concentration-dependent sediment microbial activity into methylmercury production kinetics modeling

We report in anoxic environments, anaerobic microorganisms carrying the hgcAB gene cluster can mediate the transformation of inorganic mercury (Hg(II)) to monomethylmercury (MMHg). The kinetics of Hg(II) transformation to MMHg in periphyton from East Fork Poplar Creek (EFPC) in Oak Ridge, TN have previously been modeled using a transient availability model (TAM). The TAM for Hg(II) methylation combines methylation/demethylation kinetics with kinetic expressions for processes that decrease Hg(II) and MMHg availability for methylation and demethylation (multisite sorption of Hg(II) and MMHg, Hg(II) reduction/Hg(0) oxidation). In this study, the TAM is used for the first time to describe MMHg production in sediment. We assessed MMHg production in sediment microcosms using two different sediment types from EFPC: a relatively anoxic, carbon-rich sediment with higher microbial activity (higher CO 2 production from sediment) and a relatively oxic, sandy, carbon-poor sediment with lower microbial activity (lower CO 2 production from sediment). Based on 16s rRNA sequencing, the overall microbial community structure in the two sediments was retained during the incubations. However, the hgcA containing methanogenic Euryarchaeota communities differed between sediment types and their growth followed different trajectories over the course of incubations, potentially contributing to the distinct patterns of MMHg production observed. The general TAM paradigm performed well in describing MMHg production in the sediments. However, the MMHg production and ancillary data suggested the need to revise the model structure to incorporate terms for concentration-dependent microbial activity over the course of the incubations. We modified the TAM to include Monod-type kinetics for methylation and demethylation and observed an improved fit for the carbon-rich, microbially active sediment. Overall our work shows that the TAM can be applied to describe Hg(II) methylation in sediments and that including expressions accounting for concentration-dependent microbial activity can improve the accuracy of the model description of the data in some cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Physics-Informed Neural Network Solution of Point Kinetics Equations for a Nuclear Reactor Digital Twin

A digital twin (DT) for nuclear reactor monitoring can be implemented using either a differential equations-based physics model or a data-driven machine learning model. The challenge of a physics-model-based DT consists of achieving sufficient model fidelity to represent a complex experimental system, whereas the challenge of a data-driven DT consists of extensive training requirements and a potential lack of predictive ability. We investigate the performance of a hybrid approach, which is based on physics-informed neural networks (PINNs) that encode fundamental physical laws into the loss function of the neural network. We develop a PINN model to solve the point kinetic equations (PKEs), which are time-dependent, stiff, nonlinear, ordinary differential equations that constitute a nuclear reactor reduced-order model under the approximation of ignoring spatial dependence of the neutron flux. The PINN model solution of PKEs is developed to monitor the start-up transient of Purdue University Reactor Number One (PUR-1) using experimental parameters for the reactivity feedback schedule and the neutron source. The results demonstrate strong agreement between the PINN solution and finite difference numerical solution of PKEs. We investigate PINNs performance in both data interpolation and extrapolation. For the test cases considered, the extrapolation errors are comparable to those of interpolation predictions. Extrapolation accuracy decreases with increasing time interval.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

In-situ irradiation-induced studies of grain growth kinetics of nanocrystalline UO 2

The thermal conductivity of UO 2 fuel needs to be high enough to dissipate the heat generated from the fission reaction. Since grain size affects thermal conductivity and grain size can evolve with irradiation, it is critical to understand in-reactor UO 2 grain growth. Most studies of grain growth in UO 2 are based on thermally driven processes at elevated temperatures. However, studies have shown that grain growth can occur even at cryogenic temperatures by ballistic processes. Such irradiation-induced grain growth in UO 2 is yet to be studied. Advanced in-situ Kr ion irradiation and transmission electron microscopy were systematically performed on nanocrystalline UO 2 thin films at temperatures ranging from 50 K to 1073 K; grain growth was observed at all temperatures. A combination of manual and machine learning techniques was used to measure and plot grain size evolution against irradiation fluence at various irradiation temperatures. The machine learning method has significantly improved the analysis efficiency and reduced human labors. The grain diameter data were fitted using classical grain growth and thermal spike models to describe grain growth kinetics with and without irradiation effect. Grain growth during low temperature irradiation (≤ 475 K) can be well described by the thermal spike model. Above 475 K, there were additional thermally assisted processes that further accelerate the grain growth. At the highest irradiation temperature about 1075 K, both irradiation-induced dislocation loops and cavities/bubbles were observed to form in the UO 2 . In this report, the effects of irradiation-induced defects on grain growth kinetics are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

In Situ Synchrotron X-ray Scattering Investigation of Cathodic ZIF-8 Deposition on Graphite Using 3D-Printed Cells

Electrochemical approaches provide unique means to fabricate thin films of metal organic frameworks (MOFs). However, the kinetics of electrochemical MOF deposition has not been quantitatively investigated so far. In this study, we report the first in situ measurements of electrochemical MOF growth using transmission synchrotron X-ray scattering. Electrochemical cells based on poly(lactic acid) with two windows were fabricated using fused-deposition modelling. Here, the resulting 3D-printed cells, the surface of which was coated with paraffin wax to prevent solvent percolation through the polymer material, were used to monitor the cathodic growth of zeolitic imidazolate framework-8 (ZIF-8) on graphite in a methanol solution containing ZnCl 2 and 2-methylimidazole (Hmim) at different cathodic potentials. The resulting time-resolved X-ray diffraction data showed a gradual increase in crystal size with negligible changes in crystal orientation during the cathodic ZIF-8 deposition. More importantly, the time-resolved data provided a means to quantitatively assess the kinetics of the cathodic ZIF-8 growth using the Gualtieri model, revealing that the cathodic potential and Hmim concentration affected crystal growth kinetics, but not nucleation kinetics. These ZIF-8 samples exhibited changes in X-ray diffraction pattern after being washed with methanol and dried in air, indicating that in situ measurements are essential to investigate mechanisms behind MOF electrodeposition.

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

General H2O Outgassing Model Applicable to Silica-Filled Silicones

Moisture-induced agining and degradation are continuing concerns in long-term storage and use of silicone materials. The silica fillers frequently used to provide mechanical reinforcement in silicone formulations have a high affinity for H2O and act as a source of moisture outgassing later when placed in dry or vacuum applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-exchange Resin at Initial pH ~3 – 7 and Temperature 25 – 55 °C

This work presents systematic aqueous-phase adsorption equilibrium, kinetics, and column breakthrough measurements with three key biointermediates that are common compounds in many bioprocesses. Adsorption equilibrium experiments were carried out with acetic acid, butyric acid, and lactic acid on a commercial ion-exchange resin, Amberlite IRN-78, at wide ranges of acid concentration (8-500 mmol/L), initial pH (similar to 3-7), and temperature (25-55 degrees C), simulating the effluent characteristics from different fermenter operations. The kinetics and column breakthrough experiments were conducted at an initial pH of 6 and a concentration of 200 mmol/L. The equilibrium study shows a higher loading at the initial pH < pK(a) and a lower loading at the initial pH > pK(a). Overall removal varies between 16 and 99% depending on the initial pH, temperature, and organic acid concentration and type. The study further indicates monolayer adsorption at the equilibrium pH > 10 and multilayer adsorption at the equilibrium pH < 6. The thermodynamic modeling of adsorption isotherm data was carried out using Langmuir and Freundlich isotherms. IRN-78 presents fast adsorption kinetics as the maximum loading was attained in <= 10 min and nearly the same breakthrough time for all three organic acids involved in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗