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

Fundamental Issues Identified for Thermodynamic Description of Molten Salt Systems

Thermodynamic modeling of the molten salt systems is critical in the design of advanced nuclear power plants, materials recycling, and many other important clean energy applications. As one of the most capable computational thermodynamic approaches, the CALPHAD (Calculation of Phase Diagrams) method has been widely used for multicomponent thermodynamic database development and further support more physics-based modeling. Although considerable CALPHAD modeling efforts for molten salt systems have been made, we found that no evaluation has been performed regarding the unary thermodynamic models for the molten salt systems despite the fact that more than one pure substance database is available for the CALPHAD modeling. Therefore, the thermodynamic descriptions of the pure salts should be critically evaluated to ensure a high-fidelity molten salt thermodynamic database. In this work, we comprehensively analyzed some selected molten salts using different available thermodynamic descriptions from two databases as an example. One is the SSUB database released by the SGTE (Scientific Group Thermodata Europe), and the other is the FactPS database from the FactSage. Finally, the significant discrepancies observed in the thermodynamic description between the existing CALPHAD pure substance databases and experiments call for an urgent effort to improve.

42 ENGINEERING↗
DOE OSTI2024

Rapid Design of Cost-Effective Refractory High Entropy Alloys Strengthened by Precipitation

Funded by the ARPA-E ULTIMATE program, this project aims to develop cost-effective refractory high entropy alloys (RHEAs) for gas turbine blades applications at 1300 degree Celsius and above. The alloy design focuses on carbides precipitation strengthening in a strong and ductile solid solution BCC matrix. The goal is to achieve balanced mechanical properties including room temperature ductility & fracture toughness and high temperature strength and creep resistance, while maintaining comparable oxidation resistance and densities to Ni-base superalloys. To accelerate alloy design, multiscale modeling including CALPHAD and first-principles density functional theory (DFT) calculations and machine learning are performed. Specifically, high throughput phase diagram calculations and screening are carried out using CALPHAD; DFT calculations are performed to predict intrinsic ductility, grain boundary strengthening potential, coefficient of thermal expansion (CTE) and temperature-dependent elastic constants. Temperature dependent yield stress, creep strain and CTE are also predicted using machine learning. Down-selected alloys are synthesized in small buttons of about 250g using arc melting for rapid evaluation on microstructure and mechanical properties before employing plasma arc melting and additive manufacturing for producing large ingots. Preliminary computational and experimental results of this project will be presented and future direction in RHEA development will be discussed.

high temperature materials↗
DOE OSTI2021

Sensitivity of additively manufactured AA7075 to variation in feedstock composition and print parameters

Here we developed a combined finite element and CALPHAD based model of the Laser Powder Bed Fusion (LPBF) process for AA7075 alloy that considers the effect of feedstock composition and print parameters. A single-pass of a laser on a layer of AA7075 alloy powder has been considered. Sensitivity of temperature evolution and melt pool geometry to variation in the stoichiometry of the feedstock powder and laser source characteristics have been studied. Our results indicate that deviation (up to 10%) of the feedstock composition from the AA7075 raises the maximum temperature and increases melt pool size. Excess Cu content shows the largest melt pool width and depth among all the cases. The peak temperature is higher than the standard feedstock composition in all cases, except when the Cu concentration is reduced. Increasing the scan power also results in a higher peak temperature and a larger melt pool size. Furthermore, the temperature's rise time increases by lowering the scan speed.

42 ENGINEERING↗
DOE OSTI2023

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Predictive modeling of Néel temperature in austenitic alloys using CALPHAD and data analytics

The Néel temperature is a crucial yet often overlooked parameter in calculating the stacking fault energy (SFE) of austenitic alloys. Several empirical equations have been proposed to estimate the Néel temperature of austenitic alloys, which are then used to calculate the SFE and explain deformation mechanisms. However, these empirical equations, typically derived using linear regression algorithms, are often simplistic and may fail to capture the complex interactions among multiple alloying elements that influence the Néel temperature. Moreover, their applicability is usually limited to specific compositional ranges. In this study, we propose a CALPHAD based approach and develop a surrogate decision tree based regression model capable of capturing the interactions among multiple alloying elements to predict the Néel temperature. Predictions from both the CALPHAD approach and the regression model show close agreement with experimental measurements reported in the literature. In conclusion, the implications of accurate Néel temperature predictions on the calculated SFE and deformation mechanisms are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding the microstructural stability in a y’-strengthened Ni-Fe-Cr-Al-Ti alloy

