Search NASA⌕ Search

SEARCH · Search NASA

Results for “DIFFUSION EFFECT”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiscale, mechanistic modeling of irradiation-enhanced silver diffusion in TRISO particles

Tristructural isotropic (TRISO) particles are under consideration for use in several proposed advanced nuclear reactor concepts. The silicon carbide (SiC) layer in TRISO acts as a barrier to prevent the release of the fission products. However, despite remarkable retention, silver (Ag) release has been observed from intact particles, which requires investigation since the Ag isotope ( 110m Ag) has a long half-life. Previous work focused on developing a multiscale, mechanistic model for Ag diffusion accounting for temperature and microstructure effect and has been successfully validated. In this work, we expand the previous model to account for irradiation-enhanced Ag diffusivity in SiC and improve its accuracy over a wider grain size and temperature ranges relevant for advanced reactor conditions. A temperature, grain size, and flux dependent diffusivity is therefore derived using the mesoscale code MARMOT and implemented in the fuel performance code BISON. The irradiation-enhanced Ag diffusivity in SiC is compared against experimental data and validated using BISON against Ag release measurements from the Advanced Gas Reactor Fuel Development and Qualification Program (AGR-1 and AGR-2). Herein, we quantify the impact of SiC grain size, irradiation, and temperature on Ag release. In agreement with previous studies, we find accounting for SiC grain size improves agreement between BISON predictions and experimental observations for most cases. In conclusion, we also find that accounting for irradiation improves agreement for cases where Ag release was underestimated, but the impact was less significant than accounting for microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Accelerated oxidation of epoxy thermosets with increased O 2 pressure

Polymer oxidation is usually accelerated with temperature, which is therefore applied in nearly every experimental approach dealing with predictive materials aging. Because of this simple approach, we may tend to neglect that the effective concentration of oxygen also acts as a rate multiplier for oxidation. Increasing the oxygen partial pressure in an aging environment accelerates oxidation, and while there is often a near proportional increase initially, the effect of additional oxygen usually transitions to a saturation level at some elevated pressure. This has been theoretically described in the general autoxidation scheme and is well recognized. However, for many materials the exact rate behavior under moderately increased oxygen concentration remains to be established. We therefore review epoxy oxidation and offer a broader overview of its behavior under increased O 2 partial pressure. Experimental data are given for a few thermoset materials demonstrating their rate behavior under O 2 partial pressure up to 4 atm, meaning approximately 20 times more than under standard atmospheric conditions. Confirmative evidence suggests that epoxy materials will reach saturation oxidation rates only at significantly higher O 2 partial pressure. In such a high-pressure regime it is theoretically possible to not only accelerate oxidation, but to transition into a condition where O 2 diffusion can be increased without further accelerating the oxidation rate. Finally, this can reduce diffusion limited oxidation effects under specific accelerated aging conditions as a combination of temperature and O 2 partial pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Water Dynamics of Superacid Aromatic Proton Exchange Membranes for Fuel Cell Applications

Proton exchange membranes (PEMs) with high conductivity are of critical importance for the development of fuel cells, electrolyzers, and other electrochemical technologies. In this research, poly(1,1,2,2-tetrafluoro-2-phenoxyethane-1-sulfonic acid) (PTPS) with an aromatic polymer main chain and a perfluorinated superacidic polymer side chain was synthesized. The water dynamics of PTPS were characterized across various length scales using a combination of Fourier-transform infrared spectroscopy (FTIR) and nuclear magnetic resonance (NMR) and compared with Nafion, a standard perfluorinated PEM, and sulfonated poly(ether sulfone) (SPES 40), an aromatic PEM without perfluorinated superacid side chains. The T 1 and T 2 relaxation times of water in the samples probed by NMR increase from SPES 40 to PTPS to Nafion, indicating that the local motion of the water molecules becomes faster. This trend corresponds well with the relative fraction of bulk-like water determined using FTIR. At larger length scales, the diffusion coefficient of water was characterized using pulsed-field gradient NMR (PFG-NMR). At a longer diffusion time (Δ = 100 ms), PTPS has a smaller diffusion coefficient compared with both Nafion and SPES 40, due to restricted diffusion, and this effect is also evident in the proton conductivity of the hydrated membranes. From this comparison, it is apparent that the aromatic backbone and side chain type greatly influence the water dynamics in PEMs at various length scales and the water dynamics significantly impact the bulk proton conductivity. These insights will lead to new designs for aromatic PEMs and help to identify bottlenecks in current materials.

