Search NASA⌕ Search

SEARCH · Search NASA

Results for “vacancy diffusion”

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

Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen

Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.

silicon oxidation↗

Dopant effect of yttrium and the growth and adherence of alumina on nickel-aluminum alloys

The atom superposition and electron delocalization molecular orbital theory and large cluster models have been employed to study cation vacancy diffusion in alpha-Al2O3 and the bonding of alpha-Al2O3 to nickel, aluminum, and yttrium surfaces. Al(3+) diffusion barriers in alpha-Al2O3 by the vacancy mechanism are in reasonable agreement with experiment. The barrier to Y(3+) diffusion is predicted to be much higher. Since addition of yttrium to transition metal alloys is known to reduce the growth rate and stress convolutions in protective alumina scales, this result suggests the rate-limiting step in scale growth is cation vacancy diffusion. This may partially explain the beneficial effect of yttrium dopants on scale adhesion. The theory also predicts a very strong bonding between alumina and yttrium at the surface of the alloy. This may also be important to the adhesion phenomenon. It is also found that aluminum and yttrium atoms bond very strongly to nickel because of charge transfer from their higher lying valence orbitals to the lower lying nickel s-d band.

Anderson, A. B.↗

Atomic-scale Modeling of the Structure and Dynamics of Dislocations in Complex Alloys at High Temperatures

We report on the progress made during the first year of the project. Most of the progress at this point has been on the theoretical and computational side. Here are the highlights: (1) A new code, tailored for high-end desktop computing, now combines modern Accelerated Dynamics (AD) with the well-tested Embedded Atom Method (EAM); (2) The new Accelerated Dynamics allows the study of relatively slow, thermally-activated processes, such as diffusion, which are much too slow for traditional Molecular Dynamics; (3) We have benchmarked the new AD code on a rather simple and well-known process: vacancy diffusion in copper; and (4) We have begun application of the AD code to the diffusion of vacancies in ordered intermetallics.

Daw, Murray S.↗

Tensile Creep of Polycrystalline Near-Stoichiometric NiAl

Long term tensile creep studies were conducted on binary NiAl in the temperature range 700-1200 K with the objectives of characterizing and understanding the creep mechanisms. Inverse and normal primary creep curves were observed depending on stress and temperature. It was concluded that the creep of NiAl is limited by dislocation mobility. The stress exponent for creep, n, increased from 5.5 at 1200 K to 13.9 at 700 K. The true activation energy for creep, Qc, was constant and equal to about 400 kJ per mole between 20 and 50 MPa but decreased to a constant value of 250 kJ per mole between 50 and 110 MPa. The activation energy was observed to be stress dependent above 110 MPa. The tensile creep results reported in this investigation were compared with compression creep data reported in the literature. A detailed discussion of the probable dislocation creep mechanisms governing compressive and tensile creep of NiAl is presented. It is concluded that the non-conservative motion of jogs on screw dislocations influenced the nature of the primary creep curves, where the climb of these jogs involves either the next nearest neighbor or the six-jump cycle vacancy diffusion mechanism. The probable nature of the atom vacancy exchange that occur within the core of an edge dislocation undergoing climb in NiAl are schematically examined.

Raj, Sai V.↗

Compression and Tensile Creep of Binary NiAl

Compression creep and long term tensile creep studies were conducted on cast and extruded binary NiAl in the temperature range 700-1200 K with the objectives of characterizing and understanding the creep mechanisms. Inverse and normal primary creep curves were observed in both compression and tension creep depending on stress and temperature although an asymmetrical response was observed under these two stress states. It was concluded that the primary creep of NiAl is limited by dislocation mobility. The stress exponents, n, for compression and tensile creep were similar varying between about 5 and 14. However, there were significant differences in the stress dependence of the activation energies for compression and tensile creep. The true activation energy for tensile creep, Q(sub c), was constant and equal to about 400 kJ/mol between 20 and 50 MPa but decreased to a constant value of 250 kJ/mol between 50 and 110 MPa. The activation energy was observed to be inversely stress dependent above 110 MPa. In contrast, Q(sub c) = 300 kJ/mol for compression creep was constant between 25 and 70 MPa and inversely dependent on the true stress above 70 MPa. A detailed discussion of the probable dislocation creep mechanisms governing compressive and tensile creep of NiAl is presented. It is concluded that the non-conservative motion of jogs on screw dislocations influenced the nature of the primary creep curves, where the climb of these jogs involves either the next nearest neighbor or the six-jump cycle vacancy diffusion mechanism. The probable natures of the atom-vacancy exchange that occur within the core of an edge dislocation undergoing climb in NiAl are schematically examined.

Raj, Sai V.↗

Defect distribution near the surface of electron-irradiated silicon

The surface-defect distributions of electron-irradiated n-type silicon have been investigated using a transient capacitance technique. Schottky, p-n junction, and MOS structures were used in profiling the defect distributions. Surface depletions of defects observed were attributed to the vacancy distribution, but not that of oxygen, and other capture centers' distributions. The vacancy diffusion length at 300 K was estimated to be about 3-6 microns.

Wang, K. L.↗

