Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗
DOE OSTI2017

A NEW GENERATION OF PARSEC-COLIBRI STELLAR ISOCHRONES INCLUDING THE TP-AGB PHASE

We introduce a new generation of PARSEC–COLIBRI stellar isochrones that includes a detailed treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase, covering a wide range of initial metallicities (0.0001 < Z {sub i} < 0.06). Compared to previous releases, the main novelties and improvements are use of new TP-AGB tracks and related atmosphere models and spectra for M and C-type stars; inclusion of the surface H+He+CNO abundances in the isochrone tables, accounting for the effects of diffusion, dredge-up episodes and hot-bottom burning; inclusion of complete thermal pulse cycles, with a complete description of the in-cycle changes in the stellar parameters; new pulsation models to describe the long-period variability in the fundamental and first-overtone modes; and new dust models that follow the growth of the grains during the AGB evolution, in combination with radiative transfer calculations for the reprocessing of the photospheric emission. Overall, these improvements are expected to lead to a more consistent and detailed description of properties of TP-AGB stars expected in resolved stellar populations, especially in regard to their mean photometric properties from optical to mid-infrared wavelengths. We illustrate the expected numbers of TP-AGB stars of different types in stellar populations covering a wide range of ages and initial metallicities, providing further details on the “C-star island” that appears at intermediate values of age and metallicity, and about the AGB-boosting effect that occurs at ages close to 1.6-Gyr for populations of all metallicities. The isochrones are available through a new dedicated web server.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol -1 relative to the experiment. Here, the two-point QZ/5Z and TZ/QZ CBS 1/$l^3_{max}$ extrapolation schemes are inferior to their inverse power counterpart (1/($l_{max}$ + 1/2) 4 ) by 1.3 and 4.3 kJ mol -1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol -1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double- ζ/triple- ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

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

Density functional theory study of the influence of segregation of S or Fe impurities on electromigration in nano-grained copper interconnects

The reliability of the circuits in the modern microelectronic devices remains, during last decades, one of the key topics in research and gains an attention for improving the promising candidates for conductors. Improvement of materials for such applications can be obtained by both electronic and compositional optimization. Ab initio calculations using full potential linearized augmented plane wave method in density functional theory are applied to explain the reduction in electromigration effect in the vicinity of grain boundaries (GB) in nano-structured Cu due to the segregation of some additives to the GB. Several possible mechanisms are considered. It is demonstrated that S atoms segregated to GB of nano-structured Cu lead to the growth of effective mass of the electrons. This decreases the mobility of electrons in external electric field, and, correspondingly, the momentum that they may transfer to atoms in collisions. Fe atoms segregated to GB of Cu create new empty states at the top of the valance band. These non-occupied states may stimulate the current of holes when external electric field is applied to the system, creating the “hole wind” in the direction opposite to the current of electrons. Such “hole wind” will compensate the forces generated by the electron current and therefore will reduce the total momentum transfer between charge carriers and atoms. The calculated electron density maps show that S and Fe segregating to Cu GB increases the strength of covalent bonds reducing the diffusion of Cu atoms in the vicinity of GB.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

The Mechanism of Influence of Disperse Nanoparticles on Parameters of the Martensitic Transitions in Alloys with the Shape Memory Effect

Within the diffuse martensitic transition theory based on the thermodynamic and kinetic equations and relations, the mechanism of influence of disperse nanoparticles on the parameters of martensitic transitions in alloys with the shape memory effect (SME) is analyzed. The objects of the analysis are the TiNi alloy with the varied sizes of Ti{sub 3}Ni{sub 4} particle at their constant volume concentration, and the NiMnGaTb alloy with Tb precipitate particles with constant sizes and variation of volume concentration of the precipitates. Information about these alloys is available in the literature. The analysis has shown that, due to the coherent coupling of the Ti{sub 3}Ni{sub 4} particles with the substrate, the temperature width of the transition R–B19' depends on the particle sizes, which confirms the earlier established regularity of the local interior strain influence on this parameter. Concerning the NiMnGaTb alloy, the analysis has shown that, due to the presence of interior local strains associated with the Tb particles, the temperature width of the martensitic transformations increases linearly alongside the growth of the particle concentration in the alloy. The existence of the critical value of the particle concentration, above which the temperature width of the transition becomes indefinitely large, and the martensitic transformation in the alloy is blocked, is shown.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Spatiotemporal Reaction Dynamics Control in Two‐Photon Polymerization for Enhancing Writing Characteristics

