Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse functions”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗
DOE OSTI2021

DECaNT: Simulation tool for diffusion of excitons in carbon nanotube films

In this work, we present the numerical tool DECaNT (Diffusion of Excitons in Carbon NanoTubes) that simulates exciton transport in thin films of carbon nanotubes. Through a mesh of nanotubes generated using the Bullet Physics C++ library, excitons move according to an ensemble Monte Carlo algorithm, with the scattering rates that account for tube chirality, orientation, and distance. We calculate the diffusion tensor from the position–position correlation functions and analyze its anisotropy and dependence on the film composition, morphology, and defect density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2020

Determination of Stress-Free Lattice Spacing (D₀) for Residual Stress Relaxation Measurement in Ni-based Superalloys by Neutron Diffraction

The residual stress and its relaxation in critical engineering components are key to structural materials reliability. The determination of residual stress in bulk engineering structures like Ni-based superalloys can be conducted by measuring the lattice strain with penetrating neutron diffraction. The state-of-the-art engineering diffraction beamline VULCAN at the Oak Ridge National Laboratory (ORNL) Spallation Neutron Source (SNS) provides the critical capability to evaluate the residual strain/stress relaxation in engineering component thanks to the high flux and event data features. By the means of residual stress/strain calculation, it is critical to accurately determine the stress-free lattice spacing d₀, which can be altered by the change of chemistry, unrelieved stress, and measurement scheme. Here we reported d₀ determination for determining residual stress relaxation and distribution in Alloy 718 superalloys after different quenching treatments. Selection of locations for d₀ measurement was discussed by considering the neutron path, attenuation and sample alignment. A dynamic d₀ that resulted from atom diffusion and chemistry change was estimated as function of time for the in-situ relaxation characterization. It demonstrated how the dynamic d₀ values may influence on the strain calculation in different thermally treated samples. The results highlight the importance of dynamic d₀ for in-situ relaxation involving chemical changes and provide guidance to the dynamic d₀ measurement and simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Polyelectrolyte Surface Diffusion in a Nanoslit Geometry

The surface diffusion of poly-l-lysine (PLL) in a planar nanoslit was studied using convex lens-induced confinement (CLiC) single-molecule tracking microscopy. Here, three surface chemistries were employed to understand the interplay of electrostatic and short-range interactions: an amine-functionalized silica surface, an oligo(ethylene oxide) (OEG)-modified surface, and a 1:1 mixture of the two ligands. Effective surface diffusion coefficients increased rapidly with slit height until saturating for slit heights <30 nm. While diffusion at a semi-infinite interface was significantly faster for OEG surfaces, the diffusion coefficient increased most rapidly with slit height for amine-functionalized surfaces, resulting in surface diffusion within very thin slits being nearly independent of surface chemistry. Intermittent random walks were simulated within a planar slit geometry, using experimentally measured parameters obtained from diffusion at a single interface to account for the characteristic short-range interactions between PLL and each surface chemistry, and were in good agreement with experimental measurements.

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

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗
DOE OSTI2025

Helium bubble formation in aged plutonium metal at ambient conditions

A model of helium bubble formation is developed that assumes the initial nascent bubbles are generated when two diffusing helium atoms come within a distance of about the interatomic spacing of the plutonium atoms. These nascent bubbles are assumed to be locked in place and grow as additional helium atoms diffuse across their surfaces. Calculations of the helium-bubble number density and sizes as a function of age are compared to experiment. This leads to an inferred helium diffusion coefficient of-the-order of 10 -20 m 2 /s, supports the idea that most of the helium from alpha-decay end up in bubbles, and that previous transmission electron microscopy size estimates are biased low. The energy that helium can transfer to 40-year-old plutonium, upon expansion, is determined to be <0.1 Joules per gram of metal. The model predicts differences between the bubble properties in naturally and accelerated aged plutonium.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗
DOE OSTI2020

Validation of Thermal Scattering Laws for Light Water at Elevated Temperatures with Diffusion Experiments

