A mixed interface-capturing/interface-tracking formulation for thermal multi-phase flows with emphasis on metal additive manufacturing processes

Abstract High fidelity thermal multi-phase flow simulations are in much demand to reveal the multi-scale and multi-physics phenomena in metal additive manufacturing (AM) processes, yet accurate and robust predictions remain challenging. In this paper, we present a novel computational framework by mixing interface-capturing/interface-tracking methods for simulating the thermal multi-phase flows in metal AM applications, focusing on better handling the gas-metal interface, where AM physics, such as phase transitions and laser-material interactions, mainly takes place. The framework, built on level set method and variational multi-scale formulation (VMS), features three major contributions: (1) a simple computational geometry-based re-initialization approach, which maintains excellent signed distance property on unstructured meshes, re-constructs an explicit representation of gas-metal interface from the level set, and facilitates the treatment of the multiple laser reflections during keyhole evolution in AM processes; (2) a fully coupled VMS formulation for thermal multi-phase governing equations, including Navier-Stokes, level set convection, and thermodynamics with melting, solidification, evaporation, and interfacial force models; and (3) a three-level recursive preconditioning technique to enhance the robustness of linear solvers. We first compare the geometry-based re-initialization with the Eikonal partial differential equation (PDE)-based approach on two benchmark problems on level set convection and bubble dynamics. The comparison shows the geometry-based approach attains equivalent and even better performance on key criteria than the PDE-based counterpart. We then apply the developed framework to simulate two AM experiments, which Argonne National Laboratory has recently conducted using in-situ high-speed, high-energy x-ray imaging. The proposed framework’s accuracy is assessed by thoroughly comparing the simulated results against experimental measurements on various quantities. We also report important quantities that experiments cannot measure to show the modeling capability.

[1]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[2]  Victor M. Calo,et al.  Improving stability of stabilized and multiscale formulations in flow simulations at small time steps , 2010 .

[3]  Ignacio Romero,et al.  A thermodynamically consistent numerical method for a phase field model of solidification , 2014, Commun. Nonlinear Sci. Numer. Simul..

[4]  D. Jacqmin Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[5]  Interface-Capturing Method for Free-Surface Plunging and Breaking Waves , 2019, Journal of Engineering Mechanics.

[6]  Sheldon Wu,et al.  Modulating laser intensity profile ellipticity for microstructural control during metal additive manufacturing , 2017 .

[7]  Rama Govindarajan,et al.  Dynamics of an initially spherical bubble rising in quiescent liquid , 2015, Nature Communications.

[8]  Orion L. Kafka,et al.  Multi-scale modeling of electron beam melting of functionally graded materials , 2016 .

[9]  G. Buscaglia,et al.  A geometric mass-preserving redistancing scheme for the level set function , 2009 .

[10]  Wing Kam Liu,et al.  Data-driven multi-scale multi-physics models to derive process–structure–property relationships for additive manufacturing , 2018 .

[11]  Hector Gomez,et al.  Liquid-vapor transformations with surfactants. Phase-field model and Isogeometric Analysis , 2016, J. Comput. Phys..

[12]  H. Gómez,et al.  Flow and mixing dynamics of phase-transforming multicomponent fluids , 2019, Applied Physics Letters.

[13]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[14]  A. Masud,et al.  Variationally derived interface stabilization for discrete multiphase flows and relation with the ghost-penalty method , 2021 .

[15]  A. B. Strong,et al.  Vorticity based turbulence model for thermofluids modelling of welds , 2003 .

[16]  Suck-Joo Na,et al.  A study on ray tracing method for CFD simulations of laser keyhole welding: progressive search method , 2016, Welding in the World.

[17]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[18]  Kensuke Yokoi,et al.  A density-scaled continuum surface force model within a balanced force formulation , 2014, J. Comput. Phys..

[19]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

[20]  Tayfun E. Tezduyar,et al.  Finite element stabilization parameters computed from element matrices and vectors , 2000 .

