Mechanical Properties and Microstructure of TIG and ATIG Welded 316L Austenitic Stainless Steel with Multi-Components Flux Optimization Using Mixing Design Method and Particle Swarm Optimization (PSO)

In this study, the effects of pseudo-ternary oxides on mechanical properties and microstructure of 316L stainless steel tungsten inert gas (TIG) and activating tungsten inert gas (ATIG) welded joints were investigated. The novelty in this work is introducing a metaheuristic technique called the particle swarm optimization (PSO) method to develop a mathematical model of the ultimate tensile strength (UTS) in terms of proportions of oxides flux. A constrained optimization algorithm available in Matlab 2020 optimization toolbox is used to find the optimal percentages of the selected powders that provide the maximum UTS. The study indicates that the optimal composition of flux was: 32% Cr2O3, 43% ZrO2, 8% Si2O, and 17% CaF2. The UTS was 571 MPa for conventional TIG weld and rose to 600 MPa for the optimal ATIG flux. The obtained result of hardness for the optimal ATIG was 176 HV against 175 HV for conventional TIG weld. The energy absorbed in the weld zone during the impact test was 267 J/cm2 for the optimal ATIG weld and slightly higher than that of conventional TIG weld 256 J/cm2. Fracture surface examined by scanning electron microscope (SEM) shows ductile fracture for ATIG weld with small and multiple dimples in comparison for TIG weld. Moreover, the depth of optimized flux is greater than that of TIG weld by two times. The ratio D/W was improved by 3.13 times. Energy dispersive spectroscopy (EDS) analysis shows traces of the sulfur element in the TIG weld zone.


Introduction
Austenitic stainless steels are the most common stainless steel. They are used in many fields as oil, shipbuilding, machinery, and marine applications. They are characterized by good strength, high toughness, and excellent corrosion resistance. Fusion welding is the most popular method to join workpieces in many industrial applications.
Tungsten Inert Gas (TIG) welding is a widespread process in industries, but a limited thickness can be joined in a single pass with this process. TIG welding is also very sensitive to the chemical composition of the base metal. Moreover, materials thickness greater than 3 mm requires multiple passes to achieve full penetration weld; therefore, the productivity of the process is reduced. The activated tungsten inert gas (ATIG) technic is a variant of conventional TIG welding. In ATIG welding, a thin layer of flux is deposited on the workpiece before the welding operation. ATIG offers the possibility of increasing the penetration depth using the same equipment and parameters as conventional TIG. Comparatively to TIG, ATIG has many advantages. It eliminates the need for edge preparation, increases the penetration depth, and reduces the number of weld pass [1,2]. Three mechanisms have been proposed; the increase in ATIG weld penetration can be attributed to the reversal of Marangoni convection [3,4]. The second mechanism is the arc constriction proposed by D. S. Howse et al. [5,6]. The third mechanism, proposed by Sire, consists of a constriction of arc by using a flux characterized by high melting point and high electrical resistivity [7,8]. In recent years, considerable works were carried out using Taguchi method design to optimize TIG welding process parameters to improve both depth and weld aspects [9,10]. Some studies were oriented to study the influence of welding parameters on enhancing mechanical properties [11,12]. Other works were dedicated to optimizing the composition of paste to improve the depth and weld aspects [13,14] using the mixing method. Recently, other methods like genetic algorithms, simulated annealing, particle swarm optimization [15], etc., have been used to optimize solution in many industries [16,17]. In this work, the tensile strength of weld joints is optimized as an output parameter by varying the input parameters that are the combinations of oxides powder without varying the welding parameters.
The objective of this work is to elaborate an appropriate composition of flux to improve mechanical properties of the full penetrated 316L stainless steel ATIG weld. Mixing method simplex design degree four was combined with particle swarm optimization (PSO) method to optimize combination powders. The developed mathematical model and optimization methodology can help researchers and engineers working in developing fluxes for weld improvements. This study is a contribution to expanding research dedicated to ATIG welding of 316L.

Material
The material used in this study is the austenitic stainless steel grade 316L. The chemical composition is shown in Table 1. Rectangular pieces of 6 mm thickness were cut from the received plate to perform welding of 20 cm line.

