Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion rate”

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

High-resolution measurements, line identification, and spectral modeling of K-alpha transitions in Fe XVIII-Fe XXV

A detailed analysis of the iron K-alpha emission spectrum covering the wavelength region from 1.840 to 1.940 A is presented. Measurements are made with a high-resolution Bragg crystal spectrometer on the Princeton Large Torus (PLT) tokamak for plasma conditions which closely resemble those of solar flares. A total of 40 features are identified, consisting of either single or multiple lines from eight charge states in iron, Fe XVIII - Fe XXV, and their wavelengths are determined with an accuracy of 0.1-0.4 mA. Many of these features are identified for the first time. In the interpretation of our observations we rely on model calculations that determine the ionic species abundances from electron density and temperature profiles measured independently with nonspectroscopic techniques and that incorporate theoretical collisional excitation and dielectronic recombination rates resulting in the excitation of the 1s2sr2ps configurations. The model calculations also include the effect of diffusive ion transport. Good overall agreement between the model calculations and the observations is obtained, which gives us confidence in our line identifications and spectral modeling capabilities. The results are compared with earlier analyses of the K-alpha emission from the Sun.

Beiersdorfer, P.↗
NASA NTRS1994

Active chlorine and nitric oxide formation from chemical rocket plume afterburning

Chlorine and oxides of nitrogen (NO(x)) released into the atmosphere contribute to acid rain (ground level or low-altitude sources) and ozone depletion from the stratosphere (high-altitude sources). Rocket engines have the potential for forming or activating these pollutants in the rocket plume. For instance, H2/O2 rockets can produce thermal NO(x) in their plumes. Emphasis, in the past, has been placed on determining the impact of chlorine release on the stratosphere. To date, very little, if any, information is available to understand what contribution NO(x) emissions from ground-based engine testing and actual rocket launches have on the atmosphere. The goal of this work is to estimate the afterburning emissions from chemical rocket plumes and determine their local stratospheric impact. Our study focuses on the space shuttle rocket motors, which include both the solid rocket boosters (SRB's) and the liquid propellant main engines (SSME's). Rocket plume afterburning is modeled employing a one-dimensional model incorporating two chemical kinetic systems: chemical and thermal equilibria with overlayed nitric oxide chemical kinetics (semi equilibrium) and full finite-rate chemical kinetics. Additionally, the local atmospheric impact immediately following a launch is modeled as the emissions diffuse and chemically react in the stratosphere.

Leone, D. M.↗
NASA NTRS1960

Diffusion of Sound Waves in a Turbulent Atmosphere

The directional and frequency diffusion of a plane monochromatic 2 sound wave in statistically homogeneous, isotropic, and stationary turbulence is analyzed theoretically. The treatment is based on the diffusion equation for the energy density of sound waves, using the scattering cross section derived by Kraichnan for the type of turbulence assumed here. A form for the frequency-wave number spectrum of the turbulence is adopted which contains the pertinent parameters of the flow and is adapted to ease of calculation. A new approach to the evaluation of the characteristic period of the flow is suggested. This spectrum is then related to the scattering cross section. Finally, a diffusion equation is derived as a small-angle scattering approximation to the rigorous transport equation. The rate of spread of the incident wave in frequency and direction is calculated, as well as the power spectrum and autocorrelation for the wave.

Lyon, Richard H.↗
NASA NTRS1999

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