Since 2001, 3D microfabrication based on two‐photon polymerization (TPP) has drawn extensive attention and interest in biology, optics, photonics, material science, and high‐energy physics. The in‐volume fabrication capability due to the threshold behavior of two‐photon absorption enables TPP higher flexibility compared with other nanofabrication techniques. However, as determined by the in‐volume fabrication feature as well as various reaction dynamics, the writing characteristics of TPP, such as throughput, accuracy, surface quality, and fabrication capability, are still limited. Herein, a comprehensive study is performed on the spatiotemporal behavior of reaction dynamics during TPP fabrication, mainly focusing on spatiotemporal characteristics of radical diffusion, photothermal effect, microscale mechanics, and voxel stacking process. Based on the study, a nonsequential fabrication method is established to simultaneously improve key writing characteristics of TPP and realize sharp features, high speeds, large overhang structure, and smooth surfaces. The method established in this work can be applied to improve the performance of functional devices for various fields.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A one-dimensional Fickian model to predict the Ga depth profiles in three-stage Cu(In,Ga)Se{sub 2}

We present a one-dimensional Fickian model that predicts the formation of a double Ga gradient during the fabrication of Cu(In,Ga)Se{sub 2} thin films by three-stage thermal co-evaporation. The model is based on chemical reaction equations, structural data, and effective Ga diffusivities. In the model, the Cu(In,Ga)Se{sub 2} surface is depleted from Ga during the deposition of Cu-Se in the second deposition stage, leading to an accumulation of Ga near the back contact. During the third deposition stage, where In-Ga-Se is deposited at the surface, the atomic fluxes within the growing layer are inverted. This results in the formation of a double Ga gradient within the Cu(In,Ga)Se{sub 2} layer and reproduces experimentally observed Ga distributions. The final shape of the Ga depth profile strongly depends on the temperatures, times and deposition rates used. The model is used to evaluate possible paths to flatten the marked Ga depth profile that is obtained when depositing at low substrate temperatures. We conclude that inserting Ga during the second deposition stage is an effective way to achieve this.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

A self-consistent and spatially dependent model of the multiband emission of pulsar wind nebulae

A self-consistent and spatially dependent model is presented to investigate the multiband emission of pulsar wind nebulae (PWNe). In this model, a spherically symmetric system is assumed and the dynamical evolution of the PWN is included. The processes of convection, diffusion, adiabatic loss, radiative loss, and photon–photon pair production are taken into account in the electron’s evolution equation, and the processes of synchrotron radiation, inverse Compton scattering, synchrotron self-absorption, and pair production are included for the photon’s evolution equation. Both coupled equations are simultaneously solved. The model is applied to explain observed results of the PWN in MSH 15–52. Our results show that the spectral energy distributions (SEDs) of both electrons and photons are all a function of distance. The observed photon SED of MSH 15–52 can be well reproduced in this model. With the parameters obtained by fitting the observed SED, the spatial variations of photon index and surface brightness observed in the X-ray band can also be well reproduced. Moreover, it can be derived that the present-day diffusion coefficient of MSH 15–52 at the termination shock is κ{sub 0}=6.6×10{sup 24} cm{sup 2} s{sup −1}, the spatial average has a value of κ-bar =1.4×10{sup 25} cm{sup 2} s{sup −1}, and the present-day magnetic field at the termination shock has a value of B{sub 0}=26.6 μG and the spatial averaged magnetic field is B-bar =14.9 μG. The spatial changes of the spectral index and surface brightness at different bands are predicted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗
DOE OSTI2018

Drift turbulence, particle transport, and anomalous dissipation at the reconnecting magnetopause