Ni-Fe-Cr-Al-Ti alloys, with Ni levels near 50 wt.%, have the potential to develop a microstructure consisting of a face-centered cubic ? matrix with the homogeneous precipitation of fine ordered ?’ precipitates similar to traditional Ni-based superalloys with significantly greater Ni content. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), atom probe tomography (APT), and CALPHAD -based thermodynamic modeling were employed to understand the phase stabilities and microstructural evolution in an age-hardenable Ni-27Fe-18Cr-1Co-1.6Al-3.75Ti-1.2Mo-0.03C (wt%) alloy. The primary heat-treatment of solution annealing at 1121°C for 4h and age-hardening treatment at 760 °C for 16h resulted in a microstructure consisting of fine ?’ precipitates in an austenitic matrix along with grain boundary carbides, consistent with thermodynamic calculations. Long-term aging at 900 °C for 250h resulted in the coarsening of ?’ along with a change in the morphology of the precipitates from spherical to a more cuboidal shape. In addition, ? phase formation was observed concomitant with the partial dissolution of the ?’ phase. The ability of computational thermodynamic models to predict microstructural characteristics is discussed.

Gwalani, Bharat↗
DOE OSTI2023

Thermo4PFM: Facilitating Phase-field simulations of alloys with thermodynamic driving forces

Phase-field modeling is a popular front-tracking approach used to model solidification. Its time-evolution equations are often coupled to alloy composition and/or thermal diffusion in high-resolution multiphysics approaches. Materials thermodynamic properties tabulated in CALPHAD databases can be used for phase-field modeling to parameterize bulk energies of alloys. In addition, they can be naturally integrated into models such as the Kim-Kim-Suzuki (KKS) model where driving forces depend on the differences between chemical potentials of co-existing phases. In that case, a small system of coupled nonlinear equations needs to be solved at every point in space where the phase-field order parameter is to be updated and evolved in time. Here we present Thermo4PFM, a solver for the KKS equations for binary and ternary alloys, with two or three phases, and parameterized with CALPHAD models. Thermo4PFM is open source, written in C++, and can take advantage of Graphics Processing Units (GPU) accelerators. Using OpenMP offload capabilities for C++ classes, an excellent performance is demonstrated on GPU using the LLVM compiler. CALPHAD data is read from simple JSON files using an open source parser from the boost library.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermodynamic modeling of the Pd–Zn system with uncertainty quantification and its implication to tailor catalysts

Pd–Zn intermetallic catalysts show encouraging combinations of activity and selectivity on well-defined active site ensembles. Thermodynamic description of the Pd–Zn system, delineating phase boundaries and enumerating site occupancies within intermediate alloy phases, is essential to determining the ensembles of Pd–Zn atoms as a function of composition and temperature. Combining the present extensive first-principles calculations based on density functional theory (DFT) and available experimental data, the Pd–Zn system was remodeled using the CALculation of PHAse Diagrams (CALPHAD) approach. High throughput modeling tools with uncertainty quantification, i.e., ESPEI and PyCalphad, were incorporated in the phase analysis. Here, the site occupancies across the γ phase composition region were given special attention. A four-sublattice model was used for the γ phase owing to its four Wyckoff positions, i.e., the outer tetrahedral (OT) site 8c, the inner tetrahedral (IT) site 8c, the octahedral (OH) site 12e, and the cuboctahedral (CO) site 24g. The site fractions of Pd and Zn calculated from the present thermodynamic model show the occupancy preference of Pd in the OT and OH sublattices in agreement with experimental observations. The force constants obtained from DFT-based phonon calculations further supports the tendency of Pd occupying the OH sublattice compared with the IT and CO sublattice. The catalytic ensembles changing from Pd monomers (Pd 1 ) to trimers (Pd 3 ) on the surface of γ phase are attributed to the increase of Pd occupancy in the OH sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Novel Hybrid Ab Initio Model of High-Performance Structural Alloys Design for Fossil Energy Power Plants

