Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective diffusivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

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

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗
DOE OSTI2021

Measurement of the longitudinal diffusion of ionization electrons in the MicroBooNE detector

Accurate knowledge of electron transport properties is vital to understanding the information provided by liquid argon time projection chambers (LArTPCs). Ionization electron drift-lifetime, local electric field distortions caused by positive ion accumulation, and electron diffusion can all significantly impact the measured signal waveforms. This paper presents a measurement of the effective longitudinal electron diffusion coefficient, D L , in MicroBooNE at the nominal electric field strength of 273.9 V/cm. Historically, this measurement has been made in LArTPC prototype detectors. This represents the first measurement in a large-scale (85 tonne active volume) LArTPC operating in a neutrino beam. This is the largest dataset ever used for this measurement. Using a sample of ~70,000 through-going cosmic ray muon tracks tagged with MicroBooNE's cosmic ray tagger system, we measure D L = 3.74 +0.28 -0.29 cm 2 /s.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Using θ' interfaces as templates for planar L1 2 precipitation in AlCuMnZr alloys

Controlled Mn and Zr additions to Al-Cu alloys have allowed for the improved retention of mechanical properties after extended 350°C exposures by stabilizing the main strengthening θ' (Al 2 Cu) phase. Ultimately, θ'/L1 2 (Al 3 Zr) co-precipitate formation stabilizes θ' most effectively; however, Zr diffuses sluggishly and has low solubility in aluminum castings. Increasing the Zr segregation rate would allow for faster and more effective θ'/L1 2 co-precipitation. It is demonstrated that the Zr segregation rate is faster when the Zr matrix content is higher. A much higher Zr matrix content was achieved by rapid cooling during additive manufacturing (AM) that produces θ'/L1 2 co-precipitation faster, which is shown by scanning transmission electron microscopy and atom probe tomography experiments. It was also found that Zr continuously segregates to θ' interfaces up to the most aggressive heat treatment studied such that planar L1 2 precipitates remain after the metastable θ' dissolves. In this manner, we demonstrate that θ' coherent interfaces serve as perfect templates to form stable planar L1 2 precipitates that can provide strength at higher temperatures than traditional θ' strengthened AlCu alloys. This work introduces an alloy design strategy that uses metastable precipitates to quickly nucleate and grow co-precipitates with a desired geometry that contain slow diffusing elements. These ideas can be applied to engineer more heat resistant alloys by taking advantage of high solute matrix contents enabled by rapid cooling during additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shock-driven kinetic and diffusive mix in high-Z pusher ICF designs

Revolver and Double Shell Inertial Confinement Fusion capsule designs hope to achieve a robust volumetric thermonuclear burn via the use of a high-Z pusher shell filled with a cryogenic D–T fuel. Unfortunately, mix of the pusher material into the fuel (gas) may adversely impact the burn performance. Hydrodynamic instability of the metal/gas interface as the mix source is an obvious concern, but 1D effects may also be detrimental. Such effects include plasma diffusion at material interfaces, which has been the subject of numerous theoretical, computational, and experimental investigations. However, other 1D mix mechanisms may exist, which have yet to be thoroughly explored. In particular, plasma kinetic effects may drive the mix when a shock breaks out of the metal/gas interface. Using the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov–Fokker–Planck code, iFP, we show herein that shock-driven kinetic effects can reconfigure the interface and the interfacial width subsequently grows diffusively. Finally, we consider any implications for high-Z pusher designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Diffusion of alloying and fission atoms in α-U investigated by first-principles calculation