Silicon carbide (SiC) and Si3N4 materials were tested in various turbine engine combustion environments chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high-speed aircraft. Representative chemical vapor-deposited (CVD), sintered, and composite materials were evaluated by furnace and high-pressure burner rig exposures. Although protective SiO2 scales formed in all cases, the evidence presented supports a model based on paralinear growth kinetics (i.e., parabolic growth moderated simultaneously by linear volatilization). The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were thus used to map SiO2 volatility (and SiC recession) over a range of temperatures, pressures, and velocities. The functional dependency of material recession (volatility) that emerged followed the form A[exp(-Q / RT)](P(sup x)v(sup y). These empirical relations were compared with rates predicted from the thermodynamics of volatile SiO and SiOxHy reaction products and a kinetic model of diffusion through a moving boundary layer. For typical combustion conditions, recession of 0.2 to 2 micrometers/hr is predicted at 1200 to 1400 C, far in excess of acceptable long-term limits.

Smialek, James L.↗
NASA NTRS2002

A Comparison of the Lower Stratospheric Age-Spectra Derived from a General Circulation Model and Two Data Assimilation Systems

We use kinematic and diabatic back trajectory calculations, driven by winds from a general circulation model (GCM) and two different data assimilation systems (DAS), to compute the age spectrum at three latitudes in the lower stratosphere. The age-spectra are compared to chemical transport model (CTM) calculations, and the mean ages from all of these studies are compared to observations. The age spectra computed using the GCM winds show a reasonably isolated tropics in good agreement with observations; however, the age spectra determined from the DAS differ from the GCM spectra. For the DAS diabatic trajectory calculations there is too much exchange between the tropics and mid-latitudes. The age spectrum is thus too broad and the tropical mean age is too old as a result of mixing older mid latitude air with tropical air. Likewise the mid latitude mean age is too young due to the in mixing of tropical air. The DAS kinematic trajectory calculations show excessive vertical dispersion of parcels in addition to excessive exchange between the tropics and mid latitudes. Because air is moved rapidly to the troposphere from the vertical dispersion, the age spectrum is shifted toward the young side. The excessive vertical and meridional dispersion compensate in the kinematic case giving a reasonable tropical mean age. The CTM calculation of the age spectrum using the DAS winds shows the same vertical and meridional dispersive characteristics of the kinematic trajectory calculation. These results suggest that the current DAS products will not give realistic trace gas distributions for long integrations; they also help explain why the extra tropical mean ages determined in a number of previous DAS driven CTM s are too young compared with observations. Finally, we note trajectory-generated age spectra . show significant age anomalies correlated with the seasonal cycles. These anomalies can be linked to year-to-year variations in the tropical heating rate. The anomalies are suppressed in the CTM spectra suggesting that the CTM transport scheme is too diffusive.

Schoeberl, Mark R.↗
NASA NTRS2001

Optimized suspension culture: the rotating-wall vessel

Suspension culture remains a popular modality, which manipulates mechanical culture conditions to maintain the specialized features of cultured cells. The rotating-wall vessel is a suspension culture vessel optimized to produce laminar flow and minimize the mechanical stresses on cell aggregates in culture. This review summarizes the engineering principles, which allow optimal suspension culture conditions to be established, and the boundary conditions, which limit this process. We suggest that to minimize mechanical damage and optimize differentiation of cultured cells, suspension culture should be performed in a solid-body rotation Couette-flow, zero-headspace culture vessel such as the rotating-wall vessel. This provides fluid dynamic operating principles characterized by 1) solid body rotation about a horizontal axis, characterized by colocalization of cells and aggregates of different sedimentation rates, optimally reduced fluid shear and turbulence, and three-dimensional spatial freedom; and 2) oxygenation by diffusion. Optimization of suspension culture is achieved by applying three tradeoffs. First, terminal velocity should be minimized by choosing microcarrier beads and culture media as close in density as possible. Next, rotation in the rotating-wall vessel induces both Coriolis and centrifugal forces, directly dependent on terminal velocity and minimized as terminal velocity is minimized. Last, mass transport of nutrients to a cell in suspension culture depends on both terminal velocity and diffusion of nutrients. In the transduction of mechanical culture conditions into cellular effects, several lines of evidence support a role for multiple molecular mechanisms. These include effects of shear stress, changes in cell cycle and cell death pathways, and upstream regulation of secondary messengers such as protein kinase C. The discipline of suspension culture needs a systematic analysis of the relationship between mechanical culture conditions and biological effects, emphasizing cellular processes important for the industrial production of biological pharmaceuticals and devices.

Non-NASA Center↗
DOE OSTI2020

Transport and kinetics properties of LaF 3 in FLiNaK molten salt determined by electrochemical methods

It is found that lanthanum species can exist in the form of LaF 6 3- in FLiNaK (46.5 LiF-11.5 NaF-42 KF, mol%) molten salt in previous studies. The understanding of transport and reaction kinetics parameters of LaF 6 3- ions in FLiNaK molten salt is essential for the electrochemical separation of lanthanum from FLiNaK melt to achieve the coolant clean-up in fluoride salt cooled high-temperature reactor (FHR). Through conducting chronopotentiometry tests, the diffusion coefficients of LaF 6 3- ions in FLiNaK molten salt at the temperature range of 923 K–1023 K were determined. By conducting potentiodynamic polarization measurements at different concentrations and temperatures, the exchange current density, reaction rate constant, and charge transfer coefficient were obtained. The experiment data of potentiodynamic polarization were analyzed using a non-simplified electrode kinetics equation which incorporates both mass transfer and reaction kinetics.

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

Understanding and Controlling the Formation of Nonradiative Defects in Blue Organic Triplet Emitters