The current project developed a novel HT-CALPHAD/DFT approach, which can quickly design new high-performance structural alloys for the application of FE power plants. The PI will mainly take charge of high-throughput DFT simulations and computational thermodynamics of the selected multicomponent alloy systems for the FE power plant applications. At the end of the project, a novel hybrid model based on high throughput CALPHAD/DFT simulations and computational thermodynamics will be developed to provide guidance on how to identify multi-component new high-performance structural alloys with much less computational effort needed. It will address the extensive computation time needed for DFT on the new alloys design. In addition, it will also address the well-known headache of DFT, i.e. how to make the accurate prediction of the high-temperature equilibria. The novel hybrid model the PI proposed will not only be applied to the design of high-performance structural alloys in FE power plants but in many different applications, such as nuclear reactors. This hybrid modeling approach includes four sections: 1) HT-CALPHAD modeling of Al-Co-Cr-Ni-Fe system with FCC and BCC phase. In this section, 3561 non-equiatomic compositions were randomly generated in order to investigate the phase stability of single FCC and BCC phases. Meanwhile, we proposed a data screening procedure to screen out the good candidates within these compositions, considering the temperature range, average density, and melting temperature, etc. 2) Investigation of FCC-Cr lattice stability in Fe-Cr system. We systematically assessed the reliability of FCC-Cr lattice stability derived by DFT and CALPHAD approaches. Meanwhile, the Fe-Cr binary system was chosen as a case study to verify the Cr lattice stability obtained by both approaches. 3) High-throughput DFT modeling on elastic properties of Al-Co-Cr-Ni-Fe systems. We predicted and established the FCC quinary elastic constant database of the Al-Co-Cr-Fe-Ni systems at 0K by using special quasi-random structure (SQS) approach. The predictions will start with pure elements of Al-Co-Cr-Fe-Ni system and will be continued with binaries, ternaries, quaternaries, and finally the quinary compositions. 3) Modeling of temperature-dependent elastic properties in Al-Co-Cr-Ni-Fe systems. In this part, we predict the thermal expansion coefficient and elastic stiffness coefficient as a function of temperature by applying quasiharmonic approximation. With this approach, the elastic properties of HEAs at elevate temperature can be estimated.

01 COAL, LIGNITE, AND PEAT↗
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 OSTI2023

Impact of titanium content on the thermo-mechanical and oxidation response of TiAlTa

This study aimed to analyze how the addition of titanium (Ti) to TiAlTa alloys affects their mechanical properties (Young’s modulus and hardness) and oxidation behavior. Alloys with three different Ti additions and equi-atomic Al and Ta were fabricated with nominal compositions of 33Ti-33Al-33Ta at%, 50Ti-25Al-25Ta at%, and 70Ti-15Al-15Ta at%. Phase identification of as-processed and compositionally homogenized alloys was conducted by coupling various electron microscopy, diffraction, and chemical analysis techniques. Composition and structures were modeled using CALPHAD-based phase predictions and compared with experimental results. Each of the alloys exhibited a unique set of microstructures with distinctly different ordered precipitates, dependent on cooling rate and chemical composition. Subsequently, nanoindentation tests were performed at temperatures of 25°C, 250°C, 500°C, and 750°C. Oxidation studies were conducted under static air at 750°C for up to 200 h. 33Ti-33Al-33Ta at% and 70Ti-15Al-15Ta at% alloys showed higher hardness and modulus than 50Ti-25Al-25Ta at% at room temperature. The 50Ti-25Al-25Ta at% alloy exhibited higher hardness and modulus values but lower oxidation resistance at 750°C when compared with the other specimens. The X-ray photoelectron spectroscopy experiments were conducted to further analyze the oxides present in the samples. Despite each sample displaying combinations of mixed oxide layers, TiO 2 was the preferred oxide forming in the 50Ti-25Al-25Ta composition, while the other samples preferred to form Al 2 O 3 . In conclusion, the findings from this systematic exploration of phase evolution in the TiAlTa alloy space provide insights for optimized material chemistries, processing, and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Design and development of a high-performance Ni-based superalloy WSU 150 for additive manufacturing

