Search NASA⌕ Search

SEARCH · Search NASA

Results for “triple phase boundary”

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.

51 records · Page 3

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗

Ab Initio Phase Diagram of Tungsten

The phase diagram of tungsten (W) to a pressure (P) of 2500 GPa is investigated using a comprehensive ab initio approach that includes (i) the calculation of the zero temperature (T) free energies (enthalpies) of different solid structures, (ii) the quantum molecular dynamics simulation of the melting curves of different solid structures, (iii) the derivation of the analytic form for the solid-solid phase transition boundary, and (iv) the simulations of the solidification of liquid W into the final solid states on both sides of the solid-solid phase transition boundary, in order to confirm the corresponding analytic form. There are two solid structures confirmed to be present on the phase diagram of W, the ambient body-centered cubic (bcc) and the high-pressure double hexagonal close-packed (dhcp). At T = 0, the bcc-dhcp transition occurs at 1060 GPa, and the transition boundary has a positive slope dT/dP : the bcc-dhcp-liquid triple point is at (P, T) = (1675 GPa, 23680 K).

36 MATERIALS SCIENCE↗

Atomistic and mesoscale simulations to determine effective diffusion coefficient of fission products in SiC

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) particles serves as the barrier to prevent escape of fission products produced in the fuel kernel. Knowing the diffusion coefficient of fission products through SiC is critical to determining whether fission gas can escape from the particle. It has been observed in experiments that Ag accumulated in grain boundaries and triple junctions in SiC. It is hypothesized that grain boundary diffusion is the primary pathway by which fission products penetrate the SiC layer. In this report, the effective diffusion coefficient of the fission product Ag through the grain boundary network is calculated using a combination of atomistic and phase-field methods. The grain boundary diffusion coefficient is calculated using molecular dynamics simulations. The bulk diffusion coefficient is determined using a combination of density functional theory and nudged elastic band methods. An effective diffusion coefficient is calculated, accounting for the grain structure using a phase-field method. The effective diffusion coefficient will be incorporated into Bison and fission product release calculations are compared to available experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Understanding grain refinement and intergranular gas bubble evolution in U-10Mo fuel using phase-field modeling

Monolithic U-10Mo fuel undergoes significant microstructural changes in the form of grain refinement and gas bubble formation during burnup, which degrades its mechanical properties. In this talk, I present a phase-field model for microstructure evolution in U-10Mo developed using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Simulations demonstrate that grain refinement initiates at pre-existing grain boundaries (GBs) to eliminate the lattice distortion energy caused by the accumulation of self-interstitial loops. By employing an equation of state for xenon gas, we simulate the evolution of gas bubbles in the polycrystal microstructure. Large, interconnected bubbles are found to form along the triple junctions. The effects of defect production rate, diffusivities and GB mobility on the microstructure evolution are systematically studied. Homogenization is employed on the microstructures to obtain effective elastic constants and diffusivity as a function of fission density. The simulations provide critical insights on microstructure and property degradation in U-10Mo fuel.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase-field modeling for restructuring in the dark zone of high burnup UO 2

This report summarizes the mesoscale modeling work performed in fiscal year 2024 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstructural evolution and restructuring observed in the dark regions of high burnup UO 2 nuclear fuel. This is the first attempt to realistically simulate the restructuring behavior observed in different region of a high burnup fuel. We employ a grand-potential based phase-field model to concurrently evaluate the formation of subgrains and growth of fission bubbles within the fuel. A energy-based subgrain formation criteria is introduced to simulate the restructuring process. Effect of different initial conditions and different modeling parameters are studies systematically to capture how each of these parameters influence the characteristics of the restructured fuel. It is observed that the subgrain formation begins around existing fission gas bubbles and then proceeds towards triple junctions, grain boundaries and grain interiors. It is demonstrated that restructuring is influenced by a combination of initial dislocation densities, subgrain formation rate, and temperature. Rate of restructuring increases with increase in fuel temperature. A restructuring bias is observed within the microstructure due to variation in defect accumulation among different grains. Furthermore, bubble sizes and distribution does not have a significant effect on rate of restructuring. The predicted microstructures resembles the characteristics of the restructured regions as observed in experiments. Finally, a correlation is presented that demonstrates the evolution of the restructuring volume fraction as a function of local effective burnup. This work provides a first of its kind restructuring model for darkzone that can be used by BISON for performance prediction of high burnup UO 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Electronic and magnetic phase diagrams of the Kitaev quantum spin liquid candidate Na 2 Co 2 TeO 6

