Results for “diffusion behavior” — Page 3

DOE OSTI2024

Unraveling Interdiffusion Phenomena and the Role of Nanoscale Diffusion Barriers in the Copper–Gold System

Diffusion is one of the most fundamental concepts in materials science, playing a pivotal role in materials synthesis, forming, and degradation. Of particular importance is solid state interdiffusion of metals which defines the usable parameter space for material combinations in the form of alloys. This parameter space can be explored on the macroscopic scale by using diffusion couples. However, this method reaches its limit when going to low temperatures, small scales, and when testing ultrathin diffusion barriers. Therefore, this work transfers the principle of the diffusion couples to small scales by using core–shell nanowires and in situ heating. This allows us to delve into the interdiffusion dynamics of copper and gold, revealing the interplay between diffusion and the disorder–order phase transition. Our in situ TEM experiments in combination with chemical mapping reveal the interdiffusion coefficients of Cu and Au at low temperatures and highlight the impact of ordering processes on the diffusion behavior. The formation of ordered domains within the solid-solution is examined using high-resolution imaging and nanodiffraction including strain mapping. In addition, we examine the effectiveness of ultrathin Al 2 O 3 barrier layers to control interdiffusion of the diffusion couple. Our findings indicate that a 5 nm thick layer serves as an efficient diffusion barrier. Furthermore, this research provides valuable insights into the interdiffusion behavior of Cu and Au on the nanoscale, offering potential applications in the development of miniaturized integrated circuits and nanodevices.

alloys↗

DOE OSTI2023

Diffusion of small anti-Schottky clusters in UO 2

The aggregation of irradiation-induced defect clusters in UO 2 leads to the formation of dislocation loops and cavities which contribute swelling and ultimately produce deleterious effects on the fuel. Despite their fundamental role in fuel evolution, the kinetics of these defect clusters are not currently well understood. In this work, we investigate the diffusive behavior of interstitial clusters via Molecular Dynamics (MD) simulations. In this work, our investigation considers a range of defect cluster sizes, each composed of N UO 2 units or anti-Schottky defects. We report a complex cluster size - mobility relation; increasing cluster size corresponds to an increase in mobility up until a critical size (N=4) where diffusivity reaches a maximum, after which the trend reverts and further increasing size corresponds to a decrease in mobility. The rapid migration observed correlates well with the very low barriers reported in a few historic experimental studies. Further analysis of the cluster shape and orientation reveals a strong structural preference for near-planar configurations oriented normal to the $\langle$100$\rangle$ direction. Rotational energy barriers are found to be similar in both magnitude and trend to the observed diffusion barriers. Lastly, connections are drawn to the behavior of large defect clusters including the formation of (1/3) $\langle$111$\rangle$ Frank dislocation loops.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Reaction-diffusive dynamics of number-conserving dissipative quantum state preparation

The use of dissipation for the controlled creation of nontrivial quantum many-body correlated states is of much fundamental and practical interest. What is the result of imposing number conservation, which, in closed system, gives rise to diffusive spreading? For this work, we investigate this question for a paradigmatic model of a two-band system, with dissipative dynamics aiming to empty one band and to populate the other, which had been introduced before for the dissipative stabilization of topological states. Going beyond the mean-field treatment of the dissipative dynamics, we demonstrate the emergence of a diffusive regime for the particle and hole density modes at intermediate length- and timescales, which, interestingly, can only be excited in nonlinear response to external fields. We also identify processes that limit the diffusive behavior of this mode at the longest length- and timescales. Strikingly, we find that these processes lead to a reaction-diffusion dynamics governed by the Fisher-Kolmogorov-Petrovsky-Piskunov equation, making the designed dark state unstable towards a state with a finite particle and hole density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Investigation of Molecular Diffusion at Block Copolymer Thin Films Using Maximum Entropy Method-Based Fluorescence Correlation Spectroscopy and Single Molecule Tracking

Fluorescence correlation spectroscopy (FCS) has been widely used to investigate molecular diffusion behavior in various samples. The use of the maximum entropy method (MEM) for FCS data analysis provides a unique means to determine multiple distinct diffusion coefficients without a priori assumption of their number. Comparison of the MEM-based FCS method (MEM-FCS) with another method will reveal its utility and advantage as an analytical tool to investigate diffusion dynamics. Herein, we measured diffusion of fluorescent probes doped into nanostructured thin films using MEM-FCS, and validated the results with single molecule tracking (SMT) data. The efficacy of the MEM code employed was first demonstrated by analyzing simulated FCS data for systems incorporating one and two diffusion modes with broadly distributed diffusion coefficients. The MEM analysis accurately afforded the number of distinct diffusion modes and their mean diffusion coefficients. These results contrasted with those obtained by fitting the simulated data to conventional two-component and anomalous diffusion models, which yielded inaccurate estimates of the diffusion coefficients. Subsequently, the MEM analysis was applied to FCS data acquired from hydrophilic dye molecules incorporated into microphase-separated polystyrene-block-poly(ethylene oxide) (PS-b-PEO) thin films characterized under a water-saturated nitrogen atmosphere. The MEM analysis revealed distinct fast and slow diffusion components attributable to molecules diffusing on the film surface and inside the film, respectively. SMT studies of the same materials yielded trajectories for mobile molecules that appear to follow the curved PEO microdomains. Diffusion coefficients obtained from the SMT data were consistent with those obtained for the slow diffusion component detected by MEM-FCS. Furthermore, these results highlight the utility of MEM-FCS and SMT for gaining complementary information on molecular diffusion processes in heterogeneous material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effective field theory of fluctuating wall in open systems: from a kink in Josephson junction to general domain wall

