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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2003

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.↗

NASA NTRS2013

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↗

NASA NTRS2011

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↗

NASA NTRS2015

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↗

NASA NTRS2017

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

Good, Brian S.↗

NASA NTRS2017

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.↗

NASA NTRS2017

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↗

NASA NTRS2015

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↗

NASA NTRS2011

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.↗

NASA NTRSDate not supplied

Atomistic Multi-Lattice Kinetic Monte Carlo (Kmc) Modeling of Hyperthermal Oxidation of Multi-Layer Graphene

Explore the source record for details and available documents.

Sharon Edward↗

NASA NTRS2010

Implementing the Second-Order Fermi Process in a Kinetic Monte-Carlo Simulation

Radio JOVE is an education and outreach project intended to give students and other interested individuals hands-on experience in learning radio astronomy. They can do this through building a radio telescope from a relatively inexpensive kit that includes the parts for a receiver and an antenna as well as software for a computer chart recorder emulator (Radio Skypipe) and other reference materials

Summerlin, Errol J.↗

NASA NTRSDate not supplied

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗

NASA NTRS2026

Modeling Study of Hatch Spacing’s Effect on Grain Morphology in Repairing Damaged Metal Parts with Welding and Hatch Spacing’s Potential Use in Lunar Exploration

The welding process is a potential way of repairing a damaged metal component, especially cavity damage caused by a harsh environment like the lunar environment, which is characterized by large temperature differences and reduced gravity. The adjustment of welding parameter (e.g., hatch spacing) can improve production efficiency in the repair process. Seen from the microstructural level, hatch spacing sensitivity affects the metallic grain evolution and morphology in the welding process, which can further influence a repaired part’s mechanical properties; however, the study of hatch spacing’s effect on microstructure is challenging. Traditional experimental procedures are costly and time-consuming, and any change in hatch spacing value needs roll-back of experimental procedure. A modeling study can address the above challenges in experimental observation. In this research, a modeling approach based on the Kinetic Monte Carlo (KMC) Potts theory was used to simulate grain evolution and morphology with three hatch spacings. Through quantifying and analyzing the predicted grain morphologies, the effect of hatch spacing on microstructure in a welding-fabricated part was investigated. The predicted grain morphologies were validated with an EBSD image of welding microstructure, which has been published before. The primary grain morphologies were columnar grains with a small amount of fine equiaxed grains formed in the scanning path centerline. When increasing the hatch spacing, the columnar grains become larger and more lengthy, while the effect of hatch spacing on the equiaxed grains is not obvious.

Welding for repairing↗

NASA NTRS1994

Transversely accelerated ions in the topside ionosphere

Data from the rocket campaigns Mechanism in the Auroral Region for Ion Energization (MARIE) and TOpside Probe of the Auroral Zone (TOPAZ) III, within regions of low-altitude transversely accelerated ions, are interpreted to explain the acceleration of the ions. Using the Monte Carlo kinetic technique to evaluate the ion heating produced by the simultaneously observed lower hybrid waves, we find that their observed electric field amplitudes are sufficient to explain the observed ion energies in the MARIE event. Much of the uncertainty in evaluating the efficiency of a plasma wave induced particle heating process which is dependent on a velocity resonance comes from the lack of information on the phase velocities of the waves. In the case of the MARIE observations, our modeling efforts show that features in the ion velocity distribution are consistent with the wave phase velocities inferred from interferometer measurements of wavelengths. The lower hybrid waves with which low-altitude transversely accelerated ions are associated are frequently observed to be concentrated in small-scale wave packets called 'spikelets'. We demonstrate through the scaling of the size of these wave packets that they are consistent with the theory of lower hybrid collapse. Using the Monte Carlo technique, we find that if the lower hybrid field energy is concentrated in these wave packets, it is still adequate to accelerate the ionospheric ions to the observed energies.

Retterer, John M.↗

NASA NTRS2005

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The paper discusses the following: 1. No-cost Extension. The no-cost extension is required to complete the work on the unified model codes (both hydrodynamic and kinetic Monte Carlo) as described in the initial proposal and previous annual reports. 2. Scientific Accomplishments during the Report Period. We completed a comprehensive model of Alfvtn wave reflection that spans the full distance from the photosphere to the distant heliosphere. 3. Comparison of Accomplishments with Proposed Goals. The proposal contained two specific objectives for Year 3: (1) to complete the unified model code, and (2) to apply it to various kinds of coronal holes (and polar plumes within coronal holes). Although the anticipated route toward these two final goals has changed (see accomplishments 2a and 2b above), they remain the major milestones for the extended period of performance. Accomplishments la and IC were necessary prerequisites for the derivation of "physically relevant transport and mode-coupling terms" for the unified model codes (as stated in the proposal Year 3 goals). We have fulfilled the proposed "core work" to study 4 general types of physical processes; in previous years we studied turbulence, mode coupling (Le., non-WKB reflection), and kinetic wave damping, and accomplishment lb provides the fourth topic: nonlinear steepening.

