Engineering PapersSearch

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.

72 records · Page 4

Characterization and Computational Modeling of Minor Phases in Alloy LSHR

The minor phases of powder metallurgy disk superalloy LSHR were studied. Samples were consistently heat treated at three different temperatures for long times to approach equilibrium. Additional heat treatments were also performed for shorter times, to assess minor phase kinetics in non-equilibrium conditions. Minor phases including MC carbides, M23C6 carbides, M3B2 borides, and sigma were identified. Their average sizes and total area fractions were determined. CALPHAD thermodynamics databases and PrecipiCalc(TradeMark), a computational precipitation modeling tool, were employed with Ni-base thermodynamics and diffusion databases to model and simulate the phase microstructural evolution observed in the experiments with an objective to identify the model limitations and the directions of model enhancement.

Jou, Herng-Jeng

Open-source Numerical Modeling of Solidification Cracking Susceptibility: Application to Refractory Alloy Systems

Introduction. Alloys such as aluminum, nickel-base, and austenitic stainless steels are susceptible to solidification cracking during welding and 3D printing. Compositional optimization is one method used to effectively mitigate solidification cracking of those alloy systems. With the surge in hypersonic and in-space propulsion activities, refractory metals (Nb, Mo, Ta, W, and Re) and their alloy derivatives are increasing in importance due to their extreme high melting point and retention of high-temperature strength; however, their chemistry was most typically optimized to promote ductility during mechanical operations such as drawing and forming. Welding of such alloys has been a challenge due to a number of issues including solidification cracking, atmospheric contamination (O, C, and N), as well as a shift in ductile-to-brittle transition to higher temperature following grain growth induced by welding. Compositional optimization of refractory alloys for solidification cracking resistance in particular is desirable as their usage increases with the advent of advanced manufacturing methods such as 3D printing. This work evaluates the effect of compositional variation in refractory metal systems on the solidification cracking susceptibility with the goals of optimizing existing alloys and joining process techniques, and formulating new alloys with increased solidification cracking resistance. Experimental Procedures. A python code was developed in a Jupyter notebook environment (Michael and Sowards, 2023) to facilitate the calculation of crack susceptibility index proposed by Kou (2015). Composition is entered as a single point, or as a 1-D or 2-D array. The notebook calls pycalphad (Otis and Liu, 2017 and Bocklund et al, 2020) to calculate the evolution of fraction solid as a function of temperature (under either Scheil or equilibrium assumptions) and then evaluates steepness of the fraction solid curve near the terminal stage of solidification to predict solidification cracking resistance. Open source thermodynamic databases available at online repositories are used (van de Walle). The process is setup in an automated fashion to generate plots that show variation in solidification cracking susceptibility according to composition on 1-D line plots or 2-D contour plots. The Jupyter notebook and crack susceptibility algorithm was also integrated with a widely used commercial CALPHAD code for validation and alloy exploration. Results and Discussion. The crack susceptibility model was first validated against a series of refractory alloy compositions evaluated in past work which utilized a specialized Varestraint test built inside a vacuum chamber environment (Lessman and Gold, 1971). The alloys tested in the Varestraint apparatus included T-111 (Ta-8W-2Hf), ASTAR-811C (Ta-8W-1Re-0.7Hf-0.025C), FS-85 (Nb-27Ta-10W-1Zr), T-222 (Ta-9.6W-2.4Hf-0.01C), Ta-10W, B-66 (Nb-5Mo-5V-1Zr), and SCb-291 (Nb-10W-10Ta). The initial test of the model showed a strong correlation with empirical Varestraint data, i.e., a Spearman rank correlation between model predictions and hot cracking measurements was observed to be greater than 0.8. Following the validation, a set of refractory metal binary mixtures was investigated to evaluate sensitivity of Nb, Mo, W, and Ta to C, N, and O content. A series of plots were produced that suggest ppmw ranges of C, N, and O where solidification cracking increases significantly and reaches a maximum. Also comparative ranking of each primary refractory metal to each interstitial was produced. For example C produces greater cracking response in Mo whereas O produces greater cracking response in Ta and Nb. Such compositional values have utility in setting limits on pickup of these interstitial elements during welding and printing rather than using a one-size-fits-all approach. Furthermore, the results have use in determining additive powder recycling requirements, which is especially pertinent for refractory metal powders due to their high cost compared to conventional alloys. Another application created thousands of hypothetical alloys within the nominal specified composition range of two widely used refractory alloys C103 (Nb-10Hf-1Ti) and TZM (Mo-0.5Ti-0.1Zr). The cracking index was calculated for the alloys and results were fed into machine learning regression techniques including Multiple Linear Regression, Ridge Regression, and Lasso Regression to determine relative potency each alloying element had on computed solidification cracking index. A series of linear equations were produced that relate composition of C103 and TZM to solidification cracking index. The crack susceptibility of C103 for example is described by an equation of the form: cracking index ~ O + 0.667*C + 0.635*N + 0.00037*Ta – 0.0008*Hf (in wt.%) From that equation, it is clear that O has strong propensity to induce solidification cracking. Interestingly, Hf is shown to reduce calculated cracking response. Finally, realizing the potential of this method to discover new refractory alloy formulations across the period table that have low solidification cracking sensitivity, the code was applied to new untested alloy systems including W-Zr-C, W-Ta-C, and others. Conclusions. In summary, an open source numerical method has been developed using Python code to calculate Kou’s crack susceptibility index. The method was applied to refractory metals which are inherently difficult to study from a weldability testing standpoint since inert shielding gas is not sufficient and welding is typically done in vacuum, especially in light of findings presented here where oxygen has profound influence on solidification cracking. This work revealed the effect of compositional variations on a series of refractory metals and showed the framework defined here will be useful in 1) the development of new alloys that have improved weldability and 3D printability, 2) placing compositional limits on existing alloys, and 3) ensuring adequate controls of manufacturing processes such as 3D printing where powder reuse is critical. Keywords. pycalphad; Python; refractory metals; solidification cracking. References. B. Bocklund et. al. (2020) http://doi.org/10.5281/zenodo.3630657. S. Kou. (2015) https://doi.org/10.1016/j.actamat.2015.01.034. G.G. Lessmann and R.E. Gold. Welding Journal, issue 1, pp. 1-s – 8-s (1971). F.N. Michael and J.W. Sowards. NASA/TM-20230002218 (2023). R. Otis and Z.-K. Liu. (2017) http://doi.org/10.5334/jors.140. A. Van de Wallle et. al. (2018) https://doi.org/10.1016/j.calphad.2018.04.003.

