Engineering Papers⌕ Search

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Oxidation Lifetime Modeling of 625 and 120 Foils After Long-term Exposure in Flowing Air + 10% ${\hbox {H}_2\hbox {O}}$ at 700 and 800 $^\circ$C

A model describing oxidation kinetics considering simultaneous oxidation and volatilization of external oxide scales was integrated with a CALPHAD-based coupled thermodynamic and kinetic approach to predict the oxidation-induced lifetime of foils of alloys 625 and 120. Long-term exposures in flowing air +10% H 2 O at 700 and 800 °C under two flow velocities were used to validate the modeling results. The time to a critical Cr concentration of 10 wt% at the oxide/alloy interface, t 10 , was shown to correlate with the breakdown of the Cr 2 O 3 scale for the studied experimental conditions. Finally, the modeling approach was extended to consider time-dependent oxidation and volatilization rate constants to develop a more appropriate lifetime criterion for microturbine recuperator foils.

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

Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles

The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems

Molten salt electrolysis can recover the actinides from spent nuclear fuels, and it involves a eutectic LiCl-KCl in molten form as an electrolyte. During reprocessing, the concentration of fission products such as La, Nd, Pr, etc., increases and affects the recovery efficiency of the electrolyte. In this work, thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems was carried out using the CALPHAD (Calculation of Phase Diagrams) approach for the first time. The thermodynamic functions for the pure salts were taken from the SGTE (Scientific Group Thermodata Europe) Substances (SSUB) database. The experimental thermochemical and phase equilibria data available in the literature were used as input for the assessment of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems. The model parameters for the KCl-LiCl system were adjusted to include the new Gibbs energy descriptions for the pure salts. In addition, the sublattice model for the liquid phase in the LiCl-PrCl 3 system was modified and reassessed to ensure the model compatibility for higher-order extrapolation. There is a good agreement between the experimental and calculated thermochemical and phase diagram data for all the optimized constituent binaries and the ternary system. Furthermore, this work will be beneficial for determining the solubility limit of PrCl 3 in molten LiCl-KCl electrolytes and their thermodynamic properties for improving the efficiency of the pyrochemical process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗
DOE OSTI2026

An Integrated Computational Materials Engineering (ICME) Approach to Design Nonlinear Transition Zones Between Dissimilar Metals

Current approaches to designing graded transition joints (GTJs) between dissimilar metals often rely on linear changes in both composition profiles and thickness of each sublayer. This increases fabrication cost and may not be optimal with respect to residual stress or the formation of undesirable phases. Here, in this study, GTJs between P91 ferritic/martensitic steel and 347H austenitic stainless steel were designed using Integrated Computational Materials Engineering (ICME) principles with nonlinear composition and length profiles. Guided by inputs from classical mechanics and CALPHAD predictions of carbon chemical potential, a novel transition zone consisting of five discrete compositions was proposed, with the thickness of each sublayer varying according to a brachistochrone-inspired distribution. In addition to carbon potential gradients, CALPHAD was used to predict coefficients of thermal expansion, which were incorporated into finite element models to evaluate stress evolution. The proposed nonlinear design resulted in a smoother carbon potential gradient, lower carbon depletion at the P91 interface, and a comparable residual stress under long-term thermal exposure, compared to a conventional linear design using ten sublayers with equal thickness. This work introduces a brachistochrone-inspired distribution for GTJ design, offering a general framework for optimizing graded interfaces between dissimilar metals.

Directed Energy Deposition↗
DOE OSTI2024

Boron and lithium aqueous thermochemistry to model crud deposition in pressurized water reactors