Light water is the most important neutron moderator in many reactor applications. Thermal neutron scattering kernels for H bound in H 2 O are represented by ENDF-format thermal scattering laws (TSLs) evaluated at discrete temperatures T . Nuclear data evaluations are conventionally validated utilizing a combination of critical benchmarks and experimental cross section data. Existing public critical benchmarks may be inadequate for testing water TSLs over the full range of T of interest in reactor applications, and experimental thermal scattering cross section data for water is sparse at elevated T . In this work, MC21 is used to simulate the decay of the equilibrium thermal neutron flux in light-water spheres of several radii. The fundamental-mode time eigenvalue $\mathcal{a}$ is calculated at 22 °C and 227 °C (at saturation pressure) for each sphere using three different H-H 2 O TSLs. Fitting $\mathcal{a}$ to a polynomial function of geometric buckling allows calculation of the thermal neutron diffusion length L . Extrapolation lengths are treated as a function of the transport mean free path and geometry. Experimental L data from 25 publications at 49 temperatures (from 10 °C to 295 °C), computed by several different time-dependent and space-dependent decay methods, is used to develop an empirical fit of L vs. T . The MC21-calculated L for the TSLs tested are within ±1% of the predicted values at 22 °C and 227 °C. This validation approach, which may be repeated at arbitrary T , constitutes an integral benchmark specific to and characterizing the detailed physics of the thermal scattering kernel applied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study

Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Solution-Phase Growth of Cu Nanowires with Aspect Ratios Greater Than 1000: Multiscale Theory

Penta-twinned metal nanowires are finding widespread application in existing and emerging technologies. However, little is known about their growth mechanisms. Here, we probe the origins of chloride- and alkylamine-mediated, solution-phase growth of penta-twinned Cu nanowires from first principles using multi-scale theory. Using quantum density functional theory (DFT) calculations, we characterize the binding and surface diffusion of Cu atoms on chlorine-covered Cu(100) and Cu(111) surfaces. We find stronger binding and slower diffusion of Cu atoms on chlorinated Cu(111) than on chlorinated Cu(100) – a reversal of the trend for bare Cu surfaces. We also probe inter-facet diffusion and find this proceeds faster from Cu(100) to Cu(111) than the reverse. Using the DFT rates for hopping between individual sites at Ångstrom scales, we calculate coarse-grained, inter-facet rates for nanowires of various lengths – up to hundreds of μm – and diameters in the 10-nm range. We predict nanowires with aspect ratios around 100 based on surface diffusion, alone. We also account for the influence of a self-assembled alkylamine layer that covers most of the {100} facets, but is absent or thin and disordered on the {111} facets and in an “end zone” near the {100}-{111} boundary. With an end zone, we predict a wide range of nanowire aspect ratios in the experimental ranges. Our work reveals the mechanisms by which a halide – chloride – promotes the growth of high-aspect-ratio nanowires.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First-principles studies of the concentration-dependent tritium diffusion in the zirconium hydrides with and without Sn impurity

We performed first-principles calculations to study the 3 H diffusion in the zirconium hydrides with or without Sn impurity. Our results show that the formation of Zr 3 H x becomes preferable as the 3 H concentration increases. Based on the most stable configurations of Zr 3 H x (x = 0.5, 1.0, 1.5, 2.0), we studied the interstitial 3 H diffusion with or without Sn impurity. The results show that the interstitial 3 H diffusion becomes less probable in the hydrides with a higher 3 H concentration. When introducing a substitutional Sn impurity, the diffusion barriers increase. Therefore, the substitutional Sn in the diffusion pathway hampers the 3 H diffusion in Zr hydrides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Graphite Lithiation under Fast Charging Conditions: Atomistic Modeling Insights

