Humping Formation and Suppression in High-Speed Laser Welding

Increasing welding speed can promote the productivity of laser welding. However, humping defects often occur, which limits the application of this strategy. The existing explanations for the humping formation remain vague, and mitigation and suppression methods are limited. In this research, high-speed imaging experiments and numerical simulation of the high-speed laser welding process are performed. Through careful examination, the humping phenomenon is explained. At high welding speed, the high-speed melt flow caused by recoil pressure is hindered by the solidified region in the melt pool, leading to the occurrence of a swelling. The swelling then grows, forming a valley in front of the swelling under the effect of surface tension. The solidification of the valley results in the occurrence of a second swelling. This process repeats and humping defect forms. Marangoni force and viscous force also have influence on this process. In addition, it is found that adding a Tungsten Inert Gas arc behind the laser beam can effectively suppress the humping.


Introduction
Because of the high-power density of a laser beam, generally laser welding has higher welding speed than arc welding, and thus can achieve higher productivity [1][2][3][4]. However, the productivity of laser welding cannot be improved without limit by simply increasing the welding speed and laser power proportionately. It is found that with the increase of welding speed and laser power, humping defects, which can be described as periodic undulations consisting of a series of swellings and valleys, occur easily on the weld seam [5][6][7]. The occurrence of humping severely deteriorates the product properties, and limits the further increase of welding speed.
Studies have been performed to understand the formation mechanism of humping in high-speed laser welding. Some researchers investigated the behavior of the melt pool in high-speed laser welding through high-speed imaging experiments and numerical simulations and explained the formation of humping with the Rayleigh instability model [8][9][10][11][12][13]. The Rayleigh instability model is a common explanation for the formation of humping. According to this theory, the melt flow in the melt pool is analogous to a cylindrical inviscid fluid jet freely suspended in space, and the lateral surface of the cylindrical jet becomes unstable once its length L exceeds its circumference 2πR. With the existence of tiny perturbations, the radius of the jet becomes non-uniform, so the additional pressure caused by surface tension also becomes non-uniform, leading to the periodic swelling and narrowing of the jet. This theory can explain the formation of humping to some extent, such as the periodic occurrence of swelling and narrowing in the melt pool, but it is a highly idealized model, and many important factors influencing the dynamic behavior of melt pool are ignored, such as the effects of viscosity and solidification. Hence, there is an ongoing debate about the validity of this explanation, and it needs further improvement [14].
There are also some other explanations for humping formation in high-speed laser welding. Kawahito et al. [15] observed the melt pool behavior through high-speed imaging and believed that the humping formation was dependent upon several dynamic or static factors such as melt volume, strong melt flow, solidification, and surface tension. Ai et al. [16] speculated that the collision of the backflow with the fast rearward melt flow provided the condition for the wave formation in the melt pool, which led to the occurrence of humping. Though these explanations were proposed according to the observation of the melt pool image, they were still vague. The humping formation process was still not described in detail. Therefore, the formation mechanism of humping in high-speed laser welding still needs further research.
In order to improve the productivity of laser welding while maintaining product quality, some methods have been proposed to suppress humping, but existing methods are still limited. Kawahito et al. [15] found that humping could be reduced by performing inner defocused welding, but the penetration depth would decrease in this way. Cornell et al. [17] used laser wobble welding to reduce humping, but this method only mitigated the undulation of the weld seam surface, and humping was not suppressed completely.
In this research, high-speed imaging experiment and numerical simulation of highspeed laser welding process are performed. Through the analysis of the experimental and simulation results, the humping formation process is explained. In addition, the effects of surface tension, Marangoni force, and viscous force on humping formation are studied. Finally, based on the analysis of the humping formation mechanism, a humping suppression method is proposed by adding a TIG (Tungsten Inert Gas) arc behind the laser beam, and humping is suppressed successfully in this way.