Here, using fully kinetic 3D simulations, the reconnection dynamics of asymmetric current sheets are examined at the Earth's magnetopause. The plasma parameters are selected to model MMS magnetopause diffusion region crossings with guide fields of 0.1, 0.4, and 1 of the reconnecting magnetosheath field. In each case, strong drift-wave fluctuations are observed in the lower-hybrid frequency range at the steep density gradient across the magnetospheric separatrix. These fluctuations give rise to cross-field electron particle transport. In addition, this turbulent mixing leads to significantly enhanced electron parallel heating in comparison to 2D simulations. We study three different methods of quantifying the anomalous dissipation produced by the drift fluctuations, based on spatial averaging, temporal averaging, and temporal averaging followed by integrating along magnetic field lines. A comparison of different methods reveals complications in identifying and measuring the anomalous dissipation. Nevertheless, the anomalous dissipation from short wavelength drift fluctuations appears weak for each case, and the reconnection rates observed in 3D are nearly the same as in 2D models. Finally, the 3D simulations feature a number of interesting new features that are consistent with recent MMS observations, including cold beams of magnetosheath electrons that penetrate into the hotter magnetospheric inflow, the related observation of decreasing temperature in regions of increasing total density, and an effective turbulent diffusion coefficient that agrees with predictions from quasi-linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effects of defect density on ultrathin graphene-based metal diffusion barriers

The authors investigated the effects of defect density on the performance of monolayer graphene as a barrier to metal diffusion. The defects were introduced to the graphene by controlled ultraviolet-ozone irradiation. The barrier performance of pristine graphene was found to be superior to that of defective graphene at temperatures up to 700 °C. Changes in surface morphology were more prevalent in the defective graphene-based films than in the pristine graphene-based film; the thermal stability of graphene films depends on their defect density. Defect density was found to be a determining factor in the barrier performance of graphene.

42 ENGINEERING↗
DOE OSTI2023

The role of isotope mass and transport for H-mode access in tritium containing plasmas at JET with ITER-like wall

The required heating power, $P_{\mathrm{LH}}$, to access the high confinement regime (H-mode) in tritium containing plasmas is investigated in JET with ITER-like wall at a toroidal magnetic field of $B_t = 1.8$ T and a plasma current of $I_p = 1.7$ MA. $P_{\mathrm{LH}}$, also referred to as the L-H power threshold, is determined in plasmas of pure tritium as well as mixtures of hydrogen with tritium (H-T) and mixtures of deuterium with tritium (D-T), and is compared to the L-H power threshold in plasmas of pure hydrogen and pure deuterium. It is found that, for otherwise constant parameters, $P_{\mathrm{LH}}$ is not the same in plasmas with the same effective isotope mass, $A_{\mathrm{eff}}$, when they differ in their isotope composition. Thus, $A_{\mathrm{eff}}$ is not sufficient to describe the isotope effect of $P_{\mathrm{LH}}$ in a consistent manner for all considered isotopes and isotope mixtures. The electron temperature profiles measured at the L-H transition in the outer half of the radius are very similar for all isotopes and isotope mixtures, despite the fact that the L-H power threshold varies by a factor of about six. This finding, together with the observation of an offset linear relation between the L-H power threshold, $P_{\mathrm{LH}}$, and an effective heat diffusivity, $\chi_{\mathrm{eff}}$, indicates that the composition-dependent heat transport in the low confinement mode (L-mode) determines, how much power is needed to reach the necessary electron temperatures at the edge, and hence $P_{\mathrm{LH}}$.

tritium plasmas↗
DOE OSTI2025

Validation of a comprehensive first-principles-based framework for predicting the performance of future stellarators

This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.

Gyrokinetic simulations↗
DOE OSTI2024

Topology optimization for the full-cell design of porous electrodes in electrochemical energy storage devices

