Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase behavior”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Three-dimensional phase evolution and stress-induced non-uniform Li intercalation behavior in lithium iron phosphate

A lithium iron phosphate (LiFePO 4 ) olivine cathode displays rich phase transition behavior during charging/discharging. Despite extensive studies, the evolution of the LiFePO 4 /FePO 4 two-phase coexistence structure and the underlying mechanism are still a point of debate due to the intricate and often disparate phase morphologies observed. Here we apply phase-field simulations to provide detailed insights into the three-dimensional phase evolution process in LiFePO 4 upon delithiation under a large driving force. The simulation reveals that the coherency stress arising from the LiFePO 4 /FePO 4 lattice mismatch destabilizes the initially flat delithiation front and induces filamentary growth of the FePO 4 phase along [010], which is in excellent agreement with a recent operando X-ray imaging study. The combined simulation and experimental results provide the first illustration of the phenomenon of stress-induced instability of the lithium (de)intercalation front in battery electrode compounds. We show that this phenomenon is facilitated by the surface mode of coherent spinodal decomposition in LiFePO 4 and influenced by the anisotropies of misfit strain, elasticity and lithium diffusion. It is detrimental to battery performance and life by causing stress concentration and reducing the lithium (de)intercalation kinetics and thus should be mitigated.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Performance enhancement of aqueous ionic liquids with lower critical solution temperature (LCST) behavior through ternary mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid–liquid phase separation when mixed with water and heated above a lower critical solution temperature (LCST), resulting in a water-rich (WR) and an IL-rich (ILR) phase. These binary IL–water mixtures can be employed in a variety of thermodynamic processes such as forward osmosis (FO) desalination, for which two solution properties are desirable: low phase separation temperature and high osmotic strength (osmolality). However, these two properties are interlinked, with ILs that exhibit higher osmotic strengths typically requiring higher phase separation temperatures. This behavior tends to arise from the hydrophilicity of the IL cations, which enhances osmotic strength while also elevating the phase separation temperature. In this work, we highlight a pathway to overcome this tradeoff by developing ternary IL mixtures (two ILs with varying cation hydrophilicity mixed with water), which lowers the phase separation temperature while maintaining and even enhancing the osmotic strength of the solution. We characterize the mixing behavior (osmolality, phase separation temperature, WR phase purity, and WR to ILR phase mass ratio) of four ILs as a function of their concentration in solution. We find that an enhancement of up to 81.6% in the osmolality with a concomitant reduction of up to 15.4% in the phase separation temperature can be achieved using this approach. The ternary mixture is also shown to improve the phase separation kinetics by nearly 95% compared to the binary mixture. Overall, this work highlights a new pathway to improve the performance of LCST ILs for water and energy applications.

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

Tailoring ductile-phase toughened tungsten hierarchical microstructures for plasma-facing materials

Biological materials, such as bones, nacre, etc., have been found to show attractive and unique combinations of stiffness, strength and fracture toughness. Research has been conducted in various engineering areas to mimic the naturally hierarchical microstructures or nanostructures of these materials to produce materials with optimum stiffness, high strength and fracture toughness for their intended structural applications. In the same objective, this work applies a multiscale microstructural approach recently developed to investigate the deformation and fracture behavior of ductile phase toughened tungsten (W) materials such as tungsten-nickel/iron (W-Ni/Fe) composites that possess lamellar-like and hierarchical “brick-and-mortar” (BAM) microstructures. First, the approach is used to simulate tensile loading of W-Ni/Fe specimens cut out from hot-rolled W-Ni/Fe plates which possess lamellar-like microstructures. The finite element model of the gage section of the specimen consists of a W-Ni/Fe dual-phase microstructural domain in which the constitutive behaviors of W and Ni/Fe phases are described by an elastic-plastic damage model. The predicted material stress-strain response and crack pattern development as a function of loading are then compared to the corresponding experimental results to determine the constitutive model parameters for W and Ni/Fe. Subsequently, the model parameters are used to analyze W-Ni/Fe BAM microstructures that are artificially created to investigate the effects of microstructural features and morphologies on the composite stress-strain response, damage, and fracture patterns. The modeling shows that regular BAM microstructures that experience bridging mechanism combined with crack penetration across W-phase regions exhibit significantly higher strengths than the rather random lamellar-like microstructures. More significantly, however adjusting the brick’s length-to-height ratio, the BAM microstructure can be designed to allow a more distributed damage zone which leads to increased strength, ductility and fracture energy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Defect behavior and radiation tolerance of MAB phases (MoAlB and Fe2AlB2) with comparison to MAX phases