On the Oxidation of the Silicon Bond Coat in Environmental Barrier Coatings

Deal and Grove’s formulation is applied to model the oxidation of the silicon bond coat in environmental barrier coatings. The approach is similar to the author’s previous study, except now the oxidant is assumed to flow through the coating as oxygen ions via vacancies on the oxygen sublattice. The use of an “effective permeability” of the oxidant in the coating is required under the vacancy diffusion assumption. It is found that the linear-parabolic growth equation from Deal and Grove, 2 xox + Axox = Bt , is still applicable for the oxidation of the silicon bond coat, but the equation for the parameter A must be modified to account for the coating. The expression for B is the same as that obtained by Deal and Grove and thus is unaffected by the presence of a coating. The modified expression for A is a linear function of the coating thickness and includes the ratio of the oxidant permeability in the oxide to its effective permeability in the coating

Roy M Sullivan↗

On the Oxidation of the Silicon Bond Coat in Environmental Barrier Coatings

Deal and Grove’s formulation is applied to model the oxidation of the silicon bond coat in environmental barrier coatings. The approach is similar to the author’s previous study, except now the oxidant is assumed to flow through the coating as oxygen ions via vacancies on the oxygen sublattice. The use of an “effective permeability” of the oxidant in the coating is required under the vacancy diffusion assumption. It is found that the linear-parabolic growth equation from Deal and Grove, , is still applicable for the oxidation of the silicon bond coat, but the equation for the parameter A must be modified to account for the coating. The expression for B is the same as that obtained by Deal and Grove and thus is unaffected by the presence of a coating. The modified expression for A is a linear function of the coating thickness and includes the ratio of the oxidant permeability in the oxide to its effective permeability in the coating.

silicon oxidation, ionic diffusion, permeation, li↗

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

Oxygen Tracer Diffusion in LA(z-x) SR(X) CUO(4-y) Single Crystals

The tracer diffusion of O-18 in La(2-x)Sr(x)CuO(4-y) single crystals (x = 0 to 0.12) has been measured from 400 to 700 C in 1 atm of oxygen using SIMS analysis. Evidence for diffusion by a vacancy mechanism was found at low strontium contents. Oxygen diffusivities for x greater than or = 0.07 were depressed by several orders of magnitude below the diffusivity for undoped La2CuO(4+/-y). The observed effects of strontium doping on oxygen diffusivity are discussed in terms of defect chemical models. The decreasing oxygen diffusivity with increasing strontium was attributed to the ordering of oxygen vacancies at large defect concentrations. A diffusion anisotropy D(sub ab)/D(sub c) of nearly 600 was also found at 500 C.

OXYGEN DIFFUSIVITY↗

Theoretical Model of the Effect of Crack Tip Blunting on the Ultimate Tensile Strength of Welds in 2219-T87 Aluminum

A theoretical model representing blunting of a crack tip radius through diffusion of vacancies is presented. The model serves as the basis for a computer program which calculates changes, due to successive weld heat passes, in the ultimate tensile strength of 2219-T81 aluminum. In order for the model to yield changes of the same order in the ultimate tensile strength as that observed experimentally, a crack tip radius of the order of .001 microns is required. Such sharp cracks could arise in the fusion zone of a weld from shrinkage cavities or decohered phase boundaries between dendrites and the eutectic phase, or, possibly, from plastic deformation due to thermal stresses encountered during the welding process. Microstructural observations up to X2000 (resolution of about .1 micron) did not, in the fusion zone, show structural details which changed significantly under the influence of a heat pass, with the exception of possible small changes in the configuration of the interdendritic eutectic and in porosity build-up in the remelt zone.

Beil, R. J.↗

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

Diffusivity in Alumina Scales Grown on Al-MAX Phases

Ti3AlC2, Ti2AlC, and Cr2AlC are oxidation resistant MAX phase compounds distinguished by the formation of protective Al2O3 scales with well controlled kinetics. A modified Wagner treatment was used to obtain interfacial grain boundary diffusivity, deltaD(sub gb,O,int.), from scale growth rates and corresponding grain size. It is based on the p(O2)(exp -1/6) dependency of the double charged oxygen vacancy and oxygen diffusivity, coupled with the effective diffusion constant for short circuit grain boundary paths. Data from the literature for MAX phases was analyzed accordingly, and deltaD(sub gb,O,int.) was found to nearly coincide with the Arrhenius line developed for Zr-doped FeCrAl, where: deltaD(sub gb,O,int.) = 1.8x10(exp -10) exp(-375 kJ/RT) cubic meters/s. Furthermore, this oxidation relation suggests the more general format applicable to bulk samples under ambient conditions: deltaD(sub gb,O) = 7.567x10(exp -8) exp(-544 kJ/RT) p(O2)(exp -1/6) cubic meters/[s x Pa(exp -1/6)]. Data from many other FeCrAl(X) studies were similarly assessed to show general agreement with the relation for deltaD(sub gb,O,int.). This analysis reinforces the view that protective alumina scales grow by similar mechanisms for these Al-MAX phases and oxidation resistant FeCrAl alloys.

aluminum oxide↗