Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusivity”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Short-range thermal magnon diffusion in magnetic garnet

Here, using the spin Seebeck effect (SSE), we study the propagation distance of thermally induced spin currents inside a magnetic insulator thin film in the short-range regime. We disambiguate spin currents driven by temperature and chemical potential gradients by comparing the SSE signal before and after adding a heat-sinking capping layer on the same device. We report that the measured spin decay behavior near the heat source is well accounted for by a diffusion model where the magnon diffusion length is in submicron range, in other words, two orders of magnitude smaller than previous estimates inferred from the long-range behavior. Our results highlight the caveat in applying a diffusive theory to describe thermally generated magnon transport, where a single decay length may not capture the behavior on all length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Antiphase boundary migration as a diffusion mechanism in a P3 sodium layered oxide

The electrode materials of Na-ion batteries must exhibit fast ion transport, a property that is often negatively impacted when mobile cations order. Many layered Na (and K) intercalation compounds adopt the P3 crystal structure, which tends to stabilize cation-vacancy orderings over a wide range of Na concentrations. The ordered phases have been predicted to consist of domains that are periodically separated by antiphase boundaries (APBs) to accommodate variations in Na concentration. The mechanisms and rates with which Na diffuses within such ordered phases of P3 are not understood. We have conducted a first-principles study of Na diffusion in the P3 structure, using Na x CoO 2 as a model system. We find that most simple nearest-neighbor hops in Na-vacancy ordered phases in P3 result in unstable configurations and are therefore forbidden. We identify an alternative mechanism of Na transport that is mediated by APB migration via the formation and expansion of kinks along domain boundaries. The limiting kinetic barriers for this mechanism are quite low (0.03 to 0.30 eV). Vacancy defects are not required, and furthermore have formation energies that are higher than those of the APB kink defects. Furthermore, our results suggest that APB migration is a fundamental diffusion mechanism in the P3 structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time-of-flight modulated intensity small-angle neutron scattering measurement of the self-diffusion constant of water

We report that the modulated intensity by zero effort small-angle neutron scattering (MI-SANS) technique is used to measure scattering with a high energy resolution on samples normally ill-suited for neutron resonance spin echo. The self-diffusion constant of water is measured over a q–t range of 0.01–0.2 Å -1 and 70–500 ps. In addition to demonstrating the methodology of using time-of-flight MI-SANS instruments to observe diffusion in liquids, the results support previous measurements on water performed with different methods. This polarized neutron technique simultaneously measures the intermediate scattering function for a wide range of time and length scales. Two radio frequency flippers were used in a spin-echo setup with a 100 kHz frequency difference in order to create a high-resolution time measurement. The results are compared with self-diffusion measurements made by other techniques and the general applicability of MI-SANS at a pulsed source is assessed.

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

Integral X-ray diffuse scattering for studying irradiation-induced dislocation loops in single crystals

Numerical diffuse scattering cross-section calculations are used to establish a rigorous basis for determining the concentration and size distribution of dis­location loops in irradiated single crystals from integral X-ray diffuse scattering (XRDS) measurements. Differential XRDS intensities for prismatic {111} type dislocation loops are numerically calculated as a function of loop radius R and wavevectors q relative to Bragg reflections in tungsten. The results show the well known 1/q 2 Huang scattering form at small q that transitions to a ∼1/q 4 dependence associated with the Stokes–Wilson approximation for q ≳ 1/R. More importantly, they show further that the 1/q 4 falloff is not the asymptotic large-q form of the diffuse scattering for small loops (R < 200 Å) as has often been assumed. Rather, for loop sizes as small as R ≃ 5 Å with strong curvature, the calculations show definitively that the scattering transitions to a robust 1/q 5 falloff at larger q that arises due to the local strains near the dislocation core defining the circumference of the dislocation loops. The presence of this 1/q 5 asymptotic form for both small and large loops is experimentally confirmed using an integral XRDS measurement around the 110 reflection on self-ion-irradiated tungsten combined with numerically calculated integral XRDS cross-sections. Accordingly, the historical two-region theoretical treatment of the cross-sections for integral XRDS is extended to a three-region model that has direct sensitivity to the (first-moment) dislocation line lengths of dislocation loops. In conclusion, these developments enable the use of both numerical and analytically modeled cross-sections to make accurate integral XRDS determinations of dislocation loop sizes and concentrations using modest-intensity laboratory X-ray sources.

