Impact of Wavelengths of CO₂, Disk, and Green Lasers on Fusion Zone Shape in Laser Welding of Steel

Article information

J Weld Join. 2020;38(3):235-240

doi : https://doi.org/10.5781/JWJ.2020.38.3.1

Won-Ik Cho ^*^,

, Suck-Joo Na ^**^,^***

* BIAS - Bremer Institut für angewandte Strahltechnik GmbH, 28359 Bremen, Germany

** Dept. of Mechanical Engineering, KAIST, Daejeon, 34141, Korea

*** School of Materials Science and Engineering, Xi’an Jiaotong University, Xi’an, 710049, China

^† Corresponding author : cho@bias.de

Received 2020 May 19; Revised 2020 June 10; Accepted 2020 June 11.

Abstract

Numerical simulations of the multi-kW laser welding of structural carbon steel were conducted to analyze the effects of laser wavelength on fusion zone shape. For the heat and mass transfer analysis, the basic governing equations were utilized in addition to the volume-of-fluid (VOF) conservation equation to calculate the free surface boundary. Laser welding models that considered phenomena such as multiple reflections of the laser beam at the keyhole, absorption, recoil pressure, and shear stress and heat source by metal vapor generated were used as well. Two models: a general ands a simplified model, which were both based on the Fresnel equations, were considered with the wavelengths of CO₂, disk, and green lasers for calculating of the laser absorptivity of the material. The difference found in the results of the general model and the simplified model in the range of the wavelengths considered, was negligible. In addition, the results obtained for disk laser welding and green laser welding did not exhibit any considerable difference in terms of their impact on the fusion zone shape.

Keywords: Laser welding; Numerical simulation; Laser wavelength; Fresnel reflection model

1. Introduction

Laser welding has the advantage of achieving deep penetration and high welding speed, leaving a relatively low heat affected zone in the weld thanks to its high power density. There are CO₂ lasers and lamp pumped Nd:YAG lasers that are traditionally widely used in the industry. Of them, the CO₂ laser has been widely used in thick plate applications which requires high power, yet has a disadvantage that it cannot be transmitted by optical fiber due to its 10.6 μm wavelength. In contrast, the lamp pumped Nd:YAG laser has been applied to various fields because it can be transmitted by optical fiber. However, it has some shortcomings in that not only its power is limited but the efficiency and quality are deteriorated due to the temperature rise of the laser gain medium, thermal lensing, and depolarization loss. Currently, the application of disk and fiber lasers that overcome these shortcomings is increasing. Since they are able to maintain high efficiency and quality even at high power, they have been widely applied not only in the micro field, but in the thick plate applications mainly covered by the CO₂ laser. Since these lasers operate at a wavelength of approximately 1.03 μm or 1.07 μm, they can be transmitted of course by optical fiber as with the Nd:YAG laser. With the recent increase in the demand for e-mobility, the demand for welding of materials with high electrical conductivity such as copper has also grown. However, these materials have a very low absorptivity for traditional lasers with an infrared wavelength. Thus, application of lasers with a short wavelength (e.g. a green laser with a wavelength of 532 nm) is being examined.