TIG-ATIG Welding Platform
First of all, the plates were cleaned with acetone. The oxide powders were dried in the furnace for 1 h at 180 • C to eliminate the humidity. Then, a thin layer of a mixed powder with methanol was applied using a brush to the surface subject to the welding as shown in Figure 1a. The mean coating density of flux was about 4-5 mg/cm 2 .
The tungsten inert gas welding machine was used. The electrode used has a diameter of 3.2 mm and the torch was mounted on a motorized carriage as shown in Figure 1b. The experimental parameters selected for welding are presented in Table 2.   Figure 2a shows a good external aspect for TIG weld line, also surface ATIG weld line in Figure 2b is free of slugs, which indicates that the optimal flux is well consumed with small traces of residue. Pieces were cut from the welded plates for mechanical testing, weld morphology, and microstructure study according to the schematic drawing shown in Figure 3. The cutting machine used is Hydro-Jet Eco 0515 SL-Micro waterjet Cutting (KNUTH Germany).   Figure 2a shows a good external aspect for TIG weld line, also surface ATIG weld line in Figure 2b is free of slugs, which indicates that the optimal flux is well consumed with small traces of residue.   Figure 2a shows a good external aspect for TIG weld line, also surface ATIG weld line in Figure 2b is free of slugs, which indicates that the optimal flux is well consumed with small traces of residue. Pieces were cut from the welded plates for mechanical testing, weld morphology, and microstructure study according to the schematic drawing shown in Figure 3. The cutting machine used is Hydro-Jet Eco 0515 SL-Micro waterjet Cutting (KNUTH Germany). Pieces were cut from the welded plates for mechanical testing, weld morphology, and microstructure study according to the schematic drawing shown in Figure 3. The cutting machine used is Hydro-Jet Eco 0515 SL-Micro waterjet Cutting (KNUTH Germany).

Tensile Test
Design of experiment and mathematical modelling will be applied to the tensile strength property. Optimal flux, which will give the maximum UTS resulting from the tensile test, is used for the rest of the study.
The tensile tests were performed with a computer control electrohydraulic servo universal testing machine model WAW-300E (Jinan testing equipment IE, Jinan, China), at a test rate of 0.5 mm/min, at 0.5 kN/s load rate, and at low strain rate of 1.6 × 10 −4 s −1 .

Weld Bead Aspect
The cross-sections of the weld beads for ATIG with optimal flux and conventional TIG were photographed using an optical microscope CAROLINA (CAROLINA, Burlington, NJ, USA). The morphology of the welds were checked using Motic Images plus version 2.0 software integrated with an optical microscope.

Microstructure Assessment
The microstructural characterization of the fusion zone of both TIG and ATIG welding has been analyzed. Micrographs were taken on JEOL JSM-7600F scanning electronic microscope, (SEM). The areas image processing software from Microvision Instruments (Microvision Instruments, Imager M2.m, Paris, France) was used to measure the ferrite volume proportions.

Hardness Test
Vickers hardness tests were performed by a digital hardness tester model HVS-50 (SCTMC, Shanghai, China) with a standard load of 98 N. Eight indentations in the weld bead on fusion zone (FZ) and at the heat-affected zone (HAZ) were performed on each sample, with about 0.5 mm between two indentations, as shown in Figure 4.

Tensile Test
Design of experiment and mathematical modelling will be applied to the tensile strength property. Optimal flux, which will give the maximum UTS resulting from the tensile test, is used for the rest of the study.
The tensile tests were performed with a computer control electrohydraulic servo universal testing machine model WAW-300E (Jinan testing equipment IE, Jinan, China), at a test rate of 0.5 mm/min, at 0.5 kN/s load rate, and at low strain rate of 1.6 × 10 −4 s −1 .

Weld Bead Aspect
The cross-sections of the weld beads for ATIG with optimal flux and conventional TIG were photographed using an optical microscope CAROLINA (CAROLINA, Burlington, NJ, USA). The morphology of the welds were checked using Motic Images plus version 2.0 software integrated with an optical microscope.

Microstructure Assessment
The microstructural characterization of the fusion zone of both TIG and ATIG welding has been analyzed. Micrographs were taken on JEOL JSM-7600F scanning electronic microscope, (SEM). The areas image processing software from Microvision Instruments (Microvision Instruments, Imager M2.m, Paris, France) was used to measure the ferrite volume proportions.

