Next Article in Journal
Perna canaliculus as an Ecological Material in the Removal of o-Cresol Pollutants from Soil
Previous Article in Journal
Prediction of Deflection of Reinforced Concrete Beams Considering Shear Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Laser Powder Bed Fusion (LPBF) of In718 and the Impact of Pre-Heating at 500 and 1000 °C: Operando Study

by
Asif Ur Rehman
1,2,3,*,
Fatih Pitir
1 and
Metin Uymaz Salamci
2,3,4
1
ERMAKSAN, Bursa 16065, Turkey
2
Department of Mechanical Engineering, Gazi University, Ankara 06570, Turkey
3
Additive Manufacturing Technologies Application and Research Center—EKTAM, Gazi University, Ankara 06560, Turkey
4
Manufacturing Technologies Center of Excellence—URTEMM A.S., Ankara 06980, Turkey
*
Author to whom correspondence should be addressed.
Materials 2021, 14(21), 6683; https://doi.org/10.3390/ma14216683
Submission received: 4 October 2021 / Revised: 26 October 2021 / Accepted: 1 November 2021 / Published: 5 November 2021

Abstract

:
The morphology of a melt pool has a critical role in laser powder bed fusion (LPBF). Nevertheless, directly characterizing the melt pool during LPBF is incredibly hard. Here, we present the melt pool flow of the entire melt pool in 3D using mesoscopic simulation models. The physical processes occurring within the melt pool are pinpointed. The flow patterns throughout the same are exposed and measured. Moreover, the impact of pre-heating at 500 and 1000 °C has been described. The study findings offer insights into LPBF. The findings presented here are critical for comprehending the LPBF and directing the establishment of improved metrics for process parameters optimization.

Graphical Abstract

1. Introduction

Additive manufacturing (AM) provides customized designs, reduced preparation time, and the ability to create complicated shapes. Many advanced technological applications [1] such as aerospace [2], biomedicine [3,4] and architecture [5,6] have considerable interest in it. Laser powder bed fusion (LPBF) is one of the most widely used additive manufacturing (AM) technologies because of its many advantages, which include significantly reduced structural limitations, high reproducibility, and on-time delivery [7]. In the LPBF, the metal particles are deposited layer upon layer via the blade or roller, followed by the fusing of particles by laser on particular locations to generate the desired slices, which is driven by CAD data [8]. There are a variety of flaws that are detrimental to efficiency and component quality [9,10,11] including balling, fractures, pores, and poor layer uniformity. As a result, it is necessary to better understand the deformations and the influence of input factors on the melt pool [12,13].
It has been shown that a variety of factors, including scan speed, laser power, particle sizes distribution (PSD), and layer height [14], have an impact on the melt pool and, therefore, the quality of the elements created [15,16]. Systematic attempts have been undertaken to explain the intricate melt pool dynamics [12,17,18,19,20], process parameters, and recurring defects in terms of the processing parameters that have an influence on the process. Studies by Hodge et al. and Lin et al. [21,22] looked at the effect of laser power and scan speed on the surface properties of LPBF parts. According to the research [21,22,23], irregularities, deformations, cracking, and other deficiencies on the surfaces are produced at a high scan rate, resulting in more surface defects. Studies have focused on the formation of the defect during the LPBF techniques of metallic metal powder [24]. The results of the investigation revealed that the energy density (ED) had a substantial influence on the formation of defects. The physics underpinning the dynamic interaction between the process parameters is, however, insufficiently explained by the available data. Concentrating only on the hit-and-trial strategy to extract the proper process variables throughout the LPBF experiments is both costly and time-consuming. Furthermore, the observational research [25,26] doubts the LPBF method’s ability to distinguish between the many complicated rules.
Melt stream inside the melt pool dictates heat hysteresis and mass flow in laser based 3D printing (e.g., powder bed fusion (PBF), laser metal deposition (LMD)) and therefore plays an essential part in melt pool formation [27,28], creation and progression of defects [29,30,31,32], solidification [33,34,35], and the formation of spatter [36,37,38]. As a consequence, the molten metal flow behavior has a significant impact on the surface morphology and mechanical characteristics of additively manufactured components [12,39,40,41,42]. As a result, monitoring transient characteristics during the melt is important for process optimization and also microstructure prediction models.
Given the difficulties of actually seeing the flow behavior within the melt pool, considerable scientific work has been conducted to investigate the flow dynamics under varying circumstances, i.e., with varying laser power or scan rates [43,44]. Experiments were also conducted to investigate the effect of a surface-active tracers on melt flow such as tungsten nanoparticles in LPBF of X-ray semitransparent Al based alloys. However, considering the intricacy of the underlying physics such as surface tension (surfactants itself can form high surface tension end), visualizing the associated multi-physics events is very difficult [27,45]. In fact, major visualization agents and parameters enhancements must be applied in most of the experiments [19,27,28,46,47]: certain implications are established for the purpose of visual convenience, some for the purpose of intricacy, whilst others due to the underpinning physics is not clear. Each of these inferences will have an effect on the reliability of the tests to some degree. As a consequence, despite the fact that the majority of the trials have been evaluated by evaluating the melt pool, the estimated melt flow outputs are often not accurate, and in certain cases are indeed contrary.
Flow prediction errors were discovered around the frontal depression [48,49], across the back [43,50], at the baseplate of the melt pool [51,52], over the depression’s exit [51,53], and the depression outflow [52]. Thus, comprehending the underlying melt flow characteristics inside the melt pool is essential in powder bed fusion [29,54]. The use of in-situ x-ray imaging has recently been shown for the investigation of laser additive manufacturing processes [28]. The latest research [27,55] utilized synchrotron x-ray scanning to track melt flow using tungsten particles, but only in single projection plane even though thousands of particles and tracers can be in the same plane. Given the complexity of the melt-pool in three dimensions the movement of the tracers, traced by the 2D imaging is put under question. The authors of this study recently developed the model for the flow in 3D [56,57].
Determining three-dimensional molten metal movement in PBF is complex, as the lasers scans faster, and the resultant melt pools are shallower. Owing to the difficulty, the full melt pool dynamics under LPBF settings has not really been exposed. Furthermore, part-scale thermal study shows interface heating continues to rise during deposition owing to heat buildup [58,59], and interfacial temperature fluctuation influences melt pool behavior. IR heating systems for layer pre-heating [60], scanner reheating as well as remelting methods [61,62,63], an additional unfocused laser beam [64], resistance heating of the base, [65,66] and inductive circuits on the baseplate [67], may lower the heat flow and is often used as an efficient thermal stresses reduction technique [68]. In LPBF, a substantial decrease in deformation could be seen for aluminum parts at a pre-heating degree of approximately 150 °C. Even with pre-heating level of 250 °C, distortions in the samples was no more detectable [69]. Kruth et al. [70] observed a 10 percent decrease in bending angles in the LPBF method by pre-heating up to 180 °C.
Part failure and deformation caused by residual stress may be avoided by pre-heating the chamber [65]. It was discovered that electron beam powder bed fusion (EPBF) part fabrication from pre-heated powders had an anisotropy in microstructure and mechanical properties [71]. Ali et al. showed significantly improved yielding and ductility in LPBF at 570 °C pre-heating temperatures [66]. Investigators also discovered that increasing the pre-heating temperatures allowed the breakdown of α′ martensitic in α + β stable microstructure. Recent research on the effect of pre-heating temperatures in metals printing has mostly focused on residual stress removal and deformation minimization, microstructure, and mechanical characteristics such as ductility as well as yield strength. Fully dense hard to weld materials such as alumina [72] and zirconia [64] were also printed using preheating of 800 and 1600 °C with LPBF respectively. Pre-heating can extend the range of LPBF from only weldable material processing to un-weldable materials such as ceramics.
Aviation and turbines parts frequently utilize INCONEL 718 (IN718) [73]. The slow precipitation-hardening mechanisms of IN718 make it highly weldable (ideal for LPBF) [73]. The microstructure influences the mechanical characteristics of IN718 [74]. LPBF produces a consistent microstructure with little porosity and finer dendritic grain [74,75]. The regular dendritic structure disappears and a needle-like δ phase precipitates at grain boundaries when γ′ and γ″ phases dissolve in the matrix when heat treatment is applied A needle-like δ phase precipitates along grain boundaries whenever the γ′ and γ″ phases breakdown in the matrix during heat treatment, replacing the normal dendritic structure. [74]. Heat-treated LPBF has comparable tensile as well as ductility to wrought IN718 [74].
Using laser additive manufacturing, we developed a system for detecting melt flow behavior over the whole melt pool, which we describe in this article. We uncovered the melt flow behavior of the whole melt pool and investigated the driving factors of liquid flow as well as the fundamental processes throughout the melt pool. Finally, the impact of pre-heating temperatures on melt pool formation is investigated numerically. We provide an explanation for the track’s variable melt pool size and shape. Furthermore, the flow rates is used to describe the flow caused by the Marangoni effect and the resulting matching flow.

2. Materials and Methods

2.1. Powder Bed Modeling

To simplify the process of computing the powder production and layering process, this will be divided into two stages: in the beginning, a multitude of particles are dropped directly on the surface to generate a powder stack; afterwards, the blade/recoater moves over the surface at a set speed and the particles advance into build chamber to create the layer.
The interaction technique using the nonlinear elastic equation is employed to assess the elastic real contact force [76], and the damping criterion is theoretically used to recognize the dissipation of mechanical energy [49,77,78].
It is at this point in the perpendicular planes when natural contact forces and damping forces in elastic materials occur. Relative stiffness across the plane is perpendicular, and the Young’s modulus as well as mass of the plane are both constant. In order to account for elastic contact force, no micro-slip method is used in the tangential route [76]. Table 1 shows the PSD for ERMAK-A241-IN718 supplied by ERMAKSAN, Bursa, Turkey, with D10, D50, and D90. The particle have been simulated utilizing the PSD that has been provided.
Throughout the whole study, the discrete element modeling (DEM) module from Flow Science, Santa Fe, NM, USA, has been employed to simulate the layer-by-layer deposition of SS316L stainless metal powder. A layer of powder was deposited utilizing discrete microparticles, rather than viewing the powder layer as a plate of uniform size.
The powder particles are seen clearly in Figure 1. According to the SEM image in Figure 1a, the real powder particle may be interpreted as round in a reasonable approximation. The particle sizes were known to be consistent with the experiment computation, and the particle sizes were found to be consistent with D10, D50, and D90 in the appropriate proportions (see Figure 1). Figure 1b depicts the particle that was created using the simulation.