We investigate macroscopic behaviors of fluctuating domain walls in nonequilibrium open systems with the help of the effective field theory based on symmetry. Since the domain wall in open systems breaks the translational symmetry, there appears a gapless excitation identified as the Nambu-Goldstone (NG) mode, which shows the non-propagating diffusive behavior in contrast to those in closed systems. After demonstrating the presence of the diffusive NG mode in the (2+1)-dimensional dissipative Josephson junction, we provide a symmetry-based general analysis for open systems breaking the one-dimensional translational symmetry. A general effective Lagrangian is constructed based on the Schwinger-Keldysh formalism, which supports the presence of the gapless diffusion mode in the fluctuation spectrum in the thin wall regime. Besides, we also identify a term peculiar to the open system, which possibly leads to the instability in the thick-wall regime or the nonlinear Kardar-Parisi-Zhang coupling in the thin-wall regime although it is absent in the Josephson junction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗

DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗

DOE OSTI2020

Deuterium diffusion in γ-LiAlO2 pellets irradiated with He+ and D2+ ions

Not Available

Jiang, Weilin↗

DOE OSTI2024

Underlying mechanism of structural transformation between GaSb and GaAs response to intense electronic excitation

Ion irradiation of semiconductors has emerged as a promising approach for fabricating self-organized nanosystems with high atomic precision, despite often being accompanied by undesirable phenomena. Exploring the mechanisms underlying structural transformations is crucial for assessing nanostructure array types under complex irradiation environments. By quantitatively calculating the thermodynamically driven processes and analyzing the impact of intrinsic structural parameters, distinct structural transformations in response to intense electronic excitation are systematically investigated in gallium antimonide (GaSb) and gallium arsenide (GaAs) systems. In high-energy regimes, the nanofibers layer of GaSb exhibits intriguing structural discrepancy, characterized by partial nanofibers with coherent boundaries, interspersed nanopores accompanied by antisite defects and Ga precipitates, distinguishing to a series of discontinuous latent tracks that emerged within cylindrical trajectories in GaAs. Furthermore, significant diffusion behaviors of the nanohillocks are discovered in GaAs, with higher average roughness than GaSb, driven by the gradient stress distribution influenced by the free-surface effects. The deposition energy for melting phase formation, Gibbs free energy, and Ga diffusion coefficients contribute to the distinctive structural features, evidencing relatively stable morphological configurations and higher irradiation resistance in GaAs. Consequently, special optoelectronic properties associated with structural discrepancies facilitate the design and optimization of material functionalities by irradiation technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effect of lithium incorporation on tweaking the electrocatalytic behavior of tantalum-based oxides

To study the effect of lithium incorporation on the electrochemical properties and catalytic behavior of tantalum-based oxides, a series of mesoporous lithium-contained tantalum oxides was developed adopting an inverse micelle-based sol gel synthesis method. Four various phases of amorphous Ta 2 O 5 , polycrystalline LiTaO 3 , Li 5.04 Ta 6.16 O 17.92 , and Li 3 TaO 4 were determined for samples synthesized at Li/Ta ratios of 0.5 (LTO-0.5), 1 (LTO-1), 2 (LTO-2), and 3 (LTO-3) respectively. The surface analysis revealed that LTO-2 sample is featured with lower amount of lithium atoms at the surface rather than the bulk of the sample. The electrochemical analysis of LTO samples showed that LTO-2 possesses lower charge transfer resistance extracted from impedance spectroscopy that consequently leads to a better electrocatalytic performance. The LTO-2 electrode also showed a low open circuit potential (OCP) combined with a relatively higher electroactive surface area, thus capable of running an electrooxidation reaction at a potential lower than all other LTO samples. The more negative OCP value of LTO-2 electrode in AA solution states a more facile electron transfer between the analyte and electrode compared to other synthesized samples. Mott–Schottky plots of LTO electrodes are characterized by positive slopes which indicate n-type semiconductor behavior. The charge carrier density (5.5 × 10 19 cm -3 ) in LTO-2 sample is the lowest among all other samples. The charge carrier density is dependent on the amount of applied lithium in the synthesis confirming various lithium diffusion behaviors in LTO samples. The lower carrier concentration in LTO-2 justifies its higher activity in electrooxidation providing less populated lower unoccupied molecular orbital (LUMO) for the accommodation of electrons extracted from the analyte higher occupied molecular orbital (HOMO). Therefore, LTO-2 showed an improved electrocatalytic performance due to lower charge transfer resistance (higher conductivity) offered by lower amount of lithium atom at its surface and higher oxygen vacancies in the bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

NASA NTRS2000

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗

DOE OSTI2023

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