pycalphad

An Open-Source Numerical Model for Mitigating Refractory Alloy Hot Cracking Susceptibility

Refractory alloys are susceptible to solidification cracking during welding and 3D printing. Composition control is an effective method of controlling solidification cracking. This work evaluates the effect of compositional variation in refractory metal systems on a computed solidification cracking susceptibility. A numerical model has been developed using Python code and open-source CALPHAD software to calculate Kou’s crack susceptibility index. The model is validated against past weldability studies performed on several refractory alloy systems. The approach is extended towards the development of new alloys with improved 3D printability and weldability and is shown to have utility in defining compositional limits for existing alloys and feedstocks. Furthermore, the model will aid in determining process controls for powder reuse and recycling.

pycalphad

A Combined Computational, Experimental, and Technology Development Approach to In-Space Laser Manufacturing Maturation at NASA Marshall Space Flight Center

In-space manufacturing (ISM) is emerging as a field vital to continued access and capabilities in the space environment. NASA Marshall Space Flight Center (MSFC) is advancing the frontier of in-space laser manufacturing (ISLM) techniques through work initially focused on maturing laser beam welding (LBW) and laser forming (LF) for use in space. Such techniques proffer the ability to assemble and join structures in space from sheet metal or other stock – extant satellites, in situ resource utilization of Lunar regolith, etc. – by forming to desired shapes and then joining via in-space welding (ISW). ISLM processes are useful for assembly, joining, modification, and repair of structures in free space and on the Lunar surface such as large observatories, antennas, trusses, blast/thermal/radiation shields, pressure vessels, and more. However, these techniques are not yet qualified & certified (Q&C) for regular application in space. It would be prohibitively expensive, laborious, and time-consuming to perform Q&C via traditional experimental approaches as data collection & experimentation in space is resource-intensive. As such, benchmark experiments and focused, properly instrumented technology demonstration efforts in space can collect sufficient data that – when combined with verified computational models in an integrated computational materials engineering (ICME) approach – can validate ICME tools capable of translating more readily obtained ground data to in-space, in situ, computationally informed Q&C of ISLM techniques. Several ISLM projects at MSFC are obtaining the data required to validate ICME tools through both ground and flight experiments. A parabolic flight experiment of LBW under vacuum is manifested for August 2024, including both microgravity and Lunar gravity profiles. This collaboration with the Ohio State University is investigating common aerospace alloys such as 316L stainless steel, 2219 aluminum alloy, and Ti64 titanium alloy. In situ data collection includes videography, thermography, and reference thermocouples to build a thermal model of the welds. This will elucidate the relevant physics when combined with post-flight microstructural examination and mechanical testing. MSFC is also progressing towards a suborbital flight experiment of LBW under vacuum, which could provide reams of data on ISW during sustained, high-quality reduced gravity. The effect of combined thermal (cryogenic and high-temperature) and vacuum exposure on both LBW (NASA-funded) and LF (DARPA-funded) is being investigated through ground experiments. In addition to the copious data collected during these ground experiments, ruggedization of LBW hardware will also be pursued. The datasets from these experiments will be used to validate computational models which will inform future ISLM efforts in an ICME framework. A variety of techniques across lengths scales, from CALPHAD-driven thermodynamics & kinetics to phase field modeling of solidification to kinetic Monte Carlo simulations of grain evolution at the mesoscale, will be employed to accelerate the infusion and eventual Q&C of LBW and LF for use in space. The development of data-driven surrogate models to bridge ground to flight experiments and thereby reduce the need for resource-intensive experiments in space will also be investigated. These ICME techniques, surrogate models, and datasets from ground testing can also be employed to advance manufacturing in terrestrial environments.