The 3⁢d 7 Co 2+ -based insulating magnet Na 2 ⁢Co 2 ⁢TeO 6 has recently been reported to have strong Kitaev interactions on a honeycomb lattice and is thus being considered as a Kitaev quandum spin liquid candidate. However, due to the existence of other types of interactions, a spontaneous long-range magnetic order occurs. This order is suppressed by applied magnetic fields leading to a succession of phases and ultimately saturation of the magnetic moments. The precise phase diagram, the nature of the phases, and the possibility that one of the field-induced phases is a Kitaev quantum spin liquid phase are still a matter of debate. Here, in this study, we measured an extensive set of physical properties to build the complete temperature-field phase diagrams to magnetic saturation at 10 T for magnetic fields along the a and a* axes, and a partial phase diagram up to 60 T along c. We probe the phases using magnetization, specific heat, magnetocaloric effect, magnetostriction, dielectric constant, and electric polarization, which is a symmetry-sensitive probe. With these measurements, we identify all the previously incomplete phase boundaries and find additional high-field phase boundaries. We find strong magnetoelectric coupling in the dielectric constant and moderate magnetostrictive coupling at several phase boundaries. Furthermore, we detect the symmetry of the magnetic order using electrical polarization measurements under magnetic fields. Based on our analysis, the absence of electric polarization under zero or finite magnetic field in any of the phases or after any combination of magnetic/electric field cooling suggests that a zigzag spin structure is more likely than a triple-Q spin structure at zero field. Finally, we investigate the hysteresis and first- or second-order nature of each phase transition and its entropy changes. With this information, we establish a map of the magnetic phases of this compound and its magnetic, thermodynamic, and magnetoelectric properties, and discuss where spin liquid or other phases may be sought in future studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mesoscale modeling of restructuring in high burnup UO 2 fuel

Here, this work aims to simulate the restructuring behavior observed in different regions of high burnup fuel, providing a first-of-its-kind restructuring model for the dark zone and rim region of high-burnup UO 2 fuel. We employed a grand-potential-based phase-field model to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of different initial conditions and different modeling parameters were systematically studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to begin around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Restructuring was demonstrated to be influenced by a combination of initial dislocation densities, burnup rate, subgrain formation rate, and temperature. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. Furthermore, bubble size and distribution do not significantly affect the rate of restructuring. The predicted microstructures are consistent with experimental observations of the restructured regions. Finally, we present a correlation demonstrating the evolution of the restructuring volume fraction as a function of local burnup.

UO2↗

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗

Additive Manufactured Composite Phase-Change Material for Thermal Energy Storage Applications

Phase-change materials play a critical role in industrial energy storage applications to drive efficiency improvements, thermal energy management, and carbon emissions reductions. Recently, it has been shown that rapid solidification of alloys with metastable immiscibility in the liquid phase has the potential to form unique microstructures in which a low-melting phase is uniformly distributed in a high-melting matrix. This feature can be exploited using additive manufacturing to produce components with complex geometries containing such unique phase-change microstructures. Phase-field simulations utilizing high-performance computing were used to provide a detailed description of the evolution of the active phase during service in terms of their morphology and composition in different polycrystalline matrix grain morphologies that are typically produced during additive manufacturing. Phase field simulations were performed using, MEUMAPPS-SL (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid Liquid) code that was developed in-house by the Oak Ridge National Laboratory. The simulations utilized the capabilities of the Kestrel supercomputer at the National Renewable Energy Laboratory. The simulation results were compared with experimental results generated at Siemens Energy, Inc. The results indicate that the kinetics of liquid spreading along grain boundaries is largely determined by the mobility of the triple line along the intersection of the grain boundary liquid and the grain boundary plane.