Here, a comprehensive database of boron and lithium aqueous thermochemistry has been developed for use at elevated temperatures. The Helgeson-Kirkham-Flowers (HKF) formalism provides a framework to describe thermodynamic properties over a broad range of temperatures and pressures. Accuracy at high temperatures is vital to modeling nickel oxide and ferrite fuel deposits that occur in pressurized water nuclear reactors (PWRs). CALPHAD (CALculation of PHAse Diagrams) calculations are performed at PWR crud conditions to predict the stability regions of the solid lithium metaborate (LiBO 2 ) and lithium tetraborate (Li 2 B 4 O 7 ) precipitates. In addition, similar calculations are performed using sodium and potassium instead of lithium in order to assess the thermodynamics associated with preventing the formation of such precipitates. This database contributes to understanding of crud formation and composition and will aid in the prediction of phenomena such as Crud-Induced Power Shifts (CIPS).

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

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solidification cracking of refractory alloys: a computational and machine learning study to investigate composition-dependence for improved weldability and additive manufacturability

Large-batch numerical, CALculation of PHAse Diagrams (CALPHAD)-based solidification cracking calculations are performed and then analyzed with machine learning methods to generate models that relate chemistry of refractory alloys to cracking susceptibility. Kou’s solidification cracking index is used to study the refractory alloys including O, N, C binary mixtures with Mo, Ta, Nb, and W, the molybdenum-based TZM, Niobium-based C103, and Tantalum-based T111 and Ta-10 W, as well as hypothetical refractory ternary alloys. Findings strongly validate Kou’s Crack Susceptibility Index (CSI) against Varestraint test data for Nb- and Ta-based alloys, establishing CSI thresholds where refractory alloys with CSI < 15,000 K are likely weldable, CSI > 15,000 K are prone to cracking, and CSI > 25,000 K are likely unweldable (or unprintable). Furthermore, interstitial elements C, N, and O significantly increase crack susceptibility, with some existing material specifications coinciding with peak cracking susceptibility concentrations. Finally, machine learning-derived elemental potency factors enable rapid prediction of CSI from alloy chemistry for C103, TZM, Ta-10 W, and T-111 alloys. These results provide practical guidance for feedstock selection, powder reuse limits, and alloy specification amendments for welding and additive manufacturing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental investigation and thermodynamic modeling of the binary Pt–Sn system

Here, the thermodynamics of the Pt–Sn system is revisited with attention to phase equilibria involving PtSn 4 , a topological quantum material with interesting physical properties. The phase relationship in the Sn-rich region of the Pt–Sn binary system was investigated by a combination of differential scanning calorimetry and equilibration experiments. Based on the present experiments and a critical evaluation of the available thermodynamic and phase diagram data reported in the literature, thermodynamic assessments of the entire Pt–Sn system were performed using the CALPHAD approach. Two sets of thermodynamic parameters were optimized, where the liquid phase was described using the Bragg-Williams random mixing model and the associate solution model, respectively. The modeling results are compared with the reliable experimental data and are shown to provide good representations of the phase equilibria and thermodynamic properties in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Primary solidification of ternary compounds in Al-rich Al–Ce–Mn alloys

Primary solidification of ternary compounds Al 20 Mn 2 Ce and Al 10 Mn 2 Ce were analyzed through the coupling of the thermodynamic modeling and classic nucleation theory. Thermodynamic models of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce were developed using the CALPHAD approach based on first-principles calculated enthalpy of formation and experimental data obtained from this work and the literature. The analysis suggested that despite the larger thermodynamic driving force for nucleation of Al 10 Mn 2 Ce, nucleation is preferred for the Al 20 Mn 2 Ce phase in the highly undercooled liquid due to its smaller interfacial energy. Therefore, manufacturing methods with rapid cooling rates will favor primary solidification of Al 20 Mn 2 Ce for Al-rich Al–Ce–Mn alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermodynamic re-modelling of the Cu–Nb–Sn system: Integrating the nausite phase