Phosphorescent organic light-emitting devices (PHOLEDs) suffer from destructive molecular processes due to triplet-polaron and triplet-triplet annihilation. These processes are energetically driven and hence are particularly active in decreasing the lifetime of blue PHOLEDs. It has recently been shown that increasing triplet radiative rates via the Purcell effect effectively extends the device operational lifetime by reducing the triplet radiative lifetime, thus decreasing their density and the probability of triplet-annihilation reactions. We provide an analytical framework using Marcus theory to explain the observed, approximately exponential relationship between exciton energy and device lifetime. From transient drift-diffusion dynamics, we show that the Purcell effect reduces the exciton density in the steady state and increases the photoluminescent yield, thereby reducing defect generation rates and extending the device lifetime. We control the radiative rate of excitons in microcavities, thereby connecting the exciton energy and decay rates with the observed device lifetime. The device lifetime is shown to follow a power-law dependence on the Purcell factor ( PF m ) with m = 1.5 to 2.5, dependent on the TTA-to-TPA ratio and photoluminescence quantum yield. From our analysis, a fivefold increase in PF has the potential to extend the blue PHOLED lifetime by up to 2 orders of magnitude, making the blue PHOLED lifetime comparable to that of state-of-the-art green PHOLEDs. Published by the American Physical Society 2024

Zhao, Haonan (ORCID:0000000190584422)↗
DOE OSTI2022

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2024

Experimental and modeling study of the effect of confinement on the thermal decomposition of organic materials

Simultaneous Thermal Analysis is widely used to study thermal decomposition of a variety of organic materials, and the choice of experimental conditions is very important to obtain reliable mass loss and heat flow curves. Here in this work we investigate the effect of confinement on the thermal decomposition of organic materials. As an example, we explore high-density polyethylene (HDPE) and one of its pyrolysis products, namely eicosane (C 20 H 42 ), through experiments and modeling in open crucibles and crucibles closed with a pierced lid. A new model is developed for evaporation and sublimation from open and pinhole thermal-analysis pans. It considers the kinetic resistance of evaporation at the liquid-vapor interface, inhibition of evaporation by buildup of vapor inside the pan, thermal expansion of vapor inside the pan, and diffusion of vapor from the pan. The model is validated using simultaneous heat-flow and mass-loss measurements of n-eicosane evaporation for various pinhole sizes, down to 50 µm. The importance of product evaporation inhibition for measuring polymer decomposition was demonstrated using HDPE. This study sheds light on the effect of confinement on the mass loss rate of organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Radial diffusion of low-energy plasma ions in Saturn's magnetosphere

Radial diffusion of low-energy plasma ions in Saturn's magnetosphere is investigated using a comprehensive set of equations for radial diffusion that incorporate distributed sources and sinks of ions. The results of calculations indicate that the radial-diffusion transport of low-energy O(+) ions with a source in the neutral H2O cloud of the satellites Dione and Tethys can account for Voyager observations of thermal heavy ions in Saturn's magnetosphere. The source rate was calculated to be about 10 to the 26th O(+) ions/sec, in good agreement with the sputtering calculations of Johnson et al. (1989). It is estimated that, due to fast radial diffusion, the residence time of O(+) ions in the Dione-Tethys torus is about 30 days, sufficiently short to account for the plasma density observed there. The densities of hot H(+) and N(+) resulting from the ionization and pickup of Titan's neutral clouds in the outer magnetosphere can also be accounted for within the framework of diffusive ion transport.

Barbosa, D. D.↗
DOE OSTI2023

Reversible Reactions, Mesh Size, and Segmental Dynamics Control Penetrant Diffusion in Ethylene Vitrimers

The diffusion of two aromatic dyes with nearly identical sizes was measured in ethylene vitrimers with precise linker lengths and borate ester cross-links using fluorescence recovery after photobleaching (FRAP). One dye possessed a reactive hydroxyl group, while the second was inert. The reaction of the hydroxyl group with the network is slow relative to the hopping times of the dye, resulting in a large slowdown by a factor of 50 for a reactive probe molecule. A kinetic model was fit to the fluorescence intensity data to determine rate constants for the reversible reaction of the dye from the network, which confirms the role of slow reaction kinetics. A second network cross-linker was also investigated with a substituted boronic ester showing ∼10,000 times faster exchange kinetics. In this system, the two dyes show the same diffusion coefficient, as the reaction is no longer the rate-limiting step. The role of dense meshes on small and large dyes is also discussed in the context of the existing theories. Finally, these results highlight the potential of dynamic networks to control penetrant transport through synergistic effects of the mesh size, dynamic bond kinetics, and penetrant–network interactions.

