Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALPHAD”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Design of an alumina forming coating for Nb-base refractory alloys

Refractory multi-principal element alloys (RMPEAs) promise to significantly enhance gas turbine engine efficiency, but their poor oxidation performance inhibits their implementation. Alumina-forming bond coat alloys can provide oxidation protection, but discovering suitable chemistries remains a challenge. We employed a design methodology that screens for alumina-formation capability using Al activity and phase constitution predictions from CalPhaD (Thermo-Calc). Alloy down-selection from approximately 7,800 alloys in the Nb-Si-Ti-Al-Hf system was conducted via analysis of calculated thermodynamic properties with number-density topology style maps. This approach is validated by creating and testing the composition Nb 12 Si 23 Ti 24 Al 36 Hf 5 , which forms protective alumina scales up to 1400 °C and resists pesting at 800 °C. Further, the alloy has an average coefficient of thermal expansion of ~10.1 ppm/K, making it well matched to Nb-based refractory alloys. The methodology will be useful for the design of coatings for RMPEAs, enabling their implementation and significant efficiency benefits sooner.

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 OSTI2026

Controlling homogenization length scales and microstructure in additively manufactured Ti-Ta functionally graded materials

Materials with smooth compositional gradients or functionally grade materials (FGMs) produced via additive manufacturing (AM), enables joining dissimilar materials and optimizing multiple properties in advanced engineering applications. However, as-printed AM microstructures exhibit micro-segregation and solidification defects which, when combined with controlling macroscale gradient properties, complicates necessary post-processing. Here, we use CALPHAD-informed diffusion modelling to design post-processing heat treatments for lightweight to refractory FGMs. Ti-Ta (0 to 85 at. % Ta) FGMs were fabricated using laser-based directed energy deposition AM. Post-processing heat treatments at 1000° C and 1500° C were designed to promote homogenization across specific length scales and experimentally validated. Investigation of chemical segregation and microstructures demonstrated that the length scale of homogenization is controlled as a function of time, temperature, and local composition. Ta-rich regions exhibited incomplete homogenization compared to Ti-rich layers. Unmelted Ta particles were found to completely dissolve at 1500 °C. By controlling cooling rate (200 °C/min), martensitic structures were produced between 14–36 at. % Ta, consistent with martensite-start temperatures calculations, while furnace cooling (2 °C/min) produced α+β morphologies. This work establishes a validated predictive framework for designing post-processing to tailor microstructure and chemical architecture in AM FGMs, facilitating their deployment in demanding environments.

Materials science↗
DOE OSTI2026

A robust alloy design (RAD) strategy for next-generation (IV) nuclear fission reactors

Next-generation nuclear reactors demand structural materials capable of withstanding extreme conditions, including high temperatures, intense neutron flux, and corrosive environments. Multi-Principal Element Alloys (MPEAs) have emerged as promising candidates due to their exceptional radiation tolerance, thermal stability, and compositional flexibility. This study introduces a versatile and customizable Robust Alloy Design (RAD) strategy for systematically designing MPEAs for GEN-IV reactor fuel cladding. The RAD framework integrates nuclear-relevant selection criteria, empirical parameter assessments, and high-throughput CALPHAD simulations to efficiently narrow compositional space and identify stable alloys. A unified RAD score developed for the first time, combines key performance metrics, including fuel-clad chemical interaction (FCCI), neutron absorption cross-section (NAC), valence electron configuration (VEC), and melting point factor (MPF), into a flexible ranking system adaptable to reactor-specific priorities. Among 724 candidates, V555(5Al–5Cr–5Fe–85V) emerged as the top alloy, validated experimentally with a homogeneous single-phase BCC microstructure and superior mechanical properties (nano-indentation: 3.389 ± 0.258 GPa; Vickers hardness: 240 ± 6.7 HV), significantly outperforming Zircaloy-4 and V-4Cr-4Ti. Importantly, the RAD strategy is not limited to nuclear applications; its customizable weighting system enables scalability to other extreme environments. This adaptability positions RAD strategy as a versatile tool for advanced materials design across multiple industries.

Alloy design↗
DOE OSTI2026

Response to the comments on creep-induced elemental redistribution at grain boundaries of 304L stainless steel – An experimental evidence for diffusional creep mechanisms