2.2. Modeling of Powder Bed Deposition Process

Primarily based on the theoretical model mentioned above, the ERMAK-A241-IN718 The modeling of the metal powder creation process was carried out. Layer thickness is 50 microns due to the powder particles level. Deposition process is shown in detail in the second illustration. Simulation of DEM during powder bed deposition is shown in Figure 2a. Figure 2b shows a 50-micron layer that has been deposited. In the DEM model (see Figure 2a, a 100-micron layer was deposited. It was found that when the spaces on the outer/free surface were left in place, the layer of powder had a packing density of 65 percent; however, when the gaps on the outer/free surface were eliminated, this density rose to 90 percent.

2.3. Modeling of Thermophysical Properties

Temperature related physical properties of ERMAK-A241-IN718 have been simulated utilizing Sente Software, UK, depending on the chemical composition (Table 2) of ERMAK-A241-IN718 obtained from ERMAKSAN, Bursa, Turkey, as shown in Figure 3.
The super quick melting as well as solidification that occurs throughout the LPBF process has an implication on all of the thermo-physical characteristics involved. The cornerstone of the LPBF process modeling of the metal powder is a complete material characteristics simulation.
In order to simulate the thermo-physical behavior of the ERMAK-A241-In718, we used the chemical composition specified in Table 2 of this same document. The results are shown in Figure 3. As demonstrated in the graphs in Figure 3, melting as well as solidification may have an effect on the average expansion coefficient but also density [79]. Figure 3a shows the change in solid state to liquid state with temperature, Figure 3b shows the change in density and molar volume with temperature, and Figure 3c shows the average expansion coefficient and volume change with temperature, Figure 3d shows the electrical resistivity and thermal conductivity with temperature. Similarly, Figure 3e shows total viscosity and liquid diffusivity with temperature, and Figure 3f shows Poisson’s ratio and liquid viscosity with temperature.
Because of this, although the density decreases consistently with heating, the shift in averaged expansion coefficient is not totally consistent, which is one of the most important influencing elements during the treatment with the laser irradiation. The decrease in surface tension may also be shown to be not totally continuous, which results in a non-uniform Marangoni flow characteristic. The Marangoni flow develops in the fusion zone as a result of the difference in surface tension between the heated and cooled ends [36], which causes the flow to be unstable. It is also possible to illustrate that the change in Poisson’s ratio and Youngs moduli is completely non-consistent. While the latent heat increases with temperature in a constant manner, the increase in conductivity does not, which is one of the most important influencing elements for heat dissipation inside the building platform [79].

2.4. Numerical Model

A computational fluid dynamics (CFD) framework was developed and implemented by including particular subprocesses from the FLOW-3D 11.2v CFD package as well as weld module of Flow Science, Inc., United States of America (USA). Some equations have been explained below.
Multiple parameters and assumptions are considered in the research for simplicity: (1) the melting within the melt stream is assumed incompressible Newtonian; (2) the changes in mass attributable only to metal evaporation are often not included.
In Equations (1)–(3), the corresponding equations that can be solved for mass continuity, momentum conservation, as well as energy conservation are shown, respectively:
· v = 0
v t + ( v · ) v = 1 ρ P + μ 2 v + g [ 1 α ( T T m ) ] g [ 1 α ( T T m ) ]
h t + ( v · ) h = 1 ρ ( · k T )
where v denotes the velocity profile. P denotes pressure, μ denotes viscosity, and, g denotes gravity function, h denotes specific enthalpy, ρ specifies density, and k denotes heat conductivity Volume of fluid (VOF) models employ the free surface for data acquisition [80]. A simple Equation (4) may be used to explain the VOF approach.
V F t + ( ν · V F ) = 0
where V F denotes the volume fraction of metal present within the cell. V F = 1, implies that the cell is completely filled with fluid, while V F = 0 indicates that the cell is completely devoid of fluid. Counts in the middle demonstrate that there is free space on this cell’s entire surface.
Variability in melt pools may be caused by a variety of factors, including thermophysical properties, vapor suppression, as well as penetration. Moreover, because the Rosenthal approach is re-extracted out from heat equation thus removes evaporation, convection, and even the Marangoni influence, the heat equation may be written as [81,82], the corresponding term in Equation (5) for melt pool diameter extracted from Rosenthal formula [83] to explain the role played by thermo-physical characteristics in melting pool heterogeneity in heat transfer [81]:
ω = 8 π e · P η ρ C p V ( T m T 0 )
where ω is width of the melt pool, P defines beam power, and η specifies absorptivity, C p specifies heat capacity, V specifies scanning velocity, the melting temperature is specified by T m , and the preheating level is specified by T 0 . According to the assumption of thermally independent physical qualities as well as the thermophysical conductivity utilized to determine the size of the melt pool, the Rosenthal solution is found.
Additionally taken into consideration are the effects of recoil pressure, as well as vapor suppression, on the melt pool scale [84]. Equation (6) might be used to determine each individual recoil pressure:
P S = A · exp { B ( 1 T V T ) }
the primary coefficient A can be calculated by using the following: A = β P 0 , β ∈ [0.54, 0.56]. P 0 is the atmospheric pressure. The secondary coefficient B can be computed using the following equation: in which the B = Δ H v / R T v , Δ H v stands for the vaporization heat and the gas constant (R). T v stands for the saturation temperature [84].
It is well established that the energy density of the beam follows a Gaussian distribution. Throughout scanning, the laser moves at a constant scan speed, and the ED of the laser may be described mathematically [84] as the following Equation (7):
q = 2 A p π R b 2 e x p [ 2 ( x ν t x 0 ) 2 + ( y y 0 ) 2 R b 2 ]
where A defines the absorption of its laser beam by the powder particles, p represents laser power, Rb denotes the radius of the laser, and v indicates the pace of scanning. ( x 0 ,   y 0 ) is the coordinate system. The starting position of the laser beam center [84] is represented by this value. The beam radius has been represented by R b . Convection and radiation were resolved on this free surface, but evaporation could not be disregarded on the molten pool’s surface because of the presence of water. It follows as a consequence that the energy equation for the melt pool’s surface can be written as an Equation (8) [84], which is as follows:
T n = q h C ( T 1 T 0 1 ) σ 0 ε ( T 4 T 0 4 ) q e v a p
where, h c denotes the coefficient of convective heat transfer, T 0 is the ambient temperature, and the Stefan–Boltzmann constant is denoted by σ 0 . ε denotes the emissivity of a material, and q e v a p is the heat transfer owing to evaporation, and it may be expressed [84] by the Equation (9) which reads as follows:
q e v a p = ω 0 L v = exp ( 2.52 + 6.121 18 , 836 T 0.5 l o g T ) L v
where, ω 0 denotes the evaporation rate and, L v denotes the latent heat of evaporation. In this study, a new equation was constructed explicitly for mass flow rate, as specified and determined by Equation (10) in the findings and discussion section.
m ˙ = ρ · v d A

2.5. Configuration of the Modeling Environment and Variables

Primarily In718 powder generation and deposition process was mainly modeled utilizing the above-mentioned theoretical framework. The material properties of In718 for the LPBF simulations given below were selected from the Flow-3D available material dataset. The bed layer size was kept constant at 100 μm. The power (200 W), speed (3 m/s) and spot size has been kept constant. The preheating temperature has been changed from ambient temperature to 500 and 1000 °C in three cases. The domain size has been provided in the figure within the Supplementary Data.
Mobile workstation Precision 7530 form Dell Inc. Round Rock, TX, USA, with Windows 10 Pro for Workstations 64-bit was used for the simulations. It has a processor Intel(R) Xeon(R) E-2186M (CPU @ 2.90 GHz (12 CPUs), ~2.9 GHz) with a memory of 32,768 MB RAM and available OS Memory of 32,524 MB RAM. A simulation can take from 24 h to 1 week to complete.

2.6. Experimental Procedure

To validate the mathematical method, single-track LPBF tests were also carried out. The singular melting track was produced using an ERMAKSAN 250, which employs a fiber laser manufactured by ERMAKSAN EON Photonics, Bursa, Turkey. Atomized ERMAK-A241-In718 powder with an essentially spherical shape was utilized in the experiments. The beam’s diameter has been kept at 85 μm. The laser speed and power has been kept constant as described in the simulation.
Figure 4a–c shows a printed specimen. Figure 4d depicts the ENAVISION 250 LPBF system that was used for validation. The equipment breakdown is shown in Figure 4e. General specifications for ENAVISION 250 LPBF system have been provided in Table 3.
The OM micrographs were taken using an ZEISS smartzoom digital microscope from Carl Zeiss Microscopy Deutschland GmbH, Germany was used to capture optical micrographs.

3. Results and Discussion

Figure 5a–d shows the melt pool dynamic viscosity at 200, 270, 300, and 400 µs, respectively, and Figure 5e shows the experimental comparison through the SEM of the irradiated single track.
When the powder material is exposed to heat (laser) of high intensity and is heated far beyond the melting point and even to the boiling point. A recoiling pressure is applied to the melt pool due to vaporization surface, which results in a depression zone (as shown in region I). While heat transfers from the outside to the interior of the material below via the depression zone, it also melts the layer below as well. The melt flow characteristics of LPBF form quite differently, attributable to the multifaceted underlying physics which will be discussed in detail in the following sections. Due to the formation of the depression zone the melted metal is pushed backwards as shown in the region II. This is one of the influencing factors in the melt pool dynamics apart from the others which will be discussed in the following sections.
When the powder material that is exposed to heat (laser) of high intensity and is heated far beyond the melting point and even to the boiling point. A recoiling pressure is applied to the melt pool due to vaporization surface, which results in a depression zone. While heat transfers from the outside to the interior of the material below via the depression zone, it also melts the material’s inner area as well. The melt flow characteristics of LPBF form quite differently, attributable to the multifaceted underlying physics.
Figure 6a–f presents the state of the melt pool at 50, 140, 230, 320, 410, and 500 µs, respectively, with respect to the density of the material. Because the depression zone is constantly heated, the flow from higher temperatures fluid (with lower density) continues to flow. However, as the melted metal has a greater temperature (low density) than those at the top of the melt pool, it rises up owing to buoyancy throughout backward transmission and creates a vortex at the rear edge of the melt pool, as illustrated. To make things easier to explain, the density color gradient has been presented in Figure 6.
The surface tension difference is one of the another most influencing factor in the liquid state. When the surface tension difference is generated in the liquid between the two ends of the liquid, a very strong pull force is produced from the high surface tension end to the low surface tension end called the Marangoni force.
The flow streams will be discussed in detail in the following discussion.
While Figure 7a–f presents the flow stream traces of the melt pool flow at 50, 140, 230, 320, 410, and 500 µs, respectively. When the laser starts to irradiate the powder, each powder particle starts to melt and joins the melt pool formation, melting is highly dependent on the energy density and particle size distribution (PSD), all this happens in the matter of micro-seconds as shown in Figure 7a. The region I shows the irradiating particle being melted by laser while a Marangoni flow is generated from top (high temperature end) to bottom (low temperature end) as shown in region II of Figure 7a. As the laser moves a small amount of the melted pool forms a flow in the forward direction contributed by both the Marangoni force (with the powder particle in front of the laser)and the recoil pressure as shown in Figure 7b region I. Similarly, a large amount of fluid is pulled backwards due to the strong Marangoni force and the recoil pressure as shown in Figure 7b region II. At the rear end of the melt pool a circular flow or swirl is formed due to the high surface tension at the edges of the melt pool as well as the rear end.
As the flow progresses with the movement of the irradiating laser, the recoil pressure and the Marangoni force keep contributing to pushing the flow backwards. The surface tension difference triggered by the temperature difference (between the rear end and the laser irradiation end) pulls the liquid to the rear as shown in Figure 7c,d. As the single track progresses the flow contributed by the nonuniform cooling of the track can form distorted circular flows or multiple swirls, which contribute more to the uneven width and height of the irradiated single track.
The top view of the same single track can be seen in Figure 8a–f at 40, 130, 220, 310, 400, and 490 µs, respectively, where the color gradient is showing the density of the melt pool.
A 2D cross section of the track was obtained, and the amount of liquified metal going forward as well as backward was computed using mass flow rate calculations in order to fully comprehend the mass flow rate, as shown in Figure 9a. It depicts the flow dynamics of the melt pool prior and afterwards laser irradiation. The positive mass flow rate indicates that the flow is moving forward, whereas the negative mass flow rate indicates that the flow is moving backward. A small portion of the melt pool flows forward until the laser is irradiated, but when the laser passes over that area, the melt pool is dragged backwards owing to the difference in surface tension between the two regions. X-ray synchrotron imaging of melt pool in LPBF of Al based alloys using tungsten nanoparticle markers can be seen in Figure 9c, and schematics of melt pool from the mentioned experiments can be seen in Figure 9d.
The melt pool is propelled by five main factors. the Marangoni force, which is also known as Bernoulli’s principle, impels a fluid in one direction, from the elevated temperatures zone to the reduced temperature region for such a material having a negative coefficient of surface tension (all metals and alloys, with the exception of a few special ones, will experience a negative temperature coefficient of surface tension in the range from low surface tension area to high surface tension region) [54,85,86,87,88]. Recoil pressure is exerted toward the vaporized surface imposed as a result of an inward force perpendicular to the vaporized surface [89,90]. For very fast vapor plumes (which may reach speeds of 1102–103 m s−1 [91]), friction between the liquid and gas phases produces shear stress at the interfaces [86,92]. Hydraulic pressure is capable of exchanging energy through pressure (hydrostatic force) or momentum (hydrodynamic force) [93]. Liquid is propelled down a density gradient by buoyancy force [43,93,94].
The cross section in the middle of the melt pool in lateral direction has been taken with stream traces and can be seen in Figure 10a at 50, 120, 210, 300, 390, and 480 µs, respectively, where the color gradient is showing the density of the melt pool. When the laser starts to irradiate the powder particle start to form the melt pool and with the movement of the laser the pool is pushed backwards due to the recoil and Marangoni force as shown in Figure 10a where the stream traces show the path of the melt pool from its initial point in the powder layer. Similarly, when the laser keeps moving further the flow keeps being pulled from the higher surface tension (or relatively lower temperature region) as shown in Figure 10b, and a circular flow pattern starts to appear. The pull force toward the rear end and the circular flow pattern keep on increasing as shown in Figure 10c,d. As the whole single track starts to solidify and cool down, a hump is generated which contributes to the uneven height of the melt pool. The higher concentration of the melt pool causes the rear end to cool down slowly as compared to the region with the lower concentration as seen in Figure 10d. Now, at the end of the single-track cooling, when the melt pool concentrated region has higher temperature and the less concentrated has lower temperature, the liquid metal starts to flow in the opposite direction and this also provides evidence that the Marangoni is the dominating force.
Figure 11 shows the experimental comparison in lateral direction through the OM of the irradiated single track. It is in good agreement with the simulation results where the recoil pressure and Marangoni force is causing irregular height of the melt pool.
The cross section has been taken with velocity vectors of each point within the melt pool and can be seen in Figure 12a–f at 50, 120, 210, 300, 390, 480 µs, respectively.
The isometric view of the melt pool has been taken without the solid powder (transparent region) and can be seen in Figure 13a–f at 50, 120, 230, 320, 410, and 500 µs, respectively, where the color gradient in the transparent region is again showing the density of the melt pool. Here, the solidification of the melt pool can be elaboratively described. When the laser starts to irradiate a depression is formed as shown in Figure 13a and the recoil pressure and Marangoni pushes the melt pool backwards as shown in Figure 13b, and the surface tension furthers the pull as shown in Figure 13c. After the laser irradiation has been finished the dominating or the only force that is contributing to the flow is Marangoni force as shown in Figure 13d–f.
Because of the strong vaporization that causes the depression zone, the flow profile in a depression-mode melt pool is very complicated. The flow profile around depression [the area under laser] is dominated by vaporization-related pressures (i.e., recoil pressure and vapor plumes drag) and the Marangoni force. With Marangoni force as well as the vapor plumes friction, the flow across the front portion of the depression flows upwards [52,90,95], whereas the flow towards the depression’s base flows downward owing to recoil pressure, as shown in both transverse views as well as the crosswise view in Figure 7b,c and Figure 8a–c. The flow bypassing the depression zone is primarily driven by the hydraulic differential pressure, as a highly pressurized area forms in the fluid flow front of the depression zone but the front depression wall continues to “push” the front liquid throughout its path. It is worth noting that the upward flows may momentarily vanish or increase, as a result of the interaction between the downwards recoil pressures and the upward vapor plumes friction, as well as the Marangoni forces at front depression wall [86,90,92] as shown in Figure 7b,c, Figure 8a–c and Figure 10a–c.
As the depression region progresses, the low hydraulic-pressure area develops behind the depression’s base [96], It causes drawing the adjoining fluid into the rear side of the depression zone. As for the melt, its surface flows to rear end from an elevated temperature area to a relatively low temperature (Marangoni force), emulating huge mass transfer. The stream going backwards picks up momentum as it travels backwards.
The flow keeps on flowing from elevated temperatures fluid (with lower density) because it is being continuously heated at the depression zone. However, the melted metal rises up due to buoyancy throughout backward transmission because it has higher temperature (low density) than those at the top of the melt pool and it forms a vortex at the rear side of the melt pool as shown in Figure 7b–d, Figure 10b–d and Figure 12b–d.
Apart from an upward motion, the stream in the medium level of the melt pool has a strong propensity to travel forward, drawn by the low hydraulic pressure area surrounding the depression zone caused by the high-speed flows under Bernoulli’s effect. Backward flow on the surface is produced by Marangoni forces until the “rear end” area of the melt pool, whereas the flow from the rear end starts to flow forward from the base of the rear end due to formation of vortex. When the flow comes back just at melt pool middle, it collides with the flow coming from the depression zone, momentum is transferred forming a joint flow upward and later forming two vortexes. During the collision a rear end vortex and the frontal vortex is formed.
The flows from transverse view create two more closed loops, which are likewise governed by Marangoni force (consistent with recent research [54]), following a similar process to that of the longitudinal view. The flow produces a clockwise vortex 1 front of the beam as well as a counterclockwise vortex 2 after the laser beam has irradiated, as shown in the images. The frames per micro-second in video format for (1) cross sectional stream traces at RT, (2) cross sections velocity vectors at RT, (3) velocity vectors isometric view at RT, and (4) velocity vectors top view at RT, can be obtained from Supplementary Data.
Pre-Heating at 500 °C
The state of the melt pool in Figure 14a–f at 50, 140, 230, 320 410 and 500 µs, respectively, with 500 °C pre-heating has been shown, similar to the Figure 6 at the room temperature. The difference from the one at RT can be clearly visualized in terms of melt pool and solidification in this figure while the change in flow will be discussed in the following sections. Similar to the melt pool irradiation at the room temperature, the density decreases rapidly with melting due to the change in state which in turn increases the volume of the fluid slightly. However, as the temperature of the surrounding is kept at 500 °C the solidification is much slower, even after 500 µs there is still the melt pool remaining.
Figure 15a–f presents the flow stream traces of the melt pool flow at 50, 140, 230, 320, 410, 500 µs, respectively, with 500 °C pre-heating. When the laser starts to irradiate the powder, the particles start to melt more quickly and joins the melt pool formation quicker than in RT. The region I shows the irradiating particle being melted by laser while a Marangoni flow is generated from top (high temperature end) to bottom (low temperature end) as shown in region II of Figure 15a which is stronger than the one RT. As the laser moves a small amount of the melt pool is formed in the forward direction as shown in Figure 15b region I, which is more stable when compared with the one at RT. Similarly, a large amount of fluid is pulled backwards as shown in Figure 15b region II and a circular flow or swirl is formed due to the high surface tension at the edges of the melt pool as well as the rear end. The melt pool becomes wider as compared with the one RT.
As the flow progresses with the movement of the irradiating laser, the recoil pressure and the Marangoni force keep contributing to pushing the flow backwards as shown in Figure 15c,d. However, in Figure 15e,f, a positive aspect of pre-heating can be visualized when compared with Figure 6 at RT, that there has not been the distortion in the flow pattern because of the slow solidification.
The top view of the same single track can be seen in Figure 16a–f at 40, 130, 220, 310, 400, and 490 µs, respectively, where the color gradient is showing the density of the melt pool.
A 2D cross section was taken in the single track with 500 °C pre-heating, the melt pool flowing forward and backward was calculated in Figure 17a through mass flow rate calculation to understand the mass flow rate like the one in Figure 9a. In this case, the laser irradiation at 500 °C pre-heating causes the melt pool to flow differently than it did before. The positive flowrate indicates forward movement, whereas the negative indicates backward flow. After being exposed to laser radiation, some of the melt pool flows forward; however, when the laser passes across that location, some of the melt pool is dragged backwards. However, it keeps moving forward after the laser irradiation as shown in the graph, even though it also keeps moving backwards. The cross section through which the mass flow has been described can be seen in the Figure 17b,c.
The cross section in the middle of the melt pool with 500 °C pre-heating has been shown in Figure 18a–f at 50, 120, 210, 300, 390, and 480 µs, respectively, where the color gradient is showing the density of the melt pool. Due to the chamber heating the density of the powder particle is slightly decreased. When the laser starts to irradiate the powder the melt pool is formed more rapidly as compared with the one at RT as shown in Figure 18a, where the stream traces show the path of the melt pool from its initial point in the powder layer. Similarly, when the laser keeps moving further the flow keeps being pulled from the higher surface tension as shown in Figure 18b and a circular flow pattern starts at the rear end. The melt pool is deeper due to the pre heating and the circular flow pattern is wider as shown in Figure 18c,d. The higher concentration of the melt pool causes the rear end to cool down slowly as compared to the region with the lower concentration, the liquid metal starts to flow in the opposite as shown in Figure 18e,f. Due to the pre-heating the melt pool has more time to settle down, as the flow moves backwards in Figure 18e,f and the height of the hump is decreased.
The cross section with 500 °C pre-heating has been taken with velocity vectors of each point within the melt pool and can be seen in Figure 19a–f at 50, 120, 210, 300, 390, 480 µs, respectively. Similar to the one discusses above, when the laser starts to irradiate 500 °C pre-heated chamber the powder particle starts to form the melt pool more rapidly as shown in Figure 19a. A reaction force in the deeper depression zone is generated on the laser contacting surface as shown by the velocity vectors in the top surface in Figure 19a similar to the one in RT. However, as the laser progresses, due to the surrounding air the upper region of the melt pool which is in direct contact with the air starts to cool (can be visualized with the velocity vectors on the tope surface caused by Marangoni force as well as by the density color gradient) with huge increase in the surface tension, in Figure 19b–d.
A circular flow pattern can also be visualized with the velocity vectors at the rear end of the melt pool in Figure 19b–f. However, as the front-end temperature decreases with time and the surface tension over there increase, it starts to pull the liquid in the opposite direction as seen in Figure 19e,f from the bottom of circular flow pattern.
The isometric view of the melt pool with 500 °C pre-heating has been taken without the solid powder (transparent region) and can be seen in Figure 20a–f at 50, 120, 230, 320, 410, and 500 µs, respectively, similar to the one at 500 °C pre-heating in Figure 13. Here, the melt pool flow can be visualized which is wider and deeper with the same energy density. When the laser starts to irradiate a depression is formed as shown in Figure 20a and the recoil pressure and Marangoni pushes the melt pool backwards as shown in Figure 20b, and the surface tension furthers the pull as shown in Figure 20c. After the laser irradiation has been finished the dominating or the only force that is contributing to the flow is Marangoni force as shown in Figure 20d–f. The frames per micro-second in video format for (1) cross sectional stream traces at 500 °C, (2) cross sections velocity vectors at 500 °C, (3) velocity vectors isometric view at 500 °C, and (4) velocity vectors top view at 500 °C, can be obtained from Supplementary Data.
Pre-Heating at 1000 °C
The state of the melt pool in Figure 21a–f at 50, 140, 230, 320, 410 and 500 µs, respectively, with 1000 °C pre-heating, similar to the Figure 14 and Figure 6 at the 500 °C pre-heating and RT. The difference from the previous ones can be clearly visualized in terms of melt pool size and in this figure. The width and depth of the melt pool has been significantly increases and solidification at the end of the 500 µs has not completed and much of the melt pool in in liquid state.
Figure 22a–f presents the flow stream traces of the melt pool flow at 50, 140, 230, 320, 410, 500 µs with 1000 °C pre-heating. As the powder bed is preheated at 1000 °C, the laser needs little energy to melt, the particles form the melt pool more rapidly. The region I shows the irradiating particle being melted by laser while a Marangoni flow is generated from top (high temperature end) to bottom (low temperature end) as shown in region II of Figure 22a which is stronger than the one at 500 °C pre-heating. As the laser moves a small amount of the melt pool is formed in the forward direction as shown in Figure 22b region I (it is more significant from 500 °C pre-heating as well as RT). Similarly, a large amount of fluid is pulled backwards as shown in Figure 22b region II and a circular flow or swirl is formed due to the high surface tension at the edges of the melt pool as well as the rear end. The melt pool becomes wider as compared with the one 500 °C pre-heating.
As the flow progresses with the movement of the irradiating laser, the recoil pressure and the Marangoni force keep contributing to pushing the flow backwards as shown in Figure 22c,d. However, in Figure 22e,f, a positive aspect of pre-heating can be visualized when compared with Figure 15 at RT, that there has not been the distortion in the flow pattern because of the slow solidification. Secondly, another positive aspect at 1000 °C pre-heating is that the melt pool has time to flow back towards the front end, in Figure 22e,f, when the front-end temperature drops, and the rear end temperature is high due to the accumulated melt pool.
The cross section in the middle of the melt pool with 1000 °C pre-heating has been shown in Figure 23a–f at 50, 120, 210, 300, 390, and 480 µs. Due to the pre heating the density of the powder particle and the baseplate is decreased as seen by the color gradient in comparison to previous. When the laser starts to irradiate the powder melt pool is formed easily and significantly bigger as shown in Figure 23a, where the stream traces show the path of the melt pool from its initial point in the powder layer. Similarly, when the laser keeps moving a recoil pressure reaction force pushes the liquid backwards and upward simultaneously as shown in ‘region I’ in Figure 23b, the flow keeps being pulled due to Marangoni, and a circular flow pattern starts at the rear end which is more significant in comparison. The melt pool is deeper due to the pre-heating and the circular flow pattern is wider as shown in Figure 23c,d. However, the hump is less significant when compared with the previous ones, because the surface tension difference around the melt pool is low consequent of the pre-heating. The higher concentration of the melt pool causes the rear end to cool down slowly as compared to the region with the lower concentration, the liquid metal starts to flow in the opposite as shown in Figure 23e,f. With 1000 °C pre-heating the solidification is very slow which causes the melt pool more time to stabilize. Furthermore, as the flow moves backwards in Figure 23e,f and is still moving at 500 µs.
At one point, the single track had a 2D cross section taken with 1000 °C pre-heating, the melt pool flowing forward and backward was calculated in Figure 24a through mass flow. It depicts the melt pool flow’s pattern before and after laser irradiation at 1000 °C pre-heating and compares the two. After being exposed to laser radiation, some of the melt pool flows ahead; however, when the laser passes across that location, some of the melt pool is dragged backwards. However, it keeps moving forward after the laser irradiation as shown in the graph, even though it also keeps moving backwards which is more significant than both previous ones. The cross section through which the mass flow has been described can be seen in Figure 24b,c.
The cross section with 1000 °C pre-heating has been taken with velocity vectors at each point within the melt pool and can be seen in Figure 25a–f at 50, 120, 210, 300, 390, and 480 µs, respectively. As in the example discussed above, when the laser starts to irradiate at 1000 °C the pre-heated powder particles start to form the melt pool more rapidly as shown in Figure 25a. A reaction force in the deeper depression zone is generated on the laser contacting surface as shown by the velocity vectors in the top surface in Figure 25a like the one at 500 °C pre-heating. However, as the laser progresses, due to the surrounding air the upper region of the melt pool which is in direct contact with the air starts to cool (can be visualized with the velocity vectors on the tope surface caused by Marangoni force as well as by the density color gradient) with huge increase in the surface tension, in Figure 25b–d.
A circular flow pattern can also be visualized with the velocity vectors at the rear end of the melt pool in Figure 25b–f. However, as the front-end temperature decreases with time and the surface tension over there increase, it starts to pull the liquid in the opposite direction as seen in Figure 25e,f from the bottom of the circular flow pattern.
The isometric view of the melt pool with 1000 °C pre-heating has been taken without the solid powder (transparent region) and can be seen in Figure 26 at 50, 120, 230, 320, 410, and 500 µs, respectively, similar to the one at 500 °C pre-heating in Figure 13. Here, the melt pool flow can be visualized which is wider and deeper with the same energy density. When the laser starts to irradiate a depression is formed as shown in Figure 26a and the recoil pressure and Marangoni pushes the melt pool backwards as shown in Figure 26b, and the surface tension furthers the pull as shown in Figure 26c. After the laser irradiation has been finished the dominating or the only force that is contributing to the flow is Marangoni force as shown in Figure 26d–f. The frames per micro-second in video format for (1) cross sectional stream traces at 1000 °C, (2) cross sections velocity vectors at 1000 °C, (3) velocity vectors isometric view at 1000 °C, and (4) velocity vectors top view at 1000 °C, can be obtained from Supplementary Data.
The mass flow pushed forward and pulled backward at RT, 500 °C and 1000 °C pre-heating has been shown Figure 27. Here, the change in the mass flow can be understood clearly. As we pre-heat the mass flow has more time to stabilize which leads to positive mass flow even after the laser irradiation. However, at the room temperature the forward mass flow is stopped after the irradiation.

4. Conclusions

We have exposed and measured the melt flow dynamics throughout the whole melt pool using velocity vectors at different points, stream traces of the fluid movement, and mass flow rate forward and backward movement. in the LPBF for the first time in this study. The main conclusions are discussed in depth below:
  • The melted flow patterns in every region of the entire melt pool have shown a fairly complicated flow pattern due to the interaction of several driving factors. We measured and derived equation a drop in flow velocity from the depression region to the melt pool rear region.
  • The propelling factors for various melt flow variations were studied. The Marangoni effect is responsible for the circulation flow from a low surface tension to a high surface tension region on the melted surface. The fluid movement all across depression-zone sidewalls is dominated by evaporation. Hydraulic pressure propels fluid movement from an area of high pressures to a low-pressure area. Buoyance pressure is responsible for fluid convection between low- and high-density regions.
  • Only at laser irradiating region of its melt pool, the powerful major factors are recoil liquid momentum and heat convection; just at rear of the melt pool, the dominant major factors have been determined for being high surface tension and the thermal conductivity, respectively.
  • Throughout this study, modeling shows an increase in pre-heating level across the melt pool variability (wide, height, and duration).
  • greater pre-heating temperature create a melt pool having greater depth and relatively dimension. In modeling, a front sidewall inclination of irradiating region increases with pre-heating temperatures, indicating more laser drill force. Increasing the temperature degree increases penetration.
  • As the solidification speed rises, the melting period tends to become shorter. Due to increased retribution force and fluid suppressing, the reversed melt flow from both the laser center region and Marangoni power is not possible leading to higher residual stresses. Lower pre levels typically have faster movement. At higher pre-heating levels the residual stresses can be reduced as the melt pool has more time to release the pressures.
A much more advanced in-situ monitoring system would aid in achieving more accurate heat transmission within the melting zone. The research contemplates even more comprehensive visual investigation into the development of melt pools associated variability changes using precise multi-physical simulations. This will help with the processing parameters guidelines, which are designed to prevent defects or greatly improve the productivity of the items manufactured by PBF processes. Future work may address other important processing parameter studies with respect to scanning speed and laser power with the Marangoni effects.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14216683/s1, Supplementary Data: Temperature Change, Supplementary Data: Density Change.

Author Contributions

A.U.R., conceptualization, methodology, software, validation, formal analysis, investigation, writing—original draft, writing—review and editing, project administration, funding acquisition; F.P., conceptualization, methodology, validation, formal analysis, resources, writing—review and editing, supervision, project administration, funding acquisition; M.U.S., conceptualization, methodology, validation, formal analysis, resources, writing—review and editing, supervision, project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received financial support from the European Union’s Horizon 2020 (H2020) research and innovation program under the Marie Skłodowska-Curie grant, agreement No. 764935.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Murat Bahsi, Ahmet Sever, Sibel Batu, Mücteba Burak Karakaş. Berkan Başaran, and Volkan Mesut Küçükakarsu from ERMAKSAN for their valuable inputs during this research. The kind assistance of Anil Ummet Yanalak for OM is also acknowledged. The support from the Turkish Aerospace for the CFD post-processing is also duly acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kok, Y.; Tan, X.P.; Wang, P.; Nai, M.L.S.; Loh, N.H.; Liu, E.; Tor, S.B. Anisotropy and heterogeneity of microstructure and mechanical properties in metal additive manufacturing: A critical review. Mater. Des. 2018, 139, 565–586. [Google Scholar] [CrossRef]
  2. Guo, N.; Leu, M.C. Additive manufacturing: Technology, applications and research needs. Front. Mech. Eng. 2013, 8, 215–243. [Google Scholar] [CrossRef]
  3. Tiwari, S.K.; Pande, S.; Agrawal, S.; Bobade, S.M. Selection of selective laser sintering materials for different applications. Rapid Prototyp. J. 2015, 21, 630–648. [Google Scholar] [CrossRef]
  4. Liu, F.H. Synthesis of bioceramic scaffolds for bone tissue engineering by rapid prototyping technique. J. Sol-Gel Sci. Technol. 2012, 64, 704–710. [Google Scholar] [CrossRef]
  5. Ur Rehman, A.; Sglavo, V.M. 3D printing of geopolymer-based concrete for building applications. Rapid Prototyp. J. 2020, 26, 1783–1788. [Google Scholar] [CrossRef]
  6. Ur Rehman, A.; Sglavo, V.M. 3D printing of Portland cement-containing bodies. Rapid Prototyp. J. 2021. ahead-of-print. [Google Scholar] [CrossRef]
  7. Mahmood, M.A.; Popescu, A.C.; Mihailescu, I.N. Metal matrix composites synthesized by laser-melting deposition: A review. Materials 2020, 13, 2593. [Google Scholar] [CrossRef]
  8. Mahmood, M.A.; Popescu, A.C.; Oane, M.; Channa, A.; Mihai, S.; Ristoscu, C.; Mihailescu, I.N. Bridging the analytical and artificial neural network models for keyhole formation with experimental verification in laser melting deposition: A novel approach. Results Phys. 2021, 26, 104440. [Google Scholar] [CrossRef]
  9. Körner, C. Additive manufacturing of metallic components by selective electron beam melting—A review. Int. Mater. Rev. 2016, 61, 361–377. [Google Scholar]
  10. Fayazfar, H.; Salarian, M.; Rogalsky, A.; Sarker, D.; Russo, P.; Paserin, V.; Toyserkani, E. A critical review of powder-based additive manufacturing of ferrous alloys: Process parameters, microstructure and mechanical properties. Mater. Des. 2018, 144, 98–128. [Google Scholar] [CrossRef]
  11. Ur Rehman, A. Additive Manufacturing of Ceramic Materials and Combinations with New Laser Strategies. Master’s Thesis, Nanjing University of Science and Technology, Nanjing, Jiangsu, China, 2017. [Google Scholar]
  12. Everton, S.K.; Hirsch, M.; Stavroulakis, P.I.; Leach, R.K.; Clare, A.T. Review of in-situ process monitoring and in-situ metrology for metal additive manufacturing. Mater. Des. 2016, 95, 431–445. [Google Scholar] [CrossRef]
  13. Sing, S.L.; An, J.; Yeong, W.Y.; Wiria, F.E. Laser and electron-beam powder-bed additive manufacturing of metallic implants: A review on processes, materials and designs. J. Orthop. Res. 2016, 34, 369–385. [Google Scholar] [CrossRef]
  14. Olakanmi, E.O.; Cochrane, R.F.; Dalgarno, K.W. A review on selective laser sintering/melting (SLS/SLM) of aluminium alloy powders: Processing, microstructure, and properties. Prog. Mater. Sci. 2015, 74, 401–477. [Google Scholar] [CrossRef]
  15. Ansari, P.; Rehman, A.U.; Pitir, F.; Veziroglu, S.; Mishra, Y.K.; Aktas, O.C.; Salamci, M.U. Selective Laser Melting of 316L Austenitic Stainless Steel: Detailed Process Understanding Using Multiphysics Simulation and Experimentation. Metals 2021, 11, 1076. [Google Scholar] [CrossRef]
  16. Ullah, A.; Wu, H.; Ur Rehman, A.; Zhu, Y.; Liu, T.; Zhang, K. Influence of laser parameters and Ti content on the surface morphology of L-PBF fabricated Titania. Rapid Prototyp. J. 2021, 27, 71–80. [Google Scholar] [CrossRef]
  17. Bidare, P.; Maier, R.R.J.; Beck, R.J.; Shephard, J.D.; Moore, A.J. An open-architecture metal powder bed fusion system for in-situ process measurements. Addit. Manuf. 2017, 16, 177–185. [Google Scholar] [CrossRef]
  18. Yamada, T.; Shobu, T.; Nishimura, A.; Yonemoto, Y.; Yamashita, S.; Muramatsu, T. In-situ X-ray observation of molten pool depth during laser micro welding. J. Laser Micro Nanoeng. 2012, 7, 244–248. [Google Scholar] [CrossRef] [Green Version]
  19. Scipioni Bertoli, U.; Guss, G.; Wu, S.; Matthews, M.J.; Schoenung, J.M. In-situ characterization of laser-powder interaction and cooling rates through high-speed imaging of powder bed fusion additive manufacturing. Mater. Des. 2017, 135, 385–396. [Google Scholar] [CrossRef]
  20. Repossini, G.; Laguzza, V.; Grasso, M.; Colosimo, B.M. On the use of spatter signature for in-situ monitoring of Laser Powder Bed Fusion. Addit. Manuf. 2017, 16, 35–48. [Google Scholar] [CrossRef]
  21. King, W.; Anderson, A.T.; Ferencz, R.M.; Hodge, N.E.; Kamath, C.; Khairallah, S.A. Overview of modelling and simulation of metal powder bed fusion process at Lawrence Livermore National Laboratory. Mater. Sci. Technol. 2015, 31, 957–968. [Google Scholar] [CrossRef]
  22. Le, T.-N.; Lo, Y.-L.; Lin, Z.-H. Numerical simulation and experimental validation of melting and solidification process in selective laser melting of IN718 alloy. Addit. Manuf. 2020, 36, 101519. [Google Scholar] [CrossRef]
  23. Rahman, M.S.; Schilling, P.J.; Herrington, P.D.; Chakravarty, U.K. A Comparison of the Thermo-Fluid Properties of Ti-6Al-4V Melt Pools Formed by Laser and Electron-Beam Powder-Bed Fusion Processes. J. Eng. Mater. Technol. 2021, 143. [Google Scholar] [CrossRef]
  24. Gong, H.; Rafi, K.; Gu, H.; Starr, T.; Stucker, B. Analysis of defect generation in Ti-6Al-4V parts made using powder bed fusion additive manufacturing processes. Addit. Manuf. 2014, 1, 87–98. [Google Scholar] [CrossRef]
  25. DebRoy, T.; Wei, H.L.; Zuback, J.S.; Mukherjee, T.; Elmer, J.W.; Milewski, J.O.; Beese, A.M.; Wilson-Heid, A.; De, A.; Zhang, W. Additive manufacturing of metallic components—Process, structure and properties. Prog. Mater. Sci. 2018, 92, 112–224. [Google Scholar]
  26. Frazier, W.E. Metal additive manufacturing: A review. J. Mater. Eng. Perform. 2014, 23, 1917–1928. [Google Scholar]
  27. Guo, Q.; Zhao, C.; Qu, M.; Xiong, L.; Hojjatzadeh, S.M.H.; Escano, L.I.; Parab, N.D.; Fezzaa, K.; Sun, T.; Chen, L. In-situ full-field mapping of melt flow dynamics in laser metal additive manufacturing. Addit. Manuf. 2020, 31, 100939. [Google Scholar] [CrossRef]
  28. Guo, Q.; Zhao, C.; Qu, M.; Xiong, L.; Escano, L.I.; Hojjatzadeh, S.M.H.; Parab, N.D.; Fezzaa, K.; Everhart, W.; Sun, T.; et al. In-situ characterization and quantification of melt pool variation under constant input energy density in laser powder bed fusion additive manufacturing process. Addit. Manuf. 2019, 28, 600–609. [Google Scholar] [CrossRef]
  29. Leung, C.L.A.; Marussi, S.; Atwood, R.C.; Towrie, M.; Withers, P.J.; Lee, P.D. In situ X-ray imaging of defect and molten pool dynamics in laser additive manufacturing. Nat. Commun. 2018, 9, 1355. [Google Scholar] [CrossRef] [Green Version]
  30. Kawahito, Y.; Mizutani, M.; Katayama, S. Elucidation of high-power fibre laser welding phenomena of stainless steel and effect of factors on weld geometry. J. Phys. D Appl. Phys. 2007, 40, 5854–5859. [Google Scholar] [CrossRef]
  31. Kawahito, Y.; Uemura, Y.; Doi, Y.; Mizutani, M.; Nishimoto, K.; Kawakami, H.; Tanaka, M.; Fujii, H.; Nakata, K.; Katayama, S. Elucidation of the effect of welding speed on melt flows in high-brightness and high-power laser welding of stainless steel on basis of three-dimensional X-ray transmission in situ observation. Weld. Int. 2017, 31, 206–213. [Google Scholar] [CrossRef]
  32. Olakanmi, E.O. Selective laser sintering/melting (SLS/SLM) of pure Al, Al-Mg, and Al-Si powders: Effect of processing conditions and powder properties. J. Mater. Process. Technol. 2013, 213, 1387–1405. [Google Scholar] [CrossRef]
  33. Yadroitsev, I.; Krakhmalev, P.; Yadroitsava, I.; Johansson, S.; Smurov, I. Energy input effect on morphology and microstructure of selective laser melting single track from metallic powder. J. Mater. Process. Technol. 2013, 213, 606–613. [Google Scholar] [CrossRef]
  34. Amato, K.N.; Gaytan, S.M.; Murr, L.E.; Martinez, E.; Shindo, P.W.; Hernandez, J.; Collins, S.; Medina, F. Microstructures and mechanical behavior of Inconel 718 fabricated by selective laser melting. Acta Mater. 2012, 60, 2229–2239. [Google Scholar] [CrossRef]
  35. Ansari, P.; Salamci, M.U. On the selective laser melting based additive manufacturing of AlSi10Mg: The process parameter investigation through multiphysics simulation and experimental validation. J. Alloys Compd. 2022, 890, 161873. [Google Scholar] [CrossRef]
  36. Nasab, M.H.; Gastaldi, D.; Lecis, N.F.; Vedani, M. On morphological surface features of the parts printed by selective laser melting (SLM). Addit. Manuf. 2018, 24, 373–377. [Google Scholar] [CrossRef]
  37. Simonelli, M.; Tuck, C.; Aboulkhair, N.T.; Maskery, I.; Ashcroft, I.; Wildman, R.D.; Hague, R. A Study on the Laser Spatter and the Oxidation Reactions During Selective Laser Melting of 316L Stainless Steel, Al-Si10-Mg, and Ti-6Al-4V. Metall. Mater. Trans. A Phys. Metall. Mater. Sci. 2015, 46, 3842–3851. [Google Scholar] [CrossRef]
  38. Liu, Y.; Yang, Y.; Mai, S.; Wang, D.; Song, C. Investigation into spatter behavior during selective laser melting of AISI 316L stainless steel powder. Mater. Des. 2015, 87, 797–806. [Google Scholar] [CrossRef]
  39. Herzog, D.; Seyda, V.; Wycisk, E.; Emmelmann, C. Additive manufacturing of metals. Acta Mater. 2016, 117, 371–392. [Google Scholar] [CrossRef]
  40. Hebert, R.J. Viewpoint: Metallurgical aspects of powder bed metal additive manufacturing. J. Mater. Sci. 2016, 51, 1165–1175. [Google Scholar] [CrossRef]
  41. Zhou, Q.; Hayat, M.D.; Chen, G.; Cai, S.; Qu, X.; Tang, H.; Cao, P. Selective electron beam melting of NiTi: Microstructure, phase transformation and mechanical properties. Mater. Sci. Eng. A 2019, 744, 290–298. [Google Scholar] [CrossRef]
  42. Kamath, C.; El-Dasher, B.; Gallegos, G.F.; King, W.E.; Sisto, A. Density of additively-manufactured, 316L SS parts using laser powder-bed fusion at powers up to 400 W. Int. J. Adv. Manuf. Technol. 2014, 74, 65–78. [Google Scholar] [CrossRef] [Green Version]
  43. Sohail, M.; Han, S.W.; Na, S.J.; Gumenyuk, A.; Rethmeier, M. Characteristics of weld pool behavior in laser welding with various power inputs. Weld. World 2014, 58, 269–277. [Google Scholar] [CrossRef]
  44. Matsunawa, A.; Kim, J.-D.; Seto, N.; Mizutani, M.; Katayama, S. Dynamics of keyhole and molten pool in laser welding. J. Laser Appl. 1998, 10, 247–254. [Google Scholar] [CrossRef]
  45. Rauniyar, S.K.; Chou, K. Melt Pool Analysis and Mesoscale Simulation of Laser Powder Bed Fusion Process (L-PBF) with Ti-6Al-4V Powder Particles. Jom 2019, 71, 938–945. [Google Scholar] [CrossRef]
  46. Shrestha, S.; Rauniyar, S.; Chou, K. Thermo-Fluid Modeling of Selective Laser Melting: Single-Track Formation Incorporating Metallic Powder. J. Mater. Eng. Perform. 2019, 28, 611–619. [Google Scholar] [CrossRef]
  47. Wolff, S.J.; Wu, H.; Parab, N.; Zhao, C.; Ehmann, K.F.; Sun, T.; Cao, J. In-situ high-speed X-ray imaging of piezo-driven directed energy deposition additive manufacturing. Sci. Rep. 2019, 9, 962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Gan, Z.; Liu, H.; Li, S.; He, X.; Yu, G. Modeling of thermal behavior and mass transport in multi-layer laser additive manufacturing of Ni-based alloy on cast iron. Int. J. Heat Mass Transf. 2017, 111, 709–722. [Google Scholar] [CrossRef] [Green Version]
  49. Tian, Y.; Yang, L.; Zhao, D.; Huang, Y.; Pan, J. Numerical analysis of powder bed generation and single track forming for selective laser melting of SS316L stainless steel. J. Manuf. Process. 2020, 58, 964–974. [Google Scholar] [CrossRef]
  50. Khairallah, S.A.; Anderson, A.T.; Rubenchik, A.M.; King, W.E. Laser powder-bed fusion additive manufacturing: Physics of complex melt flow and formation mechanisms of pores, spatter, and denudation zones. In Additive Manufacturing Handbook Product Development for the Defense Industry; CRC Press: Boca Raton, FL, USA, 2017; pp. 613–628. [Google Scholar] [CrossRef]
  51. Khairallah, S.A.; Anderson, A.T.; Rubenchik, A.M.; King, W.E. Laser powder-bed fusion additive manufacturing: Physics of complex melt flow and formation mechanisms of pores, spatter, and denudation zones. Acta Mater. 2016, 108, 36–45. [Google Scholar] [CrossRef] [Green Version]
  52. Peng, J.; Zhang, W.; Wang, X.; Guo, G.; Zhang, F. Numerical simulation of effect of filler wire melting and filling mode on molten pool in laser welding. Hongwai yu Jiguang Gongcheng/Infrared Laser Eng. 2018, 47, 306005. [Google Scholar] [CrossRef]
  53. Xia, M.; Gu, D.; Yu, G.; Dai, D.; Chen, H.; Shi, Q. Porosity evolution and its thermodynamic mechanism of randomly packed powder-bed during selective laser melting of Inconel 718 alloy. Int. J. Mach. Tools Manuf. 2017, 116, 96–106. [Google Scholar] [CrossRef]
  54. Aucott, L.; Dong, H.; Mirihanage, W.; Atwood, R.; Kidess, A.; Gao, S.; Wen, S.; Marsden, J.; Feng, S.; Tong, M.; et al. Revealing internal flow behaviour in arc welding and additive manufacturing of metals. Nat. Commun. 2018, 9, 5414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Hojjatzadeh, S.M.H.; Parab, N.D.; Yan, W.; Guo, Q.; Xiong, L.; Zhao, C.; Qu, M.; Escano, L.I.; Xiao, X.; Fezzaa, K.; et al. Pore elimination mechanisms during 3D printing of metals. Nat. Commun. 2019, 10, 3088. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Ur Rehman, A.; Pitir, F.; Salamci, M.U. Full-Field Mapping and Flow Quantification of Melt Pool Dynamics in Laser Powder Bed Fusion of SS316L. Materials 2021, 14, 6264. [Google Scholar] [CrossRef]
  57. Ur Rehman, A.; Mahmood, M.A.; Pitir, F.; Salamci, M.U.; Popescu, A.C.; Mihailescu, I.N. Mesoscopic Computational Fluid Dynamics Modelling for the Laser-Melting Deposition of AISI 304 Stainless Steel Single Tracks with Experimental Correlation: A Novel Study. Metals 2021, 11, 1569. [Google Scholar] [CrossRef]
  58. Denlinger, E.R.; Irwin, J.; Michaleris, P. Thermomechanical modeling of additive manufacturing large parts. J. Manuf. Sci. Eng. Trans. ASME 2014, 136, 061007. [Google Scholar] [CrossRef]
  59. Li, C.; Gouge, M.F.; Denlinger, E.R.; Irwin, J.E.; Michaleris, P. Estimation of part-to-powder heat losses as surface convection in laser powder bed fusion. Addit. Manuf. 2019, 26, 258–269. [Google Scholar] [CrossRef]
  60. Agarwala, M.; Bourell, D.; Beaman, J.; Marcus, H.; Barlow, J. Direct selective laser sintering of metals. Rapid Prototyp. J. 1995, 1, 26–36. [Google Scholar] [CrossRef]
  61. Aboulkhair, N.T.; Everitt, N.M.; Ashcroft, I.; Tuck, C. Reducing porosity in AlSi10Mg parts processed by selective laser melting. Addit. Manuf. 2014, 1, 77–86. [Google Scholar] [CrossRef]
  62. Demir, A.G.; Previtali, B. Investigation of remelting and preheating in SLM of 18Ni300 maraging steel as corrective and preventive measures for porosity reduction. Int. J. Adv. Manuf. Technol. 2017, 93, 2697–2709. [Google Scholar] [CrossRef]
  63. Colopi, M.; Demir, A.G.; Caprio, L.; Previtali, B. Limits and solutions in processing pure Cu via selective laser melting using a high-power single-mode fiber laser. Int. J. Adv. Manuf. Technol. 2019, 104, 2473–2486. [Google Scholar] [CrossRef]
  64. Wilkes, J.; Hagedorn, Y.C.; Meiners, W.; Wissenbach, K. Additive manufacturing of ZrO2-Al2O3 ceramic components by selective laser melting. Rapid Prototyp. J. 2013, 19, 51–57. [Google Scholar] [CrossRef]
  65. Kempen, K.; Vrancken, B.; Buls, S.; Thijs, L.; Van Humbeeck, J.; Kruth, J.P. Selective Laser Melting of Crack-Free High Density M2 High Speed Steel Parts by Baseplate Preheating. J. Manuf. Sci. Eng. Trans. ASME 2014, 136, 061026. [Google Scholar] [CrossRef]
  66. Ali, H.; Ma, L.; Ghadbeigi, H.; Mumtaz, K. In-situ residual stress reduction, martensitic decomposition and mechanical properties enhancement through high temperature powder bed pre-heating of Selective Laser Melted Ti6Al4V. Mater. Sci. Eng. A 2017, 695, 211–220. [Google Scholar] [CrossRef]
  67. Gussone, J.; Hagedorn, Y.C.; Gherekhloo, H.; Kasperovich, G.; Merzouk, T.; Hausmann, J. Microstructure of γ-titanium aluminide processed by selected laser melting at elevated temperatures. Intermetallics 2015, 66, 133–140. [Google Scholar] [CrossRef]
  68. Ding, C.; Cui, X.; Jiao, J.; Zhu, P. Effects of substrate preheating temperatures on the microstructure, properties, and residual stress of 12CrNi2 prepared by laser cladding deposition technique. Materials 2018, 11, 2401. [Google Scholar] [CrossRef] [Green Version]
  69. Buchbinder, D.; Meiners, W.; Pirch, N.; Wissenbach, K.; Schrage, J. Investigation on reducing distortion by preheating during manufacture of aluminum components using selective laser melting. J. Laser Appl. 2014, 26, 012004. [Google Scholar] [CrossRef]
  70. Kruth, J.P.; Deckers, J.; Yasa, E.; Wauthlé, R. Assessing and comparing influencing factors of residual stresses in selective laser melting using a novel analysis method. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2012, 226, 980–991. [Google Scholar] [CrossRef]
  71. Leung, C.L.A.; Tosi, R.; Muzangaza, E.; Nonni, S.; Withers, P.J.; Lee, P.D. Effect of preheating on the thermal, microstructural and mechanical properties of selective electron beam melted Ti-6Al-4V components. Mater. Des. 2019, 174, 107792. [Google Scholar] [CrossRef]
  72. Zhang, K.; Liu, T.; Liao, W.; Zhang, C.; Yan, Y.; Du, D. Influence of laser parameters on the surface morphology of slurry-based Al2O3 parts produced through selective laser melting. Rapid Prototyp. J. 2018, 24, 333–341. [Google Scholar] [CrossRef]
  73. Antonsson, T.; Fredriksson, H. The effect of cooling rate on the solidification of INCONEL 718. Metall. Mater. Trans. B Process Metall. Mater. Process. Sci. 2005, 36, 85–96. [Google Scholar] [CrossRef]
  74. Wang, Z.; Guan, K.; Gao, M.; Li, X.; Chen, X.; Zeng, X. The microstructure and mechanical properties of deposited-IN718 by selective laser melting. J. Alloys Compd. 2012, 513, 518–523. [Google Scholar] [CrossRef]
  75. Gruber, K.; Smolina, I.; Kasprowicz, M.; Kurzynowski, T. Evaluation of Inconel 718 Metallic Powder to Optimize the Reuse of Powder and to Improve the Performance and Sustainability of the Laser Powder Bed Fusion (LPBF) Process. Materials 2021, 14, 1538. [Google Scholar] [CrossRef]
  76. Cleary, P.W.; Sawley, M.L. DEM modelling of industrial granular flows: 3D case studies and the effect of particle shape on hopper discharge. Appl. Math. Model. 2002, 26, 89–111. [Google Scholar] [CrossRef]
  77. Parteli, E.J.R.; Pöschel, T. Particle-based simulation of powder application in additive manufacturing. Powder Technol. 2016, 288, 96–102. [Google Scholar] [CrossRef]
  78. Cao, L. Numerical simulation of the impact of laying powder on selective laser melting single-pass formation. Int. J. Heat Mass Transf. 2019, 141, 1036–1048. [Google Scholar] [CrossRef]
  79. Juan, J.V.; Peter, N.Q. Thermophysical Properties; ASM Handbook Committee: Materials Park, OH, USA, 2018; Volume 17. [Google Scholar]
  80. Lee, Y.S.; Zhang, W. Modeling of heat transfer, fluid flow and solidification microstructure of nickel-base superalloy fabricated by laser powder bed fusion. Addit. Manuf. 2016, 12, 178–188. [Google Scholar] [CrossRef] [Green Version]
  81. Tang, M.; Pistorius, P.C.; Beuth, J.L. Prediction of lack-of-fusion porosity for powder bed fusion. Addit. Manuf. 2017, 14, 39–48. [Google Scholar] [CrossRef]
  82. Promoppatum, P.; Yao, S.C.; Pistorius, P.C.; Rollett, A.D. A Comprehensive Comparison of the Analytical and Numerical Prediction of the Thermal History and Solidification Microstructure of Inconel 718 Products Made by Laser Powder-Bed Fusion. Engineering 2017, 3, 685–694. [Google Scholar] [CrossRef]
  83. Rosenthal, D. Mathematical Theory of Heat Distribution During Welding and Cutting. Weld. J. 1941, 20, 220–234. [Google Scholar]
  84. Chen, Q.; Zhao, Y.Y.; Strayer, S.; Zhao, Y.Y.; Aoyagi, K.; Koizumi, Y.; Chiba, A.; Xiong, W.; To, A.C. Elucidating the Effect of Preheating Temperature on Melt Pool Morphology Variation in Inconel 718 Laser Powder Bed Fusion via Simulation and Experiment. Available online: https://www.sciencedirect.com/science/article/pii/S2214860420310149#bb8 (accessed on 30 April 2021).
  85. Paul, A.; Debroy, T. Free surface flow and heat transfer in conduction mode laser welding. Metall. Trans. B 1988, 19, 851–858. [Google Scholar] [CrossRef]
  86. Peng, J.; Hu, S.; Wang, X.; Wang, J.; Zhang, F. Effect of Filler Metal on Three-Dimensional Transient Behavior of Keyholes and Molten Pools in Laser Welding. Zhongguo Jiguang/Chinese J. Lasers 2018, 45, 0102003. [Google Scholar] [CrossRef]
  87. Abderrazak, K.; Bannour, S.; Mhiri, H.; Lepalec, G.; Autric, M. Numerical and experimental study of molten pool formation during continuous laser welding of AZ91 magnesium alloy. Comput. Mater. Sci. 2009, 44, 858–866. [Google Scholar] [CrossRef]
  88. Tenner, F.; Berg, B.; Brock, C.; Klämpfl, F.; Schmidt, M. Experimental approach for quantification of fluid dynamics in laser metal welding. J. Laser Appl. 2015, 27, S29003. [Google Scholar] [CrossRef]
  89. Eriksson, I.; Powell, J.; Kaplan, A.F.H. Melt behavior on the keyhole front during high speed laser welding. Opt. Lasers Eng. 2013, 51, 735–740. [Google Scholar] [CrossRef] [Green Version]
  90. Nakamura, H.; Kawahito, Y.; Nishimoto, K.; Katayama, S. Elucidation of melt flows and spatter formation mechanisms during high power laser welding of pure titanium. J. Laser Appl. 2015, 27, 032012. [Google Scholar] [CrossRef]
  91. Bidare, P.; Bitharas, I.; Ward, R.M.; Attallah, M.M.; Moore, A.J. Fluid and particle dynamics in laser powder bed fusion. Acta Mater. 2018, 142, 107–120. [Google Scholar] [CrossRef]
  92. Matsunawa, A.; Seto, N.; Mizutani, M.; Katayama, S. Liquid Motion in Keyhole Laser Welding; Laser Institute of America: Orlando, FL, USA, 2018; pp. G151–G160. [Google Scholar]
  93. Chang, B.; Allen, C.; Blackburn, J.; Hilton, P.; Du, D. Fluid Flow Characteristics and Porosity Behavior in Full Penetration Laser Welding of a Titanium Alloy. Metall. Mater. Trans. B Process Metall. Mater. Process. Sci. 2015, 46, 906–918. [Google Scholar] [CrossRef]
  94. Peng, J.; Li, L.; Lin, S.; Zhang, F.; Pan, Q.; Katayama, S. High-Speed X-Ray Transmission and Numerical Study of Melt Flows inside the Molten Pool during Laser Welding of Aluminum Alloy. Math. Probl. Eng. 2016, 2016, 1409872. [Google Scholar] [CrossRef] [Green Version]
  95. Kawahito, Y.; Nakada, K.; Uemura, Y.; Mizutani, M.; Nishimoto, K.; Kawakami, H.; Katayama, S. Relationship between melt flows based on three-dimensional X-ray transmission in situ observation and spatter reduction by angle of incidence and defocussing distance in high-power laser welding of stainless steel. Weld. Int. 2018, 32, 485–496. [Google Scholar] [CrossRef]
  96. Zhang, T.; Li, H.; Liu, S.; Shen, S.; Xie, H.; Shi, W.; Zhang, G.; Shen, B.; Chen, L.; Xiao, B.; et al. Evolution of molten pool during selective laser melting of Ti-6Al-4V. J. Phys. D Appl. Phys. 2019, 52, 055302. [Google Scholar] [CrossRef]
Figure 1. (a) Powder particle SEM and (b) discrete element modeling.
Figure 1. (a) Powder particle SEM and (b) discrete element modeling.
Materials 14 06683 g001
Figure 2. (a) Powder bed deposition; (b) deposited layer.
Figure 2. (a) Powder bed deposition; (b) deposited layer.
Materials 14 06683 g002
Figure 3. Temperature-dependent properties: (a) solid state to liquid state, (b) density and molar volume, (c) average expansion coefficient and volume change, (d) electrical resistivity and thermal conductivity, (e) total viscosity and liquid diffusivity, (f) Poisson’s ratio and liquid viscosity, (g) electrical conductivity and young’s modulus, and (h) bulk modulus and shear modulus.
Figure 3. Temperature-dependent properties: (a) solid state to liquid state, (b) density and molar volume, (c) average expansion coefficient and volume change, (d) electrical resistivity and thermal conductivity, (e) total viscosity and liquid diffusivity, (f) Poisson’s ratio and liquid viscosity, (g) electrical conductivity and young’s modulus, and (h) bulk modulus and shear modulus.
Materials 14 06683 g003
Figure 4. (a) Faucet design, (b) ball within cage design, (c) lattice structure manufactured specimens, (d) ENAVISION 250 used for validation, (e) machine breakdown.
Figure 4. (a) Faucet design, (b) ball within cage design, (c) lattice structure manufactured specimens, (d) ENAVISION 250 used for validation, (e) machine breakdown.
Materials 14 06683 g004
Figure 5. Single track dynamic viscosity profile at (a) 200, (b) 270, (c) 300, (d) 400 µs and (e) shows the experimental comparison through the SEM of the irradiated single track.
Figure 5. Single track dynamic viscosity profile at (a) 200, (b) 270, (c) 300, (d) 400 µs and (e) shows the experimental comparison through the SEM of the irradiated single track.
Materials 14 06683 g005
Figure 6. Single track melt pool profile at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Figure 6. Single track melt pool profile at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Materials 14 06683 g006
Figure 7. Stream traces of single-track flow at (a) 50, (b) 140, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Figure 7. Stream traces of single-track flow at (a) 50, (b) 140, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Materials 14 06683 g007
Figure 8. Longitudinal velocity vectors of single track at (a) 40, (b) 130, (c) 220, (d) 310, (e) 400, (f) 490 µs.
Figure 8. Longitudinal velocity vectors of single track at (a) 40, (b) 130, (c) 220, (d) 310, (e) 400, (f) 490 µs.
Materials 14 06683 g008
Figure 9. (a) Mass flow rate forward and backwards, (b) flow streams at 410 µs, (c) X-ray synchrotron imaging of melt pool [27], and (d) schematics of melt pool from the same [27].
Figure 9. (a) Mass flow rate forward and backwards, (b) flow streams at 410 µs, (c) X-ray synchrotron imaging of melt pool [27], and (d) schematics of melt pool from the same [27].
Materials 14 06683 g009
Figure 10. Stream traces of the cross-section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Figure 10. Stream traces of the cross-section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Materials 14 06683 g010
Figure 11. (a) OM of the single track; (b) stream traces of the cross-section at 300.
Figure 11. (a) OM of the single track; (b) stream traces of the cross-section at 300.
Materials 14 06683 g011
Figure 12. Velocity vectors of the cross section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Figure 12. Velocity vectors of the cross section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Materials 14 06683 g012
Figure 13. The isometric view of the melt pool (a) 50, (b) 120, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Figure 13. The isometric view of the melt pool (a) 50, (b) 120, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Materials 14 06683 g013
Figure 14. Single track melt pool profile at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Figure 14. Single track melt pool profile at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Materials 14 06683 g014
Figure 15. Stream traces of single track in 3D at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Figure 15. Stream traces of single track in 3D at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Materials 14 06683 g015
Figure 16. Longitudinal velocity vectors of single track at (a) 40, (b) 130, (c) 220, (d) 310, (e) 400, (f) 490 µs.
Figure 16. Longitudinal velocity vectors of single track at (a) 40, (b) 130, (c) 220, (d) 310, (e) 400, (f) 490 µs.
Materials 14 06683 g016
Figure 17. (a) Mass flow rate forward and backwards, (b,c) melt pool passing from cross section at different time.
Figure 17. (a) Mass flow rate forward and backwards, (b,c) melt pool passing from cross section at different time.
Materials 14 06683 g017
Figure 18. Stream traces of the cross-section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Figure 18. Stream traces of the cross-section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Materials 14 06683 g018
Figure 19. Velocity vectors of the cross section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Figure 19. Velocity vectors of the cross section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Materials 14 06683 g019
Figure 20. The isometric view of the melt pool with 500 °C pre-heating at (a) 50, (b) 120, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Figure 20. The isometric view of the melt pool with 500 °C pre-heating at (a) 50, (b) 120, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Materials 14 06683 g020
Figure 21. Single track melt pool profile with 1000 °C pre-heating at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Figure 21. Single track melt pool profile with 1000 °C pre-heating at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Materials 14 06683 g021
Figure 22. Streamtraces of single track in 3D at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Figure 22. Streamtraces of single track in 3D at (a) 50, (b) 140, (c) 230, (d) 320 (e) 410 and (f) 500 µs.
Materials 14 06683 g022
Figure 23. Stream traces of the cross-section with 1000 °C pre-heating at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Figure 23. Stream traces of the cross-section with 1000 °C pre-heating at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Materials 14 06683 g023
Figure 24. (a) Mass flow rate forward and backwards with 1000 °C pre-heating, (b,c) melt pool passing from cross section at different time.
Figure 24. (a) Mass flow rate forward and backwards with 1000 °C pre-heating, (b,c) melt pool passing from cross section at different time.
Materials 14 06683 g024
Figure 25. Velocity vectors of the cross section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Figure 25. Velocity vectors of the cross section at (a) 50, (b) 120, (c) 210, (d) 300, (e) 390, (f) 480 µs.
Materials 14 06683 g025
Figure 26. The isometric view of the melt pool with 1000 °C pre-heating (a) 50, (b) 120, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Figure 26. The isometric view of the melt pool with 1000 °C pre-heating (a) 50, (b) 120, (c) 230, (d) 320, (e) 410, (f) 500 µs.
Materials 14 06683 g026
Figure 27. Mass flow rate forward and backwards from cross section at different time at RT in pink and yellow line, at 500 °C pre-heating in green and blue line and at 1000 °C pre-heating, in black and red line.
Figure 27. Mass flow rate forward and backwards from cross section at different time at RT in pink and yellow line, at 500 °C pre-heating in green and blue line and at 1000 °C pre-heating, in black and red line.
Materials 14 06683 g027
Table 1. IN718 particle size (µm).
Table 1. IN718 particle size (µm).
D10D50D90
192941
Table 2. Chemical composition of In718.
Table 2. Chemical composition of In718.
NiCrFeAlCoCuMnMoNbSiTaTiBC
52.519.517.30.50.50.150.1753.055.1250.1750.0250.90.0030.04
Table 3. Technical Specification ERMAKSAN Enavision 250.
Table 3. Technical Specification ERMAKSAN Enavision 250.
GENERAL SPECIFICATIONENAVISION 250
Production Volume (mm3)250 × 250 × 300 (9.8 × 9.8 × 11.8 inch)
Adjustable Layer Height20–100 μm (0.0007–0.004 inch)
Laser TypeFiber Laser
Laser Power500W
Scanning SpeedUp to 11 m/s (433.07 inch)
Scanning System3D Dynamic Focused Scanning System
Dimension(L × W × H) 2700 × 1440 × 2030 (106.3 × 56.7 × 79.9 inch)
Voltage400 V, 3 PH, 50/60 Hz
Current32 A
Inert GasArgon/Nitrogen
02 Level100 ppm
Vacuum PompYes
Operating SystemOperating SystemWindows 10/X
CONTROL UNIT
Control SystemBeckhoff Industrial PC
ProcessorProcessorIntel i5–i7
Operating SystemWindows 10/X
HMI15.6 inch, Touch Operated
SOFTWARE
Data Preparation SoftwareMaterilliase Magics and Modules
Data Processing SoftwareErmaksan Build Processor
Supported File TypesSTL, 3MF, AMF, DAE, FBX, VRML.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ur Rehman, A.; Pitir, F.; Salamci, M.U. Laser Powder Bed Fusion (LPBF) of In718 and the Impact of Pre-Heating at 500 and 1000 °C: Operando Study. Materials 2021, 14, 6683. https://doi.org/10.3390/ma14216683

AMA Style

Ur Rehman A, Pitir F, Salamci MU. Laser Powder Bed Fusion (LPBF) of In718 and the Impact of Pre-Heating at 500 and 1000 °C: Operando Study. Materials. 2021; 14(21):6683. https://doi.org/10.3390/ma14216683

Chicago/Turabian Style

Ur Rehman, Asif, Fatih Pitir, and Metin Uymaz Salamci. 2021. "Laser Powder Bed Fusion (LPBF) of In718 and the Impact of Pre-Heating at 500 and 1000 °C: Operando Study" Materials 14, no. 21: 6683. https://doi.org/10.3390/ma14216683

APA Style

Ur Rehman, A., Pitir, F., & Salamci, M. U. (2021). Laser Powder Bed Fusion (LPBF) of In718 and the Impact of Pre-Heating at 500 and 1000 °C: Operando Study. Materials, 14(21), 6683. https://doi.org/10.3390/ma14216683

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop