Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion processes”

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

Potential–Dependent BDAC Adsorption on Zinc Enabling Selective Suppression of Zinc Corrosion for Energy Storage Applications

Utility-scale zinc (Zn) batteries are a promising solution to address the problem of intermittency of renewable energy sources; however, Zn-metal anodes in these batteries suffer from capacity loss due to spontaneous corrosion of the Zn especially when high-surface area anode configurations are employed. Additionally, Zn dendrites are known to form during battery charging limiting the cycle-life of these batteries. Electrolyte additives have been explored that prevent aforementioned issues, but these too come at a cost, i.e., surface-blocking additives polarize the electrode surface leading to loss in the voltaic and energy efficiencies of the battery. In this contribution, a novel electrolyte additive, benzyldimethylhexadecylammonium chloride (BDAC), is investigated for its ability to suppresses corrosion of Zn in an acidic (pH = 3) electrolyte. An attribute of BDAC distinct from previously studied additives is that it selectively suppresses electrochemical activity when the Zn electrode is at its corrosion potential; however, during high-rate Zn deposition (charging) or stripping (discharging), BDAC is essentially deactivated and thus it does not appreciably polarize the electrode surface, thus minimizing voltaic efficiency losses. This selective corrosion suppression behavior is explored using slow-scan voltammetry, which reveals hysteresis implying a potential- or current-dependent BDAC adsorption mechanism in which BDAC reaches higher surface coverages when the partial currents at the Zn surface are low (e.g., at or near the corrosion potential), but BDAC coverage is reduced considerably when the Zn deposition or stripping rates are increased. Numerical simulations of the BDAC diffusion-adsorption process corroborate this mechanism. Ramifications of our approach to the selective suppression of Zn dendrites are discussed.

25 ENERGY STORAGE↗
DOE OSTI2018

Compact Heat Exchanger Design and Testing for Advanced Reactors and Advanced Power Cycles

The goal of the proposed research is to demonstrate the thermal hydraulic performance of innovative surface geometries in compact heat exchangers used as intermediate heat exchangers (IHXs) and recuperators for the supercritical carbon dioxide (s-CO 2 ) Brayton cycle. Printed-circuit heat exchangers (PCHEs) are the primary compact heat exchangers of interest. The overall objectives are: To develop optimized PCHE designs for different working fluid combinations including helium to s-CO 2 , liquid salt to s-CO 2 , sodium to s-CO 2 , and liquid salt to helium; To experimentally and numerically investigate thermal performance, thermal stress and failure mechanism of PCHEs under various transients; and To study diffusion bonding techniques for elevated-temperature alloys and examine post-test material integrity of the PCHEs. The project objectives were accomplished by defining and executing five different tasks corresponding to these specific objectives. The first task involved a thorough literature review and a selection of IHX candidates with different surface geometries as well as a summary of prototypic operational conditions. The second task involved optimization of PCHE design with numerical analyses of thermal-hydraulic performances and mechanical integrity. The subsequent task dealt with the development of testing facilities and engineering design of PCHE to be tested in s-CO 2 fluid conditions. The next task involved experimental investigation and validation of the thermal-hydraulic performances and thermal stress distribution of prototype PCHEs manufactured with particular surface geometries. The last task involved an investigation of diffusion bonding process and posttest destructive testing to validate mechanical design methods adopted in the design process. The experimental work utilized the two test facilities at The Ohio State University (OSU) including one existing High-Temperature Helium Test Facility (HTHF) and the newly developed s-CO 2 test loop (STL) facility and s-CO 2 test facility at University of Wisconsin – Madison (UW).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Evaluation of Oxygen Consumption as a Sensitive Measure of Electrical Cable Polymer Insulation Degradation