MAB phases are a new class of layered ternary materials that have already shown a number of outstanding properties. Here, we investigate defect evolution and radiation tolerance of two MAB phases, MoAlB and Fe 2 AlB 2 , using a combination of experimental characterization and first-principles calculations. We find that Fe 2 AlB 2 is more tolerant to radiation-induced amorphization than MoAlB, both at 150 °C and at 300 °C. The results can be explained by the fact that the Mo Frenkel pair is unstable in MoAlB and as a result, irradiated MoAlB is expected to have a significant concentration of Mo Al antisites, which are difficult to anneal even at 300 °C. We find that the tolerance to radiation-induced amorphization of MAB phases is lower than in MAX phases, but it is comparable to that of SiC. However, MAB phases do not show radiation-induced cracking which is observed in MAX phases under the same irradiation conditions. This research suggests that MAB phases might be a promising class of materials for applications that involve radiation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A review of the metastable omega phase in beta titanium alloys: the phase transformation mechanisms and its effect on mechanical properties

Since its discovery in 1954, the omega (ω) phase in titanium and its alloys has attracted substantial attention from researchers. The β-to-ω and ω-to-α phase transformations are central to β-titanium alloy design, but the transformation mechanisms have been a subject of debate. With new generations of aberration-corrected transmission electron microscopy and atom probe tomography, both the spatial resolution and compositional sensitivity of phase transformation analysis have been rapidly improving. Furthermore, this review provides a detailed assessment of the new understanding gained and related debates in this field enabled by advanced characterization methods. Specifically, new insights into the possibility of a coupled diffusional-displacive component in the β-to-ω transformation and key nucleation driving forces for the ω-assisted α phase formation are discussed. Additionally, the influence of ω phase on the mechanical properties of β-titanium alloys is also reviewed. Lastly, a perspective on open questions and future direction for research is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Increased Defect Resistance and Ordering in MnBi 2 (Se 1– x Te x ) 4 via Accurate Diffusion Monte Carlo

Stabilizing materials and controlling defect formation remain key challenges in materials science, particularly for theory, where small energy differences must be resolved for accurate predictions. Here, we applied state-of-the-art theoretical methods to topological materials, focusing on MnBi 2 Te 4 (MBT), which is a promising intrinsic magnetic topological insulator. Antisite defects in MBT alter its electronic structure and magnetism, degrading topological properties and causing experimental inconsistencies. Using diffusion Monte Carlo and density functional theory, we investigated the thermodynamic stability and defect formation in MBT, MnBi 2 Se 4 (MBS), and MnBi 2 (Se 1–x Te x ) 4 . We found that MnBi 2 Se 2 Te 2 can be stable at finite temperatures, with higher defect formation energies due to stronger Mn–Se bonding and reduced internal strain. Se preferentially substitutes Te near Mn instead of Te in the outer layer, encouraging long-range ordering when incorporated. For MnBi 2 (Se 1–x Te x ) 4 , cluster expansion phase diagrams revealed solid solution behavior when x <0.5 and phase separation for larger x. MBT and MBS are topological insulators; therefore, the MnBi 2 (Se 1–x Te x ) 4 family could offer tunable topological behavior and improved stability.

MnBi2Te4↗
DOE OSTI2024

Paradigm for approaching the forbidden spontaneous phase transition in the one-dimensional Ising model at a fixed finite temperature