This work proposes a design and development strategy of a new nickel-based superalloy for additive manufacturing. A new Ni-based superalloy has been developed by the application of the combinatorial alloy development technique coupled with CALPHAD based solidification modeling by effectively suppressing the precipitation kinetics. The suppression of precipitation during processing paved a way for prevention of cracks during deposition. The new alloy "WSU 150" revealed excellent room temperature mechanical properties with a yield strength of 867 MPa, an ultimate tensile strength of 1188 MPa, and an elongation of 27.9% in the as-deposited condition. The mechanical properties of the heat-treated alloy were improved significantly to 1114 MPa yield strength, 1396 MPa ultimate tensile strength, and an elongation of 16.1%. Improvement in the mechanical properties is attributed to the additional precipitation and coarsening of γ' and carbides during heat-treatment. Microstructural investigation of the alloy revealed spherical γ' with a rippled size distribution from the center to the interdendritic region. The average size of the γ' particles in the as-deposited condition was found to be around 48 nm in the interdendritic region. Heat-treatment promoted the coarsening of γ' which is explained in the paper.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High throughput structure–property relationship for additively manufactured 316L/IN625 alloy mixtures leveraging 2-step Bayesian estimation

While the fabrication of graded materials by directed energy deposition (DED) has led to accelerated materials discovery, the ability to rapidly explore sufficiently large material composition spaces is limited due to the time-intensive nature of conventional materials characterization techniques. The present study investigates the viability of small punch test (SPT) protocols for rapidly evaluating DED-fabricated alloy mixtures of stainless steel 316L (316L) and Inconel 625 (IN625). The SPT protocols evaluated in this study include both the recently established two-step Bayesian estimation framework as well as the empirical relationships established in prior literature. It is shown that these protocols are capable of reliably and quantitatively tracking the changes in the mechanical properties of the alloy mixtures studied. Enhancement of mechanical properties was observed with the addition of IN625 to 316L, which is attributed to the austenite stabilization in the matrix and the formation of fine δ - Ni3Nb precipitates. It is shown that CALPHAD-based Scheil model simulations predicted the formation of different precipitate phases for each composition. The novel protocols presented in this paper open new avenues for high throughput material explorations for additive manufacturing.

2-step Bayesian estimation↗
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 OSTI2021

Developing Practical Models of Complex Salts for Molten Salt Reactors

Molten salt reactors (MSRs) utilize salts as coolant or as the fuel and coolant together with fissile isotopes dissolved in the salt. It is necessary to therefore understand the behavior of the salts to effectively design, operate, and regulate such reactors, and thus there is a need for thermodynamic models for the salt systems. Molten salts, however, are difficult to represent as they exhibit short-range order that is dependent on both composition and temperature. A widely useful approach is the modified quasichemical model in the quadruplet approximation that provides for consideration of first- and second-nearest-neighbor coordination and interactions. Its use in the CALPHAD approach to system modeling requires fitting parameters using standard thermodynamic data such as phase equilibria, heat capacity, and others. A shortcoming of the model is its inability to directly vary coordination numbers with composition or temperature. Another issue is the difficulty in fitting model parameters using regression methods without already having very good initial values. The proposed paper will discuss these issues and note some practical methods for the effective generation of useful models.

Besmann, Theodore M. (ORCID:0000000155980550)↗
DOE OSTI2022

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tailoring critical Al concentration to form external Al 2 O 3 scale on Ni–Al alloys by computational approach

Abstract Nickel (Ni)‐based superalloys for high‐temperature applications are often designed to form a continuous and slow‐growing oxide scale by adding Al and Cr and other beneficial elements. In the present work, the critical Al concentration in Ni–Al alloys needed to establish an α‐Al 2 O 3 scale in contrast to internal oxide formation is predicted as a function of temperature by means of the CALPHAD approach coupled with models in the literature, which account for the thermodynamics and kinetics of oxidation. The present thermodynamic remodeling of the Ni–O system results in a better agreement with experimental data of oxygen solubility in Ni at high temperatures. The oxygen solubility is combined with kinetic parameters to determine oxygen permeability in Ni, and the critical Al concentration needed to establish an α‐Al 2 O 3 scale at a given exposure temperature. Good agreement is found with available experimental data for both oxygen permeability and critical Al concentration, indicating the capacity of the CALPHAD approach to tailor oxidation resistance for materials of interest using thermodynamic and kinetic knowledge.

Ross, Austin↗