Search NASASearch

SEARCH · Search NASA

Results for “Kinetic Monte Carlo”

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 19 records

Kinetic Monte Carlo Framework for Coupled Degradation and Dehydration of Anion Exchange Membranes

Kinetic Monte Carlo (kMC) simulations, augmented with temporal-acceleration schemes, can efficiently handle stiff reaction-transport networks when fast processes rapidly relax to quasi-equilibrium on a fixed lattice. However, in glassy anion-exchange membranes (AEM), rare and irreversible chemical degradation events continuously reshape the nanoscale morphology, and the associated hydration and transport degrees of freedom remain far from a well-defined local equilibrium. This combination of evolving state space and nonequilibrated fast dynamics lies outside the scope of existing kMC acceleration frameworks. Here, to address this challenge, we introduce an auxiliary-particle kinetic Monte Carlo (AP-kMC) scheme. In AP-kMC, short-lived mobile particles spawned at degradation sites execute hop, water-elimination, and decay moves, enforcing rapid local relaxation of the hydration structure while preserving the stochastic rules of kMC. Parameterized with molecular-dynamics morphologies and experimental solution degradation kinetics, AP-kMC reproduces the evolution of ion-exchange capacity, water uptake, and conductivity, and reveals a feedback loop in which poorly hydrated sites degrade first and each degradation event induces further local dehydration. The resulting thinning and fragmentation of water channels cause loss of hydrophilic percolation and abrupt conductivity collapse well before complete charge loss. AP-kMC thus reframes AEM durability as a coupled degradation–drying–percolation problem and provides a transferable strategy to simulate reactive, out-of-equilibrium polymer electrolytes where local solvation controls reactivity.

organic

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization

Development of Accelerated Kinetic Monte Carlo Code for Simulation of Helium Bubble Evolution

A mesoscale model to predict helium bubble evolution is needed for tritium applications. Such a model requires that the conventional kinetic Monte Carlo (kMC) simulations be significantly accelerated. The objective of this report is to (a) highlight the concepts and mathematical expressions of the accelerated method for defect implementation that have not been published, (b) show an example input file to run the kMC code, and (c) provide suggestions on future improvement following my retirement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

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

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Ytterbium disilicate is of interest as a potential environmental barrier coating for aerospace applications, notably for use in next generation jet turbine engines. In such applications, the transport of oxygen and water vapor through these coatings to the ceramic substrate is undesirable if high temperature oxidation is to be avoided. In an effort to understand the diffusion process in these materials, we have performed kinetic Monte Carlo simulations of vacancy-mediated and interstitial oxygen diffusion in Ytterbium disilicate. Oxygen vacancy and interstitial site energies, vacancy and interstitial formation energies, and migration barrier energies were computed using Density Functional Theory. We have found that, in the case of vacancy-mediated diffusion, many potential diffusion paths involve large barrier energies, but some paths have barrier energies smaller than one electron volt. However, computed vacancy formation energies suggest that the intrinsic vacancy concentration is small. In the case of interstitial diffusion, migration barrier energies are typically around one electron volt, but the interstitial defect formation energies are positive, with the result that the disilicate is unlikely to exhibit experience significant oxygen permeability except at very high temperature.

Coatings

Kinetic Monte Carlo Simulations of Oxygen Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Composite (CMC) materials are of interest for use in next-generation turbine engine components, offering a number of significant advantages, including reduced weight and high operating temperatures. However, in the hot environment in which such components operate, the presence of water vapor can lead to corrosion and recession, limiting the useful life of the components. Such degradation can be reduced through the use of Environmental Barrier Coatings (EBCs) that limit the amount of oxygen and water vapor reaching the component. Candidate EBC materials include Yttrium and Ytterbium silicates. In this work we present results of kinetic Monte Carlo (kMC) simulations of oxygen diffusion, via the vacancy mechanism, in Yttrium and Ytterbium disilicates, along with a brief discussion of interstitial diffusion.

Good, Brian S.

