Search NASA⌕ Search

SEARCH · Search NASA

Results for “grain boundary mobility”

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.

90 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Effects of oxygen on the tensile properties and fracture behavior of the Ta-10 W refractory alloy at quasi-static and impact strain rates

This study investigates effects of pre-charged oxygen at levels up to ∼1750 ppm on the tensile properties of the refractory alloy Ta-10 W. Tests were conducted at temperatures to 1100 °C and strain rates from quasi-static (∼10−3 s−1) to impact (∼103 s−1). At all temperatures, ductility decreased with increasing oxygen concentration and could be visually identified as a ductile-to-brittle transition (DBT). Additionally, for a given oxygen concentration, ductility was lower at impact strain rate than at quasi-static. The DBT was accompanied by a change from ductile transgranular to more brittle modes, intergranular and transgranular cleavage. The oxygen concentrations associated with the DBT depend on test temperature and strain rate. However, the transitions themselves are not sharp and cannot be defined as occurring at a specific oxygen concentration. At impact strain rates, ductility is relatively insensitive to test temperatures in the range 300–1100 °C with similar embrittling effects of oxygen; however, oxygen has a more severe embrittling effect at room temperature. This behavior is reversed at quasi-static strain rates where the embrittling effect of oxygen is more severe at elevated temperatures than at room temperature. It appears this difference is related to the kinetics of oxygen diffusion to grain boundaries and crack tips. Specimens that were not pre-charged with oxygen but tensile tested in various partial pressures of oxygen underwent dynamic embrittlement during the test, with the degree of embrittlement increasing as the partial pressure of oxygen in the environment increased. Consistent with such a mechanism, ductility increased with increasing strain rate for a given oxygen partial pressure, which would be the reverse if only intrinsic dislocation mobilities were controlling.

Miller, Roger [ORNL] (ORCID:0000000172155851)↗
DOE OSTI2025

Beyond curvature-driven grain growth: Insights from fully anisotropic Monte Carlo Potts simulations

Grain boundary (GB) motion away from the center of curvature, termed anti-curvature behavior, has recently been observed in 3D experiments but is not predicted by classical grain growth theory. In this study, we investigate this behavior using a novel, fully anisotropic Monte Carlo Potts (MCP) model that incorporates both misorientation and inclination dependencies of GB energy. We perform 3D grain growth simulations with isotropic and anisotropic GB energies to explore the relationship between GB velocity and curvature. Contrary to the classical relation that velocity is a product of reduced mobility and mean curvature, we observe no consistent correlation between velocity and mean curvature for individual GBs, even under isotropic conditions, though there is correlation between the average velocity and curvature for some cases. The 3D simulations exhibit frequent anti-curvature motion regardless of GB energy anisotropy including with isotropic GBs, though larger curvatures occur with anisotropic GB energy functions that promote low-energy GBs. In 2D simulations, anti-curvature behavior only occurs with the anisotropic functions that favor low-energy GBs. This difference between 2D and 3D results suggests that anti-curvature behavior results in part from the increased freedom of motion intrinsic in 3D GB networks. Furthermore, our results support recent experimental observations that demonstrate that simple curvature-driven models are insufficient for describing GB migration in polycrystals.

Anti-curvature↗
DOE OSTI2021

Irradiation-enhanced diffusion and diffusion-limited creep in U 3 Si 2

U 3 Si 2 is an advanced fuel candidate due to its relatively high fissile density and attractive thermal properties. Compared to standard UO 2 fuel, there are significant data gaps for the thermophysical and thermomechanical properties of U 3 Si 2 . Point defect concentrations and mobilities under irradiation govern a number of important fuel performance properties, such as creep and fission gas release. In this work, we utilized density functional theory (DFT) data to inform a cluster dynamics framework to predict point defect concentrations in U 3 Si 2 under irradiation. Molecular dynamics (MD) simulations were used to examine the contribution of atomic mixing during ballistic cascades to diffusion, as well as the diffusivity of U and Si at grain boundaries. These atomic scale models for diffusivity were then used to inform a creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusional creep, and climb-limited dislocation creep. The model compares well against available experimental data and has been implemented in the BISON fuel performance code. Finally, a demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U 3 Si 2 in-reactor, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Metasomatic alteration of coarse‐grained igneous calcium‐aluminum‐rich inclusions from CK3 carbonaceous chondrites

