Open Access

Solidification cracking in laser beam welding (LBW) is governed by local thermal conditions at the weld-pool boundary, yet the sub-micron microstructural response across the full mushy zone remains poorly understood. This study employs an automated phase-field (PF) simulation workflow, integrated within the Kadi4Mat research data management platform, to conduct a systematic, FAIR-compliant parametric study of dendritic solidification in the quaternary EN 1.4301 (Fe–Cr–Ni–C) alloy. Thermal conditions—thermal gradient (𝐺), solidification velocity (𝑉𝑠), and grain misorientation angle (𝜃𝑅)—are extracted from a thermocouple-validated ANSYS Fluent weld-pool model and used as inputs to 2D and 3D PACE3D phase-field simulations spanning 𝐺 from 100 to 900 K/mm, 𝑉𝑠 from 5 to 40 mm/s, and 𝜃𝑅 from 0◦ to 45◦. Key findings demonstrate that 𝜃𝑅 is a critical parameter—alongside 𝐺 and 𝑉𝑠—governing the cellular-to-dendritic morphological transition, secondary den drite arm formation, and the topology of inter-dendritic liquid (continuous films versus isolated pockets), each carrying distinct solidification cracking risk pathways. The Kadi4Mat workflow reduces manual pre-processing effort substantially while ensuring data reproducibility and reuse. Quantitative validation against electron probe micro-analysis confirms primary dendrite arm spacing predictions within 12% at the upper weld surface and within 4% at the mid-section.

This study presents mechanisms of the evolution of a narrowed region in the weld pool center during deep penetration laser beam welding. In numerous numerical studies presented in this study, it was also found that the local reduction of the weld pool size can cause detrimental effects on the melt flow behavior and the resulting properties of the welds. A particularly large influence of this effect was identified in three aspects. Firstly, the local variation of the solidification sequence of the weld pool causes an increase in the hot-cracking susceptibility due to a locally delayed solidification. Secondly, it was proven that a change in the local length and width of the weld pool is associated with an adverse impact on the potential flow routes of the molten material that induces stronger local variations of its solidification. Thus, the element mixing, e.g. during the welding with filler materials, is blocked. This leads to a non-homogeneous chemical composition of the final weld and can cause undesired effects on the final material properties. Finally, another observed effect is related to the reduced ability of process pores to reach the top surface. As this type of porosity is usually produced around the keyhole tip, the change of the fluid flow regime above this area plays a significant role in determining the final path of the pores until the premature solidification in the middle of the weld pool captures them. This study summarizes mainly numerical results which were supported by selected experimental validation results.

The numerical simulations of high power laser keyhole welding at different welding positions are performed by using Volume-Of-Fluid (VOF) method. The main material is SS400. The multi-physics phenomenon is considered using several models, such as the heat flux of Gaussian heat source, the recoil pressure with Clausisus-Clapeyron equation, the Marangoni flow considering temperature gradient, the buoyancy force with Boussinesq approximation, the additional shear stress and heat source due to metallic vapor ejected through keyhole entrance, the bubble formation assumed as adiabatic bubble, and the multiple-reflection by solving proper discriminant, are used. To analyze the fluid flow pattern, the concept of streamline formed by reconstructing the value of the velocity vector is applied. Partial and full penetration cases at different welding positions are considered. The welding position seems to have only a minor influence on bead formation characteristics in both cases. This is probably due to the fact that the recoil pressure has a major influence when compared to other driving forces. The flow characteristics and fluid velocity in weld pool are analyzed to compare the gravity direction effect at different welding positions. It is observed that the clockwise flow pattern is mainly formed by the recoil pressure on the keyhole surface in the case of partial penetration. The laser energy can't maintain the whole weld pool when the weld pool size becomes too large. And then the solidification starts from the middle part of weld pool and a necked weld pool shape is formed. In the full penetration welding, the weld pool flow patterns are affected by the leakage of laser power through the full penetration keyhole and also by surface tension. Furthermore, the numerical simulation of full penetration welding with AISI316L is also performed to analyze the effect of material properties. The weld bead shapes obtained by simulations were compared with the corresponding experimental results to confirm the validity of the process models adopted and the CFD simulation tool.

An accurate and efficient description of the spatial distribution of laser energy is a crucial factor for the modeling of laser material processing, e.g., laser welding, laser cutting, or laser-based additive manufacturing. In this study, a 3D heat transfer and fluid flow model coupled with the volume-of-fluid algorithm for free surface tracking is developed for the simulation of molten pool dynamics in high-power laser beam welding. The underlying laser-material interactions, i.e., the multiple reflections and Fresnel absorption, are considered by a ray-tracing method. Two strategies of free surface reconstruction used in the ray-tracing method are investigated: a typical piecewise linear interface calculation (PLIC)-based method and a novel localized Level-Set method. The PLIC-based method is discrete, resulting in non-continuous free surface reconstruction. In the localized Level-Set method, a continuous free surface is reconstructed, and thus the exact reflection points can be determined. The calculated spatial laser energy distribution and the corresponding molten pool dynamics from the two methods are analyzed and compared. The obtained numerical results are evaluated with experimental measurements to assure the validity of the proposed model. It is found that distinct patterns of the beam multiple reflections are obtained with the different free surface reconstructions, which shows significant influence not only on the molten pool behaviors but also on the localized keyhole dynamics.