Experimental and Numerical Analysis on TIG Arc Welding of Stainless Steel Using RSM Approach

: This study involves the validating of thermal analysis during TIG Arc welding of 1.4418 steel using ﬁnite element analyses (FEA) with experimental approaches. 3D heat transfer simulation of 1.4418 stainless steel TIG arc welding is implemented using ABAQUS software (6.14, ABAQUS Inc., Johnston, RI, USA), based on non-uniform Goldak’s Gaussian heat ﬂux distribution, using additional DFLUX subroutine written in the FORTRAN (Formula Translation). The inﬂuences of the arc current and welding speed on the heat ﬂux density, weld bead geometry, and temperature distribution at the transverse direction are analyzed by response surface methodology (RSM). Validating numerical simulation with experimental dimensions of weld bead geometry consists of width and depth of penetration with an average of 10% deviation has been performed. Results reveal that the suggested numerical model would be appropriate for the TIG arc welding process. According to the results, as the welding speed increases, the residence time of arc shortens correspondingly, bead width and depth of penetration decrease subsequently, whilst simultaneously, the current has the reverse effect. Finally, multi-objective optimization of the process is applied by Derringer’s desirability technique to achieve the proper weld. The optimum condition is obtained with 2.7 mm/s scanning speed and 120 A current to achieve full penetration weld with minimum fusion zone (FZ) and heat-affected zone (HAZ) width.


Introduction
Welding is a method of joining two similar or dissimilar metals with or without filler material [1]. One of the well-known conventional Arc welding processes is Tungsten Inert Gas welding known as TIG welding. This method is used due to having powerful controllable features to transmit heat to the welding line locally [2]. TIG fusion zone size is narrower than other single-pass arc welding processes comparatively [3]. In the TIG arc welding process, a non-consumable tungsten electrode is used to create an arc on the workpiece. This welding method was applied to the metals with inert gas as the shielding weld pool. Experiences have shown TIG welding method is a reliable welding process with high-quality welds. Spatters and fumes rarely occur in this process. Low productivity of this method causes to restrict in the limited industries [4]. One of the current 2.1. Initialization 2.1. 1

. Element Type
In the first step, a solid deformable body model of the workpiece is offered for heat transfer welding simulation. To reduce simulation time and improve efficiency, half of the workpiece is modeled in the part module. For the current model, the workpiece dimension was 127 mm × 38 mm × 4 mm. To compute the temperature of each node, the whole model needs to be split into tiny elements. As explained in the mesh section, an eight-node heat transfer element (DC3D8) is defined for the 80% model, and a four-node linear tetrahedron heat transfer element (DC3D4) is used as intermediate elements for thermal analysis [16].

Material Properties
In this part, due to extreme variation in temperature during TIG welding, the temperature-dependent material properties should be recorded in the FEM simulation [16]. These properties of 1.4418 stainless steel are defined to be homogeneous during heat transfer analysis. As shown in Table 1, the temperature-dependent physical properties, including specific heat and thermal conductivity, are given. Due to negligible variation in density, this factor is supposed to be constant at 7700 Kg/m 3 . Table 1. Temperature-dependent material properties [17]. The accuracy and convergence of the numerical simulation results depend on the element type and the dimensions of various regions significantly. Figure 1 represents the three-dimensional geometrical model with identical mesh used in this research. Higher thermal gradient in the fusion zone and heat-affected zone and the need for more accurate temperature computation in these regions, a finer mesh is considered [16]. This subject causes us to model part geometry with two types of elements that differ in the number of nodes. The partitioning technique is utilized to separate different welding regions whose differences in experienced temperature. In thermal analysis, because of the temperature degree of freedom, DC3D8 and DC3D4 elements are used. This DC3D8 element is an eight-node linear heat transfer brick element. Therefore, it is necessary to mention that intermediate elements are utilized to connect coarser elements to more delicate elements. The DC3D4 is defined as an intermediate element. This is a four-node tetrahedron linear heat transfer element. The mesh is composed of 22,202 3D heat transfer elements. The minimum size of each element is 0.25 mm in the fusion zone (FZ). differences in experienced temperature. In thermal analysis, because of the temperature degree of freedom, DC3D8 and DC3D4 elements are used. This DC3D8 element is an eight-node linear heat transfer brick element. Therefore, it is necessary to mention that intermediate elements are utilized to connect coarser elements to more delicate elements. The DC3D4 is defined as an intermediate element. This is a four-node tetrahedron linear heat transfer element. The mesh is composed of 22,202 3D heat transfer elements. The minimum size of each element is 0.25 mm in the fusion zone (FZ).

