A spatiotemporal two-level method for high-fidelity thermal analysis of laserpowder bed fusion

Numerical simulation of the laser powder bed fusion (LPBF) procedure for additive manufacturing (AM) is difficult due to the presence of multiple scales in both time and space, ranging from the part scale (order of millimeters/seconds) to the powder scale (order of microns/milliseconds). This difficulty is compounded by the fact that the regions of small-scale behavior are not fixed, but change in time as the geometry is produced. While much work in recent years has been focused on resolving the problem of multiple scales in space, there has been less work done on multiscale approaches for the temporal discretization of LPBF problems. In the present contribution, we extend on a previously introduced two-level method in space by combining it with a multiscale time integration method. The unique transfer of information through the transmission conditions allows for interaction between the space and time scales while reducing computational costs. At the same time, all of the advantages of the two-level method in space (namely its geometrical flexibility and the ease in which one may deploy structured, uniform meshes) remain intact. Adopting the proposed multiscale time integration scheme, we observe a computational speed-up by a factor ×2.44 compared to the same two-level approach with uniform time integration, when simulating a laser source traveling on a bare plate of nickel-based superalloy material following an alternating scan path of fifty laser tracks.

[1]  Silvia Bertoluzza,et al.  The Fat Boundary Method: Semi-Discrete Scheme and Some Numerical Experiments , 2005 .

[2]  Ernst Rank,et al.  A hierarchical computational model for moving thermal loads and phase changes with applications to selective laser melting , 2017, Comput. Math. Appl..

[3]  C. Kamath,et al.  Laser powder bed fusion additive manufacturing of metals; physics, computational, and materials challenges , 2015 .

[4]  Numerical solution of additive manufacturing problems using a two‐level method , 2021 .

[5]  T. Hughes,et al.  An improved implicit-explicit time integration method for structural dynamics , 1989 .

[6]  Frédéric Hecht,et al.  New development in freefem++ , 2012, J. Num. Math..

[7]  Li Ma,et al.  Single-Track Melt-Pool Measurements and Microstructures in Inconel 625 , 2018, 1802.05827.

[8]  B. Stucker,et al.  A Generalized Feed Forward Dynamic Adaptive Mesh Refinement and Derefinement Finite Element Framework for Metal Laser Sintering—Part I: Formulation and Algorithm Development , 2015 .

[9]  Silvia Bertoluzza,et al.  Analysis of the fully discrete fat boundary method , 2011, Numerische Mathematik.

[10]  P. Alam.  ‘W’ , 2021, Composites Engineering.

[11]  Lin Cheng,et al.  An optimally-coupled multi-time stepping method for transient heat conduction simulation for additive manufacturing , 2021 .

[12]  P. Michaleris,et al.  Numerical verification of an Octree mesh coarsening strategy for simulating additive manufacturing processes , 2019 .

[13]  R. Codina,et al.  An adaptive Finite Element strategy for the numerical simulation of additive manufacturing processes , 2020, Additive Manufacturing.

[14]  N. Hodge,et al.  Towards improved speed and accuracy of laser powder bed fusion simulations via representation of multiple time scales , 2020 .

[15]  Julia Mergheim,et al.  Thermal modelling of selective beam melting processes using heterogeneous time step sizes , 2019, Comput. Math. Appl..

[16]  K. Mills.  Recommended Values of Thermophysical Properties for Selected Commercial Alloys , 2001 .

[17]  Adrian Sandu,et al.  Multirate generalized additive Runge Kutta methods , 2016, Numerische Mathematik.

[18]  Adrian Sandu.  A Class of Multirate Infinitesimal GARK Methods , 2019, SIAM J. Numer. Anal..

[19]  Thomas J. R. Hughes,et al.  Implicit-Explicit Finite Elements in Transient Analysis: Stability Theory , 1978 .

[20]  C. Emmelmann,et al.  Additive Manufacturing of Metals , 2016 .

[21]  Bertrand Maury,et al.  A Fat Boundary Method for the Poisson Problem in a Domain with Holes , 2002, J. Sci. Comput..

[22]  Alessandro Reali,et al.  Suitably graded THB-spline refinement and coarsening: Towards an adaptive isogeometric analysis of additive manufacturing processes , 2018, Computer Methods in Applied Mechanics and Engineering.

[23]  Alessandro Reali,et al.  Accurate Prediction of Melt Pool Shapes in Laser Powder Bed Fusion by the Non-Linear Temperature Equation Including Phase Changes , 2019, Integrating Materials and Manufacturing Innovation.

[24]  Ted Belytschko,et al.  Mixed methods for time integration , 1979 .

[25]  Ted Belytschko,et al.  Partitioned rational Runge Kutta for parabolic systems , 1984 .

[26]  T. Belytschko,et al.  Stability of explicit‐implicit mesh partitions in time integration , 1978 .

[27]  Francois-Xavier Roux,et al.  Domain Decomposition Methodology with Robin Interface Matching Conditions for Solving Strongly Coupled Fluid-Structure Problems , 2009 .

[28]  Silvia Bertoluzza,et al.  A Fat boundary-type method for localized nonhomogeneous material problems , 2019, ArXiv.