Electromagnetic Modeling and Thermal Analysis of a Non-Axisymmetric System for Induction Brazing

: The electromagnetic process happening between an induction brazing coil and the brazing assembly was studied and is described in the following pages. Using the tools of numerical analysis, we managed to achieve 3-D simulation of the process, which includes presentation of the magnetic field, thermal heat transfer, and time dependence. The results presented are considered for two of commonly used materials—copper and stainless steel in the shape of pipes. The induction brazing coil is a complex C-shape to match with the real industrial practices. The results obtained from the model are in accordance with the experimental results, as in both the model and the experiment, the required temperature was reached in about 6 s. It was found that during the brazing process, the inductance of the winding increases by about 4 nH and the resistance by 1.5 m Ω , which is important for the coordination of the power supply.


Introduction
The joining of two metal materials by the melting and adding of a filler metal material is known as brazing. The process is performed at temperatures that range from 450 • C to more than 1000 • C depending on the certain case. Achieving such temperatures requires heating technologies like electrical resistive furnaces, gas-flame heated furnaces, or direct torch-flame heating of the brazing area. While the goal is to heat the metal material itself, the use of large heating devices leads to a large loss of energy due to the preheating process. Metals can be heated by technology called induction heating. Any electroconductive material placed in the alternating electromagnetic field becomes heat due to electromagnetic induction phenomena. Using such a technique, while the target materials are electroconductive, the brazing process can be achieved by induction heating. This process is called induction brazing and it was the focus of our work [1][2][3][4].
Brazing is a complex process defined by its physical parameters, such as material properties, process temperatures, time dependence, energy transfer, and technology application. Achieving stable and repeatable results by keeping the quality at high levels is a basic requirement for industrial application. Knowing the influence of each parameter and property gives a strong ability to solve the problem in the best way.
Typical applications of induction brazing are in various industries mostly related with heating, ventilation, and air conditioning (HVAC) systems. Such a process can be successfully used in automotive, home appliances, industrial heating and cooling, and many more areas of the industry manufacture. All these processes are currently performed by highly expensive, energy-consuming, and low-efficiency devices and techniques. Most of them are a difficult subject of automation. The latest achievements in robotics and automation systems correlate well to induction brazing.
The key parameters of the induction brazing process are related with the physics of the materials, their reaction to the electromagnetic energy, and the relation between material properties and the temperature change. To define the process properly, we have to consider the following basic parameters [4][5][6][7]: • Electrical resistivity-metals with relatively higher values acquire more energy from the field. Different metals also have a completely different temperature relation of this parameter.

•
Operating process temperature or brazing temperature-this value is defined by the specifics of the joined metals and the brazing alloy used. Usually, this point or area is specified by the process engineers designing the process. • Time for brazing-the requested time is usually as short as possible. However, this also has a process time for alloy flow between the two metals. So, fast heating is just one side of the coin, while achieving brazing is the other. Usually, the temperature of the joint area has to be kept in a certain range for a certain amount of time to ensure good flow of the brazing alloy.

•
Size of the joined pieces or size of the brazing zone-pipe sizes that have been brazed can vary from less than 1 mm to up to 50 mm or even more. The problem with bigger sizes is the provision of relatively equal heating to the whole of the joint. Occasionally, brazing can be made in steps from two or even more sides (directions).

•
Shape of the assembly-sometimes, the brazing area is not fully directly accessible. This is a source of challenges to define and achieve the proper process, by just performing experiments or tests.
Finding a technical solution of the induction brazing can be done by two main approaches: • Laboratory experiments and industrial test-this approach gives good results and definitely provides experience and data for post-process analysis. The problem is that the complex shapes of the coils and the joint assemblies require time for engineering and manufacturing of complex tools to perform the test. Each new iteration starts from the beginning.

