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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Diffuse auroral precipitation driven by lower-band chorus second harmonics

Diffuse aurora at the Earth’s high latitude regions is mainly caused by the low-energy (0.1–30 keV) electron precipitation which carries the major energy flux into the nightside upper atmosphere. Previous studies have demonstrated that combined scattering by the upper- and lower- band chorus waves acts as the dominant cause of diffuse auroral precipitation, but that is not necessarily the case as these two types of waves do not always occur simultaneously, with the lower-band more often. Here we report that the lower-band chorus satisfying the preferred condition can generate their second harmonics so as to trigger the diffuse auroral electron precipitation. We find that the lower-band chorus alone can only cause the precipitation of electrons greater than 4 keV, while the self-consistently generated second harmonic is weak but still able to result in the electron precipitation below 4 keV. The combined effect of those modes results in the observed pancake electron distributions and the diffuse aurora. Our results clearly demonstrate an alternative but universal mechanism of chorus-driven diffuse aurora in the Earth, which may also apply to the auroral formation in other planetary magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials

Diffusion involving atom transport from one location to another governs many important processes and behaviors such as precipitation and phase nucleation. The inherent chemical complexity in compositionally complex materials poses challenges for modeling atomic diffusion and the resulting formation of chemically ordered structures. Here, we introduce a neural network kinetics (NNK) scheme that predicts and simulates diffusion-induced chemical and structural evolution in complex concentrated chemical environments. The framework is grounded on efficient on-lattice structure and chemistry representation combined with artificial neural networks, enabling precise prediction of all path-dependent migration barriers and individual atom jumps. To demonstrate the method, we study the temperature-dependent local chemical ordering in a refractory NbMoTa alloy and reveal a critical temperature at which the B2 order reaches a maximum. The atomic jump randomness map exhibits the highest diffusion heterogeneity (multiplicity) in the vicinity of this characteristic temperature, which is closely related to chemical ordering and B2 structure formation. The scalable NNK framework provides a promising new avenue to exploring diffusion-related properties in the vast compositional space within which extraordinary properties are hidden.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Bifurcation and Pattern Symmetry Selection in Reaction-Diffusion Systems with Kinetic Anisotropy

Abstract Ordering and self-organization are critical in determining the dynamics of reaction-diffusion systems. Here we show a unique pattern formation mechanism, dictated by the coupling of thermodynamic instability and kinetic anisotropy. Intrinsically different from the physical origin of Turing instability and patterning, the ordered patterns we obtained are caused by the interplay of the instability from uphill diffusion, the symmetry breaking from anisotropic diffusion, and the reactions. To understand the formation of the void/gas bubble superlattices in crystals under irradiation, we establish a general theoretical framework to predict the symmetry selection of superlattice structures associated with anisotropic diffusion. Through analytical study and phase field simulations, we found that the symmetry of a superlattice is determined by the coupling of diffusion anisotropy and the reaction rate, which indicates a new type of bifurcation phenomenon. Our discovery suggests a means for designing target experiments to tailor different microstructural patterns.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of Sr doping and external strain on relieving bottleneck of oxygen diffusion in La 2- x Sr x CuO 4-δ

In many complex oxides, the oxygen vacancy formation is a promising route to modify the material properties such as a superconductivity and an oxygen diffusivity. Cation substitutions and external strain have been utilized to control the concentration and diffusion of oxygen vacancies, but the mechanisms behind the controls are not fully understood. Using first-principles calculations, we find how Sr doping and external strain greatly enhances the diffusivity of oxygen vacancies in La 2-x Sr x CuO 4-δ (LSCO) in the atomic level. In hole-doped case (2x > δ), the formation energy of an apical vacancy in the LaO layer is larger than its equatorial counterpart by 0.2 eV that the bottleneck of diffusion process is for oxygen vacancies to escape equatorial sites. Such an energy difference can be reduced and even reversed by either small strain (< 1.5%) or short-range attraction between Sr and oxygen vacancy, and in turn, the oxygen diffusivity is greatly enhanced. For fully compensated hole case (2x ≦ δ), the formation energy of an apical vacancy becomes too high that most oxygen vacancies cannot move but would be trapped at equatorial sites. From our electronic structure analysis, we found that the contrasting change in the formation energy by Sr doping and external strain is originated from the different localization natures of electron carrier from both types of oxygen vacancies.