Experimental Setup and Procedure
A high-speed imaging experiment is conducted to observe the dynamic behavior of the melt pool in high-speed laser welding. The welding method is bead-on-plate welding. The workpieces are 304 stainless steel sheets with the following dimensions: length 200 mm × width 50 mm × thickness 1 mm. Before the welding process, the surfaces of workpieces are polished with sandpaper and cleaned with absolute ethyl alcohol. A MFSC-4000W single-mode fiber laser (MAX, Shenzhen, China) and a YC52 laser processing head (PRECITEC, Gaggenau, Germany) are integrated to output a laser beam. The maximum power of the laser is 4000 W, and the wavelength of the output laser is 1070 ± 10 nm. The diameter of the laser spot is 0.4 mm. The laser processing head is installed on a MH24 robot (YASKAWA, Kitakyushu, Japan). Argon is used as the shielding gas, and its flowrate is 20 L/min. A Memrecam HX-6 high-speed camera (NAC, Tokyo, Japan)is used to capture images of the melt pool during the welding process, and its frame rate is set to 5 kHz. The axis of the camera forms a 45 • angle with the horizontal plane. Since the reflected laser and bright laser-induced vapor plume in the image can disturb the observation of the melt pool, a CAVILUX HF pulsed high power diode laser light source (CAVITAR, Tampere, Finland) is used as an active light source for illumination, and its central wavelength is 810 nm, and a narrow band-pass filter with the central wavelength of 810 nm is attached on the camera lens. In this way, the interference factors in the image can be weakened, and the melt pool can be clearer. The configuration of the experiment is shown in Figure 1.
Welding parameters used in the experiment are listed in Table 1. The linear energy density of the laser is kept consistent for the convenience of comparing the behavior of the melt pool with similar penetration depth. Welding parameters used in the experiment are listed in Table 1. The linear energy density of the laser is kept consistent for the convenience of comparing the behavior of the melt pool with similar penetration depth.

Assumptions and Governing Equations
In order to understand the dynamic behavior of the melt pool in high-speed laser welding in detail, a simulation model is built. In the simulation, fluids are assumed to form an incompressible laminar flow with Newtonian viscosity. Solutions are found for the governing equations of mass conservation, momentum conservation (Navier-Stokes), and energy conservation in the simulation.
Mass conservation equation: Momentum conservation equation: Energy conservation equation:

Assumptions and Governing Equations
In order to understand the dynamic behavior of the melt pool in high-speed laser welding in detail, a simulation model is built. In the simulation, fluids are assumed to form an incompressible laminar flow with Newtonian viscosity. Solutions are found for the governing equations of mass conservation, momentum conservation (Navier-Stokes), and energy conservation in the simulation.
Mass conservation equation: Momentum conservation equation: Energy conservation equation: where ρ, → v , p, µ, → S f , H, k, T, and S h are density, velocity, pressure, dynamic viscosity, external body force, enthalpy, thermal conductivity, temperature, and heat sources of fluids, respectively.
In addition, the multiphase flow problem is involved since both the gas phase and liquid phase of fluids are included in the simulation, so the VOF (Volume of Fluid) method is adopted to model the multiphase flow problem. VOF equation: where the subscript i represents the i-th phase of fluids, α i is the volume of fluid of the i-th phase, the sum of α i for all phases is 1, and . m i is the mass change of the i-th phase caused by phase transformation.

Melting and Solidification
The enthalpy-porosity technique [18] is adopted to model melting and solidification. The solid-liquid interface is not traced explicitly. A variable called liquid fraction β is introduced, and it represents the ratio of the liquid phase in the cell. The enthalpy of the material is composed of sensible enthalpy h and latent heat ∆H: where where h ref is the reference enthalpy, c p is the specific heat, L is the latent heat of the material, and β is defined as where T s and T l are the solidus temperature and the liquidus temperature of the material, respectively. The momentum source due to the solidification and melting is given as where C K is a small number to prevent division by zero, and C mush is the mushy zone constant.

Vaporization and Condensation
The condensation process is thought to have little influence on the fluid flow, so it is ignored in the simulation. The mass loss of the liquid phase due to vaporization is written as [19] q where L v is the latent heat of vaporization, T v is the boiling temperature, M is the molar mass, R is the universal gas constant, and P 0 is the ambient pressure.