•
Advanced method with combination of real test data and numerical analysis-a strong instrument adds efficiency and effectiveness to the final technical solution. Nowadays, the use of computer hardware and software tools provides a different view over this complex process.
Numerical analysis or just simulation is an engineering application of scientific methods extracted from physics, mathematics, chemistry, and computer technology. Processes like brazing may look simple the first time, but combined with an electromagnetic field, it can become real complex. The aim when performing simulation of some set items and their environment is to simplify the shapes and volumes as much as possible. This makes the task more achievable. For example, if there is an axial symmetry, then 2-D analysis can be easily used. It is less time consuming due to its quick setup and computational time, which is often in a few seconds. On the other hand, 3-D analysis requires more time and resources dedicated to it. Solutions are complex and reaching them is not an easy job. Do we really need to perform such a complex analysis to obtain our goals? The answer to this question comes up by trying to realize the volume of the simple task-the complex shape of the coil, electromagnetic field distribution, current density distribution, resultant heat sources, and some process properties like proper coil cooling.
The common case for induction brazing is to join two pipes of a similar size made of the same or different metal materials. The filler material is also an alloy of several compositions mostly on a metal base. The problem here is how to generate a proper electromagnetic field around the joining area. The use of a simple helical coil is not applicable in many cases because the pipes are part of an assembly or even of a whole device. The only option left is to access the joint from one of the sides. This raises the need of more complex coils to be put on a service.
Modeling is based on a real brazing system with a C-shape coil. The complex design of the coil provides a non-symmetric distribution of the electromagnetic field and the resulting heating pattern. This appears to be a rarely studied matter, with special attention paid to the parameter dependence of the induction coil, which is important for induction power generator load matching. To simulate the brazing process, we need several inputs as the basis of our task: • A 3-D model of the complex coil-it can be constructed by using of a software designing tool like ® Solid Works.
• Some experimental data about the electromagnetic field operating frequency, coil excitation currents, and expected process time-data taken from application test reports provided by © Ultraflex Power Technologies application Lab [4].
• Geometry of the brazing assembly, its size, and materials.

•
Software tool that is capable of performing complex 3-D numerical analysis-® COMSOL Multiphysics.

•
Hardware platform to perform this heavy task.
For the presented application of induction brazing, no studies have been performed related to 3-D modeling of an asymmetric inductor and on this basis to obtain the dependences of the change of R and L during the technological process. Determining the power losses that are related to the cooling of the inductor and the total power required are also included. All these parameters defined by us are basic for the design of power electronic devices that generate the necessary high-frequency energy.
The paper is organized into the following sections: • Description-details of the study inputs combined with the theoretical fundamentals included as methods and approaches.

Description
The equations describing the brazing process, or more generally, the processes in induction heating systems, are well known (see, e.g., [1,2]). Additionally, the numerical methods for their solution were developed years ago. Initially, many authors used relatively simple own codes for their models. This approach was successful in 1-D and 2-D models and lead to many important results being obtained that explain the basic features of the studied systems. At the same time, the simplifications in these models do not permit their application for obtaining the parameters of real systems. The next step-realistic 3-D modelling of existing systems-requires the use of commercial software like Comsol Multiphysics. Of course, the software package is only an instrument for the modelling of systems with a complicated geometry and cannot take the role of the author in the choice of the processes into account, as well as the studied parameters and in the interpretation of the results.

•
Brazing assembly pipes-simulation performed for two cases: Copper and stainless steel AISI 316.