Sandström and He reanalyzed our 304L creep data (Acta Materialia 121137, 2025) at 10 MPa/750 °C and proposed that: (i) primary dislocation creep dominates based on the ϕ -model; (ii) measured creep rates deviate from published 304H data and Coble creep predictions using literature grain boundary (GB) diffusivity; and (iii) creep-induced elemental redistribution (CIER) could be explained by dislocation emission/absorption (DEAA) and GB creep (GBC), as opposed to diffusional creep. Here, in this response, we propose that: (i) the creep behavior at 15 MPa/750 °C—for which we demonstrated CIER—is only partially captured by the ϕ -model; (ii) direct comparison between 304L and 304H creep data is inappropriate due to carbides in 304H, and Coble creep predictions estimated from CALPHAD-based bulk diffusivity span experimental rates; and (iii) CIER exhibits distinct fingerprint for diffusional creep characterized by anisotropy and stress orientation dependence whereas neither DEAA nor GBC have been mechanistically substantiated.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Determination of Site Occupancy in the M–Pd–Zn (M = Cu, Ag, and Au) γ-Brass Phase by CALculation of PHAse Diagrams Modeling and Rietveld Refinement

The Pd–Zn γ-brass phase provides exciting opportunities for synthesizing site-isolated catalysts with precisely controlled Pd active site ensembles. Introducing a third metallic element into the γ-brass lattice further perturbs the catalytic active site ensembles. Here, in this work, we introduce coinage metallic elements M (M = Cu, Ag, and Au) into the Pd–Zn γ-brass phase and investigate the site occupation factors of each element in the γ-brass lattice. The CALculation of PHAse Diagrams (CALPHAD) modeling approach supported by energetics predicted by the density functional theory and X-ray and neutron diffraction with Rietveld refinement were used to identify the SOF on each Wyckoff site for various M amounts alloyed into the Pd–Zn γ-brass phase. The present analysis unveils the strong preference for Pd occupying the outer tetrahedral (OT) site in the γ-brass lattice, while the coinage metallic elements tend to substitute for Zn on the octahedral (OH) site. The determination of site occupancy in the bulk M–Pd–Zn γ-brass phase provides opportunities to investigate and tailor potential catalytically active site ensembles in the γ-brass phase materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First principles free energy model with dynamic magnetism for δ -plutonium

We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.

dynamic↗
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 OSTI2026

ML-AMD/exa-pd

Exa-pd is a highly parallelizable workflow for constructing multi-element phase diagrams (PDs). It uses standard sampling techniques—molecular dynamics (MD) and Monte Carlo (MC)—as implemented in the LAMMPS package, to simultaneously sample multiple phases on a fine temperature–composition mesh for free-energy calculations. The workflow uses Parsl as a global controller to manage the MD/MC jobs to achieve massive parallelization with almost ideal scalability. The resulting free energies of both liquid and solid phases (including solid solutions) are then fed to CALPHAD modeling using the PYCALPHAD package for the construction of a multi-element PD.