Hardness Test
Vickers hardness tests were performed by a digital hardness tester model HVS-50 (SCTMC, Shanghai, China) with a standard load of 98 N. Eight indentations in the weld bead on fusion zone (FZ) and at the heat-affected zone (HAZ) were performed on each sample, with about 0.5 mm between two indentations, as shown in Figure 4.

Tensile Test
Design of experiment and mathematical modelling will be applied to the tensile strength property. Optimal flux, which will give the maximum UTS resulting from the tensile test, is used for the rest of the study.
The tensile tests were performed with a computer control electrohydraulic servo universal testing machine model WAW-300E (Jinan testing equipment IE, Jinan, China), at a test rate of 0.5 mm/min, at 0.5 kN/s load rate, and at low strain rate of 1.6 × 10 −4 s −1 .

Weld Bead Aspect
The cross-sections of the weld beads for ATIG with optimal flux and conventional TIG were photographed using an optical microscope CAROLINA (CAROLINA, Burlington, NJ, USA). The morphology of the welds were checked using Motic Images plus version 2.0 software integrated with an optical microscope.

Microstructure Assessment
The microstructural characterization of the fusion zone of both TIG and ATIG welding has been analyzed. Micrographs were taken on JEOL JSM-7600F scanning electronic microscope, (SEM). The areas image processing software from Microvision Instruments (Microvision Instruments, Imager M2.m, Paris, France) was used to measure the ferrite volume proportions.

Hardness Test
Vickers hardness tests were performed by a digital hardness tester model HVS-50 (SCTMC, Shanghai, China) with a standard load of 98 N. Eight indentations in the weld bead on fusion zone (FZ) and at the heat-affected zone (HAZ) were performed on each sample, with about 0.5 mm between two indentations, as shown in Figure 4.

Impact Test
Impact testing was performed on three samples for TIG and three samples for ATIG weld only in the fusion zone (FZ) with the Charpy "V" notch impact testing machine model JBS-500 (Jinan testing equipment IE, Jinan, China), specimens are shown in Figure 5.

Impact Test
Impact testing was performed on three samples for TIG and three samples for ATIG weld only in the fusion zone (FZ) with the Charpy "V" notch impact testing machine model JBS-500 (Jinan testing equipment IE, Jinan, China), specimens are shown in Figure  5.

Design of Experiment Methodology
To get the maximum amount of useful information with the least amount of experimentation, the design of experiments (DOE) has been used; mixing method in Minitab 17 software was applied. First, to compare the effect of the mono-flux oxides on UTS of welds, thirteen oxides (SiO2, TiO2, Fe2O3, MnO2, Cr2O3, ZrO2, CaO, Mn2O3, V2O5, MoO3, SrO, Co2O3, and MgO) were tested. Among these thirteen oxides, three oxides Cr2O3, TiO2, and ZrO2 that gave the best UTS were selected to be used in the mixing design method. Based on the simplex lattice degree four designs, nineteen combinations from the selected oxides have been prepared. For each combination, the three selected oxides Cr2O3, TiO2, and ZrO2 vary and 25% of (8% SiO2 + 17% CaF2) was added and kept fixed to get a pseudo ternary combination. Particle swarm optimization (PSO) method is used to establish an equation relating UTS to the proportions of the selected oxides. Finally, Matlab R2020 is carried out to obtain the optimal combination of oxides, hence, to maximize the UTS of weld through a constrained optimization algorithm.
Silica (SiO2) has been added because it increases the current carrying capacity. Consequently, silicate (SiO2) increases arc voltage and depth penetration of the weld bead [18]. On the other hand, the addition of calcium fluoride (CaF2) had several advantages in that it reduces the dissolved silicon content of weld metal, prevents the deleterious effect of silicon on hot cracking, and lowers the melting range of the flux. Moreover, fluorine gases escaped from the weld pool interact with outer arc electrons leading to constrict arc [19]. The presence of fluorine in arc welding reduces the anode spot and tends to increase the energy density of the heat source and electromagnetic force in the weld pool. As a result, relatively narrow and deep weld morphology is formed [20,21]. Nineteen weld lines were executed and three samples were cut for the tensile test from each of the nineteen combinations.

Mathematical Modelling
For the first step of work, we developed a mathematical model where the UTS is written in terms of these selected oxides' percentages. The technic used to get a mathematical model will be presented later in this paper. In the second step, the optimal combination that maximizes UTS is determined.