•
Brazing alloy-copper to copper pipes-copper (similar to CuP7). Steel to steel-silver (similar to high silver brazing alloys). The electrical conductivity of the material σ is given by its temperature dependence (Linearized resistivity): where: ρ0-reference resistivity, Ω.m; α-resistivity temperature coefficient, n/a; and Tref-reference temperature, °C.
Coil excitation can be defined in several ways depending on the goals of the simulation requirements. This part also defines the energy source of the system. In our case, the induction brazing system power source consists of a Radio Frequency (RF) generator, a matching circuit station and remote inductor station, know more as a brazing gun. All these parts of the system are suitable to be presented as a current source while during real operation, control is applied to the RF generator output current. For both brazing cases simulated, copper to copper and steel to steel, we used the following parameters (Table 1). The electrical parameters of the simulation were based on real laboratory data experiments and industrial trials.
The complex heat transfer solution for the heated parts (brazing joint and coil) is presented by (2): where: md-material solid density, kg/m 3 ; C-material solid heat capacity at constant pressure, J/(kg•K); k-material solid thermal conductivity, W/(m•K); and The electrical conductivity of the material σ is given by its temperature dependence (Linearized resistivity): where: ρ 0 -reference resistivity, Ω.m; α-resistivity temperature coefficient, n/a; and Coil excitation can be defined in several ways depending on the goals of the simulation requirements. This part also defines the energy source of the system. In our case, the induction brazing system power source consists of a Radio Frequency (RF) generator, a matching circuit station and remote inductor station, know more as a brazing gun. All these parts of the system are suitable to be presented as a current source while during real operation, control is applied to the RF generator output current. For both brazing cases simulated, copper to copper and steel to steel, we used the following parameters (Table 1). The electrical parameters of the simulation were based on real laboratory data experiments and industrial trials.
The complex heat transfer solution for the heated parts (brazing joint and coil) is presented by (2): where: In the term of the heat source: where E is the induced electric field intensity. Energy losses related to the brazing assembly during the heating process can be summarized in two ways: Convection heat flux-energy lost in surrounding air heating-and radiative heat flux.
The boundary conditions on the surface of the coil and of the workpiece are for normal energy flux determined by convection and radiation: where: h-heat transfer coefficient; ε-surface emissivity; σ-Stefan-Boltzmann constant; and T amb -ambient temperature (surrounding air in our case).
In addition, the coil is cooled by water flowing in an internal cooling channel. Following [11], the water cooling is presented as an additional (negative) heat source: where: m t -the water mass flow, kg/s; T w -water inlet temperature; and V-the volume of the water inside the coil.
Geometry meshing as a vital part of the simulation process is presented in Figure 2. Considering a phenomenon like the skin effect and proximity effect, some special attention was focused over the elements of the mesh related to both electroconductive pars-the coil and brazing assembly. In the term of the heat source: where E is the induced electric field intensity. Energy losses related to the brazing assembly during the heating process can be summarized in two ways: Convection heat flux-energy lost in surrounding air heating-and radiative heat flux.
The boundary conditions on the surface of the coil and of the workpiece are for normal energy flux determined by convection and radiation: where: h-heat transfer coefficient; ε-surface emissivity; σ-Stefan-Boltzmann constant; and Tamb-ambient temperature (surrounding air in our case).
In addition, the coil is cooled by water flowing in an internal cooling channel. Following [11], the water cooling is presented as an additional (negative) heat source: where: mt-the water mass flow, kg/s; Tw-water inlet temperature; and V-the volume of the water inside the coil.
Geometry meshing as a vital part of the simulation process is presented in Figure 2. Considering a phenomenon like the skin effect and proximity effect, some special attention was focused over the elements of the mesh related to both electroconductive pars-the coil and brazing assembly.

Validation of the Model
The shape of the coil is complicated, and its inductance cannot be calculated analytically. In order to check the reliability of the obtained results, they were compared with results for axisymmetric coils with a similar size and five windings. First, the analytical formula (6):

Validation of the Model
The shape of the coil is complicated, and its inductance cannot be calculated analytically. In order to check the reliability of the obtained results, they were compared with results for axisymmetric coils with a similar size and five windings. First, the analytical formula (6): was applied for a coil with N = 5, length l = 20 mm, and radius R = 8 mm. The result is inductance of the coil L = 3.15827 × 10 −7 H and resistance R = 0.0013 Ohm (ρ-specific electrical resistivity, S-cross section of the coil conductor). The analytical Formula (6) is valid for an infinite coil and it overestimates the inductance value. The next step is a 2-D axis-symmetric numerical model for a coil with the same parameters. The finite length of the coil leads to a decrease of the inductance to L = 2.2131×10 −7 H while the value of the resistance is the same (R = 0.0013 Ohm). These results were obtained at low frequency. Increasing the frequency to 100 kHz (in the range of operation of the real brazing system) leads to an increase of the resistance to R = 0.0048 Ohm because of the skin effect. In addition, a 3-D model of a helical coil also taking into account the leads was developed. The value of the inductance is in the same range like in the previous cases (L = 2.9154 × 10 −7 H) and because of the increased length of the conductor, the resistance is higher, changing from 0.0018 Ohm at low frequencies to 0.0074 Ohm at 100 kHz. The inductance and the resistance from the model of the brazing system without a workpiece are 2.82 × 10 −7 H and 0.0066 Ohm, respectively. Therefore, the results from the model of the brazing system are in the expected limits.
Experimental results taken from the application test records provided by UltraFlex Power Technologies are presented at Table 2. Tests were performed by the use of its handheld induction brazing system UBraze [15]. It can be noted that the two experiments were performed by using two different coil configurationsthree turns for copper and one turn for stainless steel. In order to better match to our model, we decided to decrease the coil current for the steel-to-steel study while our coil had three turns for both cases [16][17][18].