25 ENERGY STORAGE↗
DOE OSTI2017

Dislocation-pipe diffusion in nitride superlattices observed in direct atomic resolution

Device failure from diffusion short circuits in microelectronic components occurs via thermally induced migration of atoms along high-diffusivity paths: dislocations, grain boundaries, and free surfaces. Even well-annealed single-grain metallic films contain dislocation densities of about 10 14 m -2 ; hence dislocation-pipe diffusion (DPD) becomes a major contribution at working temperatures. While its theoretical concept was established already in the 1950s and its contribution is commonly measured using indirect tracer, spectroscopy, or electrical methods, no direct observation of DPD at the atomic level has been reported. We present atomically-resolved electron microscopy images of the onset and progression of diffusion along threading dislocations in sequentially annealed nitride metal/semiconductor superlattices, and show that this type of diffusion can be independent of concentration gradients in the system but governed by the reduction of strain fields in the lattice.

42 ENGINEERING↗
DOE OSTI2019

Ion diffusion across a disorder-to-order phase transition in a poly(ethylene oxide)- b -poly(silsesquioxane) block copolymer electrolyte

Nanostructured block copolymer electrolytes composed of organic and inorganic moieties have the potential to enable solid-state batteries. Practical uses of these materials, however, require an understanding of the microscopic and macroscopic ion transport properties across the microphase-separated systems. The self-diffusion of salt ions across a disorder-lamellar phase transition in a nanostructured poly(ethylene oxide)-b-poly(silsesquioxane) copolymer was studied using pulsed-field gradient NMR (PFG-NMR) and changes in the morphology were studied using small-angle X-ray scattering. The diffusion of the salt is isotropic when the polymer electrolyte is disordered and locally anisotropic when the polymer is microphase separated. The difference between the diffusion coefficient parallel to the lamellae, D ∥ , and the diffusion coefficient perpendicular to the lamellae, D ⊥ , measured using PFG-NMR, increases with increasing segregation strength. Finally, the anisotropy of diffusion parallels changes in the morphology measured by small-angle X-ray scattering.

42 ENGINEERING↗
DOE OSTI2021

Fast Na diffusion and anharmonic phonon dynamics in superionic Na 3 PS 4

The design of new solid electrolytes (SEs) hinges on identifying and tuning relevant descriptors. Phonons describe the atomic dynamics in crystalline materials and provide a basis to encode possible minimum energy pathways for ion migration, but anharmonic effects can be large in SEs. Identifying and controlling the pertinent phonon modes coupled most strongly with ionic conductivity, and assessing the role of anharmonicity, could therefore pave the way for discovering and designing new SEs via phonon engineering. In this work, we investigate phonons in Na 3 PS 4 and their coupling to fast Na diffusion, using a combination of neutron scattering, ab initio molecular dynamics (AIMD), and extended molecular dynamics based on machine-learned potentials. We identify that anharmonic soft-modes at the Brillouin zone boundary of the anharmonically stabilized cubic phase constitute key phonon modes that control the Na diffusion process in Na 3 PS 4 . We demonstrate how these strongly anharmonic phonon modes enable Na-ions to hop along the minimum energy pathways. Further, the quasi-elastic neutron scattering (QENS) measurements, supplemented with ab initio and machine-learned molecular dynamics simulations, probe the Na diffusivity and diffusion mechanism. These results offer detailed microscopic insights into the dynamic mechanism of fast Na diffusion and provide an avenue to search for further Na solid electrolytes.

42 ENGINEERING↗
DOE OSTI2021

Stoichiometric tuning of lattice flexibility and Na diffusion in NaAlSiO 4 : quasielastic neutron scattering experiment and ab initio molecular dynamics simulations