In this paper, we introduce a density-based topology optimization framework to design porous electrodes for maximum energy storage. We simulate the full cell with a model that incorporates electronic potential, ionic potential, and electrolyte concentration. The system consists of three materials, namely pure liquid electrolyte and the porous solids of the anode and cathode, for which we determine the optimal placement. We use separate electronic potentials to model each electrode, which allows interdigitated designs. As a result, a penalization is required to ensure that the anode and cathode do not touch, i.e., causing a short circuit. We compare multiple 2D designs generated for different fixed conditions, e.g. material properties. A 3D design with complex channel and interlocked structure is also created. All optimized designs are far superior to the traditional monolithic electrode design with respect to energy storage metrics. We observe up to a 750% increase in energy storage for cases with slow effective ionic diffusion within the porous electrode.

25 ENERGY STORAGE↗
DOE OSTI2015

First-principles-based kinetic Monte Carlo studies of diffusion of hydrogen in Ni–Al and Ni–Fe binary alloys

The diffusion of dilute hydrogen in fcc Ni–Al and Ni–Fe binary alloys was examined using kinetic Monte Carlo method with input kinetic parameters obtained from first-principles density functional theory. The simulation involves the implementation of computationally efficient energy barrier model that describes the configuration dependence of the hydrogen hopping. The predicted hydrogen diffusion coefficients in Ni and Ni 89.4 Fe 10.6 are compared well with the available experimental data. In Ni–Al, the model predicts lower hydrogen diffusivity compared to that in Ni. Overall, diffusion prefactors and the effective activation energies of H in Ni–Fe and Ni–Al are concentration dependent of the alloying element. Furthermore, the changes in their values are the results of the short-range order (nearest-neighbor) effect on the interstitial diffusion of hydrogen in fcc Ni-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulations of n -dodecane/oxygen/nitrogen cellular detonations

In this work, two-dimensional n-dodecane/air/nitrogen cellular detonations are simulated with various equivalence ratios (ERs). A skeletal mechanism consisting of 54 species and 269 reactions is used. The lower and upper equivalence ratio boundaries for self-sustained detonation are 0.3 and 2.2, respectively. Detonation with different regimes characterized by the detonation cell patterns is observed, which aligns well with the category based on the stability parameter, i.e., weakly and highly unstable detonations, and extinction. Further, in terms of the frontal structure, non-negligible effect of diffusion on the cellular detonation is revealed, especially in the vicinity of the leading shock front. In highly unstable and quenched detonations, the alternation in reaction pathway within the induction zone accounts for the changes of detonation dynamics, such as the absence or extended sequence of important radical formation, e.g., OH. In addition, the composition of the unburned pockets depends on both pocket location from the leading shock front and the ER in the fresh mixture, because the former determines the residence time, whilst the latter affects the pocket reaction rate.

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

Determination of mineral dissolution regimes using flow-through time-resolved analysis (FT-TRA) and numerical simulation

Flow-through time resolved analysis (FT-TRA) involves subjecting small mineral samples (< 10 mg) inserted in a miniature flow-through cell (50 μL) to controlled flows of eluent analyzed on-line by ICP-MS. In this study, FT-TRA is used to empirically determine the dissolution regimes for the two well-studied minerals forsterite and calcite, representing minerals with relatively slow and fast dissolution kinetics. A proportional increase in steady-state effluent [Mg, Si] concentrations with increasing flow-through cell eluent residence times confirms a dominantly surface-controlled dissolution regime for a powdered forsterite sample at pH 2.3, implying that transport limitations are negligible. In contrast, the relationship between flow rates and dissolution rates for single grain calcite samples at pH 2.3-4 reveals that transport limitations affect the rate of calcite dissolution. In order to provide a quantitative and process-based assessment of the effect of diffusive transport limitations, simulations of the calcite experiments were performed with a high resolution, pore-scale model that considers the geometry of the calcite grain and the FT-TRA flow-through reactor. The pore-scale model reproduces the observed effluent [Ca] concentrations for all experimental conditions using a single set of surface kinetic parameters, by accounting for the formation of a diffusive boundary layer (DBL) that varies in thickness as a function of flow rates. These results demonstrate that combining FT-TRA with pore-scale modeling makes it possible to obtain unprecedented insights not achievable by either method separately, including quantification of DBL thicknesses and the determination of transport controls as a function of pH, flow velocity and residence times.

58 GEOSCIENCES↗