confinement↗
DOE OSTI2026

Dates and rates of loess deposition and soil formation on the Snake River Plain, Idaho, USA

Loess-paleosol sequences in eolian deposits on the Snake River Plain, Idaho, western United States, preserve records of late Pleistocene to Holocene glacial-interglacial cycles. Here, we examine climate-driven changes in the timing and rate of loess accretion, soil formation, and pedogenic carbonate accumulation in the eastern and central Snake River Plain. High rates of loess deposition often correspond with dry, cold, and windy conditions and/or high sediment supply, while soil formation indicates landscape stability. Optically stimulated luminescence dating of four loess-soil sequences provides a record of the timing of loess deposition, and radiocarbon dating ( 14 C) of soil inorganic carbon provides ages for soil formation and pedogenic carbonate precipitation. On the central Snake River Plain, marine isotope stage (MIS) 3 is marked by loess deposition between ca. 50 ka and 33 ka, and pedogenic carbonate accumulation between ca. 42 ka calibrated years before present (cal. B.P.) and 31 ka cal. B.P. During MIS 2, loess deposition between ca. 21 ka and 15 ka is contemporaneous with pedogenic carbonate precipitation. MIS 2 is characterized by diffuse pedogenic carbonate formation and rapid loess accumulation. Both the eastern and central Snake River Plain study sites experienced high rates of loess deposition during MIS 2; however, the eastern Snake River Plain had almost three times faster rates of loess accumulation (0.56−0.63 m/k.y. from ca. 25 ka to 17 ka) as compared to the central Snake River Plain (0.18−0.26 m/k.y. from ca. 21 ka to 15 ka). During MIS 1, loess deposition on the central Snake River Plain ca. 11 ka was coincident with pedogenic carbonate formation ca. 10 ka cal. B.P. MIS 3 and MIS 1 are characterized by A horizons and well-developed Bk horizons, indicating soil stability and lower rates of loess accumulation during these intervals. Finally, this study highlights the role of late Quaternary−scale climatic variability on rates and amounts of soil inorganic carbon storage.

Geosciences↗
DOE OSTI2021

Structural Evolution and Transition Dynamics in Lithium Ion Battery under Fast Charging: An Operando Neutron Diffraction Investigation

Abstract Fast charging (<15 min) of lithium‐ion batteries (LIBs) for electrical vehicles (EVs) is widely seen as the key factor that will greatly stimulate the EV markets, and its realization is mainly hindered by the sluggish diffusion of Li + . To have a mechanistic understanding of Li + diffusion within LIBs, in this study, structural evolutions of electrodes for a Ni‐rich LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) || graphite cylindrical cell with high areal loading (2.78 mAh cm −2 ) are developed for operando neutron powder diffraction study at different charging rates. Via sequential Rietveld refinements, changes in structures of NMC622 and Li x C 6 are obtained during moderate and fast charging (from 0.27 C to 4.4 C). NMC622 exhibits the same structural evolution regardless of C‐rates. For phase transitions of Li x C 6 , the stage I (LiC 6 ) phase emerges earlier during the stepwise intercalation at a lower state of charge when charging rate is increased. It is also found that the stage II (LiC 12 ) → stage I (LiC 6 ) transition is the rate‐limiting step during fast charging. The LiC 12 → LiC 6 transition mechanism is further analyzed using the Johnson–Mehl–Avrami–Kolmogorov model. It is concluded as a diffusion‐controlled, 1D phase transition with decreasing nucleation kinetics under increasing chargingrates.

25 ENERGY STORAGE↗
NASA NTRS1982

Oceanic corrosion test of bare and zinc-protected aluminum alloys for seawater heat exchangers

Bare 3004 tubes, 7072 Alclad 3004 tubes, and bare and zinc diffusion treated 3003 extrusions from a brazed aluminum, plate-fin heat exchanger were exposed to 1.8 m/sec flowing seawater aboard an open ocean test facility moored 3.4 km off the southeast coast of Puerto Rico. After six months exposure, the average corrosion rates for most varieties of aluminum materials converged to a low value of 0.015 mm/yr (0.6 mils/yr). Pitting did not occur in bare 3003 and 3004 samples during the six month test. Pitting did occur to varying degrees in the Alclad and zinc diffusion treated material, but did not penetrate to the base metal. Biofouling countermeasures (intermittent chlorination and brushing) did not affect the corrosion rates to any significant extent. Intermittent chlorination at a level of 0.5 ppm for 28 minutes daily controlled microbiofouling of the samples but did not prevent the development of a macrobiofouling community in areas of the plumbing with low flow.