Abstract We report on the primary and secondary mineralogies of three coarse‐grained igneous calcium‐aluminum‐rich inclusions (CAIs) (Compact Type A [CTA], Type B [B], and forsterite‐bearing type B [FoB]) from the Northwest Africa (NWA) 5343 (CK3.7) and NWA 4964 (CK3.8) carbonaceous chondrites, compare them with the mineralogy of igneous CAIs from the Allende (CV3.6) chondrite, and discuss the nature of the alteration processes that affected the CK and CV CAIs. The primary mineralogy and mineral chemistry of the CK3 CAIs studied are similar to those from Allende; however, primary melilite and anorthite are nearly completely absent. Although the secondary minerals identified in CK CAIs (Al‐diopside, andradite, Cl‐apatite, clintonite, forsterite, ferroan olivine, Fe,Ni‐sulfides, grossular, ilmenite, magnetite, plagioclase, spinel, titanite, and wadalite) occur also in the Allende CAIs, there are several important differences: (i) In addition to melilite and anorthite, which are nearly completely replaced by secondary minerals, the alteration of CK CAIs also affected high‐Ti pyroxenes (fassaite and grossmanite) characterized by high Ti 3+ /Ti 4+ ratio and spinel. These pyroxenes are corroded and crosscut by veins of Fe‐ and Ti‐bearing grossular, Fe‐bearing Al,Ti‐diopside, titanite, and ilmenite. Spinel is corroded by Fe‐bearing Al‐diopside and grossular. (ii) The secondary mineral assemblages of grossular + monticellite and grossular + wollastonite, commonly observed in the Allende CAIs, are absent; the Fe‐bearing grossular + Fe‐bearing Al‐diopside ± Fe,Mg‐spinel, Fe‐bearing grossular + Fe,Mg‐olivine ± Fe,Mg‐spinel, and Ca,Na‐plagioclase + Fe‐bearing Al‐diopside + Fe‐bearing grossular assemblages are present instead. These mineral assemblages are often crosscut by veins of Fe‐bearing Al‐diopside, Fe,Mg‐olivine, Fe,Mg‐spinel, and Ca,Na‐plagioclase. The coarse‐grained secondary grossular and Al‐diopside often show multilayered chemical zoning with distinct compositional boundaries between the layers; the abundances of Fe and Ti typically increase toward the grain edges. (iv) Sodium‐rich secondary minerals, nepheline and sodalite, commonly observed in the peripheral portions of the Allende CAIs, are absent; Ca,Na‐plagioclase is present instead. We conclude that coarse‐grained igneous CAIs from CK3.7–3.8 s and Allende experienced an open‐system multistage metasomatic alteration in the presence of an aqueous solution–infiltration metasomatism. This process resulted in localized mobilization of all major rock‐forming elements: Si, Ca, Al, Ti, Mg, Fe, Mn, Na, K, and Cl. The metasomatic alteration of CK CAIs is more advanced and occurred under higher temperature and higher oxygen fugacity than that of the Allende CAIs.

Geochemistry & Geophysics↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2021

Pore-scale CFD simulations of clay mobilization in natural porous media due to fresh water injection

The present work investigates mechanisms of permeability impairment as a result of low-salinity fluid injection into brine-saturated porous media containing dispersible clays. Here, we present a computational fluid dynamics model at the pore-scale to simulate detachment, migration and straining of fine particles in porous media. The model uses an immersed boundary method to simulate the motion of clay fines in a fluid. In addition to the hydrodynamic forces, we model the Derjaguin-Landau-Verwey-Overbeek forces (DLVO) between clay fines and grains. Our simulations show the impact of the injected fluid's salinity and velocity on the concentration of clay fines retained on the grain surface. We see clay particles dislodging from the grain surface in clusters of up to 12 particles. Our simulation results also demonstrate clogging of narrow pore spaces by the detached particles from the grain surfaces.

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

Behavior of helium cavities in ion-irradiated W-Ni-Fe ductile-phase toughened tungsten

This study reports on the distribution of helium (He) cavities in a hot-rolled W-Ni-Fe ductile-phase toughened tungsten (DPT W) composite irradiated to a dose and a helium concentration that are comparable to those in the material after 5-year irradiation in a conceptual fusion power plant. Here, the DPT W sample consists of W particles embedded in a ductile-phase NiFeW matrix with a nominal composition of 90W-7Ni-3Fe by weight. It was hot-rolled to a thickness reduction by 87% (87R DPT W). Sequential irradiation of the material with 1.2 MeV Ni + ions to a fluence of 2.15×10 16 Ni + /cm 2 and 90 keV He + ions to 6.5×10 15 He + /cm 2 was performed at 973 K. Larger He cavities with a lower number density are observed in NiFeW than W. Helium cavities are aggregated preferentially along the NiFeW/W interphase boundary. This behavior is not observed along the W/W grain boundary under the same irradiation conditions. A similar distribution of He cavities with comparable sizes is observed in 87R DPT W irradiated with He + ions only, suggesting that the He cavity distribution is not significantly affected by pre-irradiation of Ni + ions at 973 K. The data from this study does not provide any evidence for formation of visible voids or Ni precipitates in W phase irradiated with Ni + ions only. Diffusion and clustering of He atoms and mobile vacancies during ion irradiation at 973 K are responsible for the formation and growth of the He cavities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unraveling the Unconventional Order of a High-Mobility Indacenodithiophene-Benzothiadiazole Copolymer