in-space welding

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials

List of Commercial and Advanced Developmental Niobium-Based Alloys of the Space Age

This report compiles a list of commercial and developmental niobium-based alloys developed during the Space Age (late 1950s through mid-1970s) for extreme-temperature applications, including rocket engine thrust chambers, hypersonic re-entry thermal protection systems, and space fission reactor loops. Niobium (Nb) was widely pursued because it provided the lowest density (~8.6 g/cc) among the primary refractory metals, a high melting temperature (~2470°C), exceptional low-temperature ductility, good formability, and compatibility with liquid alkali metals. An evaluation of physical metallurgy mechanisms, focusing on solid-solution strengthening via heavy refractory solutes (W, Mo, Ta), dual-purpose reactive solutes (Hf, Zr, Ti), and dispersion strengthening using carbides, nitrides, and oxides is presented. Additionally, the report compares Western and Soviet Union metallurgical approaches, explaining how supply chain factors and manufacturing infrastructure influenced element selection, interstitial chemistry, and alloy identification/naming conventions. Cataloging these historical alloy chemical compositions serves as a foundational reference for modern alloy additive manufacturing, thermodynamic CALPHAD modeling, and machine-learning discovery pipelines for next-generation extreme-environment niobium-based alloys.

Physical Metallurgy

Design, Processing, and Properties of WTaCrV-Hf Multi-principal Element Alloys

Refractory multi-principal element alloys are candidates for high-temperature structural components due, in part, to their high strength and high melting points. Single-phase materials are initially preferred for isotropic material properties as a function of time and temperature in service conditions. This work outlines a computational rank-ordering and experimental validation methodology for single-phase body-centered-cubic phase stability in WTaCrV-Hf alloys using order–disorder transition temperature. Eight compositions were fabricated by arc-melting and heat-treated at 1400 °C for 24 hrs. X-ray diffraction, energy-dispersive x-ray spectroscopy, and Vickers hardness testing showed alloys with order–disorder transition temperatures below 600 °C formed a single-phase body-centered-cubic structure during solidification and remained single-phase after heat-treatment. The sample possessing the lowest order–disorder transition temperature exhibited slip traces suggestive of room-temperature plastic deformation under Vickers indentation, with both heat-treated single-phase samples exhibiting hardnesses over 800 HV with little cracking compared to tungsten. These results establish order–disorder transition temperature as a viable predictive parameter for multi-principal element alloy phase stability. The methodology outlined in this work provides a framework for future design, fabrication, and characterization of high-temperature structural multi-principal element alloys.

CALPHAD

Computational Design of Interlayers for Thermally Stable Compositionally Graded Coatings on Nickel Alloys

To extend the service life of Ni-based superalloys, refractory metal coatings are often used. However, direct bonding between metals with dissimilar crystal structure promotes brittle intermetallic phase formation. This work presents a computational thermodynamic framework for high throughput design of functionally graded interlayers to suppress deleterious phases that may form at the interlayer. The Thermo-Calc software package was used to screen candidate metallic interlayer elements based on stability of solid-solution phases. Vanadium was identified as a promising interlayer due to its consistent suppression of intermetallic phases. Temperature-dependent phase diagram mapping between 600 and 1000 °C guided selection of a compositional pathway that significantly reduced intermetallic formation compared to directly joining the Ni-based and Nb refractory alloys. Time–temperature–transformation analysis was performed to assess whether equilibrium-predicted phases are kinetically accessible along regions of the graded path where non-solid-solution phases are not fully suppressed. The methodology was further applied to additional Ni-based alloy and coating systems, illustrating its transferability as an approach for rapid computational design of graded interlayers in dissimilar high-temperature materials.

36 MATERIALS SCIENCE

Crystallographic stability of the BCC phase during quenching of metastable beta titanium alloy Ti-5553 and comparison of structural criteria for the predictive capability of α” martensite

