Search NASASearch

SEARCH · Search NASA

Results for “Solidification”

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

Deep operator network surrogate for phase-field modeling of metal grain growth during solidification

A deep operator network (DeepONet) has been constructed that generates accurate representations of phase-field model simulations for evolving two dimensional metal grain morphology growing from melt. These representations serve as lower resolution, computationally efficient stand-ins for quick parameter space exploration of solutions to the the Allen-Cahn equations that dictate the phase-field model simulations. The experimental target for the phase-field model is a uranium casting system cooling a 434 g uranium charge from a maximum temperature of 1400° C at an average rate of 30° C / min , traversing the crystallographic phases of the pure metal. Experimental parameters inform the phase-field model, whose higher resolution computational model solutions are used to train the DeepONet in a given parameter space with the aim of developing a faster, more efficient method for predicting the solidifying metal's microstructure at different potential experimental values. The final DeepONet generates high accuracy, lower resolution predictions with cumulative relative approximation error over all timesteps of less than 0.5%, while ensuring solutions remain within physically feasible ranges. Further, these relative error values are comparable with other state-of-the-art DeepONet models for microstructure evolution, while significantly reducing the amount of training data required. Training a convolutional neural network simultaneously with the DeepONet, enforcing realistic values at the complex metal grain boundaries, and mathematically encoding boundary conditions into the structure of the DeepONet improved prediction accuracy and computational efficiency over a standard DeepONet model.

36 MATERIALS SCIENCE

Completion of lunar magma ocean solidification at 4.43 Ga

Crystallization of the lunar magma ocean yielded a chemically unique liquid residuum named KREEP. This component is expressed as a large patch on the near side of the Moon and a possible smaller patch in the northwest portion of the Moon’s South Pole-Aitken basin on the far side. Thermal models estimate that the crystallization of the lunar magma ocean (LMO) could have spanned from 10 and 200 My, while studies of radioactive decay systems have yielded inconsistent ages for the completion of LMO crystallization covering over 160 My. Here, we show that the Moon achieved >99% crystallization at 4,429 ± 76 Ma, indicating a lunar formation age of ~4,450 Ma or possibly older. Using the 176 Lu– 176 Hf decay system (t 1/2 = 37 Gy), we found that the initial 176 Hf/ 177 Hf ratios of lunar zircons with varied U–Pb ages are consistent with their crystallization from a KREEP-rich reservoir with a consistently low 176 Lu/ 177 Hf ratio of 0.0167 that emerged ~140 My after solar system formation. The previously proposed younger model age of ~4.33 Ga for the source of mare basalts (240 My after solar system formation) might reflect the timing of a large impact. Our results demonstrate that lunar magma ocean crystallization took place while the Moon was still battered by planetary embryos and planetesimals leftover from the main stage of planetary accretion. The study of Lu–Hf model ages for samples brought back from the South Pole-Aitken basin will help to assess the lateral continuity of KREEP and further understand its significance in the early history of the Moon.

Science & Technology - Other Topics

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials using MALAMUTE

The Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy aims to develop and qualify additively manufactured materials for nuclear applications. One key challenge to this is the microstructural variability observed in the additively manufactured products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-throughput experimental and modeling techniques to accelerate qualification. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the additive manufacturing process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture microstructural variability is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning models to develop a digital twin for additive manufacturing that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation. The melting and subsequent solidification that occurs during the additive process is a complex phenomenon that requires multiscale multiphysics analysis. This work package focuses on understanding the role of process variabilities on the unique microstructural characteristics of additively manufactured materials. Microstructural features at the subgrain level, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. Idaho National Laboratory’s Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for additively manufactured materials in an efficient, reliable, and cost-effective way. This work focuses on capturing the microstructural variabilities at the subgrain level that are often missing in the part-scale models. In fiscal year 2025, we significantly advanced upon our work in the last fiscal year, in terms of the predictive capabilities of the physics-based and ML models, by adding the capabilities to capture subgrain-level micro-segregation during solidification using phase-field model and to predict the time-dependent dynamics of the AM process through the MOGPAR model. The alloy solidification model in MOOSE incorporates the thermodynamic properties and free energy relevant to 316 stainless steel. The model demonstrates the Cr and Ni segregation that occurs during solidification, including that the rate of solidification. The microstructural evolution model is connected to the process conditions via the surrogate model developed in this work. This enables predictions of the final microstructure in conjunctions with the manufacturing process. This work supports AMMT's rapid qualification goals by laying the foundation for an efficient and cost-effective model establishing the PSPP correlation for AM. The generated microstructures and predicted micro-segregation can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work helps to identify the key microstructural features at the subgrain level that are significant in property and performance predictions of additively manufactured components. This work will also provide inputs to the large-scale process variability models to reevaluate and validate assumptions and simplifications made in the part-scale models. Furthermore, through active learning this work can help identify the data need from both modeling and experimental sides for development of a robust digital twin for additive manufacturing and accelerate the AMMT's qualification efforts.