We report a new class of donor–acceptor (D-A) copolymers found to produce high charge carrier mobilities competitive with amorphous silicon (>1 cm 2 V –1 s –1 ) exhibit the puzzling microstructure of substantial local order, however lacking long-range order and crystallinity previously deemed necessary for achieving high mobility. Here, we demonstrate the application of low-dose transmission electron microscopy to image and quantify the nanoscale and mesoscale organization of an archetypal D-A copolymer across areas comparable to electronic devices (≈9 μm 2 ). The local structure is spatially resolved by mapping the backbone (001) spacing reflection, revealing nanocrystallites of aligned polymer chains throughout nearly the entire film. Analysis of the nanoscale structure of its ordered domains suggests significant short- and medium-range order and preferential grain boundary orientations. Moreover, we provide insights into the rich, interconnected mesoscale organization of this new family of D-A copolymers by analysis of the local orientational spatial autocorrelations.

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

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Machine-Learning-Guided Insights into Solid-Electrolyte Interphase Conductivity: Are Amorphous Lithium Fluorophosphates the Key?

Despite decades of study, the identity of the dominant Li + -conducting phase within the inorganic SEI of Li-ion batteries remains unresolved. While the mosaic model describes LiF/Li 2 O/Li 2 CO 3 nanocrystallites within a disordered matrix, these crystalline phases inherently offer limited ionic conductivity. Growing evidence suggests that interfaces, grain boundaries, and amorphous phases may instead host the primary fast-ion pathways. Using diffusion-based generative structure prediction and machine-learning interatomic potentials (MLIPs), we investigate lithium difluorophosphate (LiPO 2 F 2 ), a key mixed-anion decomposition product of phosphorus- and fluorine-containing electrolytes. We identify a stable crystalline polymorph and demonstrate that the amorphous counterpart is conductive, with projected room-temperature σ ≈ 0.18 mS cm –1 and E a ≈ 0.40 eV. Here, this enhancement stems from structural disorder flattening the Li site-energy landscape and a low formation energy for Li-interstitial defects, which supplies additional mobile carriers. We propose amorphous mixed-anion Li-P-O-F phases as a promising conducting medium in the SEI, offering a specific target for engineering improved battery interfaces.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2020

Water-Surface Drag Coating: A New Route Toward High-Quality Conjugated Small-Molecule Thin Films with Enhanced Charge Transport Properties

Electronic properties of organic semiconductor (OSC) thin films are largely determined by their morphologies and crystallinities. However, solution-processed conjugated small-molecule OSC thin films usually exhibit abundant grain boundaries and impure grain orientations because of complex fluid dynamics during solution coating. Here, a novel methodology, water-surface drag coating, is demonstrated to fabricate high-quality OSC thin films with greatly enhanced charge transport properties. This method utilizes the water surface to alter the evaporation dynamics of solution to enlarge the grain size, and a unique drag-coating process to achieve the unidirectional growth of organic crystals. Furthermore, using 2,8-difluoro-5,11-bis(triethylsilylethynyl)anthradithiophene (Dif-TES-ADT) as an example, thin films with millimeter-sized single-crystal domains and pure crystallographic orientations are achieved, revealing a significant enhancement (4.7 times) of carrier mobility. More importantly, the resulting film can be directly transferred onto any desired flexible substrates, and flexible transistors based on the Dif-TES-ADT thin films show a mobility as high as 16.1 cm 2 V -1 s -1 , which represents the highest mobility value for the flexible transistors reported thus far. The method is general for the growth of various high-quality OSC thin films, thus opening up opportunities for high-performance organic flexible electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical and Structural Evolution of AgCu Catalysts in Electrochemical CO 2 Reduction