The Ising model describes collective behaviors such as phase transitions and critical phenomena in various physical, biological, economical, and social systems. It is well known that spontaneous phase transition at finite temperature does not exist in the Ising model with short-range interactions in one dimension. Yet, little is known about whether this forbidden phase transition can be approached arbitrarily closely—at fixed finite temperature. Here I use symmetry analysis of the transfer matrix to reveal the existence of spontaneous ultranarrow phase crossover (UNPC) at finite temperature in one class of one-dimensional Ising models on decorated two-leg ladders, in which the crossover temperature T 0 is determined solely by on-rung interactions and decorations, while the crossover width 2 δ T is independently, exponentially reduced ( δ T = 0 means a genuine phase transition) by on-leg interactions and decorations. These findings establish a simple ideal paradigm for realizing an infinite number of one-dimensional Ising systems with spontaneous UNPC at desirable T 0 , which would be characterized in routine laboratory measurements as a genuine first-order phase transition with large latent heat thanks to the ultranarrow δ T (say less than one nanokelvin), paving a way to push the limit in our understanding of phase transitions and the dynamical actions of frustration arbitrarily close to the forbidden regime. Published by the American Physical Society 2024

1-dimensional systems↗
DOE OSTI2025

Engineering Electrical Transport by Implantation‐Induced Defects in CrN Films Without Affecting Thermal Conductivity

The transport properties of CrN thin films deposited on sapphire have been tailored through structural modifications induced by cumulative argon implantation. As-grown samples experience the typical structural transition in CrN films from orthorhombic at low temperature to cubic above the Néel temperature (≈280 K) and exhibit a metallic-like conduction in both phases. With increasing implantation dose, the conduction mode shifts to a semiconductor-like behavior in both phases, albeit at different damage levels. Analysis of the results suggests that hopping conduction becomes dominant beyond a given damage threshold. The results highlight a promising correlation between defect engineering and conduction mechanisms, offering valuable insights into the versatile electrical properties of CrN films. These implantation-induced defects scatter carriers, leading to a decrease in their mobility. As the implantation dose increases, the defect landscape evolves, modifying the density of states. However, up to a dose of 0.050 dpa, no significant influence on phonon scattering is observed. This approach demonstrates that ion implantation enables precise tuning of CrN's electrical properties without affecting thermal conductivity, offering valuable insights into defect engineering in transition metal nitrides and underscoring its potential for transport properties decorrelation.

defects↗
DOE OSTI2025

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗
DOE OSTI2024

Understanding the High-Temperature Deformation Behaviors in Additively Manufactured Al6061+TiC Composites via In Situ Neutron Diffraction

Aluminum matrix composites (AMCs) are designed to enhance the performance of conventional aluminum alloys for engineering applications at both room and elevated temperatures. However, the dynamic phase-specific deformation behavior and load-sharing mechanisms of AMCs at elevated temperatures have not been extensively studied and remain unclear. Here, in situ neutron diffraction experiments are employed to reveal the phase-specific structure evolution of additively manufactured Al6061+TiC composites under compressive loading at 250 °C. It is found that the addition of a small amount of nano-size TiC significantly alters the deformation behavior and increases the strength at 250 °C in comparison to the as-printed Al6061. Unlike the two-stage behavior observed in Al6061, the Al6061+TiC composites exhibit three stages during compression triggered by changes in the interphase stress states. Further analysis of Bragg peak intensity and broadening reveals that the presence of TiC alters the dislocation activity during deformation at 250 °C by influencing dislocation slip planes and promoting dislocation accumulation. These findings provide direct experimental observations of the phase-specific dynamic process in AMCs under deformation at an elevated temperature. The revealed mechanisms provide insights for the future design and optimization of high-performance AMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Influences of Mg Ions and Al Ions on Deposition Behavior of Supersaturated Silicic Acid - 20261