The charging of lithium ion batteries in a fast and safe manner is critical for promoting the mass adoption of electric vehicles. Li intercalation in graphite electrodes is known to be one of the bottlenecks during the fast charging process. The mechanism of Li diffusion in highly polarized graphite anode at high current rates remains, however, not well understood. In this work, Density Functional Theory (DFT) calculations are used to gain insights into the Li diffusion process in graphite when it is far from equilibrium under fast charging conditions. The effect of uncompensated charges on Li mobility is determined in the highly polarized regions of the anode close to the interfaces. The extra charge was found to increase the interlayer spacing in the diffusion layer and adjacent channels, increasing the diffusivity and promoting the formation of Li clusters. A concerted diffusion mechanism at the edge of high-concentration Li domains is proposed to enhance the diffusion of Li.

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

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2022

Report summarizing demonstration of thermal conductivity measurement methods for coated zirconium cladding

Thermal diffusivity measurements on Zirconium-based cladding materials have historically been a challenge due to the difficulty to measure on specimens with curved geometries, including nuclear grade Zircaloy cladding materials. In this work, we first used Laser Flash Analysis methods to measure the thermal diffusivity of Zircaloy tubes, which show good agreement with Zircaloy plates in this work, as well as previously published data. The consistent results proved the applicability of the Laser Flash Analysis setup for investigating thermal diffusivity of Zircaloy tubes. We further investigate the hydrogen effect on thermal diffusivity of Zircaloy tubes. Hydrogen plays significant roles in temperature dependent thermal diffusivity trends of Zircaloy tube. For low hydrogen concentration, where hydride dissolved into single phase α-Zr at higher temperature, thermal diffusivity is relatively unchanged at low temperature at two phase regions (α-Zr and δ-hydride), and it increases as a function of temperature at single phase α-Zr region. For high hydrogen concentration, where phase transformation (α-Zr + δ-hydride → α-Zr + β-Zr) occurs at 567 °C, thermal diffusivity decreases as a function of temperature below 567 °C before the transformation (α-Zr and δ-hydride), while increase at higher temperature after the transformation (α-Zr and β-Zr).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Boric acid removal with polyol-functionalized polyether membranes

Poor selectivity of conventional desalination membranes for boron often necessitates multi-stage treatment trains to achieve desired boron removal for end uses such as irrigation. One approach to membrane design for improved single-pass boric acid removal is via incorporation of chelating ligands that selectively sorb boron. In this study, membranes based on epoxy-amine cross-linked poly(glycidyl glycidyl ether) (PGGE) were synthesized and functionalized with N-methyl-d-glucamine (NMDG), a polyol known to interact selectively with boron. PGGE and PGGE-NMDG membranes exhibited boron sorption isotherms that were well-described by dual mode isotherms. PGGE-NMDG sorbed 2.5 mmol B/g dry polymer from a neutral aqueous solution containing 100 mmol B/L, which was almost three times the adsorption density of a commercial boron selective resin, Amberlite IRA743. Further, the membranes were regenerated in acid without a significant loss of boron sorption capacity over four cycles. Interactions between boron and NMDG and PGGE sites (e.g., epoxides) impacted boron diffusion in both membranes. The use of ligand functionalized membranes to capture target solutes such as boron requires high loading of interacting sites to maximize uptake capacity. Establishing fundamental structure/property rules for boron selectivity could lead to new material designs with improved boron separation properties for water purification.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Selective acetylene hydrogenation over single metal atoms supported on Fe 3 O 4 (001): A first-principle study

Supported single-atom catalysts (SACs) have gained increasing attention for improved catalytic activity and selectivity for industrially-relevant reactions. In this study, we explore the hydrogenation of acetylene over single Pt, Ru, Rh, Pd, and Ir atoms supported on the Fe 3 O 4 (001) surface using density functional theory (DFT) calculations. The thermodynamic profile of H diffusion is significantly changed with the type of single metal atoms used and suggests that H spillover from the single atom dopant to the Fe 3 O 4 (001) surface is favored and will likely lead to high H coverages in the functioning catalyst. Correspondingly, as the surface H coverage increases, the important desorption step of ethylene becomes energetically competitive against the detrimental hydrogenation steps of ethylene to ethane. Overall, we show herein that the selectivity of SACs towards ethylene production can be tuned by considering the proper type of metal and controlling the redox state of the support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