Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “impurity diffusion coefficients”

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

A comprehensive diffusion mobility database comprising 23 elements for magnesium alloys

We report that reliable experimental diffusion coefficients of 10 key alloying elements in Mg obtained by the present authors together with experimental data in the literature enabled us to perform a systematic test of the reliability of diffusion coefficients obtained from DFT calculations. The computed activation energy values were found to be quite accurate (mostly within 0.2 eV) but the computed pre-factors were less reliable. Such insights allowed us to develop a practical and yet robust strategy to perform diffusion mobility assessments by adopting the computed activation energy while fitting only the pre-factor when available experimental data are limited to a narrow temperature range. The overall good agreement between the DFT data and experimental data also gave us the confidence to employ the computed data for those that were still missing or inaccessible from experimental measurements. A systematic assessment of both the measured and computed diffusion data in hcp Mg was performed using the above holistic approach to yield the most comprehensive open Mg mobility database to date, comprising 23 elements (Mg, Ag, Al, Be, Ca, Cd, Ce, Cu, Fe, Ga, Gd, In, La, Li, Mn, Nd, Ni, Pu, Sb, Sn, U, Y, Zn). This more reliable mobility database will contribute to future development of advanced Mg alloys. The holistic approach developed in this study will be very beneficial to the future establishment of reliable mobility databases for other alloy systems as well.

36 MATERIALS SCIENCE↗

Albedo and flux extinction coefficient of impure snow for diffuse shortwave radiation

Impurities enter a snowpack as a result of fallout of scavenging by falling snow crystals. Albedo and flux extinction coefficient of soot contaminated snowcovers were studied using a two stream approximation of the radiative transfer equation. The effect of soot was calculated by two methods: independent scattering by ice grains and impurities and average refractive index for ice grains. Both methods predict a qualitatively similar effect of soot; the albedo is decreased and the extinction coefficient is increased compared to that for pure snow in the visible region; the infrared properties are largely unaffected. Quantitatively, however, the effect of soot is more pronounced in the average refractive index method. Soot contamination provides a qualitative explanation for several snow observations.

Choudhury, B. J.↗

Modeling turbulent impurity transport in the SOL of DIII-D with a reduced model

A novel impurity transport model that approximates SOL turbulence as a fluctuating poloidal electric field is shown to be an acceptable replacement for the traditional approach of assigning an arbitrary radial diffusion coefficient to the impurity ions. The model is implemented in the DIVIMP impurity transport code and applied to an L-Mode tungsten divertor experiment on DIII-D. The poloidal electric field is represented as fluctuating between ±1000 V m –1 based on previous measurements. The resulting intermittent v r = E θ × B T transport causes ions to transport both into the core as well as into the far-SOL. Simultaneous agreement with estimates of the W density just inside the separatrix as well as in the far-SOL is obtained (n W ~ 10 14 m –3 and n W ~ 10 12 m –3 , respectively). Prompt re-deposition of the W ions was necessary to obtain agreement (f redep ~ 99%). Here, we conclude that simulating impurity transport using a physics-based approximation for turbulence in the SOL, versus arbitrarily assigning diffusion coefficients, may enable better reactor scale predictions of core impurity contamination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing and controlling oxygen impurity diffusion in h -BN semi-bulk crystals

Combining its unique features of ultrawide bandgap (UWBG) and two-dimensional nature, h-BN has been explored for emerging applications such as deep ultraviolet optoelectronic devices and single photon emitters. One of the unusual applications of h-BN is for solid-state neutron detectors by utilizing the property of high thermal neutron capture cross section of B-10 as well as its UWBG properties. Although a record high detection efficiency of 59% has been attained by h-BN detectors, the understanding/minimization of defects and impurities is still needed to further advance the h-BN material and detector technologies. We report metal organic chemical vapor deposition growth and oxygen (O) impurity diffusion in thick h-BN. The diffusion coefficient (D) of O impurities has been measured via the evolution of an oxygen related emission with the etching depth, providing a value of D of ∼ 2 × 10−13 cm2/s at 1450 °C and supporting the interpretation that oxygen in h-BN is a substitutional donor. A multiple-buffer-layer approach was employed to mitigate to a certain degree the issue of oxygen diffusion from sapphire substrate during growth. It was demonstrated that the performance of h-BN neutron detectors fabricated from the wafer incorporating multiple buffer layers was significantly improved, as manifested by the enhanced thermal neutron detection efficiency. The advancement of the crystal growth technology of h-BN semi-bulk crystals creates applications in optoelectronic and power electronic devices utilizing the UWBG semiconductor properties of h-BN, while high efficiency h-BN neutron detectors have the potential to supplant the traditional He-3 gas detectors in various application areas by offering the obvious advantages of UWBG semiconductor technologies.

