Numerical Simulations Based on a Meshfree Method for Nickel-Steel Welded Joint Manufactured by Micro-Jet Cooling

The article presents a numerical–experimental approach to the weldability and mechanical resistance of the joint of Alloy 59 (2.4605, nickel-chromium-molybdenum) and S355J2W (1.8965) structural steel manufactured by the MIG process with the use of micro-jet cooling. This research was considered because the standard MIG process does not guarantee the procurement of a mixed hard-rusting structural steel superalloy weld of a repeatable and acceptable quality. Welds made through the classic MIG process express cracks that result from their unfavorable metallographic microstructure, while the joint supported by micro-jet cooling does not reflect any cracks and has a high strength with good flexibility. This was achieved by the application of helium for cooling. The joining technology was also considered in the numerical stage, represented by calculations in situ. For this purpose, the fundamental solution method (FSM) for the simulation of heat transfer during the process of welding with micro-jet cooling was implemented according to the initial boundary value problem (IBVP). The problem was solved employing the method of combining the finite difference method, Picard iterations, approximation by the radial basis function, and the fundamental solution method so as to solve the IVBP. The proposed method was validated by the data and results obtained during in situ experiments. The numerical approach enabled us to obtain variations in the temperature distribution values in HAZ with its different dimensional variants, ranging between 600 °C and 1400 °C.


Introduction
Mixed steel-nickel alloy joints have become more prevalent in civil engineering in light of the technological features of the components. Inconel-grade nickel alloys (nickel alloys with chromium, molybdenum, and iron) are highly weldable in shielding gases when utilizing conventional welding technologies, such as MIG (metal inert gas) and TIG (tungsten inner gas). Their good weldability is a result of their favorable austenitic microstructure, guaranteeing their acceptable ductility.
The Alloy 59 (NiCr23Mo16Al) is a superalloy used in structures intended for operations at elevated and high temperatures, as well as in aggressive environments [1,2]. This kind of material is particularly resistant to pitting and crevice corrosion [3] and, therefore, The practical aspect of the paper aims to provide the most suitable technology and its parameters for the manufacturing of a mixed joint with the best possible microstructure and mechanical properties. Therefore, we decided to verify the welding potential of a mixed joint made of Alloy 59 (nickel-chromium-molybdenum) and the non-alloy structural steel S355J2W using the MIG process, as well as various micro-jet cooling parameters. The presented results provided the basis for our conclusions regarding the selection of parameters to produce a mixed butt welded joint. Once the mixed joint was prepared with the use of the classic MIG process and the modified MIG process (with the application of micro-jet cooling), the microstructure and properties of the joint were determined.

Details on the Material and Welding Method
The selected materials were mixed welded joints made using the MIG process (131). The joints were made from 8 mm-thick sheets. The mechanical properties of the S355J2 steel in the as-received state, Alloy 59 [37], and the NiCr23Mo16 electrode wire ( Böhler, Germany,) weld metal (according to the EN ISO 18274 standard) are presented in Table 1. Table 1. Materials used to build the rotor platforms: mechanical properties [6].

