Zhuojin Pipe Fitting Co., LtdNumerical Analysis of Heat and Mass Transfer in TIG Arc Welding
Literature Overview
The paper by Fan Ding, Chen Jianhong, and Uoura Makoto, published in the Journal of Mechanical Engineering in 1998 (Vol. 34, No. 2, pp. 39-45), represents one of the foundational Chinese contributions to computational modeling of the TIG welding arc. The work establishes a comprehensive mathematical framework for coupled heat and mass transfer phenomena within the welding arc plasma column. At a time when computational fluid dynamics (CFD) was still emerging in the Chinese welding research community, this study demonstrated that rigorous numerical simulation could reproduce arc behavior with fidelity sufficient for engineering application.
Core Methodology and Model Architecture
The authors constructed a complete set of governing equations for the TIG arc, incorporating conservation of mass, momentum, energy, and charge. The model treats the arc as a magnetohydrodynamic (MHD) fluid, where the electric current density is computed directly from Ohm's law rather than being prescribed as a boundary condition. This is a critical modeling decision, as it allows the current density distribution to self-consistently couple with the flow field and temperature field.
The turbulence closure employs the k-epsilon model, a two-equation turbulence model that has proven robust for high-speed, high-temperature gas flows. The numerical implementation was carried out in FORTRAN on a SUN workstation, which reflects the computing environment of the late 1990s. The simulation was conducted across a current range of 100 A to 300 A, covering the practical operating window for most industrial TIG welding applications on carbon steel and alloy pipe.
Key Technical Findings
The numerical results provide several physically meaningful outputs:
| Parameter | Typical Value / Observation | Engineering Significance |
|---|---|---|
| Current range simulated | 100 A to 300 A | Covers thin-wall pipe to heavy-wall fitting welding |
| Turbulence model | k-epsilon | Captures arc column turbulence and jet behavior |
| Temperature distribution | Peak near cathode, decreasing radially | Explains cathode spot erosion and tungsten tip geometry effects |
| Flow field distribution | Downward axial jet with radial expansion | Determines arc force on weld pool and penetration characteristics |
| Current density distribution | Concentrated at arc root, diverging toward anode | Governs energy deposition pattern on workpiece |
The agreement between simulated and experimental results was reported to be satisfactory, particularly for temperature distribution and flow velocity profiles. This validation is essential, as it confirms that the model assumptions and numerical discretization adequately represent the real physics.
Integration with Engineering Practice
From a pipe welding perspective, the arc temperature and flow field distributions directly influence weld penetration depth, bead width, and dilution ratio. For seamless pipe welding and pipe fitting fabrication, understanding how current density concentrates at the arc root allows engineers to predict and control the fusion zone geometry. The numerical approach also provides a means to optimize welding parameters without extensive trial-and-error on actual pipe specimens.
In practice, when welding thin-walled stainless steel pipe or alloy pipe fittings, the concentrated heat input at the arc root can lead to excessive burn-through if not properly managed. The model's prediction of temperature gradients near the cathode helps explain why tungsten electrode angle and stick-out distance are critical parameters. Similarly, the flow field analysis provides insight into gas shielding effectiveness, which is vital for preventing nitrogen pickup in austenitic stainless steel pipe welds.
Key Questions and Reflections
Several questions arise from studying this work. The k-epsilon model, while robust, may not fully capture the compressible, high-temperature behavior of arc plasma, particularly near the arc root where pressure and temperature gradients are extreme. The two-dimensional axisymmetric assumption simplifies the problem significantly but may overlook three-dimensional effects such as arc wander, which is a well-documented phenomenon in TIG welding. Furthermore, the model does not explicitly account for electrode evaporation and metal vapor transport, which are known to affect arc stability and weld composition.
Nevertheless, the study laid important groundwork for subsequent CFD-based welding research in China. The methodology can be extended to simulate arc behavior in constrained geometries such as pipe internal welding or back-side welding of pipe girth welds, where gas flow and arc stability present unique challenges.
Study Insights and Implications
The most valuable takeaway from this paper is the demonstration that a self-consistent MHD model, where current density is computed rather than imposed, yields physically realistic arc behavior. This modeling philosophy has since become standard in computational welding research. For engineering practice, the results reinforce the importance of controlling arc geometry through electrode preparation, stick-out, and travel speed. The work also highlights the value of numerical simulation as a complementary tool to experimental welding trials, particularly when the cost of trial welds on expensive alloy pipe or when testing conditions are difficult to replicate in the laboratory.