[21]  Thomas J. R. Hughes,et al.  Liquid–vapor phase transition: Thermomechanical theory, entropy stable numerical formulation, and boiling simulations , 2015 .

[22]  Wentao Yan,et al.  Evaporation Model for Keyhole Dynamics During Additive Manufacturing of Metal , 2020, Physical Review Applied.

[23]  D. Gu,et al.  Laser energy absorption behavior of powder particles using ray tracing method during selective laser melting additive manufacturing of aluminum alloy , 2018 .

[24]  Y. Bazilevs,et al.  Isogeometric analysis of multi-phase flows with surface tension and with application to dynamics of rising bubbles , 2019, Computers & Fluids.

[25]  A. Marsden,et al.  The nested block preconditioning technique for the incompressible Navier-Stokes equations with emphasis on hemodynamic simulations , 2019, Computer methods in applied mechanics and engineering.

[26]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[27]  Wentao Yan,et al.  Spattering and denudation in laser powder bed fusion process: Multiphase flow modelling , 2020 .

[28]  Moataz M. Attallah,et al.  Mesoscale modelling of selective laser melting: Thermal fluid dynamics and microstructural evolution , 2017 .

[29]  P. Guillaume,et al.  Modeling of laser beam and powder flow interaction in laser cladding using ray-tracing , 2015 .

[30]  A. Masud,et al.  Residual‐based turbulence models for moving boundary flows: hierarchical application of variational multiscale method and three‐level scale separation , 2013 .

[31]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[32]  Kamel Fezzaa,et al.  Keyhole threshold and morphology in laser melting revealed by ultrahigh-speed x-ray imaging , 2019, Science.

[33]  William E. Frazier,et al.  Metal Additive Manufacturing: A Review , 2014, Journal of Materials Engineering and Performance.

[34]  J. Strain Fast Tree-Based Redistancing for Level Set Computations , 1999 .

[35]  I. Harari,et al.  Modeling of steep layers in singularly perturbed diffusion–reaction equation via flexible fine-scale basis , 2020 .

[36]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

[37]  Ju Liu,et al.  Functional entropy variables: A new methodology for deriving thermodynamically consistent algorithms for complex fluids, with particular reference to the isothermal Navier-Stokes-Korteweg equations , 2013, J. Comput. Phys..

[38]  Gunilla Kreiss,et al.  A conservative level set method for two phase flow II , 2005, J. Comput. Phys..

[39]  Stephen Lin,et al.  A volume-conserving balanced-force level set method on unstructured meshes using a control volume finite element formulation , 2019, J. Comput. Phys..

[40]  J. Hattel,et al.  Keyhole-induced porosities in Laser-based Powder Bed Fusion (L-PBF) of Ti6Al4V: High-fidelity modelling and experimental validation , 2019, Additive Manufacturing.

[41]  Thomas J. R. Hughes,et al.  Multiscale and Stabilized Methods , 2007 .

[42]  Feng Lin,et al.  Meso-scale modeling of multiple-layer fabrication process in Selective Electron Beam Melting: Inter-layer/track voids formation , 2018 .

[43]  Aiden A. Martin,et al.  Controlling interdependent meso-nanosecond dynamics and defect generation in metal 3D printing , 2020, Science.

[44]  Stephen Lin,et al.  A fully coupled finite element formulation for liquid–solid–gas thermo-fluid flow with melting and solidification , 2018, Computer Methods in Applied Mechanics and Engineering.

[45]  Y. Shin,et al.  Investigation of keyhole plume and molten pool based on a three-dimensional dynamic model with sharp interface formulation , 2013 .

[46]  G. Tryggvason,et al.  Computations of film boiling. Part I: numerical method , 2004 .

[47]  Tayfan E. Tezduyar,et al.  Stabilized Finite Element Formulations for Incompressible Flow Computations , 1991 .

[48]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[49]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[50]  Muriel Carin,et al.  A complete model of keyhole and melt pool dynamics to analyze instabilities and collapse during laser welding , 2014 .

[51]  T. Tezduyar,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. I: The concept and the preliminary numerical tests , 1992 .