Sasscer, D. S.↗
NASA NTRS1997

High Pressure Gas Permeation and Liquid Diffusion Studies of Coflon and Tefzel Thermoplastics

The life of fluid-carrying flexible or umbilical pipes during service at elevated temperatures and pressures depends inter alia on their resistance to attack by the fluids present and the rate at which these fluids are absorbed by the pipe lining materials. The consequences of fluid ingress into the thermoplastic lining could mean a) a reduction in its mechanical strength, to increase chances of crack formation and growth and thus a loss of integrity, b) the occurrence of permeation right through the lining material, with pressure build- up in the outer pipe wall construction (of flexible pipes) or chemical attack (from a hostile permeant) on outer layers of reinforcements. Therefore it is important within this project to have relevant permeation data for Coflon and Tefzel thermoplastics: the former is plasticised, the latter is not. A previous report (CAPP/M.2) described experimental equipment and techniques used by MERL when measuring high pressure (up to 5000 psi) gas permeation and liquid diffusion through thermoplastic samples cut from extruded bar or pipe, and provided the basic theory involved. Norsk Hydro are also performing gas permeation tests on pipe sections, at up to 100 bars (1450 psi) pressure or so, and reporting separately. Some comparisons between data from Norsk Hydro and MERL have been made herein. The tests should be considered as complementary, as the Norsk Hydro test has the obvious benefit of using complete pipe sections, whilst MERL can test at much higher pressures, up to 1000 bar if necessary. The sophisticated analytical measuring equipment of Norsk Hydro can distinguish the individual components of mixed gases and hence the various permeation-linked coefficients whereas MERL, in using pressure increase at constant volume to determine permeation rate, is limited to obtaining single gas data, or apparent (or representative) coefficients for a mixed gas as a whole. Except for the initial fluid diffusion data for Tefzel described in CAPP/M.2, the present report covers all aspects of fluid permeation and diffusion for Coflon and Tefzel, including all the pen-neation data accumulated in the project to date. Test gases have mainly been methane (CH4) and carbon dioxide (CO2). More high pressure (HP) gas permeation tests have been performed since the last issue of this report, most being concerned with changes in permeation characteristics brought about by ageing in various relevant fluids. This revision supersedes previous issues.

Morgan, G. J.↗
DOE OSTI2020

In Operando Calorimetric Measurements for Activated Carbon Electrodes in Ionic Liquid Electrolytes under Large Potential Windows

This study aims to investigate the effect of the potential window on heat generation in carbon-based electrical double layer capacitors (EDLCs) with ionic-liquid (IL)-based electrolytes using in operando calorimetry. The EDLCs consisted of two identical activated-carbon electrodes with either neat 1-butyl-1-methylpyrrolidinium bis(trifluoromethane-sulfonyl)imide ([Pyr 14 ][TFSI]) electrolyte or 1.0 m [Pyr 14 ][TFSI] in propylene carbonate (PC) as electrolyte. The instantaneous heat generation rate at each electrode was measured under galvanostatic cycling for different potential windows ranging from 1 to 4 V. First, the heat generation rates at the positive and negative electrodes differed significantly in neat IL owing to the differences in the ion sizes and diffusion coefficients. However, these differences were minimized when the IL was diluted in PC. Second, for EDLC in neat [Pyr 14 ][TFSI] at high potential window (4 V), a pronounced endothermic peak was observed at the beginning of the charging step at the positive electrode owing to TFSI - intercalation in the activated carbon. On the other hand, for EDLC in 1.0 m [Pyr 14 ][TFSI] in PC at potential window above 3 V, an endothermic peak was observed only at the negative electrode owing to the decomposition of PC. Third, for both neat and diluted [Pyr 14 ][TFSI] electrolytes, the irreversible heat generation rate increased with increasing potential window and exceeded Joule heating. This was attributed to the effect of potential-dependent charge redistribution resistance. Additionally, a further increase in the irreversible heat generation rate was observed for the largest potential windows owing to the degradation of the PC solvent. Finally, for both types of electrolyte, the reversible heat generation rate increased with increasing potential window because of the increase in the amount of ion adsorbed/desorbed at the electrode/electrolyte interface.

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

Accelerated kinetic Monte Carlo: A case study; vacancy and dumbbell interstitial diffusion traps in concentrated solid solution alloys

Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.

36 MATERIALS SCIENCE↗