Stability and diffusion of small defect clusters containing alloying or fission atoms (“F” is used to represent both atoms) in an α-uranium (U) crystal are investigated using first-principles calculations. Here, results indicate that different stable defect structures are observed for an interstitial defect and a cluster (F n V m ) containing different alloying or fission atoms and vacancies. An interstitial defect is found to be most stable either at the center of a pyramid pentahedron or in a 〈010〉 mixed dumbbell configuration and the latter one is reported for the first time. Formation energy of a small F n V m cluster depends on the distribution of vacancies. Based on these stable defects, new migration energies and paths have been explored in this work. Anisotropic migration of an interstitial defect by jumping between nearest centers of two pyramid pentahedrons or jumping between two mixed U-Zr dumbbells has been suggested. Furthermore, a new migration-rotation mechanism has been explored for the first time for a FV 2 cluster, resulting in a 3D diffusion. Finally, the effect of substitutional Pu atoms on the migration of a FV 2 cluster suggests that optimizing the Pu concentration may prohibit the diffusion of some alloying or fission atoms, increasing the safe performance of metallic fuels.

3D diffusion↗
DOE OSTI2020

Phase mediated dynamics of self-propagating Co/Al nanolaminate reactions

The ignition of sputter deposited nanolaminate foils comprised of alternating Co and Al layers results in rapid, self-propagating formation reactions. The propagating waves present after ignition of 150 nm-thick foils are characterized in movie mode dynamic transmission electron microscopy where these are found to have reaction speeds and wave morphology that vary with bilayer thickness. High speed videography reveals different bilayer thickness-wave character relationships in 750 nm- and 7500 nm-thick Co/Al foils. The reaction speed dependencies on bilayer thickness are calculated for each total thickness by treating the effect of radiation loss as a perturbation from an analytical model described by the difference in the heat of reaction measured in calorimetry and the adiabatic heat of product formation. From this model, an effective activation energy, diffusion constant, and flame temperatures are obtained, which allows for an interpretation of the reaction phase variations with laminate design and their effects on propagating wave morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An evaluation of size-resolved cloud microphysics scheme numerics for use with radar observations Part II: Condensation and evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part I of this study addressed collision–coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision–coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Here, large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effect of Mg and Ni impurities on tritium diffusion in lithium ceramics through cluster dynamics simulations

This study investigates the impact of Mg and Ni doping on tritium diffusion in LiAlO 2 and LiAl 5 O 8 ceramics, that are used in tritium-producing burnable absorber rods (TPBARs). Utilizing Centipede simulations across a broad temperature range (500 K to 1250 K), we explore the interplay between defect dynamics, cluster formation, and tritium mobility. In LiAlO 2 , Mg doping significantly enhances tritium diffusivity by increasing tritium interstitial concentrations and diffusion coefficients of key species, thereby doubling the overall tritium diffusivity. Ni doping, while shifting the dominant defect to Li vacancies, maintains high tritium mobility due to the low binding energy of Li vacancy-tritium complexes, which ensures effective tritium migration. In LiAl 5 O 8 , Mg and Ni doping results in a slight reduction in the diffusion coefficients of key species, yet the dramatic increase in tritium interstitial concentrations compensates, leading to a net small increase in tritium diffusivity. In conclusion, the findings highlight the critical role of defects in tritium transport and the effect of Mg and Ni defects on the performance of these ceramics in demanding nuclear environments.

Cluster dynamics↗
DOE OSTI2021

Machine learning assisted modeling of mixing timescale for LES/PDF of high-Karlovitz turbulent premixed combustion

Accurate modeling of mixing in the transported probability density function (PDF) method remains a great challenge, especially for turbulent premixed combustion under extreme conditions such as high Karlovitz number Ka. Recently, a power-law based mixing timescale model was developed for the large-eddy simulations (LES)/PDF modeling of high-Ka number turbulent premixed flames. It is found in this work that the power-law mixing timescale model is highly sensitive to the model parameters. It is thus critically needed to develop accurate calibration of these model parameters. The empirical specification of the model parameters developed in Zhang et. al. is found to be inadequate for accurate modeling of the mixing timescale. Here, machine learning is introduced as an attractive alternative in this work for the specification of the model parameters. A high-Ka number DNS jet flame is used as the training and validation of the machine learning models. The choices of the input parameters are discussed and compared for the machine learning models. The effect of differential molecular diffusion on mixing is examined by including the effect of the Lewis number in the training of the machine learning models. The performance of different machine learning algorithms is compared for the specification of the mixing model parameters. Overall, excellent performance of the machine learning models is observed for assisting the mixing modeling. The feasibility, interpretability, applicability, generality, and portability of using machine learning are discussed in general to provide a perspective on applying data-driven machine learning for turbulent combustion modeling studies.