Silicic acid can change into various species such as deposit, colloidal, and soluble form under certain conditions: for example, gradual pH changes around geological disposal repositories. In particular, silicic acid deposits narrow groundwater flow paths, thereby restricting the flow of groundwater. Such effects may contribute to a delay in radionuclide migration. Previous studies have reported that silicic acid supersaturated in sodium and calcium ions can narrow the flow paths. This study examined the effect of magnesium (Mg) and aluminum (Al) ions on the deposition of supersaturated silicic acid on the solid phase. The deposition behavior of supersaturated silicic acid on the solid phase was investigated in a batch system while changing the initial supersaturation concentration, the coexisting ion concentration, and the amount of added solid phase. In the presence of Mg ions, silicic acid precipitation increased with increasing Mg ion concentration due to the influence of the electrolyte, and crosslinking reactions. The coexistence of Al ions significantly suppressed silicic acid deposition due to the formation of a covering on the solid surface during Al ion hydrolysis under conditions of high Al concentration (1.4 mM). On the basis of these experimental results, the apparent deposition rate constant of silicic acid was estimated to be in the range from 10{sup -11} to 10{sup -10} m/s. The Damkoehler number (the dimensionless number expressed by the ratio of the characteristic time of mass transfer and the chemical reaction constant) was sufficient to cause the clogging of the host rock by deposition. This suggests that supersaturated silicic acid deposition, regardless of the presence of Mg ions (and Al ions), contributes to improving barrier performance against radionuclide migration in the host rock. (authors)

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

Role of spin-orbit coupling on crystal-field splitting and phase-stability of rare-earth based layered intermetallic

Layered and incommensurate heterostructures have attracted much attention for the occurrence of superconductivity and charge density waves with the possibility of intercalating foreign atoms. However, lanthanide-based heterostructures where (Eu x Sm 1-x )S and TaS 2 are alternatively stacked have been scarcely investigated. In this work, we performed phase stability, bonding behavior and electronic-structure analysis of (Eu x Sm 1-x )TaS 3 using first-principles density-functional theory methods. Our phase stability analysis suggests 50 at.% solubility of Eu in (Eu x Sm 1-x )TaS 3 compared to Eu solubility in all proportions in cubic SmS. The instability of Eu beyond 50 at.% in (Eu x Sm 1-x )TaS 3 was attributed to higher density of Eu-4f states at the Fermi-level. Based on band position calculated from spin-orbit coupling effect, we constructed a qualitative schematic of possible crystal-field analysis. The local change in bond-length and bond-angle around Eu-S in (Eu x Sm 1-x )TaS 3 correlate well with our crystal-field analysis. We believe that quantum mechanical insights provided in this work will be useful to understand other complex heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LDRD conclusion poster - Synthesizing Heterometallic Uranium Single Crystals to Understand the Influence of the Secondary Metals on Uranyl Axial Bond Strength

Understanding how transition metals influenced the chemistry of lanthanide and actinide (f-element) materials is critical for advancing separation technologies, materials design, and coordination chemistry. This project examined how incorporating first-row transition metals affected the structural and spectroscopic properties of f-element coordination polymers. In uranium(VI)-based systems synthesized with 2,6-pyridinedicarboxylic acid (PDC) ligands, single-crystal X-ray diffraction and Raman spectroscopy revealed that the presence of transition metals shortened the uranyl axial bond and induced a blue shift in its symmetric stretching vibration—evidence of increased bond strength. Electronic structure analysis, including Density of States (DOS) calculations using density functional theory (DFT), revealed altered orbital overlaps and highlighted the role of transition metal d-orbitals in modulating bonding. Raman modes were modeled using truncated structural fragments in collaboration with the University of Notre Dame, and although the predicted frequencies were lower than experimental values, they remained within expected ranges. In parallel, similar experiments with cerium (Ce) in the presence of cobalt (Co) and PDC demonstrated multi-step single-crystal-to-single-crystal transformations—behavior not observed in the uranium systems. Initial products included light yellow, orange, and polycrystalline materials. Single-crystal X-ray diffraction studies, conducted in collaboration with the Colorado School of Mines, identified the yellow phase as monometallic Ce(PDC)2(H2O)2·4H2O and the orange phase as heterometallic Ce2Co(PDC)4(H2O)6. After standing in solution for one week, both phases fully transformed into a dark yellow crystalline phase, [Ce3(PDC)5(H2O)8].6(H2O). Remarkably, this transformation was reversible—disturbing the equilibrium by removing some crystals caused reversion to the initial Ce(PDC)2(H2O)2·4H2O phase, highlighting dynamic behavior. All three structures were previously unreported. Solid-state UV-visible and Raman spectroscopy further distinguished these phases, revealing ligand-to-metal charge transfer involving Ce and characteristic d–d transitions from Co(II). The precise mechanism driving these transformations remained unclear; however, pH-dependent experiments confirmed that the transformation did not occur when the pH decreased. Overall, the project demonstrated that transition metals could be employed to tune bonding interactions, structural dimensionality, and optical properties in f-element materials, establishing new pathways for designing functional heterometallic systems. The work resulted in several novel structural discoveries and fostered productive collaborations with the University of Notre Dame and the Colorado School of Mines.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2025