Design of Experiment Methodology
To get the maximum amount of useful information with the least amount of experimentation, the design of experiments (DOE) has been used; mixing method in Minitab 17 software was applied. First, to compare the effect of the mono-flux oxides on UTS of welds, thirteen oxides (SiO 2 , TiO 2 , Fe 2 O 3 , MnO 2 , Cr 2 O 3, ZrO 2 , CaO, Mn 2 O 3 , V 2 O 5 , MoO 3 , SrO, Co 2 O 3, and MgO) were tested. Among these thirteen oxides, three oxides Cr 2 O 3 , TiO 2 , and ZrO 2 that gave the best UTS were selected to be used in the mixing design method. Based on the simplex lattice degree four designs, nineteen combinations from the selected oxides have been prepared. For each combination, the three selected oxides Cr 2 O 3 , TiO 2 , and ZrO 2 vary and 25% of (8% SiO 2 + 17% CaF 2 ) was added and kept fixed to get a pseudo ternary combination. Particle swarm optimization (PSO) method is used to establish an equation relating UTS to the proportions of the selected oxides. Finally, Matlab R2020 is carried out to obtain the optimal combination of oxides, hence, to maximize the UTS of weld through a constrained optimization algorithm.
Silica (SiO 2 ) has been added because it increases the current carrying capacity. Consequently, silicate (SiO 2 ) increases arc voltage and depth penetration of the weld bead [18]. On the other hand, the addition of calcium fluoride (CaF 2 ) had several advantages in that it reduces the dissolved silicon content of weld metal, prevents the deleterious effect of silicon on hot cracking, and lowers the melting range of the flux. Moreover, fluorine gases escaped from the weld pool interact with outer arc electrons leading to constrict arc [19]. The presence of fluorine in arc welding reduces the anode spot and tends to increase the energy density of the heat source and electromagnetic force in the weld pool. As a result, relatively narrow and deep weld morphology is formed [20,21]. Nineteen weld lines were executed and three samples were cut for the tensile test from each of the nineteen combinations.

Mathematical Modelling
For the first step of work, we developed a mathematical model where the UTS is written in terms of these selected oxides' percentages. The technic used to get a mathematical model will be presented later in this paper. In the second step, the optimal combination that maximizes UTS is determined.
Finally, based on the optimal formulation, ATIG and TIG weld lines have been carried out on a single plate. Thirteen oxides were tested. The weld line was executed on the plain plate using welding parameters cited in Table 2. It can be seen in Table 3 that the highest value of UTS is 565 MPa, which was obtained for the sample welded with ZrO 2 flux, followed by the sample welded with Cr 2 O 3 flux with 559 MPa and, then by the sample welded with TiO 2 flux with 542 MPa. Therefore, the selected oxides were Cr 2 O 3 , TiO 2 , and ZrO 2 . Based on the simplex lattice degree four design, nineteen combinations have been prepared. As shown in Table 4, for each combination, the three selected oxides Cr 2 O 3 , TiO 2 , and ZrO 2 vary and 25% of (8% SiO 2 + 17% CaF 2 ) is kept fixed.  Table 5 shows the Tensile test (UTS) results as well as the standards deviation (σ). The value of σ is less than 21 MPa for UTS results, which indicates that the disparities are acceptable.

Mathematical Model
In this study, the three selected oxides (Cr 2 O 3 , TiO 2 , ZrO 2 ) percentages vary for the combined oxide with a fixed 25% of (SiO 2 + CaF 2 ) oxide to form the pseudo-ternary ATIG are used. The UTS is then expressed as a function of the percentages of the three oxides as follows: UTS = f (X1, X2, X3, X4). X1, X2, X3 denote, respectively, the proportions (in terms of percentages) of the three oxides, whereas the proportion of the fourth oxide is set to be X4 = 25%.
The chosen mathematical model that describes the effect of these percentages on the UTS is the coupling of second-order model as proposed in [22] with conventional linear regression as in [23]. This model includes three components: (i) the linear effect of the proportions, (ii) their quadratic effects, and (iii) the interactions between those proportions. The general form of the mathematical model is given below in expression 1: With a fixed X(4) = 25%.