Experimental studies account for the majority of studies on laser welding, but numerical studies have also been actively underway thanks to the advancement of computing technology. Under the assumption that the model used is verified, numerical studies have the advantage of closely observing areas that are difficult to reach because of experimental constraints. They are of great help in understanding the physical phenomenon of welding. In particular, since the molten pool is formed in fusion welding such as laser welding, molten pool analysis is widely conducted for realistic simulation. With reference to molten pool analysis for laser deep welding, it is crucial to realize the phenomenon that the keyhole is generated by the vaporization of the material due to the high power density (over 10⁶ W/cm²). The simulation of the early laser welding mainly focused on the heat transfer in conduction mode, thus was performed similarly to gas tungsten arc (GTA) welding simulation1). Kaplan predicted the 2-dimensional keyhole shape using the energy balance equation considering the multiple reflections and plasma in the laser keyhole. The energy absorbed by multiple reflections was calculated by assuming a conical keyhole. The result showed that the increase in the energy absorbed by multiple reflections had a significant effect on the keyhole shape2). The 3-dimensional keyhole shape was also predicted in a similar way3). Because it was still hard to predict accurate keyhole shapes, the observed keyhole shape using an X-ray transmission method was also used in laser welding simulation4). A study using a predicted keyhole shape to analyze the behavior of the molten pool was also conducted5). With the advancement of computing technology and algorithms, numerical simulation of the keyhole shape and the molten pool behavior began to be integrated by applying the level set method and VOF method, which are free surface tracking algorithms6,7). Because of the mass conservation error of the level set method and the complex algorithm that attempts to solve it, the VOF-based method is currently widely employed. By extracting the keyhole shape in real-time using the VOF method and applying the ray-tracing technique to the shape, the multiple reflection phenomenon of the laser beam in the keyhole was modeled8). Simulation of the laser drilling process taking into account the effect of laser polarization was also performed9). As the rapid keyhole change can be calculated in detail in this way, it is possible to simulate the occurrence of pores arising from the collapse of the laser keyhole with a bubble model10). It was also simulated how the alloy elements supplied from a filler metal due to the strong flow around the keyhole were distributed on the weld11). The path in which the laser beam was focused and diffracted through the laser optics was taken into account, and the effect of the metal vapor generated by the laser was modeled using a shear stress and a heat source. It enabled the prediction of relatively precise results even in a high power laser welding12). Recently, a numerical study considering the oscillation of the laser beam was also conducted, and the effect of oscillation frequency on the molten pool was analyzed13). In addition, the buttonhole phenomenon arising from the effect of surface tension pressure in laser welding of thin plates was also simulated14).

In this study, the 3-dimensional computational fluid dynamics (CFD) simulation of laser welding using multi-kW CO₂, disk, and virtual green lasers was performed. The focused and diffracted beam path was used, and the multiple reflections of the beam were considered. The material absorptivity of the laser beam was calculated by the two types of Fresnel reflection models to investigate the effects of laser wavelengths. Buoyancy, Marangoni force, recoil pressure, shear stress and heat source by metal vapor, and bubble model was applied. The final fusion zone shape was used to compare the results.

2. Methodology

Fig. 1 shows the computational domain and mesh for the heat and mass transfer analysis of laser welding. Of the entire domain (56 mm long, 26 mm wide, 15.2 mm high), a fine mesh for the central analysis domain (36 mm long, 6 mm wide, and 15.2 mm high) where the molten pool was formed, and a coarse mesh for the other domain was used. The size of the fine mesh was chosen to be 0.2 mm through a mesh test. The boundary condition on the six sides of the computational domain was defined for the partial derivative value to be 0 (Neumann boundary condition) for all independent variables. In order to perform 3-dimensional transient CFD simulation, the mass conservation equation, momentum equation, and energy conservation equation, which are basic governing equations, were used. Additionally, the VOF equation for free surface tracking was used. Newtonian fluid and incompressible laminar flow were assumed for molten metal. For a laser welding model, the beam path considering the beam focusing phenomenon and diffraction phenomenon by the laser optics was modeled, and multiple reflection phenomenon was simulated using the ray-tracing method. The material absorptivity of the laser beam at the reflected position was determined by the Fresnel reflection model. Along with this, the recoil pressure, Marangoni effect, buoyancy, and shear stress and heat source by metal vapor were considered, and an adiabatic bubble pressure model was applied to investigate the behavior of bubbles generated inside the molten pool caused by the collapse of the keyhole. The literature was referred for the descriptions and formulas for these models and the physical properties of structural carbon steel (A36 (ASTM) or SS400 (JIS)) used for simulation12).

Fig. 1

Schematic of the computational domain and mesh for laser welding simulation.

Two Fresnel models based on the Fresnel equations were used to perform numerical analysis considering different wavelengths of CO₂, disk, and green lasers. Fig. 2 shows the absorptivity data for the laser incidence angle of the general model (based on the Drude theory) and the simplified model (simplified using the Hagen- Rubens relation) (red: general model, blue: simplified model). Table 1 shows the laser and optics parameters used in the simulation. It was assumed that the virtual parameter values for the green laser were the same as the disk laser except for the wavelength and M². In all simulation cases, the welding speed was fixed at 1 m/min.

Fig. 2

Absorptivity vs. incident angle based on general and simplified fresnel models for different lasers

Table 1

Laser and optics parameters used

3. Results and discussion