Material
Yield The mentioned welding wire (superalloy from the alloy group 59) is used in highly aggressive environments. It is recommended for the joining of duplex and super duplex steels, stainless steels, and various nickel-based alloys [38]. Based on the data presented in Table 1, significant differences in the ultimate tensile strength and yield stress values of both joined materials can be concluded. The given properties result from the different chemical compositions of the used materials. Table 2 shows the chemical compositions of the S355J2 steel, Alloy 59, and weld metal of the electrode wire. Table 2. Materials used to build the rotor platforms: chemical composition [6]. Due to the purpose of the use of the materials outlined above and subjected to an extensive range of temperatures during welding, the thermal conductivity coefficient and specific heat capacity must be treated as a temperature function. The proper characteristics are given below: c p,Ni (T) = −5.739462 · 10 −10 T 4 + 1.59927 · 10 −6 T 3 −1.404298 · 10 −3 T 2 + 4.915666 · 10 −1 T + 40.95882 (1)
Concerning the different thermal conductivity coefficients of the two welded alloys, we decided to prepare joints with one-sided bevels of 1/2 V. The throat angle was at the level of 45 • . The beveling method of the steel sheet is presented in Figure 1. are the specific heat capacity (J and T denotes the temperature field (K). The subscripts mean the following: Nialloy, St-steel alloy, and w-weld.
Concerning the different thermal conductivity coefficients of the two welde we decided to prepare joints with one-sided bevels of 1/2V. The throat angle wa level of 45°. The beveling method of the steel sheet is presented in Figure 1.
The welding parameters of the joint were as follows: the electrode wire diam 1.2 mm, the arc voltage U = 21 V, and the welding current intensity was at the lev A ( Figure 2). The welded sheets had the dimensions of 800 × 200 × 8 mm. Arg selected to act as a shielding gas in the MIG process. The shielding gas flow rat the level of 15 l/min. The welding speed during the laying of all three stitches mm/min. The MIG welding method (131) in the downward position (PA) was cording to the requirements of the EN 15614-1 standard. The rotor joints were with a direct current with a positive polarity of the electrode.  The welding parameters of the joint were as follows: the electrode wire diameter was 1.2 mm, the arc voltage U = 21 V, and the welding current intensity was at the level of 137 A (Figure 2). The welded sheets had the dimensions of 800 × 200 × 8 mm. Argon was selected to act as a shielding gas in the MIG process. The shielding gas flow rate was at the level of 15 l/min. The welding speed during the laying of all three stitches was 240 mm/min. The MIG welding method (131) in the downward position (PA) was used according to the requirements of the EN 15614-1 standard. The rotor joints were welded with a direct current with a positive polarity of the electrode.  In this case, preheating was not used. Micro-jet cooling was applied directly welding. The weld metal deposit was prepared using the MIG welding method w without micro-jet cooling. Helium was selected as a micro-jet gas. The main par of the micro-jet cooling varied slightly, including the cooling steam diameter (50 60 µm) and gas pressure (0.4 MPa and 0.5 MPa).
The assembly of the welding head and the micro-jet injector is illustrated in F One of the obtained joints is presented in Figure 3. In this case, preheating was not used. Micro-jet cooling was applied directly after the welding. The weld metal deposit was prepared using the MIG welding method with and without micro-jet cooling. Helium was selected as a micro-jet gas. The main parameters of the micro-jet cooling varied slightly, including the cooling steam diameter (50 µm and 60 µm) and gas pressure (0.4 MPa and 0.5 MPa).
The assembly of the welding head and the micro-jet injector is illustrated in Figure 2. One of the obtained joints is presented in Figure 3.
welding. The weld metal deposit was prepared using the MIG welding method w without micro-jet cooling. Helium was selected as a micro-jet gas. The main par of the micro-jet cooling varied slightly, including the cooling steam diameter (50 60 µm) and gas pressure (0.4 MPa and 0.5 MPa).
The assembly of the welding head and the micro-jet injector is illustrated in F One of the obtained joints is presented in Figure 3.

Numerical Approach
The method of fundamental solutions is applied to simulate the heat transfer during the welding of two types of metals. The phenomenon is modeled as an boundary value problem, which is described by the system of non-linear partial d tial equations with initial conditions and non-linear boundary conditions. The no ity results from the non-linear thermal characteristics of the welded materials. The conductivity and heat capacity of the superalloy and steel are described by the n function of the temperature. It is necessary to consider these nonlinearities du range of temperatures in the welding and cooling process.

Problem Description and Mathematical Model
The transient temperature is considered in the case of the welded materials. A of the welded objects is presented in Figure 4. There are three regions: the weld w ing and the welded materials, the superalloy and steel. The cross-section shown in 5 is taken into account.
The geometry of these regions is described as follows. An area occupied by eralloy is defined as Ni Ω , the steel alloy is St Ω , and the area occupied by the w Ω , where dNi is the length of the horizontal rectangle edge of the area Ni Ω , thickness of the plate, dw is the width of the weld, and dst is the length of the low zontal edge of the steel zone.

Numerical Approach
The method of fundamental solutions is applied to simulate the heat transfer process during the welding of two types of metals. The phenomenon is modeled as an initialboundary value problem, which is described by the system of non-linear partial differential equations with initial conditions and non-linear boundary conditions. The nonlinearity results from the non-linear thermal characteristics of the welded materials. The thermal conductivity and heat capacity of the superalloy and steel are described by the nonlinear function of the temperature. It is necessary to consider these nonlinearities due to the range of temperatures in the welding and cooling process.

Problem Description and Mathematical Model
The transient temperature is considered in the case of the welded materials. A scheme of the welded objects is presented in Figure 4. There are three regions: the weld with cooling and the welded materials, the superalloy and steel. The cross-section shown in Figure 5 is taken into account.   The boundaries of the regions are indicated in Figure 6, and the boundary of t eralloy is , the boundary of the steel a 6 10 9 , and the boundary of the weld is  The geometry of these regions is described as follows. An area occupied by the superalloy is defined as Ω Ni , the steel alloy is Ω St , and the area occupied by the weld is Ω w , where d Ni is the length of the horizontal rectangle edge of the area Ω Ni , h is the thickness of the plate, d w is the width of the weld, and d st is the length of the lower horizontal edge of the steel zone.
The boundaries of the regions are indicated in Figure 6, and the boundary of the superalloy is Γ Ni = Γ 1 ∪ Γ 2 ∪ Γ 3 ∪ Γ 4 , the boundary of the steel alloy is Γ St = Γ 8 ∪ Γ 9 ∪ Γ 10 ∪ Γ 6 , and the boundary of the weld is The boundaries of the regions are indicated in Figure 6, and the bounda eralloy is , the boundary of the ste , and the boundary of the weld is 6 5 ~ ∪ Γ ∪ Γ = Γ w Figure 6. Geometry: nomenclature of areas and boundaries.
The equations of the transient heat transfer are introduced for each regi sideration. Thus, the equations are in the forms:

•
For the superalloy: For the steel alloy: For the weld: • For the superalloy: • For the steel alloy: • For the weld: where t denotes time, and x, y are geometrical coordinates. The other terminology is adequate for each considered area, noted by subscripts: Ni-for the superalloy, St-for the steel alloy, and w-for the weld. Additionally, the terminology is as follows: t Ni , t St , t wtemperature in the areas, ρ Ni , ρ St , ρ w -mass density, c p, The values of the physical parameters applied in the numerical approach at the areas examined are listed in Equations (1)- (6). The problem is defined for the chosen time interval, i.e., t ∈ [t 0 , ∞), and t 0 is the initial time. The associated initial conditions describe the initial temperature of the superalloy and steel alloys, which is the ambient temperature: t Ni t St and the temperature of weld, which is the temperature of the welding: Materials 2022, 15, 8579 8 of 28 wherex = ( x, y), f Ni ( x, y) describes the initial temperature of the superalloy, f St ( x, y) describes the initial temperature of the steel alloy, and f w ( x, y) describes the initial temperature of the welding element. The problem is determined by the boundary conditions. The free surfaces of the superalloy-steel alloys are subjected to natural convection: where a Ni and a St are the superalloy and steel convection coefficients, respectively, and t am is the ambient temperature.
Regarding the boundaries Γ 4 and Γ 9 , it is assumed that the boundaries are far from the HAZ, and there is a constant temperature equal to the ambient temperature: On the outer boundaries of the weld, the cooling process is applied. Thus, the boundary condition is given by the function, which describes the cooling as follows: It is assumed that there is perfect contact between the contacted regions. This means that the following boundary conditions should be fulfilled: where λ Ni (t Ni ), λ St (t St ), and λ w (t w ) are the thermal conductivity of the superalloy, steel, and weld, respectively. Due to the range of temperatures that appear in the welding process, the material characteristics, such as the specific heat capacity and thermal conductivity coefficient, for every material used are described by Equations (1)- (6).
To avoid numerical errors, it is convenient to consider the problem in a non-dimensional form. The dimensionless variables are introduced as follows: where t max is the maximum temperature of the weld, t am is the ambient temperature, λ 0 is the reference thermal conductivity coefficient, and c p,0 is the reference specific heat capacity. The area given in the dimensionless form and occupied by the superalloy is Ω Ni , the steel alloy is Ω St , and the dimensionless area occupied by the weld is Ω w . In the dimensionless form, the boundaries of the areas are as follows: Γ 1 , Γ 2 , Γ 3 , Γ 4 , Γ 5 , Γ 6 ,Γ 7 , Γ 8 , Γ 9 , Γ 10 .
Therefore, the boundary of the region Ω Ni is Γ Ni = Γ 1 ∪ Γ 2 ∪ Γ 3 ∪ Γ 4 , the boundary of the region Ω St is Γ St = Γ 8 ∪ Γ 9 ∪ Γ 10 ∪ Γ 6 , and the boundary of the region Ω w is described Then, the governing Equations (7)-(9) have the dimensionless form in the superalloy, steel, and the weld areas: where the dimensionless forms of the specific heat capacity and thermal conductivities are: The initial conditions have the following form.
The associated initial conditions describe the dimensionless initial temperature of the superalloy and steel alloys, which is the ambient temperature (in the dimensionless form): And the temperature of the welding: where x = (x, y) and τ 0 is the initial non-dimensional time.
The boundary conditions have the following dimensionless form: T St = T w on Γ 6 (38) To solve the initial boundary problem given by Equations (24)- (39), the numerical algorithm is proposed.

Numerical Procedure
The main assumption for the numerical procedure is that one must apply one of the mesh-free methods: the method of fundamental solutions (MFS). However, for the initial problem, an auxiliary method should be used to treat the differential temperature function according to time. This paper proposes approximating the differential concerning the time by the finite difference (FD). To implement the FD, the time interval, τ ∈ [τ 0 , τ max ], in which the solution is sought, is divided into n t subintervals. The solutions are obtained in certain time steps: where dτ = τ max −τ 0 n t and i = 1, . . . , n t . The partial derivative with respect to time is approximated by: For i = 1, the initial conditions (27)-(29) are defined as: Applying the FD, the boundary value problem (BVP) at each time step is determined. The obtained BVPs are described by non-linear equations and the boundary conditions (conditions (37) and (39) are non-linear). We propose the implementation of Picard iterations, and we iteratively rewrite the governing equations: w,i (x, y) denote the temperatures of the superalloy, steel alloys, and weld at the i-th time step and j-th iteration, and: The boundary conditions (30)-(39) are still valid and are considered in the iteration process. Therefore, for each time step and iteration, the BVP i given by Equations (46)-(48) and (30)-(39) must be solved.
The proposal of the paper is to apply one of the meshless methods, i.e., the method of fundamental solutions supported by approximation by the radial basis functions (RBF) and monomials.
It is assumed that the solutions of the BVP described by the inhomogeneous Equations (46)-(48)and the boundary conditions (30)- (39) are the sums of both the particular and homogeneous solutions: p,w,i are the particular solutions, and T h,w,i are the homogeneous solutions of the superalloy, steel, and weld, respectively, at the i-th time step and j-th iteration.
The particular solutions are obtained using approximation by the RBFs and monomials. The right-hand side of Equations (46)-(48) (inhomogeneous part) is approximated as follows: w,k (x, y) (for k = 1, . . . , Na) are the RBFs defined for the superalloy steel alloy and the weld, respectively, p k (x, y) (for k = 1, . . . , Nm) are monomials, a Ni,k , a St,k , and a w,k (for k = 1, . . . , Na + Nm) are real numbers, Na is the number of approximation points in each area Ω Ni , Ω St , and Ω w , and Nm is the number of monomials. The set of approximation points is defined in each area: w,k ⊂ Ω w (for k = 1, . . . , Na), as presented in Figure 7. The w,k (x, y) are defined as follows:

approximation points in each area
Ni Ω , St Ω , and w Ω , and Nm is the number of monomials. The set of approximation points is defined in each area: (for k = 1,…,Na), as presented in Figure 7 The next step of the proposed procedure is the calculation of the homogeneous solution, which is obtained by applying the MFS. It is assumed that the homogeneous solutions are in the form of: The system of the linear algebraic Equations (52)-(54) is written for each approximation point (the proper domain) and solved. The solutions are the coefficients a Ni,k , a St,k , and a w,k (for k = 1, . . . , Na + Nm). Thus, the particular solutions are in the form of: where ψ(r) and P k (x, y) (for k = 1, . . . , Nm) are the particular solutions of Poisson's equation with the inhomogeneous part of φ(r) and p k (x, y) (for k = 1, . . . , Nm), respectively. The next step of the proposed procedure is the calculation of the homogeneous solution, which is obtained by applying the MFS. It is assumed that the homogeneous solutions are in the form of: are the sets of source points placed outside the regions Ω Ni , Ω St , and Ω w , respectively, and Ns is the number of source points for each area, s is the distance between the boundaries Γ Ni , Γ St , and Γ w , and the fictitious boundaries with the source points c Ni,k , c St,k , and c w,k (for k = 1, . . . , Ns) are real numbers. The forms of solutions (49)-(51) are introduced to the boundary conditions (30)- (39), and the approximated form of the homogeneous solutions (58)-(60) is included. Thus, the boundaries Γ Ni , Γ St , and Γ w are discretized, i.e., the sets of boundary points are defined as: Nb is the number of boundary points at each boundary. Additionally, for each boundary point, the proper boundary condition is written. This gives the system of linear algebraic equations with the unknowns c Ni,k , c St,k , and c w,k (for k = 1, . . . , Ns). The solution of this equations system completes the procedure of finding the solution in the i-th time step and j-th iteration.
The idea of the selection of boundary and source points is presented in Figure 8. The black points are the boundary points. The black stars denote the boundary of the superalloy, the black ovals lie on the weld boundary, and the black crosses denote the boundary steel zone. The source points are also plotted in Figure 8. The blue stars denote the fictitious contour of the superalloy, the red ovals lie on the fictitious weld shape, and the green crosses denote the fictitious contour of the steel zone. The distance between the boundary and fictitious contour s is clearly indicated. Therefore, for this case the following sentence can be formulated: different patterns of the points reflects sections of single bevel joint.
, respectively, and Ns is the number of source points for each area, s is the distance between the boundaries Ni Γ , and w Γ are discretized, i.e., the sets of boundary points are defined as: Nb is the number of boundary points at each boundary. Additionally, for each boundary point, the proper boundary condition is written. This gives the system of linear algebraic equations with the unknowns The idea of the selection of boundary and source points is presented in Figure 8. The black points are the boundary points. The black stars denote the boundary of the superalloy, the black ovals lie on the weld boundary, and the black crosses denote the boundary steel zone. The source points are also plotted in Figure 8. The blue stars denote the fictitious contour of the superalloy, the red ovals lie on the fictitious weld shape, and the green crosses denote the fictitious contour of the steel zone. The distance between the boundary and fictitious contour s is clearly indicated. Therefore, for this case the following sentence can be formulated: different patterns of the points reflects sections of single bevel joint. Figure 8. The sets of boundary and source points: the black points-the boundary points, black stars-the boundary region for superalloy, the black ovals-the weld boundary, the black crossessteel boundary zone, the blue stars -the fictitious contour of the superalloy, the red ovals-the fictitious shape of the weld, the green crosses-the fictitious contour of the steel.
The procedure of the iterations at the i-th time step is stopped if the following condition is fulfilled: Figure 8. The sets of boundary and source points: the black points-the boundary points, black stars-the boundary region for superalloy, the black ovals-the weld boundary, the black crosses-steel boundary zone, the blue stars-the fictitious contour of the superalloy, the red ovals-the fictitious shape of the weld, the green crosses-the fictitious contour of the steel.
The procedure of the iterations at the i-th time step is stopped if the following condition is fulfilled: where where w,k ⊂ Ω w ∪ Γ w (for k = 1, . . . , Nt) are the sets of trial points placed in each area, and on its boundary, Nt is the number of trial points in each area and boundary, and ε is a small number.

Validation of the Numerical Procedure and Numerical Results
Computer simulations were carried out using a computer program written by the authors of the article. The numerical calculation stage was performed to validate the proposed numerical approach. The proper values of the numerical procedure parameters were chosen to obtain the demanded results. These results of the numerical simulations were compared with the real data. Then, the different values of the parameters of the cooling procedure were tested. The aim was to obtain as a heat-affected zone that was as small as possible.
The dimensions of the specimens were expressed as the following values: h = 8 mm, d Ni = 30 mm, d w = 2 mm, and d St = 30 mm.
The cooling process was performed using a micro-jet system. The values of the welding velocity and the cooling process were equal to 240 mm/min. The gas used was helium. The other parameters of the cooling system were the cooling jet velocity v (related to the pressure of the gas) and the diameter of the micro-jet injector d j . Thus, the four variants of the cooling system were tested.
These For the application of the MFS, the special function, called the fundamental solution, for the differential operator in the governing equations has to be known. The same operator, the Laplace one, is applied in all the governing Equations (46)-(48). The fundamental solution for the Laplace operator is given by the formula: where r is the distance from one of the source points and the chosen point that lies in the domain in which the equations are fulfilled.
To apply the approximation given by Equations (52)-(54), the radial basis functions must be chosen. In this numerical implementation, the thin-plate spline RBF is chosen, and seven monomials are used. The thin-plate spline is expressed as: The parameters of the numerical method were selected. The selection was based on some simulations using the trial values of the method parameters (number of boundary points, number of source points, distance between the real and source contours). This calibration procedure was performed for the cooling variant D. Ten time steps were calculated, and at each time step, some iterations were performed. The number of iterations depends on the criterion's fulfillment (61).
The validation procedure consists of simulations using different values of the method parameters and checking the mean square error applied to the boundary points. The fulfillment of the boundary conditions by the obtained solution was tested for every iteration at each time step. As an example of the values, the mean square error for the 5th time step and 3rd iteration is presented in Table 3. For every calculated time step and iteration, the error had the lowest value for the number of boundary points equal to 40, and the distance between the boundary and fictitious boundary was equal to 0.2. The parameter values were selected if the mean square error had the smallest value.
The trial calculations were also performed to select the number of approximation points in the regions (superalloy, steel, and weld). This selection is conducted by calculating a mean square error of the fulfillment of the right-hand-side functions of Equations (46)-(48) by approximations (52)-(54). The values of the mean square error of the approximation error for the 5th time step and 3rd iteration are presented in Table 4. As we can observe, the lowest value appears for the number of approximation points equal to 121. The other values of the mean square error calculated for the rest of the time steps and iterations were lowest when 121 approximation points were defined.  Mean square error 4 × 10 −4 3 × 10 −6 4 × 10 −9 5 × 10 −9 9 × 10 −8 Finally, the following parameters were chosen: the number of boundary points was 40 (in each considered region); the number of source points was 121 (in each considered area); and the distance between the real and fictitious boundary was 0.2 (in each considered area).
For this referential example, the values of (61)-(64) were also calculated for the set of trial points. Thus, it was observed that values of (61)-(64) of less than 10 −5 were achieved for a relatively low number of iterations (see Table 5). The number of calculated iterations is presented in Table 6, which was prepared for the first two time steps. As concluded in Table 6, the proposed numerical procedure requires little work.   After determining the parameters of the numerical procedure, the numerical experiment can be started. The range of the heat-affected zone of the superalloy and steel alloys was simulated. The calculations were made for the welding parameters given in the natural experiment.
Thus, the temperature distribution during the welding and cooling was calculated for all the variants of the cooling process. The distribution of the temperature with the cooling variant D is presented in Figure 9a for the superalloy and Figure 9b for the steel alloy. The graphs were created for the time step in which the maximum cooling appeared. ment can be started. The range of the heat-affected zone of the superalloy and steel alloys was simulated. The calculations were made for the welding parameters given in the natural experiment.
Thus, the temperature distribution during the welding and cooling was calculated for all the variants of the cooling process. The distribution of the temperature with the cooling variant D is presented in Figure 9a for the superalloy and Figure 9b for the steel alloy. The graphs were created for the time step in which the maximum cooling appeared. We can also observe changes in the temperature over time. The graphs presented in Figure 10 show the temperature at the chosen time steps, obtained by applying different variants of the cooling system to the steel region. The graphs ( Figure 11) show the temperature at the selected time steps obtained using another variant of the cooling system in the superalloy region. We can also observe changes in the temperature over time. The graphs presented in Figure 10 show the temperature at the chosen time steps, obtained by applying different variants of the cooling system to the steel region. The graphs ( Figure 11) show the temperature at the selected time steps obtained using another variant of the cooling system in the superalloy region.   In every graph in Figures 10 and 11, we can observe that the boundary condition at the free end is fulfilled, i.e., the temperature should be equal to the ambient temperature for x = 0 in the superalloy and for x = dNi + dw + dSt in the steel zone. At the other ends of In every graph in Figures 10 and 11, we can observe that the boundary condition at the free end is fulfilled, i.e., the temperature should be equal to the ambient temperature for x = 0 in the superalloy and for x = d Ni + d w + d St in the steel zone. At the other ends of both areas, the temperature is dependent on the cooling of the weld area, and it is also straightforward to observe that these boundary conditions are fulfilled.
From the engineering point of view, it is interesting to estimate the heat-affected zone of welded materials. The process of cooling and decreasing the HAZ has an important role, and the rate of cooling determines the welded materials' characteristics. Observing Figures 10 and 11, we found changes in the range of the HAZ in the superalloy and steel zones. The temperature limit of 773 K was considered to examine the range of the heat influence zone. The content of the HAZ for the chosen time steps is presented in Table 7. The numbers presented in Table 7 are the length of the zone with a temperature higher than 773 K divided by the vertical dimension of the specimens. We can observe that the HAZ decreases fastest with cooling variant D, and the slightest decrease in the HAZ appeared when cooling variant A was applied. Additionally, variants B and C give similar decreasing rates of the HAZ, while variant B yields a smoother cooling.
It is most interesting to simulate and compare the processes for all the four variants of the cooling after welding. The simulation results are presented in Figure 12a for the superalloy, Figure 12b for the steel alloy). It shows the temperature distribution at specific points (the points are placed in the HAZ, 0.5 mm from the welded and cooled boundary of each welded material, i.e., the superalloy and steel alloys). Table 7 are the length of the zone with a temperature higher than 773 K divided by the vertical dimension of the specimens. We can observe that the HAZ decreases fastest with cooling variant D, and the slightest decrease in the HAZ appeared when cooling variant A was applied. Additionally, variants B and C give similar decreasing rates of the HAZ, while variant B yields a smoother cooling.