42 ENGINEERING↗
DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Diffusion-assisted molecular beam epitaxy of CuCrO 2 thin films

Using molecular beam epitaxy (MBE) to grow multielemental oxides (MEOs) is generally challenging, partly due to difficulty in stoichiometry control. Occasionally, if one of the elements is volatile at the growth temperature, stoichiometry control can be greatly simplified using adsorption-controlled growth mode. Otherwise, stoichiometry control remains one of the main hurdles to achieving high-quality MEO film growths. Here, in this paper, we report another kind of self-limited growth mode, dubbed diffusion-assisted epitaxy, in which excess species diffuses into the substrate and leads to the desired stoichiometry, in a manner similar to the conventional adsorption-controlled epitaxy. Specifically, we demonstrate that using diffusion-assisted epitaxy, high-quality epitaxial CuCrO 2 films can be grown over a wide growth window without precise flux control using MBE.

36 MATERIALS SCIENCE↗
DOE OSTI2021

AIR Algebraic Multigrid for a Space-Time Hybridizable Discontinuous Galerkin Discretization of Advection(-Diffusion)

This article investigates the efficiency, robustness, and scalability of approximate ideal restriction (AIR) algebraic multigrid as a preconditioner in the all-at-once solution of a space-time hybridizable discontinuous Galerkin discretization of advection-dominated flows. The motivation for this study is that the time-dependent advection-diffusion equation can be seen as a “steady” advection-diffusion problem in $(d+1)$-dimensions and AIR has been shown to be a robust solver for steady advection-dominated problems. Numerical examples demonstrate the effectiveness of AIR as a preconditioner for advection-diffusion problems on fixed and time-dependent domains, using both slab-by-slab and all-at-once space-time discretizations, and in the context of uniform and space-time adaptive mesh refinement. A closer look at the geometric coarsening structure that arises in AIR also explains why AIR can provide robust, scalable, space-time convergence on advective and hyperbolic problems, while most multilevel parallel-in-time schemes struggle with such problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Diffusion-Model-Assisted Supervised Learning of Generative Models for Density Estimation

Here, we present a supervised learning framework of training generative models for density estimation. Generative models, including generative adversarial networks (GANs), normalizing flows, and variational auto-encoders (VAEs), are usually considered as unsupervised learning models, because labeled data are usually unavailable for training. Despite the success of the generative models, there are several issues with the unsupervised training, e.g., requirement of reversible architectures, vanishing gradients, and training instability. To enable supervised learning in generative models, we utilize the score-based diffusion model to generate labeled data. Unlike existing diffusion models that train neural networks to learn the score function, we develop a training-free score estimation method. This approach uses mini-batch-based Monte Carlo estimators to directly approximate the score function at any spatial-temporal location in solving an ordinary differential equation (ODE), corresponding to the reverse-time stochastic differential equation (SDE). This approach can offer both high accuracy and substantial time savings in neural network training. Once the labeled data are generated, we can train a simple, fully connected neural network to learn the generative model in the supervised manner. Compared with existing normalizing flow models, our method does not require the use of reversible neural networks and avoids the computation of the Jacobian matrix. Compared with existing diffusion models, our method does not need to solve the reverse-time SDE to generate new samples. As a result, the sampling efficiency is significantly improved. We demonstrate the performance of our method by applying it to a set of 2D datasets as well as real data from the University of California Irvine (UCI) repository.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method