Currently available Cu–Nb–Sn phase diagrams lack the recently discovered nausite phase (Cu,Nb)Sn 2 , which is an important intermediate in the course of thermal processing of superconducting Nb 3 Sn wires. Processing decisively determines the resulting microstructure of Nb 3 Sn and, thus, its superconducting properties. Lack of suitable and complete phase diagrams, however, obstructs rational design of such thermal processing procedures. To close this gap and to obtain valid knowledge of homogeneity and stability range of nausite, various Cu–Nb–Sn samples, which are heat-treated between 300 °C and 500 °C, are investigated. By means of energy-dispersive X-ray spectroscopy (EDX), a temperature-dependent homogeneity range of nausite is observed, which covers average mole fractions of Cu between 0.09 and 0.15. This is correlated with a change in the mean atomic volume and can be seen in the lattice parameters determined by X-ray diffraction (XRD). Additionally performed first-principles calculations on different CuSn 2 and NbSn 2 model structures confirm this trend. Furthermore, the peritectic decomposition of nausite to NbSn 2 and liquid at 586 °C is determined by means of in situ XRD and differential scanning calorimetry (DSC). By using the CALPHAD (CALculation of PHase Diagrams) approach, all these findings are used to extend a previous thermodynamic description of the Cu–Nb–Sn system by including the nausite as an additional phase. Finally, with this noteworthy integration, the updated modelling of the Cu–Nb–Sn system can be used for optimizing the multistage heat-treatment steps during processing superconducting Nb 3 Sn wires.

36 MATERIALS SCIENCE↗
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 OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal Property Modeling and Assessment of the Physical Properties of FLiNaK

Here, the thermodynamic properties of the LiF–NaF–KF system and its subsystems were re-evaluated using the CALPHAD method, focusing on the deviations from ideality in the heat capacity [C p (T)] of mixtures, temperature-dependence of mixing enthalpies (Δ mix H), vapor pressures of mixtures, and enthalpies of fusion (ΔH fus ) using reported values and measurements. The results of this work confirmed a pseudoternary eutectic at 732.9 K and 46.5LiF–11.5NaF–42KF mol % (FLiNaK) with determined Cp(T) = 69.822 + 0.000679T + 2,137,152T –2 J mol –1 K –1 (732.9 < T < 1200 K) and ΔH fus = 18,051 J mol –1 in good agreement to the experimental values. The resulting accurate thermodynamic representation can be used to find relevant excess thermochemical properties in the whole pseudoternary space and provide insights into their composition and temperature dependence. A detailed evaluation of select thermophysical properties was also conducted to determine recommended density, thermal diffusivity, thermal conductivity, vapor pressures, and boiling point with calculated uncertainties over the temperature range of interest for FLiNaK–salt energy applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamic re-modeling of the Yb-Sb system aided by first-principles calculations

Here, the thermodynamic description of the Yb-Sb binary system is developed by means of the CALculations of PHAse Diagrams (CALPHAD) method by combining experimental data in the literature and predictions from first-principles calculations based on density functional theory (DFT) in the literature and the present work. Two pseudopotentials of Yb are compared in the present DFT-based calculations with 14 and 13 f-electrons frozen in the core, i.e., 5p 6 6s 2 and 5p 6 6s 2 5 d 1 electrons as valence electrons, termed Yb_2 and Yb_3, respectively. It is shown that the phonon spectrum of the YbSb phase calculated using the Yb_3 pseudopotential does not have imaginary phonon modes and is subsequently used to predict its temperature dependent thermodynamic properties by the DFT-based quasiharmonic phonon calculations. The present thermodynamic database includes the Yb 16 Sb 11 phase in addition to five intermetallic phases that were considered in previous modeling studies, i.e., YbSb 2 , YbSb, Yb 11 Sb 10 , Yb 4 Sb 3 , and Yb 5 Sb 3 . The high temperature orthorhombic structure of the Yb 5 Sb 3 phase is not considered in the present work as it was stabilized by hydrogen. The associate solution model is used to describe the short-range ordering behavior in the liquid phase. The calculations from the present thermodynamic model show good agreement with thermochemical and phase equilibrium data from both the present work and the literature.

36 MATERIALS SCIENCE↗