Numerical Calculation of Anode Current Density in DC TIG Welding Arc - Literature Study Note
Literature Overview
The paper by Fan Honggang, Yang Jun, Shi Yaowu, and Lei Yongping, published in the Journal of Xi'an Jiaotong University (Vol. 31, No. 4, 1997, pp. 77–82), presents a numerical study of the anode current density distribution in a DC TIG (tungsten inert gas) welding arc. The authors developed an axisymmetric two-dimensional numerical model that accounts for the cathode geometry and solves the magnetohydrodynamic (MHD) equations governing the arc plasma. The model was validated by comparing the computed arc column temperature distribution with experimental measurements, and good agreement was reported. The paper then uses the validated model to analyze the effects of welding current, tungsten electrode tip cone angle, and arc length on the anode current density distribution. The classification code TG444 confirms the focus on arc welding process technology. This work is significant because the anode current density distribution directly determines the heat input profile on the workpiece, which in turn controls the weld pool geometry, penetration depth, and HAZ width.
Numerical Model and Methodology
The numerical model is based on the following assumptions and governing equations:
| Aspect | Description |
|---|---|
| Geometry | Axisymmetric 2D (cylindrical coordinates r, z) |
| Cathode | Tungsten electrode with conical tip geometry |
| Anode | Flat workpiece surface (infinite plate approximation) |
| Flow regime | Laminar, incompressible |
| Plasma behavior | Ideal MHD, local thermodynamic equilibrium (LTE) |
| Governing equations | Continuity, momentum, energy, magnetic field (Maxwell's equations), current conservation |
| Solution algorithm | SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) |
| Mesh | Structured, non-uniform, refined near electrode tips |
| Boundary conditions | Fixed temperature at cathode and anode surfaces, no-slip velocity, prescribed current density |
The SIMPLE algorithm is a well-established finite volume method for solving coupled velocity-pressure systems in fluid dynamics. Its application to arc plasma MHD is appropriate because the arc is a low-Mach-number flow where pressure and velocity are strongly coupled. The axisymmetric assumption simplifies the problem from 3D to 2D, reducing computational cost while retaining the essential physics. The non-axisymmetric effects (such as arc wandering due to magnetic field asymmetries) are not captured, but for steady-state analysis, the axisymmetric model provides a good approximation of the average behavior.
The model includes both the cathode region (near the tungsten tip) and the arc column region. The cathode region is critical because it is where the current density is concentrated and where the electric field is highest. The current density distribution at the anode is the output of interest, as it determines the heat flux profile on the workpiece.
Analysis of Key Parameters
The paper examines the effects of three key parameters on the anode current density distribution:
Welding Current
As the welding current increases, the total current density at the anode increases proportionally. However, the distribution shape also changes: at higher currents, the current density becomes more concentrated near the arc center, resulting in a narrower and deeper weld pool. This is because the increased current enhances the Lorentz force (J × B), which compresses the arc column and drives the plasma flow toward the arc axis. The peak current density can be estimated from the total current and the effective current-carrying area at the anode.
Tungsten Electrode Tip Cone Angle
The cone angle of the tungsten electrode tip significantly affects the current density distribution. A sharp tip (small cone angle, e.g., 60° included angle) concentrates the current into a smaller area, producing a higher peak current density and a narrower current distribution. This results in deeper penetration and a narrower weld pool. A blunt tip (large cone angle, e.g., 120° included angle) spreads the current over a larger area, producing a lower peak current density and a wider current distribution. This results in shallower penetration and a wider weld pool. The cone angle is therefore a critical parameter for controlling weld geometry.
Arc Length
The arc length (distance between the tungsten tip and the workpiece surface) affects the current density distribution by changing the arc column geometry and the magnetic field configuration. A short arc (1–2 mm) produces a highly concentrated current density with a sharp peak, while a long arc (4–6 mm) produces a more spread-out distribution with a lower peak. The short arc also provides better shielding gas coverage, reducing the risk of atmospheric contamination. However, a very short arc increases the risk of tungsten inclusion due to electrode-workpiece contact.
| Parameter | Effect on Anode Current Density | Effect on Weld Geometry |
|---|---|---|
| Increased welding current | Higher peak density, narrower distribution | Deeper penetration, narrower HAZ |
| Smaller cone angle | Higher peak density, more concentrated | Deeper penetration, narrower bead |
| Shorter arc length | Higher peak density, sharper peak | Deeper penetration, better gas coverage |
Validation and Experimental Comparison
The model was validated by comparing the computed arc column temperature distribution with experimental measurements obtained using optical pyrometry or spectrometric methods. The agreement was reported to be good, which provides confidence in the model's ability to predict the anode current density distribution. However, the validation is limited to the arc column temperature; the anode current density itself was not directly measured experimentally. Direct measurement of the anode current density is challenging because it requires a non-intrusive technique that does not disturb the arc, such as magnetic field mapping or infrared thermography of the workpiece surface.
Integration with Engineering Practice
The numerical model and the parameter analysis described in this paper have direct implications for welding procedure development and optimization. Engineers can use the model to predict the weld geometry for different parameter combinations, reducing the number of trial welds required during WPS development. The following practical guidelines can be derived:
- For deep penetration welds (e.g., root pass of a multi-pass weld), use a sharp tungsten tip (60° cone angle), short arc length (1–2 mm), and high current density.
- For wide, shallow welds (e.g., cap pass of a multi-pass weld), use a blunt tungsten tip (120° cone angle), longer arc length (3–4 mm), and lower current density.
- For thin-wall welding, use a smaller current and a shorter arc to minimize heat input and prevent burn-through.
- For thick-wall welding, use a larger current and a sharp tungsten tip to achieve sufficient penetration.
The model can also be used to predict the HAZ width and the peak temperature at the weld toe, which are critical for assessing the risk of distortion, cracking, and mechanical property degradation. For example, if the model predicts that the peak temperature at the weld toe exceeds the austenitization temperature (approximately 727°C for low-carbon steel), there is a risk of grain coarsening in the HAZ, which can reduce impact toughness.
Key Questions and Reflections
The model presented in this paper is a steady-state, axisymmetric analysis. In practice, the welding arc is dynamic, with fluctuations in current, arc length, and gas flow that cause the current density distribution to vary with time. A transient model would be more representative of actual welding conditions, but it would require significantly more computational resources. Additionally, the model assumes LTE, which may not hold in the near-electrode regions where the plasma is non-equilibrium. A multi-fluid model or a kinetic model would be more accurate but computationally prohibitive for routine engineering use.
Another limitation is the absence of the workpiece in the model. The anode is treated as an infinite flat plate, which does not account for the thermal conductivity of the workpiece, the effect of the weld pool on the current distribution, or the interaction between the arc and the liquid metal in the weld pool. A coupled arc-workpiece model would provide a more complete picture, but it is beyond the scope of the present study.
Study Insights and Implications
This paper demonstrates the power of numerical modeling in understanding and optimizing welding processes. The MHD model provides insights into the physics of the welding arc that are difficult to obtain from experimental measurements alone. The parameter analysis provides practical guidelines for welding procedure development. For engineers, the key takeaway is that the tungsten electrode geometry, arc length, and welding current are the primary control parameters for the anode current density distribution, and therefore for the weld geometry and quality. A deeper understanding of these relationships enables more informed decision-making in welding process design.
Zhuojin Pipe Fitting Co., Ltd