The numbers presented in
It is most interesting to simulate and compare the processes for all the four variants of the cooling after welding. The simulation results are presented in Figure 12a for the superalloy, Figure 12b for the steel alloy). It shows the temperature distribution at specific points (the points are placed in the HAZ, 0.5 mm from the welded and cooled boundary of each welded material, i.e., the superalloy and steel alloys). We can observe that the most intensive cooling is obtained by variant D for both metal regions, but the heating and cooling progress faster in the steel alloy compared to the superalloy region.
The parameter known in the literature which demonstrates the efficiency of the cooling system is the cooling rate. The computer simulations make it possible to obtain the precise weld cooling conditions of the welding process, especially in the range of 500-800 °C. The cooling rate parameter is presented in Table 8. One can observe the cooling rate according to the four variants of the nozzle diameter and gas pressure in each welded material: the steel and superalloy. The highest values are those for the cooling by the D variant. Two of them, B and C, give similar results. The slowest cooling appears for variant A. However, the cooling process should not be too intensive. It is crucial to note that the We can observe that the most intensive cooling is obtained by variant D for both metal regions, but the heating and cooling progress faster in the steel alloy compared to the superalloy region.
The parameter known in the literature which demonstrates the efficiency of the cooling system is the cooling rate. The computer simulations make it possible to obtain the precise weld cooling conditions of the welding process, especially in the range of 500-800 • C. The cooling rate parameter is presented in Table 8. One can observe the cooling rate according to the four variants of the nozzle diameter and gas pressure in each welded material: the steel and superalloy. The highest values are those for the cooling by the D variant. Two of them, B and C, give similar results. The slowest cooling appears for variant A. However, the cooling process should not be too intensive. It is crucial to note that the cooling rate of the steel region is lower than that of the superalloy zone. This means that the cooling process is more intensive in the superalloy region. Therefore, the parameters of the micro-jet cooling process should be selected based on the steel region in order to obtain the intensity of the cooling, which gives the demanded HAZ.