This work evaluates the compositional standard currently developed for Ti-5553 powder and explores the metastable β region for sensitivity to martensitic formation during rapid quenching from the melt. Ti-5553 is a beta-stabilized titanium alloy that is increasingly being used in additive manufacturing applications. A series of alloys within the Ti-5553 compositional space was processed using two-piston spat quenching to perform rapid solidification and quenching for each of these alloy compositions. Typical empirical models such as molybdenum equivalency are not able to fully separate the retained BCC and martensitic compositions. Further, the use of thermodynamic data to estimate the transformation energy needed to form martensite can differentiate the alloys and provide a metric to further develop compositional limits for metastable beta titanium alloy development.

36 MATERIALS SCIENCE

Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system

The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.

CALPHAD

New insights on the effects of chemistry and temperature on α’ precipitation during aging of FeCrAl alloys

Here, the effect of Al on α' precipitation in FeCrAl alloys was studied through thermal aging of several model binary FeCr and FeCrAl alloys with Cr content of 13 wt.%, 17 wt.%, and 25 wt.% with and without 5.5 wt.% Al. Aging was performed up to 1,000h at varying temperatures from 400-500°C. At 400°C no age hardening was observed due to slow kinetics at this temperature. For both the 17Cr and 25Cr alloys, the addition of Al shows a lowering of the miscibility gap, consistent with other reports in literature. Interestingly, however, for the 25Cr alloys the addition of Al in the FeCrAl ternary alloy accelerated α' precipitation below the miscibility gap. Such enhanced precipitation has also been predicted by our atomistic kinetic Monte Carlo (AKMC) simulations. While previous literature has often focused on Al suppressing α', here we show that while Al can lower the miscibility gap in the Fe-Cr ferritic system thermodynamically, it may also enhance the kinetics of precipitation.

36 MATERIALS SCIENCE

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Thermodynamic properties of ZrCl 4 with LiCl, NaCl, KCl, CsCl, MgCl 2 , and UCl 3 for molten salt reactor applications

Here, a set of self-consistent Gibbs energy functions has been developed to describe the thermochemical behavior of the major reactive fission product ZrCl 4 with key components in chloride molten salt reactors (MSRs): LiCl–ZrCl 4 , NaCl–ZrCl 4 , KCl–ZrCl 4 , CsCl–ZrCl 4 , MgCl 2 –ZrCl 4 and UCl 3 –ZrCl 4 . Low ZrCl 4 concentration phase equilibria most relevant to MSR applications have been confirmed via differential scanning calorimetry for NaCl–ZrCl 4 , and X-ray diffraction analysis performed on equilibrated samples of NaCl–ZrCl 4 , KCl–ZrCl 4 , MgCl 2 –ZrCl 4 , and UCl 3 –ZrCl 4 . Within the framework of the modified quasi-chemical model in the quadruplet approximation, extrapolations of pseudo-binary models were generated to represent KCl–MgCl 2 –ZrCl 4 , KCl–NaCl–ZrCl 4 , and MgCl 2 –NaCl–ZrCl 4 , which show agreement with phase equilibria data where available without the use of ternary interaction parameters. The optimized thermodynamic descriptions for these systems and others are available in the open-source compendium Molten Salt Thermal Properties Database – Thermochemical (MSTDB–TC).

CALPHAD

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE

Computational and Physics-Based Modeling for the Development of in-Space Welding Technology

Welding will be an enabling technology for in-space assembly and manufacturing (ISAM) and in-space repair operations in support of the growing Space Economy. Prospective use cases of welding include metal repair on long-duration lunar surface missions and assembly of large structures that exceed dimensions of the currently available launch vehicle envelopes such as space arrays, tall lunar transmission towers, or even large transit vehicles powered by nuclear propulsion. There is a brief history of beam welding for in-space use. Soviet and NASA astronauts performed electron beam welds (EBW) in space in the late-1960s through mid-1980s and several parabolic flights of laser beam welding (LBW) were undertaken in the USA and Canada in the early-1990s. Work has since stagnated until the recent boom of the in-space economy has demonstrated need for continued development of EBW and LBW processes, which are desirable for their compactness, low mass, low power consumption, and high wall plug efficiency. No space-based manufacturing laboratory capable of demonstrating, developing, and applying welding processes exists, and costs are high to develop and fly individual space-based trials. This drives the current approach to develop welding technology with computational and physics-based models in tandem with earth-based experiments. Several approaches are taken to complement and drive the selection of instrumentation for planned parabolic and suborbital flights. First, a thermal model of Skylab EBW experiments is established and compared to empirical data. Next, the model is used to show variation between Earth-based and in-space weld structure and properties through exploration of key variables that define the Low earth orbit and lunar surface space environments, namely temperature (40 – 400 K), pressure (10-10 to 104 Pa), and gravity (10-6 g to 1 g). Finally, physical processes associated with heat and mass transport during welding are evaluated through several case studies across the above variable ranges demonstrating need for further research.

in-space welding