We have performed quasielastic neutron scattering (QENS) experiments up to 1243 K and ab initio molecular dynamics (AIMD) simulations to investigate the Na diffusion in various phases of NaAlSiO4 (NASO), namely, low-carnegieite (L-NASO; trigonal), high-carnegieite (H-NASO; cubic) and nepheline (N-NASO; hexagonal) phases. The QENS measurements reveal Na ions localized diffusion behavior in L-NASO and N-NASO, but long-range diffusion behavior in H-NASO. Furthermore, the AIMD simulation supplemented the QENS measurements and showed that excess Na ions in H-NASO enhance the host network flexibility and activate the AlO 4 /SiO 4 tetrahedra rotational modes. These framework modes enable the long-range diffusion of Na across a pathway of interstitial sites. The simulations also show Na diffusion in Na-deficient N-NASO through vacant Na sites along the hexagonal c-axis.

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

Optical/X-ray/radio view of Abell 1213: A galaxy cluster with anomalous diffuse radio emission

Context. Abell 1213, a low-richness galaxy system, is known to host an anomalous radio halo detected in data of the Very Large Array (VLA). It is an outlier with regard to the relation between the radio halo power and the X-ray luminosity of the parent clusters. Aims. Our aim is to analyze the cluster in the optical, X-ray, and radio bands to characterize the environment of its diffuse radio emission and to shed new light on its nature. Methods. We used optical data from the Sloan Digital Sky Survey to study the internal dynamics of the cluster. We also analyzed archival XMM-Newton X-ray data to unveil the properties of its hot intracluster medium. Finally, we used recent data from the LOw Frequency ARray (LOFAR) at 144 MHz, together with VLA data at 1.4 GHz, to study the spectral behavior of the diffuse radio source. Results. Both our optical and X-ray analysis reveal that this low-mass cluster exhibits disturbed dynamics. In fact, it is composed of several galaxy groups in the peripheral regions and, in particular, in the core, where we find evidence of substructures oriented in the NE–SW direction, with hints of a merger nearly along the line of sight. The analysis of the X-ray emission adds further evidence that the cluster is in an unrelaxed dynamical state. At radio wavelengths, the LOFAR data show that the diffuse emission is ~510 kpc in size. Moreover, there are hints of low-surface-brightness emission permeating the cluster center. Conclusions. The environment of the diffuse radio emission is not what we would expect for a classical halo. The spectral index map of the radio source is compatible with a relic interpretation, possibly due to a merger in the N–S or NE–SW directions, in agreement with the substructures detected through the optical analysis. The fragmented, diffuse radio emissions at the cluster center could be attributed to the surface brightness peaks of a faint central radio halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A horseshoe-shaped ring of diffuse emission detected at 1.4 GHz

We identify a horseshoe-shaped ring (HSR) of diffuse emission in J1407+0453 from the Faint Images of Radio Sky at Twenty-cm (FIRST) survey using the Very Large Array telescope at 1.4 GHz. An optical galaxy, SDSS J140709.01+045302.1, is present near the limb of the HSR of J1407+0453, with a spectroscopic redshift ofz = 0.13360. The total extent of the source, including the diffuse emission of J1407+0453, is 65 arcsec (with a physical extent of 160 kpc), whereas the diameter of the HSR is approximately 10 arcsec (25 kpc). The flux density of the HSR is ∼47 mJy at 1400 MHz, whereas the flux densities of the whole diffuse emission of J1407+0453 at 1400 MHz and 150 MHz are 172 mJy and 763 mJy, respectively. We measure the radio luminosity of HSR J1407+0453 as 1.94 × 10 24 W Hz −1 , with a spectral index,α 150 1400 = −0.67. The black hole mass of J1407+0453 is 5.8 × 10 8 M ⊙ . We compare the HSR of diffuse emission of J1407+0453 with other discovered diffuse circular sources. The possible formation scenarios for J1407+0453 are discussed, so as to understand the nature of the source. We present a spectral index map of source J147+0453 to study the spectral properties of the source.

Astronomy & Astrophysics↗
DOE OSTI2017

Reaction rates for reaction-diffusion kinetics on unstructured meshes