Polymers used as dielectric materials in nuclear grade electrical cables, such as cross-linked polyethylene (XLPE) and ethylene-propylene rubber (EPR), are long-lived and durable. Estimating the decades-long useful service life of polymers experimentally in a feasible timeframe of weeks or months requires greatly accelerated aging in the laboratory at elevated temperatures and/or gamma dose rates. Many limitations of accelerated aging, such as due to artefacts of the process including diffusion limited oxidation and dose rate effects, have been identified that can make accurate lifetime prediction at service conditions challenging. This can be due to differences in degradation mechanisms between rapid and extended aging or to the presence of induction periods during which early aging is undetectable using common metrics of polymer insulation characterization such as tensile elongation at break. Conceptually, an ideal measure of aging would track clearly and predictably with material changes throughout the lifecycle of the evolving material state from early aging through end of useful life. With such a metric accelerated aging could be performed at much milder conditions, conditions closer to those found in long term service, and end of life calculated from the early data. Oxygen consumption has previously been proposed as a sensitive measure of insulation polymer oxidation that might meet such desirable characteristics. In this work, nuclear cable insulation samples of the two most common types, XLPE and EPR, from two of the most sourced manufacturers, RSCC and The Okonite Company, were subjected to gamma irradiation at dose rates of 100, 200, and 300 Gy/hr for up to 42 days at 26°C. The aged samples were used to evaluate oxygen consumption as a sensitive measure of aging for these materials and these conditions. For comparison and validation, the aged samples were also characterized using the more conventional tests of mass change, density, carbonyl index, and yellowness index. For both materials, oxygen consumption during irradiation was found to track consistently with applied dose—oxygen consumption rate was linear with time. Mass change and carbonyl index were also clearly proportional to exposure dose. Yellowness index, a surface sensitive measure, was observed to increase linearly with dose at low exposure levels before plateauing. Sample density was not observed to change under the conditions explored. Oxygen consumption was found to be a particularly sensitive metric to measure the degradation of polymeric cable insulation exposed to radiation. Oxygen consumption was compared to Arrhenius and time-temperature-superposition methods for activation energy calculation using literature date and found to be a useful and viable alternative for lifetime prediction. The techniques developed in this work are anticipated to be valuable in the estimation of insulation aging rates under mild conditions approaching service conditions. Direction comparison of useful life prediction by oxygen consumption with prediction via Arrhenius methodology for thermal aging and existing methods for radiation and combined thermal/radiation aging would serve to further validate this promising approach and refine our understanding of material-specific aging rates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Surface-Enhanced Raman Spectroscopy Combined with Multivariate Analysis for Fingerprinting Clinically Similar Fibromyalgia and Long COVID Syndromes

Fibromyalgia (FM) is a chronic central sensitivity syndrome characterized by augmented pain processing at diffuse body sites and presents as a multimorbid clinical condition. Long COVID (LC) is a heterogenous clinical syndrome that affects 10–20% of individuals following COVID-19 infection. FM and LC share similarities with regard to the pain and other clinical symptoms experienced, thereby posing a challenge for accurate diagnosis. This research explores the feasibility of using surface-enhanced Raman spectroscopy (SERS) combined with soft independent modelling of class analogies (SIMCAs) to develop classification models differentiating LC and FM. Venous blood samples were collected using two supports, dried bloodspot cards (DBS, n = 48 FM and n = 46 LC) and volumetric absorptive micro-sampling tips (VAMS, n = 39 FM and n = 39 LC). A semi-permeable membrane (10 kDa) was used to extract low molecular fraction (LMF) from the blood samples, and Raman spectra were acquired using SERS with gold nanoparticles (AuNPs). Soft independent modelling of class analogy (SIMCA) models developed with spectral data of blood samples collected in VAMS tips showed superior performance with a validation performance of 100% accuracy, sensitivity, and specificity, achieving an excellent classification accuracy of 0.86 area under the curve (AUC). Amide groups, aromatic and acidic amino acids were responsible for the discrimination patterns among FM and LC syndromes, emphasizing the findings from our previous studies. Overall, our results demonstrate the ability of AuNP SERS to identify unique metabolites that can be potentially used as spectral biomarkers to differentiate FM and LC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Nonlocal Convection-Diffusion Problems on Bounded Domains and Finite-Range Jump Processes

In this paper, a nonlocal convection-diffusion model is introduced for the master equation of Markov jump processes in bounded domains. With minimal assumptions on the model parameters, the nonlocal steady and unsteady state master equations are shown to be well-posed in a weak sense. Finally, then the nonlocal operator is shown to be the generator of finite-range nonsymmetric jump processes and, when certain conditions on the model parameters hold, the generators of finite and infinite activity Lévy and Lévy-type jump processes are shown to be special instances of the nonlocal operator.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Diffusion Blurring of GaAs Quantum Wells Grown at Low Temperature

The processes of the diffusion blurring of a periodic system of GaAs quantum wells separated by AlGaAs barriers are studied by photoluminescence spectroscopy. The system is grown by molecular-beam epitaxy at a low temperature (200°C) and additionally doped with Sb and P isovalent impurities. Postgrowth annealing at the temperature 750°C for 30 min induces an increase in the energy corresponding to the photoluminescence peak of the e1–hh1 exciton state in quantum wells because of blurring of the epitaxial GaAs/AlGaAs interfaces due to enhanced Al–Ga interdiffusion in the cation sublattice. For the Al concentration profile defined by linear diffusion into quantum wells, the Schrödinger equation for electrons and holes is solved. It is found that the experimentally observed energy position of the photoluminescence peak corresponds to the Al–Ga interdiffusion length 3.4 nm and to the effective diffusion coefficient 6.3 × 10{sup –17} cm{sup 2} s{sup –1} at the temperature 750°C. This value is found to be close to the corresponding value for GaAs quantum wells grown at low temperatures without additional doping with Sb and P impurities. From the results obtained in the study, it is possible to conclude that enhanced As–Sb and As–P interdiffusion in the anion sublattice only slightly influences the processes of Al–Ga interdiffusion in the cation sublattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters

Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Theoretical analysis of mixing in liquid clouds – Part IV: DSD evolution and mixing diagrams