Wagner, William↗

NASA NTRS2003

Temporal Evolution of Nanostructures in a Model Nickel-Base Superalloy: Experiments and Simulations

The temporal evolution of the nanostructure of a model Ni-base superalloy (Ni-5.2 at.% Al-14.2 at.% Cr) is studied experimentally employing three-dimensional atom-probe (3DAP) microscopy in conjunction with kinetic Monte Carlo (KMC) simulations at 600 C. It is demonstrated that not only can the mean compositions of individual gamma' (Ni3Al with the Li2 structure) precipitates be measured but the Ni, Al, and Cr concentration profiles within the precipitates can also be determined for precipitates with a mean radius ( ) as small as 0.85 nm. The three asymptotic tinie dependencies of the Lifshitz-Slyzov-Wagner (LSW) theory of coarsening (Ostwald ripening) are measured and found to deviate from its theoretical predictions: possible explanations for these discrepancies are discussed. At 0.25 hr. there is 3DAP microscope evidence for the presence of precipitates of another nickel-rich phase. ="Ni3Cr" (Ni3Cr(1-x)Al(x)), which exhibits short-range order (SRO) and that is metastable with respect to Ni3Al. This metastable phase is also found by KMC simulations and has the composition Ni3Cr(1-x)Al(x), which is Ni-2.91 at.% Al-21.98 at.% Cr at 16 hours. Our results demonstrate that the decomposition of gamma the primary gamma (FCC) phase results in the concurrent formation of an ordered phase and a disordered phase by 0.25 hours.

Sudbrack, Chantal K.↗

NASA NTRSDate not supplied

A Combined Computational, Experimental, and Technology Development Approach to In-Space Laser Manufacturing Maturation at NASA Marshall Space Flight Center

In-space manufacturing (ISM) is emerging as a field vital to continued access and capabilities in the space environment. NASA Marshall Space Flight Center (MSFC) is advancing the frontier of in-space laser manufacturing (ISLM) techniques through work initially focused on maturing laser beam welding (LBW) and laser forming (LF) for use in space. Such techniques proffer the ability to assemble and join structures in space from sheet metal or other stock – extant satellites, in situ resource utilization of Lunar regolith, etc. – by forming to desired shapes and then joining via in-space welding (ISW). ISLM processes are useful for assembly, joining, modification, and repair of structures in free space and on the Lunar surface such as large observatories, antennas, trusses, blast/thermal/radiation shields, pressure vessels, and more. However, these techniques are not yet qualified & certified (Q&C) for regular application in space. It would be prohibitively expensive, laborious, and time-consuming to perform Q&C via traditional experimental approaches as data collection & experimentation in space is resource-intensive. As such, benchmark experiments and focused, properly instrumented technology demonstration efforts in space can collect sufficient data that – when combined with verified computational models in an integrated computational materials engineering (ICME) approach – can validate ICME tools capable of translating more readily obtained ground data to in-space, in situ, computationally informed Q&C of ISLM techniques. Several ISLM projects at MSFC are obtaining the data required to validate ICME tools through both ground and flight experiments. A parabolic flight experiment of LBW under vacuum is manifested for August 2024, including both microgravity and Lunar gravity profiles. This collaboration with the Ohio State University is investigating common aerospace alloys such as 316L stainless steel, 2219 aluminum alloy, and Ti64 titanium alloy. In situ data collection includes videography, thermography, and reference thermocouples to build a thermal model of the welds. This will elucidate the relevant physics when combined with post-flight microstructural examination and mechanical testing. MSFC is also progressing towards a suborbital flight experiment of LBW under vacuum, which could provide reams of data on ISW during sustained, high-quality reduced gravity. The effect of combined thermal (cryogenic and high-temperature) and vacuum exposure on both LBW (NASA-funded) and LF (DARPA-funded) is being investigated through ground experiments. In addition to the copious data collected during these ground experiments, ruggedization of LBW hardware will also be pursued. The datasets from these experiments will be used to validate computational models which will inform future ISLM efforts in an ICME framework. A variety of techniques across lengths scales, from CALPHAD-driven thermodynamics & kinetics to phase field modeling of solidification to kinetic Monte Carlo simulations of grain evolution at the mesoscale, will be employed to accelerate the infusion and eventual Q&C of LBW and LF for use in space. The development of data-driven surrogate models to bridge ground to flight experiments and thereby reduce the need for resource-intensive experiments in space will also be investigated. These ICME techniques, surrogate models, and datasets from ground testing can also be employed to advance manufacturing in terrestrial environments.

in-space welding↗

NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