25 ENERGY STORAGE↗

Electrical control of topological 3Q state in intercalated van der Waals antiferromagnet Co x -TaS 2

Van der Waals (vdW) magnets have opened a new avenue of opportunities encompassing various interesting phases. Co 1/3 TaS 2 –an intercalated metallic vdW antiferromagnet–is one of the latest additions to this growing list of materials due to its unique topologically nontrivial triple-Q (3Q) ground state. This 3Q tetrahedral structure, which critically depends on the Co content, yields the highest-density Skyrmion lattice with scalar spin chirality, resulting in a noticeable anomalous Hall effect. In this work, we demonstrate control of this topological phase via ionic gating. Using four Co x TaS 2 devices with different Co compositions, we show that ionic gating can cover the entire 3Q topological phase and reveal the nature of the thermodynamically inaccessible phase space. Another striking finding in our data is the existence of an adiabatic discontinuity in the phase boundary between the 3Q and 1Q phases. Our work constitutes one of the first examples of electrical control of scalar spin chirality using an antiferromagnetic metal.

Magnetic properties and materials↗

Assessment of effective elastic constants of U-10Mo fuel microstructures

Monolithic U-10Mo fuel undergoes significant microstructural changes during fuel burnup which degrades its mechanical properties. In this talk, we present results form a numerical method to assess the impact of the various microstructural features--grains, intragranular and intergranular Xe gas bubbles--on the elastic stiffness tensor. Using the Multiphysics Object-Oriented Simulation Environment (MOOSE), phase-field-based microstructures are combined with asymptotic expansion homogenization method to obtain effective elastic constants as a function of porosity and fission density. The results are verified and compared against analytical homogenization models. With this approach, elastic degradation in operating nuclear fuels can be quantified when the distributions of microstructural features are known from experimental characterization or rate-theory based models. We further develop an evolution model based on the virial equation of state for Xe gas and investigate the effect of growth and coalescence of the bubbles at the grain boundary faces and triple junctions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effects of Post-sintering Annealing on (NdLa)-(FeCo)-B Magnets

A common limitation of La substitutions into (Nd, La) 2 Fe 14 B magnets is the reduction of the anisotropy field (H a ), which decreases magnet coercivity (H cj ). In this study, adding a small amount of Pr-Cu to a La-containing neo magnet may modify the grain boundary (GB), which can help recover a fraction of the degraded H cj . An optimal GB modification requires a carefully designed post-sinter heat treatment, which is the focus of this work on characterizing the effects of multi-step post-sintering annealing on the evolutions of microstructure and magnetic properties. We find that Nd and Cu concentrations at GBs and triple junctions (TJs) increase when the annealing temperatures are lowered. Annealing temperatures of 580°C and 480°C (near rare-earth-rich eutectic temperatures) enabled the development of a thick and continuous GB phase, which helps to magnetically decouple the grains, and, thus, enhance H cj . For the alloy with 25% Nd replaced by La, and with the addition of 7.5 wt.% of Pr 68 Cu 32 as GB modifier, the optimized multi-step post-sinter annealing improved the H cj from 4.8 kOe to 9 kOe.

36 MATERIALS SCIENCE↗

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

Here we present a series of simple shear numerical simulations of dynamic recrystallization of two-phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (V p ) and fast S wave (V s ) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest V p quickly aligns with the maximum c-axis orientation. At the same time, the maximum c-axis development reduces V s in this orientation. When water is present, the developed maximum c-axis orientation is less intense, which results in lower V p and V s . At high percentage of water, V p does not align with the maximum c-axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of V s while V p remains close to the value for purely solid ice. These results suggest that the decrease in V s observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c-axis-dependent orientation and the magnitude of seismic anisotropy.

58 GEOSCIENCES↗

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