Surface Tension and Recoil Pressure
The pressure difference between the liquid and gas on either side of the free surface due to surface tension is P σ = σK (11) where K is the surface curvature, and σ is the surface tension coefficient: The surface tension coefficient is calculated as [12,20,21] where σ 0 = 1.77 N·m −1 is the surface tension coefficient of pure SUS304 stainless steel at the solidus temperature, A = −0.43 × 10 −3 N·m −1 ·K −1 is the temperature coefficient of surface tension coefficient of pure SUS304 stainless steel, Γ S = 1.3 × 10 −5 mol·m −2 and Γ O = 2.03 × 10 −5 mol·m −2 are the surface excesses at saturation of sulfur and oxygen, are the equilibrium constants for segregation of sulfur and oxygen, respectively, k 1,S = 0.00318 and k 1,O = 0.0138 are the entropy factors of sulfur and oxygen, respectively, ∆H S = −1.66 × 10 5 J·mol −1 and ∆H O = −1.463 × 10 5 J·mol −1 are the standard heat of absorption of sulfur and oxygen, respectively, α S = 0.03 and α O = 0.01 are the weight percentages of soluble sulfur and oxygen, respectively. According to Equation (12), the temperature coefficient of surface tension coefficient is The dependence of σ and dσ/dT on temperature is shown in Figure 2. It can be found that the value of dσ/dT is positive when the temperature is below about 2300 K, and it turns negative with the further increase of temperature.
The pressure difference between the liquid and gas on either side of the free surface due to surface tension is where is the surface curvature, and is the surface tension coefficient: The surface tension coefficient is calculated as [12,20,21]  = 0.03 and = 0.01 are the weight percentages of soluble sulfur and oxygen, respectively. According to equation (12), the temperature coefficient of surface tension coefficient is The dependence of and d d ⁄ on temperature is shown in Figure 2. It can be found that the value of d d ⁄ is positive when the temperature is below about 2300 K, and it turns negative with the further increase of temperature. The recoil pressure is calculated as follows The recoil pressure is calculated as follows

Boundary Condition
The heat boundary condition of the top free surface is influenced by laser heat flux, the heat dissipation by convection and radiation, and the heat loss by vaporization: where q laser is the absorbed laser energy by the workpiece, h A is the convection coefficient, T ∞ is the ambient temperature, ε r is the emissivity, σ s is the Stefan-Boltzmann constant, and q vap is the heat loss due to vaporization. For the heat boundary conditions of other free surfaces, only the heat dissipation by convection and radiation are taken into consideration.

Laser Heat Source
The surface power distribution of the laser is assumed to be a Gaussian function which has a mathematical form like the following equation where Q is the total laser power, the coefficient 3 means that 95% of the laser power exists in the area of radius r 0 , r 0 = 0.2 mm denotes the focal radius of the laser beam, x 0 is the welding starting position on the x axis, and u w is the welding speed along x axis. In each step of the calculation, the free surface between metal and gas is searched, and then the laser heat is applied on corresponding cells.

Numerical Consideration
A symmetry calculation domain is built and its size is length 12 mm × width 1 mm × height 2 mm. The heights of the metal zone and the gas zone are both 1 mm, as shown in Figure 3. The calculation domain is divided into uniform cubic grids with edge length of 0.025 mm, so the total number of grids is 1,536,000. The welding direction is along x axis and the starting and ending positions are x = 1.5 mm and x = 11.5 mm, respectively. The physical property parameters of SUS304 stainless steel and other parameters used in the simulation are listed in Table 2. The commercial computational fluid dynamics software ANSYS Fluent is adopted, and the PISO (Pressure-Implicit with Splitting of Operators) scheme is chosen to solve this model. The time step is set to 1 × 10 −6 s. The simulation program runs on the High Performance Computing (HPC) cluster supported by the Center of High performance computing, Tsinghua University [22]. 56 Intel(R) Xeon(R) CPU E5-2680 v4 processors are used in parallel.     [21].