The reaction-diffusion master equation is a stochastic model often utilized in the study of biochemical reaction networks in living cells. It is applied when the spatial distribution of molecules is important to the dynamics of the system. A viable approach to resolve the complex geometry of cells accurately is to discretize space with an unstructured mesh. Diffusion is modeled as discrete jumps between nodes on the mesh, and the diffusion jump rates can be obtained through a discretization of the diffusion equation on the mesh. Reactions can occur when molecules occupy the same voxel. Here in this paper, we develop a method for computing accurate reaction rates between molecules occupying the same voxel in an unstructured mesh. For large voxels, these rates are known to be well approximated by the reaction rates derived by Collins and Kimball, but as the mesh is refined, no analytical expression for the rates exists. We reduce the problem of computing accurate reaction rates to a pure preprocessing step, depending only on the mesh and not on the model parameters, and we devise an efficient numerical scheme to estimate them to high accuracy. We show in several numerical examples that as we refine the mesh, the results obtained with the reaction-diffusion master equation approach those of a more fine-grained Smoluchowski particle-tracking model.

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

Quasilinear diffusion coefficients in a finite Larmor radius expansion for ion cyclotron heated plasmas

In this study, a reduced model of quasilinear velocity diffusion by a small Larmor radius approximation is derived to couple the Maxwell’s equations and the Fokker Planck equation self-consistently for the ion cyclotron range of frequency waves in a tokamak. The reduced model ensures the important properties of the full model by Kennel-Engelmann diffusion, such as diffusion directions, wave polarizations, and H-theorem. The kinetic energy change (Wdot ) is used to derive the reduced model diffusion coefficients for the fundamental damping (n = 1) and the second harmonic damping (n = 2) to the lowest order of the finite Larmor radius expansion. The quasilinear diffusion coefficients are implemented in a coupled code (TORIC-CQL3D) with the equivalent reduced model of the dielectric tensor. We also present the simulations of the ITER minority heating scenario, in which the reduced model is verified within the allowable errors from the full model results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

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

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2022

Evidence for superionic H 2 O and diffusive He–H 2 O at high temperature and high pressure

In this work, we present the evidence of superionic phase formed in H 2 O and, for the first time, diffusive H 2 O–He phase, based on time-resolved x-ray diffraction experiments performed on ramp-laser-heated samples in diamond anvil cells. The diffraction results signify a similar bcc-like structure of superionic H 2 O and diffusive He–H 2 O, while following different transition dynamics. Based on time and temperature evolution of the lattice parameter, the superionic H 2 O phase forms gradually in pure H 2 O over the temperature range of 1350–1400 K at 23 GPa, but the diffusive He–H 2 O phase forms abruptly at 1300 K at 26 GPa. We suggest that the faster dynamics and lower transition temperature in He–H 2 O are due to a larger diffusion coefficient of interstitial-filled He than that of more strongly bound H atoms. This conjecture is then consistent with He disordered diffusive phase predicted at lower temperatures, rather than H-disordered superionic phase in He–H 2 O.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Investigations on a fuzzy process: effect of diffusion on calibration and particle identification in Liquid Argon Time Projection Chambers

Ionization electron diffusion in Liquid Argon Time Projection Chambers (LArTPCs) has typically been considered at the detector design stage, but little attention has been given to its effects on calibration and particle identification. We use a GEANT4-based simulation to study how diffusion impacts these techniques, and give consideration to how this effect is simulated. We find that diffusion can cause a drift-dependent bias to both the median and Most Probable Value (MPV) of dQ/dx distributions. The bias is estimated to be ~2.5% (median) and ~5.0% (MPV) for typical maximum drift times in currently running LArTPCs before adding detector specific considerations such as electric field non-uniformities. This information indicates that these metrics should not be used for calibration without care, contrary to the conventional wisdom. The impact of diffusion on the ability of LArTPCs to separate muons and protons is small, and not expected to pose any problems in future detectors. Diffusion may however be a significant source of systematic uncertainty when separating particles of more similar masses (muons and pions, kaons and protons). Separation of such populations may be improved by implementation of a drift-time dependent particle identification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