Heat Transfer Analysis
Thermal analysis during TIG welding and evaluating temperature distribution along the model is applied using the Energy conservation equation. In the transient thermal analysis of the welding process, the origin of the cartesian coordinate system is supposed to locate on the model surface; Y-axis and Z-axis are defined along the welding direction and the thickness direction, respectively. The governing equation for the 3D transient heat transfer is given by Equation (1) [18]: T is the temperature, p is the material density, c is the specific heat, k is the material conductivity, and Q is the internal heat generation per unit volume. Moreover, x, y, and z are the coordinates system attached to the heat source.

Heat Source Model
The accurate prediction of the fusion zone geometry and temperature distribution has a close relationship with the utilized heat source model and its parameters. In the welding process, defining the moving heat source with precise volumetric heat flux along the thickness is vital. Previous results represent that the double ellipsoid shape provides the best distribution, especially for arc welding joints; Since the double-ellipsoidal heat power density model presented by Goldak et al., is implemented in this study to simulate the TIG arc welding process [19]. Gaussian distribution of heat source is utilized to compute and analyze the thermal distribution; subsequently, dimensions of weld bead such as bead width and depth of penetration are achieved.
In the current simulation work, two rear and front ellipsoidal shapes merged, as shown in Figure 2 include two different ellipses, i.e., one in the front of the heat source and the other in the rear section.

Heat Transfer Analysis
Thermal analysis during TIG welding and evaluating temperature distribution along the model is applied using the Energy conservation equation. In the transient thermal analysis of the welding process, the origin of the cartesian coordinate system is supposed to locate on the model surface; Y-axis and Z-axis are defined along the welding direction and the thickness direction, respectively. The governing equation for the 3D transient heat transfer is given by Equation (1) [18]: T is the temperature, p is the material density, c is the specific heat, k is the material conductivity, and Q is the internal heat generation per unit volume. Moreover, x, y, and z are the coordinates system attached to the heat source.