Zhang, Feng [Ames Laboratory (AMES), Ames, IA (Uni↗
DOE OSTI2026

Supporting Data for "Pitfalls in Parameters: Practical Process Development in Chemical Vapor Processing of SiC"

CALculation of PHAse Diagrams (CALPHAD) was used to produce equilibrium phase diagrams of the MTS/H2 system. These were calculated in FactSage 8.2 using a combination of the FactPS and SGTE Solutions 2022 databases. Equilibrium was calculated for 185 gas species, 36 pure liquids, 8 pure solids, 1 liquid solution, and 22 solid solutions. Phase diagrams were calculated with 25 °C steps in temperature from 700 to 1325 °C and 10 torr steps in pressure from 20 to 760 torr.

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

High-Entropy Alloy R&D for Accelerator Beam Window Applications

High-Entropy Alloys are a class of novel material that can offer improved resistance to beam-induced radiation damage and thermal shock. Development of these new alloys to serve as beam windows in multi-megawatt accelerator target applications is ongoing at Fermilab. Currently we are investigating AlCoCrMnTiV alloy systems of 4-6 components for service as beam windows; these compositions are predicted by CALPHAD simulation to have a low density and single-phase BCC crystal structure. The microstructures of these systems are being studied by electron microscopy techniques such as energy dispersive X-ray spectroscopy (EDS) to determine elemental homogeneity and composition, electron backscatter diffraction (EBSD) to quantify grain structure and orientation, and transmission electron microscopy (TEM) to observe defect structures and precipitate formation; nanoindentation is used to probe microstructural mechanical properties. Initial bulk property characterization utilizes differential scanning calorimetry to quantify specific heat capacity (cp), dilatometry to determine the coefficient of thermal expansion (CTE), and time-domain thermoreflectance (TDTR) to measure thermal conductivity (K). A miniature tensile testing apparatus is also being developed to test tensile properties. Post-irradiation examination is currently ongoing for several of these compositions that have been irradiated by low-energy ions to damage levels and at a temperature relevant to beam window applications at future next-generation accelerator facilities. This talk will briefly describe the alloy design and synthesis before going into more depth covering microstructural pre-characterization, and post-irradiation examination results from low-energy ion irradiated specimens. This will be followed by the plans for alloy down selection and future prototypic proton irradiations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Dissimilar metal joining of Ti-6Al-4V and Inconel 625

The dissimilar joining of Ti-6Al-4V and Inconel 625 by using a suitable interlayer shows potential design flexibility by avoiding the brittle intermetallic phases from direct welding. These phases include NiTi , NiTi 2 , and Ni 3 Ti . The objective of this work is to facilitate laser welding of Ti-6Al-4V and Inconel 625 through a vanadium interlayer to form an intermetallic free, strong, and corrosion resistant joint. Vanadium was chosen based on binary phase diagram solubilities, literature, and theoretical CalPhaD calculations in PanDat. Laser welding was systematically explored from based on beam mode, offset, power, and scan speed. Varied parameters include laser power from 325 to 650 watts, laser scan speed from 80 to 160 in/min, and beam offset up to a ¼ of the beam diameter. Parameters were chosen based on previous dissimilar joining research and those found in the literature. Metallography, including optical and scanning electron microscopy, was used to analyze the weld chemical composition, grain structure, and defects through both the fusion and heat affected zones. The results and discussion will examine the relationship between process parameters, microstructure, and observed weld defects. Our findings will not only assess the weldability of vanadium and Inconel 625 but also establish a transferable framework for studying other dissimilar material combinations, providing valuable insights for future research in this field.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Entropy Alloys for Accelerator Beam Window Applications

Development of novel high-entropy alloys (HEAs) is currently underway for potential use as beam windows in future multi-megawatt target systems at Fermilab. HEAs encompass a new class of materials with a vast design space allowing for material properties to be tailored for particular applications and to potentially offer improved resistance to beam-induced radiation damage and thermal shock effects. The alloy systems being studied consist of several compositions of AlCoCrMnTiV with 4 6 component elements. These alloys are all predicted by CALPHAD simulation to have a single-phase BCC crystal structure and low density, with some compositions displaying ordered, nanoscale precipitates. This presentation will briefly discuss alloy design and synthesis before giving a detailed description of the characterization studies of these HEAs in both the pristine state and post-irradiation by low-energy heavy ions to high damage levels. Electron microscopy techniques to quantify elemental homogeneity and composition, determine grain size, shape, and orientation, and quantify lattice parameters, defect structures and precipitate phases are all being used to study alloy microstructures. Mechanical properties of the alloys at the microscale will be reported. The evolution of these properties as a function of radiation damage will also be described. Bulk thermal characteristics of these HEAs have been tested to measure specific heat capacity and coefficient of thermal expansion as a function of temperature. To determine bulk tensile properties a miniature tensile testing apparatus is under development; it s commissioning will be covered briefly. The talk will conclude with our plans for alloy down-selection.

Burleigh, A. [Fermilab]↗
DOE OSTI2026

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↗
NASA NTRS2000

Thermodynamic Analysis of the Combustion of Metallic Materials

Two types of computer codes are available to assist in the thermodynamic analysis of metallic materials combustion. One type of code calculates phase equilibrium data and is represented by CALPHAD. The other type of code calculates chemical reaction by the Gordon-McBride code. The first has seen significant application for alloy-phase diagrams, but only recently has it been considered for oxidation systems. The Gordon-McBride code has been applied to the combustion of metallic materials. Both codes are limited by their treatment of non-ideal solutions and the fact they are limited to treating volatile and gaseous species as ideal. This paper examines the significance of these limitations for combustion of metallic materials. In addition, the applicability of linear-free energy relationships for solid-phase oxidation and their possible extension to liquid-phase systems is examined.

Wilson, D. Bruce↗
NASA NTRS2001

Calculation of Phase Equilibria in the Y2O3-Yb2O3-ZrO2 System

Rare earth oxide stabilized zirconias find a wide range of applications. An understanding of phase equilibria is essential to all applications. In this study, the available phase boundary data and thermodynamic data is collected and assessed. Calphad-type databases are developed to completely describe the Y2O3-ZrO2, Yb2O3-ZrO2, and Y2O3-Yb2O3 systems. The oxide units are treated as components and regular and subregular solution models are used. The resultant calculated phase diagrams show good agreement with the experimental data. Then the binaries are combined to form the database for the Y2O3-Yb2O3-ZrO2 psuedo-ternary.

Jacobson, Nathan S.↗
NASA NTRS2001

Thermodynamic Activity Measurements with Knudsen Cell Mass Spectrometry

Coupling the Knudsen effusion method with mass spectrometry has proven to be one of the most useful experimental techniques for studying the equilibrium between condensed phases and complex vapors. The Knudsen effusion method involves placing a condensed sample in a Knudsen cell, a small "enclosure", that is uniformly heated and held until equilibrium is attained between the condensed and vapor phases. The vapor is continuously sampled by effusion through a small orifice in the cell. A molecular beam is formed from the effusing vapor and directed into a mass spectrometer for identification and pressure measurement of the species in the vapor phase. Knudsen cell mass spectrometry (KCMS) has been used for nearly fifty years now and continues to be a leading technique for obtaining thermodynamic data. Indeed, much of the well-established vapor specie data in the JANAF tables has been obtained from this technique. This is due to the extreme versatility of the technique. All classes of materials can be studied and all constituents of the vapor phase can be measured over a wide range of pressures (approximately 10(exp -4) to 10(exp -11) bar) and temperatures (500-2800 K). The ability to selectively measure different vapor species makes KCMS a very powerful tool for the measurement of component activities in metallic and ceramic solutions. Today several groups are applying KCMS to measure thermodynamic functions in multicomponent metallic and ceramic systems. Thermodynamic functions, especially component activities, are extremely important in the development of CALPHAD (Calculation of Phase Diagrams) type thermodynamic descriptions. These descriptions, in turn, are useful for modeling materials processing and predicting reactions such as oxide formation and fiber/matrix interactions. The leading experimental methods for measuring activities are the Galvanic cell or electro-motive force (EMF) technique and the KCMS technique. Each has specific advantages, depending on material and conditions. The EMF technique is suitable for lower temperature measurements, provided a suitable cell can be constructed. KCMS is useful for higher temperature measurements in a system with volatile components. In this paper, we briefly review the KCMS technique and identify the major experimental issues that must be addressed for precise measurements. These issues include temperature measurements, cell material and cell design and absolute pressure calibration. The resolution of these issues are discussed together with some recent examples of measured thermodynamic data.

Copland, Evan H.↗
NASA NTRS2006

Compositional Pathways and Capillary Effects during Early-stage Isothermal Precipitation in a Nondilute Ni-Al-Cr Alloy

For a Ni-5.2 Al-14.2 Cr at.% alloy with moderate solute supersaturations, the compositional pathways, as measured with atom-probe tomography, during early to later stage y'(LI2)-precipitation (R = 0.45-10 nm), aged at 873 K, are discussed in light of a multi-component coarsening model. Employing nondilute thermodynamics, detailed model analyses during quasistationary coarsening of the experimental data establish that the y/y' interfacial free-energy is 22- 23+/-7 mJ/sq m. Additionally, solute diffusivities are significantly slower than model estimates. Strong quantitative evidence indicates that an observed y'-supersaturation of Al results from the Gibbs-Thomson effect, providing the first experimental verification of this phenomenon. The Gibbs-Thomson relationship, for a ternary system, as well as differences in measured phase equilibria with CALPHAD assessments, are considered in great detail.

Sudbrack, Chantal K.↗