Fig. 3 shows the mesh test results. It was found that the calculated penetration depth increased with the decrease in the mesh size, and gradually converged to the keyhole depth of the experimental results. These results are considered to be because in reality the laser beam can reach the bottom of the keyhole as the finer mesh is used. In the simulation, there should be at least an empty space larger than the mesh size in order for the laser beam to pass. Based on the test results, accuracy of simulation can be secured when a mesh size of 0.2 mm or less is used. However, given the calculation time spent, when reducing the mesh size from 0.2 mm to 0.15 mm, it took approximately 2.8 times more time. Therefore, in this study, a mesh size of 0.2 mm having both accuracy and efficiency was used. The mesh size used in the previous studies on the molten pool behavior of laser welding was different depending on the keyhole depth of the target process. 0.25 mm was used as an optimal value for the shallow keyhole15), and 0.2 mm for the deep keyhole16). This is consistent with the mesh test results of this study. Next, in order to investigate if there was an effect depending on the size of the computational domain, simulation was performed by increasing the width of the computational domain from 26 mm to 46 mm as seen in Fig. 4. The results showed that the predicted differences in fusion zone shapes were negligible, thus the computational domain used in this study was sufficient. As for the computational domain of a width of 26 mm used in this study, the temperature at the boundary showed an increase of approximately 14°C from the initial temperature during the simulation.

Fig. 3

Mesh test in terms of accuracy and efficiency of calculation

Fig. 4

Fusion zone shapes calculated in different computational domains

Lastly, the simulation results of the CO₂, disk, and virtual green lasers are shown in Fig. 5. It shows the fusion zone shapes of the YZ cross-section. All of the keyholes grow sufficiently after the start of welding, indicating the result in the steady state. In our previous study, the CO₂ laser and disk laser welding models have already been verified11,12). They are considered to be effective even in the range of welding parameters used in this study. The results of applying the two absorption models shown in Fig. 2 were compared. When using the simplified model instead of the general model, the absorptivity has a larger value. In the CO₂ laser wavelength, the two models have almost the same value, yet the difference between them increases with the decrease in wavelength. Despite the difference in the absorptivity between these models, the prediction results of the fusion zone using the two models did not show much difference. In addition, when comparing the results of the disk laser and the green laser, the difference arising from wavelength changes was also negligible. These results are considered to be because the effect of the models and the wavelengths decreases with the increase in the total absorption (over 80% in this study12), up to 90% is reported17)) by the multiple reflection phenomenon for keyhole welding. After the keyhole grew sufficiently, there was no significant change in the fusion zone shapes in the XZ cross-section as well as the YZ cross-section of Fig. 5. These results can also be found in the previous study using the disk laser12). However, in the regions where the welding starts and where the welding ends, the depth of the keyhole decreases, thus the difference in the energy absorbed is expected to increase. The difference is also expected to influence the properties of the final weld by affecting the molten pool behavior and the temperature distribution. Especially, since more welding defects are observed in those regions, there is a need for further research in the future.

Fig. 5

Fusion zone shapes calculated based on general and simplified Fresnel models for different lasers

4. Conclusion

In this study, the numerical simulation considering the wavelengths of CO₂, disk, and green lasers was performed using laser welding models, and the following results were derived.

1) The mesh test showed that the size of 0.2 mm had the accuracy of result and the efficiency of calculation in the range of welding parameters used in this study. Therefore, this mesh size can be recommended for similar laser welding simulations.
2) In terms of the final fusion zone shape, the general model and the simplified model did not show any significant difference. In addition, when comparing the analysis results of the disk laser and the green laser, there was no significant difference even at different wavelengths. This is considered to be because the effect caused by the difference in the models and the wavelengths is negligible due to the significant increase in the total absorption by multiple reflections (80-90%) for the keyhole welding of this study.

Acknowledgement

The authors received no financial support for the research, authorship, and/or publication of this article. This paper contains a part of the contents in the author´s Ph.D. dissertation18).

ORCID

ORCID: Won-Ik Cho: https://orcid.org/0000-0002-4047-4230

ORCID: Suck-Joo Na: https://orcid.org/0000-0002-1832-3997

References

1. Zacharia T, David S. A, Vitek J. M, Debroy T. Weld Pool Development During GTA and Laser Beam Welding of Type 304 Stainless Steel, Part I- Theoretical Analysis. Weld. J 68(12)1989;:499s–509s.

2. Kaplan A. A Model of Deep Penetration Laser Welding Based on Calculation of the Keyhole Profile. J. Phys. D:Appl. Phys 27(9)1994;:1805–1814. https://doi.org/10.1088/0022-3727/27/9/002.

