skip to main content


Title: A Phase-Field Model for In-Space Manufacturing of Binary Alloys
The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.  more » « less
Award ID(s):
2042683
NSF-PAR ID:
10392792
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Materials
Volume:
16
Issue:
1
ISSN:
1996-1944
Page Range / eLocation ID:
383
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Hydrogen-induced porosity formed during solidification of aluminum-based alloys has been a major issue adversely affecting the performance of solidification products such as castings, welds or additively manufactured components. A three-dimensional cellular automaton model was developed, for the first time, to predict the formation and evolution of hydrogen porosity coupled with grain growth during solidification of a ternary Al-7wt.%Si-0.3wt.%Mg alloy. The simulation results fully describe the concurrent nucleation and evolution of both alloy grains and hydrogen porosity, yielding the morphology of multiple grains as well as the porosity size and distribution. This model, successfully validated by X-ray micro-tomographic measurements and optical microscopy of a wedge die casting, provides a critical tool for minimizing/controlling porosity formation in solidification products.

     
    more » « less
  2. In fluid dynamics applications that involve flow adjacent to a porous medium, there exists some ambiguity in how to model the interface. Despite different developments, there is no agreed upon boundary condition that should be applied at the interface. We present a new analytical solution for laminar boundary layers over permeable beds driven by oscillatory free stream motion where flow in the permeable region follows Darcy's law. We study the fluid boundary layer for two different boundary conditions at the interface between the fluid and a permeable bed that was first introduced in the context of steady flows: a mixed boundary condition proposed by Beavers and Joseph [“Boundary conditions at a naturally permeable bed,” J. Fluid Mech. 30, 197–207 (1967)] and the velocity continuity condition proposed by Le Bars and Worster [“Interfacial conditions between a pure fluid and a porous medium: Implications for binary alloy solidification,” J. Fluid Mech. 550, 149–173 (2006)]. Our analytical solution based on the velocity continuity condition agrees very well with numerical results using the mixed boundary condition, suggesting that the simpler velocity boundary condition is able to accurately capture the flow physics near the interface. Furthermore, we compare our solution against experimental data in an oscillatory boundary layer generated by water waves propagating over a permeable bed and find good agreement. Our results show the existence of a transition zone below the interface, where the boundary layer flow still dominates. The depth of this transition zone scales with the grain diameter of the porous medium and is proportional to an empirical parameter that we fit to the available data. 
    more » « less
  3. The complex solidification cycles experienced by multi-principal element alloys (MPEAs) during laser-based additive manufacturing (LBAM) often lead to structural defects that affect the build quality. The underlying thermal processes and phase transformations are a function of the process parameters employed. With a moving Gaussian heat source to mimic LBAM and leveraging material thermodynamics guidelines from CALculation of PHAse Diagrams (CALPHAD), we estimate the temperature-dependent thermal properties, phase fractions, and melt pool geometry using an experimentally validated computational fluid dynamics model. The results substantiate that the peak temperatures are inversely correlated to the scan speeds, and the melt pool dimensions can assist in the predictive selection of process parameters such as hatch distance and layer thickness. A relatively low cooling rate recorded during the process is ascribed to the preheating of the substrate to ensure printability of the alloy.

     
    more » « less
  4. In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn– Hilliard–Navier–Stokes equations in the free flow region and Cahn–Hilliard–Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions. 
    more » « less
  5. The continued development of metal additive manufacturing (AM) has expanded the engineering metallic alloys for which these processes may be applied, including beta-titanium alloys with desirable strength-to-density ratios. To understand the response of beta-titanium alloys to AM processing, solidification and microstructure evolution needs to be investigated. In particular, thermal gradients (Gs) and solidification velocities (Vs) experienced during AM are needed to link processing to microstructure development, including the columnar-to-equiaxed transition (CET). In this work, in situ synchrotron X-ray radiography of the beta-titanium alloy Ti-10V-2Fe-3Al (wt.%) (Ti-1023) during simulated laser-powder bed fusion (L-PBF) was performed at the Advanced Photon Source at Argonne National Laboratory, allowing for direct determination of Vs. Two different computational modeling tools, SYSWELD and FLOW-3D, were utilized to investigate the solidification conditions of spot and raster melt scenarios. The predicted Vs obtained from both pieces of computational software exhibited good agreement with those obtained from in situ synchrotron X-ray radiography measurements. The model that accounted for fluid flow also showed the ability to predict trends unobservable in the in situ synchrotron X-ray radiography, but are known to occur during rapid solidification. A CET model for Ti-1023 was also developed using the Kurz–Giovanola–Trivedi model, which allowed modeled Gs and Vs to be compared in the context of predicted grain morphologies. Both pieces of software were in agreement for morphology predictions of spot-melts, but drastically differed for raster predictions. The discrepancy is attributable to the difference in accounting for fluid flow, resulting in magnitude-different values of Gs for similar Vs. 
    more » « less