Evolution of droplet size distribution (DSD) due to mixing between cloudy and dry volumes is investigated for different values of the cloud fraction and for different initial DSD shapes. The analysis is performed using a diffusion–evaporation model which describes time-dependent processes of turbulent diffusion and droplet evaporation within a mixing volume. Time evolution of the DSD characteristics such as droplet concentration, LWC and mean volume radii is analyzed. The mixing diagrams are plotted for the final mixing stages. It is shown that the difference between the mixing diagrams for homogeneous and inhomogeneous mixing is insignificant and decreases with an increase in the DSD width. The dependencies of the normalized cube of the mean volume radius on the cloud fraction were compared with those on normalized droplet concentration and found to be quite different. If the normalized droplet concentration is used, mixing diagrams do not show any significant dependence on relative humidity in the dry volume. The main conclusion of the study is that traditional mixing diagrams cannot serve as a reliable tool for analysis of mixing type.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Visualization, understanding, and mitigation of process-induced-membrane irregularities in gas diffusion electrode-based polymer electrolyte membrane fuel cells

Polymer electrolyte membrane fuel cells (PEMFC) show substantial promise for their application in electric vehicles. For large-scale manufacturing of PEMFCs, roll-to-roll coated gas-diffusion-electrodes (GDE) offer certain advantages over other production pathways. Procedures including hot pressing and coating an ionomer overlayer may be necessary for this manufacturing pathway to enable a suitable catalyst layer/membrane interface. The same procedures may potentially introduce membrane irregularities, especially when thin membranes are used. Limited understanding exists regarding if and to what extent such irregularities impact PEMFC performance and lifetime, and therefore be considered defects. In this study, NREL's customized fuel cell hardware that enables quasi in-situ infrared (IR) thermography studies was utilized to visualize spatial hydrogen crossover and identify membrane irregularities that originated from the GDE-based MEA fabrication process. The structure of these membrane irregularities was investigated by scanning electron microscopy (SEM) and its impact on initial H2/air performance was determined. Accelerated stress testing (AST) revealed that these irregularities develop into failure point locations. These results were validated across many MEAs with identified process-induced membrane irregularities. By selecting specific gas diffusion media properties and by fine tuning the MEA hot pressing parameters, the formation of such membrane irregularities was mitigated.

08 HYDROGEN↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2025

Carbon mineralization pathways in interfacial adsorbed water nanofilms

Carbon mineralization in humidified carbon dioxide offers a promising route to mitigate anthropogenic emissions in a world stressed by water security. Despite its technological importance, our understanding of carbonation in water-poor environments lags, as traditional dissolution-precipitation pathways struggle to explain the adsorbed water nanofilm-mediated reactivity. Here, we utilize in operando X-ray diffraction (XRD) and advanced molecular simulations to investigate nanoconfined reactions driving forsterite carbonation, the magnesium-rich olivine. By examining magnesium ion dissolution and transport in atomistic simulations of the forsterite-water-carbon dioxide interface and comparing these with the in operando XRD activation energies, we identify both processes as rate-limiting at saturation. Our simulations reveal a mechanistic view of interfacial carbonation, where dissolution and precipitation are mediated by anomalous quasi two-dimensional diffusion. The transport process involves intermittent diffusive hopping in the desorbed state, separated by crawling events that are spatially short but temporally long. This understanding transcends carbon mineralization, with implications for understanding the transport of contaminants in geosystems, the design of multifunctional materials, water desalination, and molecular recognition systems.

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

Ultrafast (1‐5 sec) Lamination of Perovskite Solar Cells With Self‐Encapsulation Using Rapid Joule Heating