25 ENERGY STORAGE↗
DOE OSTI2026

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2026

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Tackling the curse of dimensionality in fractional and tempered fractional PDEs with physics-informed neural networks

Fractional and tempered fractional partial differential equations (PDEs) are effective models of long-range interactions, anomalous diffusion, and non-local effects. Traditional numerical methods for these problems are mesh-based, thus struggling with the curse of dimensionality (CoD). Physics-informed neural networks (PINNs) offer a promising solution due to their universal approximation, generalization ability, and mesh-free training. In principle, Monte Carlo fractional PINN (MC-fPINN) estimates fractional derivatives using Monte Carlo methods and thus could lift CoD. However, this may cause significant variance and errors, hence affecting convergence; in addition, MC-fPINN is sensitive to hyperparameters. In general, numerical methods and specifically PINNs for tempered fractional PDEs are under-developed. Herein, we extend MC-fPINN to tempered fractional PDEs to address these issues, resulting in the Monte Carlo tempered fractional PINN (MC-tfPINN). To reduce possible high variance and errors from Monte Carlo sampling, we replace the one-dimensional (1D) Monte Carlo with 1D Gaussian quadrature, applicable to both MC-fPINN and MC-tfPINN. We validate our methods on various forward and inverse problems of fractional and tempered fractional PDEs, scaling up to 100,000 dimensions. Our improved MC-fPINN/MC-tfPINN using quadrature consistently outperforms the original versions in accuracy and convergence speed in very high dimensions.

42 ENGINEERING↗
DOE OSTI2025

Understanding grain refinement and intergranular gas bubble evolution in U-10Mo fuel using phase-field modeling

Monolithic U-10Mo fuel undergoes significant microstructural changes in the form of grain refinement and gas bubble formation during burnup, which degrades its mechanical properties. In this talk, I present a phase-field model for microstructure evolution in U-10Mo developed using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Simulations demonstrate that grain refinement initiates at pre-existing grain boundaries (GBs) to eliminate the lattice distortion energy caused by the accumulation of self-interstitial loops. By employing an equation of state for xenon gas, we simulate the evolution of gas bubbles in the polycrystal microstructure. Large, interconnected bubbles are found to form along the triple junctions. The effects of defect production rate, diffusivities and GB mobility on the microstructure evolution are systematically studied. Homogenization is employed on the microstructures to obtain effective elastic constants and diffusivity as a function of fission density. The simulations provide critical insights on microstructure and property degradation in U-10Mo fuel.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Suppressing proximity effects during rapid serial two-photon lithography through tuning of reaction–diffusion kinetics

The ability of two-photon lithography (TPL) to print deterministic nanoporous 3D structures is highly valuable for many applications. However, it is challenging to print such structures rapidly due to the proximity effects that cause closely spaced features to enlarge and merge. A key challenge is the limited understanding of the origins and spatiotemporal dynamics of the long-range proximity effects that extend beyond the optical focal spot. Here, we empirically investigate these proximity effects in serial TPL using custom-made acrylate-based photoresists. We demonstrate that the complex spatiotemporal dynamics of the proximity effects can be explained through the kinetics of the reaction–diffusion photopolymerization mechanisms that underlie the curing process. Specifically, we show that the proximity effects arise due to comparable timescales of reaction and diffusion of oxygen. Furthermore, we demonstrate that long-range proximity effects can be suppressed by introducing phenolic inhibitors, which reduce oxygen consumption through non-inhibiting side reactions and promote faster termination of curing reactions. These insights enabled improving the linewidths from >400 to 260 nm during printing of nanoporous 3D woodpiles at scanning speeds of 50 mm s −1 . Thus, the knowledge generated here can be applied to deterministically tune the proximity effects and enable rapid printing of high-resolution nanoporous 3D structures.

36 MATERIALS SCIENCE↗
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 OSTI2025

Extending the capillary wave model to include the effect of bending rigidity: X-ray reflectivity and diffuse scattering