Silver–copper (AgCu) bimetallic catalysts hold great potential for electrochemical carbon dioxide reduction reaction (CO 2 RR), which is a promising way to realize the goal of carbon neutrality. Although a wide variety of AgCu catalysts have been developed so far, it is relatively less explored how these AgCu catalysts evolve during CO 2 RR. The absence of insights into their stability makes the dynamic catalytic sites elusive and hampers the design of AgCu catalysts in a rational manner. Here, in this study, we synthesized intermixed and phase-separated AgCu nanoparticles on carbon paper electrodes and investigated their evolution behavior in CO 2 RR. Our time-sequential electron microscopy and elemental mapping studies show that Cu possesses high mobility in AgCu under CO 2 RR conditions, which can leach out from the catalysts by migrating to the bimetallic catalyst surface, detaching from the catalysts, and agglomerating as new particles. Besides, Ag and Cu manifest a trend to phase-separate into Cu-rich and Ag-rich grains, regardless of the starting catalyst structure. The composition of the Cu-rich and Ag-rich grains diverges during the reaction and eventually approaches thermodynamic values, i.e., Ag 0.88 Cu 0.12 and Ag 0.05 Cu 0.95 . The separation between Ag and Cu has been observed in the bulk and on the surface of the catalysts, highlighting the importance of AgCu phase boundaries for CO 2 RR. In addition, an operando high-energy-resolution X-ray absorption spectroscopy study confirms the metallic state of Cu in AgCu as the catalytically active sites during CO 2 RR. Taken together, this work provides a comprehensive understanding of the chemical and structural evolution behavior of AgCu catalysts in CO 2 RR.

Chen, Peng-Cheng↗
DOE OSTI2023

Degradation of electrical resistivity of tungsten following shielded neutron irradiation

A major challenge for heat transfer in nuclear materials is to ensure thermal mobility after high amounts of neutron irradiation. Tungsten is widely selected as a heat transfer material in fusion reactors. In metals, thermal conductivity is dominated by electrons’ ability to transfer energy. Neutron irradiation generates point defects, clusters, and solid transmutation (e.g.rhenium and osmium in tungsten), which inhibit electron motion. The purpose of this work is to quantify the irradiation-induced change in electron mobility and deconvolute transmutation and microstructural effects on observed changes to electron mobility. Single and polycrystalline tungsten were fast neutron irradiated in the High Flux Isotope Reactor at Oak Ridge National Laboratory to doses between 0.2 and 0.7 displacements per atom (dpa) and temperatures from 500 °C to 1000 °C. Grain growth was observed in all samples. Microstructure and transmutation were quantified. The geometric orientation of samples with elongated grains has been shown to affect electrical resistivity. A mathematical model was developed and used to deconvolute solid-solution transmutation, grain, and temperature-dependent lattice effects on resistivity. At ~0.4 dpa at ~590 °C, the combined resistivity degradation due to voids, vacancies, interstitials, and dislocations is estimated to be greater than the contribution from solid solution Re transmutation, which is greater than the contribution from grain boundaries. At doses of ~0.7 dpa at ~750 °C, solid solution Re contributions are greater than all other effects combined. As a result, this work establishes a basis to predict the effects of irradiation temperature and transmutation on thermal properties of tungsten and highlights the importance of irradiation temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct observation of annealing-driven recrystallization behavior in magnesium alloy at low strain condition

Grain-twin interactions are significant in texture modification under thermodynamic driving force. In this study, annealing-driven twinning/detwinning behavior, grain growth, and corresponding texture evolution in a pre-deformed AZ31 magnesium alloy were systematically tracked and investigated via in-situ heating synchrotron X-ray diffraction and quasi in-situ electron backscattered diffraction techniques. A twinning texture is generated in the pre-deformed sample due to the activation of {$10\bar{1}2$} tensile twinning. During annealing, dislocation annihilation occurs between 100 and 280 °C, and recrystallization occurs above 280 °C, manifesting as the initial residual matrix and twins being competitively swallowed by each other, forming a bimodal texture. The recrystallization process is completed by boundary movement, which depends on the energy difference across the boundary. In addition, it is found that the grain boundaries favor movement towards the side with higher stored energy, regardless of the boundary type or the boundary energy.

In-situ observation↗
DOE OSTI2024

Anisotropic weak antilocalization in thin films of the Weyl semimetal TaAs

Device applications of topological semimetals await the development of epitaxial films in the ultrathin limit. Weak antilocalization (WAL) has been extensively utilized in the understanding of surface states in topological insulators and shows promise for use in elucidating the properties of thin film topological semimetals. Here, we report insights from WAL in the surface state and interface transport properties of our recently synthesized single-crystal-like thin films of the Weyl semimetal TaAs(001) grown on GaAs(001). We observe robust, anisotropic WAL in the magnetoconductance from 2 to 20 K in films from 10 to 200 nm thick. We link the anisotropic WAL magnetoconductance to anisotropic mobility stemming from film topography. We conclude that WAL in the films likely originates from the antilocalization of topological surface states. The WAL magnetoconductance is impacted by the film thickness and topography, solidifying the useful role of WAL in the study of topological semimetal/semiconductor heterointerfaces. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The role of irradiation-enhanced interstitial diffusion in over-pressurizing fission gas bubbles in UO 2

Fission gas bubbles in UO 2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, in this work, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2022

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