Results
Obtaining useful results for understanding the brazing process was the main goal of this work. Visualization of the invisible physical phenomena and its time dependence is one of the main advantages of our approach. Such a technique gives us a strong instrument for detailed analysis of the usually hidden aspects of the complex relationship between electromagnetism and heat transfer. Using it for Energies 2020, 13, 3656 7 of 14 analysis of an industrial process like brazing can optimize more of its parameters and significantly improve the performance [18][19][20].
One of the most interesting results of the simulation model is the magnetic flux density distribution across the 3-D space. As it can be expected, there are areas with more concentrated magnetic flux because of the complex coil shape. The 2-D cross-section views are presented in Figure 3. While the magnetic flux is a vector, its direction in a particular moment of time can also be visualized (Figure 4).  The electric current is one of the key parameters of the whole process. Its surface density distribution is visualized on Figure 5. Looking at the coil figure 5b, it is easy to observe the proximity effect of the high frequency electrical current flowing though the coil conductor. The C-shape of the coil defines this complex pattern. From a practical point of view, regions with higher current density can cause potential local overheating. This information can be used for fine adjustment of the coil shape   The electric current is one of the key parameters of the whole process. Its surface density distribution is visualized on Figure 5. Looking at the coil figure 5b, it is easy to observe the proximity effect of the high frequency electrical current flowing though the coil conductor. The C-shape of the coil defines this complex pattern. From a practical point of view, regions with higher current density can cause potential local overheating. This information can be used for fine adjustment of the coil shape  The electric current is one of the key parameters of the whole process. Its surface density distribution is visualized on Figure 5. Looking at the coil Figure 5b, it is easy to observe the proximity effect of the high frequency electrical current flowing though the coil conductor. The C-shape of the coil defines this complex pattern. From a practical point of view, regions with higher current density can cause potential local overheating. This information can be used for fine adjustment of the coil shape [20-22].
Magnetic flux causes a "footprint" of electric current density over the brazing assembly presented on Figure 5a. It can be directly related to the resulting heat source (Figure 6a).
The final resulting heating of the brazing assembly is one that can be observed during real laboratory or industrial tests. However, it is not so easy to be analyzed or even measured as the process times are relatively short. The simulation approach allows more detailed study and estimation at any moment. Figure 6b presents the temperature distribution over the brazing assembly after 6 s of heating. It matches very well to the test results. The electric current is one of the key parameters of the whole process. Its surface density distribution is visualized on Figure 5. Looking at the coil figure 5b, it is easy to observe the proximity effect of the high frequency electrical current flowing though the coil conductor. The C-shape of the coil defines this complex pattern. From a practical point of view, regions with higher current density can cause potential local overheating. This information can be used for fine adjustment of the coil shape [20-22]. Magnetic flux causes a "footprint" of electric current density over the brazing assembly presented on Figure 5a. It can be directly related to the resulting heat source (Figure 6a).  The final resulting heating of the brazing assembly is one that can be observed during real laboratory or industrial tests. However, it is not so easy to be analyzed or even measured as the process times are relatively short. The simulation approach allows more detailed study and estimation at any moment. Figure 6b presents the temperature distribution over the brazing assembly after 6 s of heating. It matches very well to the test results.
The results for the more detailed analysis about our case are presented in the following figures. Using the simulation, we managed to compare two cases common for the industrial practicebrazing of copper-to-copper pipes and brazing of stainless-steel pipes.
Brazing of these two cases requires similar temperatures to be reached in order to melt the brazing alloy ring and provide good flow into the gap between the pipes. Copper pipes are brazed by using an alloy with a melting temperature in the range of 640-720 °C. The requested temperature into the brazing area is around 800 °C. Figure 7a shows that this value was achieved at around 3.5 s after the start. At this point, it is reasonable to slow down the heating (by lowering the power) and try to keep the temperature in the requested range for 2 or 3 s to provide complete melting and flow of the alloy. Stainless-steel pipes are brazed by using an alloy with a temperature range of 780-840 °C. The operating temperature here is around 900 °C. From Figure 7b, we can note that this happens 4.5 s after the start. Again, the temperature has to be kept for 2-3 s. These two figures give us a moment of time when we have to start maintaining the temperatures. The results for the more detailed analysis about our case are presented in the following figures. Using the simulation, we managed to compare two cases common for the industrial practice-brazing of copper-to-copper pipes and brazing of stainless-steel pipes.
Brazing of these two cases requires similar temperatures to be reached in order to melt the brazing alloy ring and provide good flow into the gap between the pipes. Copper pipes are brazed by using an alloy with a melting temperature in the range of 640-720 • C. The requested temperature into the brazing area is around 800 • C. Figure 7a shows that this value was achieved at around 3.5 s after the start. At this point, it is reasonable to slow down the heating (by lowering the power) and try to keep the temperature in the requested range for 2 or 3 s to provide complete melting and flow of the alloy. Stainless-steel pipes are brazed by using an alloy with a temperature range of 780-840 • C. The operating temperature here is around 900 • C. From Figure 7b, we can note that this happens 4.5 s after the start. Again, the temperature has to be kept for 2-3 s. These two figures give us a moment of time when we have to start maintaining the temperatures.
The acquisition of power from the magnetic field is strongly affected by the material properties. As it is presented in Figure 8, it is also time dependent as the temperature significantly changes during the heating. The performance of the two materials is different and it needs to be taken into consideration. Here, it is the main difference between induction heated brazing and the other contact or convection methods. The magnetic field needs to have enough power to provide stable and reliable heating. From these figures, we can estimate the required power needed to be transferred to the brazing assemblies. The main difference about the power required for the brazing of each material is caused by their electrical conductivity and temperature dependence of this parameter. after the start. At this point, it is reasonable to slow down the heating (by lowering the power) and try to keep the temperature in the requested range for 2 or 3 s to provide complete melting and flow of the alloy. Stainless-steel pipes are brazed by using an alloy with a temperature range of 780-840 °C. The operating temperature here is around 900 °C. From Figure 7b, we can note that this happens 4.5 s after the start. Again, the temperature has to be kept for 2-3 s. These two figures give us a moment of time when we have to start maintaining the temperatures. The acquisition of power from the magnetic field is strongly affected by the material properties. As it is presented in Figure 8, it is also time dependent as the temperature significantly changes during the heating. The performance of the two materials is different and it needs to be taken into consideration. Here, it is the main difference between induction heated brazing and the other contact or convection methods. The magnetic field needs to have enough power to provide stable and reliable The next two figures (Figure 9a,b) are the essential estimation of the induction brazing system. Here, we can see the required electrical power that needs to be applied to the induction coil in both cases of brazing. As expected, it is quite different: Copper requires significantly more power than stainless steel. This is mostly caused by the different coil current excitation. A similar heating effect for both materials is achieved by using significantly different coil currents-copper is heated by the use of 700 A, while for the stainless steel, we used 350 A. The heat power losses are presented in Figure 10. Only the copper-to-copper case is visualized for reference. It can be estimated that the sum of the heating losses is a significant part of the transferred energy, at around 20%. The next two figures (Figure 9a,b) are the essential estimation of the induction brazing system. Here, we can see the required electrical power that needs to be applied to the induction coil in both cases of brazing. As expected, it is quite different: Copper requires significantly more power than stainless steel. This is mostly caused by the different coil current excitation. A similar heating effect for both materials is achieved by using significantly different coil currents-copper is heated by the use of 700 A, while for the stainless steel, we used 350 A. The next two figures (Figure 9a,b) are the essential estimation of the induction brazing system. Here, we can see the required electrical power that needs to be applied to the induction coil in both cases of brazing. As expected, it is quite different: Copper requires significantly more power than stainless steel. This is mostly caused by the different coil current excitation. A similar heating effect for both materials is achieved by using significantly different coil currents-copper is heated by the use of 700 A, while for the stainless steel, we used 350 A. The heat power losses are presented in Figure 10. Only the copper-to-copper case is visualized for reference. It can be estimated that the sum of the heating losses is a significant part of the transferred energy, at around 20%. The heat power losses are presented in Figure 10. Only the copper-to-copper case is visualized for reference. It can be estimated that the sum of the heating losses is a significant part of the transferred energy, at around 20%. At the end, a heat loss dissipated into the induction coil is shown in Figure 11. Again, only the copper-to-copper case is presented. This estimation is important for induction cooling system evaluation. Figure 11. Cooling power-coil (copper).