Heat Source Model
The accurate prediction of the fusion zone geometry and temperature distribution has a close relationship with the utilized heat source model and its parameters. In the welding process, defining the moving heat source with precise volumetric heat flux along the thickness is vital. Previous results represent that the double ellipsoid shape provides the best distribution, especially for arc welding joints; Since the double-ellipsoidal heat power density model presented by Goldak et al., is implemented in this study to simulate the TIG arc welding process [19]. Gaussian distribution of heat source is utilized to compute and analyze the thermal distribution; subsequently, dimensions of weld bead such as bead width and depth of penetration are achieved.
In the current simulation work, two rear and front ellipsoidal shapes merged, as shown in Figure 2 include two different ellipses, i.e., one in the front of the heat source and the other in the rear section.
The moving Cartesian coordinate system (x, y, and z) attached to the heat flux distribution (q) presented in Equations (2) and (3): The power density distribution inside the front quadrant is: Similarly, the power density for the back of the heat source is defined as Equation (3):  The moving Cartesian coordinate system (x, y, and z) attached to the heat flux distribution (q) presented in Equations (2) and (3): The power density distribution inside the front quadrant is: Similarly, the power density for the back of the heat source is defined as Equation For the front ellipse of the heat source, the front heat absorption coefficient defined. Also, for the rear ellipse of the heat source, the rear heat absorption coefficient is defined. Elliptic lengths of the front and rear sections of the heat source are different. is defined as front elliptical length and, defined as rear elliptical length. The values of parameters a, b, c in the front and rear ellipse can differ [16]. The variables of the double ellipsoidal heat source equation used in the FEM are listed in Table 2. V is the scanning speed along the welding direction, and t is the current time. According to Equation (4), Q is the heat source of the TIG welding process and is equal to the multiple arc voltage (V), welding current (I), and process efficiency (η).  To apply a moving volumetric heat source, the DFLUX user subroutine was programmed and utilized. Non-uniform distributed heat flux as a function of time and position was determined using DFLUX as an interface during thermal analysis. For the front ellipse of the heat source, the front heat absorption coefficient f f defined. Also, for the rear ellipse of the heat source, the rear heat absorption coefficient f r is defined. Elliptic lengths of the front and rear sections of the heat source are different. a f is defined as front elliptical length and, a r defined as rear elliptical length. The values of parameters a, b, c in the front and rear ellipse can differ [16]. The variables of the double ellipsoidal heat source equation used in the FEM are listed in Table 2. V is the scanning speed along the welding direction, and t is the current time. According to Equation (4), Q is the heat source of the TIG welding process and is equal to the multiple arc voltage (V), welding current (I), and process efficiency (η). To apply a moving volumetric heat source, the DFLUX user subroutine was programmed and utilized. Non-uniform distributed heat flux as a function of time and position was determined using DFLUX as an interface during thermal analysis.

Boundary Conditions
During thermal analysis simulation of TIG arc welding, the boundary conditions were identified for the model. The initial temperature of the model is considered 25°C. Heat loss through radiation and convection is considered boundary conditions. The top surface of the welding plate is exposed to a heat source, whereas the lateral surfaces during thermal analysis experience heat losses by convection and radiation only. The governing equation used for boundary condition given as Equation (5): where K is isotropic thermal conductivity, q is heat flux (W/m 2 K), h is the convective coefficient, T is the local temperature, T ∞ is the ambient temperature, ε is the emissivity coefficient, and σ is the Stefan-Boltzmann constant which is 5.67 × 10 −8 W/m 2 K 4 for radiation [20]. The free convection coefficient of airflow is considered to be 10 W/m 2 K for stainless steel. It is used for all surfaces not exposed to the shielding gas. The welding line, which is directly below the nozzle due to the shielding gas flow, experiences forced convection heat transfer. Subsequently, the heat transfer caused by the shielding gas flow is considered as forced convection heat transfer in the simulation. The following empirical relation Equation (6) initially suggested by Gardon and Cobonpue for air and later modified by Steen for any gas was used for a part of the top surface under the nozzle.
k gas thermal conductivity of shielding gas (W /m K), NPD Nozzle-to-plate distance (m), V gas Shielding gas flow velocity (m/s), v gas Kinematic viscosity of shielding gas (m 2 /s), D Nozzle diameter (m), ρ gas The density of shielding gas (kg/m 3 ) µ gas Dynamic viscosity (kg/m/s), α gas Thermal diffusivity of shielding gas (m 2 /s) C pgas Specific heat of shielding gas (J/kgK). Re (Reynolds number) and Pr (Prandtl number) based on the properties of the shielding gas. The effect of free and forced convection and radiation as thermal boundary conditions imported to the simulation model as an effective combined temperature-dependent constant (h c ) with Equation (9) using FILM user subroutine employed for surfaces during the welding process. Subroutine FILM calculates (h c ) for every integration point at a time step.

Post-Processing
In the current section, simulation results of the 1.4418 stainless steel TIG welding process, temperature distribution, FZ dimensions consist bead width, and depth of penetration examined. Validation of the numerical model is discussed in detail by comparing the simulation and experimental results in the results section.