3. Solana P, Ocaña J. L. A Mathematical Model for Penetration Laser Welding as a Free-Boundary Problem. J. Phys. D:Appl. Phys 30(9)1997;:1300–1313. https://doi.org/10.1088/0022-3727/30/9/005.

4. Kaplan A. F. H, Mizutani M, Katayama S, Matsunawa A. Unbounded Collapse Keyhole and Bubble Formation During Pulsed Laser Interaction with Liquid, Zinc. J. Phys. D:Appl. Phys 352002;:1218–1228. https://doi.org/10.1088/0022-3727/35/11/319.

5. Mahrle A, Schmidt J. The Influence of Fluid Flow Phenomena on the Laser Beam Welding Process. Int. J. Heat Fluid Fl 232002;:288–297. https://doi.org/10.1016/S0142-727X(02)00176-5.

6. Ki H, Mazumder J, Mohanty P. S. Modeling of laser keyhole welding Part I, mathematical modeling, numerical methodology, role of recoil pressure, multiple reflections, and free surface evolution. Metall. Mater. Trans. A 33(6)2002;:1817–1830. https://doi.org/10.1007/s11661-002-0190-6.

7. Lee J. Y, Ko S. H, Farson D. F, Yoo C. D. Mechanism of keyhole formation and stability in stationary laser welding. Phys J. D:Appl. Phys 35(13)2002;:1570–1576. https://doi.org/10.1088/0022-3727/35/13/320.

8. Cho J. H, Na S. J. Implementation of Real-Time Multiple Reflection and Fresnel Absorption of Laser Beam in Keyhole. J. Phys. D:Appl. Phys 392006;:5372–5378. https://doi.org/10.1088/0022-3727/39/24/039.

9. Cho J. H, Na S. J. Theoretical Analysis of Keyhole Dynamics in Laser Drilling Considering the Polarization of Laser. J. Phys. D:Appl. Phys 402007;:7638–7647. https://doi.org/10.1088/0022-3727/40/24/007.

10. Cho W. I, Cho J. H, Cho M. H, Lee J. B, Na S. J. Numerical Simulation Bubble of and Pore Generation by Molten Metal Flow in Laser-GMA Hybrid, Welding. J. Weld. Join 26(6)2008;:67–73. https://doi.org/10.5781/KWJS.2008.26.6.067.

11. Cho W. I, Na S. J, Cho M. H, Lee J. S. Numerical Study of Alloying Element Distribution in CO₂, Hybrid Laser- GMA Welding. Comput. Mater. Sci 49(4)2010;:792–800. https://doi.org/10.1016/j.commatsci.2010.06.025.

12. Cho W. I, Na S. J, Thomy C, Vollertsen F. Numerical simulation of molten pool dynamics in high power disk laser welding. J. Mater. Process. Technol 2122012;:262–275. https://doi.org/10.1016/j.jmatprotec.2011.09.011.

13. Cho W. I, Schultz V, Woizeschke P. Numerical study of the effect of the oscillation frequency in buttonhole welding. J. Mater. Process. Technol 2612018;:202–212. https://doi.org/10.1016/j.jmatprotec.2018.05.024.

14. Cho W. I, Woizeschke P. Analysis of molten pool behavior with buttonhole formation in laser keyhole welding of sheet metal. Int. J. Heat Mass Trans 1522020;:119528. https://doi.org/10.1016/j.ijheatmasstransfer.2020.119528.

15. Cho M. H, Farson D. F. Simulation Study of a Hybrid Process for the Prevention of Weld Bead, Formation Hump. Weld. J 86(9)2007;:253s–262s.

16. Cho J. H, Na S. J. Three-Dimensional Analysis of Molten Pool in GMA-Laser Hybrid, Welding. Weld. J 88(2)2009;:35s–43s.

17. Kawahito Y, Matsumoto N, Abe Y, Katayama S. Laser absorption characteristics in high-power fibre laser welding of stainless steel. Weld. Int 27(2)2013;:129–135. https://doi.org/10.1080/09507116.2011.606151.

18. Cho W. I. Ph. D. dissertation, A Numerical Study on Molten Pool Behavior in Laser-Arc Hybrid Welding Process. KAIST Daejeon, Korea 2012;:48–111.

Article information Continued

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.