Numerical Simulation of TIG Welding Arc Based on Multi-Physics Field Coupling

Literature Overview

This research, conducted by Guo Chaobo, Cui Lulu, Tao Kai, and Wang Huimin from Henan Institute of Technology and the Henan Province Engineering Technology Research Center for Metal Material Modification, was published in 2020 and supported by the Henan Province Science and Technology Project (182102210260) and the Higher Education Key Research Project Basic Research Special Plan (20B430003). The study presents a comprehensive numerical simulation of the TIG welding arc based on multi-physics field coupling, integrating electromagnetic, fluid dynamics, thermal, and plasma physics phenomena into a unified computational model. The work represents a significant advancement in understanding the complex physics of the TIG welding arc and provides valuable insights for process optimization and defect prevention.

Core Technical Content

The TIG welding arc is a complex plasma phenomenon involving multiple coupled physical fields:

Electromagnetic field: Determines current density distribution, magnetic field generation, and Lorentz force acting on the plasma
Fluid dynamics: Governs plasma flow patterns, gas velocity fields, and momentum transfer
Thermal field: Controls temperature distribution, heat transfer mechanisms, and energy balance
Plasma physics: Describes ionization states, species transport, and radiative transfer

Traditional numerical models often treat these fields independently, solving one field at a time and using the results as boundary conditions for the next field. This sequential approach fails to capture the strong coupling between fields, leading to inaccurate predictions of arc behavior and heat transfer. The multi-physics coupled model presented in this study solves all fields simultaneously, allowing for a more accurate representation of the arc physics.

The governing equations for the coupled model include:

Physical Field	Governing Equation	Key Parameters
Electromagnetic	Maxwell's equations	Current density, magnetic field
Fluid dynamics	Navier-Stokes equations	Velocity, pressure, viscosity
Thermal	Energy equation	Temperature, heat flux, thermal conductivity
Plasma	Species transport equations	Ion density, electron density, ionization rate
Momentum	Momentum conservation	Lorentz force, pressure gradient

The numerical solution employs a finite element method with adaptive mesh refinement to capture the sharp gradients near the electrode surfaces and the arc root. The model incorporates realistic boundary conditions, including electrode heat flux distributions, gas flow injection conditions, and radiative heat transfer at the arc boundaries.

Simulation Results and Analysis

The numerical simulation reveals several important characteristics of the TIG welding arc:

Current density distribution: The current density is highly concentrated at the cathode (tungsten) tip, with peak values exceeding 10^9 amperes per square meter. The current spreads gradually along the arc axis, reaching the anode (workpiece) with a more uniform distribution. This non-uniform current distribution explains the asymmetric heat input observed in TIG welds.

Temperature distribution: The arc temperature reaches a maximum of approximately 20,000 to 25,000 Kelvin at the arc root, decreasing rapidly along the arc axis to 10,000 to 15,000 Kelvin at the cathode. The workpiece surface temperature at the arc center reaches 3,000 to 4,000 Kelvin, with a heat-affected zone extending 5 to 10 millimeters from the weld centerline.

Flow patterns: The plasma flow exhibits a complex pattern with upward convection along the arc axis and downward flow near the electrode surfaces. The flow velocity reaches maximum values of 50 to 100 meters per second near the arc root, creating significant momentum transfer to the weld pool.

Radiative heat transfer: Radiation accounts for approximately 30 to 50 percent of the total heat input to the workpiece, depending on arc current and shielding gas composition. The radiative heat flux is concentrated in a narrow zone directly below the arc, with a Gaussian-like distribution.

Engineering Practice Implications

The numerical simulation results have direct implications for TIG welding process optimization:

Heat input control: The simulation provides accurate predictions of heat input distribution, enabling engineers to optimize welding parameters for specific applications. For pressure vessel fabrication, the heat input profile can be used to predict distortion patterns and design appropriate fixturing strategies.
Defect prediction: The simulation identifies regions of high cooling rates and thermal stresses that are susceptible to cracking. By adjusting welding parameters, such as travel speed and current, the simulation can predict and prevent defects such as hot cracking, cold cracking, and porosity.
Process window determination: The simulation can be used to establish process windows for different materials and geometries. For example, the simulation can identify the maximum allowable current for a given tungsten electrode diameter to prevent electrode melting and tungsten contamination.
Shielding gas optimization: The simulation predicts the effect of shielding gas composition and flow rate on arc behavior and heat transfer. This information can be used to optimize gas selection for specific welding applications, such as welding reactive metals or thick-section materials.

The simulation also provides insights into the effects of arc length, electrode angle, and workpiece geometry on weld quality. For example, increasing the arc length from 3 to 6 millimeters reduces the peak heat flux at the workpiece surface by approximately 20 percent but increases the heat-affected zone width by 15 percent. This trade-off must be carefully managed in pressure vessel fabrication to achieve the desired balance between penetration and distortion.

Key Technical Points and Reflections

The most significant contribution of this study is the development of a comprehensive multi-physics coupled model that accurately captures the complex physics of the TIG welding arc. The model provides a powerful tool for process optimization and defect prevention, reducing the need for extensive experimental trials and enabling the design of welding processes for novel materials and geometries.

However, the study also acknowledges the limitations of numerical simulation. The model relies on several simplifying assumptions, such as local thermodynamic equilibrium (LTE) and ideal gas behavior, which may not hold under all welding conditions. The computational cost of the coupled simulation is significant, requiring high-performance computing resources and limiting its application to real-time process control. Future work should focus on developing reduced-order models that capture the essential physics while reducing computational cost, enabling integration into process monitoring and control systems.

Study Insights and Implications

This research demonstrates the power of multi-physics numerical simulation in understanding and optimizing the TIG welding process. For engineers involved in bimetal pressure vessel fabrication, the simulation provides a valuable tool for predicting weld behavior and optimizing process parameters. The ability to simulate complex welding scenarios, such as multi-pass welds, welding of dissimilar metals, and welding under constrained conditions, enables the design of welding procedures that minimize defects and maximize productivity. The key insight is that a deep understanding of arc physics, enabled by numerical simulation, is essential for achieving consistent weld quality in challenging applications.