We show that non-diffusive volume advection in two-dimensions is achieved with several benchmark problems using a newly developed high-order volume of fluids (VOF) interface reconstruction method. (1) A new VOF interface reconstruction method using circular/corner facets (linear facets are a degenerate case of arcs). We create a circular interface facet in each mixed zone by matching neighbor volume with a hybrid Newton’s-bisection method and the local solution is final. In the general case, the new VOF interface reconstruction has 3rd order accuracy and can be easily made seamless. The new method addresses intrinsic issues with Young’s method such as gaps between interface facets in the case of a curved interface, and inability to define curvature nor identify corners. (2) A non-diffusive volume advection scheme. In an ALE advection step, a well-defined interface can be carried over through a Lagrange step and used to compute volume distribution into a relaxed mesh. Then, an interface reconstruction step is performed to redefine the interface in the relaxed mesh. We must point out that the interface carried over is also a solution of interface reconstruction because all the volume fractions in the relaxed mesh are naturally matched. We provide an interface tracking method compatible with our reconstruction scheme, where it is granted to use the prior info as an initial guess to capture sub-mesh resolution features. As a result, we are able to treat multiple facets inside a single mixed cell and obtain highly accurate, non-diffusive solution for advection problems with rather coarse meshes. We show our solutions for two-dimensional incompressible flows with two materials with a) the X + O diagonal translation; b) the Zalesak rotational test; and c) the single vortex spiral test.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Magnetic-Field Diffusion Effects in Beam Position Monitors III: Application to DARHT-II Beam Data

Beam position monitors (BPMs) provide time-resolved measurements of the current and centroid position of high-current electron beams in linear induction accelerators (LIAs). One of the types of detectors used in BPMs is the B-dot loop, which generates a signal from the EMF due to the time varying magnetic flux through the loop. If some of the boundaries of the loop are composed of thick metal walls with finite conductivity, the resulting signal must be corrected for the magnetic field diffusion into the metal. The theoretically predicted flux due to diffusion is in remarkable agreement with experimental measurements. Although accurate BPM measurements of beam current require correction of magnetic field diffusion, accurate measurement of beam position requires no correction. In this note, we present experimental validation of current and position results from a prototype detector employing finite conductivity sensing areas, based on experiments on the DARHT-II LIA.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

First-Principles Studies of Tritium Species Diffusivity Across the Interfaces of Ni-NiZr Alloy-Zircaloy-4

A tritium producing assembly, also known as tritium producing burnable absorber rods (TPBARs), consists of a metal getter tube located between the cladding and γ-LiAlO2 pellets. The metal getter tube is composed of a nickel (Ni) layer coated on Zircaloy-4. The getter assembly is used to capture tritium (3H) species (mainly 3H2 and 3H2O) generated from γ-LiAlO2 pellets during irradiation. Exploring 3H species (3H2, 3H2O) dissociation on the Ni surface and diffusion in the Ni layer and across the interface of Ni-plated Zircaloy-4 getters can provide insights on tritium transport and retention in the pellets and the getter materials. . In FY25, based on the ideal interface model of Ni-Zircaloy-4 generated in FY24, we will further explore the diffusion pathways of 3H across such Ni-ZrNi-alloy layer-Zircalory-4 interfaces under different conditions, including with oxide/hydroxide clusters on the Ni layer, and impurities located in ZrNi-alloy and Zircalory-4 layers. Due to the size limitation of DFT simulation, we will separate the 3-layer system into 3 sub-systems: Ni-ZrNi interface, ZrNi alloy layer, ZrNi-alloy-Zircaloy-4 interface, and simulate the diffusion barriers for tritium.

diffusion barrier↗
DOE OSTI2023

Diffusion Probabilistic Modeling for Video Generation

Denoising diffusion probabilistic models are a promising new class of generative models that mark a milestone in high-quality image generation. This paper showcases their ability to sequentially generate video, surpassing prior methods in perceptual and probabilistic forecasting metrics. We propose an autoregressive, end-to-end optimized video diffusion model inspired by recent advances in neural video compression. The model successively generates future frames by correcting a deterministic next-frame prediction using a stochastic residual generated by an inverse diffusion process. We compare this approach against six baselines on four datasets involving natural and simulation-based videos. We find significant improvements in terms of perceptual quality and probabilistic frame forecasting ability for all datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

ZrN Phase Formation, Hardening and Nitrogen Diffusion Kinetics in Plasma Nitrided Zircaloy-4