Comparison with Another Numerical Approach
In engineering simulations, the finite element method (FEM) is often used. The considered problem is given by Equations (7)- (22), which were solved using the software based on the FEM. Computer simulations based on the FEM were carried out using a computer program written by the authors of the article. It was effortless to introduce the geometry of the specimens. It required some effort to define the temperature-dependent material characteristics and proper initial boundary conditions and coefficients in the governing equations. Using FEM software computations, the domains of the materials were divided into a large number of elements. Triangular elements with three nodes were used in the calculations. Each domain was divided into 10,000 elements. This caused the system of algebraic equations to be extensive (the number of equations was proportional to the number of 3 × 10,000 nodes in the base functions, i.e., at least 90,000 equations) compared to what we observed in the case of the MFS. Using the proposed algorithm, two systems of algebraic equations were solved. One method of linear algebraic Equations (52)-(54) was used for the calculation of the particular solutions. The number of equations was3 × (121 + 7) = 384. The second set of equations were used for the calculation of the homogeneous solutions (58)-(60), which consisted of 100 equations. Due to such a low number of equations in the system, the proposed algorithm based on MFS is less time-consuming than the calculations performed using FEM. MFS is a meshless method, and RBF does not require the construction of any mesh in the considered region. Only a cloud of points is defined in the considered areas, and only the boundaries of the areas are subject to discretization (the dimension of the problem is reduced by 1). However, in FEM, the procedure of creating the mesh, defining the element, and the basis functions must be applied. This procedure is also highly time-consuming. Moreover, the MFS approach requires less memory (for the data collection and for the processor to perform calculations) than FEM. The values of the other parameters of the procedure based on FEM, analogous to the parameters of the proposed procedure based on MFS, were assumed to be the same in both simulations. The parameters, having been adopted in this way, allowed us to obtain numerically stable and convergent procedures. Both methods deliver results in the form of continuous analytical functions. This is a very useful property, because it allows for the performance of post-processing analysis. The possibility of such an analysis is significant from an engineering point of view. The advantage of the MFS approach is that the analytical functions, which construct the solutions, have derivatives of the order of infinity. In FEM, this aspect depends on the used basis functions, but in the general approach, the basis functions possess derivatives of a few first orders.
A comparison of the results obtained using both procedures, the MFS-based and FEM-based methods, is presented in Table 9. The calculations were performed for variant C of the cooling process. The values of the temperature calculated at the point of half of the specimen height and 0.5 cm from the edge of the welded edges (for the superalloy and steel zones) are provided. As can be observed, the differences between the calculated values using the different numerical procedures are not substantial. However, the comparison of the time and memory consumption using both procedures has to be taken into account as well. In our computations, the algorithm based on MFS is more appropriate. It gives the same results as FEM but is less time-and memory-consuming. The other method used to solve initial boundary engineering problems is the finite differences method (FDM). This method partially supported the algorithm proposed in this paper. The basic form of the FDM is dedicated to linear problems. For the purpose of using the FDM to solve all the aspects of the considered problems of welding and cooling, the other numerical methods should be implemented as well. However, there is a fundamental imperfection of the FDM. This method gives results in a discrete form and not a continuous analytical function. Thus, additional procedures have to be applied to perform the engineering post-processing analysis. For this purpose, the interpolation technique is implemented. However, such an approach does not ensure the compatibility of the used interpolation function with the character of the calculated phenomenon. Due to this shortcoming of the FDM, the calculations performed by this method were omitted.

