Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibria”

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

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Experimental study of mechanistic factors influencing solvent-driven fractional crystallization of calcium sulfate

To advance dimethyl ether-driven fractional crystallization (DME-FC), a more sustainable method of water treatment and mineral recovery, a range of chemical equilibria were measured. These include varying concentrations of miscible organic solvents (MOS) used to experimentally measure the solvent-induced solid-liquid equilibrium (SLE) of calcium sulfate (CaSO 4 ) in water. Seven MOS, including dimethyl ether (DME), acetonitrile (MeCN), 1,4-dioxane, tetrahydrofuran (THF), acetone, ethanol, and diethylamine, were screened to establish trends associated with molecular volume, functional groups, and physical properties. The effect of MOS on CaSO 4 removal differed at concentrations <0.15 mol fraction MOS; MOS with greater molecular volume (THF, 1,4-dioxane, and diethylamine) induced greater CaSO 4 precipitation on a per mole basis. The solvent-induced SLE for all MOS converged between 0.15 and 0.2 mol fraction MOS, reaching a CaSO 4 concentration consistent with a water to MOS hydration ratio of 5:1 to 6:1, which may correspond to the solvent generating a solution-based pseudo-clathrate structure with continuity within the solution. Finally, solution pseudo-clathrate structures provide a mechanistic basis for DME-FC.

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

No detectable redox exchange between sulfur and iron during rapid cooling of basalts

The valence states of iron and sulfur in mantle-derived melts influence and respond to the composition of the mantle as well as the subsequent evolution of those melts in Earth's crust. Characterizing the valence states of iron and sulfur in silicate melts that have quenched to form glasses is therefore critical to understanding the composition of the mantle and the magmatic processes that take place before and during volcanic eruptions. Glasses quenched from silicate melts offer the opportunity to precisely measure sulfur and iron valence states, but whether the transition from melt to glass affects the valence state of sulfur and iron is uncertain. Here, we use the glassy margins of two pillow basalts, one from a mid-ocean ridge and one from a back-arc basin, as a natural experiment to test the effect of quench rate on sulfur and iron valence states. We report micro-X-ray Absorption Near Edge Structure (XANES) measurements of S 6+ /ΣS (S 6+ /[S 6+ +S 2- ]) and Fe 3+ /ΣFe (Fe 3+ /[Fe 2+ +Fe 3+ ]) along transects from the rapidly quenched rims of the pillows down into the slowly cooled crystal-rich interiors. The range of average quench rates estimated for our selected analysis areas range from 1.2 x 10 5 °C/sec at the pillow margin to 1.1 °C/sec in the interior and span quench rates experienced by natural and experimental glasses. On average, S6+/ΣS = 0.09 ± 0.01 and Fe 3+ /ΣFe = 0.147 ± 0.002 in the mid-ocean ridge sample and S 6+ /ΣS = 0.21 ± 0.01and Fe 3+ /ΣFe = 0.153 ± 0.004 in the back-arc sample, where the quoted uncertainty is the 1σ standard deviation of n= 9-16 analyses. We find that the S 6+ /ΣS and Fe 3+ /ΣFe values along transects in each sample show no systematic changes with depth, and by proxy, no systematic changes with quench rate. In each pillow, all measured S 6+ /ΣS and Fe 3+ /ΣFe ratios are within the absolute uncertainty of each technique. We propose that these observations are consistent with either sluggish reaction kinetics between sulfur and iron, or an insignificant temperature effect on sulfur-iron redox exchange equilibria. These observations demonstrate that sulfur and iron valence states in glasses with varying quench rates can be directly compared, and that changes in sulfur and iron valence during quench are below the limits of detection with XANES.

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

Thermodynamic Modeling of Aqueous LiCl, LiBr, LiI, and LiNO 3 Solutions

Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO 3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each electrolyte-molecule pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li + , while in the second case hydration of the lithium ion from Li + to Li(H 2 O) + is considered. Here, the calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.

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

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2022

Regimes of evaporation and mixing behaviors of nanodroplets at transcritical conditions