36 - MATERIALS SCIENCE

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE

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

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic

Effect of vanadium on phase composition and hard magnetic properties of as-solidified and heat-treated Sm–Fe–(Ti,V) alloys

Although the intrinsic magnetic properties of Ti-stabilized Sm(Fe,Co,Ti) 12 compounds exhibit potential of excellent rare-earth-lean permanent magnets, it has been much easier to realize large coercivities with the isostructural compounds stabilized by either V or by certain combinations of Ti and V. To elucidate the influence of V on the microstructure and magnetic properties, a series of Sm 8.1 Fe 78.4 (Ti 1-x V x ) 13.5 alloys was studied after arc-melting and annealing at 850–1000 °C. The alloys were found to fall into three groups. For x ≤ 0.2, solidification generates mostly the Sm(Fe,Ti,V) 12 phase, but annealing converts at least part of it into the non-magnetic Sm(Fe,Ti,V) 11 and the magnetically soft Sm 2 (Fe,Ti,V) 17 phases. For 0.2 < x < 0.6, the alloys solidify into a near-equilibrium mixture of the Sm(Fe,Ti,V) 12 , TiFe 2 and Sm-rich phases. For x ≥ 0.6, solidification generates large fractions of α-Fe solid solution and Sm-rich phases; an annealing step is necessary to complete the formation of Sm(Fe,Ti,V) 12 phase. Also, for x ≥ 0.6 the temperature below which the Sm(Fe,Ti,V) 12 phase is stable decreases with x, as does the fraction of this phase formed during solidification. Here, the differences between these three groups of alloys suggest different strategies for developing hard magnetic properties, with the likelihood of a success increasing with increasing x. For x ≥ 0.6, heat treatment alone is demonstrated to generate a microstructure of micron and submicron Sm(Fe,Ti 1-x V x ) 12 crystallites separated by a Sm-rich phase and exhibiting a coercivity with values up to 3.5 and 5.7 kOe for x = 0.8 and 1.0, respectively.

36 MATERIALS SCIENCE

LPBF Processability of NiTiHf Alloys: Systematic Modeling and Single-Track Studies

Research into the processability of NiTiHf high-temperature shape memory alloys (HTSMAs) via laser powder bed fusion (LPBF) is limited; nevertheless, these alloys show promise for applications in extreme environments. This study aims to address this limitation by investigating the printability of four NiTiHf alloys with varying Hf content (1, 2, 15, and 20 at. %) to assess their suitability for LPBF applications. Solidification cracking is one of the main limiting factors in LPBF processes, which occurs during the final stage of solidification. To investigate the effect of alloy composition on printability, this study focuses on this defect via a combination of computational modeling and experimental validation. To this end, solidification cracking susceptibility is calculated as Kou’s index and Scheil–Gulliver model, implemented in Thermo-Calc/2022a software. An innovative powder-free experimental method through laser remelting was conducted on bare NiTiHf ingots to validate the parameter impacts of the LPBF process. The result is the processability window with no cracking likelihood under diverse LPBF conditions, including laser power and scan speed. This comprehensive investigation enhances our understanding of the processability challenges and opportunities for NiTiHf HTSMAs in advanced engineering applications.

36 MATERIALS SCIENCE

Understanding the influence of boron in additively manufactured GammaPrint®-700 CoNi-based superalloy

