Impact of fume particles in the keyhole vapour

During laser material processing, the energy of the laser beam needs to be efficiently and constantly transported to the processing zone to guarantee constant processing. However, spatters or ejected particles from the keyhole can absorb and scatter the laser energy leading to inhomogeneous heat input and can initiate defect occurrence like pore formation. The impact of ejected particles from the keyhole on the energy transport of the laser beam is not completely understood since they are difficult to observe due to the small size and high speeds of the ejections. In this work, the particle characteristics were derived from a simulation of the keyhole wall movement. The behavior of the calculated particles in a side shielding gas jet was calculated to derive the height, at which the particles leave the laser beam and are not interrupting the laser energy transfer to the processing zone. Low impulse values of the particles were calculated e.g., at defocusing positions slightly underneath the material surface, where also highest melt pool sizes were found. These observations indicate that the fume particles can be one reason to limit the energy delivery. An efficiency increase can be achieved by adjusting the keyhole parameters to a more stable keyhole.


Introduction
Laser deep penetration welding is a highly complex process due to the occurrence of high temperatures leading to melting, vaporization and the interaction of solid, liquid and vaporized material. The characteristic vapor channel, the keyhole, existing during processing, is kept open during processing due to the recoil pressure induced by the vaporization of the keyhole walls due to laser irradiation. The vapor outflow is known to be unstable, e.g., pore formation requires a reversed gas flow of ambient gas into the keyhole (e.g., [1]). Fabbro et al. [2] found that depending on the process parameters the vapor outflow can show a pulsed characteristic. The keyhole itself is highly dynamic showing moving melt waves on the keyhole front wall [3] and small humps appearing and moving at frequencies up to 1 MHz [4]. Those fluctuations are thought to be the origin of typical defects, such as spatters and pores (e.g., demonstrated by Volpp [5]). Besides spatter emission, the vapor plume appears above the keyhole exit due to vaporization happening on the keyhole walls during processing. This dynamic vapor characteristics could be related to fluctuations in the keyhole [6], which means that the keyhole behavior seems to define the vapor outflow characteristics.
When welding with a CO 2 -laser (wavelength 10.6 µm) instead of solid-state lasers (wavelength around 1 µm), gases in the laser beam can be even ionized leading to plasma creation. The plasma can absorb laser beam energy, which cannot be transported to the processing zone and can interrupt or even stop the processing [7]. It has been shown that when processing with solid-state lasers, plasma effects can be neglected [8]. Kawahito et al. [9] used spectroscopic measurements when welding stainless steel with a 10-kW fiber laser and found that the induced plume mainly consists of the alloying elements from the base material and that the plume is not highly ionized to form a plasma.
However, it was observed that the plume affects the laser energy transport to the keyhole. Michalowski et al. [10] observed that more than 30% of the energy of the laser beam can be lost due to scattering and absorption in the plume. Scholz et al. [11] observed that fluctuations of the particle appearance lead to inherent variations of the energy delivery to the processing zone. This, in consequence, could in turn influence the keyhole wall fluctuations.
Shcheglov et al. [12] observed two characteristic parts of the plume when laser welding mild steel (Fig. 1). The lower plume glows bright and shows fluctuation frequencies of several kHz, while the upper part shows weaker emissions and is related to the laser beam path.
The size of the particles, which are present in the plume, were measured by different authors (Fig. 2).
The detected particle sizes vary from 10 nm (Scholz et al. [13]) to 192 µm (Michalowski et al. [10]). Due to the wide range, the particles could origin from different mechanisms. The small particles can be related to the vaporization of material from the keyhole wall and the re-condensation during their outflow on other particles (e.g., [13]). The bigger particles could be also micro-spatters detached from the keyhole walls. Particles in the laser beam, in general, influence the laser energy delivery by absorption of energy on the particle surface or scattering (e.g., [14]), which both leads to a lower energy input into the keyhole. The detection of the particles was partly possible; their origin could not be sufficiently observed yet. Modeling approaches were conducted to simulate the behavior of the keyhole and melt pool. Pointto-point calculations enabled to model keyhole geometries depending on processing parameters (e.g., [15]) or could help predict the porosity [16]. Simplified models can help to predict the keyhole depth [17], which might not fully represent the dynamic observed behavior of the keyhole (e.g., observed in ice [18] or during copper welding [19]).
It is known that a shielding gas flow can blow away the particles, which are ejected from the keyhole, which can lead to a more efficient and constant energy absorption. Especially when processing with a CO 2 laser, where plasma forming is expected, an application of shielding gas can significantly improve the energy input by blowing the plasma away from the laser beam circumference. When processing with a fiber laser, Mahrle et al. [20] found that a gas flow of 2 m/s applied by a shielding gas nozzle from a rectangular position to the straight upwards directed vapor plume is sufficient to suppress interactions between the laser beam and the plume. It seems that a reduction of the particle density in the plume reduces the energy absorption and/or scattering by the plume.
However, it is not known, where and how these particles are created and how they influence the welding results. It is known that at negative focal positions (focal point below the material surface) can be beneficial to achieve a higher heat input into the keyhole, which is assumed to be related to increased multiple reflections in the keyhole or a more stable keyhole behavior (e.g., [21]). However, it is possible that the characteristics of the vapor can affect the energy delivery as well. Knowledge about the creation and transporting of those particles could help find methods to avoid particle creation in general and thereby increase the process efficiency. In addition, knowledge about the occurrence of particle sizes and their impact on the laser beam absorption could help define critical particle sizes and open the discussion of influencing the keyhole wall characteristics to avoid creating critical particle sizes that induce process changes and defects. Therefore, in this work, the possible origin, their size distribution and characteristics of particles that are created on the dynamic keyhole wall are examined and related to experimental observations.