Welding Process
In the current study, TIG arc welding of 1.4418 stainless steel plates was carried out. Scanning speed and welding current selected as input process parameters. Also, weld bead width and depth of penetration are examined as responses. It's worth mentioning that primarily experimentally applied a range of the input parameters simulated with FEM, and their results are analyzed. Pure argon was utilized as a shielding gas on the welding zone during the tests. The specimens were manually welded by a direct current electrode negative TIG (DCENTIG) arc with a maximum current of 250 A.

Specimen Preparation
1.4418 stainless steel was used as a material workpiece in the current work. The chemical composition is given in Table 3, which is the average of three X-ray fluorescence (XRF) measurements. The welding sheets were prepared by wire EDM with the dimensions of 127 mm × 38 mm × 4 mm. Surface grinding of the specimens was performed to eliminate the probable recast layer created by the wire EDM process [21].

Metallurgical Characterization
Welded specimens were cut by a wire-EDM process with specified dimensions transversely to examine the weld pool shape. The samples were prepared using standard metallographic techniques, which include grinding and polishing. The welding geometry responses, welding depth of penetration, and weld bead width were measured using Leica MEF 4A optical microscope (Leica Microsystems, Wetzlar, Germany) at a magnification of 50× and investigated by Clemex image software.

Design of Experiments
Design of experiments (DOE) is not a series of equations, but a method utilized to organize and design tests; behind there is a theoretical basis with proper mathematical logic. Thereby, it is an essential technique for the industry and by using this method, more reliable results are obtained. Accordingly, saving more time and resources occurs. Consequently, this method assists in using a minimum number of experiments, in which several input variables varied simultaneously to reach a proper connection between them to evaluate responses. The process of modeling is generally affected by the procedure consider for designing and organizing experiments. RSM is one of these procedures, which is a series of mathematical and statistical techniques that is practical for tracking the relations between input factors and output responses [22]. When all independent variables are measurable and experiment repeated with minor errors, the RSM could be defined using Equation (10): where k is the number of independent input factors. Finding a relationship between independent input factors and the output responses is necessary. To precise the model and responses, a quadratic polynomial is utilized generally as follows Equation (11): where β 0 is the constant value or intercept, β i linear coefficient, β ii is the quadratic coefficient and β ij is the interaction coefficient, and ε is the random error of the developed regression. It's worth mentioning that a proposed regression model will not be a precise approximation of the actual relationship between variables in the space variable, but for a relatively limited area using such models is recommended regularly. In the current study, the experiments were organized based on a central composite design (CCD). Input factors are considered in a five-level RSM design. Welding scanning speed (V) of 1 to 5 mm/s, and welding current (C) of 80 to 120 A chosen as independent input factors. This experimental design contains ten tests: four points in the axial direction, four as factorial points in the corners of the cube, and two points in the center of the cube. Table 4 shows the input factors and experimental design levels demonstrated with coded and actual values. In this study, as mentioned above, the design of experiments was carried out based on the response surface methodology using Minitab 2019 software (Minitab Inc., State College, PA, USA). Given that two levels were assigned to each of the input factors, including the welding current and welding scanning speed, ten experiments were conducted. Weld bead width and weld depth of penetration are considered as responses. The results of the experiments are presented in Table 5. To obtain the mathematical models of responses and optimization, ANOVA is performed using Minitab 2019 software; F-value and adjusted R 2 parameters are used as criteria to fit a regression polynomial that can predict the responses minimizing the regression residuals.     The predicted temperature distributions of welded parts at different increment times of the cooling time after the TIG welding process depicted in Figure 4a    In the following, the subsequent section develops the comparative analysis for FEM simulation, and experimental approaches of weld pool dimensions include weld bead width and depth of penetration for the TIG welding process of 1.4418 stainless steel at different conditions. For the proposed purpose, primarily, three main domains FZ, HAZ, and BM, are identified as resultant regions after the TIG welding process. The Fusion zone In the following, the subsequent section develops the comparative analysis for FEM simulation, and experimental approaches of weld pool dimensions include weld bead width and depth of penetration for the TIG welding process of 1.4418 stainless steel at different conditions. For the proposed purpose, primarily, three main domains FZ, HAZ, and BM, are identified as resultant regions after the TIG welding process. The Fusion zone is evaluated by all nodes whose temperatures are higher than the melting point of the 1.4418 stainless steel. It is possible to estimate the geometry and dimensions of the fusion zone from the simulation considering the applied heat flux. These distances can be measured between nodes in the simulation easily. As stated above, the experimental weld pool geometries were measured by Clemex image analysis software as well. Weld pool geometries (weld bead width and weld depth of penetration of butt welds) between FEM results and experimentally welded samples at ten different cases compared and present in Table 6. Furthermore, their cross-sectional compared images were taken and depicted in Figure 5. According to reported conclusions, the average percent of relative error between experiments and predicted with the model was 9.7% placed in the permissible area. It should be said that outputs revealed a good correlation between the experimental and predicted geometries of the TIG welds at different levels; Eventually, the predicted weld pool geometries verify the proper application of the heat source model.