Kinetic Monte Carlo Simulations of Oxygen Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Composite (CMC) materials are of interest for use in next-generation turbine engine components, offering a number of significant advantages, including reduced weight and high operating temperatures. However, in the hot environment in which such components operate, the presence of water vapor can lead to corrosion and recession, limiting the useful life of the components. Such degradation can be reduced through the use of Environmental Barrier Coatings (EBCs) that limit the amount of oxygen and water vapor reaching the component. Candidate EBC materials include Yttrium and Ytterbium silicates. In this work we present results of kinetic Monte Carlo (kMC) simulations of oxygen diffusion, via the vacancy mechanism, in Yttrium and Ytterbium disilicates, along with a brief discussion of interstitial diffusion. An EBC system typically includes a bond coat located between the EBC and the component surface. Bond coat materials are generally chosen for properties other than low oxygen diffusivity, but low oxygen diffusivity is nevertheless a desirable characteristic, as the bond coat could provide some additional component protection, particularly in the case where cracks in the coating system provide a direct path from the environment to the bond coat interface. We have therefore performed similar kMC simulations of oxygen diffusion in this material.

Good, Brian S.

Kinetic Monte Carlo Simulations of Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Components (CMC) components for use in turbine engines offer a number of advantages compared with current practice. However, such components are subject to degradation through a variety of mechanisms. In particular, in the hot environment inside a turbine in operation a considerable amount of water vapor is present, and this can lead to corrosion and recession. Environmental Barrier Coating (EBC) systems that limit the amount of oxygen and water reaching the component are required to reduce this degradation and extend component life. A number of silicate-based materials are under consideration for use in such coating systems, including Yttterbium and Yttrium di- and monosilicates. In this work, we present results of kinetic Monte Carlo computer simulations of oxygen diffusion in Yttrium disilicate, and compare with previous work on Yttterbium disilicate. Coatings may also exhibit cracking, and the cracks can provide a direct path for oxygen to reach the component. There is typically a bond coat between the coating and component surface, but the bond coat material is generally chosen for properties other than low oxygen diffusivity. Nevertheless, the degree to which the bond coat can inhibit oxygen diffusion is of interest, as it may form the final defense against oxygen impingement on the component. We have therefore performed similar simulations of oxygen diffusion through HfSiO4, a proposed bond coat material.

Computer Simulation

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Silicon-based ceramic components for next-generation jet turbine engines offer potential weight savings, as well as higher operating temperatures, both of which lead to increased efficiency and lower fuel costs. Silicon carbide (SiC), in particular, offers low density, good strength at high temperatures, and good oxidation resistance in dry air. However, reaction of SiC with high-temperature water vapor, as found in the hot section of jet turbine engines in operation, can cause rapid surface recession, which limits the lifetime of such components. Environmental Barrier Coatings (EBCs) are therefore needed if long component lifetime is to be achieved. Rare earth silicates such as Yb2Si2O7 and Yb2SiO5 have been proposed for such applications; in an effort to better understand diffusion in such materials, we have performed kinetic Monte Carlo (kMC) simulations of oxygen diffusion in Ytterbium disilicate, Yb2- Si2O7. The diffusive process is assumed to take place via the thermally activated hopping of oxygen atoms among oxygen vacancy sites or among interstitial sites. Migration barrier energies are computed using density functional theory (DFT).

Computer Simulation

Kinetic Monte Carlo Investigation of the Effects of Vacancy Pairing on Oxygen Diffusivity in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia s high oxygen diffusivity and corresponding high ionic conductivity, and its structural stability over a broad range of temperatures, have made the material of interest for use in a number of applications, for example, as solid electrolytes in fuel cells. At low concentrations, the stabilizing yttria also serves to increase the oxygen diffusivity through the presence of corresponding oxygen vacancies, needed to maintain charge neutrality. At higher yttria concentration, however, diffusivity is impeded by the larger number of relatively high energy migration barriers associated with yttrium cations. In addition, there is evidence that oxygen vacancies preferentially occupy nearest-neighbor sites around either dopant or Zr cations, further affecting vacancy diffusion. We present the results of ab initio calculations that indicate that it is energetically favorable for oxygen vacancies to occupy nearest-neighbor sites adjacent to Y ions, and that the presence of vacancies near either species of cation lowers the migration barriers. Kinetic Monte Carlo results from simulations incorporating this effect are presented and compared with results from simulations in which the effect is not present.

Good, Brian S.

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH