Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Communication: Diverse nanoscale cluster dynamics: Diffusion of 2D epitaxial clusters

The dynamics of nanoscale clusters can be distinct from macroscale behavior described by continuum formalisms. For diffusion of 2D clusters of N atoms in homoepitaxial systems mediated by edge atom hopping, macroscale theory predicts simple monotonic size scaling of the diffusion coefficient, D N ~ N –β , with β = 3/2. However, modeling for nanoclusters on metal(100) surfaces reveals that slow nucleation-mediated diffusion displaying weak size scaling β < 1 occurs for “perfect” sizes N p = L 2 and L(L+1) for integer L = 3,4,… (with unique square or near-square ground state shapes), and also for N p +3, N p +4,…. In contrast, fast facile nucleation-free diffusion displaying strong size scaling β ≈ 2.5 occurs for sizes N p +1 and N p +2. D N versus N oscillates strongly between the slowest branch (for N p +3) and the fastest branch (for N p +1). All branches merge for N = O(10 2 ), but macroscale behavior is only achieved for much larger N = O(10 3 ). Here, this analysis reveals the unprecedented diversity of behavior on the nanoscale.

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

A diffusion-limited reaction model for self-propagating Al/Pt multilayers with quench limits

A diffusion-limited reaction model was calibrated for Al/Pt multilayers ignited on oxidized silicon, sapphire, and tungsten substrates, as well as for some Al/Pt multilayers ignited as free-standing foils. The model was implemented in a finite element analysis code and used to match experimental burn front velocity data collected from several years of testing at Sandia National Laboratories. Moreover, both the simulations and experiments reveal well-defined quench limits in the total Al + Pt layer (i.e., bilayer) thickness. At these limits, the heat generated from atomic diffusion is insufficient to support a self-propagating wave front on top of the substrates. Quench limits for reactive multilayers are seldom reported and are found to depend on the thermal properties of the individual layers. Here, the diffusion-limited reaction model is generalized to allow for temperature- and composition-dependent material properties, phase change, and anisotropic thermal conductivity. Utilizing this increase in model fidelity, excellent overall agreement is shown between the simulations and experimental results with a single calibrated parameter set. However, the burn front velocities of Al/Pt multilayers ignited on tungsten substrates are over-predicted. Finally, possible sources of error are discussed and a higher activation energy (from 41.9 kJ/mol.at. to 47.5 kJ/mol.at.) is shown to bring the simulations into agreement with the velocity data observed on tungsten substrates. Finally, this higher activation energy suggests an inhibited diffusion mechanism present at lower heating rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Comprehensive evaluation and parametric sensitivity of interatomic potential models for diffusion kinetics of Cr 2 O 3 in molecular dynamics

While molecular dynamics (MD) has proven to be a promising approach to investigate the diffusion properties, the grand challenge resides in evaluating potential model parameters to accurately replicate experimentally measured properties. The Buckingham potential model with Columbic interaction is widely employed in MD simulations of chromia (Cr 2 O 3 ) systems, as it allows for reasonable computational cost and accuracy. However, considering the well-known limitation of classical potential models in simultaneous reproduction of various physical phenomena, further comprehensive evaluation of the potential is required for calculation of diffusion properties. In this study, we benchmark the performance of three different Buckingham models with the experimental data by calculating structural, thermodynamic, and mechanical properties of defect-free Cr 2 O 3 , and diffusion properties of Cr 2 O 3 with vacancy defects. Available Buckingham models display limited accuracies, consolidating the necessity of retraining the potential parameters for all properties impacting the diffusion dynamics. Oversimplification in parameterization procedures is suggested to impede the universal performance in property reproduction. Furthermore, this research also demonstrates effective guidelines for choosing a proper parameter set of current Buckingham potential for MD simulation with Cr 2 O 3 depending on properties and for potential reparameterization.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Solid-phase reactive inter-diffusion of Mg/B multilayers

A better understanding of the reaction of Mg and B in the solid-phase regime is needed for the low-temperature synthesis of MgB 2 films. Here in this paper, we study the kinetics of reactive inter-diffusion of Mg and B multilayers on glassy carbon substrates in the temperature range of 400-650°C. Results show that, at these temperatures, inter-diffusion is characterized by a single activation energy of ~0.45eV. The formation of the superconducting MgB 2 phase with critical temperatures of 25–31 K occurs at reaction temperatures of 450°C and above, with the rate of inter-diffusion obeying a power law with a kinetic exponent of ~0.3. This suggests that rate-limiting processes are the nucleation and growth of MgB 2 grains rather than diffusion and interfacial reactions. Implications of these results to the low-temperature synthesis of MgB 2 films are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The impact of nanoparticle softness on its tracer diffusion coefficient in all polymer nanocomposites

The diffusion of nanoparticles in a polymer matrix is an area of current interest. However, a complete understanding is still limited as it is often difficult to quantify the much slower motion of nanoparticles in a polymer matrix. To combat this problem, we have developed a protocol to measure the diffusion coefficient of soft nanoparticles in a linear polymer matrix. Recently developed synthetic control over soft nanoparticle structures combined with this protocol provides a pathway to separately elucidate the effects of the molecular weight and nanoparticle softness on its diffusive behavior. These results indicate that the nanoparticle softness and deformability dictate its motion. Increasing the cross-linking density of the nanoparticle for all molecular weights increases its hardness and suppresses its motion in the linear matrix. Additionally, the nanoparticle molecular weight dependence deviates from the exponential dependence for star polymers suggesting that these nanoparticles benefit from the cooperative motion of the matrix to open pathways for the nanoparticle. Finally, comparison of these experimentally determined values to the Stokes–Einstein theory demonstrates that the nanoparticles diffuse much slower than a hard sphere. This is interpreted to indicate that there exist additional interactions between the nanoparticle and polymer matrix that are not captured by Stokes–Einstein, including threading or entanglement of the linear chain with the nanoparticle.

Rostom, Sahar↗
DOE OSTI2023

Diffusion of Sn donors in β-Ga 2 O 3

Diffusion of the n-type dopant Sn in β-Ga 2 O 3 is studied using secondary-ion mass spectrometry combined with hybrid functional calculations. The diffusion of Sn from a Sn-doped bulk substrate with surface orientation (001) into an epitaxial layer is observed after heat treatments in the temperature range of 1050–1250 °C. Calculated formation energies of Sn-related and intrinsic defects show that the migration of Sn is mediated by Ga vacancies ( V Ga ) through the formation and dissociation of intermittent mobile V Ga Sn Ga complexes. The evolution of the Sn concentration vs depth profiles after heat treatments can be well described by a reaction–diffusion model. Using model parameters guided by the hybrid functional calculations, we extract a V Ga Sn Ga complex migration barrier of 3.0 ± 0.4 eV with a diffusion coefficient of 2 × 10 –1 cm 2 /s. The extracted migration barrier is consistent with our theoretical predictions using the nudged elastic band method, which shows migration barriers of 3.42, 3.15, and 3.37 eV for the [100], [010], and [001] directions, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models

Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.

17 WIND ENERGY↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2019

A Multilevel in Space and Energy Solver for 3-D Multigroup Diffusion and Coarse-Mesh Finite Difference Eigenvalue Problems

The Multilevel in Space and Energy Diffusion (MSED) method accelerates the iterative convergence of multigroup diffusion eigenvalue problems by performing work on lower-order equations with only one group and/or coarser spatial grids. It consists of two primary components: (1) a grey (one-group) diffusion eigenvalue problem that is solved via Wielandt-shifted power iteration (PI) and (2) a multigrid-in-space linear solver. In previous work, the efficiency of MSED was verified using Fourier analysis and numerical results from a one-dimensional multigroup diffusion code. Since that work, MSED has been implemented as a solver for the coarse-mesh finite difference (CMFD) system in the three-dimensional Michigan Parallel Characteristics Transport (MPACT) code. In this paper, the results from the implementation of MSED in MPACT are introduced, and the changes needed to make MSED more suitable for MPACT are described. For problems without feedback, the conclusions in this paper show that MSED can reduce the CMFD run time by an order of magnitude and the overall run time by a factor of 2 to 3 compared to the default CMFD solver in MPACT [PI with the generalized minimal residual (GMRES) method]. For problems with feedback, the convergence of the outer Picard iteration scheme is worsened by the well-converged CMFD solutions produced by the standard MSED method. To overcome this unintuitive deficiency, MSED may be run with looser convergence criteria (a modified version of the MSED method called MSED-L) to circumvent the issue until the multiphysics iteration in MPACT is improved. Results show that MSED-L can reduce the CMFD run time in MPACT by an order of magnitude, without negatively impacting the outer Picard iteration scheme.

3-D↗
DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Search for a diffuse flux of photons with energies above tens of PeV at the Pierre Auger Observatory

Diffuse photons of energy above 0.1 PeV, produced through the interactions between cosmic rays and either interstellar matter or background radiation fields, are powerful tracers of the distribution of cosmic rays in the Galaxy. Furthermore, the measurement of a diffuse photon flux would be an important probe to test models of super-heavy dark matter decaying into gamma-rays. In this work, we search for a diffuse photon flux in the energy range between 50 PeV and 200 PeV using data from the Pierre Auger Observatory. For the first time, we combine the air-shower measurements from a 2 km 2 surface array consisting of 19 water-Cherenkov surface detectors, spaced at 433 m, with the muon measurements from an array of buried scintillators placed in the same area. Using 15 months of data, collected while the array was still under construction, we derive upper limits to the integral photon flux ranging from 13.3 to 13.8 km -2 sr -1 yr -1 above tens of PeV. We extend the Pierre Auger Observatory photon search program towards lower energies, covering more than three decades of cosmic-ray energy. This work lays the foundation for future diffuse photon searches: with the data from the next 10 years of operation of the Observatory, this limit is expected to improve by a factor of ∼20.

cosmic rays detectors↗
DOE OSTI2026

Effect of controlled magnetic island bifurcation on electron diffusion

Magnetic islands strongly influence cross-field electron transport in magnetized plasmas. In particular, bifurcations of the island topology modify the number and location of O-points, X-points, and separatrix boundaries, thereby altering diffusion pathways. In recent DIII-D experiments, external magnetic perturbations were used to rotate and periodically bifurcate the island on the q = 2 surface, causing a switchback between a q = 2/1-dominated structure and a narrower q = 4/2-dominated structure. To investigate how this topological change affects electron transport, we employ the field line tracing code TRIP3D with an implemented collisional operator. Thermal, tracer electrons launched from O-points, X-points, and outside separatrix boundaries reveal distinct diffusion regimes, including classical, subdiffusive, and superdiffusive behavior, depending on both the dominant island mode and launch location. These results suggest that island bifurcation can alter electron diffusion across rational surfaces, with direct implications for particle confinement. While the present work emphasizes diffusion as a general framework, the findings provide insight into the conditions under which electron trapping into an island or stochastization of the island's separatrix can enable additional mechanisms, such as the generation of energetic electrons.

DIII-D↗
DOE OSTI2017

Analysis of non-diffuse characteristics of the seismic noise field in southern California based on correlations of neighbouring frequencies

Non-diffuse characteristics of the ambient seismic noise wavefield recorded at 154 broad-band stations in southern California are analysed by computing the correlation matrix of power spectral values at neighbouring frequencies. The similarities of the derived correlation matrices are compared by choosing a different reference station at a time. The data recorded by the different stations are classified into five groups using a hierarchical clustering algorithm based on the similarity of the correlation matrices. Stations belonging to different groups are clustered spatially in the Los Angeles basin, mountains and regions farther from the coast. The deviations from diffuse wavefield in the correlation matrices of representative stations from the five groups are inverted for dominant cross-frequency components. These are used to derive power spectral ratios of non-diffuse noise and fully diffuse noise. The similarity maps for different reference stations and the inverted cross-frequency components provide information on properties of the noise sources and propagation/scattering characteristics in the different regions.

58 GEOSCIENCES↗
DOE OSTI2026

Observation of the Jet Diffusion Wake Using Dijets in Heavy-Ion Collisions

Energetic quarks and gluons traversing a hot and dense quark-gluon plasma deposit energy and momentum into the medium before hadronizing to collimated sprays of particles, known as jets. This energy-momentum deposition is expected to produce medium responses, collectively known as jet wakes, with “diffusion wake” denoting a depletion of particles in the direction opposite to the propagating jet. The diffusion wake is studied by comparing dijet-hadron correlations measured in lead-lead (PbPb) and proton-proton (pp) collisions. The analysis uses PbPb and pp data recorded at a nucleon-nucleon (NN) center-of-mass energy s NN = 5.02 TeV with the CMS detector at the CERN LHC. By exploring how the dijet-hadron correlation distributions differ for various pseudorapidity separations of the two jets in the dijet, the presence of a jet diffusion wake is firmly established. The diffusion wake signal has a significance greater than 5 standard deviations for charged particles in the transverse momentum range 1 < p T < 2 GeV . The measurements are compared with model predictions with and without jet wake effects, providing new insights into quark-gluon plasma properties and the formation of jet-induced wakes.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2017

Electric-field noise from carbon-adatom diffusion on a Au(110) surface: First-principles calculations and experiments

The decoherence of trapped-ion quantum bits due to heating of their motional modes is a fundamental science and engineering problem. This heating is attributed to electric-field noise arising from processes on the trap-electrode surfaces. In this work, we address the source of this noise by focusing on the diffusion of carbon-containing adsorbates on the surface of Au(110). We show by detailed scanned probe microscopy and density functional theory how the carbon adatom diffusion on the gold surface changes the energy landscape, and how the adatom dipole moment varies with the diffusive motion. Lastly, a simple model for the diffusion noise, which varies quadratically with the variation of the dipole moment, qualitatively reproduces the measured noise spectrum, and the estimate of the noise spectral density is in accord with measured values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Energy barriers for Pb adatom diffusion on stepped ultrathin Pb(111) quantum nanofilms: First-principles calculations

Various properties of Pb(111) nanofilms, prototypical quantum films, have been studied extensively. However, key ab initio -level energy barriers for Pb adatom diffusion on stepped Pb(111) nanofilm surfaces are still not available. Using first-principles density functional theory, we calculate these barriers for films with thicknesses of few monolayers (ML). We find that two-atom exchange is always much more favorable than direct hopping to cross A- or B-type steps. Ehrlich-Schwoebel (ES) barriers for downward transport to a higher-coordination step-edge site depend strongly on the film thickness. For such transport from 2- to 1-ML terraces, or from 4- to 3-ML terraces, there is no an ES barrier, but large ES barriers of more than 100 meV are found from 3- to 2-ML terraces. We also obtain the barriers for diffusion along the step edges and find that these step-edge barriers are significantly larger than terrace diffusion barriers. In addition, we analyze energetics for diffusion on the top flat surface of a nanofilm supported on a vicinal surface, and thus having underlying buried or ghost steps. We quantify the tilted potential energy surface in both ghost A- and B-step regions separating 2- and 3-ML (as well as 3- and 4-ML) terraces. Finally, consequences are discussed for the growth kinetics of supported Pb nanofilms, where the support does not strongly affect electronic states within the film.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Thermal diffusivity above the Mott-Ioffe-Regel limit

We present high-resolution thermal diffusivity measurements on several near optimally doped electron- and hole-doped cuprate systems in a temperature range that passes through the Mott-Ioffe-Regel limit, above which the quasiparticle picture fails. Our primary observations are that the inverse thermal diffusivity is linear in temperature and can be fitted to $D^{-1}_Q = aT + b$. The slope a is interpreted through the Planckian relaxation time $τ ≈ ℏ/k_BT$ and a thermal diffusion velocity v B, which is close, but larger than the sound velocity. The intercept b represents a crossover diffusion constant that separates coherent from incoherent quasiparticles. These observations suggest that both phonons and electrons participate in the thermal transport, while reaching the Planckian limit for relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