Physics↗

Resistive drift wave turbulence and anomalous transport of multi-species plasma

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.

Physics↗

A study of the properties of beryllium doped silicon with particular emphasis on diffusion mechanisms: Profiles of depth dependent conductivity as determined by electrical surface probes

Very large diffusion coefficients were encountered and required the determination of impurity profiles for samples approximately 1 cm thick. Since conductivity values are readily converted into concentrations of electrically active impurities, the major problem became that of accurately determining the conductivity profiles of beryllium diffused silicon samples. Four-point probe measurements on samples having depth conductivities are interpreted in terms of conductivity profiles, based on an exact solution of the problem of exponentially depth dependent conductivity. Applications include surface conductivity determination where the form of the conductivity profile is known, and conductivity profile determination from probe measurements taken as the sample surface is progressively lapped away. The application is limited to samples having conductivity monotonically decreasing with depth from the probed surface.

Franks, R. K.↗

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Oxygen diffusion in alpha-Al2O3

Oxygen self diffusion coefficients were determined in single crystal alpha-Al2O3 using the gas exchange technique. The samples were semi-infinite slabs cut from five different boules with varying background impurities. The diffusion direction was parallel to the c-axis. The tracer profiles were determined by two techniques, single spectrum proton activation and secondary ion mass spectrometry. The SIMS proved to be a more useful tool. The determined diffusion coefficients, which were insensitive to impurity levels and oxygen partial pressure, could be described by D = .00151 exp (-572kJ/RT) sq m/s. The insensitivities are discussed in terms of point defect clustering. Two independent models are consistent with the findings, the first considers the clusters as immobile point defect traps which buffer changes in the defect chemistry. The second considers clusters to be mobile and oxygen diffusion to be intrinsic behavior, the mechanism for oxygen transport involving neutral clusters of Schottky quintuplets.

Cawley, J. D.↗

Observation of impurity accumulation and its compatibility with high plasma performance in W7-X

At the W7-X stellarator, the bolometer system has measured an intensive radiation zone in the inner plasma region (at a normalized radius $ρ$ ~ 0.3–0.4) in the hydrogen plasma generated by electron cyclotron resonance heating; it differs from the normal plasma radiation distribution with an edge-localized emission zone. Spectroscopic diagnostics have recorded high-Z elements such as iron. This phenomenon happens in the plasma phases after gas supply turn-off, which results in all impurity relevant diagnostic signals increasing for several seconds. Despite the enhancement of the core radiation, the plasma energy confinement is improved. A transport analysis shows that this impurity radiation behavior is associated with a low diffusion coefficient ($D$~ 0.02 m 2 s -1 ) and a reversal of the convection around the radial position of the emission peak, which, under normal conditions, separates the zones of outward convection in the central (|$V$| ~ 0.1 m s -1 ) and inward convection in the outer region (|$V$| ~ 0.3 m s -1 ). An impurity accumulation around this radial position has been identified. The transport coefficients obtained are comparable with the theoretical predictions of collisional impurity transport. In the plasma phases studied, both impurity and energy confinement are enhanced. The mechanism responsible for the improvement is believed to be a reduction of micro-instabilities associated with the observed steepening of the density profile, initiated by a low edge plasma density (<1.0 × 10 19 m -3 ) after switching off the gas fueling. The normalized temperature and density gradients fulfil the condition for the suppression of ITG turbulence.

suppression of ITG turbulence↗