Coil Analysis
In this section, the results for the changes of the coil inductance and resistance with time are presented.
The coil and the workpiece could be presented as a transformer, where the workpiece plays the role of a one turn secondary coil ( Figure 12). Eg is the electromotive force (EMF) of the generator (more precisely, this is the voltage at the ends of the coil), R1 and L1 are the resistance and the inductance of the coil, and R2 and L2 are the resistance and the inductance of the workpiece, respectively.

R1 R2
L1 L2 At the end, a heat loss dissipated into the induction coil is shown in Figure 11. Again, only the copper-to-copper case is presented. This estimation is important for induction cooling system evaluation. At the end, a heat loss dissipated into the induction coil is shown in Figure 11. Again, only the copper-to-copper case is presented. This estimation is important for induction cooling system evaluation. Figure 11. Cooling power-coil (copper).

Coil Analysis
In this section, the results for the changes of the coil inductance and resistance with time are presented.
The coil and the workpiece could be presented as a transformer, where the workpiece plays the role of a one turn secondary coil ( Figure 12). Eg is the electromotive force (EMF) of the generator (more precisely, this is the voltage at the ends of the coil), R1 and L1 are the resistance and the inductance of the coil, and R2 and L2 are the resistance and the inductance of the workpiece, respectively.

Coil Analysis
In this section, the results for the changes of the coil inductance and resistance with time are presented. The coil and the workpiece could be presented as a transformer, where the workpiece plays the role of a one turn secondary coil ( Figure 12). Eg is the electromotive force (EMF) of the generator (more precisely, this is the voltage at the ends of the coil), R 1 and L 1 are the resistance and the inductance of the coil, and R 2 and L 2 are the resistance and the inductance of the workpiece, respectively.
Assuming time dependence of the currents and voltages are of the type exp(jωt), from Kirchhoff's second law, we obtain: where I 1 and I 2 are the currents in the coil and in the workpiece, and E 12 and E 21 are the induced EMFs: where M is a mutual inductance. Expressing I 2 from (8): The equation for E g is obtained in the form of: From (10), the resistance and the inductance of the system are: with: Therefore, introducing the workpiece into the coil leads to an increase of its resistance and decrease of its inductance. Figure 11. Cooling power-coil (copper).

