Three-Dimensional TIG Welding Molten Pool Flow and Thermal Field Numerical Analysis Under Penetration Conditions Study Note

Literature Overview

This foundational study, published in 1992 in Acta Metallurgica Sinica (金属学报) by Wu Chuansong, Cao Zhening, and Wu Lin from the Welding Department of Harbin Institute of Technology, represents one of the earliest comprehensive three-dimensional numerical analyses of TIG welding molten pool dynamics under full penetration conditions. The work was supported by the National Natural Science Foundation of China. At a time when two-dimensional modeling was the norm, this study's adoption of full three-dimensional analysis marked a significant advancement in computational welding science.

Core Technical Content

Model Formulation

The numerical model developed in this study incorporates the following governing equations:

Equation	Description	Significance
Continuity equation	∇·v = 0	Mass conservation in the molten pool
Momentum equation	ρ(Dv/Dt) = -∇P + ∇·τ + F_electromagnetic + F_buoyancy + F_surface_tension	Fluid dynamics with all major driving forces
Energy equation	ρc_p(DT/Dt) = ∇·(k∇T) + Q_arc	Heat transfer with arc heat source
Free surface boundary condition	Kinematic and dynamic conditions at the pool surface	Accurate representation of the pool surface shape

Arc Heat Source Model

The study employs a Gaussian-type heat source distribution on the pool surface:

Q(r) = (2Q₀/πr₀²) exp(-2r²/r₀²)

where Q₀ is the total heat input, r is the radial distance from the arc center, and r₀ is the heat source radius. This model was subsequently refined by various researchers, but remains a widely used baseline for TIG welding simulations.

Boundary Conditions

The model incorporates the following boundary conditions:

Pool surface: Free surface with surface tension gradient (Marangoni convection), arc pressure, and heat flux.
Solid-liquid interface: Isothermal condition at the melting temperature (T_m), with Stefan condition for latent heat absorption.
Solid boundary: Adiabatic or convective cooling, depending on the specific condition modeled.
Symmetry planes: Symmetry conditions applied to reduce computational domain where applicable.

Key Numerical Results

Thermal Field Characteristics

Parameter	Value	Description
Maximum temperature	~2200–2500 °C	At the pool surface center, near the arc
Pool depth	3–6 mm	For typical TIG parameters (150–200 A)
Pool width	8–12 mm	Surface width at the pool boundary
Cooling rate at fusion boundary	50–200 °C/s	Critical for solidification microstructure
Thermal gradient at fusion boundary	500–2000 K/mm	Drives columnar dendrite growth

Flow Field Characteristics

Parameter	Value	Description
Maximum flow velocity	1.0–2.5 m/s	Near the pool surface
Dominant flow direction	Downward at center, upward at edges	Large convection cell
Pool surface depression	0.5–2.0 mm	Due to electromagnetic force and arc pressure
Flow pattern	Single dominant convection cell	In the studied parameter range

Penetration Analysis

The study specifically addresses full penetration conditions, where the molten pool extends through the entire plate thickness. Key findings include:

The penetration depth increases with current and decreases with travel speed, following the expected relationship.
At full penetration, the pool shape transitions from a shallow ellipsoid to a deep, narrow keyhole-like profile.
The flow velocity at the penetration tip is significantly reduced compared to the pool surface, indicating a quasi-static condition at the deepest point.
The thermal field shows a steep gradient at the root of the penetration, with rapid cooling that influences solidification microstructure.

Comparison with Subsequent Research

This 1992 study laid the groundwork for subsequent numerical welding research. Key developments that followed include:

Development	Year	Significance
Two-phase flow model (VOF)	1990s	More accurate free surface tracking
Moving heat source model	1990s–2000s	Dynamic simulation of welding process
Multi-phase solidification model	2000s	Coupled thermal-fluid-solidification analysis
Keyhole model for laser welding	1990s–2000s	Extended to high-energy-density processes
CFD-based welding simulation	2000s–present	Commercial software packages

The fundamental physics captured in this 1992 study remain valid and form the basis of modern welding simulation codes. The key contribution was demonstrating that three-dimensional analysis is feasible and provides significantly more accurate predictions than two-dimensional models, particularly for penetration behavior and pool shape.

Engineering Practice Implications

Process Optimization

The numerical results provide a basis for rational process optimization:

Penetration control: Adjust current and travel speed to achieve the desired penetration depth without excessive heat input.
Weld profile control: Optimize parameters to achieve a desirable weld profile (narrow, deep penetration for structural applications; wider profile for corrosion resistance).
Cooling rate control: Manage heat input and preheat to control the cooling rate at the fusion boundary, influencing solidification microstructure and mechanical properties.
Distortion prediction: Use thermal field predictions to estimate welding distortion and design appropriate fixturing.

Quality Control

The understanding of thermal and flow fields enables more effective quality control:

Radiographic testing: Predicted pool shape helps interpret radiographic images and identify potential defect locations.
Ultrasonic testing: Understanding the solidification pattern aids in UT technique development and interpretation.
Mechanical property prediction: Cooling rate predictions enable estimation of hardness and strength distributions across the weld.

Key Questions and Reflections

Despite its foundational significance, this study has limitations that are important to recognize:

The heat source model is relatively simple (Gaussian distribution), and more sophisticated models (double-ellipsoid, moving source) have been developed since.
The model does not account for gas metal transfer or spatter, which are relevant for some welding conditions.
The material properties are assumed to be constant or linearly varying with temperature, which may not accurately represent the complex behavior of real alloys.
The study focuses on steady-state conditions, while real welding involves transient effects at the start and end of each weld.

Nevertheless, the fundamental physics captured in this study remain valid. The governing equations, the inclusion of all major flow-driving forces, and the three-dimensional approach represent the correct methodology for welding process analysis. Modern commercial codes (such as ANSYS Fluent, OpenFOAM, and specialized welding simulation packages) build upon the foundations laid by this research.

Study Insights and Implications

This 1992 study represents a milestone in computational welding science. At a time when computational resources were severely limited, the authors demonstrated that three-dimensional numerical analysis of welding processes was feasible and provided valuable insights into molten pool dynamics. The study's emphasis on full penetration conditions is particularly relevant for pressure vessel fabrication, where full penetration welds are required for structural integrity.

For today's engineers, the value of this study lies not only in its specific results but also in its methodology. The approach of coupling thermal, fluid, and free surface analysis provides a comprehensive framework for understanding and optimizing welding processes. Modern engineers can build upon this foundation using more advanced computational tools and more sophisticated models, but the fundamental principles remain the same.

The study also underscores the importance of understanding the physics of welding processes. Without a fundamental understanding of the thermal and flow fields, process optimization remains empirical and limited. Numerical simulation provides a powerful tool for process development, but its value is maximized when combined with experimental validation and a deep understanding of the underlying physics. For engineers involved in pressure vessel fabrication and quality control, this study provides a valuable foundation for understanding weld quality and optimizing welding procedures.