Diffusion, atomic transport, and ordering in Al-Zr alloys: FCC and liquid phases

Additive manufacturing of materials with controlled microstructure demands knowledge of atomic scale properties near the solid-liquid transition state. Many of these properties are not affordable by experimental techniques and computer modeling is the possible solution to the problem. In this paper, we present the results of an extended atomistic study of intrinsic atomic transport due to vacancy diffusion in FCC and L12 solid phases and diffusion in the liquid phase of Al-Zr alloys. A deceleration of the overall self-diffusion was observed when Zr was added to Al. The effect was stronger in the solid and weaker in the liquid. Additionally, the effect was strongly temperature dependent in the solid phases, but not in the liquid. Atomic transport was chemically biased: transport of Zr atoms was significantly slower than that of Al atoms, and this bias effect was stronger in the solid phases. The overall diffusion and chemical ordering processes in the liquid state were five to six orders in magnitude faster than in the solid. Chemical short-range order parameters in the liquid saturated at values close to those in the ordered L12 structure of Al3Zr. Chemical and structural ordering in the solid phases was negligible over the modeled microsecond time scale. Here, the results are discussed in view of optimizing additive manufacturing parameters for the controlled formation of metastable L1 2 precipitates.

36 MATERIALS SCIENCE↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Tritium adsorption and absorption on (100) and (001) surfaces of pure and tin defective zirconium

Zirconium alloys such as zircaloy-4 are used as tritium (T) getter materials in tritium-producing burnable absorber rods (TPBARs) due to their ability to capture T, thereby forming metal hydrides. Developing an understanding of T adsorption onto zircaloy prior to diffusion into the subsurface is relevant for rational tritium getter and TPBAR design, to improve material properties for nuclear applications. Herein, density functional theory calculations revealed the preferred binding sites for T adsorption on Zr(001) and Zr(100). The energy barriers of T transfer, along the surface and from the surface to the subsurface were computed. The adsorption properties of Zr(001) were found to be superior to those of Zr(100). Surface tin impurities were found to strongly repel T. The presence of subsurface and surface tin resulted in higher absorption energy barriers for both the forward and reverse processes. Based on the calculated energy barriers, a surface to surface T diffusion coefficient of 9.53 × 10 -10 m 2 s -1 is expected for pristine Zr(001). Finally, a surface to subsurface T diffusion coefficient on the order of 10 -13 m 2 s -1 is predicted in pristine Zr, decreasing to 10 -19 m 2 s -1 for the transfer with a subsurface tin impurity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of positive ion implantation into antireflection coating of silicon solar cells

The state of technological development of Si solar cells for highest obtained efficiency and radiation resistance is summarized. The various theoretical analyses of Si solar cells are reviewed. It is shown that factors controlling blue response are carrier diffusion length, surface recombination, impurity concentration profile in surface region, high level of surface impurity concentration (degeneracy), reflection coefficient of oxide, and absorption coefficient of Si. The theory of ion implantation of charge into the oxide antireflection coating is developed and side effects are discussed. The experimental investigations were directed at determining whether the blue response of Si solar cells could be improved by phosphorus ion charges introduced into the oxide antireflection coating.

Middleton, A. E.↗

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗

Electric field effects during disruptions

Tokamak disruptions are associated with breaking magnetic surfaces, which makes magnetic field lines chaotic in large regions of the plasma. The enforcement of quasi-neutrality in a region of chaotic field lines requires an electric potential that has both short and long correlation distances across the magnetic field lines. The short correlation distances produce a Bohm-like diffusion coefficient ∼Te/eB and the long correlation distances aT produce a large scale flow ∼Te/eBaT. This cross-field diffusion and flow are important for sweeping impurities into the core of a disrupting tokamak. The analysis separates the electric field in a plasma into the sum of a divergence-free, E→B, and a curl-free, E→q, part, a Helmholtz decomposition. The divergence-free part of E→ determines the evolution of the magnetic field. The curl-free part enforces quasi-neutrality, E→q=−∇→Φq. Magnetic helicity evolution gives the required boundary condition for a unique Helmholtz decomposition and an unfortunate constraint on steady-state tokamak maintenance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