TIG Welding Experimental Results
In the current section, the weld bead width and weld depth of penetration are exam-

TIG Welding Experimental Results
In the current section, the weld bead width and weld depth of penetration are examined as process parameters. An analysis of variance (ANOVA) was conducted to investigate and specify "significantly effective parameters" on the TIG welding process and evaluate their influences on output responses. These analyses were carried out by full quadratic polynomial functions by using Minitab 2019 statistical software. Statistical data about the weld bead width and weld depth of penetration as a function of variables can be extracted and analyzed with Minitab using reported results illustrated in Table 5.

Weld Bead Width
According to the statistical report of an ANOVA performed on the weld bead width, all main input factors are significant, as shown in Table 7. Among quadratic terms and interactions, none are effective. Equation (12) shows the regression model with coded values for the weld bead width considering significant parameters. The ANOVA, table of weld bead width shows that the regression model output fits the weld bead width with a good estimation. However, the lack of fit is insignificant. It should be recalled that in ANOVA if the regression model is significant whereas the lack-of-fit is insignificant simultaneously, the analysis is performed correctly.
Weld bead width (W) = −0.704 − 0.3008 × scanning speed (V) + 0.05992 × Current (C) (12) Figure 6 represents the main effect plots for the weld bead width response. The method used to distinguish the effect of parameters in the main effect plots is investigating the slope of the line. If the calculated line for a parameter is near horizontal; this parameter has no significant impact on response. However, a parameter with the highest inclined line will have the most significant effect. It is apparent from the plots that factor C as Current is the most influential variable. This obtained result was realized from the physical phenomenon of the TIG welding process. Figures 7 and 8 illustrate the contour and the response surface plots for weld bead width based on the input parameters. As can be seen in the contour plots, no interaction between input variables existed. There is a relation between scanning speed and current on the weld bead width, as shown in these figures. It is important to note that, the physical properties of the material extensively affect its behavior during the welding process. The 1.4418 stainless steel has low thermal conductivity among martensitic stainless steels. Low thermal conductivity causes lower heat dissipation and hence, reduces the cooling step after welding. 'In this metal, FZ experiences a higher local temperature than others, due to a low heat dissipation rate. For this reason, changes in input parameters influence weld bead geometries severely. Furthermore, from a metallurgical point of view, these physical properties result in coarser grain size in the fusion zone and low hardness in this area. Therefore, with increasing current, weld bead width grows. It can be seen that overheating of metal in the fusion zone occurred with increasing the welding current at a constant scanning speed. Extending in weld bead width was due to a physical phenomenon. This is widely referred to as the Marangoni effect. Differences in surface tension cause the Marangoni effect. Fluid flow for any fluid system is always from zones with higher surface tension towards lower. The surface tension depends on the temperature as the weld center has a higher temperature compared to edges, the surface tension is lower at the edges compared to the center. Thus, fluid in the molten state was drifted from the center to the edges. It should be noted that increasing the welding current causes to rise heat input in FZ and consequently, the temperature gradient, is increased. Subsequently, the FZ and HAZ regions shorten with the reduction of the welding current.
Scanning speed(V) R-Squared (Adj) = 88.92% Figure 6 represents the main effect plots for the weld bead width response. The method used to distinguish the effect of parameters in the main effect plots is investigating the slope of the line. If the calculated line for a parameter is near horizontal; this parameter has no significant impact on response. However, a parameter with the highest inclined line will have the most significant effect. It is apparent from the plots that factor C as Current is the most influential variable. This obtained result was realized from the physical phenomenon of the TIG welding process.  Figures 7 and 8 illustrate the contour and the response surface plots for weld bead width based on the input parameters. As can be seen in the contour plots, no interaction between input variables existed. There is a relation between scanning speed and current on the weld bead width, as shown in these figures. It is important to note that, the physical properties of the material extensively affect its behavior during the welding process. The 1.4418 stainless steel has low thermal conductivity among martensitic stainless steels. Low thermal conductivity causes lower heat dissipation and hence, reduces the cooling step after welding. `In this metal, FZ experiences a higher local temperature than others, due to a low heat dissipation rate. For this reason, changes in input parameters influence weld bead geometries severely. Furthermore, from a metallurgical point of view, these physical properties result in coarser grain size in the fusion zone and low hardness in this area. Therefore, with increasing current, weld bead width grows. It can be seen that overheating of metal in the fusion zone occurred with increasing the welding current at a constant scanning speed. Extending in weld bead width was due to a physical phenomenon. This is widely referred to as the Marangoni effect. Differences in surface tension cause the Marangoni effect. Fluid flow for any fluid system is always from zones with higher surface tension towards lower. The surface tension depends on the temperature as the weld center

1, 11, x FOR PEER REVIEW
has a higher temperature compared to edges, the surface tension is lower at the compared to the center. Thus, fluid in the molten state was drifted from the center edges. It should be noted that increasing the welding current causes to rise heat in FZ and consequently, the temperature gradient, is increased. Subsequently, the F HAZ regions shorten with the reduction of the welding current.
Furthermore, reducing the scanning speed has the same effect as increasing the ing current. In such cases, different parts in the welding joints experience the same on the microstructure. As described above, the heat input has a direct relation wi welding current but its relationship with scanning speed is in reverse. By reduci scanning speed, the residence time of the arc extends, and input power density -the energy transmitted to the welding zone in the form of heat within a unit of distan cends, in turn, it causes to increase the melted material volume. Thus, the weld be ometry, especially the weld bead width, is extended. Defects like undercutting oc the weld zone when the welding current is too high [23].  Furthermore, reducing the scanning speed has the same effect as increasing the welding current. In such cases, different parts in the welding joints experience the same effect on the microstructure. As described above, the heat input has a direct relation with the welding current but its relationship with scanning speed is in reverse. By reducing the scanning speed, the residence time of the arc extends, and input power density-the input energy transmitted to the welding zone in the form of heat within a unit of distanceascends, in turn, it causes to increase the melted material volume. Thus, the weld bead geometry, especially the weld bead width, is extended. Defects like undercutting occur in the weld zone when the welding current is too high [23].

Weld Depth of Penetration
According to the statistical report of an ANOVA performed on the weld bead width, all main input factors are significant, as shown in Table 8. Among second-order terms and interactions, none are effective. Equation (13) shows the regression model with coded values for the weld bead width considering significant parameters. The ANOVA table of weld depth of penetration indicates that the regression model output fits the weld bead

Weld Depth of Penetration
According to the statistical report of an ANOVA performed on the weld bead width, all main input factors are significant, as shown in Table 8. Among second-order terms and interactions, none are effective. Equation (13) shows the regression model with coded values for the weld bead width considering significant parameters. The ANOVA table of weld depth of penetration indicates that the regression model output fits the weld bead width with a good estimation. However, the lack of fit is insignificant. It should be recalled that in ANOVA if the regression model is significant but simultaneously the lack-of-fit is insignificant, the analysis is performed correctly. Weld depth of penetration (P) = −0.776 − 0.1017 × Scanning speed (V) + 0.02383 × Current (C) (13) Figure 9 represents the main effect plots for the weld depth of penetration response. The method used to distinguish the effect of parameters in the main effect plots is investigating the slope of the line. On the one hand, the calculated line for a parameter is near horizontal; but then this parameter has no significant impact on the response. On the other hand, a parameter with the highest inclined line will have the most significant effect. It is clear from the plots that factor C as Current is the most effective variable. This obtained result could be realized from the physical phenomenon of the TIG welding process. Moreover, it should be mentioned that current and scanning speed has direct and reverse impacts on the weld depth of penetration, respectively. tigating the slope of the line. On the one hand, the calculated line for a parameter is near horizontal; but then this parameter has no significant impact on the response. On the other hand, a parameter with the highest inclined line will have the most significant effect. It is clear from the plots that factor C as Current is the most effective variable. This obtained result could be realized from the physical phenomenon of the TIG welding process. Moreover, it should be mentioned that current and scanning speed has direct and reverse impacts on the weld depth of penetration, respectively.  Figures 10 and 11 illustrate the contour and the response surface plots for weld bead width based on the input parameters. As can be seen in contour plots, no interaction between input variables existed. The relation between scanning speed and current on the weld depth of penetration is shown in these figures. With increasing current, penetration grows. We conclude that Buoyancy force is negligible in this case. The density difference of the molten metal in the weld region caused Buoyancy force. With a rise in temperature, molten metal density decreases. Due to the Buoyancy effect fluid tends to flow from high  Figures 10 and 11 illustrate the contour and the response surface plots for weld bead width based on the input parameters. As can be seen in contour plots, no interaction between input variables existed. The relation between scanning speed and current on the weld depth of penetration is shown in these figures. With increasing current, penetration grows. We conclude that Buoyancy force is negligible in this case. The density difference of the molten metal in the weld region caused Buoyancy force. With a rise in temperature, molten metal density decreases. Due to the Buoyancy effect fluid tends to flow from high density to low density. In this case, it always tends to flow from the bottom to the surface of the weld. However, buoyancy forces are less potent as compared to other forces during TIG welding. As described above, overheating of metal in the fusion zone occurred with the increase of welding current at a constant scanning speed. It should be noted that increasing the welding current causes a rise in heat input in FZ and consequently, the temperature gradient decreases. Subsequently, the FZ and HAZ regions shorten with the reduction of the welding current. Hence, the welding current makes the FZ and HAZ regions deeper.
etals 2021, 11, x FOR PEER REVIEW density to low density. In this case, it always tends to flow from the bottom of the weld. However, buoyancy forces are less potent as compared to other TIG welding. As described above, overheating of metal in the fusion zone o the increase of welding current at a constant scanning speed. It should be n creasing the welding current causes a rise in heat input in FZ and consequen perature gradient decreases. Subsequently, the FZ and HAZ regions shorten duction of the welding current. Hence, the welding current makes the FZ gions deeper.
Furthermore, as described above, the heat input has a direct relation with current but its relationship with scanning speed is in reverse. By reducing speed, the residence time of the arc extends, and input power density -the transmitted to the welding zone in the form of heat within a unit of distance turn, it causes an increase of the melted material volume. Thus, the weld be especially the weld depth of penetration, is shortened. Therefore, in such c of the welded joint can remain unwelded, leading to the low tensile strength tion of the welded joints [24].

Optimization
After statistical analysis of data obtained from the experimental study, the regression equations are given by Equations (5)-(7), developed for TIG welding process modeling These equations explain logical relationships between input variables and responses. The Furthermore, as described above, the heat input has a direct relation with the welding current but its relationship with scanning speed is in reverse. By reducing the scanning speed, the residence time of the arc extends, and input power density-the input energy transmitted to the welding zone in the form of heat within a unit of distance-ascends, in turn, it causes an increase of the melted material volume. Thus, the weld bead geometry, especially the weld depth of penetration, is shortened. Therefore, in such cases, the root of the welded joint can remain unwelded, leading to the low tensile strength and elongation of the welded joints [24].

Optimization
After statistical analysis of data obtained from the experimental study, the regression equations are given by Equations (5)- (7), developed for TIG welding process modeling. These equations explain logical relationships between input variables and responses. The research carried out the optimization of the TIG welding process of 1.4418 stainless steel to access optimum welding settings, which minimize weld bead width and makes full penetration. Table 9 offers the criteria for the optimization of input factors. In Table 9, the importance values and composite desirability of all three responses are mentioned, and as one may see, weld penetration is more important than the other ones. The optimization was performed according to the criteria mentioned in Table 9. The optimum setting was achieved by using Minitab software, which is presented in Table 10. In addition to maximizing depth penetration and minimizing weld bead width, it should be noted that the low productivity of TIG welding should be covered and optimized [25].  In criteria 1 of the optimum conditions, partial penetration occurs. Criteria 3 is chosen as the best condition of this process because of the better quality and the higher scanning speed compared to criteria 2. Optimal settings of experimental tests carried out at the optimized predicted results due to verifying statistical analysis. Table 10 depicts the predicted results of criteria 2. Figure 12 describes Derringer's desirability with graphs.   Figure 13 illustrates an overlay plot composed of the contour plots from each response stacked on top of each other to create the combined plot. The developed white area reveals the final optimal factor settings of experimental tests. On the other hand, the undesirable area grayed out.
The over plot suggests an accepted area to achieve the desired condition for responses [26].

Conclusions
In this investigation, the simulation model is primarily developed for the TIG welding process on 1.4418 stainless steel sheets followed by a mathematical model with Abaqus and validated with experimental approaches. In the following, the effect of the scanning speed and the current are investigated in this process by using Minitab software. According to numerical, experimental, and statistical works, the following results concluded: Figure 13. Overlay contour plot.

Conclusions
In this investigation, the simulation model is primarily developed for the TIG welding process on 1.4418 stainless steel sheets followed by a mathematical model with Abaqus and validated with experimental approaches. In the following, the effect of the scanning speed and the current are investigated in this process by using Minitab software. According to numerical, experimental, and statistical works, the following results concluded:

1.
Among the input variables of 1.4418 stainless steel TIG welding, scanning speed is the most significant factor with a reverse effect on all weld bead dimensions.

2.
Results revealed that increasing the current and decreasing scanning speed, causes the weld bead width to widen and the depth of penetration extended, due to increased input energy transmitted to the material and it conducts to an increase in the melted volume.

3.
By applying a multi-objective optimization, through the desirability approach, the optimum settings of the TIG welding process are scanning speed (V) = 2.69 mm/s, Current (C) = 120 A. 4.
In the highest desired condition, the weld bead width of the welded part diminishes to 5.7 mm, and the depth of penetration increases to 1.8 mm.

5.
By conducting the thermal FEM analysis tests for 10 different cases and comparing them with experimental approaches, the predicted model has an average 9.7% relative error. As a result, this model shows the robustness of the present FEM model utilized for TIG Arc welding of 1.4418 stainless steel. 6.
With obtained results in this research, we can choose optimum input factors to reach a convenient weld seam with full penetration weld and wide weld bead. Moreover, with these input factors, optimum samples with high mechanical strength can be reached.
Although the results of this study are promising, the research can be extended to analyze the presented method on a wide range of materials. Future studies can also include mechanical analyses such as residual stress and surface deformation of the weld pool. Finally, microstructural analyses of the welding area can be considered to describe the heat-affected zone (HAZ).