Coil Analysis
In this section, the results for the changes of the coil inductance and resistance with time are presented.
The coil and the workpiece could be presented as a transformer, where the workpiece plays the role of a one turn secondary coil ( Figure 12). Eg is the electromotive force (EMF) of the generator (more precisely, this is the voltage at the ends of the coil), R1 and L1 are the resistance and the inductance of the coil, and R2 and L2 are the resistance and the inductance of the workpiece, respectively. The time dependence of the resistance and of the inductance of the system (Figure 4) is determined by the dependence of the material conductivity with respect to the temperature. Because of the water cooling, the increase of R1 with time is relatively weak. The behavior of ∆R is more complicated (12). At R 2 = 0 and R 2 = ∞, its value is zero. The maximum of ∆R is at R 2 = ωL 2 . At the conditions of our model, 0 < R 2 < ωL 2 and ∆R increases with the resistance R 2 ( Figure 5).
With the temperature increases and the conductivity decreases, the skin depth in the coil and in the workpiece increases and the electric current penetrates deeper in the metal. For the coil, this means an increase of its cross section and a relatively small increase of the inductance L 1 (Figure 13). In opposite, for the workpiece, the effective radius becomes smaller and L 2 decreases. In combination with the decrease of the mutual inductance M = √ L 1 L 2 and increase of R 2 , this leads to lower ∆L ( Figure 14) and a rise of L ( Figure 13) with the temperature.
The changes of the coil parameters during the brazing process are of significant practical interest as they affect the matching of the generator and the coil and the effectiveness of the power transfer. The changes of the coil parameters during the brazing process are of significant practical interest as they affect the matching of the generator and the coil and the effectiveness of the power transfer.

