Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic simulations”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Self-Learning Kinetic Monte Carlo Simulations of Radiation Damage in Nuclear Fuels

Understanding how irradiation affects the thermo-physical and mechanical properties of nuclear materials, such as thermal conductivity degradation in fuels and embrittlement of structural components, is critical to the safety and efficiency of nuclear reactors. These effects are largely governed by the formation and evolution of atomic-scale point defects and defect clusters. Due to their small sizes, however, these defects are invisible under high-resolution scanning transmission electron microscopy. This project aims to fill this experimental knowledge gap by integrating density functional theory (DFT), machine learning interatomic potential (MLIP), and kinetic Monte Carlo (KMC) techniques to predict longtime evolution of irradiation-induced defects in nuclear fuels.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal and dynamic partition of dumbbell interstitials in complex concentrated alloys

Complex concentrated alloys (CCAs) are promising candidates for applications in extreme conditions, such as irradiation where interstitial mediated diffusion is important. In CCAs with N principal elements, $\frac{N(N+1)}{2}$ types of dumbbell interstitials exist. Currently, there is no way to predict the thermal partition (fractional concentration at equilibrium) and the dynamic partition (fractional time an interstitial spends during diffusion) of each type of dumbbell interstitial. To mitigate this issue, this work proposes a theoretical model for computing the equilibrium concentrations and thermal partition of dumbbell interstitials in CCAs and validates the model using grand canonical Monte Carlo simulations. Lattice kinetic Monte Carlo simulations show that the thermal partition is equivalent to the dynamic partition, and both are governed by composition and formation energies of dumbbells. As a result, the model proposed provides a foundation for understanding radiation enhanced diffusion and induced segregation in CCAs under irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tailoring Hierarchical Structure and Rare Earth Affinity of Compositionally Identical Polymers via Sequence Control

Macromolecule sequence, structure, and function are inherently intertwined. While well-established relationships exist in proteins, they are more challenging to define for synthetic polymer nanoparticles due to their molecular weight, sequence, and conformational dispersities. Furthermore, to explore the impact of sequence on nanoparticle structure, we synthesized a set of 16 compositionally identical, sequence-controlled polymers with distinct monomer patterning of dimethyl acrylamide and a bioinspired, structure-driving di(phenylalanine) acrylamide (FF). Sequence control was achieved through multiblock polymerizations, yielding unique ensembles of polymer sequences which were simulated by kinetic Monte Carlo simulations. Systematic analysis of the global (tertiary- and quaternary-like) structure in this amphiphilic copolymer series revealed the effect of multiple sequence descriptors: the number of domains, the hydropathy of terminal domains, and the patchiness (density) of FF within a domain, each of which impacted both chain collapse and the distribution of single- and multichain assemblies. Furthermore, both the conformational freedom of chain segments and local-scale, β-sheet-like interactions were sensitive to the patchiness of FF. To connect sequence, structure, and target function, we evaluated an additional series of nine sequence-controlled copolymers as sequestrants for rare earth elements (REEs) by incorporating a functional acrylic acid monomer into select polymer scaffolds. We identified key sequence variables that influence the binding affinity, capacity, and selectivity of the polymers for REEs. Collectively, these results highlight the potential of and boundaries of sequence control via multiblock polymerizations to drive primary sequence ensembles hierarchical structures, and ultimately the functionality of compositionally identical polymeric materials.

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

Anomalous phase separation in a correlated electron system: Machine-learning–enabled large-scale kinetic Monte Carlo simulations

Significance Phase separation is crucial to the functionalities of many correlated electron materials with notable examples including colossal magnetoresistance in manganites and high- T c superconductivity in cuprates. However, the nonequilibrium phase-separation dynamics in such systems are poorly understood theoretically, partly because the required multiscale modeling is computationally very demanding. With the aid of machine-learning methods, we have achieved large-scale dynamical simulations in a representative correlated electron system. We observe an unusual relaxation process that is beyond the framework of classical phase-ordering theories. We also uncover a correlation-induced freezing behavior, which could be a generic feature of phase separation in correlated electron systems.

Zhang, Sheng↗
DOE OSTI2022

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2021

Assessment of a Detailed Biomass Pyrolysis Kinetic Scheme in Multiscale Simulations of a Single-Particle Pyrolyzer and a Pilot-Scale Entrained Flow Pyrolyzer

A detailed biomass pyrolysis kinetic scheme was assessed in the multiscale simulations of a single-particle pyrolyzer with slow pyrolysis and a pilot-scale entrained flow pyrolyzer with fast pyrolysis. The detailed kinetic scheme of biomass pyrolysis developed by the CRECK group consists of 32 reactions and 58 species. A multiscale simulation model was developed, where the CRECK kinetics was employed to simulate biomass pyrolysis reactions, a one-dimensional particle model was utilized to simulate the intraparticle transport phenomena, and the particle-in-cell (PIC) model was employed to simulate the hydrodynamics. The multiscale model was first applied to simulate a single-particle pyrolysis experiment. The simulation with nonisothermal particles matched the experimental data better than the simulation with isothermal particles. Then the multiscale model was applied to simulate the pilot-scale entrained flow pyrolyzer. In this case, the simulation with isothermal particles matched the experimental data better than the simulation with nonisothermal particles. The reason for this difference might be that the kinetics itself already partially included the intraparticle transport effect as it was fitted using both TGA data (slow pyrolysis of small size biomass) and fluidized bed data (fast pyrolysis of relatively large size biomass). This study provides some insights into biomass pyrolysis kinetics development and pyrolyzer multiscale simulation for a future study.

09 BIOMASS FUELS↗
DOE OSTI2025

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

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

Global gyrokinetic simulations of electrostatic microturbulent transport using kinetic electrons in LHD stellarator

Global gyrokinetic simulations of ion temperature gradient (ITG) and trapped electron mode (TEM) in the LHD stellarator are carried out using the gyrokinetic toroidal code (GTC) with kinetic electrons. ITG simulations show that kinetic electron effects increase the growth rate by more than 50% and more than double the turbulent transport levels compared with simulations using adiabatic electrons. Zonal flow dominates the saturation mechanism in the ITG turbulence. Nonlinear simulations of the TEM turbulence show that the main saturation mechanism is not the zonal flow but the inverse cascade of high to low toroidal harmonics. Further nonlinear simulations with various pressure profiles indicate that the ITG turbulence is more effective in driving heat conductivity whereas the TEM turbulence is more effective for particle diffusivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Catalytic resonance theory for the kinetic signatures of multiple wavelength photocatalysis

Multi-step surface chemistry promoted with one or two wavelengths of light exhibited distinct kinetic signatures indicative of the light-adsorbate interactions and the number of photon-sensitive elementary steps. In this work, kinetic Monte Carlo simulations identified the kinetic response of general surface mechanisms to variation in the per-site photon flux of one or two wavelengths of incident light. Photocatalytic rates were described via a non-dimensional photon flux, identifying multiple kinetic regimes with unique degrees of rate control. Under photon-controlled kinetic conditions, maximum quantum yield and turnover frequency were obtained under constant illumination, while pulsing of one or more light sources was shown to exhibit slower rates and less efficient photocatalytic promotion due to the inherent dynamic nature of light, which comprises a stream of photons. Simulations provided distinct kinetic signatures in catalytic rates and Arrhenius kinetics for specific light-surface-adsorbate interactions corresponding to photocatalytic promotion of specific steps in surface chemistry.

catalytic resonance↗