Perovskite solar cells (PSCs) are traditionally fabricated using sequential layer‐by‐layer deposition, in which each layer of the device is processed on top of the preceding layer. This constrains the processing techniques and selection of transport layer materials that can be used in the solar cell. To overcome these challenges, two half‐cells can be processed independently and then diffusion‐bonded through a lamination process. However, current lamination processes for perovskite solar cells suffer from relatively long process times, which can limit throughput when moving toward high‐volume manufacturing. In this study, a custom platform was designed for rapid‐joule heating of perovskite materials and devices. This enabled more than a 99% reduction in lamination time from 26 min to 1 s. Perovskite samples that were laminated in 1 s exhibited comparable values of percent bonded area, interfacial toughness, grain domain size, and X‐ray diffraction spectra to those laminated in greater than 10 min. As a proof‐of‐concept, 18.3% efficient devices were successfully laminated in 5 s. A transient heat transfer model was developed to describe the relationship between the perovskite temperature and the electrical power supplied to the heaters, establishing a baseline for predicting processing conditions in large‐scale manufacturing systems. Ultra‐fast lamination provides a pathway toward scalable roll‐to‐roll or sheet‐to‐sheet manufacturing of PSCs.

heat transfer↗
DOE OSTI2022

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. While simulations have been able to capture key aspects of these processes to reproduce film morphology, these models typically ignore any subsurface diffusion that occurs. The results of these diffusion and intermixing processes are often observed experimentally, but cannot be explained by invoking bulk lattice diffusion at the deposition temperature. In this work, we present quantitative experimental data that reveals anion and cation intermixing over long length scales during the deposition of epitaxial Fe 2 O 3 and Cr 2 O 3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18 O and/or 57 Fe, and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. These reveal that adatoms on the film surface act to “pull up” subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. A simple model incorporating these events qualitatively agrees with the experimental intermixing results. In addition to film deposition, these intermixing mechanisms may to be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Mapping the Temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, IR-IR photodissociation mass spectrometry to the D3O+?(HDO)(D2O)19 isotopologue of the “magic number” protonated water cluster, H+?(H2O)21. The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90-120 K). Starting at 100 K, large “jumps” are observed between two OH frequencies separated by ~300 cm-1 indicating migration of the OH group from the bound OH site at 3350 cm-1 to the free position at 3686 cm-1. Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Yang, Nan↗
DOE OSTI2018

Spectroscopic and Computational Analysis of Chemical Pollutants Interactions on Steel - 18409

There is limited understanding of the scientific interaction of the chemical contaminants to DOE facility materials which can alter the stability of the materials, affecting the general safety. The resulting studies will have a direct impact towards developing an enhanced cost-effective decontamination method and understanding the negative effects caused by the chemical waste on the structural integrity of the solid surface in contact. First order kinematic models provide a better fit to the experimental data, implying that one layer of Hg adsorbs to the steel surface. Energy dispersive X-ray analysis confirmed that stainless steel samples are mainly composed of Fe, Ni, and Cr, and have a low carbon content. Objective: To determine the physiochemical interactions of Mercury (Hg) and Cesium (Cs) to steel and concrete using spectroscopic, adsorption and computational studies. To develop a multiscale simulation that shows the effects of adsorption and diffusion of alkali metals on the mechanical performance of steel alloys used in nuclear power plants. Approach: Molecular Dynamics (MD) simulations were used to compute the stress gradient generated in steel due to Hg adsorption. Finite element models were developed for simulating the effects of temperature and pressure in a diffusion analysis. Based on experimental data, one layer of Hg was used for the MD simulations. Material properties at the atomic level were determined from the MD simulations The stress gradient from the MD simulations is used as an input parameter in the FE simulation. Conclusions and Future Work: Computational methods such as MD and FE simulations were used to determine the effects of mercury adsorption in stainless steel samples. The adsorption process induces a stress gradient in the facility materials, which can affect the mechanical performance of facility materials under service. Tensile stresses help the diffusion process while compressive have the opposite effect. Temperature has a positive effect in the mass diffusion process. Future work will evaluate the effects of temperature on the material properties at the molecular level using both experimental and MD methods. Future work will evaluate the mechanical performance of facility materials at a wide range of environmental variables.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Water-Mediated Surface Diffusion Mechanism Enables the Cold Sintering Process: A Combined Computational and Experimental Study

The cold sintering process (CSP) densifies ceramics at much lower temperatures than conventional sintering processes. Several ceramics and composite systems have been successfully densified under cold sintering. For the grain growth kinetics of zinc oxide, reduced activation energies are shown, and yet the mechanism behind this growth is unknown. In this work, we investigate these mechanisms in more detail with experiments and ReaxFF molecular dynamics simulations. We investigated the recrystallization of zinc cations under various acidic conditions and found that their adsorption to the surface can be a rate-limiting factor for cold sintering. In conclusion, our studies show that surface hydroxylation in CSP does not inhibit crystallization; in contrast, by creating a surface complex, it creates an orders of magnitude acceleration in surface diffusion, and in turn, accelerates recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