[52]  T. Hughes,et al.  Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows , 2007 .

[53]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[54]  I. Akkerman,et al.  Isogeometric analysis of free-surface flow , 2011, J. Comput. Phys..

[55]  A. Rubenchik,et al.  Laser powder-bed fusion additive manufacturing: Physics of complex melt flow and formation mechanisms of pores, spatter, and denudation zones , 2015, 1512.02593.

[56]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[57]  Thomas J. R. Hughes,et al.  The multiscale formulation of large eddy simulation: Decay of homogeneous isotropic turbulence , 2001 .

[58]  Feng Lin,et al.  Multi-physics modeling of single/multiple-track defect mechanisms in electron beam selective melting , 2017 .

[59]  Taehun Lee,et al.  Single bubble rising dynamics for moderate Reynolds number using Lattice Boltzmann Method , 2010 .

[60]  Pierre Saramito,et al.  Improving the mass conservation of the level set method in a finite element context , 2010 .

[61]  A. Korobenko,et al.  Computational free-surface fluid–structure interaction with application to floating offshore wind turbines , 2016 .

[62]  Eugenio Oñate,et al.  Surface tension problems solved with the Particle Finite Element Method using large time-steps , 2016 .

[63]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[64]  E. Shirani,et al.  Interface pressure calculation based on conservation of momentum for front capturing methods , 2005 .

[65]  Arif Masud,et al.  A variational multiscale stabilized formulation for the incompressible Navier–Stokes equations , 2009 .

[66]  A. Rollett,et al.  Critical instability at moving keyhole tip generates porosity in laser melting , 2020, Science.

[67]  Ng Niels Deen,et al.  Numerical simulation of gas bubbles behaviour using a three-dimensional volume of fluid method , 2005 .

[68]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[69]  Ming-Chen Hsu,et al.  An immersogeometric formulation for free-surface flows with application to marine engineering problems , 2020 .

[70]  Stephen Lin,et al.  A conservative level set method on unstructured meshes for modeling multiphase thermo-fluid flow in additive manufacturing processes , 2020 .

[71]  S. H. Lo,et al.  A new insertion sequence for incremental Delaunay triangulation , 2013 .

[72]  Liping Lei,et al.  A new ray tracing heat source model for mesoscale CFD simulation of selective laser melting (SLM) , 2020 .

[73]  Jinhui Yan,et al.  Variational multi-scale modeling of interfacial flows with a balanced-force surface tension model , 2020 .

[74]  Marcus Herrmann,et al.  A balanced force refined level set grid method for two-phase flows on unstructured flow solver grids , 2008, J. Comput. Phys..

[75]  Hanif Montazeri,et al.  A balanced-force algorithm for two-phase flows , 2014, J. Comput. Phys..

[76]  P. Masson,et al.  A new approach to compute multi-reflections of laser beam in a keyhole for heat transfer and fluid flow modelling in laser welding , 2013 .

[77]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[78]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[79]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[80]  R. Lahey,et al.  Computation of incompressible bubble dynamics with a stabilized finite element level set method , 2005 .

[81]  Carolin Körner,et al.  Evaporation model for beam based additive manufacturing using free surface lattice Boltzmann methods , 2014 .

[82]  Gerald L. Knapp,et al.  Mechanistic models for additive manufacturing of metallic components , 2021 .

[83]  T. Sun,et al.  Revealing transient powder-gas interaction in laser powder bed fusion process through multi-physics modeling and high-speed synchrotron x-ray imaging , 2020 .

[84]  T. Tezduyar,et al.  Parallel finite element computation of free-surface flows , 1999 .

[85]  James J. Feng,et al.  A diffuse-interface method for simulating two-phase flows of complex fluids , 2004, Journal of Fluid Mechanics.

[86]  Arif Masud,et al.  Residual-based turbulence models and arbitrary Lagrangian–Eulerian framework for free surface flows , 2015 .

[87]  J. Best The formation of toroidal bubbles upon the collapse of transient cavities , 1993, Journal of Fluid Mechanics.