Optimization Process Details
The first step objective was to find the optimal parameters that minimize the quadratic error between the measured UTS and the UTS provided by the model. Thus, the modeling problem is converted into an optimization problem having the model coefficients as decision variables and the quadratic error as the objective function to be minimized. Since this problem may present many irregularities such as non-convexity of the criterion and many feasibility constraints, we used a metaheuristic technique called particle swarm optimization (PSO) to solve this problem. In PSO, a set of candidate solutions (coefficients of the model) are initialized randomly within the search-space limits and then "flown" progressively toward a sub-optimal solution while combining three components: (i) follow their current velocities, (ii) go back to their best positions visited so far, and (iii) go to the position of the best neighbor. Each particle in the group (called swarm) is assigned four vectors: the position including the model parameters, a velocity, the best personal position, and the global best position. The move equations of each particle at the kth iteration of the optimization process are provided below [24,25]: As provided in many papers such as [24,25], the inertia weight decreases linearly from 0.9 to 0.4. c 1 = c 2 = 0.75 are the cognitive and social factors. r 1 and r 2 are two random numbers generated randomly between 0 and 1. The model parameters' search limits are set, respectively, to −1 and +1 for the overall model stability [25]. The optimization process is stopped after a certain number of iterations (k max = 5000). After this global search step, a local search is conducted to find better solutions around the global sub-optimal solution yielded by the PSO algorithm. After running the optimization process many times, the selected UTS model as expressed by X1, X2, and X3 is given below: Three performance metrics are used to evaluate the effectiveness of the developed models as follows [25]: Mean Absolute Percentage Error (MAPE (%)): • Coefficient of determination (R 2 (%)) • Root Mean Square Error (RMSE (MPa)) With N 2 = 19 (number of measurements) The accuracy of this model as measured by the above cited measurement indicators are found to be as follows: MAPE = 0.7952%, R 2 = 72.04%, and RMSE = 5.5955 MPa.
The second step of the optimization process consists of finding the optimal percentages of the three oxides that provide the maximum UTS. To achieve this goal, a constrained optimization algorithm is used. This algorithm has provided the following optimal combination shown in Table 6. The predicted value is UTS = 588 MPa as shown in Table 7.

Experimental Validation
The validation test was the last step in the experimental process. A confirmation test was performed according to the optimum flux composition. Table 8 shows that UTS for ATIG weld reached 600 MPa, which is higher than that of conventional TIG welding (571 MPa). Moreover, the UTS of ATIG weld is greater than the expected value calculated by the mathematical model (588 MPa). The UTS of the optimal flux (ATIG) weld (600 MPa) is close to that of base metal (624 MPa) shown in Table 9. The ATIG welding can reduce the heat input per unit length in welds and the residual stress of the weld can be reduced [26,27]. On the other hand, the increase of UTS about 29 MPA in favor of ATIG comparatively to conventional TIG weld can be attributed to the increase of the retained ferrite volume proportions ATIG weld. Moreover, high heat input in TIG weld exhibits a coarse ferrite distribution and lower strength than the lower heat input in ATIG weld, which has a relatively fine ferrite distribution [27]. The value of standard deviation (σ) is less than 5 MPa.  Table 9. UTS comparison between mathematical model, ATIG, TIG, and base metal.

Tensile Break Zone Investigation
Micrographs of the break zone in the tensile test were conducted with SEM. The images show the same profile with the formation of multiple dimples, which demonstrates that the fracture is in the ductile mode for both cases of TIG and ATIG weld as shown in Figure 6a We notice that the fracture location for both welds occurs in base metal as shown in Figure 7. The results show the elongation percentages very close to 22.63% for ATIG specimen against 20.52% for TIG weld.

Weld Bead Morphology
The values obtained for depth D and D/W ratio of the optimal flux ATIG weld are listed in Table 10. The results show that the depth of optimized flux is greater than TIG weld by two times. The ratio D/W was improved by 3.13 times. The weld of optimal flux is a full-penetrated weld (6.8 mm). The D/W value obtained for ATIG weld is 0.72. We notice that the fracture location for both welds occurs in base metal as shown in Figure 7. The results show the elongation percentages very close to 22.63% for ATIG specimen against 20.52% for TIG weld. We notice that the fracture location for both welds occurs in base metal as shown in Figure 7. The results show the elongation percentages very close to 22.63% for ATIG specimen against 20.52% for TIG weld.

Weld Bead Morphology
The values obtained for depth D and D/W ratio of the optimal flux ATIG weld are listed in Table 10. The results show that the depth of optimized flux is greater than TIG weld by two times. The ratio D/W was improved by 3.13 times. The weld of optimal flux is a full-penetrated weld (6.8 mm). The D/W value obtained for ATIG weld is 0.72.

