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 OSTI2024

Heat affected zone liquation cracking evaluation on FeMnAl alloys

FeMnAl steels are currently generating a lot of interest with potential applications for structural parts in armored vehicles due to their lower density and outstanding mechanical properties. Despite the extensive mechanical performance and heat treatment exploration of this alloy class, further weldability investigation is required for future large-scale deployment. In the present study, the liquation cracking of four heats of cast FeMnAl alloys was investigated by the spot-Varestraint technique. The study focuses primarily on the effect of the major elements of the FeMnAl system: C, Mn and Al. Optical and electron microscopy were employed to investigate the microstructural features, and CALPHAD was employed to aid the discussion regarding the alloy’s composition differences and their effect on the liquation cracking susceptibility. The study was able to identify that compositions with the higher Mn, C, and lower Al presented the highest liquation cracking susceptibility. Conversely, composition presenting lower Mn, C, and Al showed the most resistant behavior. Furthermore, lower Al content promoted a fully-γ microstructure at low temperatures, which encouraged the appearance of longer cracks as a γ-matrix is more susceptible to HAZ cracking than a fully ferritic (α) or duplex (α + γ) microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational design of high entropy alloy coating for hydrogen turbine applications

This project aims to develop novel high entropy alloy (HEA)-based coatings to protect critical components in hydrogen-fueled turbine power system. The HEA-coatings will demonstrate superior performance in hydrogen combustion environment to commercial NiCoCrAlY coating in current natural gas turbine system. The HEA coating facilitates the formation of a protective scale of alpha-alumina to slow down the inward diffusion of oxidizing species and the outward diffusion of metal elements, and possesses ultrahigh corrosion and spallation resistance to prolong the service lifetime of critical components in hydrogen turbine power system. Aimed to accelerate the discovery of novel HEA coating compositions, high throughput computational modeling including CALPHAD and density functional theory and machine learning are performed to predict phase stability, oxygen permeability, oxidation rate constant, coefficient of thermal expansion, and mechanical properties. Based on the modeling and machine learning prediction, experimental validation is performed. Preliminary results will be presented and approaches to minimize oxidation will be discussed.

alloy design↗
DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗
DOE OSTI2024

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrodeposition of near-equiatomic CoCuFeNi multi-principal element alloys from an acidic glycine-citrate-triton X100 aqueous electrolyte

Understanding the composition and morphology control of electrodeposited CoCuFeNi is the first step to finding a general strategy for developing electrodeposition processes of unconventional alloys with large redox potential differences and complicated deposition mechanisms. Here, in this work, we have successfully synthesized the near-equiatomic (<5 at% error) CoCuFeNi films with ~200nm thickness by electrodeposition from glycine-citrate-Triton X-100 acidic electrolytes. This system generally follows the Principle II in Brenner’s paradigm on alloy composition control in electrodeposition. X-ray diffraction (XRD) profiles show that the films only consist of one crystalline phase, different from the deposit from ammonia-citrate-boric acid electrolytes and the equilibrium phases predicted by CALPHAD. The near-equiatomic deposits at Ru substrates were successfully annealed at 400 ° C without significant intermixing between the substrate and the deposit, in which no phase separation of the crystalline phase was observed in its XRD profile. Xray photoelectron spectroscopy (XPS) reveals that some depositing elements (mainly Fe) exist as metal oxides. Cu is dissolved in the crystalline phases, stabilized by small crystalline domain size and the amorphous metal oxides inside the CoCuFeNi deposits.

CoCuFeNi↗
DOE OSTI2025

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗
DOE OSTI2025

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