Experiment conduction
A single-mode fiber laser (IPG YLR-1000SM) was used to produce the laser beam, which was collimated (160 mm focal length) and focused (560 mm focal length). The resulting focal beam diameter was 52 µm at a Rayleigh length of 1.4 mm. The focal spatial laser intensity distribution shows a Gaussian beam profile. Reference parameters were defined (Table 1).
Based on these parameters, variations of the laser power (from 800 to 1000 W), the welding speed (from 0.5 to 2 m/min) and the focal position (− 1 to + 1 mm) were conducted. No shielding gas was used. Base material was stainless steel 1.4301, while specimens were of the size 40 mm x 6 mm x 100 mm. 80 mm long weld seams were produced.

Weld seam evaluation
Cross-section polishes were produced to evaluate the weld seam dimension. The weld seam area was derived as visualized in Fig. 3.

Keyhole modeling
A simplified semi-analytical model was developed to calculate the keyhole properties, which consist of two parts. An analytical model calculates the quasi-static conditions of the keyhole, while a subsequent time-dependent calculation simulates the dynamic behavior of the keyhole wall movement. The model was described in detail in [22].
The quasi-static model is based on a segmentation of the keyhole into ten cylindrical elements (Fig. 4). Initially, the depth of the keyhole must be given and was derived from experimental measurements. Based on the laser beam caustic, the energy input in every element is calculated assuming Fresnel absorption on the inclined keyhole walls (Fig. 4).
The spatial laser intensity distribution of the laser beam was considered a Gaussian distribution i gs and was calculated in every element depending on the actual beam diameter and was calculated using the following equation: with the laser power P L , the beam radius r and the focal laser beam radius r L0 .
The energy balance equation was considered including the contribution of the energy input from the laser beam and energy losses from heat conduction into the surrounding material and due to vaporization effects. The solution of the energy equation gives the actual temperature in each keyhole element. The temperature values are needed to calculate the pressure balance equation. In all keyhole elements, the vapor pressure was calculated based on the actual temperature in the element. The vapor pressure keeps the keyhole element open against pressure of the surrounding melt pool due to surface tension. The pressure balance calculation gives the actual radius of the cylindrical keyhole element. Therefore, the quasi-static keyhole shape was derived, which was used as boundary condition for the temporal calculation. The dynamic model used the quasi-static geometrical keyhole data as starting conditions of the movement of the keyhole walls of the single elements. To describe the keyhole elements as a dynamic system, a set of differential equations was used. The pressure differential equation contains the impact of evaporation and the change of the element volume at element expansion or shrinkage. Evaporation was based on the actual energy input by the laser beam, considering Fresnel absorption effects and multiple reflections inside the keyhole using an iterative routine. The second differential equation described the radius development of the elements, including the effects of the ablation pressure, which increases the radius and the counteracting surface tension pressure initiated by the melt pool and the bulging of the melt pool at the surface. This description gives the characteristics (frequencies, amplitudes) of each single keyhole element. The solution of the equations (solver: ode45), calculated in Matlab Simulink (Version R2009a), gives the temporal development of the dynamic properties of   Table 2 shows the parameters used for the calculation.
The calculated keyhole wall movements were compared to acoustic emissions recorded by a microphone during the processing (Fig. 5).
The average of the dynamic characteristic properties was evaluated at different parameter sets to be compared to the experimental evaluations. The dynamic values were used to calculate the impulse I of the particle as with the amplitude A, the density met and the particle speed v part . The particle speed was calculated assuming the maximum wall speed being transferred onto the particle by