Weld Bead Morphology
The values obtained for depth D and D/W ratio of the optimal flux ATIG weld are listed in Table 10. The results show that the depth of optimized flux is greater than TIG weld by two times. The ratio D/W was improved by 3.13 times. The weld of optimal flux is a full-penetrated weld (6.8 mm). The D/W value obtained for ATIG weld is 0.72. This can be explained by the fact that Arc weld constriction and reversed Marangoni convection in ATIG occur when: -Fluorine from flux migrates to the arc weld and contributes to enhancing the energy density according to the constriction arc mechanism explained before, as mentioned in several works [19,28]. -Oxygen liberated from oxides as surfactant element affects the surface tension of the molten metal resulting in a centripetal movement, the metal moves from the edges to the center as cited in related works [29,30]. A full penetration weld is performed in a single pass without edge preparation or the use of filler metal, which meets the needs of industries as shown in Figure 8b.
However, in TIG, Marangoni convection occurs. A molten metal moves from the center of the weld pool to the edges as pure metal leading to a wide and shallow weld bead as shown in Figure 8a. This can be explained by the fact that Arc weld constriction and reversed Marangoni convection in ATIG occur when: -Fluorine from flux migrates to the arc weld and contributes to enhancing the energy density according to the constriction arc mechanism explained before, as mentioned in several works [19,28]. -Oxygen liberated from oxides as surfactant element affects the surface tension of the molten metal resulting in a centripetal movement, the metal moves from the edges to the center as cited in related works [29,30]. A full penetration weld is performed in a single pass without edge preparation or the use of filler metal, which meets the needs of industries as shown in Figure 8b.
However, in TIG, Marangoni convection occurs. A molten metal moves from the center of the weld pool to the edges as pure metal leading to a wide and shallow weld bead as shown in Figure 8a.

Microstructural Assessment
The differences in microstructure between TIG weld bead and ATIG weld are shown in Figures 9a,b and 10a,b respectively. In both cases, the ferrite morphology has a discontinuous skeletal network of delta ferrite (δ) structures in a predominant austenite matrix. Alloy 316L with Creq/Nieq that is 1.67 solidifies in ferritic austenitic mode [31]. The delta ferrite is located mainly at the dendrite axes. The residual primary ferrite results from incomplete δ → γ transformation during solidification.

Microstructural Assessment
The differences in microstructure between TIG weld bead and ATIG weld are shown in Figure 9a,b and Figure 10a,b respectively. In both cases, the ferrite morphology has a discontinuous skeletal network of delta ferrite (δ) structures in a predominant austenite matrix. Alloy 316L with Creq/Nieq that is 1.67 solidifies in ferritic austenitic mode [31]. The delta ferrite is located mainly at the dendrite axes. The residual primary ferrite results from incomplete δ → γ transformation during solidification. The delta ferrite (δ) in the TIG weld zone is coarser than that of the ATIG weld as shown in Figures 9b and 10b, respectively. This difference in delta ferrite (δ) size is related to the heat provided by the weld bead, which is lesser in the case of ATIG weld.  The delta ferrite (δ) in the TIG weld zone is coarser than that of the ATIG weld as shown in Figures 9b and 10b, respectively. This difference in delta ferrite (δ) size is related to the heat provided by the weld bead, which is lesser in the case of ATIG weld. In ATIG and TIG weld zone, the ferrite is the result of the incomplete primary delta ferrite-austenite transformation. The proportions of delta ferrite diminish from weld zone to base metal is as shown in Figure 11a,b.
In the ATIG heat-affected zone, austenitic grains crossed by parallel elongated inclined stringers of delta ferrite as shown in Figure 11a. In TIG heat-affected zone, the globular delta ferrite is randomly distributed and in the form of small slats of reduced size in the same rolling direction in a matrix of austenite as shown in Figure 11b. The delta ferrite (δ) in the TIG weld zone is coarser than that of the ATIG weld as shown in Figure 9b and Figure 10b, respectively. This difference in delta ferrite (δ) size is related to the heat provided by the weld bead, which is lesser in the case of ATIG weld.
In ATIG and TIG weld zone, the ferrite is the result of the incomplete primary delta ferrite-austenite transformation. The proportions of delta ferrite diminish from weld zone to base metal is as shown in Figure 11a Figure 12a,b represent the measurements of ferrite volume proportions in the austenite matrix. The advisable proper amount of δ-ferrite in austenitic stainless steel welds is less than 10% volume to ensure a better ductility, toughness, and corrosion resistance, and no less than 5% to avoid solidification cracking [32]. In the ATIG heat-affected zone, austenitic grains crossed by parallel elongated inclined stringers of delta ferrite as shown in Figure 11a. In TIG heat-affected zone, the globular delta ferrite is randomly distributed and in the form of small slats of reduced size in the same rolling direction in a matrix of austenite as shown in Figure 11b. Figure 12a,b represent the measurements of ferrite volume proportions in the austenite matrix. The advisable proper amount of δ-ferrite in austenitic stainless steel welds is less than 10% volume to ensure a better ductility, toughness, and corrosion resistance, and no less than 5% to avoid solidification cracking [32]. input resulted in a slow cooling rate and a further transformation from ferrite phase to the austenite phase. Consequently, a ferrite volume proportion is reduced (5.56%). The ATIG welding was associated with a low heat input resulting in higher ferrite content (8.63%), which is in good agreement with several works [35,36]. On the other hand, the high heat input leads to coarse substructure during solidification in the case of TIG weld as shown in Figure 9a,b, which results in a more widely spaced ferrite network. But in the case of ATIG weld, heat provided is low; a fast cooling rate occurs and consequently a finer skeletal ferrite is formed as shown in Figure 10a,b [37].

Hardness Test
The hardness values are shown in Table 12. It is visible that the hardness of ATIG weld and TIG weld in both FZ and HAZ are very close. The delta-ferrite volume proportions in the weld metals are in both TIG and ATIG welds increased, and have a beneficial effect in increasing the hardness of as received 316L stainless steel welds (160 HV).
In ATIG weld, the property of hardness is not affected by the optimal flux used. On the other hand, the standard deviation is less than 9 HV, which attests to the small disparities in the obtained hardness values of the maximum and minimum. This result indicates good hardness homogeneities in the joints.

Impact Test
The impact tests were carried out on the fusion zone in ATIG and TIG welds. The experimental values obtained for the impact tests are shown in Table 13. The energy absorbed in the fusion zone in the case of ATIG weld (267 J/cm 2 ) is slightly higher than that of TIG weld (256 J/cm 2 ) by 11 J/cm 2 . The standard deviation is less than 15 J/cm 2 . The results reported in Table 11 show the measurements taken in seven different locations in the weld zone and the mean value of ferrite volume. The proportions are up to 5.56% in the conventional TIG process and 8.63% for ATIG weld. The higher retained ferrite may be attributed to the rate of cooling of the welds. Elements from optimal flux can constrict arc weld. The constriction of the arc will increase the temperature at the anode due to the increase in current density and arc voltage as reported in many studies [33,34].
The ATIG welding with optimal flux increased the energy density of the heat source that leads to low heat input. However, the arc heat of TIG welding without flux has a lower energy density; therefore, high heat input is provided to the workpiece. High heat input resulted in a slow cooling rate and a further transformation from ferrite phase to the austenite phase. Consequently, a ferrite volume proportion is reduced (5.56%). The ATIG welding was associated with a low heat input resulting in higher ferrite content (8.63%), which is in good agreement with several works [35,36]. On the other hand, the high heat input leads to coarse substructure during solidification in the case of TIG weld as shown in Figure 9a,b, which results in a more widely spaced ferrite network. But in the case of ATIG weld, heat provided is low; a fast cooling rate occurs and consequently a finer skeletal ferrite is formed as shown in Figure 10a,b [37].

Hardness Test
The hardness values are shown in Table 12. It is visible that the hardness of ATIG weld and TIG weld in both FZ and HAZ are very close. The delta-ferrite volume proportions in the weld metals are in both TIG and ATIG welds increased, and have a beneficial effect in increasing the hardness of as received 316L stainless steel welds (160 HV). In ATIG weld, the property of hardness is not affected by the optimal flux used. On the other hand, the standard deviation is less than 9 HV, which attests to the small disparities in the obtained hardness values of the maximum and minimum. This result indicates good hardness homogeneities in the joints.