Uptake and binding of La 3+ and Cr 3+ ions by alkali metal substituted alpha-zirconium phosphate

The conversion of alpha-zirconium phosphate, Zr(HPO 4 ) 2 ·H 2 O, to the K-, Rb-, and Cs-phases and the subsequent ion exchange behavior of these alkali metal phases, A-ZrP, with La 3+ and Cr 3+ have been investigated. The conversion to the A-ZrP phases was achieved by reaction with metal chloride, metal hydroxide solution, and confirmed through various techniques, including X-ray powder diffraction, thermogravimetry analysis, IR spectroscopy, scanning and transmission electron microscopy, and X-ray photoelectron spectroscopy. The effect of material hydration, which increased the interlayer distance, was also examined. The ion exchange behavior of the A-ZrP showed strong affinity for both La 3+ , a representative for trivalent lanthanide metals, and Cr 3+ , a representative for the trivalent transition metals, with a rapid, near quantitative removal of the M 3+ ions at a pH of 3 when the ion concentration was ≤25% of the ion exchange capacity of the materials. Additionally, the ion affinity was shown to be pH and concentration-dependent, decreasing with a decrease in pH or an increase in ion concentration. Lastly, the binding strength of La- and Cr-loaded ZrP materials was examined through a series of leaching experiments in a carbonate buffer, phosphate buffer, HEPES buffer, and a series of EDTA solutions at various concentrations. No observable leaching occurred in the carbonate buffer, phosphate buffer, HEPES buffer, or when the EDTA concentration was below 0.1 mM. These results highlight the potential for the A-ZrP materials to provide a platform for trivalent transition metal and trivalent lanthanide radionuclide capture for various applications.

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

Generic character of charge and spin density waves in superconducting cuprates

Charge density waves (CDWs) have been observed in nearly all families of copper-oxide superconductors. But the behavior of these phases across different families has been perplexing. In La-based cuprates, the CDW wavevector is an increasing function of doping, exhibiting the so-called Yamada behavior, while in Y- and Bi-based materials the behavior is the opposite. Here, we report a combined resonant soft X-ray scattering (RSXS) and neutron scattering study of charge and spin density waves in isotopically enriched La 1.8–x Eu 0.2 Sr x CuO 4 over a range of doping 0.07 ≤ x ≤ 0.20. In this work, we find that the CDW amplitude is temperature independent and develops well above experimentally accessible temperatures. Further, the CDW wavevector shows a nonmonotonic temperature dependence, exhibiting Yamada behavior at low temperature with a sudden change occurring near the spin ordering temperature. We describe these observations using a Landau–Ginzburg theory for an incommensurate CDW in a metallic system with a finite charge compressibility and spin-CDW coupling. Extrapolating to high temperature, where the CDW amplitude is small and spin order is absent, our analysis predicts a decreasing wavevector with doping, similar to Y and Bi cuprates. Our study suggests that CDW order in all families of cuprates forms by a common mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