Particle flow modeling
An additional model was developed to analyze the metal particle flow paths that were ejected from the keyhole walls after exiting the keyhole at a defined side gas application. The particle paths were calculated at the following assumptions: • all particles exit the keyhole vertically to the material surface, which can be assumed due to the strong vapor outflow from the keyhole, • all particles are perfect spheres due to the surface tension.
When exiting the keyhole, the ejected particle experiences the flow of a horizontally applied gas jet with a gas flow speed of 2 m/min (Fig. 6).
This side gas flow induces a pressure p on the particle of with air being the density of the gas and v air being the speed of the air. The pressure applied on the projected area of a particleB (4) p = air 2 ⋅ air 2 , induces a forceF on the particle with the massm leading to an acceleration a of the particle According to the relations at accelerated movement, the time till leaving the laser beam  when assuming a laser beam with a constant diameter 2 • r L0 , which can be assumed since the calculation is done close to the focal plane. When considering the initial speed of the particle assumed in vertical upward direction when exiting the keyhole v part , the height h above the material surface was determined, at which the particle leaves the laser beam and is not an obstacle for the transfer of laser energy to the processing zone anymore to

Weld seam characteristics
The weld seam areas at different parameter sets were evaluated as described in the Methods for comparison of different processing parameters. Figure 7 shows the measured weld seam areas derived from cross-section polishes.
At a higher heat input (lower velocity or higher laser power), the weld seam area increases. At varied focal position, the highest value was found at −1 mm defocusing.

Keyhole characteristics
The keyhole simulation calculated the characteristics of the temporal keyhole wall development in the single keyhole elements. The average of all amplitudes and frequencies calculated in all keyhole elements including the range of occurring data points are visualized in Figs. 8 and 9.
The occurring amplitudes (Fig. 8) show the lowest values at defocusing underneath the material surface, while the frequencies of the keyhole wall movement are comparably high (Fig. 9). These results describe the dynamic movement of the keyhole wall elements depending on the processing parameters.
From the calculated frequencies and amplitudes of the keyhole elements, the impulse was calculated using Eq. 2 (Fig. 10). The impulse results from the keyhole wall movement and shows lowest values at high velocities and defocusing underneath the material surface.

Particle paths
The calculation of the particle path of the material exiting the keyhole with the vapor when applying a side gas jet was conducted. This simple approach makes it possible to determine the height above the material surface, at which a particle with the diameter of the keyhole wall amplitude in the keyhole at its calculated impulse (Fig. 10) leaves the laser beam (Fig. 11) and can give indications how much material is inside the vapor plume that can affect the laser energy transmission.
At lower welding speeds, particles can be removed from the laser beam by shielding gas more efficiently, while the  laser power showed minor impact on the particle paths. At a focal position slightly below the material surface (negative values), the particle height was calculated in average the lowest, which indicates the setting with the most efficient removal of particles from the laser beam.

Discussion
In the literature, it was often observed that when positioning the laser beam focal position slightly underneath the material surface, a higher energy input efficiency into the keyhole can be achieved (e.g., [21]). In the actual work, this phenomenon could be also seen in the measured weld seam areas (Fig. 7). Possible reasons for this phenomenon could be: • an increased energy absorption due to increased amounts of multiple reflection, • a more stable keyhole behavior, leading to a more constant absorption on the keyhole walls, or, • the effect of different ejected plume particles and their impact on energy absorption and scattering.
Especially, the plume content and its impact on the processes are barely considered, since the measuring of the particle characteristics can be difficult due to high temperatures and pressures and the fast processes, as well as the limited access to their origin.
The proposed model in this paper calculated the expected ranges of particles that can be ejected from the keyhole wall at its dynamic movement. The keyhole wall movement creates the material bulging into the keyhole area, while the created vapor gas flow can be able to detach the small particles from the keyhole wall [5].
Since in other works (Fig. 2) a big range of particle sizes were found in the plume, it can be assumed that different origins are possible. Small particles are most likely re-condensated from vapor content, the slightly larger particles could also origin from keyhole wall detachments as micro-spatters as derived from the results also in this work.
The simulation results in this work show that different characteristics (especially the impulse) of detached particles can occur depending on the process parameters. Since usually shielding gas is used to blow away the plume, a standard shielding gas flow was considered that blows away the plume including the particles that are ejected from the keyhole. In general, a higher impulse leads to a later leaving of the laser beam and can therefore lead to more energy absorption or scattering in the plume and therefore to a lower laser energy transport to the keyhole walls.
Tendencies show that at a focal position underneath the material surface (Fig. 10), the particles show a comparably low impulse, which leads to a fast removal of the particles by the shielding gas. This observation shows that the particle characteristics can possibly be an impacting factor on the energy absorption efficiency and can partly explain the increased absorption at negative focal positions. This observation shows that the particle scattering can have an impact on the laser energy absorption leading to different keyhole characteristics and absorption mechanisms.

Conclusions
Comparing experimental results of weld seam area dimensions, which are a measure of energy absorption, and particle detachment simulations from the keyhole wall, it can be concluded that the particles that are ejected from the keyhole wall can affect the laser beam energy transport to the processing zone. It seems that the ejected particles are part of the mechanism that leads to the increased energy absorption when e.g., the focal position is positioned below the material surface.