A thermodynamic approach to modeling multicomponent FeTi-based alloys for hydrogen storage
Not Available
08 HYDROGEN↗
SEARCH · Engineering Papers
Results for “CALPHAD modeling”
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.
Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds
The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors
Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Aging heat treatment design for Haynes 282 made by wire-feed additive manufacturing using high-throughput experiments and interpretable machine learning
Wire-feed additive manufacturing (WFAM) produces superalloys with complex thermal cycles and unique microstructures, often requiring optimized heat treatments. To address this challenge, we present a hybrid approach that combines high-throughput experiments, precipitation simulation, and machine learning to design effective aging conditions for the WFAM Haynes 282 superalloy. Our results demonstrate that the γ’ radius is the critical microstructural feature for strengthening Haynes 282 during post-heat treatment compared with the matrix composition and γ’ volume fraction. New aging conditions at 770°C for 50 hours and 730°C for 200 hours were discovered based on the machine learning model and were applied to enhance yield strength, bringing it on par with the wrought counterpart. This approach has significant implications for future AM alloy production, enabling more efficient and effective heat treatment design to achieve desired properties.
CALPHAD↗
Advanced data science toolkit for non-data scientists – A user guide
Emerging modern data analytics attracts much attention in materials research and shows great potential for enabling data-driven design. Data populated from the high-throughput CALPHAD approach enables researchers to better understand underlying mechanisms and to facilitate novel hypotheses generation, but the increasing volume of data makes the analysis extremely challenging. Here in this paper, we introduce an easy-to-use, versatile, and open-source data analytics frontend, ASCENDS (Advanced data SCiENce toolkit for Non-Data Scientists), designed with the intent of accelerating data-driven materials research and development. The toolkit is also of value beyond materials science as it can analyze the correlation between input features and target values, train machine learning models, and make predictions from the trained surrogate models of any scientific dataset. Various algorithms implemented in ASCENDS allow users performing quantified correlation analyses and supervised machine learning to explore any datasets of interest without extensive computing and data science background. The detailed usage of ASCENDS is introduced with an example of experimental high-temperature alloy data.
96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
High-throughput design of high-performance lightweight high-entropy alloys
Developing affordable and light high-temperature materials alternative to Ni-base superalloys has significantly increased the efforts in designing advanced ferritic superalloys. However, currently developed ferritic superalloys still exhibit low high-temperature strengths, which limits their usage. Here we use a CALPHAD-based high-throughput computational method to design light, strong, and low-cost high-entropy alloys for elevated-temperature applications. Through the high-throughput screening, precipitation-strengthened lightweight high-entropy alloys are discovered from thousands of initial compositions, which exhibit enhanced strengths compared to other counterparts at room and elevated temperatures. The experimental and theoretical understanding of both successful and failed cases in their strengthening mechanisms and order-disorder transitions further improves the accuracy of the thermodynamic database of the discovered alloy system. This study shows that integrating high-throughput screening, multiscale modeling, and experimental validation proves to be efficient and useful in accelerating the discovery of advanced precipitation-strengthened structural materials tuned by the high-entropy alloy concept.
36 MATERIALS SCIENCE↗
Ab Initio Modeling on The Thermodynamic and Temperature-Dependent Elastic Properties of Subsystems of The FCC FeNiCoCr Medium Entropy Alloys (MEAs)
The stability, phonon spectra, thermodynamic, and temperature-dependent elasticity of subsystems of the FCC FeNiCoCr MEAs are systematically studied by the ab initio approach. Especially, the quasi-harmonic approximation (QHA) and the innovative Zentropy theory were utilized to predict the thermodynamic properties and elastic properties of FeNi, NiCo, FeNiCo, and FeNiCoCr MEAs with the consideration of magnetic transition. With the ensemble of the partition function based on the multiplicity of each magnetic configuration, the current work successfully reproduced the Curie temperature and the Schottky anomaly of heat capacity of these four MEAs purely based on the ab initio input. Meanwhile, the elastic properties of these alloys at finite temperatures are also successfully predicted with the consideration of magnetic transition. The overall results agree well with the available experimental data and CALPHAD prediction.
36 MATERIALS SCIENCE↗
Solid-State Transformation of an Additive Manufactured Inconel 625 Alloy at 700 °C
Inconel 625, a nickel-based superalloy, has drawn much attention in the emerging field of additive manufacturing (AM) because of its excellent weldability and resistance to hot cracking. The extreme processing condition of AM often introduces enormous residual stress (hundreds of MPa to GPa) in the as-fabricated parts, which requires stress-relief heat treatment to remove or reduce the internal stresses. Typical residual stress heat treatment for AM Inconel 625, conducted at 800 °C or 870 °C, introduces a substantial precipitation of the δ phase, a deleterious intermetallic phase. In this work, we used synchrotron-based in situ scattering and diffraction methods and ex situ electron microscopy to investigate the solid-state transformation of an AM Inconel 625 at 700 °C. Our results show that while the δ phase still precipitates from the matrix at this temperature, its precipitation rate and size at a given time are both smaller when compared with their counterparts during typical heat treatment temperatures of 800 °C and 870 °C. A comparison with thermodynamic modeling predictions elucidates these experimental findings. Our work provides the rigorous microstructural kinetics data required to explore the feasibility of a promising lower-temperature stress-relief heat treatment for AM Inconel 625. The combined methodology is readily extendable to investigate the solid-state transformation of other AM alloys.
36 MATERIALS SCIENCE↗
Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system
The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.
CALPHAD↗
Acceleration of Thermochemistry Solves in MOOSE and Pronghorn
This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.
36 - MATERIALS SCIENCE↗
Crystallographic stability of the BCC phase during quenching of metastable beta titanium alloy Ti-5553 and comparison of structural criteria for the predictive capability of α” martensite
This work evaluates the compositional standard currently developed for Ti-5553 powder and explores the metastable β region for sensitivity to martensitic formation during rapid quenching from the melt. Ti-5553 is a beta-stabilized titanium alloy that is increasingly being used in additive manufacturing applications. A series of alloys within the Ti-5553 compositional space was processed using two-piston spat quenching to perform rapid solidification and quenching for each of these alloy compositions. Typical empirical models such as molybdenum equivalency are not able to fully separate the retained BCC and martensitic compositions. Further, the use of thermodynamic data to estimate the transformation energy needed to form martensite can differentiate the alloys and provide a metric to further develop compositional limits for metastable beta titanium alloy development.
36 MATERIALS SCIENCE↗
Data Mining and Visualization of High-Dimensional ICME Data for Additive Manufacturing
Integrated computational materials engineering (ICME) methods combining CALPHAD with process-based simulations can produce rich, high-dimensional data for alloy and process design. In ICME methods for metallurgical applications, the visualization and interpretation of such high-dimensional data has previously been through heat maps represented in 2 or 3 dimensions. While such an approach is ideal when one variable is varied at a time, in the case of high-dimensional data with multiple variables varied simultaneously, as is the case in additive manufacturing, interpreting the trends through two- or three-dimensional heat maps becomes challenging. Here, we propose a strategy of mixed visual data mining and quantitative analysis for high-dimensional metallurgical and process data using high-throughput thermodynamic calculations. Two case studies show the application of the proposed approach. The first case study investigated the effects of feedstock chemistry on the δ ferrite formation in 316L stainless steel powders used for binder jet additive manufacturing. The second case study linked Scheil–Gulliver calculations to a process model for dissimilar joining of aluminum alloys 5356 and 6111 during laser hot-wire additive manufacturing. Both cases contained thousands of calculated data points, showcasing the utility of visual data analysis through parallel coordinate plotting, Pearson correlation coefficient matrices, and scatter matrices compared to traditional process maps. These visualization techniques can be extended to many additive manufacturing problems to capture process–structure–property relationships for additively manufactured components.
36 MATERIALS SCIENCE↗
Understanding the effect of refractory metal chemistry on the stacking fault energy and mechanical property of Cantor-based multi-principal element alloys
Multi-principal-element alloys (MPEAs) based on 3d-transition metals show remarkable mechanical properties. In this study, the stacking fault energy (SFE) in face-centered cubic (fcc) alloys is a critical property that controls underlying deformation mechanisms and mechanical response. Here, we present an exhaustive density-functional theory study on refractory- and copper-reinforced Cantor-based systems to ascertain the effects of refractory metal chemistry on SFE. We find that even a small percent change in refractory metal composition significantly changes SFEs, which correlates favorably with features like electronegativity variance, size effect, and heat of fusion. For fcc MPEAs, we also detail the changes in mechanical properties, such as bulk, Young's, and shear moduli, as well as yield strength. A Labusch-type solute-solution-strengthening model was used to evaluate the temperature-dependent yield strength, which, combined with SFE, provides a design guide for high-performance alloys. We also analyzed the electronic structures of two down-selected alloys to reveal the underlying origin of optimal SFE and strength range in refractory-reinforced fcc MPEAs. These new insights on tuning SFEs and modifying composition-structure-property correlation in refractory- and copper-reinforced MPEAs by chemical disorder, provide a chemical route to tune twinning- and transformation-induced plasticity behavior in fcc MPEAs.
36 MATERIALS SCIENCE↗
Applications of Thermochemical Modeling in Molten Salt Reactors
The extensively evaluated and consistent thermodynamic database, the Molten Salt Thermal Properties Database—Thermochemical (MSTDB-TC), was used along with additional thermodynamic values from other sources as examples of ways to examine molten salt reactor (MSR) fuel behavior. Relative stability with respect to halide potential and temperature for likely fuel and fission product components were mapped in Ellingham diagrams for the chloride and fluoride systems. The Ellingham diagrams provide a rich, visual means for identifying halide-forming components in proposed fuel/solvent salt systems. Thermochemical models and values from MSTDB-TC and ancillary sources were used in global equilibrium calculations to provide compositions for a close analysis of the behavior of a possible Molten Chloride Salt Fast Reactor and a Molten Salt Reactor Experiment-type system at high burnup (100 GWd/t). The results illustrated the oxidative nature of burnup in MSRs and provided information about redox behavior and possible control.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Determination of Location-Specific Solidification Cracking Susceptibility for a Mixed Dissimilar Alloy Processed by Wire-Arc Additive Manufacturing
Solidification cracking is a major obstacle when joining dissimilar alloys using additive manufacturing. In this work, location-specific solidification cracking susceptibility has been investigated using an integrated computational materials engineering (ICME) approach for a graded alloy formed by mixing P91 steel and Inconel 740H superalloy. An alloy mixture of 26 wt.% P91 and 74 wt.% Inconel 740H, with high configurational and total entropy, was fabricated using wire arc additive manufacturing. Microstructure characterization revealed intergranular solidification cracks in the FCC matrix, which increased in length along with the enrichment of Nb (~27 to 56 wt.%) and Cu (~87 wt.%) in the middle and top regions. DICTRA simulations to model location-specific solidification cracking susceptibility showed that the top region with the highest cooling rate (270 K/s) has the highest solidification cracking susceptibility in comparison with the middle and bottom regions. This is in good agreement with the experimentally observed varying crack length. From Scheil simulations, it was deduced that enrichment of Nb and Cu affected the solidification range as high as ~77%, in comparison with the matrix composition. The overall solidification cracking susceptibility and freezing range was highest for the 26 wt.% P91 alloy amongst the mixed compositions between P91 steel and 740H superalloy, proving that solidification characteristics play a major role in alloy design for additive manufacturing.
20 FOSSIL-FUELED POWER PLANTS↗
Computational thermodynamics and its applications
Thermodynamics is a science concerning the state of a system, whether it is stable, metastable or un- stable, when interacting with the surroundings. In this overview, fundamentals of thermodynamics are briefly reviewed through the combination of first and second laws of thermodynamics for open and nonequilibrium systems to demonstrate that the reversible equilibrium and irreversible nonequilibrium thermodynamics can be integrated to enhance the power and utilities of thermodynamics. The recent progresses in computational thermodynamics, the remaining challenges, and potential impacts in broad scientific fields are discussed here. It is shown that computational thermodynamics enables the modeling of thermodynamics of a state as a function of both external and internal variables and quantitative calculations of a broad range of properties of a multicomponent system in terms of first and second derivatives of energy, including not only equilibrium states when there are no driving forces for any internal processes and but also non-equilibrium states with driving forces for internal processes. Consequently, external constraints such as fixed strain and internal degree of freedoms such as ordering and defects can be described in a coherent framework and applied to materials design. Furthermore, two important but largely overlooked aspects in thermodynamics will be discussed, i.e. the rigorous application of statistical thermodynamics with the probability of configurations and their contributions to system properties, and the applications of second derivatives of energy with respect to either two extensive variables or two potentials or a mixture of them in terms of understanding and predicting emergent behaviors, critical phenomena, kinetic coefficients, and mechanical properties.
36 MATERIALS SCIENCE↗
Degradation science: Integrating modeling and experiments to predict localized corrosion processes (Annual Progress Report)
Additively manufactured (AM) eutectic high-entropy alloys (EHEAs), such as nano-lamellar AlCoCrFeNi 2.1 , have excellent strength, ductility, and wear resistance even at elevated temperatures, but their corrosion behavior in aggressive acids at different length scales remain poorly understood. This work investigates the corrosion behavior of laser powder bed–fused (L-PBF) AlCoCrFeNi 2.1 as a function of annealing temperatures, probing degradation mechanisms from nanoscopic to macroscopic length scales. The alloy is dual phase consisting of a ductile FCC L1 2 phase and a high-strength BCC B2 phase. Rapid solidification during L-PBF produces a far-from-equilibrium nano-lamellar structure with nearly homogeneous elemental distribution, which tends to evolve upon annealing toward Cr/Co/Fe-enriched FCC and Al/Ni-enriched B2. Three conditions were studied: as-printed, 600 °C/5 h, and 1000 °C/1 h, over which B2 lamellae coarsen, lamellar spacing increases, and elemental segregation becomes more prominent. Microstructure and chemistry were characterized by scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS), while in-situ electrochemical atomic force microscopy (EC-AFM) was used to link early (<5 h) local dissolution to microstructure after exposure in sulfuric acid. EC-AFM highlights preferential dissolution of the BCC/B2 phase where the surrounding matrix is Cr-depleted and directly quantifies the dissolution rates within individual phases, tracks the transition from early nano-scale attack to partial repassivation, to correlate height differences with current and impedance responses. To monitor longer-term behavior (up to 96 h), ex-situ AFM, SEM, and confocal imaging were combined with conventional bulk electrochemical tests, bridging nanoscale observations to micro/meso-scale damage morphologies. At the meso-scale, the deepest dissolution channels align with the build-direction lamellae and melt-pool boundaries, indicating that printing directionality guides the propagation of these localized corrosion sites. Annealing modifies corrosion by restructuring BCC/FCC phase fractions, lamellar spacing, and Cr/Al segregation, thereby changing the cathode/anode ratio and passive film stability. The results clarify how as-printed and annealed nano-lamellar architectures differ in their susceptibility to selective dissolution; how elemental segregation competes with residual stresses along the build direction. With these insights, future work will use CALPHAD-guided alloy modification to stabilize higher Cr contents in the B2 phase while retaining a dominant FCC+B2/BCC microstructure, with the goal of designing mechanically robust, corrosion-resistant EHEAs for safety-critical applications to leverage the LLNL’s broader national and global security mission.
36 MATERIALS SCIENCE↗
Thermodynamic properties of ZrCl 4 with LiCl, NaCl, KCl, CsCl, MgCl 2 , and UCl 3 for molten salt reactor applications
Here, a set of self-consistent Gibbs energy functions has been developed to describe the thermochemical behavior of the major reactive fission product ZrCl 4 with key components in chloride molten salt reactors (MSRs): LiCl–ZrCl 4 , NaCl–ZrCl 4 , KCl–ZrCl 4 , CsCl–ZrCl 4 , MgCl 2 –ZrCl 4 and UCl 3 –ZrCl 4 . Low ZrCl 4 concentration phase equilibria most relevant to MSR applications have been confirmed via differential scanning calorimetry for NaCl–ZrCl 4 , and X-ray diffraction analysis performed on equilibrated samples of NaCl–ZrCl 4 , KCl–ZrCl 4 , MgCl 2 –ZrCl 4 , and UCl 3 –ZrCl 4 . Within the framework of the modified quasi-chemical model in the quadruplet approximation, extrapolations of pseudo-binary models were generated to represent KCl–MgCl 2 –ZrCl 4 , KCl–NaCl–ZrCl 4 , and MgCl 2 –NaCl–ZrCl 4 , which show agreement with phase equilibria data where available without the use of ternary interaction parameters. The optimized thermodynamic descriptions for these systems and others are available in the open-source compendium Molten Salt Thermal Properties Database – Thermochemical (MSTDB–TC).
CALPHAD↗