Validation of the Simulation Model
Appearances of the weld seam and melt pool at different welding parameters in the experiments are shown in Figure 4. It can be found that at a welding speed of 16 m/min, no humping occurs on the weld seam, and the melt pool is relatively short and flat. At a welding speed of 24 m/min, humping occurs on the weld seam, and the melt pool is long, and obvious swellings appear in it. Two simulation cases are performed using welding parameters in Table 1, and appearances of the weld seam in the simulation results are shown in Figure 5. It can be found that no humping occurs at a welding speed of 16 m/min, and humping occurs at a welding speed of 24 m/min, which agrees with the experimental results. In addition, cross sections of the weld seams in experiments and simulations are given in Figure 6, and width and depth of the weld seams are listed in Table 3, and the simulation results are close to the experiment results. Therefore, the validity of the simulation model is verified. simulation results are close to the experiment results. Therefore, the validity of the simulation model is verified.

Melt Pool Behavior and Humping Formation
The behavior of melt pool during humping formation process at welding speed of 24 m/min is observed through high-speed imaging, as shown in Figure 7. In the image at t 0 , it can be observed that a long and narrow groove exists behind the laser irradiation region. Fabbro [10] also noticed this phenomenon and referred to this groove as "elongated keyhole". In addition, the melt pool is also elongated and high-speed melt flow can be noticed. A swelling of molten metal (swelling 1) exists at the tail of the melt pool. The melt flow in front of this swelling is thin and forms a valley (valley 1). Most of the melt here has solidified, and the flow channel has narrowed obviously. At t 0 + 4 ms, molten metal accumulates in front of valley 1 and forms a new swelling (swelling 2). At t 0 + 11 ms, a new valley (valley 2) appears in front of swelling 2, and the melt here solidifies earlier than swelling 2. At t 0 + 15 ms, molten metal accumulates in front of valley 2 and forms a new swelling (swelling 3). Swellings and valleys appear alternately, and then, humping forms after the melt pool solidifies.  First the cause of the high-speed melt flow is analyzed. The velocity fields on the longitudinal section of the melt pool at different welding speeds are compared in Figure  9. It can be observed that at welding speed of 24 m/min the melt flow has higher velocity, and the maximum velocity can reach 6 m/s or higher. In addition, shapes of the melt pool are much different between these two simulation cases. At a welding speed of 16 m/min, the melt pool is short and thick, but with a welding speed of 24 m/min, the melt pool is long and thin. The differences of melt flow velocity and melt pool shape can be attributed First the cause of the high-speed melt flow is analyzed. The velocity fields on the longitudinal section of the melt pool at different welding speeds are compared in Figure 9. It can be observed that at welding speed of 24 m/min the melt flow has higher velocity, and the maximum velocity can reach 6 m/s or higher. In addition, shapes of the melt pool are much different between these two simulation cases. At a welding speed of 16 m/min, the melt pool is short and thick, but with a welding speed of 24 m/min, the melt pool is long and thin. The differences of melt flow velocity and melt pool shape can be attributed to two reasons. The first reason is the different recoil pressure. The recoil pressure is the main driving force of the melt flow and is influenced by the temperature of the keyhole front wall as Equation (14) shows. The maximum temperatures of the keyhole front wall in these two simulation cases are 3547 K and 3805 K for welding speeds of 16 m/min and 24 m/min, respectively. With the same heat input, the keyhole front wall has less cooling time under higher welding speed, and thus, it has higher temperature. Therefore, it has larger recoil pressure, and therefore, the melt flow has a higher initial velocity when leaving the keyhole front wall. The other reason is the tilt angle α of the melt pool rear boundary. As shown in Figure 10, α can be calculated as where v s is the moving speed of the melt pool rear boundary, and v w is the welding speed. When the welding speed is high, α is small. That is to say that the melt pool rear boundary has a gentle slope, which can also be observed in Figure 9. When the melt flow reaches the melt pool rear boundary, a smaller α can help it to retain a high x-direction velocity component and flow rearward, and thus, the high-speed melt flow can be maintained.  (14) shows. The maximum temperatures of the keyhole front wall in these two simulation cases are 3547 K and 3805 K for welding speeds of 16 m/min and 24 m/min, respectively. With the same heat input, the keyhole front wall has less cooling time under higher welding speed, and thus, it has higher temperature. Therefore, it has larger recoil pressure, and therefore, the melt flow has a higher initial velocity when leaving the keyhole front wall. The other reason is the tilt angle α of the melt pool rear boundary. As shown in Figure 10, α can be calculated as where vs is the moving speed of the melt pool rear boundary, and vw is the welding speed. When the welding speed is high, α is small. That is to say that the melt pool rear boundary has a gentle slope, which can also be observed in Figure 9. When the melt flow reaches the melt pool rear boundary, a smaller α can help it to retain a high x-direction velocity component and flow rearward, and thus, the high-speed melt flow can be maintained.  Next, the cause of the alternate appearance of swellings and valleys is analyzed. In Figure 8, at the beginning of the welding process, the molten metal flows out from the  (14) shows. The maximum temperatures of the keyhole front wall in these two simulation cases are 3547 K and 3805 K for welding speeds of 16 m/min and 24 m/min, respectively. With the same heat input, the keyhole front wall has less cooling time under higher welding speed, and thus, it has higher temperature. Therefore, it has larger recoil pressure, and therefore, the melt flow has a higher initial velocity when leaving the keyhole front wall. The other reason is the tilt angle α of the melt pool rear boundary. As shown in Figure 10, α can be calculated as where vs is the moving speed of the melt pool rear boundary, and vw is the welding speed. When the welding speed is high, α is small. That is to say that the melt pool rear boundary has a gentle slope, which can also be observed in Figure 9. When the melt flow reaches the melt pool rear boundary, a smaller α can help it to retain a high x-direction velocity component and flow rearward, and thus, the high-speed melt flow can be maintained.  Next, the cause of the alternate appearance of swellings and valleys is analyzed. In Figure 8, at the beginning of the welding process, the molten metal flows out from the melt pool to the surface of the workpiece and then cools down and solidifies rapidly. Dur- Next, the cause of the alternate appearance of swellings and valleys is analyzed. In Figure 8, at the beginning of the welding process, the molten metal flows out from the melt pool to the surface of the workpiece and then cools down and solidifies rapidly. During the solidification, the viscosity of the molten metal increases sharply [23], so the molten metal accumulates and swelling 1 forms. At this time, the melt pool has just formed, so the formation process of swelling 1 lacks generality. Therefore the formation of the subsequent swellings and valleys is the main focus below. In order to analyze the flow of the melt pool quantitatively, for every swelling except swelling 1 in Figure 8, three measuring cross sections are picked in the melt pool at the moment 1 ms before the swelling starts to form, and the flowrates along the minus x direction on these cross sections are measured. For these three cross sections, one is 1.5 mm behind the laser beam center and denoted as cross section A, one is at the tail of the swelling when it starts to form and denoted as cross section C, and one is 1 mm in front of cross section C and denoted as cross section B. For example, as shown in Figure 11, for swelling 2, t = 16 ms is regarded as the moment it starts to form, so the flowrates on cross sections A, B, and C at t = 15 ms are measured. Then, the changes of flowrate per unit length ∆q/∆x between two adjacent cross sections are calculated. The positions of measuring cross sections and the measurement results of flowrates for different swellings are listed in Table 4. It can be noticed that at the moment 1 ms before swellings start to form, the absolute value of ∆q/∆x between B and C is much larger than that between A and B. The average value of ∆q/∆x is −15.6 mm 2 /s between B and C, while −5.8 mm 2 /s between A and B. This indicates that the melt flow decelerates more sharply between B and C, leading to the accumulation of molten metal, so the swelling starts to form. The sharp deceleration of the melt flow is caused by the solidification of the melt pool. As shown in Figure 11, before the swelling starts to form, the melt flow at C is thin, so the heat capacity is small, and the melt pool solidifies rapidly here. The subsequent melt flow is hindered by the sharply increasing viscosity and decelerates when reaching the solidified region near C, so molten metal accumulates, and the swelling starts to form here. This process also agrees with the observed process in the experiment, in which a valley appears and solidifies rapidly in the melt pool and then a swelling appears in front of the valley. ing the solidification, the viscosity of the molten metal increases sharply [23], so the molten metal accumulates and swelling 1 forms. At this time, the melt pool has just formed, so the formation process of swelling 1 lacks generality. Therefore the formation of the subsequent swellings and valleys is the main focus below. In order to analyze the flow of the melt pool quantitatively, for every swelling except swelling 1 in Figure 8, three measuring cross sections are picked in the melt pool at the moment 1 ms before the swelling starts to form, and the flowrates along the minus x direction on these cross sections are measured. For these three cross sections, one is 1.5 mm behind the laser beam center and denoted as cross section A, one is at the tail of the swelling when it starts to form and denoted as cross section C, and one is 1 mm in front of cross section C and denoted as cross section B. For example, as shown in Figure 11, for swelling 2, t = 16 ms is regarded as the moment it starts to form, so the flowrates on cross sections A, B, and C at t = 15 ms are measured. Then, the changes of flowrate per unit length ∆ ∆ ⁄ between two adjacent cross sections are calculated. The positions of measuring cross sections and the measurement results of flowrates for different swellings are listed in Table 4. It can be noticed that at the moment 1 ms before swellings start to form, the absolute value of ∆ ∆ ⁄ between B and C is much larger than that between A and B. The average value of ∆ ∆ ⁄ is −15.6 mm 2 /s between B and C, while −5.8 mm 2 /s between A and B. This indicates that the melt flow decelerates more sharply between B and C, leading to the accumulation of molten metal, so the swelling starts to form. The sharp deceleration of the melt flow is caused by the solidification of the melt pool. As shown in Figure 11, before the swelling starts to form, the melt flow at C is thin, so the heat capacity is small, and the melt pool solidifies rapidly here. The subsequent melt flow is hindered by the sharply increasing viscosity and decelerates when reaching the solidified region near C, so molten metal accumulates, and the swelling starts to form here. This process also agrees with the observed process in the experiment, in which a valley appears and solidifies rapidly in the melt pool and then a swelling appears in front of the valley. Figure 11. Illustration of measuring cross sections (taking swelling 2, for example). Figure 11. Illustration of measuring cross sections (taking swelling 2, for example).
When a swelling starts to form, the radius of the melt pool at the swelling increases. Because the melt pool is thin and long, it is unstable under the effect of surface tension, as described in the Rayleigh instability model. The additional pressure due to surface tension decreases at the swelling. Therefore, the molten metal in front of the swelling flows into the swelling further, resulting in the further growing up of the swelling. Meanwhile, the melt flow in front of the swelling narrows, leading to the formation of a valley. Then, the valley solidifies and induces the formation of the next swelling. It should be pointed out that the proposed explanation here has similarity to the traditional explanation of humping formation by using Rayleigh instability model [8][9][10][11][12][13]. The difference is that, in the traditional explanation, the origin of the non-uniform radius of the melt pool is attributed to perturbations, which is vague. In the proposed explanation here, the origin of the non-uniform radius of the melt pool is made clear. The solidified region in the melt pool hinders the subsequent melt flow, leading to the accumulation of molten metal, so the radius of the melt pool becomes non-uniform.