The objective of this paper is to examine the fundamental mechanisms responsible for the transition between subcritical evaporation and supercritical dense-fluid-mixing in the absence of convection effects, specifically focusing on the liquid–vapor interfacial dynamics. To isolate the dynamics of this transition process, we characterize the different physical behaviors exhibited by an -dodecane nanoscale droplet placed in different nitrogen ambient conditions across the fuel’s critical point. We employ a continuum-based interface-resolving diffuse-interface method to explore the underlying phase-exchange mechanisms that bring about such distinct dynamics. Following the comparison against molecular dynamics simulations and experiments of evaporating droplets and experimental data for vapor–liquid equilibria, a parametric study at various ambient conditions and droplet sizing is performed to identify four regimes of evaporation/mixing behaviors: sub- and supercritical droplet evaporation, and sub- and supercritical dense-fluid-mixing. It is shown that the distinction in the phase-exchange mechanisms in these four regimes are brought about by the different thermodynamic phases the droplet center can exhibit during the evaporation/mixing process: subcritical liquid, supercritical liquid-like, subcritical gaseous, and supercritical gas-like, respectively. It is shown that the subcritical dense-fluid-mixing behavior is a direct result of nanoconfinement of the liquid–vapor interfacial structure and thus is not present for large droplet sizes. Finally, the present study also shows that the supercritical phase-exchange dynamics can follow two different pathways: supercritical droplet-like evaporation and supercritical dense-fluid-mixing. Furthermore, promoting the early transition to supercritical dense-fluid-mixing can significantly expedite the phase-exchange process through the disintegration of the liquid-like droplet core.

42 ENGINEERING↗
DOE OSTI2025

Shadow masks predictions in SPARC tokamak plasma-facing components using HEAT code and machine learning methods

Here, this work uses machine learning (ML) to complement HEAT (Heat flux Engineering Analysis Toolkit) by developing 3-D footprint surrogate models for fast and accurate heat load calculations in the divertor of the SPARC tokamak. The focus is on shadowed regions, or magnetic shadows, caused by the 3-D geometry of plasma-facing components (PFCs). ML classifiers are employed to create a surrogate model for HEAT generated shadow masks, predicting these shadow masks and divertor heat flux profiles based on a diverse range of equilibria and only the plasma current, safety factor(q95) at the edge, and magnetic flux angles as input parameters. The ultimate goal is to integrate the model for real-time control and future operational decisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Disruption modelling for engineering and physics design of ST-E1 fusion power plant

Plasma disruptions represent a critical challenge for high-performance tokamak operations, as they can compromise machine integrity and reduce operational availability. Although future fusion devices essentially need to incorporate strategies to minimise disruption occurrence, complete avoidance remains unattainable. Consequently, assessing and characterising unmitigated disruption consequences is fundamental for the design and qualification of next-generation fusion power plants. This work supports the pre-conceptual design of ST-E1, a low aspect-ratio Tokamak Fusion Power Plant developed by Tokamak Energy Ltd., by presenting a comprehensive disruption modelling approach applied across different design stages. The methodology integrates both physics and engineering considerations to evaluate the impact of disruptions on machine performance and structural integrity. From an engineering perspective, several ST-E1 layout options were analysed to investigate the electromagnetic response of key components under disruption-induced loads, enabling comparison between alternative design solutions. On the physics side, a broad set of disruption scenarios was explored, scanning operational space parameters, plasma-material interactions, and associated thermal loads. Furthermore, the study examined variations in disruption behaviour arising from different reference equilibria, focusing on a range starting from Double Null to Single Null configurations, reflecting the increasing up-down asymmetry consequences. The results reveal significant contrasts in plasma dynamics and structures electromagnetic behaviour between configurations, highlighting the importance of disruption modelling in guiding design choices. These analyses have proven instrumental in shaping ST-E1 development, offering critical insights for mitigating risks and optimising future fusion power plant designs.

Borowiec, Katarzyna [ORNL] (ORCID:0000000335911739↗
DOE OSTI2023

Dual carbonate clumped isotopes (Δ 47 -Δ 48 ) constrains kinetic effects and timescales in peridotite-associated springs at the Cedars, Northern California

The Cedars is an area in Northern California with a chain of highly alkaline springs resulting from CO 2 -charged meteorological water interacting with a peridotite body. Serpentinization resulting from this interaction at depth leads to the sequestration of various carbonate minerals into veins accompanied by a release of Ca 2+ and OH – enriched water to the surface, creating an environment which promotes rapid precipitation of CaCO 3 at surface springs. This environment enables us to apply the recently developed Δ 47 -Δ 48 dual clumped isotope analysis to probe kinetic isotope effects (KIEs) and timescales of CO 2 transformation in a region with the potential for geological CO 2 sequestration. Here, we analyzed CaCO 3 recovered from various localities and identified significant kinetic fractionations associated with CO 2 absorption in a majority of samples, characterized by enrichment in Δ 47 values and depletion in Δ 48 values relative to equilibrium. Surface floes exhibited the largest KIEs (ΔΔ 47 : 0.163‰, ΔΔ 48 : -0.761‰). Surface floe samples begin to precipitate out of solution within the first hour of CO 2 absorption, and the dissolved inorganic carbon (DIC) pool requires a residence time of >100 h to achieve isotopic equilibria. The Δ 48 /Δ 47 slope of samples from the Cedars (-3.223 ± 0.519) is within the range of published theoretical values designed to constrain CO 2 hydrolysis-related kinetic fractionation (-1.724 to -8.330). The Δ 47 /δ 18 O slope (-0.009 ± 0.001) and Δ 47 /δ 13 C slope (-0.009 ± 0.001) are roughly consistent with literature values reported from a peridotite in Oman of -0.006 ± 0.002 and -0.005 ± 0.002, respectively. The consistency of slopes in the multi-isotope space suggests the Δ 47 -Δ 48 dual carbonate clumped isotope framework can be applied to study CO 2 -absorption processes in applied systems, including sites of interest for geological sequestration.

58 GEOSCIENCES↗
DOE OSTI2026

Hydrogen release from aqueous based LOHCs: The role of water in the potassium formate/bicarbonate cycle

Here, in this study, we evaluate the chemical and physical properties governing the potential round-trip efficiency of an aqueous formate/bicarbonate cycle for hydrogen storage, with a particular focus on potassium formate as a potential liquid hydrogen carrier. Using thermodynamic parameters, we predict the conversion of formate to bicarbonate across varying temperatures and pressures, revealing that hydrogen release is relatively insensitive to temperature but highly dependent on pressure. Our findings indicate that hydrogen uptake is highly efficient, whereas hydrogen release poses a greater challenge, necessitating detailed optimization to enhance round-trip efficiency. We calculate the solubility limits of bicarbonate salts influenced by the common ion effect to enable the prediction of target conversion ranges that prevent the precipitation of potassium bicarbonate in a reactor. Furthermore, the energy efficiency of hydrogen release was assessed based on the heating requirements of the aqueous solutions. This analysis maximizes the round-trip efficiency by balancing the solubility limits of bicarbonate, the heat capacity of aqueous formate solutions and conversion based on thermodynamic equilibria. Taking these factors into consideration we suggest reaction conditions that could be utilized in a systems analysis to calculate the levelized cost of storage using the bicarbonate/formate cycle at commercial scales.

common ion effect↗
DOE OSTI2020

Critical evaluation and thermodynamic modeling of the Pd–Sn system

The single crystalline material PdSn 4 , a homologue of the Dirac nodal arc semimetal PtSn 4 , is a promising candidate to search for new topological states with fascinating quantum physical properties. In the present work, the thermochemistry and phase diagram information of the Pd–Sn system published in the literature was collected and critically reviewed, two thermodynamic reassessments of the Pd–Sn system were carried out in the frame of the CALPHAD approach to optimize the Gibbs free energy of each phase presented in the system. The liquid phase was described using respectively the Bragg-Williams random mixing model and the associate solution model. All the reliable thermodynamic property and phase equilibria data can be generally described using the presently obtained thermodynamic descriptions of the Pd–Sn system. Issues related to the further improvement of the thermodynamic descriptions of the Pd–Sn system are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2025

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. In conclusion, we present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Unmasking the Conformational Stability and Inhibitor Binding to SARS-CoV-2 Main Protease Active Site Mutants and Miniprecursor

We recently demonstrated that inhibitor binding reorganizes the oxyanion loop of a monomeric catalytic domain of SARS CoV-2 main protease (MPro) from an unwound (E) to a wound (active, E*) conformation, independent of dimerization. Here we assess the effect of the flanking N-terminal residues, to imitate the MPro precursor prior to its autoprocessing, on conformational equilibria rendering stability and inhibitor binding. Thermal denaturation (T m ) of C145A mutant, unlike H41A, increases by 6.8 °C, relative to wild-type mature dimer. An inactivating H41A mutation to maintain a miniprecursor containing TSAVL[Q or E] of the flanking nsp4 sequence in an intact form [ (-6) MPro H41A and (-6*) MPro H41A , respectively], and its corresponding mature MPro H41A were systematically examined. While the H41A mutation exerts negligible effect on T m and dimer dissociation constant (K dimer ) of MPro H41A , relative to the wild type MPro, both miniprecursors show a 4–5 °C decrease in T m and > 85-fold increase in K dimer as compared to MPro H41A . The K d for the binding of the covalent inhibitor GC373 to (-6*) MPro H41A increases ~12-fold, relative to MPro H41A , concomitant with its dimerization. While the inhibitor-free dimer exhibits a state in transit from E to E* with a conformational asymmetry of the protomers’ oxyanion loops and helical domains, inhibitor binding restores the asymmetry to mature-like oxyanion loop conformations (E*) but not of the helical domains. Finally, disorder of the terminal residues 1–2 and 302–306 observed in both structures suggest that N-terminal autoprocessing is tightly coupled to the E-E* equilibrium and stable dimer formation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Effect of alloying with Sc, Nb and Zr on reduction-diffusion synthesis of magnetically hard Sm(Fe,Co,Ti) 12 -based monocrystalline powders

Powders of Sm(Fe,Co) 11.2 Ti 0.8 alloys modified with Sc, Nb and Zr, as well as with additional Ti were prepared by reducing mechanically activated raw oxides with Ca metal in the furnace preheated to 990–1250 °C. Expansion of the crystal lattice upon introduction of Nb or additional Ti implies that atoms of these elements replace the smaller Fe atoms in the tetragonal ThMn 12 -type structure. On the other hand, contraction of the lattice upon introduction of Sc or Zr was smaller than what was expected for replacement of the Sm atoms, which suggests that the Sc and Zr atoms replace both the Sm and Fe atoms. Washing away the reduction byproducts expands the crystal lattice of the 1:12 particles and increases their coercivity. The lattice expansion associated with the washing is believed to be caused by interstitial H atoms; more research, however, is needed to establish the mechanism(s) of the washing effect on the coercivity. The earlier reported development of a high coercivity in zirconium-modified monocrystalline particles achieved by increasing the reduction annealing temperature to ≈1200 °C was similarly characteristic of the particles modified with Sc (the coercivity reaches 11.5 kOe) and Nb (8.1 kOe), but not for the particles prepared with additional Ti where the maximum coercivity of 8.3 kOe develops for a lower annealing temperature. Furthermore, it is concluded that Sc, Nb and Zr modify the high-temperature phase equilibria of the Sm(Fe,Co) 11.2 Ti 0.8 alloys allowing for an effective high-temperature processing, whereas the alloy coercivity increases with the synthesis temperature through a different, still unknown mechanism which may involve suppression of the defects specific to the 1:12 crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of fission product inclusion on phase development in U 3 Si 2 fuel

Due to its high thermal conductivity and uranium density, U 3 S i2 has been considered as a candidate for use as an accident tolerant fuel (ATF). In order to fully assess its suitability and performance as a fuel, the impact of fission products (FPs) on the stability and performance of U 3 Si 2 must be investigated. The interactions of FPs and U 3 Si 2 have had relatively little study until now and require experimental and computational examination. U 3 Si 2 was doped with individual FPs to explore U-Si-FP interactions and phase equilibria that may impact the performance of the ATF during irradiation. Elemental Ce, Mo, Y, or Zr were used to individually dope U 3 Si 2 at a concentration of 5 wt% FP. A diffusion couple of a 1:1 Mo:Zr alloy and U 3 Si 2 was heated to 1200 °C in order to consider the impacts of multiple FPs on the stability and structure of the fuel. Samples were characterized for FP solubility and secondary phase formation using electron microscopy, energy dispersive spectroscopy, and x-ray diffraction. First principles density functional theory calculations complemented the experimental effort to understand FP behavior. Experimental and computational findings were used in the development of a thermodynamic database containing 8 major FPs and their associated silicide phases. Finally, fuel compositions generated from depletion calculations were used to thermodynamically model the equilibrium phases of the fuel undergoing burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