42 ENGINEERING↗
DOE OSTI2023

Late-time post-merger modeling of a compact binary: effects of relativity, r-process heating, and treatment of transport

Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first 10 s after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10 000 km, so that thermal wind velocities only asymptote beyond this radius and at higher values than often reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Active and passive defects in tetragonal tungsten bronze relaxor ferroelectrics

Tetragonal tungsten bronze (TTB) based oxides constitute a large family of dielectric materials which are known to exhibit complex distortions producing incommensurately modulated superstructures as well as significant local deviations from their average symmetry. The local deviations produce diffuse scattering in diffraction experiments. The structure as well as the charge dynamics of these materials are anticipated to be sensitive to defects, such as cation or oxygen vacancies. In this work, in an effort to understand how the structural and charge dynamical properties respond to these two types of vacancy defects, we have performed measurements of dielectric susceptibilities and single crystal diffraction experiments of two types of TTB materials with both ‘filled’ (Ba 2 NdFeNb 4 O 15 and Ba 2 PrFeNb 4 O 15 ) and ‘unfilled’ (Sr 0.5 Ba 0.5 Nb 2 O 6 ) cation sublattices. We also perform these measurements before and after oxygen annealing, which alters the oxygen vacancy concentrations. Surprisingly, we find that many of the diffuse scattering features that are present in the unfilled structure are also present in the filled structure, suggesting that the random fields and disorder that are characteristic of the unfilled structure are not responsible for many of the local structural features that are reflected in the diffuse scattering. Furthermore, oxygen annealing clearly affected both color and dielectric properties, consistent with a diminishment of the oxygen vacancy concentration, but had little effect on observed diffuse patterns.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Study of quasi-collisional effects in laboratory and astrophysical plasmas

High-amplitude turbulence excited in plasmas at small-scales can leave imprint on the radiation produced by (accelerated) plasma particles -- mostly electrons in electron-ion plasmas and both electrons and positrons in lepton pair plasmas. Furthermore, turbulence is known to introduce "effective collisional" effects -- anomalous resistivity, dissipation, diffusion, etc. -- in the otherwise collisionless plasmas. These effective collisions affect radiation transfer and modify optical and magneto-optic effects -- the transmittance and reflectivity coefficients, Faraday rotation, etc. Note that the effective collisions considered in this project are not resonant wave-particle interactions, but much less conventional randomization of particles’ paths. Therefore, we colloquially refer them to as "quasi-collisions". The study of quasi-collisions on transport, radiative and optical properties of plasmas is the main focus of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Measuring thermal diffusivity and gap conductance in uranium nitride and Zircaloy relevant for microreactor applications

Heat transfer across nuclear fuels and structural interfaces is an important factor for evaluating the performance of nuclear power systems. Specifically, heat generated as nuclear fuel fissions must be transported through the cladding material and through the reactor to reach the steam turbine for power generation. As new microreactor designs emerge, maximizing the efficiency of this heat transfer process becomes crucial to make them commercially viable. This article examines thermal diffusivity and gap conductance in uranium nitride (UN) fuel and Zircaloy-4 (Zry4) cladding using light flash analysis (LFA). Thermal diffusivity measurements were made on monolithic UN pellets and Zry4 exposed to carbon at peak operating temperatures of microreactors and show that carbon ingress has a minimal effect on thermal diffusivity when compared with identical materials not exposed to carbon. Evaluation of gap conductance at the UN-Zry4 interface was done using one-dimensional two-layer thermal transport models as a function of applied pressure. Here the results show that increasing pressure on the UN-Zry4 interface leads to gains in gap conductance per unit area in fuel-cladding assemblies at microreactor operating temperatures. While many other variables are expected to influence UN-Zry4 interfacial gap conductance (e.g. contact surface roughness, porosity, localized heating, environmental gas pressure), the work offers a demonstration of using a conventional LFA apparatus to determine this parameter at elevated temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

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

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