Effect of Surface Tension and Marangoni Force
From the above analysis, it can be noticed that the surface tension plays an important role in the humping formation process, so it is necessary to investigate the effect of surface tension. In addition, the temperature gradient of surface tension produces a Marangoni force, so the effect of the Marangoni force also needs investigation. The welding parameters used here are still 3000 W and 24 m/min. In the simulation model, the surface tension coefficient σ is determined by σ 0 and dσ/dT, as shown in Equation (12). First σ 0 is changed from its original value 1.77 N·m −1 to 1.27 N·m −1 and 2.27 N·m −1 , respectively, while dσ/dT remains unchanged. The temperatures of most regions in the melt pool are below 2000 K except the region near the laser beam, so dσ/dT can be regarded as positive. The longitudinal sections and appearances of the simulated weld seams are shown in Figure 12. Next dσ/dT is fixed to 0. That is to say, σ is set to fixed values. σ is set to 1.26 N·m −1 and 1.76 N·m −1 , respectively, and the longitudinal sections and appearances of the simulated weld seams are shown in Figure 13.
From Figures 12 and 13, it can be found that with the increase of σ, the humping becomes more severe. This is due to the fact that when σ increases, the difference of additional pressure between the swelling and the melt flow in front of it increases, so it is easier for the swelling to grow and for a valley in front of the swelling to form and induce the formation of the next swelling.
In addition, in the simulation case of Figure 5b, when T < 2000 K, it can be calculated that σ < 1.16 N·m −1 . That is to say that σ in Figure 5b is smaller than that in Figure 13a where σ = 1.16 N·m −1 . According to the above analysis, a smaller σ can help to reduce humping. However, it can be found that humping in Figure 5b is more severe than that in Figure 13a. This indicates that a positive dσ/dT can contribute to the formation of humping. When dσ/dT > 0, the Marangoni forces pointing from low temperature region to high temperature region exist in the melt pool. The center of the swelling has a large heat capacity, so the temperature here is higher than in the surrounding region. On the longitudinal section of the melt pool, the Marangoni force points from the front to the swelling, as shown in Figure 14a. On the cross section of the swelling, the Marangoni force points from sides to the center, as shown in Figure 14b. The Marangoni forces on these two sections can both contribute to the growing up of the swelling and thus make humping more severe.     capacity, so the temperature here is higher than in the surrounding region. On the longitudinal section of the melt pool, the Marangoni force points from the front to the swelling, as shown in Figure 14a. On the cross section of the swelling, the Marangoni force points from sides to the center, as shown in Figure 14b. The Marangoni forces on these two sections can both contribute to the growing up of the swelling and thus make humping more severe.