Conclusions
There are FEM brazing models, but they are usually on axial-symmetric coils. We presented that this modeling method can also be used for windings with a more complex shape. We presented, for the first time, the change of the parameters of the winding during the process (inductance and resistance) as well as brazing of different materials. These results are important for determining the parameters in the design of the power electronic converter. The main goal was to realize different combinations of brazing materials with one converter. Modeling helped us determine the limits of change in the power and frequency so that we could design the converter.
During our study, we managed to achieve a successful 3-D simulation of the induction brazing process. Graphical and numerical representation of the key parameters and their time dependence provided a better look over the complex relation between the electromagnetic field and the resulting material heating. Estimation of the required electrical power directly refers to the design consideration of the induction brazing system. Complex coil design was estimated, and the further analysis and optimizations were possible based on the presented modeling. Power heat losses were estimated and visualized for both main components of the brazing assembly with its convection and radiation emissions and the induction coil with the heat dissipated into the cooling water.
The presented approach and the obtained results can be useful for application and process development engineers in the area of industrial manufacturing. Using the tools of the simulation, induction brazing system components can be evaluated in detail. The model adds abilities to solve practical problems, such as: The changes of the coil parameters during the brazing process are of significant practical interest as they affect the matching of the generator and the coil and the effectiveness of the power transfer.

Conclusions
There are FEM brazing models, but they are usually on axial-symmetric coils. We presented that this modeling method can also be used for windings with a more complex shape. We presented, for the first time, the change of the parameters of the winding during the process (inductance and resistance) as well as brazing of different materials. These results are important for determining the parameters in the design of the power electronic converter. The main goal was to realize different combinations of brazing materials with one converter. Modeling helped us determine the limits of change in the power and frequency so that we could design the converter.
During our study, we managed to achieve a successful 3-D simulation of the induction brazing process. Graphical and numerical representation of the key parameters and their time dependence provided a better look over the complex relation between the electromagnetic field and the resulting material heating. Estimation of the required electrical power directly refers to the design consideration of the induction brazing system. Complex coil design was estimated, and the further analysis and optimizations were possible based on the presented modeling. Power heat losses were estimated and visualized for both main components of the brazing assembly with its convection and radiation emissions and the induction coil with the heat dissipated into the cooling water.
The presented approach and the obtained results can be useful for application and process development engineers in the area of industrial manufacturing. Using the tools of the simulation, induction brazing system components can be evaluated in detail. The model adds abilities to solve practical problems, such as:

Conclusions
There are FEM brazing models, but they are usually on axial-symmetric coils. We presented that this modeling method can also be used for windings with a more complex shape. We presented, for the first time, the change of the parameters of the winding during the process (inductance and resistance) as well as brazing of different materials. These results are important for determining the parameters in the design of the power electronic converter. The main goal was to realize different combinations of brazing materials with one converter. Modeling helped us determine the limits of change in the power and frequency so that we could design the converter.
During our study, we managed to achieve a successful 3-D simulation of the induction brazing process. Graphical and numerical representation of the key parameters and their time dependence provided a better look over the complex relation between the electromagnetic field and the resulting material heating. Estimation of the required electrical power directly refers to the design consideration of the induction brazing system. Complex coil design was estimated, and the further analysis and optimizations were possible based on the presented modeling. Power heat losses were estimated and visualized for both main components of the brazing assembly with its convection and radiation emissions and the induction coil with the heat dissipated into the cooling water.
The presented approach and the obtained results can be useful for application and process development engineers in the area of industrial manufacturing. Using the tools of the simulation, induction brazing system components can be evaluated in detail. The model adds abilities to solve practical problems, such as: Coil shape and size estimation and basic optimization; Power estimation of induction generator; Process time calculation; and Cooling system power and size estimation.
Author Contributions: D.G., K.T. and N.H. were involved in the full process of producing this paper including conceptualization, methodology, modelling, validation, visualization, and preparing the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research was funded by the European Regional Development Fund within the Operational Program "Science and Education for Smart Growth 2014-2020" under the Project CoE "National center of mechatronics and clean technologies "BG05M2OP001-1.001-0008".

Conflicts of Interest:
The authors declare no conflict of interest.