Impact Test
The impact tests were carried out on the fusion zone in ATIG and TIG welds. The experimental values obtained for the impact tests are shown in Table 13. The energy absorbed in the fusion zone in the case of ATIG weld (267 J/cm 2 ) is slightly higher than that of TIG weld (256 J/cm 2 ) by 11 J/cm 2 . The standard deviation is less than 15 J/cm 2 .  Figure 12 represents the fractographs of the impact Charpy "V" notch test. The images show the formation of multiple dimples, which demonstrate that the fracture is a ductile mode in both cases of TIG and ATIG as shown in Figure 13a,b. However, in ATIG weld the dimples are finer with the presence of voids, which can explain the slightly high value of absorbed energy comparatively to TIG weld. The ductile fracture mode leads to good resistance to sudden impact loads. Multiple dimples attest for high impact energy withstand as reported by several authors [38,39].   Figure 12 represents the fractographs of the impact Charpy "V" notch test. The images show the formation of multiple dimples, which demonstrate that the fracture is a ductile mode in both cases of TIG and ATIG as shown in Figure 13a,b. However, in ATIG weld the dimples are finer with the presence of voids, which can explain the slightly high value of absorbed energy comparatively to TIG weld. The ductile fracture mode leads to good resistance to sudden impact loads. Multiple dimples attest for high impact energy withstand as reported by several authors [38,39]. The results of EDS/SEM analysis in Figure 14a and Table 14 are shown in the case of ATIG weld, with the same level of silicon as in base metal. There are no trace of unwanted elements in the case of ATIG welds. These results can explain the good resistance to impact test. The results of EDS/SEM analysis in Figure 14a and Table 14 are shown in the case of ATIG weld, with the same level of silicon as in base metal. There are no trace of unwanted elements in the case of ATIG welds. These results can explain the good resistance to impact test.
(a) (b) Figure 13. Fractograph of Austenitic stainless steel 316L impact Charpy "V" notch for (a) TIG Welded Zone (500×) and (b) ATIG Welded Zone (500×) The results of EDS/SEM analysis in Figure 14a and Table 14 are shown in the case of ATIG weld, with the same level of silicon as in base metal. There are no trace of unwanted elements in the case of ATIG welds. These results can explain the good resistance to impact test.    The results of EDS/SEM analysis in Figure 14b and in Table 14 show the presence of oxygen in TIG weld, which is probably ascribed to insufficient protection of the weld pool that affects mechanical properties, particularly the toughness. The presence of sulfur leads to the formation of low-melting eutectics with iron, chromium, and nickel, which can alter the mechanical properties of TIG weld beads [40]. The obtained results with EDS/SEM related to TIG weld may explain the decrease in resistance to impact test comparatively to the ATIG weld.

Conclusions
In the present work, ATIG weld has been investigated and compared to conventional TIG weld. The starting point of this work is elaborating the optimal flux to maximize the UTS ATIG weld. The mixing design of the experiment combined to the particle swarm optimization (PSO) method is used to minimize the number of trials, which reduces the cost of materials as well as time. Based on the obtained results in this investigation, the following conclusions can be drawn: -Mixing design of the experiment combined to the particle swarm optimization (PSO) method is among the novelties of this work. Optimal flux was composed by 32% Cr 2 O 3 , 43% ZrO 2 , 8% SiO 2 , and 17% CaF 2 . The optimal flux raised from this method improves the mechanical properties in comparison to conventional TIG weld bead. Transversal tensile testing of the produced ATIG weld has a UTS value (600 MPa) close to parent metal (624 MPa). On the other hand, the UTS of TIG weld fell to 571 MPa. ATIG welding is more resistant than conventional TIG welding to sudden impact loads. The ATIG weld hardness results were close to those of conventional TIG weld. - The ATIG depth weld bead reached 6.80 mm and the weld aspect ratio increased 3.13. The ATIG depth weld was increased by two times in comparison to the conventional TIG weld bead. The fully penetrated ATIG weld is ascribed to two mechanisms cited earlier. The reversal Marangoni mechanism owing to oxygen liberated from flux and the constriction of ATIG arc weld related to the migration of fluorine from the flux to the arc. - The microstructure of welds in both welds is composed of matrix austenite interspersed by skeleton δ-ferrite. The ferrite volume proportions in ATIG weld around 8.63% and decreased to 5.56% in conventional TIG welding. The δ-ferrite in TIG weld is coarser comparatively to that of ATIG δ-ferrite.