Effect of Viscous Force
The viscous force can influence the deceleration of the melt flow, so it shall also influence the formation of humping. In order to investigate the effect of the viscous force, the viscosity of the liquid metal is changed from its original value of 0.006 to 0 and 0.02, respectively. The longitudinal sections and appearances of the simulated weld seams are shown in Figure 15. With the increase of the viscosity, the appearance of the weld seam becomes flat, and humping becomes less severe. According to the previous analysis of the humping formation process, the high-speed melt flow is hindered by the solidified region in the melt pool, leading to the accumulation of molten metal and the formation of the swelling, which is the origin of the humping formation process. With large viscosity, the melt flow has decelerated significantly under the effect of viscous force before reaching the solidified region in the melt pool, making it harder for the swelling to form, so humping becomes less severe. This result can also be a reference to humping suppression. If some methods can be taken to influence the behavior of the melt pool and make the highspeed melt flow decelerate in advance, then humping can be suppressed.

Effect of Viscous Force
The viscous force can influence the deceleration of the melt flow, so it shall also influence the formation of humping. In order to investigate the effect of the viscous force, the viscosity of the liquid metal µ is changed from its original value of 0.006 to 0 and 0.02, respectively. The longitudinal sections and appearances of the simulated weld seams are shown in Figure 15. With the increase of the viscosity, the appearance of the weld seam becomes flat, and humping becomes less severe. According to the previous analysis of the humping formation process, the high-speed melt flow is hindered by the solidified region in the melt pool, leading to the accumulation of molten metal and the formation of the swelling, which is the origin of the humping formation process. With large viscosity, the melt flow has decelerated significantly under the effect of viscous force before reaching the solidified region in the melt pool, making it harder for the swelling to form, so humping becomes less severe. This result can also be a reference to humping suppression. If some methods can be taken to influence the behavior of the melt pool and make the high-speed melt flow decelerate in advance, then humping can be suppressed.