Plasma nitridation was conducted to modify the surfaces of Zircaloy-4. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), and Raman analysis were used to characterize microstructures and phases. Surface indentation and cross-sectional indentation were performed to evaluate mechanical property changes. Nitridation forms a thin layer of ZrN phase, followed by a much deeper layer affected by nitrogen diffusion. The ZrN phase is confirmed by both TEM and Raman characterization. The Raman peaks of ZrN phase show a temperature dependence. The intensity increases with increasing nitridation temperatures, reaches a maximum at 700 °C, and then decreases at higher temperatures. The ZrN layer appears as continuous small columnar grains. The surface polycrystalline ZrN phase is harder than the bulk by a factor of ~8, and the nitrogen diffusion layer is harder by a factor of ~2–5. The activation energy of nitrogen diffusion was measured to be 2.88 eV. The thickness of the nitrogen-hardened layer is controllable by changing the nitridation temperature and duration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulations of neutron noise in the research reactor AKR-2: comparison between a discrete ordinates and a diffusion-based method

A diffusion-based and a discrete ordinates method are used to simulate a neutron noise experiment in the research reactor AKR-2 at the Technical University in Dresden, Germany. The AKR-2 reactor provides an interesting case for the comparison between the two methods because it is characterized by large heterogeneities and regions with low macroscopic neutron cross-sections. For the calculations, the same spatial discretization and the same set of two-energy macroscopic neutron cross-sections with isotropic scattering are used. Significant discrepancies between the diffusion-based and discrete ordinates methods are found in regions of the systems where the diffusion approximation is expected to be inaccurate in reproducing characteristics of the static neutron flux and neutron noise. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

ELECTRIC FIELD-ASSISTED DIFFUSION WELDING TO FABRICATE ALLOY 617 COMPACT HEAT EXCHANGERS

Compact heat exchangers are of interest for a number of applications including advanced reactors. Alloy 617 is one of the top candidate materials for the gas-cooled reactor intermediate heat exchanger. Previous endeavors to diffusion weld Alloy 617 utilized hot pressing (HP). It was reported that grain boundary migration across the interface was hindered by extensive precipitation. Bonds of this nature have been observed to reduce the elevated-temperature mechanical properties compared to the wrought-product form. It was hypothesized that the electric current applied during electric-field-assisted sintering (EFAS) can overcome these challenges, resulting in improved diffusion welding (DW). This study investigated DW of Alloy 617 via EFAS. Stacks composed of three sheets that were 20 mm in diameter were welded using EFAS. Specimens were welded with an applied electric current, a pressure of 30 MPa, hold time of 30 min, and temperatures of 1050°C, 1100°C, and 1150°C. DW using HP as the zero-current analog of EFAS was also performed at the most promising EFAS conditions. Results revealed that both the applied electric current and temperature played a key role in precipitation and grain boundary migration in diffusion-welded Alloy 617. Precipitates were observed at the interface of the hot-pressed samples which limited grain boundary migration. Electric current was found to prevent precipitate formation along the interface at 1150°C. The electric current coupled with a temperature of 1150°C during EFAS resulted in significant grain boundary migration across the interface.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

A Guideline to Mitigate Interfacial Degradation Processes in Solid‐State Batteries Caused by Cross Diffusion

Abstract Diffusion of transition metals across the cathode–electrolyte interface is identified as a key challenge for the practical realization of solid‐state batteries. This is related to the formation of highly resistive interphases impeding the charge transport across the materials. Herein, the hypothesis that formation of interphases is associated with the incorporation of Co into the Li 7 La 3 Zr 2 O 12 lattice representing the starting point of a cascade of degradation processes is investigated. It is shown that Co incorporates into the garnet structure preferably four‐fold coordinated as Co 2+ or Co 3+ depending on oxygen fugacity. The solubility limit of Co is determined to be around 0.16 per formula unit, whereby concentrations beyond this limit causes a cubic‐to‐tetragonal phase transition. Moreover, the temperature‐dependent Co diffusion coefficient is determined, for example, D 700 °C = 9.46 × 10 −14 cm 2 s −1 and an activation energy E a = 1.65 eV, suggesting that detrimental cross diffusion will take place at any relevant process condition. Additionally, the optimal protective Al 2 O 3 coating thickness for relevant temperatures is studied, which allows to create a process diagram to mitigate any degradation with a minimum compromise on electrochemical performance. This study provides a tool to optimize processing conditions toward developing high energy density solid‐state batteries.

Chemistry↗