The Application of the MFS-Based Procedure
In the previous section, it was validated that the proposed MFS-based numerical algorithm is a good tool for simulating the time-dependent temperature field in welded materials. Thus, the method may be applied to estimate the cooling process parameters so as to obtain workpieces with the demanded characteristics (according to the analysis of the HAZ or cooling rate). It is also possible to simulate the time-dependent welding-cooling process for other materials. The procedure was elaborated for the welding of two different materials with non-linear characteristics dependent on the temperature.
The simulation test was performed using materials with thermal conductivity coefficients and a specific heat capacity different from those of the materials defined in the validation section.
Thus, it is assumed that the superalloy is replaced by materials with a thermal conductivity coefficient and specific heat capacity equal to the characteristics of the alloy, multiplied by 0.8, 0.9, 1.1, and 1.2. This means that: where the first subscript, which is equal to I, refers to the material replacing the superalloy, and the second one describes the version of the material.
The parallel change is performed for the steel, replaced by the material with notation II. Additionally, the characteristics of the variant of material II are noted as: The parameters of the applied numerical procedure (including the parameters related to the welding-cooling process) are not changed. They stay the same as in the validation procedure.
For such modified materials, the values of the cooling rate parameter are calculated, and they are included in Table 10.  1  47  65  67  160  2  49  67  69  166  3  52  72  73  180  4  55  75  77  189   Material II   1  43  54  55  135  2  45  57  58  137  3  50  61  62  152  4  53  64  65  152 The cooling rate parameter has values in the range of 50-80 (which is most popularly demanded by engineers) for cooling variants B and C. The materials' characteristics are changed slightly compared to the input ones (given by Formulas (1), (2), (4), and (5)). Thus, we can also conclude that cooling schemes B and C are the most advantageous. More concerning is the fact that the cooling rate of variant D becomes very high (higher than 180). This may indicate that cracks could appear in the welded material. In Section 3.3, a high cooling rate appeared, and in chapter 2, it was pointed out that there were no cracks. However, in the case of this section, regarding the cooling, the rate is much higher, and there is no reliable assurance that imperfections (such as cracks) will not appear.

Validation In Situ and Discussion
To validate the results obtained by the numerical simulations, some mechanical tests were performed.
Once the joints were welded with various micro-jet cooling parameters, non-destructive tests were carried out, including visual, penetration, and ultrasonic tests.
The following results were collected: • Minor cracks appeared in the joints made using the MIG process, without micro-jet cooling (Figure 13a,b); • No defects or non-conformities for level B (according to the PN-EN ISO 5817:2005 standard) appeared in the joints made using the MIG process, where helium micro-jet cooling was applied.
changed slightly compared to the input ones (given by Formulas (1), (2), (4), and (5)). Thus, we can also conclude that cooling schemes B and C are the most advantageous. More concerning is the fact that the cooling rate of variant D becomes very high (higher than 180). This may indicate that cracks could appear in the welded material. In Section 3.3, a high cooling rate appeared, and in chapter 2, it was pointed out that there were no cracks. However, in the case of this section, regarding the cooling, the rate is much higher, and there is no reliable assurance that imperfections (such as cracks) will not appear.

Validation In Situ and Discussion
To validate the results obtained by the numerical simulations, some mechanical tests were performed.
Once the joints were welded with various micro-jet cooling parameters, non-destructive tests were carried out, including visual, penetration, and ultrasonic tests.
The following results were collected: • Minor cracks appeared in the joints made using the MIG process, without micro-jet cooling (Figure 13a  The next step for the joint quality assessment was performed by applying destructive tests (bending, tensile, hardness, microscopic observations). In the case of all the specimens, the bending test was performed. Five bending tests were conducted from the face and the root side of the weld. In the joints made using the classic MIG process, cracks in the weld were observed at a bending angle above 140 • (Figure 14a). The next step for the joint quality assessment was performed by applying destructive tests (bending, tensile, hardness, microscopic observations). In the case of all the specimens, the bending test was performed. Five bending tests were conducted from the face and the root side of the weld. In the joints made using the classic MIG process, cracks in the weld were observed at a bending angle above 140° (Figure 14a).
In the case of the application of micro-jet cooling in the MIG process of connecting Alloy 59 with the S355J2 hard-rusting steel, after the welding, no cracks or other incompatibilities were detected in the tested specimens ( Figure 14b) The bending test results showed that the joints were welded correctly, reaching the required quality.
(a) (b) Figure 14. View of the welding specimens after the bending test: (a) the specimen made without micro-jet cooling, (b) the specimen made with micro-jet cooling.
Next, tensile tests were carried out. The obtained results for the mixed joint made of Alloy 59 and hard-rusting S355J2 steel prepared using MIG and different micro-jet cooling parameters (average from three measurements) are presented in Table 11, Figure 15.  In the case of the application of micro-jet cooling in the MIG process of connecting Alloy 59 with the S355J2 hard-rusting steel, after the welding, no cracks or other incompatibilities were detected in the tested specimens (Figure 14b) The bending test results showed that the joints were welded correctly, reaching the required quality.
Next, tensile tests were carried out. The obtained results for the mixed joint made of Alloy 59 and hard-rusting S355J2 steel prepared using MIG and different micro-jet cooling parameters (average from three measurements) are presented in Table 11, Figure 15.   Table 11.
The data analysis (Table 11, Figure 15) shows that all the joints made with the use of cooling had comparable mechanical properties (a yield stress above the desired value of 300 MPa and ultimate tensile strength above 500 MPa). These data follow the mechanical parameters of the S355J2W defined by the material manufacturers [39]. In the connections for which micro-jet cooling was not introduced, worse mechanical properties were noted. These could have been influenced by micro-cracks that did not exist in the joints made using micro-jet cooling. The best mechanical properties of the mixed (Alloy 59 with S355J2 steel) joint were obtained using the following micro-jet cooling parameters: for the helium cooling gas, a gas pressure of 0.5 MPa and a stream diameter of 50 µm.
Hardness tests of the Alloy 59 joint with S355J2 hard-rusting steel were also used to extend our knowledge of the joint quality. The hardness of the steel was 185 MPa, while in the case of the superalloy, this parameter reached 375 MPa. Table 12 and Figure 16 show the results of the hardness tests of the heat-affected zone (HAZ) for both the welded sides and the weld hardness of the tested joints (using different micro-jet cooling parameters). Brinell hardness data calculated based on the Vickers hardness results for the HAZ were utilized to assess the values of the ultimate tensile strength for the three areas considered, as shown in Table 12, Figure 17. For this purpose, the factor equal to 3.53 as a proportion between the UTS and the hardness value was determined. In this case, the smallest value of the mechanical parameter was equal to 745 MPa for the joint region represented by the S355J2 steel. The highest value was obtained for the superalloy side, expressing 1197 MPa. Following these data, it can be determined that the hardness values become lower from the superalloy side and through the joint material up to the steel side. From a practical point of view, the steel and superalloy regions can be called the strongest and weakest ones, respectively. In regard to inspection, this means that the steel region is representative of the observation of cracks and deformation.   Table 11.
The data analysis (Table 11, Figure 15) shows that all the joints made with the use of cooling had comparable mechanical properties (a yield stress above the desired value of 300 MPa and ultimate tensile strength above 500 MPa). These data follow the mechanical parameters of the S355J2W defined by the material manufacturers [39]. In the connections for which micro-jet cooling was not introduced, worse mechanical properties were noted. These could have been influenced by micro-cracks that did not exist in the joints made using micro-jet cooling. The best mechanical properties of the mixed (Alloy 59 with S355J2 steel) joint were obtained using the following micro-jet cooling parameters: for the helium cooling gas, a gas pressure of 0.5 MPa and a stream diameter of 50 µm.
Hardness tests of the Alloy 59 joint with S355J2 hard-rusting steel were also used to extend our knowledge of the joint quality. The hardness of the steel was 185 MPa, while in the case of the superalloy, this parameter reached 375 MPa. Table 12 and Figure 16 show the results of the hardness tests of the heat-affected zone (HAZ) for both the welded sides and the weld hardness of the tested joints (using different micro-jet cooling parameters). Brinell hardness data calculated based on the Vickers hardness results for the HAZ were utilized to assess the values of the ultimate tensile strength for the three areas considered, as shown in Table 12, Figure 17. For this purpose, the factor equal to 3.53 as a proportion between the UTS and the hardness value was determined. In this case, the smallest value of the mechanical parameter was equal to 745 MPa for the joint region represented by the S355J2 steel. The highest value was obtained for the superalloy side, expressing 1197 MPa. Following these data, it can be determined that the hardness values become lower from the superalloy side and through the joint material up to the steel side. From a practical point of view, the steel and superalloy regions can be called the strongest and weakest ones,   Table 12.  Table 12.
The data presented in the table show that the hardness values of all the joints made using micro-jet cooling are comparable. The joints made without micro-jet cooling are characterized by the slightly lower hardness of the joint in both heat-affected zones. Figure 18 illustrates the microstructure of the mixed weld made using the MIG process with the application of micro-jet cooling.
The weld manufactured with the application of micro-jet cooling has an austenite dendritic cell microstructure. This guarantees better plastic properties, derived from the column microstructure. No cracks were observed in the junction with the dendritic cellular structure. The influence of the cooling rate on the structure of the welded material was considered in [40]. This is confirmed by the results obtained from the tensile test, observing the number of data collected upon fracture, the ductility value (Figure 19a,b, respectively), and the fracture zone (Figure 20b,c).  Table 12. Figure 17. Variations in the ultimate tensile strength as an effect of the parameters of the micro-jet cooling based on the data presented in Table 12.  Table 12. The data presented in the table show that the hardness values of all the joints made using micro-jet cooling are comparable. The joints made without micro-jet cooling are characterized by the slightly lower hardness of the joint in both heat-affected zones. Figure 18 illustrates the microstructure of the mixed weld made using the MIG process with the application of micro-jet cooling.
The weld manufactured with the application of micro-jet cooling has an austenite dendritic cell microstructure. This guarantees better plastic properties, derived from the column microstructure. No cracks were observed in the junction with the dendritic cellular structure. The influence of the cooling rate on the structure of the welded material was considered in [40]. This is confirmed by the results obtained from the tensile test, observing the number of data collected upon fracture, the ductility value (Figure 19a,b, respectively), and the fracture zone (Figure 20b,c).   The data presented in the table show that the hardness values of all the joints made using micro-jet cooling are comparable. The joints made without micro-jet cooling are characterized by the slightly lower hardness of the joint in both heat-affected zones. Figure 18 illustrates the microstructure of the mixed weld made using the MIG process with the application of micro-jet cooling.
The weld manufactured with the application of micro-jet cooling has an austenite dendritic cell microstructure. This guarantees better plastic properties, derived from the column microstructure. No cracks were observed in the junction with the dendritic cellular structure. The influence of the cooling rate on the structure of the welded material was considered in [40]. This is confirmed by the results obtained from the tensile test, observing the number of data collected upon fracture, the ductility value (Figure 19a,b, respectively), and the fracture zone (Figure 20b,c).  These results directly demonstrate that the weld's fracture under monotonic te occurred upon no dynamic reaction of the material examined, compared to the stru steel. It can be concluded the MIG process, when applied to the Alloy 59 and S35 joint, is very beneficial, as there is no brittle cracking and no rapid development of c before the separation.

Summary and Conclusions
The welding of two types of structural materials (Alloy 59 and S355J2W stru steel) is a challenging task, but in an in situ experiment using a micro-jet cooling pr the problem can be solved. Here, the parameters of the welding and cooling process discussed, and the results of the in situ experiments were presented and analyzed data and the experimental results are the basic steps for the numerical stage. Based o conception of the numerical approach, the method for the calculation can be deem innovative. The proper combination of the FDM, Picard iterations, approximation by and MFS was implemented to solve the heat transfer problem during the welding-co process. The validation of the proposed numerical method was performed accord the results of the in situ experiments. The theoretical research, performed based o proposed numerical algorithm, allows for the estimation of the characteristics of t sultant specimen. It is also possible to select welding-cooling process parameters to o workpieces with the required values of the mechanical properties.
Analyzing the results collected in the non-destructive and destructive tests, th lowing conclusions can be formulated: • Micro-jet cooling used in the MIG process enables one to obtain a correctly m joint, such as that of the S355J2 steel and Alloy 59. • Excellent mechanical properties of the joint were achieved, i.e., a high yield (248 MPa) and ultimate tensile strength (518 MPa, ductility 10%). These results directly demonstrate that the weld's fracture under monotonic tension occurred upon no dynamic reaction of the material examined, compared to the structural steel. It can be concluded the MIG process, when applied to the Alloy 59 and S355J2W joint, is very beneficial, as there is no brittle cracking and no rapid development of cracks before the separation.

Summary and Conclusions
The welding of two types of structural materials (Alloy 59 and S355J2W structural steel) is a challenging task, but in an in situ experiment using a micro-jet cooling process, the problem can be solved. Here, the parameters of the welding and cooling process were discussed, and the results of the in situ experiments were presented and analyzed. The data and the experimental results are the basic steps for the numerical stage. Based on the conception of the numerical approach, the method for the calculation can be deemed as innovative. The proper combination of the FDM, Picard iterations, approximation by RBS, and MFS was implemented to solve the heat transfer problem during the welding-cooling process. The validation of the proposed numerical method was performed according to the results of the in situ experiments. The theoretical research, performed based on the proposed numerical algorithm, allows for the estimation of the characteristics of the resultant specimen. It is also possible to select welding-cooling process parameters to obtain workpieces with the required values of the mechanical properties.
Analyzing the results collected in the non-destructive and destructive tests, the following conclusions can be formulated:

•
Micro-jet cooling used in the MIG process enables one to obtain a correctly mixed joint, such as that of the S355J2 steel and Alloy 59. • Excellent mechanical properties of the joint were achieved, i.e., a high yield stress (248 MPa) and ultimate tensile strength (518 MPa, ductility 10%).

•
The inspection confirmed no welding defects or incompatibilities, such as sticking, fusion, cracks, or porosity.

•
The in situ experiments showed that the best properties of the welded material can be obtained by welding with cooling variant C. The values of the cooling rate (60 for the superalloy and 71 for the steel zone) of variant C (presented in Table 9) are the most beneficial from the engineering point of view. It is known that, for such values of the cooling rate, the welding process gives the best-quality joints. There are no cracks, and the obtained materials have good mechanical and plastic properties. This was confirmed by our in situ experiments.

•
The numerical algorithm based on the MFS is sufficient, because it requires less time and PC memory than the other methods (i.e., using meshing and the discrete approach).

•
The method can also be used for other types of joints made of metallic structural materials, including hard-welded ones, such as high-strength steels and their advanced grades.

•
The algorithm can be used in the welding calibration process when two different structural materials are taken for joining. It is also very suitable for training courses for welding engineers.
The applied calculation method is offers new possibilities for the use of equations based on a free approach and is closely related to the physics of the phenomena associated with the manufacturing of joints by welding using the micro-jet cooling technique as well.
Observing the results of the numerical experiment, it is acceptable to state the proposed method based on the method of fundamental solutions is the proper tool for solving the problem of multi-material welding with micro-jet cooling.
The proposed method of joining steel-nickel alloy can be tested in situ and numerically, showing that this material can be used for the critical elements of power plant devices.