Humping Suppression by TIG Arc
According to the above analysis, a humping suppression method is proposed by adding a TIG (Tungsten Inert Gas) arc behind the laser beam. An experiment is conducted to investigate the influence of adding a TIG arc behind the laser beam on humping formation

Humping Suppression by TIG Arc
According to the above analysis, a humping suppression method is proposed by adding a TIG (Tungsten Inert Gas) arc behind the laser beam. An experiment is conducted to investigate the influence of adding a TIG arc behind the laser beam on humping formation in high-speed laser welding. The experimental setup described in Section 2 is used to perform laser welding and high-speed imaging. The laser power is 3000 W, and the welding speed is 24 m/min during the experiment. In addition, a Fronius TransTig 4000 TIG power source is used to output the TIG arc, and the direct current electrode negative (DCEN) setup is used. The arc current is 100 A. The TIG torch is fixed on the laser processing head, and the relative position relationship between the torch and the laser beam is shown in Figure 16. The flowrate of the shielding gas Argon delivered out of the TIG nozzle is 8 L/min. The active light source is not used during this experiment since its light can be blocked by the TIG torch. Appearances of the weld seam and melt pool with TIG arc added are shown in 17. No humping occurs on the weld seam, indicating that humping in high-spee welding can be suppressed by adding a TIG arc behind the laser beam. The arc an beam act on a common melt pool, and the melt pool becomes wider after going t the arc, and the surface of the melt pool is flat. Due to this, the arc force and arc h influence the behavior of the melt pool. Because the distance between the laser bea the TIG electrode is relatively large, the coupling of the arc with the laser-induced plume does not occur. The laser irradiation region does not get heat from the arc melt flow is not strengthened by the arc, which is different from the typical laserbrid welding [24,25]. Due to the inclination angle of the TIG torch, the arc force b the resistance of the high-speed melt flow, and therefore, the melt flow decelerates icantly before reaching the solidified region in the melt pool, as shown in Figure  addition, under the effect of the arc heat, the melt pool becomes wider after going t the arc, as shown in Figure 18b. The increase of the melt pool width can make th flow decelerate further in case of flow conservation and can also reduce the instab the melt flow. Therefore, when the melt flow reaches the solidified region of the me its sharp deceleration is avoided, so the swelling does not occur. As a result, hump be suppressed. Appearances of the weld seam and melt pool with TIG arc added are shown in Figure 17. No humping occurs on the weld seam, indicating that humping in high-speed laser welding can be suppressed by adding a TIG arc behind the laser beam. The arc and laser beam act on a common melt pool, and the melt pool becomes wider after going through the arc, and the surface of the melt pool is flat. Due to this, the arc force and arc heat can influence the behavior of the melt pool. Because the distance between the laser beam and the TIG electrode is relatively large, the coupling of the arc with the laserinduced vapor plume does not occur. The laser irradiation region does not get heat from the arc, so the melt flow is not strengthened by the arc, which is different from the typical laser-arc hybrid welding [24,25]. Due to the inclination angle of the TIG torch, the arc force becomes the resistance of the high-speed melt flow, and therefore, the melt flow decelerates significantly before reaching the solidified region in the melt pool, as shown in Figure 18a. In addition, under the effect of the arc heat, the melt pool becomes wider after going through the arc, as shown in Figure 18b. The increase of the melt pool width can make the melt flow decelerate further in case of flow conservation and can also reduce the instability of the melt flow. Therefore, when the melt flow reaches the solidified region of the melt pool, its sharp deceleration is avoided, so the swelling does not occur. As a result, humping can be suppressed.
addition, under the effect of the arc heat, the melt pool becomes wider after going the arc, as shown in Figure 18b. The increase of the melt pool width can make flow decelerate further in case of flow conservation and can also reduce the insta the melt flow. Therefore, when the melt flow reaches the solidified region of the m its sharp deceleration is avoided, so the swelling does not occur. As a result, hum be suppressed.

Conclusions
In this research, a detailed study was conducted to investigate the humping formation and suppression in high-speed laser welding through high-speed imaging and numerical modeling. The main conclusions are as follows: (1) At high welding speed, the high-speed rearward melt flow caused by recoil pressure is hindered by the solidified region and accumulates, causing the formation of swelling. Then, under the effect of surface tension, the swelling grows. This process repeats leading to humping. (2) The increase in surface tension and Marangoni force can lead more molten metal to flow into the swelling and promote the humping formation process.

Conclusions
In this research, a detailed study was conducted to investigate the humping formation and suppression in high-speed laser welding through high-speed imaging and numerical modeling. The main conclusions are as follows: (1) At high welding speed, the high-speed rearward melt flow caused by recoil pressure is hindered by the solidified region and accumulates, causing the formation of swelling.
Then, under the effect of surface tension, the swelling grows. This process repeats leading to humping. (2) The increase in surface tension and Marangoni force can lead more molten metal to flow into the swelling and promote the humping formation process.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.