Boron is commonly added to superalloys in small amounts to enhance creep resistance, but can lead to cracking at high concentrations, especially during the additive manufacturing process. Two variants of CoNi-based GammaPrint®-700 superalloy with different B contents (0.08 at% vs 0.16 at%) were printed via laser powder bed fusion (LPBF) with the same printing parameters, with only the high B alloy exhibiting solidification cracking. Atom probe tomography (APT) revealed stronger segregation behaviors in the high B alloy compared to the low B alloy at both the inter-dendritic regions and grain boundaries (GBs). The segregation behavior at inter-dendritic regions was well captured with Scheil simulation and can correlate with the existing cracking susceptibility index (CSI) on cracking tendencies, although high angle GBs are where cracking occurs according to electron backscatter diffraction (EBSD) measurements. Additionally, the extent of GB segregation was compared between the high B and low B alloy. Higher B additions led to significantly more GB B segregation in the high B alloy compared to the low B alloy. Further, for the high B alloy, the cracked region of one GB exhibited higher levels of B compared to the uncracked region of the same GB. However, much higher B contents were also found in two other uncracked GBs in the high B alloy, which demonstrates that higher GB B concentrations are not fully responsible for the cracking. A much larger variance in GB B segregation content was found in the high B alloy compared to the low B alloy. These phenomena were explained with a solidification model with the GB segregation content expressed explicitly by a modified Langmuir-McLean equation. This model linked the GB segregation content with solidification undercooling, which can be used as quantitative cracking criteria for future builds.

36 MATERIALS SCIENCE

A high strength Al-2Ni-0.5Zr conductor alloy fabricated via laser powder bed fusion

There is a current need for new aluminum alloy design strategies to target applications requiring high strength and conductivity with reductions in mass. A new lightweight Al-2Ni-0.5Zr (wt. %) conductor alloy was fabricated using laser powder bed fusion. A design of experiments probed the alloy's solidification cracking susceptibility. It was observed that solidification cracking was generally reduced with fast scan speeds, above 1500 mm/s, and smaller hatch spacings. The different cooling rates throughout the melt pool produced a heterogeneous distribution of cellular and equiaxed Al 3 Ni precipitates in the as-printed alloy. Additionally, the rapid solidification characteristic of laser powder bed fusion created a super-saturated Zr solid solution. An aging heat treatment at 375 °C for 24 h imparted strengthening through the precipitation of L1 2 -Al 3 Zr nanoprecipitates, which counteracted the softening caused by the fragmentation and coarsening of Al 3 Ni precipitates. The yield strength increased from 138 MPa in the as-printed condition to 168 MPa after aging, while the ductility remained constant at ∼21%. The aging treatment simultaneously increased the electrical conductivity from 40.8% IACS (International Annealed Copper Standard) to 53.5% IACS. Modeling of the strengthening mechanisms and electrical conductivity contributions rationalized the simultaneous increase in strength and conductivity upon aging. Furthermore, the strengthening efficacy of the Al 3 Ni and L1 2 -Al 3 Zr precipitates, combined with the low Ni and Zr solubility in the FCC Al matrix, facilitated both high strength and electrical conductivity. Overall, the combination of strength and electrical conductivity positions this alloy as a suitable choice for additively manufactured lightweight conductors.

Additive manufacturing

Ubiquitous short-range order in multi-principal element alloys

Recent research in multi-principal element alloys (MPEAs) has increasingly focused on the role of short-range order (SRO) on material performance. However, the mechanisms of SRO formation and its precise control remain elusive, limiting the progress of SRO engineering. Here, leveraging advanced additive manufacturing techniques that produce samples with a wide range of cooling rates (up to 10 7 K s –1 ) and an enhanced semi-quantitative electron microscopy method, we characterize SRO in three CoCrNi-based face-centered-cubic (FCC) MPEAs. Surprisingly, irrespective of the processing and thermal treatment history, all samples exhibit similar levels of SRO. Atomistic simulations reveal that during solidification, prevalent local chemical order arises in the liquid-solid interface (solidification front) even under the extreme cooling rate of 10 11 K s –1 . This phenomenon stems from the swift atomic diffusion in the supercooled liquid, which matches or even surpasses the rate of solidification. Therefore, SRO is an inherent characteristic of most FCC MPEAs, insensitive to variations in cooling rates and even annealing treatments typically available in experiments.

36 MATERIALS SCIENCE

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science

Modelling of cracking during beam oscillation laser welding of a creep-resistant nickel alloy

During high power laser welding of Inconel 740H, horizontal solidification cracks form at locations between approximately 70–80% of the weld depth. With the integration of circular beam oscillation, the laser energy density distribution and underlying thermo-mechanical conditions were altered. By coupling three-dimensional heat transfer, fluid flow, and stress modelling tools, the role that circular beam oscillation plays in the formation of the strain rates and stresses driving this cracking phenomenon was identified. While the addition of a circular oscillation pattern appeared to lower the stress levels along the solidification front, it did not eliminate the appearance of horizontal cracking. Only when the beam amplitudes reached 1.6 mm did solidification cracking disappear, but the weld depths were also significantly reduced.

Gao, Mingze