The surface roughness of a thin film at a liquid interface exhibits contributions of thermally excited fluctuations. This thermal roughness depends on temperature (𝑇), surface tension (𝛾), and elastic material properties, specifically the bending modulus (𝜅) of the film. A nonzero 𝜅 suppresses the thermal roughness at small length scales compared to an interface with zero 𝜅, as expressed by the power spectral density of the thermal roughness. The description of the x-ray scattering of the standard capillary wave model (CWM), which is valid for zero 𝜅, is extended to include the effect of 𝜅. The extended CWM (eCWM) provides a single analytical form for both the specular x-ray reflectivity (XRR) and the diffuse scattering around the specular reflection, and recovers the expression of the CWM at its zero 𝜅 limit. This theoretical approach enables the use of single-shot grazing incidence x-ray off-specular scattering (GIXOS) measurements for characterizing the structure of thin films on a liquid surface. The eCWM analysis approach decouples the thermal roughness factor from the surface scattering signal, providing direct access to the intrinsic surface-normal structure of the film and its bending modulus. Moreover, the eCWM facilitates the calculation of reflectivity at any desired resolution (pseudo-XRR approach). The transformation into pseudo-XRR provides the benefit of using widely available XRR software to perform GIXOS analysis. The extended range of the vertical scattering vector (𝑄 𝑧 ) available with the GIXOS pseudo-XRR approach allows for a higher spatial resolution than with conventional XRR. Experimental results are presented for various lipid systems, showing strong agreement between conventional specular XRR and pseudo-XRR methods. This agreement validates the proposed approach and highlights its utility for analyzing soft, thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
DOE OSTI2025

Water reactive polyurethane grouting for deep vadose zone contaminant immobilization

Contaminants in the deep vadose zone (DVZ) pose a long-term threat to groundwater, human health, and the environment. Polyurethane grouting is a type of chemical grouting that is potentially advantageous in immobilizing contaminants in DVZs. Polyurethane resin has low viscosity that makes it feasible to penetrate and fill the pore space between fine particles. Here, in this study a series of laboratory pressure grouting and leaching tests were conducted to simulate and assess the effectiveness of polyurethane grouting for immobilizing contaminants in DVZs. Soil with fine particles was prepared with 127 I (as iodide) that served as a non-radioactive surrogate for radioactive 129 I. After grouting and curing, leaching tests were used to measure and compare contaminant diffusivity and leachability index values. Additionally, changes in porosity and saturated hydraulic conductivity of grouted soil were measured. X-Ray Computer Tomography (XCT) results showed that the cured polyurethane was distributed nearly homogeneously and approximately half of the voids were filled with cured polyurethane. Grouting reduced the saturated hydraulic conductivity of soil by 37 %. The effective diffusivity decreased by more than 80 % as compared with the ungrouted soil. The leachability index of the grouted soil was 6.5; meeting the criteria established by the U.S. Nuclear Regulatory Commission (NRC) standard. The results obtained in this study provide a valuable assessment of polyurethane grouting for iodide immobilization in the DVZ and indicate this approach may be a viable method for contaminant remediation in DVZ soils.

Contaminant Diffusion↗
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
DOE OSTI2024

Numerical validation of scaling laws for stratified turbulence

Recent theoretical progress using multiscale asymptotic analysis has revealed various possible regimes of stratified turbulence. Notably, buoyancy transport can either be dominated by advection or diffusion, depending on the effective Péclet number of the flow. Two types of asymptotic models have been proposed, which yield measurably different predictions for the characteristic vertical velocity and length scale of the turbulent eddies in both diffusive and non-diffusive regimes. The first, termed a ‘single-scale model’, is designed to describe flow structures having large horizontal and small vertical scales, while the second, termed a ‘multiscale model’, additionally incorporates flow features with small horizontal scales, and reduces to the single-scale model in their absence. By comparing predicted vertical velocity scaling laws with direct numerical simulation data, we show that the multiscale model correctly captures the properties of strongly stratified turbulence within regions dominated by small-scale isotropic motions, whose volume fraction decreases as the stratification increases. Meanwhile its single-scale reduction accurately describes the more orderly, layer-like, quiescent flow outside those regions.

Mechanics↗